The Coupled Thermo-Hydro-Mechanical Numerical Analysis for the Estimation of Active Zone and Free-Field Heave in Expansive Natural Clay and Bedrock

Expansive soils are widespread in the Southern, Western, and Rocky Mountain areas of the United States. Expansive natural clays and bedrock can expand upon wetting and shrink upon drying. These soils can cause damage to structures, roadways, and pipelines. Accurate estimation of volumetric deformation of expansive soils is critical to minimizing damage to the nation’s infrastructure. Free-field heave is the amount of movement the ground subsurface experiences due to wetting/drying of the soils within the active zone, with no surface surcharge load (e.g., foundation load) applied. Active zone is the depth below the ground where the soil volume deformation can still occur due to moisture/energy changes in unsaturated soils. Several methods for free-field heave estimation have been developed over the years including empirical methods and using oedometer swell test results and soil suction values. Further, in the past decade, most numerical studies have been performed to determine the free-field heave by incorporating hydro-mechanical framework which disregards the temperature effect and usually considers one-way coupling between the moisture flow and swelling/shrinkage deformation. The main objective of the current study is utilizing a strongly coupled thermo-hydro-mechanical (THM) model for the estimation of active zone depth and free-field heave when the soil is subjected to real-time climate conditions. The THM model, input thermal, and hydraulic/mechanical data are implemented in COMSOL Multiphysics software. The results are compared with experimental data from a case study in the Colorado Front Range area. The Soil-Water Characteristic Curve (SWCC) and Hydraulic Conductivity Function (HCF) are based on the Van Genuchten (VG) model and the mechanical constitutive model, which governs swelling deformation, is based on the Extended Barcelona Model. The THM model predicts the active zone reasonably well. Moreover, the simulation of free-field heave is compared with the analytical calculation produced by VOLFLO software.