Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reaction”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reaction Layer Formation on MgO in the Presence of Humidity

Mineralization by MgO is an attractive potential strategy for direct air capture (DAC) of CO 2 due to its tendency to form carbonate phases upon exposure to water and CO 2 . Hydration of MgO during this process is typically assumed to not be rate limiting, even at ambient temperatures. However, surface passivation by hydrated phases likely reduces the CO 2 capture capacity. Here, we examine the initial hydration reactions that occur on MgO(100) surfaces to determine whether they could potentially impact CO 2 uptake. We first used atomic force microscopy (AFM) to explore changes in reaction layers in water (pH = 6 and 12) and MgO-saturated solution (pH = 11) and found the reaction layers on MgO are heterogeneous and nonuniform. To determine how relative humidity (R.H.) affects reactivity, we reacted samples at room temperature in nominally dry N 2 (~11–12% R.H.) for up to 12 h, in humid (>95% R.H.) N 2 for 5, 10, and 15 min, and in air at 33 and 75% R.H. for 8 days. X-ray reflectivity and electron microscopy analysis of the samples reveal that hydrated phases form rapidly upon exposure to humid air, but the growth of the hydrated reaction layer slows after its initial formation. Reaction layer thickness is strongly correlated with R.H., with denser reaction layers forming in 75% R.H. compared with 33% R.H. or nominally dry N 2 . The reaction layers are likely amorphous or poorly crystalline based on grazing incidence X-ray diffraction measurements. After exposure to 75% R.H. in air for 8 days, the reaction layer increases in density as compared to the sample reacted in humid N 2 for 5–15 min. Further, this may represent an initial step toward the crystallization of the reaction layer. Overall, high R.H. favors the formation of a hydrated, disordered layer on MgO. Based on our results, DAC in a location with a higher R.H. will be favorable, but growth may slow significantly from initial rates even on short timescales, presumably due to surface passivation.

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

Chemical Differences between Phenolic Secondary Organic Aerosol Formed through Gas-Phase and Aqueous-Phase Reactions

Phenolic compounds, which are significant emissions from biomass burning (BB), undergo rapid photochemical reactions in both gas and aqueous phases to form secondary organic aerosol, namely, gasSOA and aqSOA, respectively. The formation of gasSOA and aqSOA involves different reaction mechanisms, leading to different product distributions. In this study, we investigate the gaseous and aqueous reactions of guaiacol a representative BB phenol to elucidate the compositional differences between phenolic aqSOA and gasSOA. Aqueous-phase reactions of guaiacol produce higher SOA yields than gas-phase reactions (e.g., roughly 60 vs 30% at one half-life of guaiacol). These aqueous reactions involve more complex reaction mechanisms and exhibit a more gradual SOA evolution than their gaseous counterparts. Initially, gasSOA forms with high oxidation levels (O/C > 0.82), while aqSOA starts with lower O/C (0.55–0.75). However, prolonged aqueous-phase reactions substantially increase the oxidation state of aqSOA, making its bulk chemical composition closer to that of gasSOA. Additionally, aqueous reactions form a greater abundance of oligomers and high-molecular-weight compounds, alongside a more sustained production of carboxylic acids. AMS spectral signatures representative of phenolic gasSOA have been identified, which, together with tracer ions of aqSOA, can aid in the interpretation of field observation data on aerosol aging within BB smoke. The notable chemical differences between phenolic gasSOA and aqSOA highlighted in this study also underscore the importance of accurately representing both pathways in atmospheric models to better predict the aerosol properties and their environmental impacts.

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

Chemical reaction enhanced graph learning for molecule representation

Abstract Motivation Molecular representation learning (MRL) models molecules with low-dimensional vectors to support biological and chemical applications. Current methods primarily rely on intrinsic molecular information to learn molecular representations, but they often overlook effectively integrating domain knowledge into MRL. Results In this article, we develop a reaction-enhanced graph learning (RXGL) framework for MRL, utilizing chemical reactions as domain knowledge. RXGL introduces dual graph learning modules to model molecule representation. One module employs graph convolutions on molecular graphs to capture molecule structures. The other module constructs a reaction-aware graph from chemical reactions and designs a novel graph attention network on this graph to integrate reaction-level relations into molecular modeling. To refine molecule representations, we design a reaction-based relation learning task, which considers the relations between the reactant and product sides in reactions. In addition, we introduce a cross-view contrastive task to strengthen the cooperative associations between molecular and reaction-aware graph learning. Experiment results show that our RXGL achieves strong performance in various downstream tasks, including product prediction, reaction classification, and molecular property prediction. Availability and implementation The code is publicly available at https://github.com/coder-ACAC/RLM.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Carbon-Carbon Aeroshell Chemical Reaction Model in an Open RPS in Martian Environment

