Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock wave experiments”

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

The Computational Analysis of a Generic Hypersonic Store Separation

Many current hypersonic vehicles involve a main vehicle from which components separate during flight and presents simulation challenges such as shock-shock interactions, highly nonlinear interactions, and shock wave boundary layer interactions. These aerodynamic phenomena influence the attitude and trajectory of each object, which needs to be understood before large-scale experiments can be run. To gain a generic understanding of the process, high-fidelity Reynolds-averaged Navier-Stokes solutions are performed on a conical vehicle geometry passing through an oblique shock wave that is representative of different configurations and scenarios. Trajectory and applied forces are tracked and show the vehicle’s dynamics are predominantly the result of differential flow incidence angles causing a strong shock; this leads to a large pressure increase over a fraction of the vehicle, which influences pitch. The vehicle appears to follow conventional stability theory with detached eddy simulation and small disturbances in initial attitude shown to have minimal influence on the scenario.

Savery, Ryan↗
DOE OSTI2021

Time-resolved particle image velocimetry measurements of the turbulent Richtmyer–Meshkov instability

Experiments are presented on the Richtmyer–Meshkov instability (RMI) with a three-dimensional, multi-mode initial perturbation. The experiments use a vertical shock tube, where a stably stratified interface is formed between air and sulphur hexafluoride (SF 6 ) via counterflow. A perturbation is imposed at the interface by vertical oscillation of the gas column, forming Faraday waves. The interface is accelerated by a Mach 1.17 (in air) shock wave, and the development of the mixing region between the gases is investigated using particle image velocimetry. Following shock acceleration, a reflected shock wave from the bottom of the shock tube interacts with the mixing layer a second time (reshock). The experiment is initialized with both high and low amplitude perturbations to examine the effect of the perturbation amplitude on measured quantities. The instability growth exponent (θ) is determined from the kinetic energy field using the width of the mixing layer and the decay of kinetic energy, which are found to be in agreement when the flow is most strongly excited. A growth exponent of θ ≈ 0.5 is found for all cases except the high-amplitude reshocked regime (where θ ≈ 0.33). High-amplitude experiments exhibit the transitional outer Reynolds number (Re ≡ hh˙/ν >10 4 ) required for mixing transition following the incident shock, and both experiments are elevated well above this threshold following reshock. However, neither set of experiments meet the more stringent requirements proposed by Zhou et al. (Phys. Rev.E, vol. 67, issue 5, 2003) which include the time dependent aspect of the RMI, an observation which is also made when examining the spectra.

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

Shocked Transport: Experiments to Study Dynamic Particle Behavior under Varying Volume Fraction Conditions

Shocks interacting with densely-packed porous media is common in volcanic eruptions or in engineering applications such as ordinance. Large-scale, exploratory experiments have been performed by Frost and collaborators examining the behavior of a bed of particles subjected to a detonation wave. The detonation replaces the central explosive with a high-density, high-temperature gas that propagates outward. A shock wave travels through the particle bed, followed by a contact discontinuity between the shocked air and the explosive products. In the early times, the detonation wave causes many of the particles to collide and compact with each other. As the particles disperse, the experiments demonstrate a series of large-scale jets or fingers that form in the advancing particle front.

42 ENGINEERING↗
DOE OSTI2025

Experimental Pathways for Detecting Double Superionicity in Planetary Ices

The ice giant planets Uranus and Neptune are assumed to contain large amounts of planetary ices such as water, methane, and ammonia. The properties of mixtures of such ices at the extreme pressures and temperatures of planetary interiors are not yet well understood. Ab initio computer simulations have predicted that a number of ices exhibit a hydrogen superionic state and a doubly superionic state. Since the latter state has not yet been generated with experiments, we outline here two possible pathways for reaching and detecting such a state with dynamic compression experiments. Here, we suggest X-ray diffraction as the principal tool for detecting when the material becomes doubly superionic and the sublattice of one of the heavy nuclei melts. That would require a temperature of ~3500 K and pressures greater than ~200 GPa for H 3 NO 4 , which we use as an example material here. Such conditions can be reached with experiments that employ an initial shock that is followed by a ramp compression wave. Alternatively, one may use triple-shock compression because a single shock does not yield sufficiently high densities.

