Extensions of the two-phase double-point material point method to simulate the landslide-induced surge process

https://doi.org/10.1016/j.enganabound.2021.09.020Get rights and content

Abstract

Landslide-induced surges, as a complex type of fluid-solid coupling problem, are widespread in mountain regions and cause disastrous consequences. In order to reproduce the entire process of landslide-induced surge, by taking full advantage of material point method (MPM) in simulating large deformation of soil and gravity-driven flow of water, the two-phase double-point material point method (TPDP-MPM) was extended: a new algorithm handling the boundary without setting boundary particles was proposed. The classical submarine block-induced surge test and sand column collapse experiment were numerically simulated by the MPM and verified the effectiveness of the program. After that, the multiphase coupling process of Lituya Bay landslide-induced surge was simulated and verified the reliability of the TPDP-MPM code, which have achieved parallel computing by making use of OpenMP model. By comparing the results of numerical simulation with the model experiment, the two-phase double-point material point method can well simulate the starting, propagation, run-up and reflux stage of landslide-induced surge. Moreover, the extension mechanism of landslide-surge disaster chain in both time and space is revealed from the perspective of energy. This study provides a reliable tool for the analysis and assessment of landslide-induced surge disasters.

Introduction

Obtaining comprehensive hazard mechanisms, such as avalanches, landslides, and landslide-induced dammed lakes, has been and is still a major challenge for hazard assessment [13], [42], [46]. Rock avalanches and landslides can reach extremely high velocities along subaerial slopes; as a consequence, the remarkable kinetic energy of the sliding masses can be the origination of large waves or surges in water bodies. The landslide-induced surge process is an essential link of disaster chains in near-water regions, which establishes a connection between the landslide mass and water bodies, such as rivers, lakes, fjords and reservoirs, and directly affects the safety of infrastructure and residents both upstream and downstream. Therefore, the analysis and assessment of the geological disaster process is the top priority in the prevention and control work of geological disaster chain.

With the continuous improvement of computer computing efficiency, numerical calculation has become an important mean of geological disaster analysis and assessment. Large deformation continuous numerical methods such as Smoothed Particle Hydrodynamics (SPH) [2], [19], [36], Material Point Method (MPM) [3], [6], [8], [9], [12], [21], [29], [33], [34], [40], [41], [43] and discontinuous numerical methods such as DDA [15], [16] , NMM [37], [38], [39] and DEM [23] are both widely used in this field. With regard to the large deformation of solid and flow of fluid problems in the landslide-surge disaster chain, the material point method is adopted in this paper since the meshless characteristic and the higher computational efficiency.

In addition to large deformation of landslides and flow of water, the landslide-surge disaster chain also involves a complex solid-liquid coupling process, which is difficult to simulate precisely. More recently, efforts have concentrated on simulating coupled hydromechanical processes in saturated or unsaturated soils with MPM. The interaction between solid and fluid phases is formulated in two different approaches in the material point method: adopting either one set of MPs (i.e., single point manner) [20], [43] or two separate sets of MPs (i.e., double point manner) [7], [22], [39], [45]. For the landslide-induced impulsive surge process, it should be mentioned that the water is not only liquid within the pores but also a free liquid in this issue. Therefore, we extended the two-phase double point material point method (TPDP-MPM) program to simulate and investigate this typical fluid-solid coupling process.

This paper is organized as follows: First, the theory of the two-phase double-point approach of the MPM is introduced. Subsequently, by comparing two numerical examples with the benchmark model test results, the material point method is shown to accurately reproduce the gravity-driven flow of water and the large deformations of the soil. Finally, the Lituya Bay landslide-induced surge process, including the impact stage of the landslide mass, and the propagation and run-up stage of impulsive waves [10], is reproduced by using the two-phase double-point MPM.

Section snippets

Efficient two-phase double-point material point method

In this section we will give a brief introduction to the two-phase double-point material point method; for more details, refer to [31]. The governing equation for describing the motion of saturated soils is the v–w formulation, where velocities or displacements of the solid and fluid phases are considered as the primary variables. The v–w formulation was extended to the MPM by Jassim et al., [20] where the velocities of each phase become fundamental unknowns of the coupled system.

Benchmark example 1

An experimental study on impulsive waves generated by a submarine block was conducted by P. Heinrich [17]. The submarine triangular block was allowed to slide freely along a 45° slope, and a buffer was set to stop the block on the bottom. The initial state of this model and other information are described in Fig. 4(a).