A new RPS is being designed to be open to the vacuum of space to improve its thermal efficiency during deep space missions. This design, however, could lead to some potential complications if it is ever used in the low-pressure gas environment of mission destinations like Mars. This report documents the development of a chemical kinetics reaction model of the FWPF carbon-carbon aeroshell in the Martian gaseous environment. This model uses an “infinite sink” assumption where the RPS is open to the Martian environment with open exchange of flowing gases, so the gas chemistry over the aeroshell does not change with the reaction and represents the worst case scenario of the chemical reaction of the FWPF in the Martian environment. The model combines the kinetics equation of both oxidation reactions that would occur on Mars: the Langmuir-Hinshelwood kinetics of the Boudouard reaction of CO 2 with the chemical kinetics of O 2 oxidation of the carbon-carbon. The FWPF chemical reaction model was benchmarked using the Opila’s experimental chemical reaction test data (1) in a simulated Martian environment to increase the validity of the model’s predicted mass loss rate. The output of the benchmarked chemical reaction model showed a trade-off between the maximum temperature of the aeroshell to have less than 10% reaction at the mission life. For example, an RPS with a five-year intended mission life on Mars has a maximum aeroshell temperature of 835 K / 562 °C, increasing the mission life to ten years reduces the maximum operating temperature to 810 K / 527 °C. If the aeroshell operating temperature in the new RPS on the surface of Mars is 960 Kelvin, the model calculates that the life is limited to 0.22 years / 81 days. Increasing the aeroshell temperature to 1343 Kelvin, further reduces the life to 0.375 days / 9 hrs. This illustrates the thermally activated nature of the chemical kinetics model in limiting the life of the aeroshell.

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

Reaction Rate Sensitivity of the Production of γ-Ray Emitting Isotopes in Core-collapse Supernovae

Radioactive isotopes produced in core-collapse supernovae (CCSNe) provide useful insights into the underlying processes driving the collapse mechanism and the origins of elemental abundances. Their study generates a confluence of major physics research, including experimental measurements of nuclear reaction rates, astrophysical modeling, and γ-ray observations. Here we identify the key nuclear reaction rates to the nucleosynthesis of observable radioactive isotopes in explosive silicon burning during CCSNe. Using the nuclear reaction network calculator SkyNet and current REACLIB reaction rates, we evolve temperature–density–time profiles of the innermost 0.45 M⊙ ejecta from the core collapse and explosion of a 12 M⊙ star. Individually varying 3403 reaction rates by factors of 100, we identify 141 reactions that cause significant differences in the isotopes of interest, namely, 43 K, 47 Ca, 44,47 Sc, 44 Ti, 48,51 Cr, 48,49 V, 52,53 Mn, 55,59 Fe, 56,57 Co, and 56,57,59 Ni. For each of these reactions, we present a novel method to extract the temperature range pertinent to the nucleosynthesis of the relevant isotope; the resulting temperatures lie within the range T = 0.47–6.15 GK. Limiting the variations to within 1σ of STARLIB reaction rate uncertainties further reduces the identified reactions to 48 key rates, which can be used to guide future experimental research. Complete results are presented in tabular form.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Matrix isolation as a tool for studying interstellar chemical reactions

Since the identification of the OH radical as an interstellar species, over 50 molecular species were identified as interstellar denizens. While identification of new species appears straightforward, an explanation for their mechanisms of formation is not. Most astronomers concede that large bodies like interstellar dust grains are necessary for adsorption of molecules and their energies of reactions, but many of the mechanistic steps are unknown and speculative. It is proposed that data from matrix isolation experiments involving the reactions of refractory materials (especially C, Si, and Fe atoms and clusters) with small molecules (mainly H2, H2O, CO, CO2) are particularly applicable to explaining mechanistic details of likely interstellar chemical reactions. In many cases, matrix isolation techniques are the sole method of studying such reactions; also in many cases, complexations and bond rearrangements yield molecules never before observed. The study of these reactions thus provides a logical basis for the mechanisms of interstellar reactions. A list of reactions is presented that would simulate interstellar chemical reactions. These reactions were studied using FTIR-matrix isolation techniques.

