Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mechanistic and kinetic relevance of hydrogen and water in CO 2 hydrogenation on Cu-based catalysts

Here, we ally steady-state kinetics, kinetic isotope effects, and density functional theory (DFT) calculations to illustrate that Cu-based catalysts remain saturated by H-adatoms (H*) and molecular formic acid (HCOOH**) during CO 2 hydrogenation. High H* coverage under methanol synthesis conditions is evidenced by reverse water-gas shift (RWGS) rates that exhibit positive H 2 reaction orders only at P H2 ≲ 0.5 bar, above which methanol synthesis and RWGS rates exhibit first and zeroth order dependence on P H2 , respectively. HCOOH** also accumulates on the surface with increasing P CO2 as informed by the Langmuir-type dependence on P CO2 (0.25-23 bar) for both methanol synthesis and RWGS. As both HCOOH** and H* have one H-atom per site occupied, the two species share the same P H2 dependence and give rise to CO 2 reaction orders that are independent of P H2 . Surface coverages determined based on kinetic analyses are further corroborated with DFT-derived adsorption energies that show favorable HCOOH** adsorbate-adsorbate interactions as well as repulsive interactions for bidentate formate (HCOO**) on H*-saturated surfaces. Methanol selectivity remains invariant with P CO2 and P CO despite CO inhibiting reaction rates, thereby demonstrating methanol synthesis and RWGS occur on the same active site. In contrast, water preferentially inhibits methanol synthesis rates, increases methanol synthesis H 2 reaction order from 1.0 to 1.5, and alters the methanol synthesis H 2 /D 2 kinetic isotope effect; the inhibitory effect of H 2 O thus cannot be attributed to competitive adsorption alone and instead reflects a change in the rate-determining step for methanol synthesis. The disparate kinetics of methanol synthesis and RWGS evince a branching pathway where methanol is formed from formates and CO is formed from carboxylates. The presented work thus identifies the relevant surface species, underscores the distinct catalytic role of water in branching methanol synthesis and RWGS pathways, and, in doing so, details a mechanistic picture that yields predictable rates and reaction orders for both methanol synthesis and RWGS on Cu-based CO 2 hydrogenation catalysts.

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

Kinetics of the reaction of ferrous ions with hydroxyl radicals in the temperature range 25–300 °C

Here, the kinetics and mechanism of the reaction between OH radicals and ferrous ions in the temperature range 25–300 °C were studied using pulse radiolysis. At temperatures <150 °C the rate of reaction is essentially independent of temperature, while at temperatures >150 °C the activation energy is 45.8 ± 3.0 kJ mol —1 . The change in activation energy is attributed to a change in the dominant mechanism from hydrogen atom transfer (HAT) to dissociative ligand interchange. The kinetic isotope effect (KIE) was measured by repeating experiments in heavy water. A value of 2.9 was measured at room temperature where HAT is the dominant mechanism. The KIE decreases to zero at temperatures > 150 °C as ligand interchange becomes dominant and the O–H bond is no longer involved in the reaction.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Mass transfer coefficients, kL, and air-CO2 ingassing rates in 3.4 m2 and 1-acre raceway ponds.

The US DOE 2016 Billion Ton Update (Langholtz 2016) projected that CO2 in flue-gas from power plants or similar sources would limit U.S. algal biofuel potential to under 5 billion gallons gasoline equivalent (gge) per year, due to lack of sufficient land, water and other requirements near such flue gas sources. 2nd generation carbon capture technologies are proposed to overcome this constraint (Davis 2017), and could expand this resource potential nearly 10-fold, assuming CO2 costs near the future, 2025 DOE NETL, goal of $40/metric ton CO2. An alternative pathway is the direct uptake of air CO2 into the algal ponds through an increase in the air-CO2 transfer rate resulting from the reaction of CO2 with hydroxide ions. However, such a ‘chemical enhancement’ in mass transfer depends on a high pH in the culture, and, by extension, the ability of microalgae to maintain high rates of carbon assimilation under such conditions. Abiotic experiments characterized the air-CO2 mass-transfer rate in 1-acre and 3.4 m2 raceway ponds as a function of pH, with rates approaching 10 g C/m2-day at pH 12 in the brackish water tested, equivalent to a biomass productivity potential of near 20 g AFDW/m2-day (@ 0.47 g C/g AFDW). The mass transfer rate was found to be nearly independent of the mass-transfer coefficient, indicating that at turbulence levels achievable in typical raceway ponds, the system is reaction rate, rather than diffusion limited. In biotic trials with an unknown green microalga, biomass productivity in ponds receiving CO2 from only direct surface air-CO2 exchange averaged 5 +/- 0.5 g AFDW/m2-day in August - September central California conditions, about half-that of experimental controls fertilized with supplemental CO2. pH in the air-only treatment reached a maximum of 10.5, supporting a model predicted air-CO2 ingassing rate between 2-3 g C/m2-day, consistent with the observed biomass carbon assimilation rate. Results suggest that unique alkaliphilic strains are required if reliant on direct in-pond air-CO2 transfer for algal inorganic carbon supply, and that reaching an economically viable biomass productivity will require strains that thrive at pH 11 and above.

