Search NASA⌕ Search

SEARCH · Search NASA

Results for “detonation velocity”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

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

A Geometrical Feature of a Travelling Concave Detonation Front

Experiments show that when two HE materials placed adjacently, a fast-propagating detonation in one of the HE regions can sometimes initiate a detonation in the other HE. Direct numerical simulations support this phenomenon. Let A and B stand for the two HE materials. In a typical case, the driving detonation in HE A has a propagation speed faster than the CJ velocity in HE B, and the induced detonation in HE B is necessarily concave.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The effects of electrically exploding gold bridgewires into inert and explosive powder beds

Abstract The particle velocity created in beds of both low-density inert sugar and explosive PETN as a function of distance from an exploding bridgewire was measured using optical velocimetry and a silvered PMMA window. As expected, more violent bridge-bursts (from a greater-stored-energy capacitive discharge unit) resulted in greater particle velocities and a better supported compaction wave in sugar. In all cases, ramp waves, not shocks, were observed in the inert sugar. Large window velocities were observed for very powerful bursts (up to 270 m/s), but bursts required for stochastic detonator operation conditions resulted in sugar/PMMA window velocities of only 8–10 m/s 0.85 mm from the bridge location. In contrast, after a distance of only 0.65 mm, a building shock wave was observed in PETN under both threshold and reliable firing conditions. Subsequently a hot-spot-driven shock-to-detonation (SDT) process was observed prior to full detonation. The measured buildup process accounts for $$\approx $$ ≈ 66% of the so-called excess transit time (ETT) between the observed and theoretical total function time for the particular exploding-bridge-wire (EBW) detonator studied. The remainder must occur in the powerful output pellet region. In contrast to a common understanding, the ETT is found to be a weak function of the discharge energy. Thus, the operation of the detonator after a bridge-burst energy-to-powder reaction transition process is found to be hot-spot-driven SDT in both the low- and high-density pellets.

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

Near-Field Imaging of Shallow Chemical Explosions in Granite Using Change Detection Methods with Surface and Borehole Seismic Data

Explosions detonated in geologic media damage it in various ways via processes that include vaporization, fracturing, crushing of interstitial pores, etc. Seismic waves interact with the altered media in ways that could be important to the discrimination, characterization, and location of the explosions. As part of the Source Physics Experiment, we acquired multiple pre- and post-explosion near-field seismic datasets and analyzed changes to seismic P-wave velocity. Our results indicate that the first explosion detonated in an intact media can cause fracturing and, consequently, a decrease in P-wave velocity. After the first explosion, subsequent detonations in the pre-damaged media have limited discernible effects. We hypothesize this is due to the stress-relief provided by a now pre-existing network of fractures into which gasses produced by the explosion migrate. We also see an overall increase in velocity of the damaged region over time, either due to a slow healing process or closing of the fractures by subsequent explosions.

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

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
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 OSTI2025

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Effect of Heterogeneities on the Reaction-zone of a Propagating Detonation Wave

