Engineering Papers⌕ Search

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.

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 OSTI2020

Response of PETN detonators to elevated temperatures

In the present study, commercially available detonators with pentaerythritol tetranitrate (PETN) were subjected to elevated temperatures. The detonators were thermally ignited over a range of heating rates to measure ignition delay time and assess detonator violence. The violence of the detonator was quantified by measuring the velocity of the detonator closure disc (or “flyer”). The maximum flyer velocity of a thermally ignited detonator was comparable in magnitude to that obtained by initiating a room temperature pristine detonator with an exploding bridge wire (under the same confinement); however, the high flyer velocity was not an indication of deflagration to detonation transition (DDT) in the thermally ignited detonator. The detonator responded more violently than a thermally ignited detonator when initiated at 95% of the ignition delay time. Furthermore, inoperability thresholds were also measured by varying the detonator temperature and the threshold was found to be sensitive at detonator temperatures below the melting point of PETN.

42 ENGINEERING↗
DOE OSTI2021

A reactive flow model for the 3,3'-diamino-4,4'-azoxyfurazan based plastic bonded explosive (PBX 9701)

Here, this paper describes the calibration of the Arrhenius Wescott–Stewart–Davis (AWSD) reactive flow model for the recently developed high explosive PBX 9701, which consists of 97% 3,3'-diamino-4,4'-azoxyfurazan (DAAF) and 3% FK-800 binder by weight. DAAF-based explosive formulations have several desirable qualities as they are relatively insensitive to non-shock insults but have higher performance than triaminotrinitrobenzene based formulations. Equations of state for the explosive reactants and products are calibrated using a combination of existing and new experimental data and theoretical calculations. The AWSD rate law calibration utilizes both one-dimensional shock-to-detonation and multidimensional rate stick experiments to capture the shock initiation and propagating detonation regimes. Validation of the calibrated model is demonstrated through comparison with recent gas-gun experiments.

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

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

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

Operability thresholds for thermally damaged EBW detonators

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (~406 K). Furthermore, the temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

42 ENGINEERING↗
DOE OSTI2020

Understanding Microstructure Variability in Vapor-Deposited Energetic Materials by Using Phase-Field Methods

Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CISME report on experiments with HMX powder, PBX 9501 prills, and heated PBX 9501

The ability to predict post-ignition explosive response is relevant to explosives safety, particularly weapons assembly/disassembly operations that must operate under a sufficiently conservative safety basis to meet risk mitigation criteria for inadvertent nuclear detonation (IND). In particular, hazard determinations hinge on whether the PBX 9501 main charge can ultimately attain detonation from an initial, non-shock, ignition stimulus. The CISME experiment uses a spherical explosive charge, which is thermally ignited in the center using a fiber-delivered infrared laser pulse. This ignition mechanism - i.e. a localized thermal hotspot - is representative of the ignition-causing insults typical in accident scenarios, including low-velocity impact, drop, skid, and punch events. In fiscal year 2018 CISME experiments were conducted with pristine PBX 9501 spheres, confined and unconfined. Sub-detonative, low-violence reactions were observed. Here we report experiments conducted since 2018. The explosive charge was varied: experiments were conducted with room temperature pour-density HMX powder and prills of PBX 9501, and heated, thermally damaged PBX 9501. Experiments were conducted in a confined configuration, with two 4 mm-thick stainless-steel hemispherical shells bolted around the sphere. Detonation, via DDT (deflagration-to-detonation transition) was observed in the 3-inch diameter powdered HMX and the 6-inch diameter heated and thermally-damaged PBX 9501 tests. No transition to DDT was observed in smaller diameters for those configurations, and no DDT was observed for pour-density prills of PBX 9501 of any diameter. The original test series demonstrated that the response of pristine, undamaged PBX 9501 to a central ignition is relatively benign. This new test series shows that increased porosity, lack of polymer binders and thermal damage can greatly increase the violence of the post-ignition response, up to and including detonation. These results bridge the explosive response of PBX 9501 from sub-detonative through detonation as a function of diameter and provide valuable test cases for modeling efforts that strive to capture the DDT threshold. By exploring explosive response in a more representative 3D geometry, we can also benchmark these new results to more plentiful DDT datasets generated from testing in 1D tube confinement (the historical test configuration)–the configuration upon which almost all DDT models were developed.

