A mesh-less smoothed particle hydrodynamics (SPH) model for bed-load transport on erosional dam-break floods is presented. This mixture model describes both the liquid phase and the solid granular material. The model is validated on the results from several experiments on erosional dam breaks. A comparison between the present model and a 2-phase SPH model for geotechnical applications (Gadget Soil; TUHH) is performed. A demonstrative 3D erosional dam break on complex topography is investigated. The present 3D mixture model is characterised by: no tuning parameter for the mixture viscosity; consistency with the Kinetic Theory of Granular Flow; ability to reproduce the evolution of the free surface and the bed-load transport layer; applicability to practical problems in civil engineering. The numerical developments of this study are represented by a new SPH scheme for bed-load transport, which is implemented in the SPH code SPHERA v.8.0 (RSE SpA), distributed as FOSS on GitHub.

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.

This work illustrated an application of the FOSS code SPHERA v.8.0 (RSE SpA, Milano, Italy) to the simulation of landslide hazard at the slope of a water basin. SPHERA is based on the weakly compressible SPH method (WCSPH) and holds a mixture model, consistent with the packing limit of the Kinetic Theory of Granular Flow (KTGF), which was previously tested for simulating two-phase free-surface rapid flows involving water-sediment interaction. In this study a limiting viscosity parameter was implemented in the previous formulation of the mixture model to limit the growth of the apparent viscosity, thus saving computational time while preserving the solution accuracy. This approach is consistent with the experimental behavior of high polymer solutions for which an almost constant value of viscosity may be approached at very low deformation rates near the transition zone of elastic–plastic regime. In this application, the limiting viscosity was used as a numerical parameter for optimization of the computation. Some preliminary tests were performed by simulating a 2D erosional dam break, proving that a proper selection of the limiting viscosity leads to a considerable drop of the computational time without altering significantly the numerical solution. SPHERA was then validated by simulating a 2D scale experiment reproducing the early phase of the Vajont landslide when a tsunami wave was generated that climbed the opposite mountain side with a maximum run-up of about 270 m. The obtained maximum run-up was very close to the experimental result. Influence of saturation of the landslide material below the still water level was also accounted, showing that the landslide dynamics can be better represented and the wave run-up can be properly estimated.