SEARCH · Engineering Papers

Results for “ignition”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2024

The role of stereochemistry in combustion processes

Stereochemical effects significantly influence chemical processes, yet it is not well understood if they are a leading source of uncertainty in combustion modeling. Stereochemistry influences a combustion model (i) at the earliest stage of its construction when mapping the reaction network, (ii) in the computation of individual thermochemical and rate parameters, and (iii) in the prediction of combustion observables. The present work reviews the importance of enumerating stereochemical species and reactions at each of these steps. Further, it analyzes the separate influence of several types of stereochemistry, including geometric, optical, and fleeting transition state diastereomers. Three reaction networks serve to examine which stages of low-temperature oxidation are most affected by stereochemistry, including the first and second oxidation of n-butane, the third oxidation of n-pentane, and the early stages of pyrolysis of 1- and 2-pentene. The 149 reactions in the n-butane mechanism are expanded to 183 reactions when accounting for diastereomerism. Each of these 183 reactions is parameterized with ab initio kinetics computations to determine that, for the n-butane mechanism, the median factor of diastereomeric deviation is 3.5 at 360 K for rate constants and as high as 1.6 for mechanism reactivity, in terms of ignition delay times, as opposed to a mechanism without stereochemical expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗

DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, \(l_T\) , and its relative ratio with the most energetic length scale of turbulence, \(l_T/l_e\) , together with two different levels of the turbulent velocity fluctuation, \(u'\) . It is found that \(l_T\) / \(l_e\) and the ratio of ignition delay time to eddy-turnover time, \(\tau _{ig}/\tau _t\) , are the key parameters that control the detonation development. An increase in either \(l_T\) or \(l_e\) enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, \(l_T\) , is significantly modified by high turbulence intensity achieved by a large \(l_e\) and \(u'\) . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small \(l_e\) and a large \(u'\) , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either \(\tau _{ig}/\tau _t\) less than unity or \(l_e = l_T\) even with a large \(u'\) is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗

DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗

DOE OSTI2022

Quantifying Firebrand Production and Transport Using the Acoustic Analysis of In-Fire Cameras

Firebrand travel and ignition of spot fires is a major concern in the Wildland-Urban Interface and in wildfire operations overall. Firebrands allow for the efficient breaching across fuel-free barriers such as roads, rivers and constructed fuel breaks. Existing observation-based knowledge on medium-distance firebrand travel is often based on single tree experiments that do not replicate the intensity and convective updraft of a continuous crown fire. Recent advances in acoustic analysis, specifically pattern detection, has enabled the quantification of the rate at which firebrands are observed in the audio recordings of in-fire cameras housed within fire-proof steel boxes that have been deployed on experimental fires. The audio pattern being detected is the sound created by a flying firebrand hitting the steel box of the camera. This technique allows for the number of firebrands per second to be quantified and can be related to the fire's location at that same time interval (using a detailed rate of spread reconstruction) in order to determine the firebrand travel distance. A proof of concept is given for an experimental crown fire that shows the viability of this technique. When related to the fire's location, key areas of medium-distance spotting are observed that correspond to regions of peak fire intensity. Trends on the number of firebrands landing per square metre as the fire approaches are readily quantified using low-cost instrumentation.

58 GEOSCIENCES↗

DOE OSTI2021

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗

DOE OSTI2020

Some Observations About the Drop-weight Explosive Sensitivity Test

The apparent simplicity of the drop-weight apparatus for explosive sensitivity testing hides the reality that it is actually a complex integrated test of both ignition and growth of deflagration. Although the drop-weight test is undeniably a useful screening test for explosive properties, a misunderstanding of the technique’s limitations has blinded many researchers to its limited wider applicability. This monograph discusses how the test actually works, the significant engineering difficulties with standardization between machines, which types of explosives are suited to the test and which are not, and finally offers a few suggestions for alternatives when a more quantified understanding of a material’s response is required for other applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Multiple injection for improving knock, gaseous and particulate matter emissions in direct injection SI engines

Advances in fuel injector technology have enabled research and development on a large variety of direct injection spark ignition (DISI) engine fueling strategies targeted to improve engine performance and reduce engine-out emissions. This study explores the effect of multiple injections on knock, engine efficiency and stability, and particulate number and gaseous emissions on a single-cylinder DISI research engine. Work to date on multiple injections in the literature is reviewed, and then two aspects of multiple injection strategies are experimentally investigated: the number of injections (up to five times in a cycle); and the timing of the injections (classified relative to the timing of intake valve opening and closing). A boosted, single-cylinder research engine equipped with a state-of-the-art piezoelectric hollow cone spray direct injector and research-grade E10 gasoline was used for the study. The results show multiple injections maintain torque and combustion stability compared with single injection and slightly increase the knock limits and indicated thermal efficiencies (maximum 0.7% absolute, 2.8% relative compared with single-injection baseline) due to improved heat release phasing, especially with an additional late injection during the intake valve closed (compression stroke) period. The gaseous pollutant emissions including nitrogen oxides and unburned hydrocarbons were reduced by 25% with multiple injections, particularly with injection during the compression stroke. In contrast, carbon monoxide emissions increased with multiple injections for all conditions. Increasing the number of injection events significantly reduced particulate number emissions (half with each additional injection), and the decrease in particulate number was not sensitive to the injection timing.

