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

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Novel method to control explosive shock sensitivity: Formulations, experiments and modeling

The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Re-calibration of PBX9501 SURF model

PBX 9501 is a plastic bonded explosive composed of 95 wt % HMX and a binder; see [Gibbs and Popolato, 1980, pp. 109–119]. SURF is a reactive burn model for shock initiation and propagation of detonation waves. It has previously been calibrated for PBX 9501. Here the SURF model is recalibrated for PBX 9501; specifically, lot 730-010 at ρ = 1.837 g/cc 3 . The new calibration uses the Davis reactants and products EOS calibrated for the AWSD model [Aslam et al., 2020]. The burn rate in the shock initiation regime is fit to the Pop plot from 5 embedded gauge shock-to-detonation transition (SDT) experiments from [Gustavsen et al., 1999, see fig 12 and table 5]. In the propagation regime, the burn rate is fit to curvature effect data (detonation speed as function of front curvature); see [Aslam, 2007]. Also the burn parameters are adjusted to fit the gap-stick experiment [Hill et al., 2018]. Simulating the detonation wave speed in this experiment requires a model that is accurate for initiation with complex shock loading; in particular, a pressure decreasing gradient behind a curved lead shock. This is more difficult than calibrating to the standard SDT experiments which are 1-D and driven by a sustained shock. Simulations of the gap-stick experiment for PBX 9501 with the SURF model will be discussed in a subsequent report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Structural Transformation and Chemical Stability of a Shock-Compressed Insensitive High Explosive Single Crystal: Time-Resolved Raman Spectroscopy

Despite the considerable interest in insensitive high explosives (IHE) as a safer alternative to conventional high explosives, a good understanding of the low sensitivity of IHEs to shock initiation is lacking. In particular, real-time measurements to directly probe the molecular-level response of shock compressed IHE single crystals constitute an important need. To address this need, plate impact experiments were conducted to determine time-resolved changes in the Raman spectra of 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shock compressed up to 20 GPa longitudinal stress. The Raman measurements examined vibrational frequencies from 800 cm-1 to 1500 cm-1 with 15 nanosecond time resolution and were conducted at several peak stresses. At 4 – 6 GPa, two new Raman peaks appeared, in addition to the original peaks, consistent with onset of the α' – ε structural transformation reported previously in static compression work. The measured spectra indicated completion of the transformation at 10 GPa. Raman data to 20 GPa showed neither additional transformations nor any indication of chemical decomposition. This finding, though consistent with recent continuum measurements, is in marked contrast to the chemical decomposition observed at lower stresses in shock compressed conventional high explosive single crystals. Our Raman results support the previous suggestion that strengthening of intra- and inter-molecular bonds, due to the α' – ε structural transformation, plays a significant role in the insensitivity of FOX-7 single crystals to shock initiation. In conclusion, the present work, in conjunction with previous static compression studies, provides the first experimental insight into the molecular-level response of a shock compressed IHE single crystal and can serve as a bench mark for theoretical studies.

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

A Data Set for Intercomparing the Transient Behavior of Dynamical Model-Based Subseasonal to Decadal Climate Predictions

Climate predictions using coupled models in different time scales, from intraseasonal to decadal, are usually affected by initial shocks, drifts, and biases, which reduce the prediction skill. These arise from inconsistencies between different components of the coupled models and from the tendency of the model state to evolve from the prescribed initial conditions toward its own climatology over the course of the prediction. Aiming to provide tools and further insight into the mechanisms responsible for initial shocks, drifts, and biases, this paper presents a novel data set developed within the Long Range Forecast Transient Intercomparison Project, LRFTIP. This data set has been constructed by averaging hindcasts over available prediction years and ensemble members to form a hindcast climatology, that is a function of spatial variables and lead time, and thus results in a useful tool for characterizing and assessing the evolution of errors as well as the physical mechanisms responsible for them. A discussion on such errors at the different time scales is provided along with plausible ways forward in the field of climate predictions

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2026

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Role of Sea‐Surface Salinity in Simulating Historical Decadal Variations of Atlantic Meridional Overturning Circulation in a Coupled Climate Model

Abstract Extending climate prediction from seasonal to decadal timescales requires realistic initialization of not only upper ocean heat content but also the Atlantic meridional overturning circulation (AMOC). However, it remains a major challenge to realistically initialize AMOC in a coupled system while also maintaining a balanced atmosphere‐ocean initial state. This study demonstrates the feasibility of generating fully coupled historical states with realistic AMOC variability. Employing a forced ocean—sea‐ice (FOSI) model simulation as the “truth,” we show reproducibility of key features of historical AMOC decadal variability in a fully coupled model by restoring sea‐surface salinity, in addition to sea‐surface temperature restoring widely used in seasonal prediction. The atmospheric state of the restored coupled model solution is much closer to that of the free coupled simulation than to the observations used in FOSI, pointing to potential advantages of using this approach for initializing decadal predictions with reduced ocean initialization shock.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Imaging perturbed shock propagation in powders

A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Cascading economic losses from port disruptions under capacity constrained multimodal freight networks

