Abstract:The extremely strong Wenchuan earthquake triggered thousands of landslides in Sichuan Province, China. Flow-like landslides, such as the Tangjiashan, Wangjiayan, and Donghekou landslides were among the most destructive, causing many casualties and serious economic damage. It is therefore important to identify the flow mechanisms and investigate the specific characteristics of these flow-like landslides. This paper presents a three-dimensional model based on smoothed particle hydrodynamics (SPH) to simulate rapid landslide motion across 3D terrain. The Navier–Stokes equations in a CFD framework are used as governing equations, and artificial viscosity is incorporated into the pressure terms in the momentum equation to dissipate energy to avoid numerical oscillation and particle penetration, thus improving the stability of the numerical results. A non-Newtonian fluid model, the Bingham model, has proven suitable for describing the relationship between the shear strain rate and the shear stress in highly deformed soil materials. In the proposed model it is used as the constitutive model to describe the fluidization characteristics of flow-like landslides combined with the Mohr–Coulomb yield criterion. The model incorporates a no-slip boundary condition, to consider the effect of a solid boundary on slope movement. Ghost particles are created and assigned an artificial velocity. The viscous force caused by the solid boundary is calculated using the relative velocities between the fluid and the boundary particles. Open Multiprocessing (OpenMP), an API for multi-platform shared-memory parallel programming, is used to improve the efficiency of the SPH code running. To show the validity of the proposed approach, a benchmark problem of 3D dam break was simulated. The calculated distances of the surge front at different times agree well with the test results. Numerical modeling of the propagation of the Tangjiashan, Wangjiayan, and Donghekou landslides was performed by the application of SPH models to real flow-like landslides. The whole flow processes of these flow-like landslides across the 3D terrain are represented. The landslides change direction, split or join in, and spread or contract in their flow path in response to the local topography. Time-history curves of the velocity and displacement were obtained to analyze the movement characteristics of the landslide mass. The shapes of the deposition zones after slide occurrence were investigated. Comparisons of the SPH simulated geometry and the surveyed landslide configurations for the Tangjiashan, Wangjiayan and Donghekou landslides were conducted, and show a high degree of similarity. This indicates that the proposed 3D SPH model can accurately represent the evolution of the final slide shape. The prediction of the fluidization characteristics of earthquake-induced flow-like landslides can notably reduce sudden loss of life, as it provides a means for mapping hazardous areas, for estimating the hazard intensity, and for identification and design of appropriate protective measures.

Three-Dimensional Numerical Simulation of Large-Scale LandslideGenerated Surging Waves with a GPU‒Accelerated Soil‒Water Coupled SPH Model

Numerical simulation of the large‐scale Huangtian (China) landslide‐generated impulse waves by a GPU‐accelerated three‐dimensional soil‒water coupled SPH model

Three-dimensional simulations of large-scale long run-out landslides with a GPU-accelerated elasto-plastic SPH model

DEM–SPH Coupling Method for Landslide Surge Based on a GPU Parallel Acceleration Technique

SPH-based Simulation of Flow Process of a Landslide at Hongao Landfill in China

Numerical simulation of submarine landslide tsunamis based on the smoothed particle hydrodynamics model

Numerical Simulation of Landslide-Generated Waves Using a Soil–water Coupling Smoothed Particle Hydrodynamics Model

On three-dimensional SPH modelling of large-scale landslides

Three-Dimensional Numerical Method for Simulating Large-Scale Free Water Surface by Massive Parallel Computing on a GPU

Numerical modeling of landslide-generated impulse waves in mountain reservoirs using a coupled DEM-SPH method

Modeling Near-Field Impulsive Waves Generated by Deformable Landslide Using the HBP Model Based on the SPH Method

SPH-DEM Coupling Method Based on GPU and Its Application to the Landslide Tsunami. Part I: Method and Validation

Numerical Simulation of Impulsive Wave Generated by Landslide on Reservoir Bank Using Two-dimensional Smoothed Particle Hydrodynamics Method

Three-dimensional modeling of dynamic interaction between soils and multibody system with GPU acceleration

Development of a Coupled DDA–SPH Method and Its Application to Dynamic Simulation of Landslides Involving Solid–Fluid Interaction

3D Large-Scale SPH Modeling of Vehicle Wading with GPU Acceleration

3D numerical modeling using smoothed particle hydrodynamics of flow-like landslide propagation triggered by the 2008 Wenchuan earthquake

Analysis of landslide-generated impulsive waves using a coupled DDA-SPH method

GIS Enabled SPH-Soil Modeling for the Post-Failure Flow of Landslides under Seismic Loadings

Study on landslide-induced wave disasters using a 3D coupled SPH-DEM method

3D numerical simulation of debris-flow motion using SPH method incorporating non-Newtonian fluid behavior