42 ENGINEERING↗

DOE OSTI2021

A novel machine learning-based optimization algorithm (ActivO) for accelerating simulation-driven engine design

A novel design optimization approach (ActivO) that employs an ensemble of machine learning algorithms is presented. The proposed approach is a surrogate-based scheme, where the predictions of a weak leaner and a strong learner are utilized within an active learning loop. The weak learner is used to identify promising regions within the design space to explore, while the strong learner is used to determine the exact location of the optimum within promising regions. For each design iteration, exploration is done by randomly selecting evaluation points within regions where the weak learner-predicted fitness is high. The global optimum obtained by using the strong learner as a surrogate is also evaluated to enable rapid convergence once the most promising region has been identified. First, the performance of ActivO was compared against five other optimizers on a cosine mixture function with 25 local optima and one global optimum. In the second problem, the objective was to minimize indicated specific fuel consumption of a compression-ignition internal combustion (IC) engine while adhering to desired constraints associated with in-cylinder pressure and emissions. In this work, the efficacy of the proposed approach is compared to that of a genetic algorithm, which is widely used within the internal combustion engine community for engine optimization, showing that ActivO reduces the number of function evaluations needed to reach the global optimum, and thereby time-to-design by 80%. Furthermore, the optimization of engine design parameters leads to savings of around 1.9% in energy consumption, while maintaining operability and acceptable pollutant emissions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A survey study on arsenic speciation in coal fly ash and insights into the role of coal combustion conditions

Coal fly ashes (CFAs) are the low-density byproducts of the coal combustion process. Improper or uncontrolled CFA disposal poses significant environmental and health concerns due to the potential leaching of toxic heavy metals such as arsenic (As). Previous studies have investigated the content and speciation of As in different CFA samples, yet systematic information on As speciation in CFA with representative coal source and combustion conditions is still missing. Based on a recent survey study on the typical coal sources and combustion conditions across the U.S., this study selected 19 representative CFA samples to systematically investigate As speciation and potential correlations with these parameters. The composition, morphology, mineralogy, and As speciation of these CFA samples were characterized by complementary analytical, microscopic, and spectroscopic techniques. Synchrotron X-ray spectroscopy and microscopy analyses revealed the dominant As oxidation state to be As(V) and with strong associations to Ca, with the exception of 3 samples that had 19–51% As(III), likely due to the use of selective catalytic reduction (SCR) process. Principal component analysis was conducted to identify potential correlations of As concentration and oxidation state with parameters such as major element content, loss on ignition (LOI), average particle size, coal source, and combustion condition. Al 2 O 3 and FeO content were found to capture a majority of the variability. Further, results from this study provide fundamental basis for understanding the correlations between coal source, combustion conditions, CFA characteristics, and As speciation, and providing insights for downstream beneficial utilization or disposal management.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗

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 OSTI2021

Finite-rate entrainment effects on nitrogen oxide (NO x ) emissions in staged combustors

Axially staged combustors offer the potential for achieving low NO x emissions from gas turbines at the elevated temperatures required to achieve at least 65% combined cycle efficiencies. While there is evidence that premixing between the main burner products and secondary stream reduces NO x , it is important to characterize the specific requirements needed to achieve low emissions at gas turbine conditions. This study examines the sensitivity of NO x to finite-rate, large-scale entrainment of the main and secondary streams under the simplifying assumption of infinitely-fast small-scale mixing. Essentially, this allows us to isolate the effects of limited large-scale entrainment rates by using a homogeneous/uniform condition for the entrained and reacting gases. We use a reduced-order reactor network model to examine a generic staged-combustor whose inputs are physical time scales that embody the finite-rate entrainment characteristics. With simulations conducted over a large parameter space at typical operating conditions (25 atm, 650 K air), the results show that, even when small-scale mixing is infinitely fast and the reaction zone is uniform, entrainment significantly affects NO x emissions due to its influence on the equivalence ratio — and thus the temperature and time — at which the entrained mixture ignites. Fuel-air staging is shown to be vital to NO x reduction in most practical cases where entrainment times exceed roughly 1 ms and the secondary stream finishes entraining before the main burner products. A constrained NO x minimization indicates that using the leanest possible main burner and re-routing as much air as possible to the secondary stage is key to low NO x with finite-rate entrainment; for example, less than 10 ppm NO x can be achieved if the secondary stage contains 20% of the total mass flow.

