Computational fluid dynamics-based modeling of liquefied soils

The residual shear strength of liquefied soil is a key parameter in evaluating liquefaction flow failures. Results from a series of dynamic centrifuge experiments where the shear strength of liquefied soil was inferred by measuring the force required to pull a thin metal plate (coupon) horizontally through the liquefied soil are assessed here using a computational fluid dynamics (CFD) based model. Viscosity is a key parameter for the Newtonian fluid constitutive model used in the simulations, and apparent viscosities of liquefied soil in the range of about 5,800 – 13,300 Pa·s were obtained when the CFD model was calibrated against coupons pulled through liquefied soil in dynamic centrifuge tests. These computational values agree reasonably with apparent viscosities of liquefied soil reported in the literature when the Reynold’s numbers exceeded 1.0. Importantly, the CFD simulations illustrated that in cases where Reynold’s numbers are < 1.0, apparent viscosities of liquefied soil back-calculated using simplistic closed-form solutions commonly applied in geotechnical literature are several orders of magnitude too large; and therefore, such closed-form solutions should not be used for these cases.