Engineering topics

Mulvihill, Clayton R.

Publications and source records attributed to Mulvihill, Clayton R..

12 records

DOE OSTI2024

Systematic exploration of the thermochemistry for a set of peroxy hydroperoxy-alkyl radicals

Here, the thermochemistry of peroxy hydroperoxy-alkyl ($\overline{O}$OQOOH) radicals has a significant influence on the reactivity of fuels and on the formation of highly oxygenated molecules (HOMs) in the atmosphere. Theoretical characterization of these radicals can be arduous due to their molecular size and complex fundamental interactions, such as hydrogen-bonding and torsional anharmonicity, and difficult to validate in the absence of any direct experimental thermochemical data. In this work, we systematically explore the thermochemistry of a set of $\overline{O}$OQOOH radicals with five increasingly affordable approaches with considerations for these interactions. In doing so, we present a novel conformer selection approach suited to predict properties at combustion temperatures. As a corollary, we also highlight the shortcomings in the standard choice of the ground conformer as the reference. The set of molecules is comprised of 149 C 2 -C 8 $\overline{O}$OQOOH radicals, selected to encompass a wide variety of branching and substitution patterns. Comparisons amongst the approaches help quantify the errors arising from various simplifying assumptions. For the largest of these radicals, and with the most affordable approach, final predictions of Gibbs energies are assigned a 2 σ uncertainty of 4 kcal mol -1 in the negative temperature coefficient region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2023

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Automated identification and calculation of prompt effects in kinetic mechanisms using statistical models

The kinetics of prompt dissociation involves rovibrationally excited species (generally formed by exothermic reactions) which may dissociate or isomerize prior to thermalization via collisions with the bath gas. Treating such rovibrationally excited species (so-called "hot" species) with standard kinetic phenomenology may result in incorrect macroscopic representation of their reactivity. Here this work presents the first fully automated methodology for the calculation of prompt effects of a chosen species in a kinetic mechanism, including (i) reaction selection; (ii) theoretical calculation of rate constants and prompt branching fractions; and (iii) final rate constant fitting. The energy partition between hot fragments is estimated using a variety of statistical models, including a new physically sound microcanonical statistical model based on the rovibrational density of states of the fragments. The methodology is validated against literature data for the prompt dissociations of HCO and C 3 H 7 radicals. The microcanonical statistical model is in better agreement with trajectory simulations for larger species and is thus applicable for practical systems that typically involve large molecules, for which direct dynamics calculations are impractical. The automated workflow is applied to the evaluation of the effects of prompt dissociation for two isomeric radicals C 4 H 7 1-3 (1-methylallyl) and C 4 H 7 1-4 (3-buten-1-yl). Twelve H-atom abstraction reactions are selected and the corresponding rate constants are computed with first principles theory. The microcanonical statistical model predicts that prompt dissociations of C 4 H 7 1-3 and C 4 H 7 1-4 are already significant at 1000K, resulting in differences of up to an order of magnitude at 2000K with respect to the phenomenological thermal rate constants. To illustrate the effects of prompt dissociation on simulations of experimental data, the calculated prompt rate constants are implemented in both CRECK and C3MechV3.3 kinetic mechanisms. Simulations of experimental flame data illustrate the noticeable impact of prompt dissociation kinetics on the high-temperature combustion reactivity of C 4 H 8 -1 and C 4 H 8 -2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Theoretical Kinetics Predictions for Reactions on the NH 2 O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH 2 O potential energy surface, including the NH 2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH 3 oxidation 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 2σ uncertainty) 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 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 constants as well as the branching ratios. We find that the NO + H 2 channel contributes 10% of the total NH 2 + O flux at combustion temperatures, although this channel is not included in modern NH 3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH 3 oxidation and NO x formation. Finally, the present results for NH 2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH 3 and thermal DeNO x .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Quantification of Key Peroxy and Hydroperoxide Intermediates in the Low-Temperature Oxidation of Dimethyl Ether

