Geomechanics and Geotechnics

Unsaturated Soil Evolution

In order to determine the soil-water retention (SWR) behavior of a particulate medium, the invading phase pressure in the inter-particle level is correlated to the governing effective pore area in the wetting and drying paths. In a three-phase medium that consists of air, wetting fluid and solids, the invading phase on the drying path is air, whereas on the wetting path the wetting fluid advances into the cavities. On a drying path where the area of a cavity is minimum, the air-entry pressure (AEP) of a pore throat is determined by numerically solving the Young-Laplace curvature equation. This can be done using the finite difference method and Newton-Raphson (Jacobian) approximation technique. Next, a relation between the pore area and the value of AEP is developed by varying the distance between solids around the pore throat. Similarly, the water-entry pressure (WEP) is correlated to a maximum pore area of cavity. After packing the particulate domain with the given particle size distribution (PSD) and void ratio values, the primary/main drying and wetting paths of the wetting fluid are simulated and the effect of hysteresis in SWR is shown. It is considered that the total suction equals to matric suction value and the water bridges between two adjacent particles are formed in the form of pendular rings. In this study, the considered material is non-plastic and the shrinkage and swelling during the wetting and drying phases or any change in pore structure are neglected*.




The packing algorithm and the meniscus in interparticel scale*
The hysteresis in sandy sample during initial drying, wetting and main drying conditions*

Download files


  • (*) Sattari, A.S., Siebelist, M., Toker, N.K., Wuttke, F. (2020). Simulation of hysteresis in soil–water retention with a correlation between the invading pressure to the governing effective flow area. Granular Matter 22, 44. Springer
  • Sattari, A. S. and Toker, N.K. (2016). Obtaining Soil-Water Characteristic Curves by Numerical Modeling of Drainage in Particulate Media. Computers and Geotechnics Journal, Vol. 74, 196-210. ScienceDirect