Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive burn”

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 OSTI2022

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗
DOE OSTI2023

Post-detonation fireball modeling: Validation of freeze out approximations

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

AWSD reactive flow model for PBX 9404

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive flow model for high explosive PBX 9404 is developed. We specifically calibrate an AWSD model for PBX 9404 by fitting equations of state for reactants and detonation products to the results of thermochemical calculations and to experimental data from multiple sources. The calibrated equations of state are then coupled with an Arrhenius rate law based on shock temperature that describes the reaction progress during PBX 9404 detonation. The parameter values in the rate law are calibrated to experimental gas-gun data and diameter effect results. The results of the calibrated AWSD model are in strong agreement with available experimental data for PBX 9404. A similar level of agreement between predicted and experimental results is observed when the calibrated model is validated on data that were not used in the model parameterization procedure. Our results illustrate that the AWSD model is capable of accurately describing the many important properties and observables in the reactive burn of PBX 9404. Because of the historical significance of PBX 9404 in high explosives research and its current use in aging studies, this work provides an important model of a legacy material, which can be used to make comparisons to new high explosive formulations.

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

Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives

A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2022

Optimization of artificial viscosity in production codes based on Gaussian Regression surrogate models

To accurately model flows with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity has to be introduced to convert kinetic energy into internal energy, thereby increasing the entropy across shocks. Determining the appropriate strength of the artificial viscosity is an art and strongly depends on the particular problem and experience of the researcher. The objective of this study is to pose the problem of finding the appropriate strength of artificial viscosity as an optimization problem and solve this problem using machine learning (ML) tools, specifically using surrogate models based on Gaussian Process regression and Bayesian analysis. We describe the optimization method and discuss various practical details of its implementation. The shock-containing problems for which we apply this method all have been implemented in the LANL code FLAG. First, we apply ML to find optimal values to isolated shock problems of different strengths. Second, we apply ML to optimize viscosity for a 1D propagating detonation problem based on Zel’dovich-von Neumann-Doring (ZND) detonation theory using a reactive burn model. We compare results for default (currently used values in FLAG) and optimized values of artificial viscosity for these problems demonstrating the potential for significant improvement in the accuracy of computations.

42 ENGINEERING↗
DOE OSTI2023

Pinwheel Experiment Confirmatory Shot 1 Pre-Shot Report

The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into an aluminum sandwiched PBX-9501. The follow-on experiment will drive the aluminum cylinder into a tantalum sandwiched PBX-9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.

42 ENGINEERING↗
DOE OSTI2023

Notes on Davis EOS: Historical Perspective, Derivations & Physical Considerations

Davis reactant and product semi-empirical EOS have been introduced in a series of publications by W. C. Davis. These EOS are convenient to use in conjunction with reactive burn models to parameterize the thermodynamic response of quiescent high explosives (HE) and the mixture of HE decomposition products. There is a history of development of these EOS with different formulae and functional forms found in different papers. The goal of the present notes is to track that history, clarify the physical motivation behind various modifications to Davis EOS over time, and collect all the relevant formulae in one place.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Pinwheel Experiment Tantalum Shot 1 Pre-Shot Report

The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into a tantalum sandwiched PBX‐9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Burnable absorbers in nuclear reactors – A review

Burnable absorbers can benefit nuclear reactors of virtually any design by providing reactivity control for extended fuel cycles, tritium production, burning of long-lived radionuclides, and reactor safety. When selecting the ideal burnable absorber type and its design, one must consider the resulting impact on the reactor’s fuel cycle design and cost, reactivity, thermal hydraulics, manufacturing, and radiation response. These selection criteria, as well as neutronic and thermophysical material property requirements, may be vastly different depending on whether the burnable absorber is intended for use in a commercial water-cooled reactor, a research reactor, or a next-generation advanced reactor system. A recent integration, synthesis, description of past and present technologies, and identification of existing gaps and areas of future research is lacking on these important topics. Here, this paper includes a fundamental description of the use of burnable absorbers and their impacts on reactivity, absorber depletion, self-shielding, basic thermophysical properties, and the use of burnable absorbers in next-generation nuclear applications..

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Combustion of metal powder with dinitrogen tetroxide

