Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shocks”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing shocks interacting with radiation waves with the Radishock experiment

Both radiation flows and shocks have been extensively studied in the laboratory in the past few decades due to their critical roles in many astrophysical and high-energy density physics processes. In the Radishock experiment, a halfraum-powered radiation wave is driven into a low-density foam and interacts with an ablatively driven, counter-propagating shock. The interacting waves produce a spike in energy density with a temperature greater than the local temperature of the individual waves. As in the successful predecessor experiment, COAX, the primary diagnostic uses absorption spectroscopy at many locations down the cylindrical target, enabling a spatial temperature inference of the radiation wave and its interactions with the shock. Combined with a radiography diagnostic that is capable of imaging the shock and interaction features, we are able to study and inform model predictions of the interaction spike phenomenon. We describe the underlying physics behind the shock interactions with the radiation front and the implications of this experimental study for a broad range of astrophysical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

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 OSTI2024

Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach

Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Quadrilateral particle arrangement within shocks in a two-dimensional dusty plasma

The microscopic structure within a two-dimensional shock was studied using data from a dusty plasma experiment. A single layer of charged microparticles, levitated in a glow-discharge plasma, was perturbed by an electrically floating wire that was moved at a steady supersonic speed to excite a compressional shock. A rearrangement of particles was observed, from a hexagonal lattice in the preshock into a quadrilateral microstructure on the front side of the shock. This quadrilateral structure would not be stable in a monolayer of identical repulsive particles, under equilibrium conditions. Glaser-Clark polygon analysis of the microstructure helped in identifying quadrilaterals. Voronoi analysis was used to characterize the defect fraction behind the shock, as an indication of shock-induced melting. Published by the American Physical Society 2025

Kananovich, Anton (ORCID:0000000225374952)↗
DOE OSTI2024

Flat Spectra of Energetic Particles in Interplanetary Shock Precursors

The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extending upstream of the shock for up to seven flux e-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy as ϵ -1 , consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivity κ on the particle flux F (nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ 3/2 that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Sound speed determination in copper shock compressed to 190 GPa

Sound speed measurements in shock compressed solids have long been valuable for the development of equations of state at extreme conditions, shock-induced phase transformations, and a comprehensive characterization of the thermophysical response of high-pressure standards. We present results from plate impact experiments to 190 GPa to determine the longitudinal sound speed in copper—an important high-pressure standard. Surprisingly, the sound speeds determined using the two most common experimental techniques—the front surface impact (FSI) approach and the release wave overtake (RWO) approach—diverge significantly for stresses greater than ~100 GPa. Further analyses, including numerical simulations, show that the FSI experiments provide the correct sound speeds and that fundamental assumptions underlying the RWO method are likely violated due to the complex release response of shock compressed copper. The sound speeds determined using the FSI approach provide for a more accurate high-pressure description of copper in dynamic compression experiments. The present findings are in contrast to the results for shock compressed silver [Wallace et al., Phys. Rev. B 104, 214106 (2021)], where both methods provided consistent sound speed results. Thus, the findings presented here demonstrate the need to experimentally verify the validity of the RWO method on a case-by-case basis. Finally, we note that even at the high stresses in the present work, the copper unloading response shows a time-dependent, quasielastic response often observed in metals at lower stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
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 OSTI2025

Biermann-Battery-Driven Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

This Letter reports the first complete observation of magnetized collisionless shock precursors formed through the compression of Biermann-battery magnetic fields in laser produced plasmas. At OMEGA, lasers produce a supersonic CH plasma flow which is magnetized with Biermann-battery magnetic fields. The plasma flow collides with an unmagnetized hydrogen gas jet plasma to create a magnetized shock precursor. The situation where the flowing plasma carries the magnetic field is similar to the Venusian bow shock. Imaging 2⁢𝜔 Thomson scattering confirms that the interaction is collisionless and shows density and temperature jumps. Proton radiographs have regions of strong deflections and FLASH magnetohydrodynamic (MHD) simulations show the presence of Biermann fields in the Thomson scattering region. Electrons are accelerated to energies of up to 100 keV in a power-law spectrum. Here, OSIRIS particle-in-cell (PIC) simulations, initialized with measured parameters, show the formation of a magnetized shock precursor and corroborate the experimental observables.

High-energy-density plasmas↗