Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Spectrally accurate, reverse-mode differentiable bounce-averaging algorithm and its applications

We present a fast, spectrally (exponentially) accurate, automatically differentiable bounce-averaging algorithm that is used to simplify kinetic models. Using this algorithm, implemented in the DESC stellarator optimisation suite, we can perform efficient optimisation of many objectives to improve stellarator performance, such as the effective ripple 𝜖 eff metric for the neoclassical transport coefficient in the low collisionality regime and proxies for energetic particle confinement. For the first time, we optimise a finite-beta stellarator to directly reduce neoclassical ripple transport using reverse-mode differentiation. This ensures the computational cost of differentiation is independent of the number of controllable parameters.

fusion plasma↗
DOE OSTI2018

Improved yield and control of spectra from high-intensity laser-generated neutron beams

Abstract Kinetic modeling of laser-ion beam generation from the “break-out afterburner” (BOA) has been modeled for several deuteron-rich solid-density target foils. Modeling the transport of these beams in a beryllium converter shows as much as a fourfold increase in neutron yield over the present state of the art through the use of alternative target materials. Additionally, species-separation dynamics during the BOA can be exploited to control the hardness of the neutron spectra, of interest for, for example, enhancing penetrability in shielded material in active neutron interrogation settings.

laser based neutron sources, active interrogation↗
DOE OSTI2022

Prospects and Limitations of Predicting Fuel Ignition Properties from Low-Temperature Speciation Data

Using chemical kinetic modeling and statistical analysis, we investigate the possibility of correlating key chemical “markers”–typically small molecules–formed during very lean (φ ~ 0.001) oxidation experiments with near-stoichiometric (φ ~ 1) fuel ignition properties. One goal of this work is to evaluate the feasibility of designing a fuel-screening platform, based on small laboratory reactors that operate at low temperatures and use minimal fuel volume. Buras et al. [Combust. Flame2020,216, 472–484] have shown that convolutional neural net (CNN) fitting can be used to correlate first-stage ignition delay times (IDTs) with OH/HO2 measurements during very lean oxidation in low-T flow reactors with better than factor-of-2 accuracy. In this work, we test the limits of applying this correlation-based approach to predict the low-temperature heat release (LTHR) and total IDT, including the sensitivity of total IDT to the equivalence ratio, φ. We demonstrate that first-stage IDT can be reliably correlated with very lean oxidation measurements using compressed sensing (CS), which is simpler to implement than CNN fitting. LTHR can also be predicted via CS analysis, although the correlation quality is somewhat lower than for first-stage IDT. In contrast, the accuracy of total IDT prediction at φ = 1 is significantly lower (within a factor of 4 or worse). Furthermore, these results can be rationalized by the fact that the first-stage IDT and LTHR are primarily determined by low-temperature chemistry, whereas total IDT depends on low-, intermediate-, and high-temperature chemistry. Oxidation reactions are most important at low temperatures, and therefore, measurements of universal molecular markers of oxidation do not capture the full chemical complexity required to accurately predict the total IDT even at a single equivalence ratio. As a result, we find that φ-sensitivity of ignition delay cannot be predicted at all using solely correlation with lean low-T chemical speciation measurements.

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

The Influence of Intake Pressure and Ethanol Addition to Gasoline on Single- and Dual-Stage Autoignition in an HCCI Engine

Autoignition in HCCI engines is known to be controlled by the combustion kinetics of the in-cylinder fuel/air mixture which is highly influenced by the amount of low-temperature and intermediate-temperature heat release (LTHR and ITHR) that occurs. At lower intake pressures (typically <1.4 bar absolute), it has been observed that gasoline behaves as a single-stage heat release fuel, while at higher intake pressures (typically >1.8 bar absolute) gasoline behaves as a two-stage heat release fuel. Furthermore, ethanol blending into gasoline strongly affects heat release characteristics, and this warrants further investigation. This paper experimentally investigates the conditions under which gasoline transitions from a single-stage heat release fuel to a two-stage heat release fuel as intake pressure is increased. Experiments were performed in single-cylinder HCCI engine fueled with two research-grade gasolines, FACE A and FACE C. These gasolines were tested neat, and with 10% and 20% (by volume) ethanol addition. In addition, these results were compared to results previously obtained for PRF 85, and new results for PRF 84 with 10% and 20% ethanol addition. Moreover, the engine experiments were supported by rapid compression machine (RCM) ignition delay data for the same fuels. The engine experiments revealed that there were minimal differences between the heat release profiles of the two gasolines, FACE A and FACE C, a result which was supported by the RCM experiments that showed similar ignition delay times for the two FACE fuels and PRF 84. On the other hand, with ethanol addition to these gasolines and PRF 84, the occurrence of LTHR shifted to higher intake pressures compared to ethanol-free cases, from 1.4 bar intake pressure for neat fuel to 2.2 bar with 20% ethanol. Consequently, the intake temperatures required to achieve constant combustion phasing for all mixtures were drastically altered. Simulations using a detailed chemical kinetic model were utilized to understand the effects of ethanol blending on the ignition characteristics of PRF 84. The addition of ethanol was found to act as a radical sink where it inhibits the radical pool formation during the low (<850 K) and intermediate (850–1050 K) temperature chemistry regimes resulting in lower reactivity. These results help explain ethanol’s significant antiknock qualities under boosted conditions in spark-ignition engines.

