login name
ITA english

Article RSE 15002151

A 3D SPH integrated model for granular flows with the transport of solid bodies: model coupling and enhanced boundary treatment based on surface wall elements


10th International SPHERIC Workshop 16-18, Giugno-2015.

Ask material

A. Amicarelli (RSE SpA), G. Agate (RSE SpA)

GRID RESILIENCE 2015 - Security and vulnerabilities of the electrical system

This study preliminarily integrates several numerical schemes into a 3D SPH code: the SPH model of Amicarelli et al. (2013, IJNME) with discrete surface wall elements, a 3D SPH model for bed-load transport (Amicarelli et al., submitted) and a 3D SPH scheme for the transport of solid bodies in fluid flows (Amicarelli et al., 2015, CAF).aries (DB-SPH).

This study aims at integrating three SPH numerical models into a 3D SPH code. The first model simulates bed-loadtransport phenomena. The second model represents the transport of solid bodies in free surface flows. The third modelis based on a boundary treatment with discrete surface wall elements (“DB-SPH” method, just for simplicity ofnotation).

So far, the model couplings mainly involve the following features: a sequential coupling of the OpenFOAMmesh generator (SnappyHexMesh) with the SPH code to define the initial distribution and conditions of the surface wallelements and the fluid semi-particles (DB-SPH); the extension of the DB-SPH boundary treatment to any 3D frontiershape and translational kinematics; a simplified procedure to initialize the integral Shepard coefficient for any 3D initialconditions. Preliminary validations are obtained by comparison with available measures and SPH model intercomparisonsare also reported.

Further, the model couplings have allowed the simulation of a 3D demonstrative dambreach (earth-filled dam break driven by a breach channel) on complex topography, with simultaneous transport of fivetree trunks and granular material. The model couplings have also permitted a preliminary implementation of themolecular viscosity term (in the momentum equation) at boundaries (DB-SPH treatment). This would allowrepresenting the bottom drag at low Reynolds numbers, which are often detected in the bed-load transport layer.

Related Links