Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constant Rate Production”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electronic and Electrochemical Control of Isostructural Ruthenium Hydricities and the Implications for Catalytic Overpotentials

Electronic tuning of metal hydrides enables precise control over potentials, mechanisms, selectivity, and rates of electrocatalytic reactions by regulating bond dissociation free energies such as the hydricity (Δ$G$ H- ° ) and $pK$ a of the catalyst. Here, we investigate a series of electronically tuned ruthenium hydrido complexes that are isostructural at the metal center: [Ru(4,4'-R 2 -bpy) 2 (CO)H] + (R = CF 3 , Cl, H, CH 3 , and CH 3 O; bpy = 2,2'-bipyridine) (denoted as (R)Ru-H+). A substantial 22 kcal mol -1 hydricity range is available across five complexes in three stable oxidation states: (R)Ru-H + , (R)Ru-H 0 , and (R)Ru-H - . Thermodynamic and mechanistic predictions of electrocatalytic proton reduction were tested experimentally by reducing protons from weak acids to H 2 . Two mechanisms are observed, depending on the acid strength and the catalyst hydricity. The rate constants for hydride transfer and protonation of the catalyst were, in some cases, extracted from the analysis of cyclic voltammetry data. A key finding is a 400 mV decrease in the catalytic overpotential for H 2 production by using a doubly reduced electron-poor metal hydride instead of a singly reduced electron-rich metal hydride. In conclusion, the former also exhibits a higher rate constant for hydride transfer, representing a strategy to disconnect rate and free energy relationships.

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

Electrolyte-Assisted Hydrogen Cycling in Lithium and Sodium Alanates at Low Pressures and Temperatures

An investigation of electrolyte-assisted hydrogen storage reactions in complex aluminum hydrides (LiAlH4 and NaAlH4) reveals significantly reduced reaction times for hydrogen desorption and uptake in the presence of an electrolyte. LiAlH4 evolves ~7.8 wt% H2 over ~3 h in the presence of a Li-KBH4 eutectic at 130 °C compared to ~25 h for the same material without the electrolyte. Similarly, NaAlH4 exhibits 4.8 wt% H2 evolution over ~4 h in the presence of a diglyme electrolyte at 150 °C compared to 4.4 wt% in ~15 h for the same material without the electrolyte. These reduced reaction times are composed of two effects, an increase in reaction rates and a change in the reaction kinetics. While typical solid state dehydrogenation reactions exhibit kinetics with rates that continuously decrease with the extent of reaction, we find that the addition of an electrolyte results in rates that are relatively constant over the full desorption window. Fitting the kinetics to an Avrami-Erofe’ev model supports these observations. The desorption rate coefficients increase in the presence of an electrolyte, suggesting an increase in the velocities of the reactant-product interfaces. In addition, including an electrolyte increases the growth parameters, primarily for the second desorption steps, resulting in the observed relatively constant reaction rates. Similar effects occur upon hydrogen uptake in NaH/Al where the presence of an electrolyte enables hydrogenation under more practical low temperature (75 °C) and pressure (50 bar H2) conditions.

25 ENERGY STORAGE↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

Thermal Stability of a Eutectic Mixture of Bis(2,2-dinitropropyl) Acetal and Formal: Part C. Kinetic Compensation Effect

Here, the aging behavior of a eutectic mixture of bis(2,2-dinitropropyl) acetal and formal [called NP here] has been studied in various atmospheres [dry (air or nitrogen) versus wet] at temperatures 70 °C and below. The properties of aged samples were analyzed using Fourier transform infrared (FTIR) spectroscopy, Karl Fischer (KF) titration, liquid chromatography/mass spectrometry (LC/MS), and thermogravimetric analysis (TGA) over a period of three years. The results indicate that at aging temperatures up to 55 °C, the initial rates of water production from nitrous acid (HONO) formation and decomposition into the water, NO, and NO 2 follows a 1st order rate law and the rate constants follow an Arrhenius law as a function of temperature. The activation energies and pre-factors for water and volatiles production yield a single linear kinetic compensation plot, suggesting a common degradation pathway between NP and the various combinations of its constituents. Within a narrow temperature range, around 55 °C, a trace amount of water in NP stabilizes its properties by preventing HONO elimination. When the aging temperature is substantially higher than 55 °C, the nature of the degradation mechanism changes. It is suspected that the degradation products of NO x , water, and HNO 3 serve as catalysts to auto-catalyze (kinetics beyond the 1st order) and further degrade NP. The effect of headspace volume on this auto-catalytic process will be discussed.

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

Sensitivity and history match analysis of a carbon dioxide “huff-and-puff” injection test in a horizontal shale gas well in Tennessee

Due to improvements in well development, shale gas production has gained importance recently, especially in the United States. To improve gas production and develop a better understanding of shale reservoirs, researchers are conducting field tests to monitor CO 2 storage and enhanced gas recovery. Reservoir simulations can then utilize these field results for sensitivity studies, uncertainty analysis, history matching, gas production forecasts, and CO 2 storage capacity estimations. One such field test was performed for the Chattanooga Shale formation in Morgan County, Tennessee. Approximately 463 tonnes (510 tons) of CO 2 was successfully injected into a hydraulically fractured horizontal well over a thirteen-day period in March 2014. The injection test was achieved in four steps: pre-injection, injection, soaking, and flowback. In this paper, those steps were modeled with a reservoir simulator to match the historic production and injection rates with gas component compositions and forecast the production for five years. The reservoir fluid was modeled as a multi-component gas, including CH 4 , C 2 H 6 , C 3 H 8 , N 2 , and CO 2 . Langmuir constants, reservoir pressure, and fracture network volume were adjusted to match simulation results with observed production rates and gas composition. Results showed that, under the same reservoir conditions, each gas component behaves differently by way of compositional changes in production. CH4 behaves like N2, while CO 2 behaves like heavier hydrocarbons such as C 2 H 6 and C 3 H 8 . CO 2 plume results showed that, after injection, the produced CO 2 is mostly derived from fractured limestone (i.e., Fort Payne Formation) because the injected CO 2 is not adsorbed in limestone but rather by shale formations.

03 NATURAL GAS↗
DOE OSTI2023

Dihydroxyterephthalate—A Trojan Horse PET Counit for Facile Chemical Recycling

Abstract Here, low‐energy poly(ethylene terephthalate) (PET) chemical recycling in water: PET copolymers with diethyl 2,5‐dihydroxyterephthalate (DHTE) undergo selective hydrolysis at DHTE sites, autocatalyzed by neighboring group participation, is demonstrated. Liberated oligomeric subchains further hydrolyze until only small molecules remain. Poly(ethylene terephthalate‐ stat ‐2,5‐dihydroxyterephthalate) copolymers were synthesized via melt polycondensation and then hydrolyzed in 150–200 °C water with 0–1 wt% ZnCl 2 , or alternatively in simulated sea water. Degradation progress follows pseudo‐first order kinetics. With increasing DHTE loading, the rate constant increases monotonically while the thermal activation barrier decreases. The depolymerization products are ethylene glycol, terephthalic acid, 2,5‐dihydroxyterephthalic acid, and bis(2‐hydroxyethyl) terephthalate dimer, which could be used to regenerate virgin polymer. Composition‐optimized copolymers show a decrease of nearly 50% in the Arrhenius activation energy, suggesting a 6‐order reduction in depolymerization time under ambient conditions compared to that of PET homopolymer. This study provides new insight to the design of polymers for end‐of‐life while maintaining key properties like service temperature and mechanical properties. Moreover, this chemical recycling procedure is more environmentally friendly compared to traditional approaches since water is the only needed material, which is green, sustainable, and cheap.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hydrophobic pockets built in polymer micelles enhance the reactivity of Cu 2+ ions

We report the hydrophobicity-enhanced reactivity of Cu2+ ions as an ester hydrolase. Using a dipicolylamine (DPA) containing reversible addition–fragmentation chain transfer agent, the synthetic sequence, either hydrophobic or hydrophilic first in amphiphilic block copolymers of polystyrene-block-poly(N'N-dimethylacrylamide) (PS-b-PDMA), was varied to control the location of the binding motif, DPA, in the hydrophobic core or on the hydrated corona of polymer micelles. The hydrophobicity of Cu 2+ sites showed a significant impact (as large as 60 times more activity) on their catalytic efficiency towards ester hydrolase. With two different kinetic modes, including Michaelis–Menten and the reverse saturation kinetics models, the binding constant Kb of the substrates to Cu 2+ sites were quantitatively analyzed and we demonstrate that hydrophobicity favors the binding of the substrates to Cu2+ sites at polymer micelles with smaller sizes, however, K b decays exponentially with micellar diameters. Despite the diffusion barrier, hydrophobicity shows a profound impact on the catalytic rate constant kc that measures the single conversion rate of bound substrates to products. There is a 16–20 times kinetic enhancement in the hydrolase activity, completely endowed by the hydrophobic microenvironment of Cu 2+ sites compared to micelles with similar sizes. Our results indicate how the hydrophobicity of Cu 2+ -containing micelles can impact the catalytic efficiency and potentially illustrate a promising way toward the design of bioinspired catalysts.

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

Effects of moisture and aging upon decomposition of methyl iodide by reduced silver mordenite

Abstract Reduced silver mordenite has been considered as a sorbent for the capture of organic iodides, especially methyl iodide, from off-gases produced by aqueous used nuclear fuel reprocessing operations. The adsorption capacity of this material has been unpredictable especially when NO x and water are present. Previous work has found that a catalytic decomposition reaction is occurring on the surface but few determinations have been made of the kinetics of this reaction. The work presented tested the adsorption behavior and apparent catalytic reaction rate in humid conditions and compared those to dry conditions testing. Both experiments observed a first order reaction with rate constants of 0.0847 L/g sorbent/s and 0.1202 L/g sorbent/s respectively. Such a reduction in apparent rate constant is possibly due to either water obstructing methyl iodide adsorption or product desorption limitation. Changes in the adsorption profile were also apparent between these two, with the humid conditions experiment reaching saturation sooner than the dry conditions experiment. Additionally, an experiment into the effects of sorbent storage in a controlled laboratory environment was performed. The performance of the sorbent materials that were stored with silver in the zerovalent state was slightly inferior to those materials that were stored in ionic form (Ag + ) and reduced to zerovalent silver immediately prior to subjecting them to sorption test. The materials stored with silver in the ionic form (and reduced just prior to application) behaved essentially similarly to the freshly synthesized (and reduced) sorbents in the sorption tests. This suggests that zerovalent silver experiences some oxidation resulting in deactivation of some sites.

Goettsche, Heinrik↗
DOE OSTI2024

Bioelectrochemical reactor to manage anthropogenic sulfate pollution for freshwater ecosystems: Mathematical modeling and experimental validation

Anthropogenic sulfate loading into otherwise low-sulfate freshwater systems can cause significant ecological consequences as a biogeochemical stressor. To address this challenge, in situ bioremediation technologies have been developed to leverage naturally occurring microorganisms that transform sulfate into sulfide rather than implementing resource-intensive physio-chemical processes. However, bioremediation technologies often require the supply of electron donors to facilitate biological sulfate reduction. Bioelectrochemical systems (BES) can be an alternative approach for supplying molecular hydrogen as an electron donor for sulfate-reducing bacteria through water electrolysis. Although the fundamental mechanisms behind BESs have been studied, limited research has evaluated the design and operational parameters of treatment systems when developing BESs on a scale relevant to environmental systems. Further, this study aimed to develop an application-based mathematical model to evaluate the performance of BESs across a range of reactor configurations and operational modes. The model was based on sulfate transformation by hydrogenotrophic sulfate-reducing bacteria coupled with the recovery of solid iron sulfide species formed by the oxidative dissolution of dissolved ferrous iron from a stainless steel anode. Sulfate removal closely corresponded to the rate of electrolytic hydrogen production and hydraulic residence time but was less sensitive to specific microbial rate constants. The mathematical model results were compared to experimental data from a pilot-scale BES tested with nonacidic mine drainage as a case study. The close agreement between the mathematical model and the pilot-scale BES experiment highlights the efficacy of using a mathematical model as a tool to develop a conceptual design of a scaled-up treatment system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Re-Examination of the N 2 O + O Reaction

The reaction of N 2 O with O is a key step in consumption of nitrous oxide in thermal processes. It has two product channels, NO + NO (R2) and N 2 + O 2 (R3). The rate constant for R2 has been measured both in the forward and the reverse direction at elevated temperature and is well established. However, the rate constant for the N 2 + O 2 channel (R3) has been difficult to quantify and has significant error limits. Here, the direct reaction on the triplet surface has a barrier of around 40 kcal mol –1 , and it is too slow for the N 2 + O 2 channel to have any practical significance. Recently, Pham and Lin (2022) suggested an alternative low activation energy reaction path that involves intersystem crossing and reaction on the singlet surface. In the present work, we re-examined a wide range of experiments relevant for the N 2 O + O reaction through kinetic modeling, paying attention to the impact of artifacts such as impurities and surface reactions. Experimental results from shock tubes and batch reactors on the final NO yield in N 2 O decomposition, covering temperatures of 973–2200 K and pressures of 0.013–11.5 atm, support k 3 ~ 0, consistent with the high activation energy for reaction on the triplet surface and a low probability of ISC.

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

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

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↗
DOE OSTI2020

Valuation and cost reduction of behind-the-meter hydrogen production in Hawaii

A 250kW hydrogen electrolysis facility was recently installed at the Natural Energy Laboratory of Hawaii Authority's (NELHA's) campus. This facility that will begin operation in 2020 to produce hydrogen for fuel cell buses on the island to demonstrate of the application of hydrogen to decarbonize transportation. Given the size of the electrolysis station, it has the potential to significantly increase electricity costs for the campus, which is subject to energy and peak demand charges from the local utility. In this paper, we analyze the cost of hydrogen production at NELHA given the rate structure options available from the utility. Production costs are estimated using optimal versus constant scheduling of the facility to meet the buses’ demand. A model of the electrolysis station is used to capture changes in production efficiency over the power range in the optimization routine. The effects of combining the station and campus load versus standalone operation and increasing solar generation are also explored. The analyses surrounding this scenario show the importance of multiple factors on the potential profitability of hydrogen production in behind-the-meter applications and show trends that could have implications for other similar installations.

08 HYDROGEN↗
DOE OSTI2021

Temperature controlled friction extrusion process

Processing temperature is decisive to the properties and performance of the product in friction extrusion. However, the constantly varying thermal boundary conditions result in an unstable processing temperature even fixed die rotational speed and feed rate were implemented, which extensively disturb the stability of the process and the quality of the product. It is key to maintain the process temperature to ensure the properties of the product is uniform. What's more, a controlled processing temperature enables the properties of the product to be manipulatable between parts and between different locations of one part

Li, Xiao↗
DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

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

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