Ball, David W.↗
NASA NTRS2007

Calculational and Experimental Investigations of the Pressure Effects on Radical - Radical Cross Combinations Reactions: C2H5 + C2H3

Pressure-dependent product yields have been experimentally determined for the cross-radical reaction C2H5 + C2H3. These results have been extended by calculations. It is shown that the chemically activated combination adduct, 1-C4H8*, is either stabilized by bimolecular collisions or subject to a variety of unimolecular reactions including cyclizations and decompositions. Therefore the "apparent" combination/disproportionation ratio exhibits a complex pressure dependence. The experimental studies were performed at 298 K and at selected pressures between about 4 Torr (0.5 kPa) and 760 Torr (101 kPa). Ethyl and vinyl radicals were simultaneously produced by 193 nm excimer laser photolysis of C2H5COC2H3 or photolysis of C2H3Br and C2H5COC2H5. Gas chromatograph/mass spectrometry/flame ionization detection (GC/MS/FID) were used to identify and quantify the final reaction products. The major combination reactions at pressures between 500 (66.5 kPa) and 760 Torr are (1c) C2H5 + C2H3 yields 1-butene, (2c) C2H5 + C2H5 yields n-butane, and (3c) C2H3 + C2H3 yields 1,3-butadiene. The major products of the disproportionation reactions are ethane, ethylene, and acetylene. At moderate and lower pressures, secondary products, including propene, propane, isobutene, 2-butene (cis and trans), 1-pentene, 1,4-pentadiene, and 1,5-hexadiene are also observed. Two isomers of C4H6, cyclobutene and/or 1,2-butadiene, were also among the likely products. The pressure-dependent yield of the cross-combination product, 1-butene, was compared to the yield of n-butane, the combination product of reaction (2c), which was found to be independent of pressure over the range of this study. The [ 1-C4H8]/[C4H10] ratio was reduced from approx.1.2 at 760 Torr (101 kPa) to approx.0.5 at 100 Torr (13.3 kPa) and approx.0.1 at pressures lower than about 5 Torr (approx.0.7 kPa). Electronic structure and RRKM calculations were used to simulate both unimolecular and bimolecular processes. The relative importance of C-C and C-H bond ruptures, cyclization, decyclization, and complex decompositions are discussed in terms of energetics and structural properties. The pressure dependence of the product yields were computed and dominant reaction paths in this chemically activated system were determined. Both modeling and experiment suggest that the observed pressure dependence of [1-C4H8]/[C4H10] is due to decomposition of the chemically activated combination adduct 1-C4H8* in which the weaker allylic C-C bond is broken: H2C=CHCH2CH3 yields C3H5 + CH3. This reaction occurs even at moderate pressures of approx.200 Torr (26 kPa) and becomes more significant at lower pressures. The additional products detected at lower pressures are formed from secondary radical-radical reactions involving allyl, methyl, ethyl, and vinyl radicals. The modeling studies have extended the predictions of product distributions to different temperatures (200-700 K) and a wider range of pressures (10(exp -3) - 10(exp 5) Torr). These calculations indicate that the high-pressure [1-C4H8]/[C4H10] yield ratio is 1.3 +/- 0.1.

Fahr, Askar↗
DOE OSTI2020

Carbon-Carbon Aeroshell Chemical Reaction Model in an Open RPS in Martian Environment

A new RPS is being designed to be open to the vacuum of space to improve its thermal efficiency during deep space missions. This design, however, could lead to some potential complications if it is ever used in the low-pressure gas environment of mission destinations like Mars. This report documents the development of a chemical kinetics reaction model of the FWPF carbon-carbon aeroshell in the Martian gaseous environment. This model uses an “infinite sink” assumption where the RPS is open to the Martian environment with open exchange of flowing gases, so the gas chemistry over the aeroshell does not change with the reaction and represents the worst case scenario of the chemical reaction of the FWPF in the Martian environment. The model combines the kinetics equation of both oxidation reactions that would occur on Mars: the Langmuir-Hinshelwood kinetics of the Boudouard reaction of CO2 with the chemical kinetics of O2 oxidation of the carbon-carbon. The FWPF chemical reaction model was benchmarked using the Opila’s experimental chemical reaction test data (1) in a simulated Martian environment to increase the validity of the model’s predicted mass loss rate. The output of the benchmarked chemical reaction model showed a trade-off between the maximum temperature of the aeroshell to have less than 10% reaction at the mission life. For example, an RPS with a five-year intended mission life on Mars has a maximum aeroshell temperature of 835 K / 562 °C, increasing the mission life to ten years reduces the maximum operating temperature to 810 K / 527 °C. If the aeroshell operating temperature in the new RPS on the surface of Mars is 960 Kelvin, the model calculates that the life is limited to 0.22 years / 81 days. Increasing the aeroshell temperature to 1343 Kelvin, further reduces the life to 0.375 days / 9 hrs. This illustrates the thermally activated nature of the chemical kinetics model in limiting the life of the aeroshell.

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

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Practical Considerations for Understanding Surface Reaction Mechanisms Involved in Heterogeneous Catalysis

