Abstract：Solid-liquid two-phase flow generally occurs during natural disasters，such as debris flows and avalanches. To address the shortcomings of current numerical methods，this study developed a two-layer two-phase smoothed particle hydrodynamics(SPH) method by combining soil constitutive equations and considering the complex interaction mechanism between the solid and liquid phases. The simulation of soil-water interaction with large deformation is performed using unified governing equations under the framework of the continuous medium theory，revealing its practical feasibility for researching large-scale debris flows. The accuracy and reliability of the developed numerical method for modeling the liquid and soil motion，as well as the interaction of soil-water are verified through a numerical example of saturated soil collapse in air，which approximates actual debris flows. The sensitivity analysis of fluid viscosity，effective grain size of soil mass，initial void ratio，and other parameters was conducted. Meanwhile，different percolation models for the soil-water coupling are compared. The results suggest that the proposed two-layer two-phase SPH method in this paper is reasonable and effective for simulating the soil-water coupling with large deformation. Notably，the environmental liquid and soil parameters significantly affect the results. When the Reynolds number of the particles is large，percolation models are essential for correctly simulating the interaction between solid and liquid and capturing the interface between the two phases. The investigated results in this paper are of significance for understanding the solid-liquid interaction mechanism in debris flows as well as disaster prevention and mitigation engineering.
[1] DE LANGE S I，SANTA N，PUDASAINI S P，et al. Debris-flow generated tsunamis and their dependence on debris-flow dynamics[J]. Coastal Engineering，2020，157：1–18.
[2] WALDER J S，WATTS P，WAYTHOMAS C F. Case study：mapping tsunami hazards associated with debris flow into a reservoir[J]. Journal of Hydraulic Engineering，2006，132(1)：1–11.
[3] DAWSON A G，STEWART I. Tsunami deposits in the geological record[J]. Sedimentary Geology，2007，200(3/4)：166–183.
[4] 文 强，胡卸文，刘 波，等. 四川丹巴梅龙沟“6•17”泥石流成灾机理分析[J]. 中国地质灾害与防治学报，2022，33(3)：23–30. (WEN Qiang，HU Xiewen，LIU Bo，et al. Analysis on the mechanism of debris flow in Meilong valley in Danba County on June 17，2020[J]. The Chinese Journal of Geological Hazard and Control，2022，33(3)：23–30.(in Chinese))
[5] 黄 琦，谢 涛，陈欢欢，等. 2020年8月甘肃武都群发泥石流基本特征及减灾建议[J]. 公路，2022，67(9)：30–37.(HUANG Qi，XIE Tao，CHEN Huanhuan，et al. Basic Characteristics and disaster reduction suggestions of cluster debris flow in Wudu，Gansu in August 2020[J]. Highway，2022，67(9)：30–37.(in Chinese))
[6] MENG X N，JOHNSON C G，GRAY J M N T. Formation of dry granular fronts and watery tails in debris flows[J]. Journal of Fluid Mechanics，2022，943(A19)：1–43.
[7] TAYLOR-NOONAN A M，BOWMAN E，MCARDELL B，et al. Influence of pore fluid on grain-scale interactions and mobility of granular flows of differing volume[J]. Earth and Space Science Open Archive，2022，127(12)：1–57.
[8] TAYLOR-NOONAN A M，BOWMAN E，MCARDELL B，et al. Effect of pore fluid on grain-scale interactions and mobility of granular flows：Comparison of distal reach for increasing source volume for dry and saturated granular flows[J]. Earth and Space Science Open Archive，2022，127(2)：1–54.
[9] TRUJILLO-VELA M G，RAMOS-CAÑÓN A M，ESCOBAR-VARGAS J A，et al. An overview of debris-flow mathematical modelling[J]. Earth-Science Reviews，2022，232：104135.
[10] HIRT C W，NICHOLS B D. Volume of fluid (VOF) method for the dynamics of free boundaries[J]. Journal of Computational Physics，1981，39(1)：201–225.
[11] TOMÉ M F，GROSSI L，CASTELO A，et al. A numerical method for solving three-dimensional generalized Newtonian free surface flows[J]. Journal of Non-Newtonian Fluid Mechanics，2004，123(2)：85–103.
[12] OSHER S，SETHIAN J A. Fronts propagating with curvature-dependent speed：Algorithms based on Hamilton-Jacobi formulations[J]. Journal of Computational Physics，1988，79(1)：12–49.
[13] GINGOLD R A，MONAGHAN J J. Smoothed particle hydrodynamics：theory and application to non-spherical stars[J]. Monthly Notices of the Royal Astronomical Society，1977，181(2)：375–389.
[14] MONAGHAN J J. Smoothed particle hydrodynamics and its diverse applications[J]. Annual Review of Fluid Mechanics，2012，44(1)：323–346.
[15] LIU C，YU Z X，ZHAO S C. A coupled SPH-DEM-FEM model for fluid-particle-structure interaction and a case study of Wenjia gully debris flow impact estimation[J]. Landslides，2021，18(7)：2 403– 2 425.
[16] SUN X S，SAKAI M，YAMADA Y. Three-dimensional simulation of a solid-liquid flow by the DEM-SPH method[J]. Journal of Computational Physics，2013，248(1)：147–176.
[17] MARKAUSKAS D，KRUGGEL-EMDEN H，SCHERER V. Numerical analysis of wet plastic particle separation using a coupled DEM-SPH method[J]. Powder Technology，2018，325：218–227.
[18] WU K，YANG D，WRIGHT N，et al. An integrated particle model for fluid-particle-structure interaction problems with free-surface flow and structural failure[J]. Journal of Fluids and Structures，2018，76：166–184.
[19] HE K，SHI H B，YU X P. Effects of interstitial water on collapses of partially immersed granular columns[J]. Physics of Fluids，2022，34：23306.
[20] HE Y，BAYLY A E，HASSANPOUR A，et al. A GPU-based coupled SPH-DEM method for particle-fluid flow with free surfaces[J]. Powder Technology，2018，338：548–562.
[21] ZENG Q Y，ZHENG M X，H D. Numerical simulation of impact and entrainment behaviors of debris flow by using SPH-DEM-FEM coupling method[J]. Open Geosciences，2022，14(1)：1 020–1 047.
[22] LI J J，QIU L C，TIAN L，et al. Modeling 3D non-Newtonian solid-liquid flows with a free-surface using DEM-MPS[J]. Engineering Analysis with Boundary Elements，2019，105：70–77.
[23] XIE F Z，MENG Q J，WAN D C. Numerical simulations of liquid-solid flows in a vertical pipe by MPS-DEM coupling method[J]. China Ocean Engineering，2022，36(4)：542–552.
[24] XIE F Z，ZHAO W W，WAN D C. Numerical simulations of liquid-solid flows with free surface by coupling IMPS and DEM[J]. Applied Ocean Research，2021，114：102771.
[25] CECCATO F，LEONARDI A，GIRARDI V，et al. Numerical and experimental investigation of saturated granular column collapse in air[J]. Soils and Foundations，2020，60(3)：683–696.
[26] BATCHELOR G K. An introduction to fluid dynamics[M]. Cambridge：Cambridge University Press，2000：53–69.
[27] 揭天琪. 碎屑流运动过程固–流转化机制数值实验研究[博士学位论文][D]. 重庆：重庆大学，2018.(JIE Tianqi，Numerical experiment research of solid-fluid transform criterion of debris flow sliding process[Ph. D. Thesis][D]. Chongqing：Chongqing University，2018.(in Chinese))
[28] BUI H H，FUKAGAWA R，SAKO K，et al. Lagrangian meshfree particles method(SPH) for large deformation and failure flows of geomaterial using elastic-plastic soil constitutive model[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2008，32(12)：1 537–1 570.
[29] CHEN W F，MIZUNO E. Nonlinear analysis in soil mechanics：theory and implementation[M]. Amsterdam：Elsevier，1990：398–409.
[30] LAMBE T W，WHITMAN R V. Soil mechanics[M]. New York：John Wiley and Sons，1969：241–294.
[31] GHOLAMI K M，GALINDO-TORRES S A，SCHEUERMANN A，et al. Smoothed Particle Hydrodynamics for investigating hydraulic and mechanical behaviour of an embankment under action of flooding and overburden loads[J]. Computers and Geotechnics，2018，94：31–45.
[32] GENT M R. Stationary and oscillatory flow through coarse porous media[R]. Delft：Delft University of Technology，1993.
[33] MONAGHAN J J. Simulating Free surface flows with SPH[J]. Journal of Computational Physics，1994，110(2)：399–406.
[34] DALRYMPLE R A，KNIO O. SPH modeling of water waves[C]// Coastal Dynamics 2001. Lund：American Society of Civil Engineers，2001：779–787.
[35] MORRIS J P，FOX P J，ZHU Y. Modeling low reynolds number incompressible flows using SPH[J]. Journal of Computational Physics，1997，136(1)：214–226.
2023, Vol. 42 
