Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock wave”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗
DOE OSTI2021

Pulse-burst spontaneous Raman thermometry of unsteady wave phenomena in a shock tube

A high-speed temperature diagnostic based on spontaneous Raman scattering (SRS) was demonstrated using a pulse-burst laser. The technique was first benchmarked in near-adiabatic H 2 -air flames at a data-acquisition rate of 5 kHz using an integrated pulse energy of 1.0 J per realization. Both the measurement precision and accuracy in the flame were within 3% of adiabatic predictions. Furthermore, this technique was then evaluated in a challenging free-piston shock tube environment operated at a shock Mach number of 3.5. SRS thermometry resolved the temperature in post-incident and post-reflected shock flows at a repetition rate of 3 kHz and clearly showed cooling associated with driver expansion waves. Collectively, this Letter represents a major advancement for SRS in impulsive facilities, which had previously been limited to steady state regions or single-shot acquisition.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2022

Shocked silica aerogel radiance transition

Silica (SiO 2 ) aerogel is widely used in high-energy-density shock experiments due to its low and adjustable density. Reported here are measurements of the shock velocity, optical radiance, and reflectivity of shocked SiO 2 aerogel with initial densities of 0.1, 0.2, and 0.3 g/cm 3 . These results are compared with similar data from three solid polymorphs of SiO 2 , silica, quartz, and stishovite with initial densities 2.2, 2.65, and 4.3 g/cm 3 , respectively. Interestingly, below a brightness temperature of T bright ≈ 35,000 K, the slope of the radiance vs shock velocity is the same for each of the SiO 2 aerogels and solid polymorphs. At T bright ≈ 35000 K, there is an abrupt change in the radiance vs shock velocity slope for aerogels, but not seen in the solid polymorphs over the pressures and temperatures explored here. An empirical model of shock front radiance as a function of SiO 2 density and laser drive parameters is reported to aid in the design of experiments requiring maximum shock front radiance.

42 ENGINEERING↗
DOE OSTI2025

Rippled shock propagation in a laser-driven target at multimegabar pressures

The evolution of non-uniform shocks produced by modulated laser irradiation or surface perturbations is relevant to studies of inertial confinement fusion and material properties at high-energy-density conditions. We present results from an experiment conducted at the OMEGA EP laser facility, where a 300 GPa shock was driven into a fused silica sample with pre-fabricated single-mode surface modulations. Using time-resolved optical velocimetry, we captured the continuous evolution of rippled shock motion, enabling a comprehensive mapping of the spatial amplitude history from formation to phase reversal in a single experiment. Initially, the ablation-driven shock inherits a fraction of the surface modulation amplitude from the sample, which subsequently grows before decaying, ultimately leading to the flattening of the rippled shock and a phase reversal. We find that two-dimensional inviscid hydrodynamic simulation of the experiment is able to qualitatively capture many aspects of the rippled shock evolution but over-predicts the initial amplitude growth. This experimental platform, capable of accommodating varying ripple wavelengths, lays the groundwork for a potential viscometry method at extreme pressures, where viscous effects manifest as differences in shock flattening times between rippled shocks of two distinct wavelengths propagating through the sample.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Detonation-induced transformation of graphite to hexagonal diamond

In this work, we explore the structural evolution of highly oriented pyrolytic graphite (HOPG) under detonation-induced shock conditions using in situ synchrotron x-ray diffraction in the ns timescale. We observe the formation of hexagonal diamond (lonsdaleite) at pressures above 50 GPa, in qualitative agreement with recent gas gun experiments. First-principles density functional calculations reveal that under uniaxial compression, the energy barrier for the transition toward hexagonal diamond is lower than that for cubic diamond. Finally, no indication of cubic diamond formation was observed up to >70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Polysulfone shock compressed above the decomposition threshold: Velocimetry and modeling of two-wave structures

Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

High-Fidelity Analysis of ULF Wave Mode Structure Following Interplanetary Shock Compression of the Dayside Magnetopause Using MMS Multi-Point Observations

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