Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetics”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Effect of NO on DME-Methanol HCCI Combustion Using a Reduced Chemical Kinetics Mechanism

Methanol is an attractive fuel for the maritime sector due to its wide availability. Its direct use as a fuel, however, is accompanied by challenges such as high latent heat of vaporization and low cetane number. A potential solution to overcome the ignition properties of methanol could be through on-board generation of dimethyl ether (DME) via catalytic dehydration of methanol. The resulting mixture from dehydration can be mixed in with the intake air to generate a homogenous charge compression ignition (HCCI) preburn for subsequent direct injection (DI) mixing controlled compression ignition (MCCI) of neat methanol. Within that context, complementary experimental work found that the influence of combustion residuals on the heat release rate (HRR) was significant, specifically for residual NO. This finding motivated the present computational and kinetic evaluation of the effects of NO on the low (LTHR) and high (HTHR) temperature heat release rates. The strong influence of small quantities of NO on the combustion process of a DME/methanol/H2O mixture (low catalyst or reactor efficiency) necessitated a kinetics-based investigation into this phenomenon. A mechanism sourced from the existing literature with NO had 172 species and 1375 reactions, making it computationally expensive for use. Hence, a mechanism reduction effort was implemented, and a rate constant (k) tuning effort based on sensitivity analysis was needed to validate experimental results using a zero-dimensional engine model in Cantera. The reduced mechanism was able to successfully capture the negligible influence of NO addition on DME HCCI combustion, whereas an advancement in LTHR and HTHR for a DME/methanol/H2O mixture was kinetically confirmed. Reaction pathway analysis showed that addition of NO chemically counteracted the OH sink created by alcohols like methanol, increasing the effectiveness of DME ignition.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗
DOE OSTI2021

Evaluating the benefits of kinetic Monte Carlo and microkinetic modeling for catalyst design studies in the presence of lateral interactions

Popular computational catalyst design strategies rely on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends. The main benefit of this approach is related to the inherent dimensionality reduction of the large material space to just a few catalyst descriptors. Conversely, it is well documented that a small set of descriptors is insufficient to capture the intricacies and complexities of a real catalytic system. The inclusion of coverage effects through lateral adsorbate-adsorbate interactions can narrow the gap between simplified descriptor predictions and real systems, but mean-field MKMs cannot properly account for local coverage effects. This shortcoming of the mean-field approximation can be rectified by switching to a lattice-based kinetic Monte Carlo (kMC) method using cluster expansion representation of adsorbate–adsorbate lateral interactions. Using the prototypical CO oxidation reaction as an example, we critically evaluate the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. Here, the mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, the MKM approach does provide a better entry point for computational catalyst design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rate expressions and kinetic parameters for metal ferrites in relation to applications of fossil fuel conversion to hydrogen: Part 1 of 2

Here, the goal of the present work was to provide the necessary reaction emulation information to enable detailed process simulation of a chemical looping H 2 production system from fossil fuels using CaFe 2 O 4 . This specifically pertained to the necessary kinetic data, reaction model development, and model rate parameters required for reaction emulation in both reducing and oxidizing environments. A logical methodology was defined, which included discretization of the reaction network, establishing a core model for reaction emulation that could be adapted based on the system phenomena, and development of a rate parameter regression tool designed around the core model. An extensive array of data sets was acquired by which parametric regressions were performed. The work presented and tabulated a comprehensive set of rate parameters for the reduction and oxidation reactions of CaFe 2 O 4 and descendent phases of Ca 2 Fe 2 O 5 , FeO, Fe 3 O 4 , Fe, and CaO to emulate reaction behavior in a looping-based process environment. This included direct reduction using CH 4 , H 2 , and CO, and direct oxidation reactions with steam, CO 2 and O 2 . Dynamic equilibrium was quantified for reactions that could utilize H 2 O and CO 2 as soft oxidants to re-saturate lattice oxygen in the depleted structure/phases. The kinetics associated with the oxidative mechanisms with the soft oxidants were quantified and compared to those of the reducing counterparts. The analysis provided critical insight to emulate reactions for a process that seeks to use natural gas (NG) or other fossil fuels as a direct reductant for the end goal of H 2 production.