The surface of the submarine block, the two inclined slopes, and the bottom were frictionless during the experiment. The elevation data of the free surface at 0.5 s, 1.0 s, and

Lituya bay landslide-induced surge analysis

On July 8, 1958, the Lituya Bay landslide was induced by an 8.3 Richter magnitude earthquake, resulting in an estimated 3 × 107 m3 of sliding mass in the Gilbert Inlet. Fig. 9 shows an aerial photograph of Lituya Bay with a graphical representation of the landslide area and the wave run-up area. There is a maximum bay depth of 122 m and a maximum run-up of 524 m above sea level on the opposite shore, which is partially overtopped. According to site conditions, the surge height of Lituya Bay is

Conclusion

In this study, the two-phase double-point material point method was used to reproduce this typical fluid-solid coupling process. By using OpenMP, the program was parallelized to analysis the dynamic process of surge initiation, wave propagation and the strong dynamic coupling between solids and fluids. The main conclusions resulting from this study are described as follows:

  • (1)

    The material point method was calibrated by reconstructing and simulating two model experiments regarding to submarine

Declaration of Competing Interest

The authors declare no conflict of interest.

Acknowledgments

The work reported in this paper is financially supported by the Youth Innovation Promotion Association CAS (no. 2021325), the International Partnership Program of Chinese Academy of Sciences Grant no. 131551KYSB20180042, and the National Natural Science Foundation of China (no. 51779250; no. 52179117).

A special acknowledgement should be expressed to China-Pakistan Joint Research Center on Earth Sciences that supported the implementation of this study.

References (46)

  • Y. Yang et al.

    Stability analysis of soil-rock-mixture slopes using the numerical manifold method[J]

    Eng Anal Bound Elem

    (2019)
  • A. Yerro et al.

    Run-out of landslides in brittle soils

    Comput Geotech

    (2016)
  • J. Yu et al.

    Experimental study and discrete element method modeling of compression and permeability behaviors of weakly anisotropic sandstones[J]

    Int J Rock Mech Min Sci

    (2020)
  • X.P. Zhou et al.

    Multiscale numerical modeling of propagation and coalescence of multiple cracks in rock masses

    Int J Rock Mech Min Sci

    (2012)
  • K. Abe et al.

    Material point method for coupled hydromechanical problems

    J Geotech Geoenviron Eng

    (2014)
  • A.M. Xenakis et al.

    Landslides and tsunamis predicted by incompressible smoothed particle hydrodynamics (SPH) with application to the 1958 Lituya Bay event and idealized experiment.Proc

    R. Soc. A.

    (2017)
  • E.E. Alonso et al.

    Recent developments of the material point method for the simualtion of landslides

  • M. Becker et al.
  • F. Ceccato et al.

    Developing and testing multiphase MPM approaches for the stability of dams and river embankments

    AIMETA

    (2019)
  • G.B. Crosta et al.

    Interaction of landslide mass and water resulting in impulse waves

    (2013)
  • J. Fern et al.

    Experimental and numerical investigations of dyke failures involving soft materials

    Eng Geol

    (2017)
  • FrancoA. et al.

    The 1958 lituya bay tsunami -pre-event bathymetry reconstruction and 3d numerical modelling utilising the computational fluid dynamics software flow-3d

    Nat. Hazards. Earth. Syst. Sci.

    (2020)
  • H.M. Fritz et al.

    Lituya Bay case: rockslide impact and wave run-up

    Sci Tsunami Hazard

    (2001)
  • Cited by (15)

    • Two-phase two-layer SNS-PFEM for hydromechanical geotechnical large deformation problems

      2024, Computer Methods in Applied Mechanics and Engineering
    • Multiphysics coupling material point method for modelling frost heave of rock slope

      2022, Computers and Geotechnics
      Citation Excerpt :

      Thus, this method is not applicable for phase transitions, soil liquefaction and consolidation, erosion and deposition, etc. Subsequently, the two-phase double-point approach was proposed (Abe et al., 2014; Bandara and Soga, 2015; Bolognin et al., 2017; Martinelli et al., 2017; Zhao et al., 2017; Du et al., 2021), which uses two sets of material points to describe the liquid and solid phases, respectively. It assumes that each area in space is occupied by particles of both phases, each representing a component in the mixture, and that the two sets of particles follow their respective laws of motion within the constraints of interactions.

    View all citing articles on Scopus
    View full text