09 BIOMASS FUELS↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

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

Solution and Interface Structure and Dynamics in Geochemistry: Gateway to Link Elementary Processes to Mineral Nucleation and Growth

Our understanding of the process of mineral nucleation is starting to evolve from the solely thermodynamical aspects in classical nucleation theory toward a mechanistic study of reactions influenced by solvation structures of ions and their dynamics. A knowledge of how these atomic and molecular interactions respond to changing solution composition and/or the nature of an interface will form a fundamental foundation to understand rate-limiting reactions and processes for nucleation and growth. Here, we review recent advances from multidisciplinary experimental and computational approaches of the structure and dynamics of aqueous solutions and mineral–liquid interfaces that will influence rates of nucleation (and growth behavior) both homogeneously in aqueous solutions and heterogeneously at mineral–liquid interfaces. Furthermore, these processes include structure (and extent) of solute clustering, solvent degrees of freedom, and the presence of noninnocent charge-balancing ions, all of which can drive the selectivity and exert control over nucleation mechanisms. Similar controlling reactions can also occur at interfaces with the added complexities of how solute ions interact with and assemble on surface sites and use different mechanisms from bulk solutes, such as lateral solvent exchanges between distinct surface sites. To gain a predictive understanding of nucleation chemistry, we emphasize several key knowledge gaps that will need to be addressed in this field and highlight the importance of combining computations and experiments.

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

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

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

Effect of a redox-mediating ligand shell on photocatalysis by CdS quantum dots

Semiconductor quantum dots (QDs) are efficient organic photoredox catalysts due to their high extinction coefficients and easily tunable band edge potentials. Despite the majority of the surface being covered by ligands, our understanding of the effect of the ligand shell on organic photocatalysis is limited to steric effects. We hypothesize that we can increase the activity of QD photocatalysts by designing a ligand shell with targeted electronic properties, namely, redox-mediating ligands. Herein, we functionalize our QDs with hole-mediating ferrocene (Fc) derivative ligands and perform a reaction where the slow step is hole transfer from QD to substrate. Surprisingly, we find that a hole-shuttling Fc inhibits catalysis, but confers much greater stability to the catalyst by preventing a build-up of destructive holes. We also find that dynamically bound Fc ligands can promote catalysis by surface exchange and creation of a more permeable ligand shell. Finally, we find that trapping the electron on a ligand dramatically increases the rate of reaction. Furthermore, these results have major implications for understanding the rate-limiting processes for charge transfer from QDs and the role of the ligand shell in modulating it.

Chemistry↗
DOE OSTI2023

Stochastic optimization of a uranium oxide reaction mechanism using plasma flow reactor measurements

Abstract In this work, a coupled Monte Carlo Genetic Algorithm (MCGA) approach is used to optimize a gas phase uranium oxide reaction mechanism based on plasma flow reactor (PFR) measurements. The PFR produces a steady Ar plasma containing U, O, H, and N species with high temperature regions (3000–5000 K) relevant to observing UO formation via optical emission spectroscopy. A global kinetic treatment is used to model the chemical evolution in the PFR and to produce synthetic emission signals for direct comparison with experiments. The parameter space of a uranium oxide reaction mechanism is then explored via Monte Carlo sampling using objective functions to quantify the model-experiment agreement. The Monte Carlo results are subsequently refined using a genetic algorithm to obtain an experimentally corroborated set of reaction pathways and rate coefficients. Out of 12 reaction channels targeted for optimization, four channels are found to be well constrained across all optimization runs while another three channels are constrained in select cases. The optimized channels highlight the importance of the OH radical in oxidizing uranium in the PFR. This study comprises a first step toward producing a comprehensive experimentally validated reaction mechanism for gas phase uranium molecular species formation.

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