calcium ferrite oxygen carriers↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Flame Chemistry Workshop: a perspective on challenges and strategic actions in combustion experiments and chemical kinetics modeling

Continued progress in the development of predictive models for combustion chemistry—including ignition and flame behavior, species evolution, and combustion system performance—relies on overcoming enduring and emerging challenges in experimental measurements, theoretical formulations, and chemical kinetics mechanism construction. As combustion science continues to coincide with advances in sustainable fuels development, plasma technologies, and automated modeling capabilities, the need for coordinated, community-driven strategies is essential. The Flame Chemistry Workshop (FCWS), held biennially before the International Symposium on Combustion, serves as a dedicated platform to identify, consolidate, and address these challenges in a structured and collaborative manner. This perspective arises from discussions at the 7th FCWS in Milan, Italy (2024), and presents a collective view of the critical barriers currently limiting progress. Across the five technical domains discussed during the 7th FCWS – sustainable fuels combustion, advanced diagnostics for combustion measurements, experiments and modeling in plasma combustion, artificial intelligence and automated methods for theory and mechanisms generation, and chemical kinetic models—a series of persistent and emerging scientific challenges were identified, highlighting the need for deeper integration between three areas: theory, experiments, and modeling. In conclusion, the present article concisely describes present challenges that were identified in each of the technical domains in an effort to streamline and coordinate solutions to accelerate progress in combustion science.

Chemical kinetics↗
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 OSTI2025

Pre-steady-state kinetics of nanocrystal:molybdenum nitrogenase biohybrids reveals hole-scavenging efficiency is critical to N 2 reduction

Molybdenum (Mo) nitrogenase is a two-component enzyme complex that catalyzes the reduction of dinitrogen to ammonia and protons to hydrogen gas. We have shown that electrons for dinitrogen reduction can be delivered photochemically to the catalytic MoFe protein component by cadmium sulfide (CdS) nanocrystals. In this study, we used electron paramagnetic resonance spectroscopy to measure the transient populations of catalytic intermediates. We fit the populations with a pre-steady-state kinetic model, which allowed us to distinguish between productive and non-productive reaction pathways and extract the rate constants for the reaction. Our results demonstrated that the rate of catalytic electron delivery into MoFe protein increased with the concentration of the sacrificial electron donor. This enabled electron delivery to exceed the rate of hydride protonation, a relaxation pathway that competes with N 2 binding. Thus, managing the balance between electron transfer and hole transfer reactions is required to achieve a kinetic regime that favors N 2 reduction.

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

The Role of Evolving Interfacial Substrate Properties on Heterogeneous Cellulose Hydrolysis Kinetics

Interfacial enzyme reactions require formation of an enzyme–substrate complex at the surface of a heterogeneous substrate, but often multiple modes of enzyme binding and types of binding sites complicate analysis of their kinetics. Excess heterogeneous substrate is often used as a justification to model the substrate as unchanging, but using the study of the enzymatic hydrolysis of insoluble cellulose as an example, we argue that reaction rates are dependent on evolving substrate interfacial properties. We hypothesize that the relative abundance of binding sites on cellulose where hydrolysis can occur (productive binding sites) and binding sites where hydrolysis cannot be initiated or is inhibited (nonproductive binding sites) contribute to rate limitations. We show that the initial total number of productive binding sites (the productive binding capacity) determines the magnitude of the initial burst phase of cellulose hydrolysis, while productive binding site depletion explains overall hydrolysis kinetics. Furthermore, we show that irreversibly bound surface enzymes contribute to the depletion of productive binding sites. Our model shows that increasing the ratio of productive to nonproductive binding sites promotes hydrolysis, while maintaining an elevated productive binding capacity throughout conversion is key to preventing hydrolysis slowdown.

09 BIOMASS FUELS↗
DOE OSTI2021

Dispersion relation for gauge-free electromagnetic drift kinetics