42 ENGINEERING↗
DOE OSTI2020

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

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

Engineering the Microstructure and Morphology of Explosive Films via Control of Interfacial Energy

Physical vapor deposition of organic explosives enables growth of polycrystalline films with a unique microstructure and morphology compared to the bulk material. This study demonstrates the ability to control crystal orientation and porosity in pentaerythritol tetranitrate films by varying the interfacial energy between the substrate and the vapor-deposited explosive. Additionally, variation in density, porosity, surface roughness, and optical properties is achieved in the explosive film, with significant implications for initiation sensitivity and detonation performance of the explosive material. Various surface science techniques, including angle-resolved X-ray photoelectron spectroscopy and multiliquid contact angle analysis, are utilized to characterize interfacial characteristics between the substrate and explosive film. Optical microscopy and scanning electron microscopy of pentaerythritol tetranitrate surfaces and fracture cross sections illustrate the difference in morphology evolution and the microstructure achieved through surface energy modification. X-ray diffraction studies with the Tilt-A-Whirl three-dimensional pole figure rendering and texture analysis software suite reveal that high surface energy substrates result in a preferred (110) out-of-plane orientation of pentaerythritol tetranitrate crystallites and denser films. Low surface energy substrates create more randomly textured pentaerythritol tetranitrate and lead to nanoscale porosity and lower density films. This work furthers the scientific basis for interfacial engineering of polycrystalline organic explosive films through control of surface energy, enabling future study of dynamic and reactive detonative phenomena at the microscale. Results of this study also have potential applications to active pharmaceutical ingredients, stimuli-responsive polymer films, organic thin film transistors, and other areas.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2020

Development of a Safe, Reliable, and Versatile Explosive Demolition Initiation System

Demolition contractors have a need to simultaneously detonate many demolition charges safely, cheaply, and reliably. An explosive distribution company has contracted the Initiation Systems Group at Lawrence Livermore Laboratory as consultants. The task is to develop a new, safer initiation system to detonate anywhere from 1 to 2000 demolition charges. The primary user of the system will be the explosive demolition contractors.

42 ENGINEERING↗
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 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

Numerical investigation of mass flow rate effects on multiplicity of detonation waves within a H 2 /Air rotating detonation combustor

This paper presents results from numerical simulations of a non-premixed hydrogen-air rotating detonation combustor with radial injection. The fuel and air mass flow rate are varied in order to hold a unity global equivalence ratio. The calculations show that multiple detonation waves co-exist when the mass flow rate is increased. Conditional statistics of the detonation structure and combustion processes suggest similarities across co-existing waves. Quantification of the injection response to the rotation of a detonation indicates that at higher flow rate the refill time is short enough to allow for a quick and well mixed composition prior to the new front passage. Details of the combustion characteristics are analyzed. Finally, the results elucidate the correlation between initial injection conditions and detonation multiplicity on the overall physics within the combustor.

08 HYDROGEN↗
DOE OSTI2022

Investigating growth to detonation in vapor-deposited hexanitrostilbene and pentaerythritol tetranitrate films using high-throughput methods

In this work, a high-throughput experimental setup was used to characterize initiation threshold and growth to detonation in the explosives hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN). The experiment sequentially launched an array of laser-driven flyers to shock samples arranged in a 96-well microplate geometry, with photonic Doppler velocimetry diagnostics to characterize flyer velocity and particle velocity at the explosive–substrate interface. Vapor-deposited films of HNS and PETN were used to provide numerous samples with various thicknesses, enabling characterization of the evolution of growth to detonation. One-dimensional hydrocode simulations were performed with reactions disabled to illustrate where the experimental data deviate from the predicted inert response. Prompt initiation was observed in 144 μm thick HNS films at flyer velocities near 3000 m/s and in 125 μm thick PETN films at flyer velocities near 2400 m/s. This experimental setup enables rapid quantification of the growth of reactions in explosive materials that can reach detonation at sub-millimeter length scales. These data can subsequently be used for parameterizing reactive burn models in hydrocode simulations, as discussed in Paper II [D. E. Kittell, R. Knepper, and A. S. Tappan, J. Appl. Phys. 131, 154902 (2022)].

36 MATERIALS SCIENCE↗
DOE OSTI2022

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