Numerical modelling of suction diffusion in unsaturated expansive soils under climatic fluctuations

Abstract

The behaviour of expansive clay soils is highly influenced by climatic variations that govern the hydric exchanges between the atmosphere and the ground surface. These fluctuations cause changes in soil suction, leading to significant volumetric variations in expansive soils and resulting in differential settlements of lightweight structures, which can induce structural damage. Despite the critical role of these hydroclimatic processes, few models are able to describe the evolution of soil suction under realistic boundary conditions. In this study, a two-dimensional diffusion model is proposed to simulate the temporal and spatial evolution of suction beneath a structure, considered as impermeable, under the influence of climatic fluctuations and over selected time steps. The model is based on the numerical solution of the diffusion equation in unsaturated porous medium, formulated in terms of suction. The boundary conditions incorporate a surface water balance accounting for infiltration and evaporation, using meteorological data from the Météo-France SIM2 model. Hydraulic soil parameters, such as unsaturated permeability and water retention curve expressing suction as a function of water content, are defined from empirical relationships reported in the literature. The model is applied to a clay soil whose parameters are derived from the Chaingy experimental site (France) using daily climatic data for the year 2025. The simulations highlight a strong seasonal variability of suction concentrated within the upper soil layers, while the presence of an impermeable surface zone representing the structure significantly attenuates suction fluctuations beneath the foundation and induces marked horizontal gradients during dry periods.