Heat of Mixing During Fast Charge/Discharge of a Li-Ion Cell: A Study on NMC523 Cathode

Predicting temperature rise accurately during fast charge/discharge of a Li-ion cell is essential to avoid thermal runaway and extend battery life. While modeling the temperature rise, it is necessary to model the heat generation correctly. Heat of mixing, one of the four main sources of heat generation in a Li-ion battery, has often been considered insignificant and therefore excluded from modeling. When included, it is modeled using the expression from a Taylor expansion approximation. In this work, we have shown the conditions when including the heat of mixing becomes important and quantified the error associated with using the Taylor expansion, especially under high charge/discharge conditions. Consequently, we carry out the calculation of the rate of enthalpy change and the heat generation rate from the most fundamental equation for the rate of the total enthalpy change without simplifying assumptions or approximations. The heat generation rate calculated doing so naturally includes irreversible, reversible and mixing heat. We then exclusively separate out the heat of mixing by subtracting the rate of enthalpy change by reaction from the rate of the total enthalpy change. Results show that the contribution of heat of mixing in the total heat generated increases with the charge/discharge rate and is as large as 23% for a 6 C discharge. This result suggests that while modeling heat generation for fast charge/discharge, it is necessary to include the heat of mixing and avoid calculating it using the Taylor expansion approximation.

25 ENERGY STORAGE↗
DOE OSTI2024

Non-equilibrium rate theory for polariton relaxation dynamics

We derive an analytic expression of the non-equilibrium Fermi’s golden rule (NE-FGR) expression for a Holstein–Tavis–Cumming Hamiltonian, a universal model for many molecules collectively coupled to the optical cavity. These NE-FGR expressions capture the full-time-dependent behavior of the rate constant for transitions from polariton states to dark states. The rate is shown to be reduced to the well-known frequency domain-based equilibrium Fermi’s golden rule (E-FGR) expression in the equilibrium and collective limit and is shown to retain the same scaling with the number of sites in non-equilibrium and non-collective cases. We use these NE-FGR to perform population dynamics with a time-non-local and time-local quantum master equation and obtain accurate population dynamics from the initially occupied upper or lower polariton states. Furthermore, NE-FGR significantly improves the accuracy of the population dynamics when starting from the lower polariton compared to the E-FGR theory, highlighting the importance of the non-Markovian behavior and the short-time transient behavior in the transition rate constant.

Chemical dynamics↗
DOE OSTI2023

Ligand effect on switching the rate-determining step of water oxidation in atomically precise metal nanoclusters

The ligand effects of atomically precise metal nanoclusters on electrocatalysis kinetics have been rarely revealed. Herein, we employ atomically precise Au 25 nanoclusters with different ligands (i.e., para-mercaptobenzoic acid, 6-mercaptohexanoic acid, and homocysteine) as paradigm electrocatalysts to demonstrate oxygen evolution reaction rate-determining step switching through ligand engineering. Au 25 nanoclusters capped by para-mercaptobenzoic acid exhibit a better performance with nearly 4 times higher than that of Au 25 NCs capped by other two ligands. We deduce that para-mercaptobenzoic acid with a stronger electron-withdrawing ability establishes more partial positive charges on Au(I) (i.e., active sites) for facilitating feasible adsorption of OH – in alkaline media. X-ray photo-electron spectroscopy and theoretical study indicate a profound electron transfer from Au(I) to para-mercaptobenzoic acid. The Tafel slope and in situ Raman spectroscopy suggest different ligands trigger different rate-determining step for these Au 25 nanoclusters. The mechanistic insights reported here can add to the acceptance of atomically precise metal nanoclusters as effective electrocatalysts.

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

Mechanistic Analysis of Metallaphotoredox C–N Coupling: Photocatalysis Initiates and Perpetuates Ni(I)/Ni(III) Coupling Activity

In this work, the combined use of reaction kinetic analysis, ultrafast spectroscopy, and stoichiometric organometallic studies has enabled the elucidation of the mechanistic underpinnings to a photocatalytic C–N cross-coupling reaction. Steady-state and ultrafast spectroscopic techniques were used to track the excited-state evolution of the employed iridium photocatalyst, determine the resting states of both iridium and nickel catalysts, and uncover the photochemical mechanism for reductive activation of the nickel cocatalyst. Stoichiometric organometallic studies along with a comprehensive kinetic study of the reaction, including rate–driving force analysis, unveiled the crucial role of photocatalysis in both initiating and sustaining a Ni(I)/Ni(III) cross-coupling mechanism. The insights gleaned from this study further enabled the discovery of a new photocatalyst providing a >30-fold rate increase.

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

