Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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

Shock compression of a low-density carbon foam

Porous polymer foams are found in numerous applications where they may be subjected to dynamic compression. The shockwave properties of polymer foams have several features that make shock Hugoniot measurements and interpretation challenging. Here, we present a comparison of dynamic X-ray phase contrast imaging and traditional optical velocimetry approaches to measuring the shocked states for a low-density carbon foam. Overall, dynamic X-ray phase contrast imaging was found to reduce errors associated with the measured Hugoniot states. In this study, the foam exhibits anomalous compression as expected due to its high initial porosity arising from shock heating.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Compensatory reserve detects subclinical shock with more expeditious prediction for need of life‐saving interventions compared to systolic blood pressure and blood lactate

Abstract Introduction We conducted a prospective observational study on 205 trauma patients at a level I trauma facility to test the hypothesis that a compensatory reserve measurement (CRM) would identify higher risk for progression to shock and/or need a life‐saving interventions (LSIs) earlier than systolic blood pressure (SBP) and blood lactate (LAC). Methods A composite outcome metric included blood transfusion, procedural LSI, and mortality. Discrete measures assessed as abnormal (ab) were SBP <90 mmHg, CRM <60%, and LAC >2.0. A graded categorization of shock was defined as: no shock (normal [n] SBP [n‐SBP], n‐CRM, n‐LAC); sub‐clinical shock (ab‐CRM, n‐SBP, n‐LAC); occult shock (n‐SBP, ab‐CRM, ab‐LAC); or overt shock (ab‐SBP, ab‐CRM, ab‐LAC). Results Three patients displayed overt shock, 53 displayed sub‐clinical shock, and 149 displayed no shock. After incorporating lactate into the analysis, 86 patients demonstrated no shock, 25 were classified as sub‐clinical shock, 91 were classified as occult shock, and 3 were characterized as overt shock. Each shock subcategory revealed a graded increase requiring LSI and transfusion. Initial CRM was associated with progression to shock (odds ratio = 0.97; p < .001) at an earlier time than SBP or LAC. Conclusions Initial CRM uncovers a clinically relevant subset of patients who are not detected by SBP and LAC. Our results suggest CRM could be used to more expeditiously identify injured patients likely to deteriorate to shock, with requirements for blood transfusion or procedural LSI.

Convertino, Victor A.↗
DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Image Distinguishability Analysis Testing Through Principal Components and Its Application to Hot Spot Scale Invariance

Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.

organic↗
DOE OSTI2026

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗
DOE OSTI2023

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Invited article: X-ray phase contrast imaging in inertial confinement fusion and high energy density research

X-ray phase contrast imaging (XPCI) provides enhanced image contrast beyond absorption-based x-ray imaging alone due to refraction and diffraction from gradients in the object material density. It is sensitive to small variations in density, such as internal voids, cracks, grains, defects, and material flow, as well as to stronger density variations such as from a shock wave. Beyond its initial use in biology and materials science, XPCI is now routinely used in inertial confinement fusion (ICF) and high energy density (HED) research, first to characterize ICF capsules and targets, and later applied in dynamic experiments, where coherent x-ray sources, ultrafast x-ray pulses, and high temporal and spatial resolution are required. In this Review article, XPCI image formation theory is presented, its diverse use in ICF and HED research is discussed, the unique requirements for ultrafast XPCI imaging are given, as well as current challenges and issues in its use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

AT2025ulz and S250818k: Zooming in with the Hubble Space Telescope

AT2025ulz is an optical/near-infrared transient discovered during follow-up of the candidate gravitational wave (GW) event S250818k. Its young age (≲1 day), rapid decline, and strong color evolution over the first 48 hr classify it as a potential kilonova candidate. In this work, we present the results of our observing campaign, carried out with the Gran Telescopio Canarias and the Hubble Space Telescope (HST). Although the early time evolution of AT2025ulz resembles some aspects of a kilonova, its rapid onset (∼3 hr after the GW trigger) and luminosity (a factor of ∼5 brighter than AT2017gfo in g band) are difficult to reproduce. Only a small subset of our kilonova models matches its multicolor light curve, and the inferred ejecta mass is uncomfortably large given the low chirp mass (≲0.87 M ⊙ ) of the GW candidate. HST observations place the transient within a nearby (z = 0.08489) spiral galaxy with ongoing star formation and measure a color (F336W − F160W ≈ 1.4 mag) that is too blue to match with a kilonova. Our data support the classification of AT2025ulz as a supernova, initially undergoing a shock-cooling phase and later entering its photospheric phase and spectroscopically identified via its broad absorption features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Experimental measurement of energy release from an initiating layer in an insensitive explosive

When subjected to a shock of insufficient strength to trigger prompt reaction, heterogenous condensed phase explosives can form regions where significant amounts of the explosive remain unreacted for times much greater than the reaction time of the detonating explosive. This phenomena is observed for the explosive PBX 9502 (95 wt% TATB) both for planar and oblique input shocks. Here, we build on previous results by performing cylinder expansion (CYLEX) tests where the explosive charge is comprised of a faster core of PBX 9501 (95 wt% HMX) inside a slower annulus of PBX 9502. The detonation in the faster PBX 9501 drives an oblique shock into the adjacent PBX 9502, and an annular transverse initiating layer (IL) results. In the test geometry, the IL travels steadily down the length of the test after a short run distance. At radial positions beyond the IL, an annular region of detonating PBX 9502 is observed. Using standard CYLEX test diagnostics, we infer the total energy release of this experiment. By making the assumptions that (1) the combined energy release is comprised of contributions from detonating PBX 9501, detonating PBX 9502, and the IL in the PBX 9502 and (2) mass-specific energy release for the detonating explosives is approximately the same as typically observed for each explosive, the IL energy release and reaction efficiency can be computed. Results are compared to prior results for a similar geometry, and indicate that while shock deadened PBX 9502 does not detonate promptly, it does eventually release a significant portion of its chemical potential energy over longer timescales on the order of 10 µs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗
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

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