In this work, we present analysis of a novel reactive material system that employs dinitrogen tetroxide (N 2 O 4 ) as a liquid oxidizer with metal powder fuels. The oxidizer was added to micron scale aluminum and zirconium powders by a remote injection system. When ignited with a high voltage spark, the mixtures were observed to possess reactivity comparable to nanocomposite reactive materials, with open-tube flame expansion velocities from 500 to 1400 m/s depending on fuel/oxidizer ratio and tube diameter. Temperatures were observed to range from 3000 to 3500 K as measured with gray body fits to 16-channel time-resolved pyrometer and time-integrated spectrometer data. These values were significantly below calculated adiabatic flame temperatures, which we attribute to local deviations from stoichiometry and kinetic/energetic limitations similar to those observed in studies of particles burning in high pressures gaseous oxidizers. Al/N 2 O 4 reactivity was found to be most likely limited by the vaporization of the metal from the particle surface and Zr/N 2 O 4 was limited by slower burning and complex interactions involving the solubility of nitrogen and oxygen in the molten Zr. We also discuss the potential for these materials to be used to create an “on/off” reactive material, since N2O4 can be added remotely and driven to evaporate via vacuum or purge of inert gases to return it to a safe condition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurement report: Extreme heat and wildfire emissions enhance volatile organic compounds in a temperate forest

Climate extremes are projected to cause unprecedented deviations in the emission and transformation of volatile organic compounds (VOCs), which trigger feedback mechanisms that will impact the atmospheric oxidation and formation of aerosols and clouds. However, the response of VOCs to future conditions such as extreme heat and wildfire events is still uncertain. This study explored the modification of the mixing ratio and distribution of several anthropogenic and biogenic VOCs in a temperate oak–hickory–juniper forest as a response to increased temperature and transported biomass burning plumes. A chemical ionization mass spectrometer was deployed on a tower at a height of 32 m in rural central Missouri, United States, for the continuous and in situ measurement of VOCs from June to August of 2023. The maximum observed temperature in the region was 38 °C, and during multiple episodes the temperature remained above 32 °C for several hours. Biogenic VOCs such as isoprene and monoterpene followed closely the temperature daily profile but at varying rates, whereas anthropogenic VOCs were insensitive to elevated temperature. During the measurement period, wildfire emissions were transported to the site and substantially increased the mixing ratios of acetonitrile and benzene, which are produced from burning of biomass. An in-depth analysis of the mass spectra revealed more than 250 minor compounds, such as formamide and methylglyoxal. Extreme heat and presence of wildfire plumes modified the overall volatility, reactivity, O : C, and H : C ratios of the extended list of VOCs. The calculated OH reactivities during extreme temperature condition and transport of biomass burning plumes were 106.37±4.27 and 106.22±5.15 s −1 , respectively, which are substantially higher than background level of 98.78±1.16 s −1 . Multivariate analysis also clustered the compounds into five factors, which highlighted the sources of the unaccounted-for VOCs. Ultimately, results here underscore the effect of extreme heat and wildfire emissions on the overall chemical properties VOC in a temperate forest.