Reactive burn models for propagating detonation waves in an explosive are based on the ZND (Zel'dovich-Von Neumann-Doring) theory. A key property of the ZND theory, which leads to self-sustaining detonation waves, is a sonic point relative to the front in the reaction zone. Moreover, burn models assume the explosive is a homogeneous material. This assumption is used for plastic-bonded explosives, even though they are heterogeneous materials. The heterogeneities are accounted for only with an empirical burn rate. Furthermore, the reaction-zone width can be less than the length scale of heterogeneities. This raises the question of the effect of heterogeneities on the reaction zone. In particular, whether the reaction zone can be steady and have a well defined sonic point. To examine this question we focus on PBX 9501 since there is data on the grain-binder heterogeneities and on the reaction-zone velocity time history of a steady planar detonation wave. Using pseudo-mesoscale 1-D simulations, we examine the effect of a detonation wave transversing explosive and binder segments of random length. The simulations show that the reaction zone is quasi-steady with only short wavelength small amplitude variations in the detonation wave pressure and speed. The main effect of the heterogeneities is to introduce noise in the pressure and velocity fields that propagates into the Taylor wave behind the detonation. This is consistent with the scatter in the reaction-zone measurements for PBX 9501.

42 ENGINEERING↗
DOE OSTI2020

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

3D optical diagnostics for explosively driven deformation and fragmentation

High-speed, optical imaging diagnostics are presented for three-dimensional (3D) quantification of explosively driven metal fragmentation. At early times after detonation, Digital Image Correlation (DIC) provides non-contact measures of 3D case velocities, strains, and strain rates, while a proposed stereo imaging configuration quantifies in-flight fragment masses and velocities at later times. Experiments are performed using commercially obtained RP-80 detonators from Teledyne RISI, which are shown to create a reproducible fragment field at the benchtop scale. DIC measurements are compared with 3D simulations, which have been ‘leveled’ to match the spatial resolution of DIC. Overall, results demonstrate improved ability to identify predicted quantities-of-interest that fall outside of measurement uncertainty and shot-to-shot variability. Similarly, video measures of fragment trajectories and masses allow rapid experimental repetition and provide correlated fragment size-velocity measurements. Measured and simulated fragment mass distributions are shown to agree within confidence bounds, while some statistically meaningful differences are observed between the measured and predicted conditionally averaged fragment velocities. Together these techniques demonstrate new opportunities to improve future model validation.

42 ENGINEERING↗
DOE OSTI2020

Scaling Law for Cloud-Rise Velocity vs. Scaled Height of Burst (SHOB)

The shockwave produced by a nuclear detonation leaves a low-density region in the vicinity of air zero (i.e., the point of detonation). As the shockwave degenerates into an acoustic wave, the buoyant force takes over and causes the low-density region to begin to rise at speeds that can reach as high as ~200-300 mph. The objective of this project was to find out if this initial cloud rise followed any scaling law, by stimulating a cloud rise for a 1 kiloton (kt) detonation and determining if a correlation existed between the scaled-height-of-burst (SHOB) and the scaled-cloud-rise velocity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Laser Wavelength Shift from Splitters and Switches

The Detonation Science and Technology Group (Q-6) operates several photon Doppler velocimetry (PDV) systems as a diagnostic tool for measuring velocities in the 1–10 km/s range during high explosive (HE) detonation experiments. Q-6 uses “heterodyne” PDV systems that are produced by Traverse Solutions and Instrumentation, LLC (TS&I). These heterodyne systems combine wavelength (WL) tunable source and “reference” lasers with a WL offset of 10–100 pm from each other. Properly setting the source and reference WLs for an experiment can be critical for collecting the desired experimental data. The tunable lasers report WL readings in their control software, but these readings are based on a lookup table from a temperature reading and are not able to be calibrated. Therefore, recent test fire programs that require PDV as a diagnostic have also implemented requirements for monitoring the WL of the tunable lasers using a precision WL meter. To monitor the WL of several PDV lasers without disconnecting any patch fibers form the experimental setup, Q-6 developed a system splitters with a coupling ratio (%) of 99:1 and an optical switch with the switch output going to a single WL meter. A brief experiment was developed to determine if the splitters and/or switch introduced any WL bias into the readings. The purpose of this report is to document the system, the experiment, and its results.

42 ENGINEERING↗
DOE OSTI2021

Mass flow rate effect on a rotating detonation combustor with an axial air injection

Changing the flow rate of reactants being injected into a rotating detonation combustor (RDC) results in interesting behavior of the system. Prior studies have found that an increase in mass flow rate gradually increases the detonation wave speed before splitting the wave into multiple fronts. Here, the focus of this study is in understanding the physics of such behavior through a combination of experiments and numerical simulations. For this purpose, the axial air inlet-based RDC geometry is used. In the experiments, the wave velocity increased similar to prior studies. Corresponding full-scale simulations show that increase in mass flow rate by increasing pressure of the feed plenums shortens the recovery time of the injectors. This causes a more uniform fuel–air mixture to form prior to the arrival of the detonation wave. As a result, a more ideal detonation is observed, which leads to an increased wave velocity. Details of the detonation structure and variation with mass flow rate are analyzed. The presence and variation of the different deflagration zones in such non-premixed RDCs are also discussed.

42 ENGINEERING↗
DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

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

Design and validation of a nanosecond short-wave infrared streak camera using gaseous detonation experiments

A custom short-wave infrared (SWIR) galvanometric streak camera was developed to provide nanosecond- scale, time-resolved optical diagnostics. The system was designed as a compact, mechanically driven alternative to traditional tube-based streak cameras, using a galvanometer mirror and multiple-reflection optical path to achieve high temporal resolution while maintaining electronic synchronization for trigger sequencing. Temporal resolution achieved 139.5 ns/px, allowing direct conversion of image-slope gradients to physical wave velocities. The instrument was validated using hydrogen-oxygen-argon detonations within a modular detonation tube, capturing transient emission structures and post-detonation decay behavior with microsecond precision. Although signal levels in filtered configurations were limited by lens transmission and detector sensitivity, the system successfully resolved intensity rise and decay characteristics across 10-60 % argon mixtures. The observed trends demonstrate the feasibility of compact galvanometer-based streak systems for high-speed SWIR imaging in reactive environments!

42 ENGINEERING↗