Acquiring useful knowledge about the active site(s) of a catalyst, nature of reactant–catalyst interactions, nature of reactive intermediates, rate-determining step, reaction rate orders that affect various process parameters, and reaction mechanism as a whole is exceedingly challenging. This is especially true in the case of heterogeneous catalysts due to the complexity of the nature of surface active sites and their nonstatic behavior. Here, we present our perspective on differentiating between various surface reaction mechanisms in light of pioneering studies by leaders in the field, with the aim of clarifying some of the confusion associated with these complex mechanisms, especially the Eley–Rideal mechanism. Using bibliometric analysis, we identify and discuss the following four reactions that most commonly invoke the Eley– Rideal mechanism: H 2 activation, CO oxidation, esterification of alcohols by acids, and selective catalytic reduction (SCR) of NO x with NH 3 . Our analysis of studies utilizing well-suited experimental and computational methodologies for differentiating surface reaction mechanisms suggests that the above-mentioned four reactions do not occur via the Eley–Rideal mechanism. Instead, each reaction occurs via the Langmuir–Hinshelwood mechanism with nonidealities present. Lastly, we highlight practical considerations regarding select experimental (characterization methods and differential kinetics) and computational modeling that we believe can provide useful insights to accurately discern between the various possible reaction mechanisms in heterogeneous catalysis.

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

Low-energy measurement of the 25 Mg ⁢(𝛼,𝑛)⁢ 28 Si reaction via neutron spectroscopy

During core helium and carbon burning in massive stars, neutrons are produced mainly by the 22 Ne⁢ (𝛼,𝑛) ⁢25 Mg reaction. Some of these released neutrons are captured by heavy seed nuclei from previous nucleosynthesis events, resulting in the slow production of many of the elements between masses 60 ≤ 𝐴 ≤ 90 via the weak 𝑠 process. Determining the overall neutron flux available in this environment is one of the main challenges in modeling its contributions to heavy element production. Not only must the reaction rate be well defined for the primary neutron source, but the rates of neutron poison and recycling reactions must also be well determined. One critical reaction in the simulation network is 25 Mg ⁢(𝛼,𝑛) ⁢ 28 Si. This reaction, together with 25 Mg ⁢(𝑛,𝛾)⁢ 26 Mg and 25 Mg⁢ (𝛼,𝛾) ⁢ 29 Si, determine how 25 Mg nuclei affect the available neutron flux. Past attempts to measure the 25 Mg ⁢(𝛼,𝑛)⁢ 28 Si cross section using neutron moderator counters have been greatly hindered by lower-𝑍 background reactions. Here, in the present work, neutron spectroscopy with deuterated liquid scintillator detectors has been used. The experimental spectra have been analyzed by applying spectrum unfolding techniques to achieve improved background discrimination for the 25 Mg⁢ (𝛼,𝑛)⁢ 28 Si reaction at low energies, down to 𝐸 𝛼 = 1.75 MeV. The separation of the different background contributions gives further insight into the results of previous moderator counter measurements and the measurements lead to a revised and more reliable determination of the reaction rate.

hydrostatic stellar nucleosynthesis↗
DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Quenching of single-particle strengths in direct reactions