42 ENGINEERING↗
DOE OSTI2025

Coupling of Lipid Peroxidation and Criegee Intermediate Mediated Autoxidation in the Heterogeneous Oxidation of Linoleic Acid Aerosols

Autoxidation is an established mechanism for the degradation of organic molecules, which is relevant in the atmosphere, combustion processes, the environment, and the rancidification of lipids (commonly known as lipid peroxidation). Autoxidation proceeds via radical chain reactions involving hydroxyl (•OH), peroxy (RO2•), and alkoxy radicals, which are also prominent oxidants in the atmosphere. Recent reports have provided evidence for an alternative autoxidation mechanism driven instead by Criegee intermediates (CIs), which are produced from the reaction of β-hydroxy peroxy radicals (β-OH-RO2•). This work evaluates the contributions of these two mechanisms in the •OH initiated heterogeneous oxidation of linoleic acid (LA) aerosols. Reaction kinetics and product distributions are monitored using a vacuum ultraviolet photoionization aerosol time-of-flight mass spectrometer. To explain the observed kinetics, a kinetic model is developed that incorporates both the conventional peroxidation and alternative CI-mediated autoxidation mechanisms. We observe that the CI-mediated autoxidation pathways enhance the heterogeneous autoxidation rate, while the peroxidation reactions, although present, contributes less to the overall oxidation rate. α-Acyloxyalkyl hydroperoxides (AAHPs) are identified as key indicators for bimolecular reactions of CI with LA, highlighting the role of LA as a CI scavenger. Moreover, the measured functionalized LA products with hydroxyl or carbonyl group(s), serve as markers for the peroxidation reactions. In summary, this work presents a quantitative framework to understand the coupled reaction network of •OH, RO2•, β-OH-RO2• radicals, and CI in driving heterogeneous autoxidation, which is crucial for understanding degradation mechanisms of organic molecules in the environment and atmosphere.

Criegee intermediate↗
DOE OSTI2026

Dramatic Slowdown of Bimolecular Reactions of Criegee Intermediates in Organic Aerosols

Although Criegee intermediates (CIs) play a central role in driving the transformation of organic molecules, their condensed-phase reaction kinetics remains poorly understood. While CIs undergo unimolecular decomposition to produce OH radicals, stabilized CIs participate in a host of other bimolecular reactions forming larger molecular weight products, including secondary ozonides (kSOZ), alkoxyalkyl hydroperoxides (kAlAHP), and acyloxyalkyl hydroperoxides (kAcAHP). Here, we demonstrate, using prior literature and new measurements of heterogeneous organic aerosol reactions using O3 or OH radicals, that these condensed-phase bimolecular reaction rate coefficients kSOZ (10-20 cm3·molecule-1·s-1), kAlAHP, and kAcAHP (10-19 cm3·molecule-1·s-1) are 7-10 orders of magnitude slower than analogous gas-phase reactions (∼ 10-10-10-13 cm3·molecule-1·s-1). This contrasts with current evidence that the rate coefficients for unimolecular steps in the condensed phase are similar in magnitude to those in the gas phase. Kinetic models using gas-phase-like bimolecular coefficients consistently overestimate the importance of bimolecular relative to unimolecular pathways while underpredicting the aerosol reactivity by up to 30-fold. The slowing of bimolecular reactivity in aerosols is likely due to solvation effects and the increased steric hindrance of larger molecular weight CIs and their scavengers in the condensed phase. These results challenge the conventional view that the high reactant density of an organic aerosol should favor bimolecular over unimolecular reactions, underscoring the need for more accurate CI condensed-phase rate constants for reliable modeling of oxidative processes across a diverse set of atmospheric and environmental systems.

