Abstract:Abstract Climate warming accelerates permafrost thawing, causing warming‐driven disasters like ground collapse and retrogressive thaw slump (RTS). These phenomena, involving intricate multiphysics interactions, phase transitions, nonlinear mechanical responses, and fluid‐like deformations, and pose increasing risks to geo‐infrastructures in cold regions. This study develops a thermo‐hydro‐mechanical (THM) coupled single‐point three‐phase material point method (MPM) to simulate the time‐dependent phase transition and large deformation behavior arising from the thawing or freezing of ice/water in porous media. The mathematical framework is established based on the multiphase mixture theory in which the ice phase is treated as a solid constituent playing the role of skeleton together with soil grains. The additional strength due to ice cementation is characterized via an ice saturation‐dependent Mohr–Coulomb model. The coupled formulations are solved using a fractional‐step‐based semi‐implicit integration algorithm, which can offer both satisfactory numerical stability and computational efficiency when dealing with nearly incompressible fluids and extremely low permeability conditions in frozen porous media. Two hydro‐thermal coupling cases, that is, frozen inclusion thaw and Talik closure/opening, are first benchmarked to show the method can correctly simulate both conduction‐ and convection‐dominated thermal regimes in frozen porous systems. The fully THM responses are further validated by simulating a 1D thaw consolidation and a 2D rock freezing example. Good agreements with experimental results are achieved, and the impact of hydro‐thermal variations on the mechanical responses, including thaw settlement and frost heave, are successfully captured. Finally, the predictive capability of the multiphysics MPM framework in simulating thawing‐triggered large deformation and failure is demonstrated by modeling an RTS and the settlement of a strip footing on thawing ground.

A TOUGH-FEMM based cryogenic THM coupled model and its application to cold-region tunnels

THM model of rock tunnels in cold regions and numerical simulation

A fully coupled thermo-hydro-mechanical model for fractured rock masses in cold regions

A coupled cryogenic thermo-hydro-mechanical model for frozen medium: Theory and implementation in FDEM

Evolution Patterns of Frost-Heaving Pressure with Partial Ponding in Cold Region Tunnels

Thermal-hydraulic-mechanical coupling behavior and frost heave mitigation in freezing soil

An Analytical Solution of Frost Heaving Pressure for Cold-region Tunnel Considering Freeze-thaw Cycles and in-situ Stress

A Frost Heaving Prediction Approach for Ground Uplift Simulation Due to Freeze-Sealing Pipe Roof Method

Model test on frost heaving pressure induced by frozen partial ponding behind tunnel lining in cold region

A Fully Thermo-Hydro-Mechanical Coupling Simulation at Low Temperature in Underground LNG Storage

Numerical simulation of frost heave of saturated soil considering thermo‐hydro‐mechanical coupling

Meso-scale mechanical model for mortar deformation under freeze thaw cycles

Mesoscale Simulation of Concrete Behavior with Non-Uniform Frost Damage with Verification by CT Imaging

A Thermo-Poromechanical Model for Simulating Freeze–thaw Actions in Unsaturated Soils

An Analytical Solution for the Frost Heaving Force considering the Freeze-Thaw Damage and Transversely Isotropic Characteristics of the Surrounding Rock in Cold-Region Tunnels

A New Analytical Solution on the Frost Heaving Force of Circular Tunnel in Cold Regions

Experimental analysis of the thermo-hydro-mechanical (THM) coupling in freezing vertical shafts of unsaturated sandy soil

Thermo‐hydro‐mechanical coupled material point method for modeling freezing and thawing of porous media

Mechanical analysis of frost heaving of weakly cemented soft rock tunnels in cold regions

Experimental study and numerical modeling of the thermo-hydro-mechanical processes in soil freezing with different frost penetration directions

Preliminary identification of the failure mode of shallow tunnels in soil subjected to frost heave: Model test and numerical simulation