Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ignition kinetics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An experimental and modeling study on auto-ignition kinetics of ammonia/methanol mixtures at intermediate temperature and high pressure

A rapid compression machine (RCM) has been applied to measure the ignition delay times of NH 3 /CH 3 OH mixtures covering pressures of 20 and 40 bar, equivalence ratios of 0.5, 1.0 and 2.0, and temperatures between 845 and 1100 K. Here the measurements show that the NH 3 /CH 3 OH mixtures become more reactive with increasing methanol addition. Addition of merely 1% (molar basis) of CH 3 OH to NH 3 lowers the ignition temperature around 100 K at 40 bar in comparison to pure NH 3 . The ignition delay is a complex function of fuel mixture and stoichiometry. For the 1% CH 3 OH mixture, the leaner mixtures are more reactive, while the reverse trend is found for mixtures with 5%, 20% and pure CH 3 OH. Analysis of the pressure profiles shows three distinct ignition modes for NH 3/ CH 3 OH mixtures, facilitated by the pre-ignition heat release from NH 3 consumption. A detailed mechanism for ignition of NH 3 /CH 3 OH fuel blends has been developed, capable of reproducing the ignition behavior of mixtures with reasonable accuracy. A subset for amine / methanol interactions was established, with rate constants for the key reaction between NH 2 and CH 3 OH calculated from ab initio theory. A sensitivity analysis indicates that the critical reactions during the auto-ignition process vary with the CH 3 OH mole fraction in the fuel. The ammonia chemistry, namely NH 2 + NO, NH2 + NO 2 and NH 2 + HO 2 , is dominant for the mixture with 1% CH 3 OH, while the reactions related to CH 3 OH and H 2 O 2 are more important for the 20% CH 3 OH mixture. The interaction between ammonia and methanol shows a more prominent effect on auto-ignition for mixtures with 5% CH 3 OH in fuel as compared to those with 1% and 20% CH 3 OH. According to the modeling results, methanol addition is found to enrich the O/H radical pool, consuming ammonia and promoting auto-ignition through different reaction pathways.

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

A co-kurtosis PCA based dimensionality reduction with nonlinear reconstruction using neural networks

For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Finally, results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.

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

Experimental and modeling study of C2–C4 alcohol autoignition at intermediate temperature conditions

C2–C4 alcohols are advantageous blendstocks identified by many research groups, including the U.S. Department of Energy Co-Optima Initiative, towards enabling efficient, boosted Spark-Ignition (SI) engines. Their use in advanced engine applications requires a comprehensive understanding of their intermediate-temperature autoignition behavior. This work reports an experimental and modeling study covering their fundamental autoignition characteristics in a twin-piston rapid compression machine at pressures of 20 and 40 bar, intermediate temperatures from 750 to 980 K, and two fuel loading conditions representative of boosted SI engines. Direct comparison between these alcohols is made, where the order of reactivity is established across different thermodynamic and fuel loading conditions. Changes in preliminary exothermicity (or intermediate-temperature heat release) displayed in single-stage autoignition across different alcohols and conditions are also quantified. This provides insight into fuel-to-fuel differences, and how these could affect advanced combustion concepts such as spark-assisted compression ignition. Kinetic models are used to simulate the experiments, and reasonable agreement is obtained. Furthermore, the sensitivity analysis results demonstrate the importance of accurately capturing the autoignition kinetics, particularly H-abstraction reactions on the parent fuels by OH and HO 2 , and the branching ratio associated with these.

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

Exploring the potential benefits of high-efficiency dual-fuel combustion on a heavy-duty multi-cylinder engine for SuperTruck I