A discrepancy in the asymmetry dependence of spectroscopic factors extracted with different reaction probes calls into question whether the corresponding reaction models are properly understood. In this work, we present extracted spectroscopic factors from the Ar 46 , 34 ( p , d ) Ar 45 , 33 transfer reactions in inverse kinematics at a beam energy of 70 MeV/nucleon. The results are consistent with previous measurements of these reactions at a lower beam energy [Lee et al. , Phys. Rev. Lett. 104 , 112701 (2010) ], indicating that the transfer reaction is a reliable probe for the nuclear structure of exotic nuclei across a wide energy range. Results from a large body of transfer reaction measurements, ( p , p N ) measurements, and theoretical nuclear structure studies make a compelling case for much weaker asymmetry dependence than what is observed with single-nucleon knockout reactions on beryllium or carbon targets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

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

Predicting the potential for mineral scale precipitation in unconventional reservoirs due to fluid-rock and fluid mixing geochemical reactions

Mineral precipitation within hydraulically fractured shale may affect fluid flow pathways and impact longterm hydrocarbon production. The ability to predict geochemical reactions that lead to problematic mineral precipitation will lead to active reservoir management strategies to improve overall production. Using the Marcellus Shale as a case study, a combination of laboratory experiments and reaction path modeling was applied to determine which reactions are likely to occur upon introduction of hydraulic fracturing fluid into the shale reservoir. Experimental results indicate that contact between fracturing fluid and shale will result in dissolution of primary minerals (quartz, feldspars, kaolinite, chlorite, pyrite) and secondary mineral precipitation over time periods of less than one week. Precipitation of barite, Fe-oxides, feldspars, amorphous silica and clay is likely to occur within the reservoir during shut in and early flowback due to mixing between fracturing fluid and reservoir brine as based on modeling saturation indices using experimental fluid data. Reaction path modeling of the experimental scenarios corroborates the specific dissolution and precipitation reactions observed experimentally. Comparison of the results to injected and produced waters from a Marcellus Shale well pad in Greene County, PA, USA, shows that the mineral reactions occur during the hydraulic fracturing, shut in, and early flowback periods. The results presented in this paper demonstrate the value in applying experimental approaches to identify mineral precipitation/dissolution reactions that may significantly impact reservoir performance. The good agreement between geochemical models and experimental results provides confidence that numerical models can be applied to screen the potential fluid-mineral and fluid-mixing reactions in unconventional reservoirs that result in undesired mineral scale precipitation.

58 GEOSCIENCES↗
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 OSTI2021

Unraveling the Intermediate Reaction Complexes and Critical Role of Support-Derived Oxygen Atoms in CO Oxidation on Single-Atom Pt/CeO 2

CeO 2 -supported Pt single-atom catalysts have been extensively studied due to their relevance in automobile emission control and for the fundamental understanding of CeO 2 -based catalysts. Though CeO 2 -supported Pt nanoparticles are often more active than their single-atom counterparts, the former could easily redisperse to Pt single atom under oxidizing diesel conditions. Therefore, to maximize the reactivity of every Pt atom, it is important to fully understand the reaction mechanism of CeO 2 -supported Pt single atoms. Here, we report a CO oxidation study on a Pt/CeO 2 single-atom catalyst, where we can account for all of the neighbors using in situ and operando spectroscopy techniques and microcalorimetric measurements. Coupled with density functional theory calculations, we present a comprehensive picture of the dynamics of the surface species, the role of surface intermediates, and explain the observed reaction kinetics. We started with a catalyst containing exclusively single atoms and used in situ/operando spectroscopy to provide evidence for their stability during the reaction and to identify the Pt 1 complexes before and during the reaction and their binding to CO. The results reveal that in the precatalyst, Pt is present as Pt(O) 4 on the CeO 2 (111) step edge sites, but during CO oxidation, we find that two Pt 1 complexes coexist, representing two states of the same active site in the reaction cycle. The dominant state/complex remains Pt(O) 4 , which adsorbs CO very weakly as shown by CO microcalorimetry. The second, minority state/complex, Pt(CO)(O) 3 is generated through the reaction of Pt(O) 4 with CO, and CO is bound strongly to Pt 1 . Labile oxygen adatoms from the CeO 2 surface play a major role in the regeneration of Pt(O) 4 either directly from Pt(O) 3 or by reaction with the strongly adsorbed CO in Pt(CO)(O) 3 . We show that the formation of an oxygen vacancy and generation of a labile O* are not barrierless, which explains the long lifetime of Pt(CO)(O) 3 and its detectability despite being a minority complex. The results help to develop a comprehensive view of the dynamic evolution of Pt 1 complexes along the reaction cycle and provide mechanistic insights to guide the design of Pt-based single-atom catalysts.

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

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