SEARCH · Engineering Papers

Results for “detonation 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Shock initiation of low density polymer bonded explosive LX-14: A study of two morphologies

A series of six shock initiation experiments have been carried out on low density LX-14 powder in order to simulate a shock insult on the heavily damage polymer bonded explosive, LX-14. Two distinct morphologies were studied (tap density molding powder and machined swarf), both at the same density of 0.942 g/cm 3 , or 50.1% theoretical maximum density. The purpose of these experiments was to provide data to help make an assessment of the effects that damage has on the material sensitivity to a planar shock. This was achieved primarily by providing shock sensitivity data in the form of a Pop plot, and also reactants equation of state data, which aids in the determination of input conditions for both the experiments performed in this work and also future experiments of this material type. The experiments were of a cut-back format, consisting of four sample heights on each shot and diagnosed with optical velocimetry. The experiments were carried out at the Technical Area 40 Chamber 9 gas gun facility at the Los Alamos National Laboratory, where the LX-14 targets were subjected to Al 6061 and Oxygen Free High Conductivity copper impactors launched to velocities up to 2.14 km/s. Time corrected reactive growth wave profiles are presented in this paper along with the derivation of the following Hugoniot parameters for this explosive, where the molding powder and machine swarf are represented by the following linear equations, respectively, U s = 1.64 (±0.64)u p + 2.25 (±0.49) and U s = 0.87 (±.0.60)u p + 4.08 (±0.49). The results show that the steady increase in shock sensitivity with increasing void fraction reaches an inflection point beyond which the shock sensitivity begins to decease. This inflection point lies between 65% and 50% of the theoretical maximum density of the LX-14. In both cases, the damaged LX-14 was not as sensitive as expected, with a relative shock sensitivity of the molding powder being less than the pressed LX-14, and the machined swarf having a shock sensitivity that is comparable to pressed LX-14.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Opto-thermal laser detonator

An opto-thermal laser detonator uses resonantly absorptive tuned nano-material associated with secondary explosives for optical absorption and initiation by an integral laser diode. The opto-thermal laser detonator includes main explosive material; resonantly absorptive tuned nano-material; secondary explosive material, wherein the resonantly absorptive tuned nano-material and the secondary explosive material are associated to form associated material made of the resonantly absorptive tuned nano-material and the secondary explosive material; and a laser diode operatively connected to the associated material, wherein the laser diode initiates the associated material which in turn initiates the main explosive material.

Wilkins, Paul R.↗

DOE OSTI2023

Opto-thermal laser detonator

Wilkins, Paul R.↗

DOE OSTI2023

Effect of thermal conditioning on the initiation threshold of secondary high‐explosives

While most performance metrics of high–explosive (HE) based devices like detonation velocity, detonation pressure, and energy output are expected to degrade over time, the evolution of initiation threshold appears less clear, with claims of both increasing and decreasing trends in threshold having been made in the literature. This work analyzes D–optimally designed sequential binary test data for a few thermally conditioned porous–powder and polymer–bonded HE initiator systems using a Bayesian likelihood method employing the probit regression model. Here we find that in most cases the initiation threshold decreases (i.e., sensitivity increases) upon accelerated thermal conditioning. However, such results are nuanced and influenced by factors like the contact area of initiating stimulus, HE characteristics like density and specific surface area, as well as possible thermally induced changes to other materials and interfaces involved.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. In this work, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N–N and axial NO2 stretching modes with the long-lived, excited phonon bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

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↗

DOE OSTI2023

Technical Report on Waveform Fit Metrics for Global Models

