Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK WAVES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2020

Shocking Cerium into a New Phase

A material's equation of state (EOS) - which indicates its state under various conditions of pressure, temperature, volume, and energy, and relates them to various phases - is a vital part of understanding the material’s behavior and response in different environments and applications. However, determining the EOS and complete phase diagrams of some materials can pose daunting challenges. One of these materials is cerium, for which some regions of the phase diagram have proven elusive. Investigators from the U.S. Department of Energy’s (DOE’s) Los Alamos National Laboratory (LANL) employed the U.S. DOE’s Advanced Photon Source (APS) to probe the high-pressure solid phase of cerium through shock-wave experiments that provided a detailed look at the cerium's transition from the α-ϵ phases. Their work, which provides the first evidence that an α-ϵ phase transition can be shock-induced in cerium, was published in the Journal of Applied Physics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Collapse of helium-filled voids in extreme deformation: Dislocation mechanisms

The mechanisms responsible for the collapse of helium-filled bubbles during the passage of shock waves in monocrystalline copper are revealed. Both internal pressure (caused by pre-existing helium atoms) and bubble size are varied in molecular dynamics simulations to understand the atomistic scale deformation as they are subjected to shock compression at pressures of 48, 123, and 170 GPa, corresponding to particle velocities of 1.0, 2.0, and 2.5 km/s. Both empty and helium filled bubbles serve as dislocation sources, generating intense, localized plastic regions. There are distinct differences in the collapse of empty voids compared to He-filled bubbles, the former requiring less stress and generating a greater density of dislocations for a given shock strength. Additionally, a generalized model for dislocation emission is proposed, where the inclusion of shear stress generated by the helium bubble increases the critical stress to generate dislocations at the defect surface, demonstrating the change in plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ultrafast structural response of shock‐compressed plagioclase

Abstract Meteor impacts can induce unique pressure‐dependent structural changes in minerals due to the propagation of shock waves. Plagioclase—ubiquitous throughout the Earth’s crust, extraterrestrial bodies, and meteorites—is commonly used for reconstructing the impact history and conditions of the parent bodies. However, there have been unresolved inconsistencies in the interpretation of shock transformations across previous studies: The pressure at which amorphization begins and the process by which it occurs is the subject of ongoing debate. Here, we utilize time‐resolved in situ X‐ray diffraction (XRD) to probe the phase transformation pathway of plagioclase during shock compression at a sub‐nanosecond timescale. Direct amorphization begins at pressures much lower than what was previously assumed, just above the Hugoniot elastic limit of 5 GPa, with full amorphization to a high‐density amorphous phase, observed at 32(10) GPa and 20 ns. Upon release, the material partially recrystallizes back into the original structure, demonstrating a memory effect.

58 GEOSCIENCES↗
DOE OSTI2026

A wide-range thermodynamically consistent constitutive model for dynamic loading of partially saturated porous media

SUMMARY This paper presents a novel thermodynamically consistent constitutive model for partially saturated porous rocks across a wide range of conditions. The material states generated behind the shock wave from an explosive source can vary significantly, ranging from crushed and melted rock near the source to a poroelastic medium in the far field. In the model, rock strength is determined by the effective pressure, which is calculated using two independent equations of state: one for the solid rock and another for the pore fluid. The model accounts for shock-induced liquefaction resulting from fluid pressure buildup in the pore spaces near the explosive source. Simultaneously, it describes the increase in wave propagation speed due to elastic pore contraction in both dry and partially saturated rocks. This model is applied to investigate how fluid saturation affects the amplitude and shape of the generated waves, as well as the residual stress surrounding the cavity formed by spherical explosions.

