skip to main content

This content will become publicly available on March 14, 2025

Title: A MPM Lagrangian‐Eulerian hydrocode for simulating buried explosions in transversely isotropic geomaterials

Shock waves in geological materials are characterized by a sudden release of rapidly expanding gas, liquid, and solid particles. These shock waves may occur due to explosive volcanic eruptions or be artificially triggered. In fact, underground explosions have often been used as an engineering solution for large‐scale excavation, stimulating oil and gas recovery, creating cavities for underground waste storage, and even extinguishing gas field fires. As such, hydrocodes capable of simulating the rapid and significant deformation under extreme conditions can be a valuable tool for ensuring the safety of the explosions. Nevertheless, as most of the hydrocodes are often formulated in an Eulerian grid, this setting makes it non‐trivial to track the deformation configuration of the materials without a level set. The objective of this paper is to propose the use of the material point method equipped with appropriate equation of state (EOS) models as a hydrocode suitable to simulate underground explosions of transverse isotropic geomaterials. To capture the anisotropic effect of the common layered soil deposits, we introduce a new MPM hydrocode where an anisotropic version of the Mie‐Gruneisen EOS is coupled with a frictional Drucker‐Prager plasticity model to replicate the high‐strain‐rate constitutive responses of soil. By leveraging the Lagrangian nature of material points to capture the historical dependence and the Eulerian calculation of internal force, the resultant model is capable of simulating the rapid evolution of geometry of the soil as well as the high‐strain‐rate soil mechanics of anisotropic materials.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
International Journal for Numerical and Analytical Methods in Geomechanics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)

    Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with a nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion.

    more » « less
  2. Abstract

    The core of this research is separated into three domains, the ultrahigh strain rate response of elastomeric polymers, laser‐induced shock waves , and terahertz time‐domain spectroscopy (THz‐TDS). Elastomers, e.g., polyurea, constitute an advance class of materials suitable for many applications, specifically in high impact loading scenarios, thus, a laser‐induced shock wave (LSW) experimental technique is used to investigate the mechanical response of shock‐loaded polyurea. LSW can submit samples to a strain rate exceeding 106s−1at low strains, enabling determination of material intrinsic failure modes. The large deformation induced during shock loading may alter the macromolecule structure, which can only be detected spectroscopically. Therefore, this research incorporated terahertz bulk spectroscopy to detect and report molecular conformational changes. Microscopy techniques were also used to elucidate changes in the microscale properties, morphology, and topography. The interpretation of the results explicated brittle failure in terms of partial and total spallation and, remarkably, ductile failure leading to plastic deformation, including plastic bulging and adiabatic shearing, not previously associated with LSW technique. Furthermore, spectral changes found in the terahertz regime substantiated the validity of terahertz spectroscopy in elucidating the underlying mechanism associated with the impact mitigating properties of dynamically loaded polyurea.

    more » « less
  3. A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [ O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ∼30% lateral expansion and ∼25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The results revealed that lighter particles such as aluminum gained significant momentum, reaching up to ∼96% of the shocked CH's speed, compared to ∼29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models. 
    more » « less

    The hadron-quark phase transition in quantum chromodynamics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\tt {-}100\, \text{M}_\odot$. We find explosions only for two low-compactness models (14 and $16\, \text{M}_\odot$) with the DD2F_SF EoS, both with low explosion energies of ${\sim }10^{50}\, \mathrm{erg}$. These weak explosions are characterized by a neutrino signal with several minibursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition-driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as 90Zr for the $14\hbox{-} \text{M}_\odot$ zero-metallicity model and 94Zr for the $16\hbox{-}\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the protocompact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to 4kB/baryon and may have a distinctive gravitational-wave signature, also for a second-order phase transition.

    more » « less

    The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones of reacted shale.


    Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.

    more » « less