Numerical Simulations of Coupled Vapor, Water, and Heat Flow in Unsaturated Deformable Soils During Freezing and Thawing

Freezing and thawing cycles significantly affect the mechanical and hydraulic behavior of soils, posing detrimental challenges for infrastructures in cold climates. This study develops and validates a coupled Thermal–Hydraulic–Mechanical (THM) model using COMSOL Multiphysics (Version 6.3) to demonstrate the complexities of vapor and water flux, heat transport, frost heave, and vertical stress build-up in unsaturated soils. The analysis focuses on fine sand, sandy clay, and silty clay by examining their varying susceptibilities to frost action. Silty clay generated the highest amount of frost heave and steepest vertical stress gradients due to its high-water retention and strong capillary forces. Fine sand, on the other hand, produced a minimal amount of frost heave and a polarized vertical stress distribution. The study also revealed that vapor flux is more noticeable in freezing fine sand, while silty clay produces the greatest water flux between the frozen and unfrozen zones. The study also assesses the impact of soil properties including the saturated hydraulic conductivity, the particle thermal conductivity, and particle heat capacity on the frost-induced phenomena. Findings show that reducing the saturated hydraulic conductivity has a greater impact on mitigating frost heave than other variations in thermal properties. Silty clay is most affected by these changes, particularly near the soil surface, while fine sand shows less noticeable responses.