SPH modeling of water-related natural hazards

This paper collects some recent smoothed particle hydrodynamic (SPH) applications in the field of natural hazards connected to rapidly varied flows of both water and dense granular mixtures including sediment erosion and bed load transport. The paper gathers together and outlines the basic aspects of some relevant works dealing with flooding on complex topography, sediment scouring, fast landslide dynamics, and induced surge wave. Additionally, the preliminary results of a new study regarding the post-failure dynamics of rainfall-induced shallow landslide are presented. The paper also shows the latest advances in the use of high performance computing (HPC) techniques to accelerate computational fluid dynamic (CFD) codes through the efficient use of current computational resources. This aspect is extremely important when simulating complex three-dimensional problems that require a high computational cost and are generally involved in the modeling of water-related natural hazards of practical interest. The paper provides an overview of some widespread SPH free open source software (FOSS) codes applied to multiphase problems of theoretical and practical interest in the field of hydraulic engineering. The paper aims to provide insight into the SPH modeling of some relevant physical aspects involved in water-related natural hazards (e.g., sediment erosion and non-Newtonian rheology). The future perspectives of SPH in this application field are finally pointed out

SPH Simulation of an Air-Assisted Atomizer Operating at High Pressure : Influence of Non-Newtonian Effects

A twin-fluid atomizer configuration is predicted by means of the 2D weakly-compressible Smooth Particle Hydrodynamics (SPH) method and compared to experiments. The setup consists of an axial liquid jet fragmented by a co-flowing high-speed air stream (Ug ~ 60 m/s) in a pressurized atmosphere up to 11 bar (abs.). Two types of liquid are investigated: a viscous Newtonian liquid (µ = 200 mPa.s) obtained with a glycerol/water mixture and a viscous non-Newtonian liquid (µ ~ 150 mPa.s) obtained with a carboxymethyl cellulose (CMC) solution. 3D effects are taken into account in the 2D code by introducing (i) a surface tension term, (ii) a cylindrical viscosity operator and (iii) a modified velocity accounting for the divergence of the volume in the radial direction. The numerical results at high pressure show a good qualitative agreement with experiment, i.e. a correct transition of the atomization regimes with regard to the pressure, and similar dynamics and length scales of the generated ligaments. The predicted frequency of the Kelvin-Helmholtz instability needs a correction factor of 2 to be globally well recovered with the Newtonian liquid. The simulation of the non-Newtonian liquid at high pressure shows a similar breakup regime with finer droplets compared to Newtonian liquids while the simulation at atmospheric pressure shows an apparent viscosity similar to the experiment

Bingham fluid simulations using a physically consistent particle method

The Bingham fluid simulation model was constructed and validated using a physically consistent particle method, i.e., the Moving Particle Hydrodynamics (MPH) method. When a discrete particle system satisfies the fundamental laws of physics, the method is asserted as physically consistent. Since Bingham fluids sometimes show solid-like behaviors, linear and angular momentum conservation is especially important. These features are naturally satisfied in the MPH method. To model the Bingham feature, the viscosity of the fluid was varied to express the stress-strain rate relation. Since the solid-like part, where the stress does not exceed the yield stress, was modeled with very large viscosity, the implicit velocity calculation was introduced so as to avoid the restriction of the time step width with respect to the diffusion number. As a result, the present model could express the stopping and solid-like behaviors, which are characteristics of Bingham fluids. The proposed method was verified and validated, and its capability was demonstrated through calculations of the two-dimensional Poiseuille flow of a Bingham plastic fluid and the three-dimensional dam-break flow of a Bingham pseudoplastic fluid by comparing those computed results to theory and experiment

DualSPHysics: from fluid dynamics to multiphysics problems

DualSPHysics is a weakly compressible smoothed particle hydrodynamics (SPH) Navier–Stokes solver initially conceived to deal with coastal engineering problems, especially those related to wave impact with coastal structures. Since the first release back in 2011, DualSPHysics has shown to be robust and accurate for simulating extreme wave events along with a continuous improvement in efficiency thanks to the exploitation of hardware such as graphics processing units for scientific computing or the coupling with wave propagating models such as SWASH and OceanWave3D. Numerous additional functionalities have also been included in the DualSPHysics package over the last few years which allow the simulation of fluid-driven objects. The use of the discrete element method has allowed the solver to simulate the interaction among different bodies (sliding rocks, for example), which provides a unique tool to analyse debris flows. In addition, the recent coupling with other solvers like Project Chrono or MoorDyn has been a milestone in the development of the solver. Project Chrono allows the simulation of articulated structures with joints, hinges, sliders and springs and MoorDyn allows simulating moored structures. Both functionalities make DualSPHysics especially suited for the simulation of offshore energy harvesting devices. Lately, the present state of maturity of the solver goes beyond single-phase simulations, allowing multi-phase simulations with gas–liquid and a combination of Newtonian and non-Newtonian models expanding further the capabilities and range of applications for the DualSPHysics solver. These advances and functionalities make DualSPHysics an advanced meshless solver with emphasis on free-surface flow modelling

Numerical simulation of single droplet dynamics in three-phase flows using ISPH

In this study, a new surface tension formulation for modeling incompressible, immiscible three-phase fluid flows in the context of incompressible smoothed particle hydrodynamics (ISPH) in two dimensions has been proposed. A continuum surface force model is employed to transform local surface tension force to a volumetric force while physical surface tension coefficients are expressed as the sum of phase specific surface tension coefficients, facilitating the implementation of the proposed method at triple junctions where all three phases are present. Smoothed color functions at fluid interfaces along with artificial particle displacement throughout the computational domain are combined to increase accuracy and robustness of the model. In order to illustrate the effectiveness of the proposed method, several numerical simulations have been carried out and results are compared to analytical data available in literature. Results obtained by simulations are compatible with analytical data, demonstrating that the ISPH scheme proposed here is capable of handling three-phase flows accurately