A fully resolved SPH-DEM method for heterogeneous suspensions with arbitrary particle shape

Peng, Chong; Zhan, Ling; Wu, Wei; Zhang, Bingyin

doi:10.1016/j.powtec.2021.04.044

Cited by 54 publications

(20 citation statements)

References 111 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…In fact, SPH-solver coupling with other methods has been defined as one of the grand challenges for the SPH community by SPHERIC in 2020 [102]. To date, the SPH-DEM solver for fluid-granule problems (see, e.g., [250][251][252][253][254]) and SPH-FE/TLSPH solvers (see, e.g., [200][201][202][203][204][210][211][212][213][214]) for the fluid-solid problem are relatively mature, whereas for fluid simulations, the "SPH+X" coupling solver is still an open topic that has been attracting more and more attention from the SPH community. As discussed in Section 1, the potential-flow theory and the RANS approach are two popular mesh-based methods for simulating OEDs; it is thus naturally the best method to couple these two schemes with the SPH method.…”

Section: Sph Coupling With Other Efficient Numerical Methodsmentioning

confidence: 99%

A Review of SPH Techniques for Hydrodynamic Simulations of Ocean Energy Devices

et al. 2022

View full text Add to dashboard Cite

This article is dedicated to providing a detailed review concerning the SPH-based hydrodynamic simulations for ocean energy devices (OEDs). Attention is particularly focused on three topics that are tightly related to the concerning field, covering (1) SPH-based numerical fluid tanks, (2) multi-physics SPH techniques towards simulating OEDs, and finally (3) computational efficiency and capacity. In addition, the striking challenges of the SPH method with respect to simulating OEDs are elaborated, and the future prospects of the SPH method for the concerning topics are also provided.

show abstract

Section: Sph Coupling With Other Efficient Numerical Methodsmentioning

confidence: 99%

A Review of SPH Techniques for Hydrodynamic Simulations of Ocean Energy Devices

et al. 2022

View full text Add to dashboard Cite

show abstract

“…First, the particle trajectories are usually simulated by the Discrete Element Method (DEM). Then, the fluid part is solved by some CFD methods using the Finite Volume Method [22], Finite Difference Method [23], and meshfree particle methods [24][25][26] or Lattice-Boltzmann method [27][28][29]. These approaches for fluid-particle two-phase flows are referred to as CFD-DEM, and these methods can be classified into resolved and unresolved CFD-DEM, as shown in Fig.…”

Section: Porous Media and Particulate Flows Simulations With The Darc...mentioning

confidence: 99%

Seepage failure prediction of breakwater using an unresolved ISPH-DEM coupling method incoprporated with Terzaghi's critical hydraulic gradient

Tsuji

Asai²,

Kasama

2022

Preprint

View full text Add to dashboard Cite

This study develops a new numerical simulation model for rubble mound failure prediction caused by piping destruction under seepage flows. The piping has been pointed out as a significant cause of breakwater failure during tsunamis. Once boiling and heaving occur on the mound surface, the piping suddenly propagates in the opposite direction of seepage flow. For the seepage failure prediction, a coupled fluid-soil-structure simulator is developed by combining the ISPH for fluid and the DEM for rubble mounds and caisson blocks. The ISPH, a Lagrangian particle method for incompressible fluids, can simulate seepage and violent flows such as tsunamis. The DEM has been applied for discrete particle and rigid body simulations that include discontinuous deformation as in the rubble mounds failure and large displacement of the caisson block. ISPH-DEM coupling simulations have already been proposed as a technique for multi-phase flows. Still, the technique cannot reproduce the sudden onset of piping from a stable mound. Two simple assumptions are applied to reduce the numerical cost for the fluid-soil-structure simulators of a breakwater structure composed of rubble mound and the caisson block. Firstly, each rubble is modeled as an idealized spherical DEM particle with the mean diameter of rubble. The ISPH particle size is assumed to be the same size as the DEM particle. Under these assumptions, the unresolved coupling model between rubble mound particle and fluid, which obtains the interaction through empirical drag force, should be applied. Whereas the interaction between the fluid and the caisson block is fully resolved with the spatial resolution with the ISPH and DEM particle size. Our new contribution in this paper is how to model the interaction as an unresolved coupling between seepage flow simulated by ISPH and rubble mound particle modeled with DEM. Our original seepage failure experiment is simulated using the proposed ISPH-DEM coupling simulator. We recognized the conventional drag force models as the unresolved coupling model are insufficient to initiate the boiling and piping observed in the experiment. It may be due in one part to excessive averaging of flow velocities caused under unresolved coupling. Therefore, Terzaghi's critical hydraulic gradient is introduced to initiate the boiling and heaving. Unstable DEM particles, judged by Terzaghi's critical hydraulic gradient, gradually lose their mass to represent unresolved suspended fine rubble mound particles. Our models qualitatively reproduce the sand boiling and backward erosion in the opposite direction of the seepage flow, as shown in the experiment.

