skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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 Tsunamis are rare, extreme events and cause significant damage to coastal infrastructure, which is often exacerbated by soil instability surrounding the structures. Simulating tsunamis in a laboratory setting is important to further understand soil instability induced by tsunami inundation processes. Laboratory simulations are difficult because the scale of such processes is very large, hence dynamic similitude cannot be achieved for small-scale models in traditional water-wave-tank facilities. The ability to control the body force in a centrifuge environment considerably reduces the mismatch in dynamic similitude. We review dynamic similitudes under a centrifuge condition for a fluid domain and a soil domain. A novel centrifuge apparatus specifically designed for exploring the physics of a tsunami-like flow on a soil bed is used to perform experiments. The present 1:40 model represents the equivalent geometric scale of a prototype soil field of 9.6 m deep, 21 m long, and 14.6 m wide. A laboratory facility capable of creating such conditions under the normal gravitational condition does not exist. With the use of a centrifuge, we are now able to simulate and measure tsunami-like loading with sufficiently high water pressure and flow velocities. The pressures and flow velocities in the model are identical to those of the prototype yielding realistic conditions of flow-soil interaction. 
    more » « less
  3. Brehm, Christoph; Pandya, Shishir (Ed.)
    Computational fluid dynamics (CFD) and its uncertainty quantification are computationally expensive. We use Gaussian Process (GP) methods to demonstrate that machine learning can build efficient and accurate surrogate models to replace CFD simulations with significantly reduced computational cost without compromising the physical accuracy. We also demonstrate that both epistemic uncertainty (machine learning model uncertainty) and aleatory uncertainty (randomness in the inputs of CFD) can be accommodated when the machine learning model is used to reveal fluid dynamics. The demonstration is performed by applying simulation of Hagen-Poiseuille and Womersley flows that involve spatial and spatial-tempo responses, respectively. Training points are generated by using the analytical solutions with evenly discretized spatial or spatial-temporal variables. Then GP surrogate models are built using supervised machine learning regression. The error of the GP model is quantified by the estimated epistemic uncertainty. The results are compared with those from GPU-accelerated volumetric lattice Boltzmann simulations. The results indicate that surrogate models can produce accurate fluid dynamics (without CFD simulations) with quantified uncertainty when both epistemic and aleatory uncertainties exist. 
    more » « less
  4. 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
  5. Computational fluid dynamics (CFD) simulations are broadly used in many engineering and physics fields. CFD requires the solution of the Navier–Stokes (N-S) equations under complex flow and boundary conditions. However, applications of CFD simulations are computationally limited by the availability, speed, and parallelism of high-performance computing. To address this, machine learning techniques have been employed to create data-driven approximations for CFD to accelerate computational efficiency. Unfortunately, these methods predominantly depend on large labeled CFD datasets, which are costly to procure at the scale required for robust model development. In response, we introduce a weakly supervised approach that, through a multichannel input capturing boundary and geometric conditions, solves steady-state N-S equations. Our method achieves state-of-the-art results without relying on labeled simulation data, instead using a custom data-driven and physics-informed loss function and small-scale solutions to prime the model for solving the N-S equations. By training stacked models, we enhance resolution and predictability, yielding high-quality numerical solutions to N-S equations without hefty computational demands. Remarkably, our model, being highly adaptable, produces solutions on a 512 × 512 domain in a swift 7 ms, outpacing traditional CFD solvers by a factor of 1,000. This paves the way for real-time predictions on consumer hardware and Internet of Things devices, thereby boosting the scope, speed, and cost-efficiency of solving boundary-value fluid problems. 
    more » « less