In support of the Daimler SuperTruck I team’s 55% brake thermal efficiency (BTE) pathway goal, researchers at Oak Ridge National Laboratory performed an experimental investigation of the potential efficiency and emissions benefits of dual-fuel advanced combustion approaches on a modified heavy-duty 15-L Detroit™ DD15 engine. For this work, a natural gas port fuel injection system with an independent injection control for each cylinder was added to the DD15 engine. For the dual-fuel strategies investigated, 65%–90% of the total fuel energy was supplied through the added port fuel injection natural gas (NG) fueling system. The remaining fuel energy was supplied by one or more direct injections of diesel fuel using the production high pressure diesel fueling system. The production DD15 air handling system and combustion geometry were unmodified for this study. Efficiency and emissions with dual-fuel strategies including both low temperature combustion (LTC) and non-LTC approaches such as dual fuel direct-injection were investigated along with control authority over combustion phasing. Parametric studies of dual-fuel NG/diesel advanced combustion were conducted in order to experimentally investigate the potential of high-efficiency, dual-fuel combustion strategies to improve BTE in a multi-cylinder engine, understand the potential reductions in engine-out emissions, and characterize the range of combustion phasing controllability. Characterization of mode transitions from mixing-controlled diesel pilot ignition to kinetically controlled ignition is presented. Key findings from this study included a reproducible demonstration of BTE approaching 48% at up to a 13-bar brake mean effective pressure with significant reductions in engine-out NOx and soot emissions. Additional results from investigating load transients in dual-fuel mode and initial characterization of particle size distribution during dual-fuel operation are presented.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Kinetic simulations of ignited mode cesium vapor thermionic converters

Cesium vapor thermionic converters are an attractive method of converting high-temperature heat directly to electricity, but theoretical descriptions of the systems have been difficult due to the multi-step ionization of Cs through inelastic electron–neutral collisions. This work presents particle-in-cell simulations of these converters, using a direct simulation Monte Carlo collision model to track 52 excited states of Cs. Here, these simulations show the dominant role of multi-step ionization, which also varies significantly based on both the applied voltage bias and pressure. The electron energy distribution functions are shown to be highly non-Maxwellian in the cases analyzed here. A comparison with previous approaches is presented, and large differences are found in ionization rates due especially to the fact that previous approaches have assumed Maxwellian electron distributions. Finally, an open question regarding the nature of the plasma sheaths in the obstructed regime is discussed. The one-dimensional simulations did not produce stable obstructed regime operation and thereby do not support the double-sheath hypothesis.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Misr: Miniature Ignition Screening Rapid Compression Machine for Kinetic Measurements of Novel Fuels (Final Report)

High throughput ignition screening would be useful for evaluating fuels before scaled production. One method for doing this involves generating a map of ignition delay times over different conditions for each fuel blend. Constant condition ignition delay time maps can be used to estimate ignition for engines in prescribed thermodynamic trajectories. A Miniaturized Ignition Screening Rapid Compression Machine (MISR) has been developed and built at the University of Illinois Chicago (UIC) intending to generate those constant condition ignition delay time maps for low temperature chemistry relevant conditions. Some design features include a bellows-actuated driver section, which is lighter compared to other RCM designs, an adjustable brake section, and full optical access. The setup was designed to achieve reaction conditions (T > 600K and P < 100 bar with repetition rates of about 8 seconds per experiment). These conditions were meant to be relevant to low-to-intermediate temperature auto-ignition and for generation of data for developing relevant chemical kinetics models or ignition maps. Because of the compression characteristics and the heat loss characteristics of the MISR, the apparatus currently tests higher temperatures over shorter test times than what would be desirable, however there are some pathways to improved compression characteristics as detailed in this report. A detailed description of the MISR design and its operating mechanism is discussed. The team faced several challenges and setbacks in the construction and commissioning of the MISR, and those are highlighted, along with the modifications that were performed to get the MISR operational. To demonstrate the performance characteristics of the MISR, limited ignition delay experiments were conducted for ethanol fuel at different equivalence ratios (Φ), and those results were usable, but again test a higher temperature and shorter ignition delay time than initially intended.

02 PETROLEUM↗
DOE OSTI2023

Hydrogen and Ammonia Blending With Natural Gas: Ignition Delay Times and Chemical Kinetic Model Validation at Gas Turbine Relevant Conditions