show abstract

“…In CFD-DEM, the particle trajectories are usually simulated by the Discrete Element Method (DEM). Then, the fluid part is solved by some CFD methods using the Finite Volume Method(FVM) [19], Finite Difference Method(FDM) [20], and meshfree particle methods [21][22][23] or Lattice-Boltzmann method(LBM) [24][25][26]. Whether the fluid around the particles is resolved, CFD-DEM can be classified into two types: resolved coupling model and unresolved coupling model, as shown in Fig.…”

Section: Introductionmentioning

confidence: 99%

Seepage failure prediction of breakwater using an unresolved ISPH-DEM coupling method enriched with Terzaghi’s critical hydraulic gradient

Tsuji

Asai

Kasama

2023

Adv. Model. and Simul. in Eng. Sci.

View full text Add to dashboard Cite

This study develops a new numerical simulation model for rubble mound failure prediction caused by piping destruction under seepage flows. The piping has been pointed out as a significant cause of breakwater failure during tsunamis. Once boiling and heaving occur on the mound surface, the piping suddenly propagates in the opposite direction of seepage flow. For the seepage failure prediction, a coupled fluid-soil-structure simulator is developed by combining the ISPH for fluid and the DEM for rubble mounds and caisson blocks. The ISPH, a Lagrangian particle method for incompressible fluids, can simulate seepage and violent flows such as tsunamis. The DEM has been applied for discrete particle and rigid body simulations that include discontinuous deformation, as in the rubble mounds failure and large displacement of the caisson block. ISPH-DEM coupling simulations have already been proposed as a technique for multi-phase flows. Still, the technique cannot reproduce the sudden onset of piping from a stable mound. Two simple assumptions are applied to reduce the numerical cost for the fluid-soil-structure simulators of a breakwater structure composed of a rubble mound and the caisson block. Firstly, each rubble is modeled as an idealized spherical DEM particle with the mean diameter of the rubble. The ISPH particle size is assumed to be the same size as the DEM particle. Under these assumptions, the unresolved coupling model between rubble mound particles and fluid, which obtains the interaction through empirical drag force, should be applied. At the same time, the interaction between the fluid and the caisson block is fully resolved with the spatial resolution with the ISPH and DEM particle size. Our new contribution in this paper is how to model the interaction as an unresolved coupling between seepage flow simulated by ISPH and rubble mound particle modeled with DEM. Our original seepage failure experiment is simulated using the proposed ISPH-DEM coupling simulator. We identified the conventional drag force models as the unresolved coupling model are insufficient to initiate the boiling and piping observed in the experiment. It may be due in one part to excessive averaging of flow velocities caused by unresolved coupling. Therefore, Terzaghi’s critical hydraulic gradient is introduced to initiate the boiling and heaving. Unstable DEM particles, judged by Terzaghi’s critical hydraulic gradient, gradually lose their mass to represent unresolved suspended fine rubble mound particles. Our models qualitatively reproduce the sand boiling and backward erosion in the opposite direction of the seepage flow, as shown in the experiment.

show abstract

A fully resolved SPH-DEM method for heterogeneous suspensions with arbitrary particle shape

Cited by 54 publications

References 111 publications

A Review of SPH Techniques for Hydrodynamic Simulations of Ocean Energy Devices

A Review of SPH Techniques for Hydrodynamic Simulations of Ocean Energy Devices

Seepage failure prediction of breakwater using an unresolved ISPH-DEM coupling method incoprporated with Terzaghi's critical hydraulic gradient

Seepage failure prediction of breakwater using an unresolved ISPH-DEM coupling method enriched with Terzaghi’s critical hydraulic gradient

Contact Info

Product

Resources

About