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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock generated vorticity in spark discharges

Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 μs), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the shock wave for a range of electrical energies to measure the shock properties and estimate shock velocity and curvature. These measurements are combined with a model to calculate the vorticity generated. The measurements show that the highest vorticity is generated near the peak shock curvature location, and the shock curvature and strength increase with electrical energy deposited. A comparison of the vorticity estimated from the model to vorticity measurements from stereoscopic particle image velocimetry shows the results to be statistically equivalent. Furthermore, this suggests that the shock curvature and velocity contribute to the vortex rings induced by spark discharges.

42 ENGINEERING↗
DOE OSTI2024

A generalization of the shock invariant relationship

Shock invariant relationship, which was conceived for inert shock waves to derive the 4th power relationship between shock pressure and maximum strain rate, is generalized for reactive shock waves such as Chapman–Jouget detonation and shock-induced vaporization. The generalization, based on the first-order reaction models, is a power function relationship between overall dissipated energy (Δe dis ) and reaction time Δτ such that Δe dis Δτ 1/α = constant, where the power coefficient α is found to be in the range of 2/3–4. Experimental data, though scarce, are consistent with the generalization. Implication of the generalization for inert shocks is also considered and suggests a broad range of the 4th power coefficient including an inequality equation that constrains the shock and particle velocity relationship.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

On various mechanisms of shear wave generation from underground chemical explosions in hard rocks

Significant shear waves observed during underground explosions can make them look similar to natural earthquakes. In the near field, they can be caused by the source geometry and emplacement conditions, heterogeneous and anisotropic properties of the rock mass, relaxation of a non-uniform in situ stress around the source after the departure of the shock wave as well as by activation of nearby faults. Here, the stress relaxation mechanism is important at depth when the difference in main principal stresses becomes significant. Another mechanism of shear motion generation is due to explosive gas flow from the source cavity into the cracks initiated by the shock wave and joints connected to the cavity. Also, the shock wave from the explosion can cause weakening of the rock either due to permanent damage or shock-induced liquefaction at saturated joints. Numerical modelling of underground chemical explosions in hard rock at various confinement stresses has been conducted to compare various mechanisms of shear wave generation and their seismic characteristics with that of a simple fault rupture process observed during an earthquake.

58 GEOSCIENCES↗
DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗
DOE OSTI2022

Multi-frame, ultrafast, x-ray microscope for imaging shockwave dynamics

Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Weak shock compaction on granular salt

This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Shock Compression of Fluorapatite to 120 GPa

Apatite is a phosphate mineral relevant to shock metamorphism in planetary materials. Here, we report on the response of natural fluorapatite from Durango, Mexico under shock wave loading between 14.5 and 119.5 GPa. Wave profile measurements were obtained in plate-impact experiments conducted on [0001]-oriented fluorapatite single crystals. To GPa peak stresses, we observed a two-wave structure indicating an elastic-inelastic response with elastic wave amplitudes of 10.5–13.1 GPa. Between 39.1 and 62.1 GPa, a complex wave structure was observed involving the propagation of three waves. At and above 73.7 GPa, only a single shock wave was observed. The data above 73.7 GPa provided the following linear shock velocity – particle velocity relationship: U s = 6.5(2) + 0.78(6) u p , (mm/μs). Above 80 GPa, the densities in the shocked state exceed both the extrapolated 300-K density of fluorapatite and the predicted 300-K density for a mixture of the high-pressure assemblage, tuite and CaF 2 . Furthermore, this result indicates that fluorapatite undergoes a transition to a denser structure under shock loading at these conditions. The shock response of fluorapatite is observed to be similar to enstatite but stiffer than quartz and albite at the stresses examined in this work.

58 GEOSCIENCES↗
DOE OSTI2021

Exceptionally high spallation strength for a high-entropy alloy demonstrated by experiments and simulations