Mechanistic understanding of pH effects on the oxygen evolution reaction

The oxygen-evolution reaction (OER) is pivotal in many energy-conversion technologies as it is an important counter reaction to others that convert stable chemicals to higher-value products using electrochemistry. The local microenvironment and pH for the anode OER can vary from acidic to neutral to alkaline depending on the system being explored, making definitive mechanistic insights difficult. In this paper, we couple experiments, first-principles calculations based on density functional theory, microkinetics, and transport modeling to explore the entire pH range of the OER. At low current densities, neutral pH values unexpectedly perform better than the acidic and alkaline conditions, and this trend is reversed at higher current densities (> 20 mA cm -2 ). Using multiscale modeling, this switch is rationalized by a change from a dual-reaction mechanism to a single rate-determining step. The model also shows how the alkaline reaction rates dominate in the middle to high pH range. Furthermore, we explore that the local pH for near-neutral conditions is much different (e.g., 2.4 at the reaction surface vs. 9 in the bulk) than the pH extremes, demonstrating the criticality that transport phenomena plays in kinetic activity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Atom-at-a-Time Radioactive Molecule Identification: Looking toward Studies of Superheavy Elements

The chemical behavior of superheavy elements (SHEs, Z > 103) remains poorly understood. Their chemical properties are expected to deviate from established trends, challenging the predictive power of the periodic table. To investigate these elements experimentally, they must first be synthesized through nuclear reactions and then quickly subjected to chemical studies before they decay. Given the low production rates of these reactions and the need for measurements on an atom-at-a-time basis, innovative techniques are needed. Here, to address these challenges, a novel gas-phase chemistry method has been developed at Lawrence Berkeley National Laboratory, utilizing the Berkeley Gas-filled Separator and FIONA. This technique enables the production, identification, and study of molecular species formed by SHEs. As a proof of concept, we present measurements on the formation and identification of 151,152 HoO + molecules, demonstrating the capability to study the production of radioactive molecules under controlled conditions and directly identify them via their mass-to-charge ratio. These measurements validate the effectiveness of this technique for low-statistics SHE studies, highlighting the potential of this approach to ignite the next generation of experimental SHE chemistry research, offering a path to re-evaluate SHE placement on the periodic table.

Chemistry↗
DOE OSTI2026

Activating magnetite ores for aqueous ironmaking at high current densities

Low-temperature electrochemical cells reducing iron oxides to metal in alkaline electrolytes can support fully electrified steelmaking processes. Previous studies on these cells have primarily focused on high-surface-area hematite, Fe 2 O 3 , reactants whereas attempts to reduce suspensions of magnetite, Fe 3 O 4 —one of the two feedstocks for existing ironmaking reactors—have generally been limited to low rates of reaction (<30 mA cm −2 ). Here, in this study, we control the crystalline domain size of Fe 2 O 3 and Fe 3 O 4 particles in 10 M NaOH electrolytes to study how the nanoscale morphology of oxides controls the rate of electrochemical ironmaking. Rotating-ring disk electrode measurements of Fe 2+ , in situ Raman spectroscopy of the electrode surface, and ex situ electron microscopy were consistent with a hypothesized passivation process at Fe 3 O 4 surfaces that may prevent the continuous formation of soluble intermediates. Sufficiently small (<100 nm diameter) oxide particles yielded Fe partial current densities >160 mA cm −2 , a fivefold increase relative to previously reported rates for Fe 3 O 4 suspensions and comparable to active Fe 2 O 3 . Electron microscopy revealed that electrodeposited films were composed of micron-scale crystalline Fe domains with a porous film of Fe 3 O 4 nanoparticles and supports a model where Fe is grown primarily from soluble Fe 2+ intermediates. Based on these insights, inactive blast-furnace-grade iron-oxides were transformed into high surface area nanoparticles (1.6 to 229.2 m 2 g −1 ) via reprecipitation, leading to a ninefold enhancement in faradaic efficiency and an Fe partial current density of 120 mA cm −2 . When integrated with chlor-iron cells producing reagents for reprecipitation, this approach could lead to a cost-competitive process for electrochemical ironmaking from industrially relevant feedstocks.

TEM↗
DOE OSTI2020

Solvent influence on non-adiabatic interfacial electron transfer at conductive oxide electrolyte interfaces