Recently, a new approach to gyrokinetics, invariant under electromagnetic gauge transformations, was developed. The gyrocenter equations of motion are now expressed in terms of the perturbed fields instead of the potentials, in a form suitable for numerical simulations and analytic studies. In this paper, we verify that the long-wavelength limit, i.e., the drift-kinetic limit of the new gyrokinetic theory, is in line with existing work, providing a solid foundation for simulations. Here, we compute the dispersion relation of the new drift-kinetic theory in slab geometry and find agreement with a long-wavelength limit of the full Vlasov–Maxwell model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

In situ thermal conductivity measurement revealing kinetics of thermochemical reactions

Utilizing thermochemical reactions for thermal energy storage and solar fuel production has been an emerging research topic. Thermal transport properties of the materials are an important parameter that can determine the kinetics and efficiency of thermochemical reactions. With the increasing number of new thermochemical materials (TCMs); however, there is a lack of reliable techniques to monitor the thermal transport property of the materials and their changes as a function of reactions in real time. In this work, we report the in situ monitoring of thermochemical reactions using modulated photothermal radiometry (MPR). The thermal conductivities of two TCMs, namely, calcium hydroxide (Ca(OH) 2 ) and Ba 0.15 Sr 0.85 FeO 3–δ (BSF1585), were measured as a function of temperature and time using the MPR technique. The measured thermal conductivities were correlated to the reaction. The work has two significant contributions to the research communities. First, it provides a non-invasive diagnostic tool for monitoring the thermal transport properties of TCMs that can potentially be a high-throughput measurement technique conducive to optimizing TCMs, reactors, and related thermal systems. Second, for TCMs that show observable changes in thermal transport properties, a correlation between the measured thermal conductivity and the conversion fraction of the reaction can be established for monitoring the reaction kinetics based on thermal characterization.

14 SOLAR ENERGY↗
DOE OSTI2023

Gas-phase microactuation using kinetically controlled surface states of ultrathin catalytic sheets

Biological systems convert chemical energy into mechanical work by using protein catalysts that assume kinetically controlled conformational states. Synthetic chemomechanical systems using chemical catalysis have been reported, but they are slow, require high temperatures to operate, or indirectly perform work by harnessing reaction products in liquids (e.g., heat or protons). Here, we introduce a bioinspired chemical strategy for gas-phase chemomechanical transduction that sequences the elementary steps of catalytic reactions on ultrathin (<10 nm) platinum sheets to generate surface stresses that directly drive microactuation (bending radii of 700 nm) at ambient conditions (T = 20 °C; P total = 1 atm). When fueled by hydrogen gas and either oxygen or ozone gas, we show how kinetically controlled surface states of the catalyst can be exploited to achieve fast actuation (600 ms/cycle) at 20 °C. Finally, we also show that the approach can integrate photochemically controlled reactions and can be used to drive the reconfiguration of microhinges and complex origami- and kirigami-based microstructures.

42 ENGINEERING↗
DOE OSTI2024

Relativistic Guiding-Center Motion: Action Principle, Kinetic Theory and Hydrodynamics

We treat the guiding-center dynamics in a varying external Maxwell field using a relativistically covariant action principle which reproduces the known Vandervoort expression for the drift velocity and extends it to curved spacetime. We derive the corresponding kinetic theory and ideal hydrodynamic theory. In contrast to conventional five-equation hydrodynamics, the guiding-center hydrodynamics needs only three equations due to a constraint on the motion across magnetic field. Furthermore, we argue that such a hydrodynamics is applicable to strongly coupled plasmas where kinetic theory fails.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A Theoretical Study of NH 2 Radical Reactions with Propane and Its Kinetic Implications in NH 3 -Propane Blends’ Oxidation