Abstract Ignition delay times from undiluted mixtures of natural gas (NG)/H2/Air and NG/NH3/Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000–1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data is compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH3 fuel addition made no noticeable difference in ignition time. NG/NH3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. The reaction CH3 + O2 = CH3O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.

Pierro, Michael↗
DOE OSTI2026

Methanol mixing−controlled compression ignition enabled via homogeneous charge compression ignition of dimethyl ether through catalytic decomposition of methanol

Methanol is a potentially attractive fuel for marine and off−road engines owing to its availability at bunkering and global distribution locations. Although methanol is well−distributed worldwide, its fuel chemistry and ignition properties make it poorly suited as a direct drop−in replacement for diesel fuel in compression−ignition engines. However, industrial processes are regularly used to convert methanol, via catalytic dehydration, to dimethyl ether (DME) over nonprecious metal catalysts. This chemical conversion can occur at relatively low pressures, temperatures, and catalyst space velocities, highlighting a potential opportunity to generate DME via onboard catalytic dehydration of methanol. DME’s fuel kinetic and ignition properties for compression ignition are much more favorable than those of methanol or even diesel fuel, but DME is more challenging than diesel fuel or methanol to pump, store, and deliver through conventional diesel fueling injection hardware. Thus, a potential opportunity exists to use the ignition and kinetic properties of DME, with the transportation and delivery advantages of methanol, in a methanol−fueled mixing−controlled compression−ignition engine. The present work explores performance, combustion behavior, and emissions reduction opportunities for methanol mixing−controlled combustion, enabled by a HCCI of DME that represents a small fraction of the total fuel energy that can be generated onboard via catalytic dehydration of methanol.

Splitter, Derek [ORNL] (ORCID:0000000174044047)↗
DOE OSTI2024

Ignition Delay Times and Chemical Kinetic Model Validation for Hydrogen and Ammonia Blending With Natural Gas at Gas Turbine Relevant Conditions

Ignition delay times from undiluted mixtures of natural gas (NG)/H 2 /Air and NG/NH 3 /Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000 and 1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data are compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H 2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH 3 fuel addition made no noticeable difference in ignition time. NG/NH 3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. Here, the reaction CH 3 + O 2 = CH 3 O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Validation of a Pressure Dependent PBX 9501 Cookoff Model

A pressure dependent cookoff model for PBX 9501 was developed. This cookoff model was implemented in the finite element (FE) heat transfer code Aria and was used to simulate a set of cookoff experiments. This set is called the Large Scale Annular Cookoff (LSAC) experiments. Three dimensional and axisymmetric heat conduction FE models were implemented to simulate these experiments. The predictive ability of the PBX 9501 pressure dependent cookoff model was evaluated using Latin Hypercube sampling (LHS). Predictions of the thermal times to ignition and temperatures are compared with the experiments.

42 ENGINEERING↗
DOE OSTI2021

Data-driven chemical kinetic reaction mechanism for F-24 jet fuel ignition

A data-driven chemical kinetic mechanism for the military version of Jet A, F-24, is developed for numerical simulations of the ignition process. The main purpose of this study is to obtain a practical F-24 mechanism across wide temperature and equivalence ratio ranges, with a particular focus on the negative temperature coefficient and low temperature regions. The new mechanism (ARLMech-HC-F24) is based on the HyChem model of a similar fuel and optimized using a micro-genetic algorithm against an experimental ignition delay data set of the target fuel. The development and optimization processes include reaction selection, population creation, shuffled tournament implementation based on a merit function, and child-individual creation for the next generation. Several techniques and parameters are proposed to generate an accurate mechanism through an efficient process. The newly introduced data-driven mechanism based on these techniques shows better merit value convergence and represents the ignition behavior more accurately than that without the techniques. This practical mechanism is suitable for the numerical simulations of the F-24 or Jet A ignition problem, and the suggested strategies can be employed in similar problems of rate coefficients estimation.

42 ENGINEERING↗
DOE OSTI2026

