Numerical Modeling of the Effect of Desaturation on Liquefaction Hazard Mitigation

Abstract

Earthquake-induced liquefaction is always a concern when the soil near the surface of a site is composed of relatively loose saturated sand. One of liquefaction mitigation methods is to induce gas bubbles into the deposit to reduce the degree of saturation. A coupled pore-scale model is presented herein to investigate liquefaction resistance of desaturated granular materials. The multiphase fluid, which mimics the behavior of air and water, is modeled using the multiphase single component lattice Boltzmann method. The solid phase is modeled using the discrete element method. The coupled framework was utilized to study the behavior of a soil deposit with the different degrees of saturation of 100%, 92%, and 82% during an earthquake loading. Based on the results of the simulations performed, liquefaction occurred in the fully saturated granular deposit and was not observed anywhere at depth in the desaturated deposits. It has also been found that reducing the saturation level from 100% to 92% significantly affects behavior. In desaturated deposits, higher average coordination number, lower pore pressure buildup, and slower effective stress decay were observed compared to fully saturated deposits. However, it turned out that a further reduction in the degree of saturation from 92% to 82% does not have a significant impact on the calculated parameters.