Modeling and observation of detonation front formation with microwave interferometry during the deflagration-to-detonation transition
Explore the source record for details and available documents.
Chemistry - Chemical explosives↗
SEARCH · Engineering Papers
Results for “Detonations”
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.
Chemical timescale effects on detonation convergence
Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.
Adaptive mesh refinement↗
Exploding Bridgewire (EBW) Detonators: An Example of Synergistic Multiphysics
Exploding bridgewire (EBW) detonators are highly temporally reproducible explosive devices that require the rapid discharge of a high‐voltage capacitance to operate and so are immune to most of the accidental hazards associated with traditional electric detonators. They have been demonstrated to be safe enough for use in high‐consequence explosive applications. Despite continued use for over 82 years, understanding the exact mechanism of operation has remained elusive. Various researchers have ascribed either deflagration‐to‐detonation (DDT) or shock‐to‐detonation (SDT) phenomena observed in other explosive events to explain the science behind the successful engineering; however, a rigorous justification has been absent. Previously, we have demonstrated a complex interaction in EBW detonators between large electrical currents, non‐equilibrium thermodynamic material states, plasma physics, powder compaction phenomena, shock physics, photochemistry, and rapid conventional explosive chemical reaction processes. Specifically, we have made progress in understanding the complex multiphysics that operates in these detonators and demonstrating that it is a serendipitous synergy between UV light emitted from the arc formed as the bridge is electrically exploded and the accompanying short‐duration shock transmitted into the explosive powder bed that allows these devices to function at practical capacitor sizes and charge voltages. This insight not only places the topic on a firmer scientific footing but potentially enables new approaches to safe detonator design.
36 MATERIALS SCIENCE↗
Additive manufacturing of high explosives with inert dilution for wave shaping and detonation velocity grading
The performance of a high explosive charge is highly influenced by the detonation velocity and any inherent internal structure. Additively manufactured (AM) high explosive charges allow for a precisely engineered internal structure and the selective placement of discrete or graded volumes of a tailored material, which are unachievable by conventional means. The dilution of a solid explosive by an inert additive for detonation velocity reduction is herein evaluated for a method of material extrusion called a direct-ink-write AM, where an 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane-based ultraviolet -curable energetic-binder paste explosive was diluted by the inert density mock 5-iodo-2′-deoxyuridine for use in the manufacture of multi-material explosive charges. The explosive ink was characterized for printing by rheological and particle size measurements, and the detonation velocities for 0–40 wt. % inert diluent were found to have a continuous linear reduction in detonation velocity down to 14.3% difference. A two-component 3D plane wave explosive lens with 8.038 km/s “fast” and 7.491 km/s “slow” diluted ink in an internal cone demonstrated a reduced time of breakout at the output face from 361 to 53 ns. A linearly graded “fast–slow–fast” cylinder was successfully printed and fired, showing engineered control along the length of the charge as the detonation velocity dropped from 7.851 to 6.700 km/s and then recovered to 7.241 km/s. The extreme control over the internal composition of a high explosive charge via a dual-material extruder/mixer challenges conventional manufacturing approaches and opens new avenues for engineered detonation wavefronts, high explosive charge performance, and new applications.
36 MATERIALS SCIENCE↗
Mathematical modelling of the concave front in the adjacent high explosive detonation problem
This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.
acceleration↗
Single-particle chemical analysis reveals organic-rich detonation soot products
The detonation of high explosives produces a wide variety of particulate matter (PM) with distinct properties, not all of which are traditionally studied for chemical composition, formation processes, and forensic applications. We report particle-resolved measurements of Composition B detonation soot using soot particle aerosol mass spectrometry (SP-AMS), identifying carbonaceous species and metals not previously characterized on a single-particle basis. Results are combined with scanning electron microscopy with energy dispersive X-ray spectroscopy (SEM-EDS) to enhance source-dependent signatures. Black carbon, including graphitic carbon and detonation nanodiamonds, contributed 50.5-71.4% of PM 2.5 mass, while 22.5-43.4% was non-refractory organic carbon, a previously overlooked component that exhibited a complex and varying composition. Detonations were performed with and without PMMA confinement and under steady and overdriven conditions. Distinct particles enriched with polycyclic aromatic hydrocarbons (PAHs) were observed in experiments utilizing confinement, with quantities dependent on manufacturing method. SEM-EDS validated SP-AMS findings of metals internally mixed with carbonaceous species and extended the particle size range to 100 µm. This work makes detecting non-refractory organics using single-particle techniques more feasible for detonation forensics and understanding high-energy soot formation. While this analysis uses offline sample collection, SP-AMS could be deployed for in-situ measurements of detonation PM 2.5 transported in the atmosphere.
36 MATERIALS SCIENCE↗
Equation of state measurement of detonation carbon condensates using optical microscopy and interferometry
Thermochemical models of detonation that estimate performance (e.g., detonation velocity, energy delivery, etc.), are based on assumptions that carbon condensates (soot) formed during detonation is largely similar to bulk carbon. However, soot constituents can range from amorphous carbon to nanodiamond and include other material phases. Since thermodynamic properties of the soot such as compressibility are imperative for accurate thermochemical modeling of detonation reaction chemistry, experimental measurements of the equation of state (EOS) which determine the compressibility are vital. Due to the mixed-phase nature of detonation soot, typical methods to measure the EOS (e.g., x-ray diffraction) are untenable. In this study, the high-pressure EOS up to 20 GPa was determined for detonation soot collected from PBX 9502, Composition B (Comp B), Hexanitrostilbene (HNS), and LX-21 high explosives by employing a direct volume technique using optical microscopy and interferometry in a diamond anvil cell. Comp B soot was determined to be the least compressible [K 0 = 57.9(17) GPa] with HNS soot [K 0 = 53.7(15) GPa], LX-21 soot [K 0 = 45.8(51) GPa], and PBX 9502 soot [K 0 = 28.2(27) GPa] being more compressible, likely due to differences in nanodiamond content as compared to amorphous carbon and graphite content.
Amorphous materials↗
Hypersonic Jets of Detonation Products in the Hydrodynamic Collapse of Macroscopic Voids
Localizing the energetic output from detonation waves has been a long-standing challenge in applied detonation physics. Here, energy localization is achieved via machined millimeter scale voids in pressed samples of PBX 9501, an HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane)-based plastic bonded explosive. A main mechanism of energy localization in these systems, the formation of hydrodynamic jets of dense product gases, is characterized experimentally using a semicylindrical geometry in witness plate impact experiments and streak imaging of the jet propagating into the air. The distance at which the jet is optimally developed is identified and the supersonic flow structure in the vicinity of this feature is explored using hydrocode simulations. This analysis found that most of the kinetic energy of the hydrodynamic jet arises from pressure gradients induced by geometrically mediated squeeze flow lateral to the direction of detonation propagation. This work presents a new development in the control of energetic output from detonation waves and applications to detonation wave shaping are discussed.
42 ENGINEERING↗
Detonation Waves in High Explosives
A material at high temperature can react or decompose. For an energetic material, the reaction is exothermic and releases chemical energy that would further increase the temperature. Under some circumstances, when a reaction is triggered, such a reaction can propagate and the material rapidly releases a large amount of energy giving rise to an explosion. Examples of such materials are aerosols, suspensions of solid particles or liquid droplets in a gas; such as coal dust, grain dust and fuel-air explosions. Frequently, explosions are due to accidents. A spectacularly destructive example is the recent explosion of a large quantity of ammonium nitrate (thousands of tons) in Beirut, Lebanon (August 2020); see for example Beirut explosion. Ammonium nitrate is used as a fertilizer. It and the aerosols are not considered to be explosives due to the limited conditions for which an explosion can occur. An aerosol gets the oxidizer from the surrounding air. Burning requires diffusion of the oxidizer to the particle surface where the reaction occurs. A large density of small particles is required for a fast enough reaction to support an explosion. In contrast, an explosive is an energetic material with both fuel and oxidizer mixed on a molecular scale (either premixed gases or within molecules of a solid). This allows fast enough reactions over a wide range of conditions to support a self-propagating reactive wave known as a detonation wave. A detonation wave can be controlled and an explosive used for useful purposes such as in mining, construction, demolition, explosive welding, argon flash lamp, pulsed power using a magnetic flux generator [see also Goforth et al., 2015], jet cutter with shaped charge, explosive art, and generating conditions to study the response of materials at high strain rates and high pressures [see for example, Marsh, 1980]. Explosives are also used in conventional munitions and nuclear weapons. The focus of this book is on the theory and phenomenology of solid high explosives (HEs); in particular, plastic-bonded explosives (PBXs). Some aspects of detonation wave theory are needed to interpret explosive data. Hence, the theory is presented before the detonation wave phenomenology. A familiarity with fluid flow, specifically the notion of shock waves and the shock loci are assumed. In the remainder of this chapter we give a brief overview on the basic properties of detonation waves and PBXs.
36 MATERIALS SCIENCE↗
Effect of a collapsing gas bubble on the shock-to-detonation transition in liquid nitromethane
We studied the shock-induced collapse of butane gas bubbles in the homogeneous explosive nitromethane (NM) to investigate the effects of hot spot formation on the detonation process. A butane bubble was injected into a sample of NM, and a shock wave from a flat plate impactor compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and bubble collapse processes. A multiband optical fiber pyrometer measured the time-resolved thermal radiance, and we used the results and emissivity values extracted from spectral fits to estimate temperatures. We measured the characteristics of the shock-to-detonation transition in NM with and without a bubble. All experiments were performed at shock pressures near 8 GPa, where neat NM can detonate. A single bubble in this system was shown to sensitize NM, leading to a reduced run-to-detonation time. We used hydrodynamic modeling to predict shock wave propagation, the extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. We used a temperature-dependent Arrhenius burn model for simulations, and it yielded much better results than reactive burn models that depend only on pressure and density.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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↗
Concept of Dynamic Heat Insulation for Rotating Detonation Engines
This work introduces a new class of materials concept to dynamically reduce instantaneous heat fluxes in Rotating Detonation Engine (RDE) combustor chamber walls. The high-frequency and high-amplitude surface heat fluxes observed in RDEs arise from large instantaneous temperature differences between the detonation shockwave and chamber wall surface. These temperature gradients drive substantial energy loss and reduce the chamber gas pressure, ultimately limiting the cycle’s thermodynamic efficiency. This work introduces a concept for dynamically insulating the combustion chamber surfaces using surface layers or coatings with low thermal time scale. With such coatings, the surface temperature may follow the fluctuations of the cyclic detonation wave temperature, thus reducing the instantaneous heat flux therefore cycle-mean heat flux. To analyze these cyclic thermal phenomena, a one-dimensional analytical conduction solver was utilized with the capability to handle multilayered structures. Parametric modeling was performed using transient heat flux boundary conditions representative of a hydrogen–air RDE across a broad range of coating thermal properties and engine conditions. The coating effectiveness scaled with the product of thermal time constant and detonation wave frequency and the results were non-dimensionalized to guide future materials development. This strategy may offer benefits in increasing material survivability, reducing cooling requirements, and enhancing pressure gain.
heat transfer↗
Laser absorption measurements of temperature, pressure, CO, and CO 2 at near-MHz rates in post-detonation fireballs with comparison to synthetic measurements
A laser absorption spectroscopy (LAS) diagnostic was used to obtain measurements of temperature, pressure, CO, and CO 2 at 500 kHz or 1 MHz in post-detonation fireballs produced by hemispherical samples of pentaerythritol tetranitrate (PETN). A quantum-cascade laser was scanned over multiple CO absorption transitions near 2008.5 cm −1 at 1 MHz, while an interband-cascade laser was scanned over a CO 2 absorption transition near 2394.8 cm −1 at 500 kHz. Light from each laser was combined onto a single path and passed through a detonation chamber approximately 83 mm above the 12-mm diameter hemispherical PETN charge. The CO and CO 2 absorption signals were post-processed to obtain time histories of temperature, pressure, species column pressures (P CO L, P CO2 L), and species column mole fractions (X CO L, X CO2 L). Additionally, schlieren imaging was performed simultaneously at 500 kHz to aid interpretation of the LAS measurements. Experimental and synthetic (i.e., CFD based) LAS measurements were compared to evaluate the accuracy of the CFD model and its ability to model the turbulent afterburning of the detonation products in air. In general, the experimental measurements exhibit reasonable agreement with the synthetic measurements at early times; thereby supporting the accuracy of the CFD model. Periods of disagreement between experimental and synthetic measurements at later times are most likely due to a reflected shock and detonator cavity jetting, which are not accounted for in the CFD model.
Schwartz, Charles J. [Purdue Univ., West Lafayette↗
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↗
Deflagration to Detonation Transition Update: XDDT Code Modularization
A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Stability and convergence of nuclear detonations in white dwarf collisions
We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.
79 ASTRONOMY AND ASTROPHYSICS↗
Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors
This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.
gas turbines↗
Z Machine Detonator Holder Testing
Current packaging for the detonators used on the Z Machine is dirty, not well suited for reusability, and labor intensive to use. A new packaging design is developed, explosively tested, and approved by the numerous stakeholders resulting in more efficient operations and lowered costs.
42 ENGINEERING↗