Dimethyl ether (DME) oxidation is a model chemical system with a small number of prototypical reaction intermediates that also has practical importance for low-carbon transportation. Although it has been studied experimentally and theoretically, ambiguity remains in the relative importance of competing DME oxidation pathways in the low-temperature autoignition regime. To focus on the primary reactions in DME autoignition, we measured the time-resolved concentration of five intermediates, CH 3 OCH 2 OO (ROO), OOCH 2 OCH 2 OOH (OOQOOH), HOOCH 2 OCHO (hydroperoxymethyl formate, HPMF), CH 2 O, and CH 3 OCHO (methyl formate, MF), from photolytically initiated experiments. We performed these studies at P = 10 bar and T = 450-575 K, using a high-pressure photolysis reactor coupled to a time-of-flight mass spectrometer with tunable vacuum-ultraviolet synchrotron ionization at the Advanced Light Source. Our measurements reveal that the timescale of ROO decay and product formation is much shorter than predicted by current DME combustion models. The models also strongly underpredict the observed yields of CH 2 O and MF and do not capture the temperature dependence of OOQOOH and HPMF yields. Adding the ROO + OH $\rightarrow$ RO + HO 2 reaction to the chemical mechanism (with a rate coefficient approximated from similar reactions) improves the prediction of MF. Increasing the rate coefficients of ROO ↔ QOOH and QOOH + O 2 ↔ OOQOOH reactions brings the model predictions closer to experimental observations for OOQOOH and HPMF, while increasing the rate coefficient for the QOOH $\rightarrow$ 2 CH 2 O + OH reaction is needed to improve the predictions of formaldehyde. In conclusion, to aid future quantification of DME oxidation intermediates by photoionization mass spectrometry, we report experimentally determined ionization cross-sections for ROO, OOQOOH, and HPMF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Low- and intermediate-temperature oxidation of dimethyl ether up to 100 atm in a supercritical pressure jet-stirred reactor

Understanding the low- and intermediate-temperature oxidation chemistry of oxygenated fuels like dimethyl ether (DME) at high pressure is paramount to the development of advanced engines with low carbon emissions. The supercritical pressure jet-stirred reactor (SP-JSR) recently developed at Princeton provides a new platform for conducting kinetic studies at low and intermediate temperatures at extremely high pressures with a uniform temperature distribution and a short flow residence time. This paper uses the SP-JSR to investigate DME oxidation at equivalence ratios of 0.175, 1.0, and 1.72, for pressures of 10 and 100 atm, and temperatures ranging from 400 to 900 K. The results demonstrate weakened NTC behavior at 100 atm relative to 10 atm due to increased flux through QOOH + O 2 = O 2 QOOH relative to QOOH = 2 CH 2 O + OH at 100 atm. Furthermore, the intermediate-temperature oxidation window is shifted to lower temperatures at 100 atm. The experimental data are compared with several chemical kinetic models from the literature. The existing models are seen to agree quite well with the experimental data at 10 atm. However, the models fail to properly capture the NTC behavior at 100 atm. Reaction pathway analyses indicate that both the low- and intermediate-temperature chemistries are controlled by RO 2 consumption pathways. The reaction rates for several of the important reactions, such as DME + OH = CH 3 OCH 2 + H 2 O, H 2 O 2 (+M) = 2 OH (+M), and 2 HO 2 = 2 OH + O 2 are updated in this work. Furthermore, the updated model improves the predictability for all key species compared to the original model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Diastereomers and Low-Temperature Oxidation

Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. Furthermore, these results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Isopropanol dehydration reaction rate kinetics measurement using H 2 O time histories

H 2 O formation during thermal decomposition of isopropanol was measured behind reflected shock waves at temperatures ranging from 1127 to 1621 K at an average pressure of 1.42 atm using a laser absorption technique. Of the five modern chemical kinetics models used to compare the H 2 O time histories, the model from Li et al (Combust. Flame 2019;207:171-185) showed the best overall agreement. Sensitivity and rate of production analyses using the Li et al model (as well as those from AramcoMech 3.0, CRECK, and Togbe et al (Energy Fuel. 2011;25:676-683)) showed unimolecular dehydration of isopropanol, iC 3 H 7 OH ⇌ H 2 O + C 3 H 6 (R1), is nearly the sole reaction controlling H2O production at early times, allowing for an a priori measurement of the forward rate constant k 1 . The Arrhenius expression k1 (s –1 ) = 2.60 × 10 13 exp(-31 120 K/T) was determined to represent best the data from this study. According to the models and previous experimental investigations, the pressure investigated is well within the high-pressure limit (HPL) for this reaction. Additional, higher-pressure experiments also confirmed the HPL assumption. An uncertainty analysis was performed by varying secondary reactions within their uncertainties and examining their effect on the overall prediction, establishing an uncertainty within ±20% for all but the highest temperature cases, which have a maximum uncertainty of ±40%. Experiments conducted with a radical trapper, toluene, showed little influence from radical chemistry, suggesting this estimated uncertainty is fairly conservative. Experimental data from Heyne et al (Z Phys Chem. 2015;229:881-907) were found to be in good agreement with the rate measurements from this study and, therefore, a second Arrhenius expression, k 1 (s –1 ) = 2.11 × 10 13 exp(-30 820 K/T), was found to represent both datasets well. This second expression has a larger temperature range of 976-1621 K. Here the present study provides the first high-temperature data collected for this reaction, adding to the limited data available for isopropanol in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