Vorobiev, Oleg [Lawrence Livermore National Labora↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
DOE OSTI2023

Electron–Ion Temperature Ratio in Astrophysical Shocks

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T e /T p inferred from the Hα profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T e /T p for supernova remnant shocks using Hα profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the Hα determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks, and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence, and varying preshock conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shock compression of molybdenum single crystals to 110 GPa: Elastic–plastic deformation and crystal anisotropy

To explore the role of crystal anisotropy on the elastic-plastic deformation of BCC single crystals at high shock stresses, molybdenum (Mo) single crystals were shock compressed along [100], [111], and [110] orientations at elastic impact stresses between 20 and 110 GPa. Laser interferometry was used to measure shock wave velocities and particle velocity histories. Along the [100] and [111] orientations, elastic-plastic waves (two wave structure) were observed up to 110 GPa. Along the [110] orientation, the two wave structure was observed only up to 90 GPa. The measured elastic wave amplitudes were analyzed to determine crystal anisotropy effects, impact stress dependence, and the activated slip systems on the Hugoniot elastic limit. The results from our work have provided insight into the role of crystal anisotropy on the elastic-plastic deformation under shock compression at high stresses.

42 ENGINEERING↗
DOE OSTI2020

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Turbulence, Structures, and Diffusive Shock Acceleration

Diffusive shock acceleration (DSA) is thought to be the primary acceleration at shock waves. While this acceleration mechanism has been very successful in explaining many observations of energetic particles in astrophysics and heliophysics, recent observations suggest that DSA, at least its simplest version, may not be adequate in explaining observations of acceleration of particles at the termination shock and some interplanetary shocks. Meanwhile, turbulence and its structures like flux ropes and current sheets have been proposed to have important effects on accelerating particles into high energies. The goal of this project is to carefully examine and extend theories of particle acceleration including the DSA and recently developed reconnection acceleration associated with turbulence structures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2020

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Visualization of shocked material instabilities using a fast-framing camera and XFEL four-pulse train

Many questions regarding dynamic materials could be answered by using time-resolved ultra-fast imaging techniques to characterize the physical and chemical behavior of materials in extreme conditions and their evolution on the nanosecond scale. In this work, we perform multi-frame phase-contrast imaging (PCI) of micro-voids in low density polymers under laser-driven shock compression. At the Matter in Extreme Conditions (MEC) Instrument at the Linac Coherent Light Source (LCLS), we used a train of four x-ray free electron laser (XFEL) pulses to probe the evolution of the samples. To visualize the void and shock wave interaction, here, we deployed the Icarus V2 detector to record up to four XFEL pulses, separated by 1-3 nanoseconds. In this work, we image elastic waves interacting with the micro-voids at a pressure of several GPa. Monitoring how the material’s heterogeneities, like micro-voids, dictate its response to a compressive wave is important for benchmarking the performances of inertial confinement fusion energy materials. For the first time in a single sample, we have combined an ultrafast x-ray framing camera and four XFEL pulse train to create an ultrafast movie of micro-void evolution under laser-driven shock compression. Eventually, we hope this technique will resolve the material density as it evolves dynamically under laser shock compression.

fusion↗
DOE OSTI2020

Shock-induced kinetics and cellular structures of liquid nitromethane detonation

Using a combination of high-speed diagnostics, optical pyrometry, velocimetry and video photography, we examine the spatiotemporal reaction kinetics of a prototypical high-energy explosive, liquid nitromethane (NM), as we drive this high explosive into detonation. The detonations were initiated by powerful shock waves whose durations (4 ns) were shorter than the characteristic time associated with the reaction (7 ns). Simple optical spectroscopy alone cannot characterize the kinetics with high time resolution because the reactants and products are flowing at a high velocity of ~6 km/s (6 µm/ns). Additionally, every detection volume behind the shock front contains material initiated over a range of times. High spatial resolution can be obtained by probing interfaces where the shock enters or breaks out of the NM. In addition, we obtain cellular patterns imprinted on the luminous shock front by the two-stage explosion in NM. These spatiotemporal patterns arise naturally as a result of the asymmetry produced by the moving shock front. In this way the shock front serves as a thin moving intrinsic optical gauge that reports and characterizes the two-stage NM explosion behind the front.

42 ENGINEERING↗
DOE OSTI2022

The Coupling of an EUV Coronal Wave and Ion Acceleration in a Fermi-LAT Behind-the-Limb Solar Flare

We present the Fermi-LAT observations of the behind-the-limb (BTL) flare of 2021 July 17 and the joint detection of this flare by STIX on board the Solar Orbiter. The separation between Earth and the Solar Orbiter was $99^{\circ}_{\centerdot}2$ at 05:00 UT, allowing STIX to have a front view of the flare. The location of the flare was S20E140 in Stonyhurst heliographic coordinates, making this the most distant behind-the-limb flare ever detected in >100 MeV gamma-rays. The LAT detection lasted for ~16 minutes, the peak flux was 3.6 ± 0.8 (10 –5 ) ph cm –2 s –1 with a significance >15σ. A coronal wave was observed from both STEREO-A and SDO in extreme ultraviolet (EUV), with an onset on the visible disk in coincidence with the LAT onset. A complex type II radio burst was observed by GLOSS also in coincidence with the onset of the LAT emission, indicating the presence of a shock wave. We discuss the relation between the time derivative of the EUV wave intensity profile at 193 Å as observed by STEREO-A and the LAT flux to show that the appearance of the coronal wave at the visible disk and the acceleration of protons as traced by the observed >100 MeV gamma-ray emission are coupled. We also report how this coupling is present in the data from three other BTL flares detected by Fermi-LAT, suggesting that the protons driving the gamma-ray emission of BTL solar flares and the coronal wave share a common origin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Use of Acoustic Cavitation to Recovery High-Assay Low-Enriched Uranium (HALEU) form TRISO Fuel Particles

Recent interest in advanced nuclear reactor concepts such as small modular reactor (SMR), micro reactor or versatile test reactor (VTR) has increased the demand for high-assay, low-enriched uranium (HALEU). The sustainable management of this resource is integral to successful resurgence of nuclear energy in the United States. The High Temperature Gas-Cooled Reactor (HTGR) concept is a prime candidate for deployment of the new fleet of SMRs, and it will become very important to have a simple and cost-effective technology to recover HALEU from short-burned, defective, or off-specification TRISO fuel particles. We have attempted to harness the high energetics of sonochemistry to penetrate the TRISO fuel particles and recover HALEU via acidic leaching. While the TRISO fuel particle design is inherently “bulletproof,” (due to the silicon carbide (SiC) layer), ultrasonic irradiation may hold the key to convenient access to its valuable contents. Sonochemistry is a field of chemistry based on acoustic cavitation, which is the formation, growth, and collapse of bubbles in liquid media.1,2 The oscillating bubble formation is produced by irradiation of a liquid media with sound waves. Literature reports indicate that collapsing bubbles induced by cavitation produces intense local heating, high pressures, and short lifetimes.1 These localized hot spots reach temperatures of ˜5000 K, pressures approaching 500 atm, and heating and cooling rates exceeding 100 K/s.1 The large temperature and pressure differentials deliver high-energy heating and microscopic explosive shock waves to a liquid media or liquid/solid interface. Cavitation at the surface of a solid in solution induces a deformation in the bubble cavity upon collapse. This deformation reinforces the bubble structure and sends a fast-moving stream of liquid through the cavity at the surface with velocities greater than 100 m/s.3 These energetic impacts have demonstrated an ability to penetrate or simply destroy the SiC shell surrounding the uranium fuel kernel during prolonged exposures to the high-power acoustic waves. Some preliminary results demonstrating HALEU recovery using sonochemistry techniques will be reported.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Characterizing IHE Response to Multiple Shock Loading

The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