Numerical simulation of unstable suction transients in unsaturated soils: the role of wetting collapse

Fluid injection is one of the major triggering factors of instability in unsaturated soils. Robust simulation tools are therefore required to examine hydrologic transients in such a class of inelastic materials. In this paper, the governing equations imposing the balance of mass and momentum in deformable unsaturated porous media are inspected from an analytical perspective. Our main goal is to examine the role of volume collapse on the stability of suction transients. Distinct scenarios are considered in terms of material behavior, distinguishing the case of elastoplastic and viscoplastic materials. It shows that in elastoplastic materials the mathematical ill-posedness of the partial differential equations leads to a diffusive instability of the pore pressure field, signaled by the loss of stress controllability under water injection. By contrast, the introduction of viscosity restores the mathematical well-posedness and ensures the positive diffusivity. Numerical implications of these findings have been explored through simulations of water injection with a 1D finite element solver. It shows that an increase of the soil hydraulic sensitivity (i.e., a marked deterioration of the yield stress upon suction removal) can deteriorate the diffusive stability and prevent the computation of suction changes across a collapsible soil layer. While the incorporation of minor amounts of viscosity can restore numerical stability by suppressing the constitutive singularity. Such numerical results have been corroborated by the computation of local stability indicators based on controllability theory, which proved to be robust diagnostic tools to identify the culprit of runaway instability emerging from coupled hydro-mechanical processes.