Zeng, Meirong↗
DOE OSTI2017

Production of Hydrogen Peroxide in Groundwater at Rifle, Colorado

The commonly held assumption that photodependent processes dominate H 2 O 2 production in natural waters has been recently questioned. Here, we present evidence for the unrecognized and light-independent generation of H 2 O 2 in groundwater of an alluvial aquifer adjacent to the Colorado River near Rifle, CO. In situ detection using a sensitive chemiluminescent method suggests H 2 O 2 concentrations ranging from lower than the detection limit (<1 nM) to 54 nM along the vertical profiles obtained at various locations across the aquifer. Our results also suggest dark formation of H 2 O 2 is more likely to occur in transitional redox environments where reduced elements (e.g., reduced metals and NOM) meet oxygen, such as oxic-anoxic interfaces. A simplified kinetic model involving interactions among iron, reduced NOM, and oxygen was able to reproduce roughly many, but not all, of the features in our detected H 2 O 2 profiles, and therefore there are other minor biological and/or chemical controls on H 2 O 2 steady-state concentrations in such aquifer. Because of its transient nature, the widespread presence of H 2 O 2 in groundwater suggests the existence of a balance between H 2 O 2 sources and sinks, which potentially involves a cascade of various biogeochemically important processes that could have significant impacts on metal/nutrient cycling in groundwater-dependent ecosystems, such as wetlands and springs. More importantly, our results demonstrate that reactive oxygen species are not only widespread in oceanic and atmospheric systems but also in the subsurface domain, possibly the least understood component of biogeochemical cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Kinetics of Direct Olefin Synthesis from Syngas over Mixed Beds of Zn–Zr Oxides and SAPO-34

A packed bed containing a physical mixture of both Zn–Zr mixed oxide catalyst and SAPO-34 converts syngas directly into a mixture of C 2 –C 5 olefins and paraffins. Specifically, the mixed oxide catalyst is responsible for intermediate oxygenate synthesis from syngas while the molecular sieve catalyzes olefin synthesis from the oxygenate intermediates. Kinetic measurements with cofed propylene over each catalyst independently confirm olefin hydrogenation activity over both components of the composite bed. The addition of either water or CO to the feed drops the activity of propylene hydrogenation over the Zn–Zr oxide. In sum, under reaction conditions of syngas feed and produced water, olefin hydrogenation predominantly occurs over the SAPO-34 catalyst, rather than over the catalyst responsible for hydrogenating CO into oxygenate intermediates. Here, a developed kinetic model consistent with this conclusion describes measurements at differing feed compositions, temperatures, space velocities, and bed catalyst mixing ratios. Technoeconomic analysis of the process indicates that the olefin-to-paraffin ratio is a key performance metric for commercial scale syngas conversion and highlights the importance of considering olefin hydrogenation rates over the molecular sieve component.

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

Group Additivity Determination for Oxygenates, Oxonium Ions, and Oxygen-Containing Carbenium Ions

Bio-oil produced from biomass fast pyrolysis often requires catalytic upgrading to remove oxygen and acidic species over zeolite catalysts. The elementary reactions in the mechanism for this process involve carbenium and oxonium ions. In order to develop a detailed kinetic model for the catalytic upgrading of biomass, rate constants are required for these elementary reactions. The parameters in the Arrhenius equation can be related to thermodynamic properties through structure-reactivity relationships, such as the Evans-Polanyi relationship. For this relationship, enthalpies of formation of each species are required, which can be reasonably estimated using group additivity. However, the literature previously lacked group additivity values for oxygenates, oxonium ions, and oxygen-containing carbenium ions. In this work, 71 group additivity values for these types of groups were regressed, 65 of which had not been reported previously and six of which were newly estimated based on regression in the context of the 65 new groups. Heats of formation based on atomization enthalpy calculations for a set of reference molecules and isodesmic reactions for a small set of larger species for which experimental data was available were used to demonstrate the accuracy of the Gaussian-4 quantum mechanical method in estimating enthalpies of formation for species involving the moieties of interest. Isodesmic reactions for a total of 195 species were constructed from the reference molecules to calculate enthalpies of formation that were used to regress the group additivity values. The results showed an average deviation of 1.95 kcal/mol between the values calculated from Gaussian-4 and isodesmic reactions versus those calculated from the group additivity values that were newly regressed. Importantly, the new groups enhance the database for group additivity values, especially those involving oxonium ions.

