Numerical Simulation of Water Migration during Soil Freezing and Its Resulting Characterization

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.