Atomistic Simulations of Pore Collapse Initiation and Deflagration in HMX
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Not Available
This report describes an improvement to the XDDT code previously described in SAND2026-19856. The improvements include a conversion from CGS units to SI units, a simplified mesh refinement capability, and the addition of tabular gas-phase equation of state (EOS) tables.
The report describes a pipeline to connect a pre-thermal ignition code called XCHEM to a post-thermal ignition code called XDDT. The coupled code enables end-to-end simulation from first exposure to an abnormal thermal environment such as a fire, through thermal ignition, to the final reactive event—a benign pressure burst or a catastrophic detonation. Several demonstration calculations are given for slab and cylindrical geometries. Suggested future work is also discussed.
UL Solutions produced this report as a summary of a test series executed at Sandia National Laboratories that they funded. This report was approved for release by UL Solutions in October 2025 and released in 2026 on their website at https://www.ul.com/insights/methods-and-recommendations-large-scale-deflagration-testing-battery-energy-storage-system.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
We have studied the laser-initiated deflagration of nitromethane (CH 3 NO 2 ) under high static pressure in the diamond anvil cell. Time-resolved images of the deflagrations have been obtained using intensified CCDs (ICCD). In contrast with previous work, we rely on spontaneous emission from the reaction, rather than changes in the speckle pattern produced by artificial illumination of the sample. Furthermore, as opposed to the 1D records obtained previously with streak cameras, ICCDs permit imaging of the burn in two dimensions providing the ability to directly observe anisotropic deflagration behavior. Here, we report several examples of this behavior and discuss its possible origins. We have also investigated the products of the reaction using Raman spectroscopy. At pressures below 25 GPa, the burn product is observed to be opaque and has a Raman spectrum consistent with a carbonaceous soot. At pressures above 25 GPa, the burn product is observed to be transparent (consistent with earlier reports) and its Raman spectrum reveals the presence of molecular N 2 and a nitrogen-carbonate species. The latter species can be recovered to ambient pressure and has a Raman spectrum consistent with ammonium carbonate or ammonium bicarbonate. To aid interpretation of the Raman spectra, measurements have also been made on the isotopologues CH 3 15 NO 2 and 13 CH 3 NO 2 . This work establishes the importance of advanced time-resolved imaging to reveal the details of deflagration under high pressure and further advances our understanding of the corresponding chemistry.
Direct Numerical Simulations (DNS) data of Moderate or Intense Low-oxygen Dilution (MILD) combustion are analysed in this paper to identify the contributions of the autoignition and flame modes. This is performed using an extended Chemical Explosive Mode Analysis (CEMA) which accounts for diffusion effects allowing it to discriminate between deflagration and autoignition. This analysis indicates that in premixed MILD combustion conditions, the main combustion mode is ignition for all dilution and turbulence levels and for the two reactant temperature conditions considered. In non-premixed conditions, the preponderance of the ignition mode was observed to depend on the axial location and mixture fraction stratification. With a large mixture fraction lengthscale, ignition is more preponderant in the early part of the domain while the deflagrative mode increases further downstream. On the other hand, when the mixture fraction lengthscale is small, sequential autoignition is observed. Finally, the various combustion modes are observed to correlate strongly with mixture fraction where lean mixtures are more likely to autoignite while stoichiometric and rich mixtures are more likely to react as deflagrative structures.
The Fusion Z Pinch Experiment (FuZE) investigates sheared-flow stabilization of classic m = 0 and m = 1 instabilities in Z pinches with an embedded axial flow. FuZE consists of a 100 cm coaxial plasma accelerator, where neutral gas is ionized and accelerated in a pulsed electrical discharge, followed by a 50 cm assembly region, where pinches are formed. Maintaining the pinch requires continuous plasma injection provided by a deflagration mode in the coaxial accelerator. Two discharge modes, with and without deflagration, are investigated on FuZE. Pinch formation is observed with deflagration only. Plasma velocities in the assembly region are found to match the E × B velocity estimated in the accelerator based on a 1D circuit model, indicating that a 1D MHD approximation may offer a valid description of the plasma in the accelerator channel. The velocity of magnetic field propagation is found to agree with the snowplow model based on momentum conservation, and the lifetime of the pinch is shown to be in agreement with constraints imposed by mass conservation.
In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.
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.
Lawrence Livermore National Laboratory (LLNL) is collaborating with Colorado School of Mines (Mines), to field the first dynamic experiments on the 3-ring flash radiography system (3-ring). One set of experiments for 3-ring concerns the deflagration to detonation transition (DDT), within a sealed column of high-explosives. The high-explosive material is ignited at one end with a BK-80, an igniter made from Boron Potassium-Nitrate (BKNO 3 ). The initial deflagration reaction will increase the pressure in the sealed vessel, resulting in the formation of a dense "plug" of material ahead of the deflagration front. The transition to detonation is thought to be related to the plug formation. Previous experiments have been able to resolve the plug geometry, but only in a two dimensional radio-graph. 3-ring offers a novel diagnostic tool which could provide three dimensional density fields via few-view reconstruction algorithms and computed tomography.
On April 11, 2018, four drums containing transuranic waste at the Idaho National Laboratory underwent over-pressurization, ejecting their lids and spreading radiological waste within a facility. An investigation has found that waste in the drums generated methane gas, which contributed to the event. Subsequent to the investigation, the potential for drums to have methane and other flammable analytes whose concentrations could approach or exceed the lower flammability limit (LFL) and the adequacy of the controls to prevent or mitigate a possible deflagration was evaluated. An extensive review of the historical records was performed to determine how many drums exceed the LFL. The historical record identified a small quantity of drums that exceeded the LFL for xylene, hydrogen, and methane. The primary codified applicable code or standard for handling drums is 29 CFR 1910.120 (j) 'Handling Drums and Containers.' It is used throughout the commercial and government sector. This code is integrated throughout the Fluor Idaho safety management programs and procedures to control the Standard Industrial Hazards (SIHs) associated with drum and container handling. Example requirements include: - Drums and containers are inspected, and their integrity assured prior to being moved. - Site operations are organized to minimize the amount of drum or container movement. - When there is a reasonable possibility of flammable atmospheres being present, material handling equipment and hand tools are of the type to prevent sources of ignition. Drums and containers under pressure, as evidenced by bulging or swelling, are not moved until the cause for excess pressure is determined and appropriate containment procedures have been implemented to protect employees from explosive relief of the drum. Drums which exceed the WIPP FGA limit have NCRs associated with them. Once a drum fails flammable gas sampling, an NCR is initiated. The NCR drives placement into NCR dense pack rows. A standing order drives placement of the FGA failures to a single planer segregation location. These drums are not allowed for further processing or shipment until the NCR has been cleared. Once in the single planer rows, follow-up FGA testing is performed by CCP. The NCR is dispositioned by follow-up sampling occurring after adequate time has been provided for the gas to diffuse through the filter assembly, or, in some cases, the drums are repackaged into additional drums to reduce the source term. To date, the ICP project has performed approximately 1.5 million drums moves without a deflagration event due to combustible gas generation. Flammable gas generation in a transuranic waste drum is not unique to the ICP but is common across the DOE complex. Based on the experience at Idaho, application of industry standards is sufficient to control the risk of drum deflagration due to drum movements. (authors)
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.
While low-temperature combustion (LTC) strategies have been found to mitigate nitrogen oxides and particulate matter emissions in diesel engines, studies have also associated LTC with an increase in unburnt hydrocarbons. With more recent studies on diesel after end-of-injection (AEOI), combustion recession is identified as a phenomenon where at near nozzle region, high-temperature ignition (HTI) combustion can propagate back to the nozzle tip consuming the unburnt hydrocarbons AEOI. Current literature has suggested that combustion recession is controlled by auto-ignition. However, high-fidelity simulations and detailed analysis of such a mechanism are missing. In this study, comprehensive Large Eddy Simulations of a reacting spray at “Spray A” conditions are performed, where detailed analysis of combustion recession concerning flame morphology and propagation modes are included. In particular, this study demonstrated for the first time that while combustion recession is mainly auto-ignition dominated (consistent with the literature), a cool flame was found to deflagrate towards the richer regions of the mixture, promoting mixing and increasing the mixture temperature. This leads to HTI kernels, which then grow and develop as deflagrative waves, therefore sustaining the combustion recession process. The study also detailed the extinction mechanism of combustion: the entrainment wave will overlean the near-nozzle mixtures, rendering it unable to support HTI, which leads to the extinction of the upstream flame AEOI in lower reactivity mixtures. Combustion recession is also observed to be contingent on the chemical and diffusion processes, even at low scalar dissipation rates. Finally, a new criterion for combustion recession based on chemical explosive mode is proposed and validated with previous combustion recession index to quantify the extent of HTI in near-nozzle mixtures AEOI. The newly developed metric combined with a previous experimentally-based metric can provide simple but valuable measurements of the degree and propensity of the upstream flame AEOI.