The reaction of NH 2 radicals with C 3 H 8 is crucial for understanding the combustion behavior of NH 3 /C 3 H 8 blends. In this study, we investigated the temperature dependence of the rate coefficients for the hydrogen abstraction reactions of C 3 H 8 by NH 2 radicals using high-level theoretical approaches. The potential energy surface was constructed at the CCSD(T)/cc-pV(T, Q)//M06-2X/aug-cc-pVTZ level of theory, and the rate coefficients were computed using conventional transition state theory, incorporating the corrections for quantum tunneling and hindered internal rotors (HIR). The computed rate coefficients showed a strong curvature in the Arrhenius behavior, capturing the experimental literature data well at low temperatures. However, at T > 1500 K, the theory severely overpredicted the experimental data. The available theoretical studies did not align with the experiment at high temperatures, and the possible reasons for this discrepancy are discussed. At 300 K, the reaction of NH 2 with C 3 H 8 predominantly occurs at the secondary C-H site, which accounts for approximately 95% of the total reaction flux. However, the hydrogen abstraction reaction at the primary C-H site becomes the dominant reaction above 1700 K. A composite kinetic model was built, which incorporated the computed rate coefficients for NH 2 + C 3 H 8 reactions. The importance of NH 2 + C 3 H 8 reactions in predicting the combustion behavior of NH 3 /C 3 H 8 blends was demonstrated by kinetic modeling.

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

Operando Neutron Imaging of Lithium Flux and Gradient Cathode Design for Enhanced Kinetics in High‐Loading All‐Solid‐State Li─S Batteries

All-solid-state Li–sulfur batteries (ASSLSBs) are considered promising candidates for next-generation energy storage owing to their inherent safety, high energy density, and abundant sulfur resources. However, slow redox kinetics greatly limit sulfur utilization during solid-solid sulfur reactions, leading to significant challenges to achieve efficient performance in high-mass-loading ASSLSBs. Here, operando neutron image is employed to directly visualize, for the first time, that sluggish Li + transport kinetics and the uneven distribution of Li + during cathodic reactions are critical factors limiting sulfur conversion. To address this issue, gradient cathode architectures comprising three and five layers are designed, in which catholyte concentrations are strategically varied to optimize Li-ion flux and enhance ionic conductivity of the whole composite cathode electrode. Operando neutron imaging distinctly visualizes and confirms that three-layer gradient approach significantly enhances Li-ion mobility, resulting in more uniform redox reactions and greatly improved sulfur utilization compared to traditional non-gradient structures. Consequently, the three-layer gradient cathode achieves superior rate performance and reduced electrode polarization at high sulfur mass loadings of 4.5 and 6.0 mg cm −2 . Furthermore, the applicability and scalability of this design are demonstrated in a five-layer gradient cathode architecture, achieving an impressive discharge specific capacity increase from 656 mAh g −1 (three-layer gradient) to 1232 mAh g −1 at 1/20 C for ultra-high sulfur loading of 7.5 mg cm −2 . In conclusion, this innovative gradient cathode design offers substantial advancements in understanding and overcoming Li-ion transport limitations, paving the way toward practical, high-energy-density ASSLSBs.

25 ENERGY STORAGE↗
DOE OSTI2020

Toward Higher Voltage Solid-State Batteries by Metastability and Kinetic Stability Design

The energy density of battery systems is limited largely by the electrochemical window of the electrolyte. Here, the combined thermodynamic and kinetic effects of mechanically induced metastability are shown to greatly widen the operational voltage window of solid-state batteries based on ceramic-sulfide electrolytes. Solid electrolyte voltage stability up to 10 V is achieved with minimal degradation, far beyond the capability of organic liquid electrolytes. Furthermore, combined experiment, ab initio computation, and theoretical modeling identify the nature of mechanically constrained Li 10 GeP 2 S 12 decomposition both within the bulk and at interfaces with cathode materials at very high voltages. Previously unclear kinetic processes are identified that, when properly implemented, can potentially allow solid-state full cells with remarkably high operational voltages.

25 ENERGY STORAGE↗
DOE OSTI2025

Contrasting Chemical Kinetic Parameters for Near In‐Service‐Temperature Aging and High‐Temperature Thermal Decomposition of Triaminotrinitrobenzene Formulations

