skip to main content


Title: 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.  more » « less
Award ID(s):
1728172
NSF-PAR ID:
10104124
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
7th International Conference on Earthquake Geotechnical Engineering
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. One of the key factors in simulating realistic wall-bounded flows at high Reynolds numbers is the selection of an appropriate turbulence model for the steady Reynolds Averaged Navier–Stokes equations (RANS) equations. In this investigation, the performance of several turbulence models was explored for the simulation of steady, compressible, turbulent flow on complex geometries (concave and convex surface curvatures) and unstructured grids. The turbulence models considered were the Spalart–Allmaras model, the Wilcox k- ω model and the Menter shear stress transport (SST) model. The FLITE3D flow solver was employed, which utilizes a stabilized finite volume method with discontinuity capturing. A numerical benchmarking of the different models was performed for classical Computational Fluid Dynamic (CFD) cases, such as supersonic flow over an isothermal flat plate, transonic flow over the RAE2822 airfoil, the ONERA M6 wing and a generic F15 aircraft configuration. Validation was performed by means of available experimental data from the literature as well as high spatial/temporal resolution Direct Numerical Simulation (DNS). For attached or mildly separated flows, the performance of all turbulence models was consistent. However, the contrary was observed in separated flows with recirculation zones. Particularly, the Menter SST model showed the best compromise between accurately describing the physics of the flow and numerical stability. 
    more » « less
  2. Abstract

    The susceptibility of a granular soil to suffusion is strongly dependent on its grain size distribution (GSD) and the mechanical and hydraulic conditions it is subjected to. This study investigates the onset of suffusion considering the effect of confining pressure and stress anisotropy using a fully resolved computational fluid dynamics and discrete element method (CFD–DEM). Three benchmarks, including the sedimentations of single and two adjacent spheres and the classic one‐dimensional (1D) consolidation are performed to demonstrate the capability of this method for high‐fidelity particle‐fluid simulations. A modified hydraulic criterion for the onset of suffusion considering stress anisotropy is presented. The microstructural changes of soil specimens before and during global suffusion are inspected, with emphasis on the evolutions of particle kinetic energy and displacements, force chain networks, and stress anisotropy. We found that the critical hydraulic gradient is negatively correlated with the confining pressure and the degree of stress anisotropy. Fine particles in the soil matrix are locally detached at small hydraulic gradients before the apparent global suffusion, as manifested by the variation of particle kinetic energy and coordination numbers. The roles of different contact types on force transmission and stress anisotropy in eroded specimens are also examined.

     
    more » « less
  3. The standard of practice when assessing the liquefaction susceptibility of geosystems uses an empirical case history database that was primarily developed for clean, poorly graded sands. However, many geosystems in the built environment are either constructed with or founded on well graded soils, creating a disconnect between the sand encountered in practice and the sand used as the basis of knowledge. Using the 9-m centrifuge at the University of California Davis’s Center for Geotechnical Modeling a centrifuge experiment was designed to test the dynamic response of embankments constructed poorly graded and well graded sands at the system level scale. The experiment consisted of two 10-degree slopes, one constructed with a poorly graded sand and the other with a well graded sand positioned side by side in the same model container. Each slope was dry pluviated to the same relative density of Dr=63%, while the absolute densities were different. The slopes were instrumented with dense arrays of pore pressure transducers and accelerometers in the level ground at the head of the slope. The stress-strain behavior between accelerometers was calculated using inverse analysis techniques, providing a 1-D shear-beam soil response at the sensor array location. Liquefaction was triggered, as defined by an excess porewater pressure ratio (ru) of 1.0, but the shear strains at triggering in the well graded sand were significantly less than the strains in the poorly graded sand. During cyclic mobility, strain accumulation in the well graded sand occurred at a slower rate. This study demonstrates that liquefaction triggering and the post-triggering response for saturated sands needs to consider gradation characteristics and clean poorly graded sands cannot act as a single predictor of dynamic response for all sand gradations. 
    more » « less
  4. null (Ed.)
    The dynamic stall phenomenon produces adverse aerodynamic loading, which negatively affects the structural strength and life of aerodynamic systems. Aerodynamic shape optimization (ASO) provides a practical approach for delaying and mitigating dynamic stall characteristics without the addition of an auxiliary system. A typical ASO investigation requires multiple evaluations of accurate but time-consuming computational fluid dynamics (CFD) simulations. In the case of dynamic stall, unsteady CFD simulations are required for airfoil shape evaluation; combining it with high-dimensions of airfoil shape parameterization renders the ASO investigation computationally costly. In this study, metamodel-based optimization (MBO) is proposed using the multifidelity modeling (MFM) technique to efficiently conduct ASO investigation for computationally expensive dynamic stall cases. MFM methods combine data from accurate high-fidelity (HF) simulations and fast low-fidelity (LF) simulations to provide accurate and fast predictions. In particular, Cokriging regression is used for approximating the objective and constraint functions. The airfoil shape is parameterized using six PARSEC parameters. The objective and constraint functions are evaluated for a sinusoidally oscillating airfoil with the unsteady Reynolds-averaged Navier-Stokes equations at a Reynolds number of 135,000, Mach number of 0.1, and reduced frequency of 0.05. The initial metamodel is generated using 220 LF and 20 HF samples. The metamodel is then sequentially refined using the expected improvement infill criteria and validated with the normalized root mean square error. The refined metamodel is utilized for finding the optimal design. The optimal airfoil shape shows higher thickness, larger leading-edge radius, and an aft camber compared to baseline (NACA 0012). The optimal shape delays the dynamic stall occurrence by 3 degrees and reduces the peak aerodynamic coefficients. The performance of the MFM method is also compared with the single-fidelity metamodeling method using HF samples. Both the approaches produced similar optimal shapes; however, the optimal shape from MFM achieved a minimum objective function value while more closely satisfying the constraint at a computational cost saving of around 41%. 
    more » « less
  5. We examine linear and nonlinear shear and extensional rheological properties using a “micelle-slip-spring model” [T. Sato et al., J. Rheol. 64, 1045–1061 (2020)] that incorporates breakage and rejoining events into the slip-spring model originally developed by Likhtman [Macromolecules 38, 6128–6139 (2005)] for unbreakable polymers. We here employ the Fraenkel potential for main chain springs and slip-springs to address the effect of finite extensibility. Moreover, to improve extensional properties under a strong extensional flow, stress-induced micelle breakage (SIMB) is incorporated into the micelle-slip-spring model. Thus, this model is the first model that includes the entanglement constraint, Rouse modes, finite extensibility, breakage and rejoining events, and stress-induced micelle breakage. Computational expense currently limits the model to micellar solutions with moderate numbers of entanglements ([Formula: see text]), but for such solutions, nearly quantitative agreement is attained for the start-up of the shearing flow. The model in the extensional flow cannot yet be tested owing to the lack of data for this entanglement level. The transient and steady shear properties predicted by the micelle-slip-spring model for a moderate shear rate region without significant chain stretch are fit well by the Giesekus model but not by the Phan–Thien/Tanner (PTT) model, which is consistent with the ability of the Giesekus model to match experimental shear data. The extensional viscosities obtained by the micelle-slip-spring model with SIMB show thickening followed by thinning, which is in qualitative agreement with experimental trends. Additionally, the extensional rheological properties of the micelle-slip-spring model with or without SIMB are poorly predicted by both the Giesekus and the PTT models using a single nonlinear parameter. Thus, future work should seek a constitutive model able to capture the behavior of the slip-spring model in shear and extensional flows and so provide an accurate, efficient model of micellar solution rheology. 
    more » « less