This study quantifies how throughput disruptions at major seaports cascade through capacity-constrained multimodal freight networks and interregional production systems. We couple an agent-based model (ABM) multimodal freight simulation that resolves rerouting, terminal queueing, and inventory drawdown under binding modal and facility capacities with a multiregional output loss input-output (MRIIM) model that propagates realized delivery shortfalls across regions and sectors. The framework is demonstrated for the Port of Los Angeles using Freight Analysis Framework flows and Bureau of Economic Analysis input-output accounts and is evaluated over a 52-week horizon under deterministic sector targeted shocks and stochastic disruption realizations with uncertain severity and duration. Results indicate nonlinear amplification: realized national losses concentrate in manufacturing and transportation/warehousing even when exogenous port shocks are dispersed, suggesting that congestion spillback and limited short-run substitution can dominate the initial shock allocation. We further evaluate a tabular reinforcement-learning (Q-learning) intervention layer that selects among a small set of implementable system level levers (truck-to-rail and truck-to-barge shift settings) without overriding shipper routing, finding that such interventions reduce total losses for moderate disruptions but yield diminishing returns once substitute modes approach capacity. By linking operational freight behavior to system wide impacts under uncertainty, the proposed ABM-MRIIM pipeline provides a reusable workflow for port disruption stress testing, identification of structurally critical sectors/corridors, and evaluation of resilience interventions under realistic capacity limits.

42 ENGINEERING↗
DOE OSTI2021

Near-optimal combination of high performance and insensitivity in a shock compressed high explosive single crystal

Achieving the desired combination of superior detonation performance and insensitivity to shock initiation has been a long-standing goal in high explosive (HE) science and technology. Having previously established the shock insensitivity of 1,1-diamino-2,2-dinitroethene (also known as DADNE or FOX-7) single crystals to 20 GPa (extended to 25 GPa in this work), the FOX-7 detonation response was determined through wave profile measurements in ~250 μm thick single crystals shock compressed to 64 GPa. Quite unexpectedly, FOX-7 demonstrated the classic Chapman–Jouguet (C–J) detonation response—reaction completion in the detonation front (<0.7 ns) at pressures of 44 GPa and higher—not observed in other insensitive high explosives. Here, the experimentally determined C–J pressure (35 GPa), detonation wave velocities and the detonation products equation of state—together with shock insensitivity to 25 GPa—demonstrate that FOX-7 single crystals display a near-optimal combination of high performance and shock insensitivity, not observed in another HE crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

SURFplus model calibration for PBX 9012

PBX 9012 is a plastic bonded explosive composed of 90.2 wt % HMX (cyclo-tetramethylene-tetranitramine) and 9.8 % Viton A (vinylidene-hexafluoropropylene copolymer). It has nearly the same composition as LX-07 (90% HMX and 10% Viton A) but with a slightly lower density. Hence it has a higher porosity and is slightly more shock sensitive. Here we describe the SURFplus model calibration for PBX 9012 using data from five shock-to-detonation transition (SDT) gas gun experiments by Burns and Chiquete, detonation front shape data from rate stick experiments and cylinder test experiment reported in Anderson et al., Jackson et al. The SDT experiments provide Pop plot data points (run distance-to-detonation) which characterize shock initiation, Lagrangian velocity time histories from embedded magnetic velocity gauges which give information on the reaction behind the lead shock, and Hugoniot data points for the reactants EOS. The rate stick experiments provide data on the curvature effect which characterize propagating detonation waves. The cylinder test experiments provide data on the detonation release isentrope used to calibrate the products EOS. The HE model uses a Davis reactants EOS with parameters that Burns and Chiquete [2020] fit to Hugoniot data, and Davis products EOS with parameters that Anderson et al. [2020] fit to cylinder test data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Femtosecond Reflectance Spectroscopy for Energetic Material Diagnostics

Understanding the fundamental mechanisms underpinning shock initiation is critical to predicting energetic material (EM) safety and performance. Currently, the timescales and pathways by which shock-excited lattice modes transfer energy into specific chemical bonds remains an open question. Towards understanding these mechanisms, our group has previously measured the vibrational energy transfer (VET) pathways in several energetic thin films using broadband, femtosecond transient absorption spectroscopy. However, new technologies are needed to move beyond these thin film surrogates and measure broadband VET pathways in realistic EM morphologies. Herein, we describe a new broadband, femtosecond, attenuated total reflectance spectroscopy apparatus. Performance of the system is benchmarked against published data and the first VET results from a pressed EM pellet are presented. This technology enables fundamental studies of VET dynamics across sample configurations and environments (pressure, temperature, etc .) and supports the potential use of VET studies in the non-destructive surveillance of EM components.

42 ENGINEERING↗
DOE OSTI2021

Equation of State Measurements on Iron Near the Melting Curve at Planetary Core Conditions by Shock and Ramp Compressions

The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. In this work, we use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.

58 GEOSCIENCES↗
DOE OSTI2021

Shock-reshock of zirconium to explore the effect of microstructure on the alpha-omega phase transformation

Phase transformations play an important role in the mechanical behavior of materials subjected to extreme loading conditions. A series of shock-reshock experiments were fielded to determine whether the phase transitions in materials are significantly enhanced or inhibited by preexisting microstructural features. Polycrystalline zirconium samples were shock loaded using gas-gun plate impact and soft recovered to examine microstructure using electron backscatter diffraction (EBSD). Drive conditions were varied to study the (hcp) alpha to (hexagonal) omega solidsolid phase transformation. Recovered samples were then subjected to a second shock loading event to determine the change in material behavior as a function of pre-shock microstructure. Crystallography of phase fragments in the final microstructure showed that prior twinning (formed during the shock to a peak stress below the transition threshold) appeared to suppress omega formation/retention after reshock. Conversely, when a material was initially shocked into the omega phase field, retained-omega was not found to have a large impact on subsequent omega formation during reshock. This suggests that nucleation and growth of omega phase are important processes, and the relative activity of nucleation vs. growth processes is modified by a pre-existing substructure. Additionally, orientation relationships reveal a reverse transformation pathway (omega to alpha) dominates the final microstructure, suggesting significant grain growth in the omega phase field is possible even for dynamic timescales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

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

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