This computational study focuses on the thermo-hydro-mechanical simulations of the behaviors of freezing soils used for artificial ground freezing (AGF) in a metro project. Leveraging the experimental and field data available in the literature, we simulate the sequential freezing and excavation of a twin tunneling that occurred in months during the actual construction of the tunnel. A thermo-hydro-mechanical model is developed to capture the multi-physical rate-dependent behaviors triggered by phase transitions, as well as the creeping and secondary consolidation of the soil skeleton and the ice crystals. We then calibrate the material models and establish the THM finite element model coupled with the rate-dependent multi-physical models, which may accurately predict the surface heave induced by ground freezing throughout the project. To showcase the potential of using simulations to guide the AGF, we simulate the scenario where a simultaneous freezing scheme is employed as an alternative to the actual sequential scheme design. We then compared the simulated performance with the recorded results obtained from the sequential scheme. Finally, parametric studies on the effect of ground temperature, the porosity of the frozen soil, and the intrinsic elastic modulus of the solid skeleton are conducted. The maximum surface heave is inferred from finite element simulations to quantify the sensitivity and the impact on the safety of AGF operations.

Accurately describing the solid-like and fluid-like behaviors of granular media is crucial in geotechnical engineering. While the unified frictional-collisional model, integrating rate-independent frictional and ratedependent collisional stresses, is widely used for solid-fluid phase transitions, an effective model is still under investigation, and comprehensive analyses are lacking. This study addresses these gaps by developing an enhanced elastoplasticity-based frictional-collisional model. The frictional stress is modeled using a critical-statebased elastoplasticity approach, and the collisional stress is formulated through an enhanced kinetic theory incorporating particle stiffness. Subsequently, comprehensive element simulations are conducted to explore the effects of concentration, particle stiffness, and strain rate paths on the model. The proposed model's effectiveness is also validated against experimental data. Finally, a detailed comparison with the typical mu(I) rheology model and a state-equation-based phase transition model is conducted. Our analyses show that the developed model effectively captures strain rate path and particle stiffness through the collisional stress component, while concentration-dependent characteristics are captured through both frictional and collisional stress components. Through comparative analyses, we also found that both the state-equation-based and elastoplasticity-based models depict solid-like behavior and replicate the rheology of granular media in a fluid-like state, similar to the mu(I) model. However, they differ in implementing critical state theory: the state-equation-based model acts as a partial-range phase transition model, describing stress evolution from the critical state to the fluid-like state, while the proposed elastoplasticity-based model serves as a full-range phase transition model, covering stress evolution from the initial to the fluid-like state.

Accurate continuum modelling of granular flows is essential for predicting geohazards such as flow-like landslides and debris flows. Achieving such precision necessitates both a robust constitutive model for granular media and a numerical solver capable of handling large deformations. In this work, a novel unified phase transition constitutive model for granular media is proposed that follows a generalized Maxwell framework. The stress is divided into an elastoplastic part and a viscous part. The former utilizes a critical-state-based elastoplasticity model, while the latter employs a strain acceleration-based mu(I) rheology model. Key characteristics such as nonlinear elasticity, nonlinear plastic hardening, stress dilatancy, and critical state concept are incorporated into the elastoplasticity model, and the non-Newtonian mu(I) rheology model considers strain rate and strain acceleration (i.e., a higher-order derivative of strain) to capture changes in accelerated and decelerated flow conditions. A series of element tests is simulated using the proposed unified phase transition model, demonstrating that the novel theory effectively describes the transition of granular media from solid-like to fluid-like states in a unified manner. The proposed unified model is then implemented within the material point method (MPM) framework to simulate 2D and 3D granular flows. The results show remarkable consistency with results from experiments and other numerical methods, demonstrating the model's accuracy in capturing solid-like behaviour during inception and deposition, as well as liquid-like behaviour during propagation.

The enzyme-induced calcium carbonate precipitation (EICP) method has been utilized for curing low-permeability clay by directly mixing the reaction solution with soil. The added reaction solution quantity is limited by the optimal water content, producing insufficient calcium carbonate. Herein, the high-activity urease and high-concentration cementation solution efficacy in treating dispersive soils was evaluated. Phase transitions and structural modifications in EICP-cured soils were investigated through oscillatory amplitude scanning. The soil gradation influence on the EICP treatment effectiveness was assessed. The fluidized EICP-cured soil cementation and rupture mechanisms were investigated by viscosity measurements, electron microscopy, and zeta potential evaluations. A 3 M cementation solution, coupled with 500g/L of soybean urease, significantly enhanced the soil shear resistance, increasing it by 339% to 1807%. The EICP-cured soil gradually transitioned from a fluid to a paste and eventually to a solid within 168 h. High-clay-particle-content soils exhibited pronounced increases in shear resistance after EICP treatment. Under dynamic loading, three shear crack types emerged in EICP-cured soils, emphasizing the importance of soybean protein viscosity and calcium carbonate crystal filling-bonding capability in enhancing soil structural stability. The fluid solidification effectiveness in treating fine-grained soils utilizing EICP was validated through erosion trenches in fluid-solidified check dams, validating its potential.

Under cyclic loading, sand will undergo a solid-liquid phase transition during the liquefaction. This study utilizes discrete element method (DEM) to investigate the stage characteristics of sand macroscopic stress-strain response during the solid-liquid phase transition. The microscopic mechanism of sand solid-liquid phase transition is elucidated from the perspective of contact network. The results indicate that based on the sand flowability, the liquefaction process can be divided into solid phase, solid-liquid transition phase, and liquid phase stages. The strong contact network within the sand is the primary contributor to its effective stress, and the degradation of the originally well-connected strong contact network are the reasons for the sand solid-liquid phase transition. A parameter xi c has been proposed to measure the connectivity of the strong contact network. The weak contacts between particles dominates the sliding and rolling between particles, which is the reason for the macroscopic deformation and flow of sand.

