Two-Phase Two-Layer Depth-Integrated SPH-FD Model: Application to Lahars and Debris Flows

The complex nature of debris flows suggests that the pore-water pressure evolution and dewatering of a flowing mass caused by the high permeability of soil or terrain could play an essential role in the dynamics behavior of fast landslides. Dewatering causes desaturation, reducing the pore-water pressure and improving the shear strength of liquefied soils. A new approach to landslide propagation modeling considering the dewatering of a mass debris flow has drawn research attention. The problem is characterized by a transition from saturated to unsaturated soil. This paper aims to address this scientific gap. A depth-integrated model was developed to analyze the dewatering of landslides, in which, desaturation plays an important role in the dynamics behavior of the propagation. This study adopted an SPH numerical method to model landslide propagation consisting of pore-water and a soil skeleton in fully or partially saturated soils. In a two-phase model, the soil–water mixture was discretized and represented by two sets of SPH nodes carrying all field variables, such as velocity, displacement, and basal pore-water pressure. The pore-water was described by an additional set of balance equations to take into account its velocity. In the developed two-layer model, an upper desaturated layer and a lower saturated layer were considered to enhance the description of dewatering. This is the so-called two-phase two-layer formulation, which is capable of simulating the entire process of landslides propagation, including the large deformation of soils and corresponding pore-water pressure evolutions, where the effect of the dewatering in saturated soils is also taken into account. A dam-break problem was analyzed through the new and previously developed model. A flume test performed at Trondheim was also used to validate the proposed model by comparing the numerical results with measurements obtained from the experiment. Finally, the model was applied to simulate a real case lahar, which is an appropriate benchmark case used to examine the applicability of the developed model. The simulation results demonstrated that taking into account the effects of dewatering and the vital parameter of relative height is essential for the landslide propagation modeling of a desaturated flowing mass.