An experimental and kinetic modeling study of the ignition of 2-methyl decane

Ignition delay times (IDTs) of 2-methyl decane (C 11 H 24 -2) are measured in a high-pressure shock tube and in a rapid compression machine at equivalence ratios in the range 0.5–2.0 at 90 % dilution, at temperatures in the range 600–1430 K and at pressures of 15 and 30 bar. To clarify the effect of the branched methyl group on fuel oxidation, IDTs of n-undecane (nC 11 H 24 ) are also measured at similar conditions to those measured for C 11 H 24 -2. A new chemical kinetic mechanism, using C3MechV4.0.1 as the core chemistry, is developed and validated against the new experimental data. The thermodynamic properties of the fuel (RH), alkyl (Ṙ), alkyl peroxy (RȮ 2 ), hydroperoxy-alkyl (Q̇OOH), and peroxy hydroperoxy alkyl (Ȯ 2 QOOH) radicals are updated in both the C 11 H 24 -2 and nC 11 H24 models using THERM25. A reaction path flux analysis for C 11 H 24 -2 at different temperatures was conducted. Compared to nC 11 H 24 , C 11 H 24 -2 shows slower reactivity. At low and intermediate temperatures, the chain propagation pathway Q̇OOH ↔ C 11 cyclic ether + ȮH is favored for C 11 H 24 -2, while the chain branching pathway Q̇OOH ↔ Ȯ 2 QOOH ↔ C 11 carbonyl hydroperoxide (KHP) + ȮH is suppressed, leading to lower reactivity compared to nC 11 H 24 . At high temperatures, the presence of the branched methyl group inhibits the direct decomposition of the fuel, resulting in reduced C 2 H 4 formation, which in turn suppresses the reactivity of the fuel.

2-methyl decane↗
DOE OSTI2021

New insights into fuel blending effects: Intermolecular chemical kinetic interactions affecting autoignition times and intermediate-temperature heat release

Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.

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

Experimental evidence of disordered crystalline premixing in sputter-deposited Ni(V)/Al multilayers

The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.

Crystal structure↗
DOE OSTI2022

Computational Fluid Dynamics Modeling of Low Temperature Ignition Processes From a Nanosecond Pulsed Discharge at Quiescent Conditions

Recent interest in nonequilibrium plasma discharges as sources of ignition for the automotive industry has not yet been accompanied by the availability of dedicated models to perform this task in computational fluid dynamics (CFD) engine simulations. The need for a low-temperature plasma (LTP) ignition model has motivated much work in simulating these discharges from first principles. Most ignition models assume that an equilibrium plasma comprises the bulk of discharge kernels. LTP discharges, however, exhibit highly nonequilibrium behavior. In this work, a method to determine a consistent initialization of LTP discharge kernels for use in engine CFD codes like CONVERGE is proposed. The method utilizes first principles discharge simulations. Such an LTP kernel is introduced in a flammable mixture of air and fuel, and the subsequent plasma expansion and ignition simulation is carried out using a reacting flow solver with detailed chemistry. Finally, the proposed numerical approach is shown to produce results that agree with experimental observations regarding the ignitability of methane-air and ethylene-air mixtures by LTP discharges.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Simulation of the PBX 9501 Self-Ignited Heavily Confined Cookoff Test

The PBX 9501 self-ignited thermally damaged explosive (SITX) test was carried out by Holmes et al. A 6-inch diameter sphere made of the PBX 9501 was heated inside a steel confinement vessel until it underwent a thermal explosion. This test provided measured temperature data that was used to define a boundary condition to carry out post-experiment finite element (FE) thermal simulations and comparisons of the FE thermal response with the experimental time to thermal ignition and internal temperature measurements at the center of the PBX 9501 sphere. Finite element models of this experiment were constructed using the heat transfer code Aria. A cookoff model for PBX 9501 was used to simulate the thermal decomposition of PBX 9501 and to predict the time to thermal ignition. An inverse heat conduction scheme was designed and applied to estimate a 1D temperature boundary condition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