High-entropy alloys are materials with an increasing number of technological applications. Amongst them, the Cantor alloy, FeMnCoCrNi, shows desirable mechanical properties at normal loading conditions. In this study we focus on the performance of the Cantor alloy at the ultra-high deformation rates of shock waves. We study shock-induced spallation using both experiments and atomistic simulations. Experimental loading is achieved using high power laser, with VISAR to obtain velocity profiles and spall strength, followed by transmission electron microscopy of the recovered samples. Molecular Dynamics (MD) simulations of shock-induced spallation are compared with experiments. Both experiments and simulations show a high spall strength which would be beneficial for certain applications, with experiments giving ~8 GPa at ~10 7 s -1 and MD giving almost ~30 GPa at ~10 9 s -1 . The difference between experiments and simulations can be explained by the difference in strain rate. Post-mortem analysis of the experimental samples shows nanotwins near the spall plane, while MD simulations show a highly disordered region giving rise to void nucleation and spall during loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics

To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase transitions in Zr at sub-nanosecond time scales

Here solid-solid phase transitions are investigated in Zr thin films using shock compression induced by a short laser pulse (<1 ns). Shock wave profiles are measured at free surfaces for films of different thicknesses (a few micrometers) using chirped-pulse line velocimetry with 10 ps time resolution. Experiments are performed at pressures up to ~50 GPa, which is sufficient to reach the ω and β phases under equilibrium conditions. The shock wave structures are analyzed using a general Lagrangian analysis method, which allows for the calculation of stress–strain paths and assessments of phase transition behavior. In agreement with recent short laser pulse experiments using ultra-fast x-ray diffraction, we do not find any clear evidence of the α–ω transition, though this would be expected from the phase diagram. Instead, we infer a direct transformation to a metastable β-phase at lower shock pressures (<20 GPa) and equilibrium β at higher pressures. Through the velocimetry analysis, we find α–β transformation onset times of less than ~100 ps and completion times of less than ~200 ps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Non-local entanglement and fast scrambling in de-Sitter holography

We study holographic entanglement and information scrambling in de-Sitter (dS) space in the context of the DS/dS correspondence. We find that our previously identified non-local entanglement structure of dS vacua can be extended out of the time-reflection symmetric slice. We extend the geometry to a two-sided configuration and calculate the zero-time mutual information between two intervals on different sides when there is a localized shock wave in the bulk. Interestingly, we find that the information scrambling time saturates the fast scrambler bound proposed by Sekino and Susskind and that the shock wave renders a wormhole to be traversable. Furthermore, we calculate a two-sided out-of-time-ordered correlator (OTOC) in the late time regime and we see that, before scrambling, it exponentially grows with an exponent whose value saturates the maximal bound of chaos proposed by Maldacena, Shenker and Stanford. At the end, we provide an explanation why the exponential growing of the late-time OTOC with the maximal bound of chaos saturated and the traversability of the wormhole are simple results of the non-local entanglement structure and we point out that this is a realization of the ER=EPR proposal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Shock Damage Analysis in Serial Femtosecond Crystallography Data Collected at MHz X-ray Free-Electron Lasers

Serial femtosecond crystallography (SFX) data were recorded at the European X-ray free-electron laser facility (EuXFEL) with protein microcrystals delivered via a microscopic liquid jet. An XFEL beam striking such a jet may launch supersonic shock waves up the jet, compromising the oncoming sample. To investigate this efficiently, we employed a novel XFEL pulse pattern to nominally expose the sample to between zero and four shock waves before being probed. Analyzing hit rate, indexing rate, and resolution for diffraction data recorded at MHz pulse rates, we found no evidence of damage. Notably, however, this conclusion could only be drawn after careful identification and assimilation of numerous interrelated experimental factors, which we describe in detail. Failure to do so would have led to an erroneous conclusion. Femtosecond photography of the sample-carrying jet revealed critically different jet behavior from that of all homogeneous liquid jets studied to date in this manner.

36 MATERIALS SCIENCE↗