09 BIOMASS FUELS↗
DOE OSTI2022

Thermal and Catalytic Decomposition of 2-Hydroxyethylhydrazine and 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid

To develop chemical kinetics models for the combustion of ionic liquid-based monopropellants, identification of the elementary steps in the thermal and catalytic decomposition of components such as 2-hydroxyethylhydrazinium nitrate (HEHN) is needed but is currently not well understood. The first decomposition step in protic ionic liquids such as HEHN is typically the proton transfer from the cation to the anion, resulting in the formation of 2-hydroxyethylhydrazine (HEH) and HNO 3 . In the first part of this investigation, the high-temperature thermal decomposition of HEH is probed with flash pyrolysis (<1400 K) and vacuum ultraviolet (10.45 eV) photoionization time-of-flight mass spectrometry (VUV-PI-TOFMS). Next, the investigation into the thermal and catalytic decomposition of HEHN includes two mass spectrometric techniques: (1) tunable VUV-PI-TOFMS (7.4-15 eV) and (2) ambient ionization mass spectrometry utilizing both plasma and laser ionization techniques whereby HEHN is introduced onto a heated inert or iridium catalytic surface and the products are probed. The products can be identified by their masses, their ionization energies, and their collision-induced fragmentation patterns. Further, formation of product species indicates that catalytic surface recombination is an important reaction process in the decomposition mechanism of HEHN. The products and their possible elementary reaction mechanisms are discussed.

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

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

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

Resolving Discrepancies between State-of-the-Art Theory and Experiment for HO 2 + HO 2 via Multiscale Informatics

Recent high-level theoretical calculations predict a mild temperature dependence for HO 2 + HO 2 inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO 2 + HO 2 rate constants nearly identical to theoretical predictions---implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Here, similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate---identified by recent theory but absent from earlier interpretations---yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

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

The Identity and Chemistry of C 7 H 7 Radicals Observed during Soot Formation

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C 7 H 7 radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H 7 , the results clearly demonstrate that the isomeric composition of C 7 H 7 strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C 7 H 7 isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C 7 H 7 reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C 7 H 7 isomers and predict negligible contributions from other C 7 H 7 isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

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

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗
DOE OSTI2016

Quantification of the Keto-Hydroperoxide (HOOCH 2 OCHO) and Other Elusive Intermediates during Low-Temperature Oxidation of Dimethyl Ether

This paper provides new temperature-dependent mole fractions of elusive intermediates relevant to the low-temperature oxidation of dimethyl ether (DME). It extends the previous study of Moshammer et al. [ J. Phys. Chem. A 2015, 119, 7361–7374] in which a combination of a jet-stirred reactor and molecular beam mass spectrometry with single-photon ionization via tunable synchrotron-generated vacuum-ultraviolet radiation was used to identify (but not quantify) several highly oxygenated species. Here, temperature-dependent concentration profiles of 17 components were determined in the range of 450–1000 K and compared to up-to-date kinetic modeling results. Special emphasis is paid toward the validation and application of a theoretical method for predicting photoionization cross sections that are hard to obtain experimentally but essential to turn mass spectral data into mole fraction profiles. The presented approach enabled the quantification of the hydroperoxymethyl formate (HOOCH 2 OCH 2 O), which is a key intermediate in the low-temperature oxidation of DME. Finally, the quantification of this keto-hydroperoxide together with the temperature-dependent concentration profiles of other intermediates including H 2 O 2 , HCOOH, CH 3 OCHO, and CH 3 OOH reveals new opportunities for the development of a next-generation DME combustion chemistry mechanism.

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

Aerosol Fragmentation Driven by Coupling of Acid–Base and Free-Radical Chemistry in the Heterogeneous Oxidation of Aqueous Citric Acid by OH Radicals