Neptune↗
DOE OSTI2021

Low mode implosion symmetry sensitivity in low gas-fill NIF cylindrical hohlraums

Achieving an efficient capsule implosion in National Ignition Facility indirect-drive target experiments requires symmetric hohlraum x-ray drive for the duration of the laser pulse. This is commonly achieved using two-sided two-cone laser irradiation of cylindrical hohlraums that, in principle, can zero the time average of all spherical harmonic asymmetry modes <6 as well as the time dependence of the usually dominant mode 2. In practice, experimental evidence indicates that maintaining symmetric drive becomes limited late in the pulse due to the inward expansion of the hohlraum wall and outward expansion of the capsule ablator plasmas impairing the propagation of the inner-cone laser beams. This effect is enhanced in hohlraums employing low gas-fill, now used almost exclusively as these provide the highest performing implosions and reduce Stimulated Brillouin and Raman backscatter losses, since the gas plasma provides less back pressure to limit blow-in of the hohlraum wall and capsule ablator plasmas. In order to understand this dynamic behavior, we combined multi-keV X-ray imaging of the wall and imploded fuel plasmas as we changed a single parameter at a time: hohlraum gas-fill, laser outer cone picket energy, radius of high density carbon capsules used, and laser beam polar and azimuthal pointing geometry. In this work, we developed a physics-based multi-parameter experimental scaling to explain the results that extend prior scalings and compare those to radiation hydrodynamic simulations to develop a more complete picture of how hohlraum, capsule, and laser parameters affect pole vs equator drive symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Measurements of the inert Hugoniot and observation of mechanical ignition in Ni(V) + Al reactive multilayers via x-ray diffraction

This study investigates the inert Hugoniot response, mechanical ignition, and reaction dynamics of Ni(V)+Al multilayers during longitudinal, laser-driven shock compression experiments. Ni(V)+Al multilayers, known for their self-propagating exothermic reactions, were subjected to longitudinal stresses exceeding 50 GPa using the laser shock facility within the Dynamic Compression Sector (DCS) at the Advanced Photon Source (APS). In situ x-ray diffraction (XRD) revealed that Ni(V) and Al were not in equilibrium during compression, with stress discrepancies attributed to twinning, grain structure effects, and/or dislocation density. However, the measured inert Hugoniot closely matched prior experimental and computational studies, confirming the utility of XRD for measuring the equation of state of thin, complex materials. Additionally, reaction was observed at significantly higher stresses than reported previously using laser-launched flyers. This discrepancy suggests a strong influence of externally imposed shear stress on reaction thresholds, which likely arose from deviations in flyer planarity during past experiments. Full reaction of the multilayer occurred within 40 ns after shock-wave passage, evidenced by complete melting of the constituents. Eulerian hydrocode simulations replicated experimental conditions, providing insights into equilibrium dynamics and experimental artifacts. The results highlight how even small shear forces facilitate ignition in Ni(V)+Al multilayers at lower stresses.

Crystallography↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2025

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗
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 OSTI2022

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2020

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2020

Mixing in ICF implosions on the National Ignition Facility caused by the fill-tube

The micrometer-scale tube that fills capsules with thermonuclear fuel in inertial confinement fusion experiments at the National Ignition Facility is also one of the implosion's main degradation sources. It seeds a perturbation that injects the ablator material into the center, radiating away some of the hot-spot energy. This paper discusses how the perturbation arises in experiments using high-density carbon ablators and how the ablator mix interacts once it enters the hot-spot. Both modeling and experiments show an in-flight areal-density perturbation and localized x-ray emission at stagnation from the fill-tube. Simulations suggest that the fill-tube is degrading an otherwise 1D implosion by ~2×, but when other degradation sources are present, the yield reduction is closer to 20%. Characteristics of the fill-tube assembly, such as the through-hole size and the glue mass, alter the dynamics and magnitude of the degradation. Overall, these aspects point the way toward improvements in the design, some of which (smaller diameter fill-tube) have already shown improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