The frost damage of rock mass poses a serious threat to the safety and stability of tunnels in cold regions, and the related thermo-hydro-mechanical (THM) coupling model under low-temperature conditions has been a key focus of research. This paper proposed a cryogenic THM coupled model (TOUGH-FEMM) to study the frost heave behavior of cold-region tunnels. Key issues including heat transfer, thermal stress, water-ice phase transition, unfrozen water, frost heave deformation, and ice-rock interaction are systematically addressed in the proposed model. Specifically, frost pressure in pores and cracks is derived separately to better simulate the ice expansion effect in rock masses. The proposed model is first validated against an experimental test and then applied to a practical cold-region tunnel to reveal the evolution of temperature, frost pressure and frost heave fields, as well as the tunnel stability. Moreover, the effects of cracks and frost damage on tunnel stability under freeze-thaw cycles are discussed. The work detailed herein provides an efficient tool for the THM coupled process in cold- region tunnels.

This manuscript investigated the thermal stability, crystal reconstruction and microstructure evolution of graphite tailing cement mortar subjected to high temperature. Simultaneously, a computational model for heat transfer and degradation considering chemical transformations had been developed by combining multiscale mechanics with the laws of thermodynamics. The results show that 20% graphite tailings can increase the tobermorite crystal content under high temperature and inhibit its transformation into disordered form. Furthermore, the dormant active SiOx in graphite tailing is gradually activated under the action of high temperature, which catalyzes and induces the formation of more belite crystal in graphite tailing cement mortar. Additionally, 40% graphite tailing can promote the generation of anorthite by the induction of high temperature. Finally, a new multi-scale model considering the chemical transformation is established to calculate the hightemperature degradation process of graphite tailing cement mortar.

Integrating liquid CO2 phase transition blasting (LCPTB) technology with hydraulic fracturing (HF) methods can help reduce wellbore damage, create multiple radial fractures, and establish a complex fracture network. This approach significantly increases the recovery efficiency of low-permeability oil and gas fields. Accurately calculating the number of fractures caused by LCPTB is necessary to predict production enhancement effects and optimize subsequent HF designs. However, few studies are reported on large-scale physical model experiments in terms of a method for calculating the fracture number. This study analyzed the initiation and propagation of cracks under LCPTB, derived a calculation formula for crack propagation radius under stress waves, and then proposed a new, fast, and accurate method for calculating the fracture number using the principle of mass conservation. Through ten rock-breaking tests using LCPTB, the study confirmed the effectiveness of the proposed calculation approach and elucidated the variation rule of explosion pressure, rock-breaking scenario, and the impact of varying parameters on fracture number. The results show that the new calculation method is suitable for fracturing technologies with high pressure rates. Recommendations include enlarging the diameter of the fracturing tube and increasing the liquid CO2 mass in the tube to enhance fracture effectiveness. Moreover, the method can be applied to other fracturing technologies, such as explosive fracturing (EF) within HF formations, indicating its broader applicability and potential impact on optimizing unconventional resource extraction technologies. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Water-filled capillary tubes are a kind of standard component in both life science (e.g., blood vessels, interstitial pores, and plant vessels) and engineering (e.g., MEMS microchannel resonators, heat pipe wicks, and watersaturated soils). Under sufficiently low temperatures, water in capillary tubes undergoes phase transition and exhibits frost heave, which can cause deformation, damage, and even fracture of tube wall. However, the thermomechanical analysis of freezing water-filled capillary tubes remains obscure, particularly regarding the rapid change in water temperature due to thermal transient effects. We develop a thermal model of freezing in a waterfilled capillary tube that is suddenly exposed to cold air flow, with the time domain divided into two regimes, separated by the thermal penetration time tp. The effect of thermal penetration on temperature distribution is solved. Then, a distinction is made between freezing occurring before thermal penetration and those occurring after thermal penetration. We next analyze transient mechanical stresses acting at tube wall, with interfacial tension and frost heave effect accounted for. Results obtained are not only useful for preventing frost heave failure but also provide theoretical guidance for tailoring the freezing resistance of microfluidic devices used in MEMS.

Water migration behavior is the main cause of engineering disasters in cold regions, making it essential to understand its mechanisms and the resulting mechanical characteristics for engineering protection. This study examined the water migration process during soil freezing through both experimental and numerical simulations, focusing on the key mechanical outcomes such as deformation and pore water pressure. Initially, a series of controlled unidirectional freezing experiments were performed on artificial kaolin soil under various freezing conditions to observe the water migration process. Subsequently, a numerical model of water migration was formulated by integrating the partial differential equations of heat and mass transfer. The model's boundary conditions and relevant parameters were derived from both the experimental processes and existing literature. The findings indicate that at lower clay water content, the experimental results align closely with those of the model. Conversely, at higher water content, the modeled results of frost heaving were less pronounced than the experimental outcomes, and the freezing front advanced more slowly. This discrepancy is attributed to the inability of unfrozen water to penetrate once ice lenses form, causing migrating water to accumulate and freeze at the warmest ice lens front. This results in a higher ice content in the freezing zone than predicted by the model, leading to more significant freezing expansion. Additionally, the experimental observations of pore water pressure under freeze-thaw conditions corresponded well with the trends and peaks projected by the simulation results.