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):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
7th International Conference on Earthquake Geotechnical Engineering
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

    Seismic design of water retaining structures relies heavily on the response of the retained water to shaking. The water dynamic response has been evaluated by means of analytical, numerical, and experimental approaches. In practice, it is common to use simplified code‐based methods to evaluate the added demands imposed by water sloshing. Yet, such methods were developed with an inherent set of assumptions that might limit their application. Alternatively, numerical modeling methods offer a more accurate way of quantifying the water response and have been commonly validated using 1 g shake table experiments. In this study, a unique series of five centrifuge tests was conducted with the goal of investigating the hydrodynamic behavior of water by varying its height and length. Moreover, sine wave and earthquake motions were applied to examine the water response at different types and levels of excitation. Arbitrary Lagrangian‐Eulerian finite element models were then developed to reproduce 1 g shake table experiments available in the literature in addition to the centrifuge tests conducted in this study. The results of the numerical simulations as well as the simplified and analytical methods were compared to the experimental measurements, in terms of free surface elevation and hydrodynamic pressures, to evaluate their applicability and limitations. The comparison showed that the numerical models were able to reasonably capture the water response of all configurations both under earthquake and sine wave motions. The analytical solutions performed well except for cases with resonance under harmonic motions. As for the simplified methods, they provided acceptable results for the peak responses under earthquake motions. However, under sine wave motions, where convective sloshing is significant, they underpredict the response. Also, beyond peak ground accelerations of 0.5 g., a mild nonlinear increase in peak dynamic pressures was measured which deviates from assumed linear response in the simplified methods. The study confirmed the reliability of numerical models in capturing water dynamic responses, demonstrating their broad applicability for use in complex problems of fluid‐structure‐soil interaction.

    more » « less
  3. 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
  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