03 NATURAL GAS↗

DOE OSTI2021

Probing intermediate temperature heat release in autoignition of C3-C4 iso-alcohol/gasoline blends

In this work, we report an experimental and modeling study on the blending effects of C3–C4 iso-alcohols, namely iso-propanol and iso-butanol, on the characteristics of intermediate temperature heat release (ITHR) for a research-grade gasoline (FACE-F) in a rapid compression machine operating at diluted/stoichiometric fuel loading, compressed pressure of 43 bar, and compressed temperatures from 700 to 970 K. Changes in ITHR behavior are characterized through ITHR extent and evolution at 0 to 30 vol% iso-alcohol blending levels. Experimental observations reveal the strong promoting and decelerating effects on ITHR extent and evolution, respectively, within the low-temperature regime for both iso-alcohols, with stronger effects seen for iso-propanol; within the intermediate-temperature regime, the influence of both iso-alcohols diminishes greatly. Comprehensive chemical kinetic modeling is undertaken using a recently developed gasoline surrogate/alcohol model in conjunction with a five-component gasoline surrogate (FGF-LLNL), with good agreement obtained with the experiments for all the blends. Sensitivity, rate of production and flux analyses highlight the importance of chemical kinetic interactions via both fuel-specific and non-fuel-specific reactions. Disabling all the chemical kinetic interactions between iso-alcohol and FGF-LLNL in carbonated/non-fuel-specific species leads to somewhat extended ITHR duration and reduced ignition reactivity. Direct comparison between FGF-LLNL/iso-propanol and FGF-LLNL/iso-butanol blends is also made, where iso-propanol and iso-butanol are found to influence ITHR characteristics via significantly different chemical kinetic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Repetitively pulsed nanosecond discharge plasma decay in propane–oxygen gas mixture in the presence of a heating electric field

Plasma decay in the afterglow of a repetitively pulsed nanosecond discharge in a stoichiometric propane–oxygen mixture was experimentally investigated when a weak heating DC electric field was applied and in its absence. The discharge was ignited at room gas temperature and a pressure of 1–2 Torr and was characterized by low specific energy inputs (<0.004 eV per molecule in one pulse). Using microwave interferometry, the temporal evolution of the electron density during plasma decay was studied, and the effective recombination coefficients were obtained from data processing. It was shown that the rate of plasma decay behaved in a non-monotonic manner with increasing degree of propane oxidation; at first the decay rate grew, then passed through a maximum, fell and saturated in the limit of a large (~2000) number of pulses. In this limit, the effect of the heating DC electric field on the plasma decay decreased with approaching chemical equilibrium. Numerical simulation of the observed effects was performed for low and high oxidation degrees of propane taking into account changes in the composition of positive ions in the plasma. Good agreement was obtained between measurements and calculations of the electron density during plasma decay in these cases. Here it was shown that the formation of cluster ions in the discharge afterglow plays a fundamental role. The plasma decay was controlled by electron recombination with hydrocarbon cluster ion at low oxidation degree of propane and with water cluster (hydrated) ions at high oxidation degree. A hypothesis was proposed to explain the observed nonmonotonic behavior of the plasma decay rate with an increase in the propane oxidation in the discharge, based on the formation of hydrated hydrocarbon ions C x H y + (H 2 O) k at moderate oxidation degrees.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Skeletal model reduction with forced optimally time dependent modes

Skeletal model reduction based on local sensitivity analysis of time dependent systems is presented in which sensitivities are modeled by forced optimally time dependent (f-OTD) modes. The f-OTD factorizes the sensitivity coefficient matrix into a compressed format as the product of two skinny matrices, i.e. f-OTD modes and f-OTD coefficients. The modes create a low-dimensional, time dependent, orthonormal basis which capture the directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active species, and reaction paths. Evolution equations for the f-OTD modes and coefficients are derived, and the implementation of f-OTD for skeletal reduction is described. For demonstration, skeletal reduction is conducted of the constant pressure ethylene-air burning in a zero-dimensional reactor, and new reduced models are generated. The laminar flame speed, the ignition delay, and the extinction curve as predicted by the models are compared against some existing skeletal models in literature for the same detailed model. Finally, the results demonstrate the capability of f-OTD to eliminate unimportant reactions and species in a systematic, efficient and accurate manner.

42 ENGINEERING↗

DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