A key uncertainty in the heterogeneous oxidation of carboxylic acids by hydroxyl radicals (OH) in aqueous-phase aerosol is how the free-radical reaction pathways might be altered by acid-base chemistry. In particular, if acid-base reactions occur concurrently with acyloxy radical formation and unimolecular decomposition of alkoxy radicals, there is a possibility that differences in reaction pathways impact the partitioning of organic carbon between the gas and aqueous phases. To examine these questions, a kinetic model is developed for the OH-initiated oxidation of citric acid aerosol at high relative humidity. The reaction scheme, containing both free-radical and acid-base elementary reaction steps with physically validated rate coefficients, accurately predicts the experimentally observed molecular composition, particle size, and average elemental composition of the aerosol upon oxidation. The difference between the two reaction channels centers on the reactivity of carboxylic acid groups. Free-radical reactions mainly add functional groups to the carbon skeleton of neutral citric acid, because carboxylic acid moieties deactivate the unimolecular fragmentation of alkoxy radicals. In contrast, the conjugate carboxylate groups originating from acid-base equilibria activate both acyloxy radical formation and carbon-carbon bond scission of alkoxy radicals, leading to the formation of low molecular weight, highly oxidized products such as oxalic and mesoxalic acid. Subsequent hydration of carbonyl groups in the oxidized products increases the aerosol hygroscopicity and accelerates the substantial water uptake and volume growth that accompany oxidation. These results frame the oxidative lifecycle of atmospheric aerosol: it is governed by feedbacks between reactions that first increase the particle oxidation state, then eventually promote water uptake and acid-base chemistry. When coupled to free-radical reactions, acid-base channels lead to formation of low molecular weight gas-phase reaction products and decreasing particle size.

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

Molecular Dynamics Study of Combustion Reactions in a Supercritical Environment. Part 2: Boxed MD Study of CO + OH → CO 2 + H Reaction Kinetics

Oxy-fuel combustion technology holds a great promise in both increasing the efficiency of the energy conversion and reducing environmental impact. Yet, effects of the higher pressures and replacement of the nitrogen with carbon dioxide diluent are not well understood at present. The title reaction is one of the most important processes in combustion. Despite numerous studies, the effects of supercritical carbon dioxide environment did not receive much attention in the past. In this work, we report the results of boxed molecular dynamics simulations of these effects at QM/MM theory level with periodical boundary conditions. The free energy barriers for HOCO intermediate formation and decomposition were tabulated in a wide range of pressures (1–1000 atm) and temperatures (400–1600 K). Pressure dependence of calculated rate constants for these reaction steps and overall reaction were analyzed. We discovered that the CO 2 environment may increase these rate constants up to a factor of 25, at near critical conditions. At higher temperatures, this effect weakens significantly. Numerical values for parameters of extended Arrhenius equation, suitable for combustion kinetic modeling are reported.

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

Large Intermediates in Hydrazine Decomposition: A Theoretical Study of the N 3 H 5 and N 4 H 6 Potential Energy Surfaces

Large complex formation involved in the thermal decomposition of hydrazine (N 2 H 4 ) is studied here using transition state theory based theoretical kinetics. A comprehensive analysis of the N 3 H 5 and N 4 H 6 potential energy surfaces was performed at the CCSD(T)-F12a/aug-cc-pVTZ//ωB97x-D3/6-311++G(3df,3pd) level of theory, and pressure-dependent rate coefficients were determined. There are no low-barrier unimolecular decomposition pathways for triazane (n-N 3 H 5 ), and its formation becomes more significant as the pressure increases; it is the primary product of N 2 H 3 + NH 2 below 550, 800, 1150, and 1600 K at 0.1, 1, 10, and 100 bar, respectively. The N 4 H 6 surface has two important entry channels, N 2 H 4 + H 2 NN and N 2 H 3 + N 2 H 3 , each with different primary products. Interestingly, N 2 H 4 + H 2 NN primarily form N 2 H 3 + N 2 H 3 , while disproportionation of N 2 H 3 + N 3 H 3 predominantly leads to the other N 2 H 2 isomer, HNNH. Stabilized tetrazane (n-N 4 H 6 ) formation from N 2 H 3 + N 2 H 3 becomes significant only at relatively high pressures and low temperatures due to fall-off back into N 2 H 3 + N 2 H 3 . Pressure-dependent rate coefficients for all considered reactions as well as thermodynamic properties of triazane and tetrazane, which should be considered for kinetic modeling of chemical processes involving nitrogen- and hydrogen-containing species, are reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