The new WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes an increased emphasis on the development of Earth models and methods to predict entire seismic and acoustic waveforms more accurately. In general, this increased emphasis is predicated on the need to better characterize seismic events and provide improved model-based discrimination between event types including earthquakes and explosions. More specifically, while current moment tensor inversion methods tend to work well for larger events (M>~4) using tuned 1-D Earth models, the development of state-of-the-art 3-D models and methods is required for the prediction of shorter period waves over large areas for discrimination of smaller events. There is no standard metric for model-based waveform prediction accuracy used in the waveform modeling/inversion community. However, there are several popular waveform misfit definitions; and minimizing the corresponding objective functions is the goal of waveform inversion. Some example misfit definitions employed for adjoint waveform tomography include measures of simple travel time differences (e.g. Tape et al., 2010), cross-correlation travel time differences (e.g. Luo and Schuster, 1991), multi-taper frequency dependent methods (e.g. Lei et al., 2020), time-frequency phase misfit functions (e.g. Fichtner 2010; Rodgers et al., 2022), normalized cross-correlation methods (e.g. Tao et al., 2018), and others. In some cases, these misfit definitions also involve complicated weighting schemes and summations over multiple frequency bands making it difficult to duplicate the misfit measurement with alternative models and datasets. Although each of the misfit definitions mentioned above are useful for developing waveform models, the actual misfit values are not usually meaningful outside of a given project, model, and/or dataset. Therefore, it is difficult to understand and communicate model performance for predicting waveforms and comparing to other models and/or new model iterations with a different dataset. Therefore, there is a need for a generalized method for evaluating overall model performance that is independent from the specific misfit chosen to develop the waveform models that is also intuitive and meaningful. In this report, we describe a new metric we refer to as ‘Percent of Correlated Signal’. The following sections describe and demonstrate the metric with a case study event and a more rigorous test using a random selection of globally distributed events. While the focus here is on global tomography models, the metric is meant to applicable to regional ‘wiggle-for-wiggle’ waveform models/studies as well.

58 GEOSCIENCES↗

DOE OSTI2024

Report on the sensitivity kernel construction and updating the SPiRaL model with waveform data

The WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes research leading towards the prediction of entire seismic and acoustic waveforms produced by natural and manmade events, including explosions. A key aspect of this research is the development of Earth (seismic) models at multiple scales including crustal, regional, and global. The work described here pertains to the global-scale seismic tomography effort led by LLNL.

58 GEOSCIENCES↗

DOE OSTI2021

PBX9012 Tracker Gauge Data

Five shock-to-detonation transition (SDT) gas gun experiments for PBX 9012 were performed by Burns and Chiquete [2020]. The data from the experiments is available on the small-scale database. Here we focus on the tracker gauge data for the lead shock trajectory, x(t). In order to determine the shock speed U s (t) and the detonation transition point, a least squared fit to the shock trajectory is used based on a pair of ODEs with a fitting form for the shock acceleration developed by Hill and Gustavsen [2002].

42 ENGINEERING↗

DOE OSTI2026

User’s Manual for RESRAD-RDD&IND Code Version 2: Vol. 2—User’s Guide for RESRAD-RDD&IND Code