Salvador, Christian Mark [Oak Ridge National Labor↗
DOE OSTI2024

Reactivity of net-zero carbon alcohol fuels and their corresponding aldehyde intermediates on PGM-based commercial oxidation catalysts for lean-burn emissions control

Alcohol net-zero carbon fuels will play a significant role in decarbonization of the hard-to-electrify transportation sectors. However, the combustion process in alcohol-fueled engines generate toxic aldehydes. For commercialization of net-zero carbon fuels on engines, compliance with the stringent U.S. EPA regulations is necessary. Here, this contribution focuses on reactivity of alcohols and aldehydes on commercial diesel oxidation catalysts (DOC) to inform selection of future net-zero carbon fuels. Alcohol (and aldehyde) light-off temperatures were measured on aged DOCs under full synthetic engine-exhaust conditions. Methanol was the most reactive, while ethanol completely oxidized at higher temperatures. Alcohols typically formed less reactive aldehyde intermediates. DRIFTS revealed strongly adsorbing surface formates and acetates during methanol and ethanol oxidation, respectively, resulting in inhibition effects. Preliminary alcohol oxidation mechanisms combining flow reactor and DRIFTS observations are presented. Aldehydes inhibited CO oxidation. Minimal N 2 O formation was observed for the alcohols investigated.

Alcohol oxidation↗
DOE OSTI2022

Path to Increasing p-B11 Reactivity via ps and ns Lasers

The Lawson criterion for proton-boron (p-11B) thermonuclear fusion is substantially higher than that for deuterium-tritium (DT) because the fusion cross section is lower and peaks at higher ion energies. The Maxwellian averaged p-11B reactivity peaks at several hundred keV, where bremsstrahlung radiation emission may dominate over fusion reactions if electrons and ions are in thermal equilibrium and the losses are unrestricted. Nonequilibrium burn has often been suggested to realize the benefits of this aneutronic reaction, but the predominance of elastic scattering over fusion reactivity makes this difficult to achieve. The development of ultrashort pulse lasers (USPL) has opened new possibilities for initiating nonequilibrium thermonuclear burns and significant numbers of p-11B alpha particles have been reported from several experiments. We present an analysis that shows that these significant alpha yields are the result of beam fusion reactions that do not scale to net energy gain. We further find that the yields can be explained by experimental parameters and recently updated cross sections such that a postulated avalanche mechanism is not required. We use this analysis to understand the underlying physics of USPL-driven nonequilibrium fusion reactions and whether they can be used to initiate fusion burns. We conclude by outlining a path to increasing the p-11B reactivity towards the goal of achieving ignition and describing the design principles that we will use to develop a computational point design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Plasma-assisted deflagration to detonation transition in a microchannel with fast-frame imaging and hybrid fs/ps coherent anti-Stokes Raman scattering measurements

Our study examines kinetic enhancement by nanosecond dielectric barrier discharge (ns-DBD) plasma on fuellean dimethyl ether (DME), oxygen (O 2 ), and argon (Ar) premixtures during deflagration to detonation transition (DDT) experiments in a microchannel. Non-equilibrium plasma produces active species and radicals and creates fast and slow heating of a mixture to promote ignition due to electronic and vibrational excitation. Experiments have been conducted to examine the influence of the plasma discharge on the premixture and on the resultant deflagration to detonation transition (DDT) onset time and distance through the use of high speed imaging and one-dimensional, two-beam, femtosecond/picosecond, coherent anti-Stokes Raman scattering (CARS). A highspeed camera is used to trace the time histories of flame front position and velocity and to identify the dynamics and onset of DDT. The results show that plasma discharge can nonlinearly affect the onset time and distance of DDT. It is shown that a small number of plasma discharge pulses prior to ignition result in reduced DDT onset time and distance by 60% and 40%, respectively, when compared to the results without pre-excitation by ns discharges. The results also show that an increase of plasma discharge pulses results in an extended DDT onset time and distance of 224% and 94%, respectively. Time history of the deflagration wave speed of DME and the analysis of ignition timescale under the choking condition of the deflagration front suggest low temperature ignition may play a role for DME near the isobaric choking condition of the burned gas and the DDT. Plasmainduced conversion of the reactive mixture was assessed via the O 2 to CO 2 ratio as measured through fs/ps CARS during the DBD discharges. The present experiments demonstrate the ability of non-equilibrium plasma to alter the chemistry of DME/O 2 /Ar premixtures in order to control DDT for applications in advanced propulsion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