Triaminotrinitrobenzene-based high explosives such as LX-17 offer high energy density and exceptional safety, yet their long-term aging behavior at low temperatures remains poorly understood. In this study, several legacy and new production lots of LX-17 were subjected to accelerated aging experiments below 100°C, during which the formation rate of the initial degradation product, monofurazan (F1), was monitored. Kinetic analysis was performed using a sample-age-aware computational approach, yielding activation energy estimates of 82–91 kJ mol −1 for the low-temperature initiation step—markedly lower than the ∼200 kJ mol −1 associated with high-temperature thermal decomposition. Extrapolation from established cookoff models supports the conclusion that the dominant degradation mechanism at low temperatures differs from that at high temperatures. In conclusion, our results provide a unified kinetic framework that bridges the gap between in-service conditions and high-temperature damage.

Chemistry - Chemical explosives↗
DOE OSTI2024

Enhancing the Reaction Kinetics and Stability of Co‐Free Li‐Rich Cathode Materials via a Multifunctional Strategy

Co-free Li-rich layered oxides (CFLLOs) with anionic redox activity are among the most promising cathode materials for high-energy-density and low-cost lithium-ion batteries (LIBs). However, irreversible oxygen release often causes severe structural deterioration, electrolyte decomposition, and the formation of unstable cathode-electrolyte interface (CEI) film with high impedance. Additionally, the elimination of cobalt elements further deteriorates the reaction kinetics, leading to reduced capacity and poor rate performance. Here, in this study, a multifunctional strategy is proposed, incorporating Li 2 MnO 3 phase content regulation, micro-nano structure design, and heteroatom substitution. The increased content of Li 2 MnO 3 phase enhances the capacity through oxygen redox. The smaller nanoscale primary particles induce greater tensile strain and introduce more grain boundaries, thereby improving the reaction kinetics and reactivity, while the larger micron-sized secondary particles help to reduce interfacial side reactions. Furthermore, Na⁺ doping modulates the local coordination environment of oxygen, stabilizing both the anion framework and the crystal structure. As a result, the designed cathode exhibits enhanced rate performance, delivering a capacity of 158 mAh g⁻¹ at 5.0 C and improved cyclic stability, with a high capacity retention of 99% after 400 cycles at 1.0 C. This multifunctional strategy holds great promise for advancing the practical application of CFLLOs in next-generation LIBs.

Co-free Li-rich layered oxide↗
DOE OSTI2019

Kinetic Study on Alloying Element Transfer During an Electroslag Remelting Process

Experimental and theoretical studies have been carried out to investigate the effects of slag on the alloying elements in ingots during the electroslag remelting (ESR) process with a focus on developing a mass-transfer model to understand the mechanism of slag-metal reaction. Stainless steel 1Cr21Ni5Ti was used as the electrode and remelted with two different kinds of slags using a 50-kg ESR furnace. The contents of sulfur, aluminum, titanium and silicon along the axial direction of the produced ingots were analyzed. On the basis of the penetration and film theories, the theoretical model developed in this work elucidates the kinetics of the slag-metal reaction revealing the mechanism of alloying element transfer during the ESR process. The calculation results obtained from the model agree well with the experimental results. The model indicates that the resultant [O] coming from the desulfurization reaction of (O{sup 2−}) + [S] = (S{sup 2−}) + [O] causes the oxidation of alloying elements in steel by [M] + [O] = (MO). The distribution ratio of sulfur L{sub S} decreases with the increase of slag temperature in the first slag-temperature-rising period, and the concentration of sulfur in the ingot at the beginning of the ESR process is lower than in the rest of the process because of the combination of the large distribution ratio of sulfur L{sub S} and excellent kinetic conditions. The concentration of aluminum along the height of the ingot has an increasing trend in the first slag-temperature-rising period, while it has a decreasing trend in the rest of the process. Two methods can solve this problem: one is starting up the ESR furnace by high temperature molten slag technology and the other is continually adding extra titania into the molten slag in the first slag-temperature-rising period.

36 MATERIALS SCIENCE↗