Version 2.0 of the RESRAD-RDD&IND computer code is designed to support the implementation of protective action guides (PAGs) after a nuclear emergency incident including a radiological dispersal device (RDD) and/or an improvised nuclear device (IND) incident (EPA 2017). Eight different group types, addressing various decisions, are available for selection. The RESRAD-RDD&IND code calculates radiological doses, stay times, etc., for the selected group that the user wishes to focus on. (That is, the results for all the groups are not calculated simultaneously, and the input for those other groups do not matter, although some parameter values are shared between groups.) Version 2.0 has a user-friendly interface so that the RESRAD-RDD&IND code can be used with minimal training. For example, the user can select the major characteristics of the problem-event type, source term, and decision type from the left side of the interface and then calculate the results with the default assumptions for the exposure scenarios. More in-depth analysis would include specifying site-specific exposure scenario characteristics in the right side of the interface. The procedures for data entry and results viewing are self-explanatory. This is because common window maneuvering features and text instructions were incorporated in the interface design. General and context-specific help are available to aid users entering parameter values, as well. The RESRAD-RDD&IND computer code gives the user the option to select either an RDD or IND incident for analysis. For an RDD event analysis, 11 radionuclides (Am-241, Cf-252, Cm-244, Co-60, Cs-137, Ir-192, Po-210, Pu-238, Pu-239, Ra-226, and Sr-90) are included. These 11 radionuclides are the radionuclides most likely used for an RDD. More than 90 radionuclides can be selected for an IND event analysis. Initial default concentrations are provided for 44 radionuclides for a uranium-fueled IND event. These 44 radionuclides are those that would contribute significantly to the radiation dose associated with a uranium-fueled bomb detonation. The radionuclides generated from ingrowth of these 44 initial radionuclides are also automatically included in the analysis. Pu-239, Cs-134m, Ru-105, and Rb-89 and their progeny can be selected for analysis if they are detected and their concentrations are determined. This user’s guide, which is Volume 2 of the User’s Manual for RESRAD-RDD&IND Code Version 2, provides instructions to users on how to install the RESRAD-RDD&IND code, navigate the interface, and use the various features, including those discussed above, to set up an analysis and view/print the results in text outputs. Volume 1 of the User’s Manual for RESRAD-RDD&IND Code Version 2 (Yu et al. 2026), which contains descriptions of the methodology and theoretical basis for dose modeling and the mathematical equations implemented in the code, can be accessed and viewed through the Help menu in the code or can be downloaded from the RESRAD website (https://resrad.evs.anl.gov).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Single-shot, double-pulse determination of the detonation energy in nanosecond-laser ablation using the blast model

We demonstrate a novel single-shot method to determine the detonation energy of laser-induced plasma and investigate its performance. This approach can be used in cases where there are significant shot-to-shot variations in ablation conditions, such as laser fluctuations, target inhomogeneity, or multiple filamentation with ultrashort pulses. The Sedov blast model is used to fit two time-delayed shadowgrams measured with a double-pulse laser. We find that the reconstruction of detonation parameters is insensitive to the choice of interpulse delay in double-pulse shadowgraphy. In contrast, the initial assumption of expansion dimensionality has a large impact on the reconstructed detonation energy. The method allows for a reduction in the uncertainties of blast wave energy measurements as a diagnostic technique employed in various laser ablation applications.

42 ENGINEERING↗

DOE OSTI2021

Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE

A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

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 OSTI2020

Shock Initiation of PBX 9502: Material Lots 006 and 007 heated to 76 Celsius

Two shock initiation experiments have been carried out on PBX 9502 material Lots 006 and 007 heated to 76 Celsius. The purpose of these shots is to assess the reactant equation of state and shock-todetonation transition behavior relative to previously characterized lots of heated PBX 9502. PBX 9502 consists of 95% tri-amino-tri-nitro-benzene (TATB) and 5% Kel-F 800 (a copolymer derived from the monomers chlorotrifluorethylene (CTFE) and vinylidine fluoride) by weight. The experiments carried out at the TA-40 Chamber 9 gas gun facility were subject to shock pressures of 12.01 and 11.44 GPa for shots 2s-1129 and 2s-1130 respectively. These pressures resulted in run-to-detonation transitions of 7.91 and 8.73 mm from the impact plane. Unfortunately, due to a series of problems with the heating system the Lot 007 shot was subject to temperature cycling twice before the experiment was fired. The results do not indicate that this had any effect on the sensitivity of the explosive. In terms of shock sensitivity the Lot 006 and 007 PBX 9502 was found to be much less sensitive than the recycled lot 004 when heated to 76 Celsius. The data actually lies closer to the suit of ambient data on the Pop Plot.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Carbon and Hydrocarbon Particle Seeding in Air-Breathing Rotating Detonation Engine

Within the power generation community, the rotating detonation engine (RDE) is only growing in popularity with its increased performance, simple mechanism, and operation. Although significant testing is underway to characterize the RDE for integration with conventional gas turbines, this entire system is still at a relatively low technology readiness level. In the midst of RDE research, there is an initiative to understand solid particle seeding effects in the detonation performance. Under investigation at the University of Central Florida is a Department of Energy (DOE) 15.24 cm (6 in.) RDE, with a solid particle seeder in parallel with its H2 and air flow lines. Previous work on this system involved carbon particle detonation; however, the tested particles were taken one step further to include more sustainable, greener hydrocarbon particles. Testing of powdered sugar, peanut flour, and cornstarch, along with previous carbon black tests have shown not only successful detonability, but a noticeable effect on the detonation wave dynamics. Side-by-side with a particle burning model being developed, an operational map can be determined for the hydrocarbon particles particularly, which can be tuned with the local flow conditions to achieve peak operability while replacing fuels with sustainable alternatives that could even be grown.

Engineering↗

DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, $l_T$ , and its relative ratio with the most energetic length scale of turbulence, $l_T/l_e$ , together with two different levels of the turbulent velocity fluctuation, $u'$ . It is found that $l_T$ / $l_e$ and the ratio of ignition delay time to eddy-turnover time, $\tau _{ig}/\tau _t$ , are the key parameters that control the detonation development. An increase in either $l_T$ or $l_e$ enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, $l_T$ , is significantly modified by high turbulence intensity achieved by a large $l_e$ and $u'$ . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small $l_e$ and a large $u'$ , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either $\tau _{ig}/\tau _t$ less than unity or $l_e = l_T$ even with a large $u'$ is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗

DOE OSTI2022

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗

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↗