The kinetics for interfacial electron transfer (ET) from a transparent conductive oxide (tin-doped indium oxide, ITO, Sn:In 2 O 3 ) to molecular acceptors 4-[N,N-di(p-tolyl)amino]benzylphosphonic acid, TPA, and [Ru II (bpy) 2 (4,4' -(PO 3 H 2 ) 2 -bpy)] 2+ , RuP, positioned at variable distances within and beyond the electric double layer (EDL), were quantified in benzonitrile and methanol by nanosecond absorption spectroscopy as a function of the thermodynamic driving force, -ΔG ° . Relevant ET parameters such as the rate constant, k et , reorganization energy, λ, and electronic coupling, Hab, were extracted from the kinetic data. Overall, k et increased as the distance between the molecular acceptor and the conductor decreased. For redox active molecules within the Helmholtz planes of the EDL, k et was nearly independent of -ΔG° , consistent with a negligibly small λ value. Rips–Jortner analysis revealed a non-adiabatic electron transfer mechanism consistent with H ab < 1 cm -1 . The data indicate that the barrier for electron transfer is greatly diminished at the conductor–electrolyte interface.

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

Reorganization Energies for Interfacial Proton-Coupled Electron Transfer to a Water Oxidation Catalyst

The reorganization energy (λ) for interfacial electron transfer (ET) and proton-coupled electron transfer (PCET) from a conductive metal oxide (In 2 O 3 :Sn, ITO) to a surface-bound water oxidation catalyst were extracted from kinetic data measured as a function of the thermodynamic driving force. Visible light excitation resulted in rapid excited-state injection (k inj > 10 8 s -1 ) to the ITO, which photo-initiated the two interfacial reactions of interest. The rate constants for both reactions increased with the driving force, −ΔG°, to a saturating limit, k max , with rate constants consistently larger for ET than for PCET. Marcus-Gerischer analysis of the kinetic data provided the reorganization energy for interfacial PCET (0.90 ± 0.02 eV) and ET (0.40 ± 0.02 eV), respectively. The magnitude of k max for PCET was found to decrease with pH, behavior that was absent for ET. Both the decrease in k max and the larger reorganization energy for an unwanted competing PCET reaction from the ITO to the oxidized catalyst showcases a significant kinetic advantage for driving solar water oxidation at high pH. Computational analysis revealed a larger inner-sphere reorganization energy contribution for PCET than for ET arising from a more significant change in the Ru–O bond length for PCET. Extending Marcus-Gerischer theory to PCET by including the excited electron-proton vibronic states and the proton donor-acceptor motion provided an apparent reorganization energy of 1.01 eV. Furthermore, this study demonstrates that the Marcus-Gerischer theory initially developed for ET can be reliably extended to PCET for quantifying and interpreting reorganization energies observed experimentally.

Catalysts Kinetic parameters↗
DOE OSTI2024

Implementation of new mixture rules has a substantial impact on combustion predictions for H 2 and NH 3

Complex-forming reactions comprise a substantial fraction of all important combustion reactions and are central to combustion behavior. Despite being often called “pressure-dependent” reactions, their rate constants depend on not only the pressure but also the composition. While modern combustion codes allow arbitrarily high accuracy in treating pressure dependence, recent work has consistently demonstrated dramatic failures of essentially all available treatments of mixture dependence. In situations where mixture dependence is treated at all, it is inevitably treated through specification of pressure-dependent rate constants for a set of pure bath gases, which are then combined to estimate the rate constant in a mixture via a “mixture rule.” While there had been a generally unquestioning confidence in these mixture rules, they had, in reality, been scarcely tested until the last decade, when comparisons against master equation calculations revealed order-of-magnitude errors for important pressure-dependent reactions. New mixture rules, based on the reduced pressure, have recently been proposed and shown to reproduce master equation calculations for broad classes of complex-forming reactions very accurately. Here, in this work, we present an implementation of one such new mixture rule (“LMR-R”) in Cantera and then use it to enable simulations that use new high-accuracy ab initio data for individual bath gases (for the first time, since codes previously could not accommodate the complex bath gas dependence). Demonstrations focus on combustion of H 2 and NH 3 , where (1) high-accuracy ab initio data are available and (2) the impact is expected to be large due to the high fractions of efficient colliders (e.g., H 2 O and NH 3 ) in the burned and unburned gases. Indeed, we find the impact of this treatment to be substantial and may explain previous modeling difficulties for these important carbon-free fuels, particularly for NH 3 , whose extraordinarily high third-body efficiency (~20) is often omitted from kinetic models.

Ammonia↗