Search NASA⌕ Search

SEARCH · Search NASA

Results for “deflagration-to-detonation transition”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Observation of asymmetric explosive density evolution in the deflagration-to-detonation transition for porous explosives

We report deflagration-to-detonation transition phenomena in cylindrical columns of porous explosives based on octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX). The HMX powder and LX-14 (95.5 wt. % HMX and 4.5 wt. % estane polyurethane) prills were observed with high-speed cameras and flash x rays under conditions that allow for precise scaling and registration of the recorded images. These experimental results test the conventional assumption of axisymmetric, effective one-dimensional burning. Polycarbonate confinement tubes allowed this direct imaging with steel spheres employed as radio-opaque fiducials. Spheres embedded in the explosive column revealed the displacement of the explosive prior to deflagration. X-ray and fast camera images show that a dense plug of compact explosive develops ahead of the deflagration front. The internal fiducials register both internal changes in the plug and changes in the position between successive x-ray images. Tantalum witness strips placed on the inner wall of the explosive channel revealed radial expansion of the tube bore at successive x-ray exposures, indicating the internal pressure of about 100 MPa just before detonation. Images indicate that the deflagration is not only asymmetric but also involves newly observed transient phenomena just before detonation. Furthermore, we typically see a brief but continuous increase in deflagration speed and the formation of a deflagration channel that bypasses the plug just prior to detonation. We hypothesize that these phenomena play an important role in the deflagration-to-detonation transition itself.

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

The Deflagration-to-Detonation Transition in Two Dimensions

The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figures 1 and 2 for an overview of the proceeding steps. It has been assumed in the past that the same mechanism which drives DDT in the one dimensional case will apply to both two and three dimensions. However, this has never been experimentally verified and it is possible that an entirely new mechanism drives DDT higher dimensions. In order to test this a series of two dimensional DDT tests were performed to help provide empirical evidence of the mechanism in higher dimensions.

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

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Type Ia supernovae deflagration-to-detonation transition explosions powered by the Zel’dovich reactivity gradient mechanism

ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.

Brooker, E.↗
DOE OSTI2026

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↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computational diagnostics for flame acceleration and transition to detonation in a hydrogen/air mixture

A new computational diagnostic method for pressure-induced compressibility is proposed by projecting its local contribution to the chemical explosive mode (CEM) in the chemical explosive mode analysis (CEMA) framework. The new method is validated for the study of detonation development during the deflagration-to-detonation transition (DDT) process. The flame characteristics are identified through the quantification of individual CEM contributions of chemical reaction, diffusion, and pressure-induced compressibility. Numerical simulations are performed to investigate the DDT processes in a stoichiometric hydrogen-air mixture. A Godunov algorithm, fifth-order in space, and third-order in time are used to solve the fully compressible Navier-Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model (CDM) described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The new diagnostic method is first applied to onedimensional (1D) canonical flame configurations followed by two-dimensional (2D) simulations of DDT in an obstructed channel where different detonation initiation scenarios are examined using the new CEMA projection formulation. Detailed examinations of the idealized configuration of detonation initiation through shock focusing mechanism at a flame front are also studied using the new formulation. A comparison of the currently proposed CEMA projection and the original formulation by the authors suggests that including the pressure-induced compressibility is essential for the use of CEMA in DDT process. The results also show that the new formulation of CEMA projection can successively capture the detonation initiation through either a gradient mechanism or a direct initiation mechanism, and therefore can be used as an effective local analytical tool for the computational diagnostics of detonation initiation in a DDT process. It was found that detonation development is characterized by a strong contribution of chemistry role to the CEM which is pivotal to the initiation of detonation. The role of compressibility is found enhanced at the edge of the detonation front where diffusion was found to have minimal effects on detonation development.

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

On The Deflagration to Detonation Transition in 1.65 g cm -3 Pentaerythritol Tetranitrate

The LANL Weapons Response Group is charged with making determinations for a wide range of scenarios regarding the probability that inadvertent ignition of the high explosive (HE) components might cause nuclear yield, or dispersal of special nuclear material, during assembly/disassembly operations at the Pantex Plant. These scenarios are categorized based on the nature of the ignition causing stimulus, i.e. impact, electrostatic discharge, pressure-shear, etc. However, for all of these scenarios, regardless of how ignition is achieved, there must also be a credible path for the burning HE to transition to detonation for yield or dispersal to occur. If it can be determined that there is no feasible path for deflagration-to-detonation transition (DDT), the scenario space that requires mitigative action is greatly reduced. The strategy to support Weapons Response determinations is to use experimentation and observation to evaluate HE response in severe “overtest” configurations. To increase with confidence in the margins of HE response, these overtests are designed, by controlling parameters known to enhance HE response violence, to be more severe than what is ever encountered in the actual operations environment at the Plant. This study is focused on DDT in pentaerythritol tetranitrate (PETN), a high explosive that is used in LANL’s detonators. PETN is known to be able to DDT at certain densities, and this represents a safety concern. The PETN in LANL’s detonators is commonly a density of 1.65 g cm -3 . However, DDT has never been observed at this density in PETN. In this work, we studied the DDT reaction in 1.65 g cm -3 PETN to determine whether DDT is possible at this density, and if so, what length of PETN is required for the transition to occur. In any of the tests conducted, we did not observe DDT in 1.65 g cm -3 PETN. It is possible that the DDT length for 1.65 g cm -3 PETN is greater than 6” (the available run length in the experiments), but this length is simply not available in a detonator, and is therefore not a concern.

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

IHE Material and IHE Subassembly Qualification Test Description and Criteria (v. 14.7)

The purpose of the deflagration-to-detonation test is to demonstrate that an IHE material will not undergo deflagration-to-detonation transition (DDT) under stockpile relevant conditions of scale, confinement, and material condition. Inherent in this test design is the assumption that ignition does occur, with onset of deflagration. The test design will incorporate large margins and replicates to account for the stochastic nature of DDT events.

36 MATERIALS SCIENCE↗
DOE OSTI2020

HMX detonation calculation using multi-reaction chain

Detonation of HMX has been simulated using a multi-reaction chain and species equation-of-states. A reacting flow code capable of handling arbitrary number of reactions and species has been used. Computational results are compared to the results produced by a pressure-based single reaction model (AJEX) and a temperature-based single-reaction model. Results show that the multi-reaction approach is capable of producing detonation calculations comparable to the AJEX model calculations, demonstrating that the multi-reaction approach possesses the properties required for the detonation modeling, a prerequisite of deflagration-to-detonation transition modeling.

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

Structural Effects of Postulated Hydrogen Explosions in Process Piping and Vessels - 20274

At Savannah River Site (SRS), preparations to startup new liquid waste operations and recent advances in understanding thermolytic hydrogen generation have led to several recent studies on the structural effects of postulated hydrogen explosions. This paper reports on seven representative piping and vessel cases that are evaluated for deflagrations, detonations, and deflagration-to-detonation transition (DDT) assuming the initial mixture is stoichiometric hydrogen and air. The evaluations are based on the approach used by the Department of Defense for analyzing structures subjected to accidental explosions and include inelasticity, dynamic structural amplification, and material strain-rate effects. Evaluation criteria consider the potential for repeated explosions causing fatigue but also consider that some degree of permanent plastic deformation is acceptable as long as rupture is prevented. Specific conclusions depend on the assumed pre-ignition pressure. If the pre-ignition pressure is atmospheric and if the components are designed to a code from the American Society of Mechanical Engineers, the components evaluated are expected to plastically deform but not rupture under the explosion scenarios considered. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