SEARCH · Engineering Papers

Results for “kinetic model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Kinetic Modeling of Secondary Organic Aerosol in a Weather-Chemistry Model: Parameterizations, Processes, and Predictions for GOAmazon

Secondary organic aerosol (SOA) forms and evolves in the atmosphere through many pathways and processes, over diverse spatial and time scales. Hence, there is a need to represent these widely-varying kinetic processes in large-scale atmospheric models to allow for accurate predictions of the abundance, properties, and impacts of SOA. In this work, we integrated a kinetic, process-level model (simpleSOM-MOSAIC) into a weather-chemistry model (WRF-Chem) to simulate the oxidation chemistry and microphysics of atmospheric SOA. simpleSOM-MOSAIC simulates multigenerational gas-phase chemistry, autoxidation reactions, heterogeneous oxidation, oligomerization, and phase-state-influenced gas/particle partitioning of SOA. As a case study, the integrated WRF-Chem-simpleSOM-MOSAIC (WC-SSM) model was used to simulate the photochemical evolution downwind of a large city (Manaus, Brazil) in the Amazon and, in turn, study the anthropogenic and biogenic interactions in an otherwise pristine environment. Consistent with previous work, we found that OA was enhanced by up to a factor of four in the urban plume due to elevated hydroxyl radical (OH) concentrations, relative to the background, and that this OA was dominated by SOA from biogenic precursors (80%). Further, in addition to accurately simulating the OA enhancement in the urban plume, the model reproduced the magnitude of the OA oxygen-to-carbon (O:C) ratio and broadly tracked the evolution of the aerosol size distribution. Our work highlights the importance of including an integrated, kinetic representation of SOA processes in an atmospheric model

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Kinetic Model of Photochemical Nitrogen Reduction [SWR-25-97]

The code in this repository was used to model the kinetics of photochemical nitrogen reduction. Each directory represents a manuscript. The directory titled Dahl_CRPS_2025 contains the code used to fit the kinetic models presented in the manuscript titled, "Pre-steady-state kinetics of nanocrystal:molybdenum nitrogenase biohybrids reveals hole-scavenging efficiency is critical to achieving N2 reduction" by Peter J. Dahl, Lauren M. Pellows, Zhi-Yong Yang, Lance C. Seefeldt, John W. Peters, Gordana Dukovic, David W. Mulder, and Paul W. King. 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 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 sacrificial electron donor. This enabled electron delivery to exceed the rate of hydride protonation, a relaxation pathway which competes with N2 binding. Thus, managing the balance between electron transfer and hole transfer reactions is required to achieve a kinetic regime that favors N2 reduction.

Dahl, Peter [National Renewable Energy Laboratory ↗

DOE OSTI2025

Measurements and kinetic modeling of O 2 vibrational kinetics in O 2 –Ar mixtures partially dissociated by a Ns pulse discharge

Vibrational kinetics of O 2 is studied during the O atom recombination in an O 2 –Ar mixture, partially dissociated by a burst of ns discharge pulses in a heated plasma flow reactor. The time-resolved temperature in the discharge afterglow is determined by Rayleigh scattering. Time-resolved O atom number density is measured by ps Two-Photon absorption Laser Induced Fluorescence, calibrated in xenon. Time-resolved vibrational level populations of molecular oxygen, O 2 (v= 8–20), are measured by ps Laser Induced Fluorescence (LIF), with the absolute calibration by NO LIF. Time-resolved ozone number density is monitored by broadband UV absorption. The results are compared with the predictions of a state-specific kinetic model. The experimental data indicate a rapid initial decay of O 2 (v) populations generated by electron impact in the discharge, due to the vibration-translation (V–T) relaxation by O atoms. This is followed by a slower population reduction, on the time scale much longer compared to that for V–T relaxation or vibration-vibration (V–V) exchange. Both O atoms and the O 2 (v) populations decay on the same time scale, indicating that chemical reactions initiated by the O atom recombination result in the generation of vibrationally excited O 2 molecules. These trends are reproduced by the kinetic model, which shows that the reaction of O atoms with ozone is the dominant pathway of O 2 (v) generation at the present conditions. The predicted relative O 2 (v) populations are close to the experimental results, but absolute number densities differ from the experimental data. This is likely due to uncertainties in the absolute calibration of LIF measurements and in the spectroscopic model used in the data reduction. The present work demonstrates the capability for the absolute, time-resolved measurements of vibrationally excited O 2 in recombining gas flows, to quantify the energy partition in the recombination reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Data for KETCHUP: Parameterizing of Large-Scale Kinetic Models Using Multiple Datasets with Different Reference States

Repository for Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP can use K-FIT [2] input files. Example K-FIT input files are located in the K-FIT repository at https://github.com/maranasgroup/K-FIT.

Metabolomics↗

DOE OSTI2025

Toward a Unified Kinetic Model of Nitrogenase Catalysis

The microbial enzyme nitrogenase catalyzes the MgATP-dependent reduction of N 2 to 2NH 3 , a transformation central to the global nitrogen cycle. While the canonical Thorneley−Lowe (TL) kinetic model has long served as a mechanistic framework, it does not incorporate several recent insights. Here, we present an updated kinetic model for Monitrogenase that incorporates these new findings. A significant insight is that electron transfer (ET) from the reduced Fe protein to the FeMo-cofactor is gated by MgATP-dependent conformational transitions and can be described as a probabilistic event that is dependent on the ligand bound to the active-site metallocofactor. The updated kinetic model quantitatively reproduces steady-state product formation rates across a broad range of experimental conditions, yielding revised estimates for key rate constants. It is demonstrated that under N 2 turnover, the probability of productive ET to the active site decreases by ∼60%, resulting in a significant fraction of Fe protein cycles that are unproductive for electron delivery. This mechanistic feature explains the observed rate limitation in N 2 reduction and implies a revised minimum energetic cost of approximately 25 MgATP per N 2 reduced. Integrating these new features into the revised kinetic model provides a more complete and usable foundation for understanding nitrogenase catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗

DOE OSTI2026

Flow Reactor Study and Kinetic Model Development of HEFA-SPK and its Surrogate

In this work, we formulate a two-component surrogate for HEFA-SPK, incorporating aromatic or cycloalkane components, and develop reduced kinetic models for the surrogates to be used in high-performance computing simulations. The HEFA-SPK surrogate was selected and optimized based on the fuel's physical and combustion properties, including ignition delay times and flame speeds. The resulting surrogate consists of 40% n-undecane and 60% 2-methylnonane. The surrogate was confirmed by flow reactor experiments for both the HEFA-SPK fuel and the suggested two-component surrogates, where excellent agreement was observed. To meet aromatic requirements, 1,2,4-trimethylbenzene (8%) was selected, and we determined that incorporating 30% propylcyclohexane into the HEFA-SPK will achieve a volume swell equivalent to 8% aromatics. The properties of the surrogates were measured, and a new reduced kinetic model was developed based on the semi-decoupling methodology, using a reduced CH4 chemistry from NUIG 1.0 as base chemistry. Kinetic models will be employed in combustor simulations to enable a comprehensive understanding of the effect of SAF fuel properties on aviation combustor performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗

DOE OSTI2025

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗

DOE OSTI2025

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

kinetics↗

DOE OSTI2025

Challenges in integrating dissolved organic matter chemodiversity into kinetic models of soil respiration

The chemodiversity of dissolved organic matter (DOM) in soil has been proposed to influence the microbial metabolism and fate of belowground organic carbon (C). However, integrating DOM chemistry into soil C cycle models to improve predictions of C stocks and fluxes—beyond simply considering DOM pool size—remains a challenge. While recent research suggests that incorporating DOM chemodiversity into models can improve predictions of microbial respiration, there is still a lack of mechanistic understanding describing how DOM chemodiversity affects microbial metabolism and soil respiration. Here, we evaluated whether DOM chemodiversity was a determinant of soil respiration using paired measurements of high-resolution DOM chemistry, obtained from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), and potential soil respiration rates from across the United States (U.S.), all data provided by the Molecular Observation Network. Our objectives were to (1) assess statistical relationships between DOM chemodiversity and microbial respiration, and (2) evaluate the ability of kinetic models to leverage DOM chemistry to explain empirical relationships found in statistical models. Statistical regressions revealed that DOM chemodiversity (alpha diversity) was nonlinearly related to potential soil respiration rates, both independently and through its interactions with DOM and total C concentrations. In soils with relatively high DOM but low total C concentrations, potential soil respiration rates were negatively correlated with DOM alpha diversity, whereas in soils with relatively low DOM and high total C concentrations showed the opposite trend. However, when metabolic transition theory kinetic models were modified to include chemodiversity, their performance was comparable to traditional Monod kinetics approaches, which simulate respiration rates as a function of DOM concentration. The inability to account for nonlinearities in DOM chemodiversity–respiration relationships highlight an opportunity to advance substrate uptake kinetics by establishing causal links between DOM chemodiversity, microbial metabolism trade-offs, and potential interactions under varied environmental conditions.

Bioenergetic model↗

DOE OSTI2026

Kinetic model development for single step ethanol to butene rich olefin process over Cu-Y/Beta catalysts

Here, this study presents the first intrinsic kinetic model for the single-step conversion of ethanol to butene-rich olefins over bifunctional Cu-Y/Beta catalysts, addressing a critical gap in the design and scale-up of Sustainable Aviation Fuel (SAF) processes. The reaction network comprises ten global steps involving dehydrogenation, aldol condensation, hydrogenation, and dehydration reactions, distributed across Cu and rare-earth (Y) active sites. The model incorporates dual-site functionality (Cu and Y site) and explicitly accounts for key intermediates such as crotonaldehyde and butanal. Reaction rates are formulated using Langmuir–Hinshelwood–Hougen–Watson (LHHW) kinetics. Kinetic parameters are extracted by fitting the model to lab-scale packed-bed reactor data across a wide range of temperatures and space velocities, demonstrating strong agreement in ethanol conversion and product selectivity. The reaction kinetics developed in this work provide a foundational basis for constructing reactor models that enable process optimization and scale-up of ethanol-to-jet fuel technologies.

Cu-Y/Beta catalyst↗

DOE OSTI2025

Density functional theory-based surrogate kinetic models for heterogeneous reactions of hydrocarbon intermediates on silicon carbide

The increasing demand for high-performance materials in advanced technologies highlights the importance of achieving a fundamental understanding and potential control of silicon carbide (SiC) deposition processes. However, existing models often lack sufficient theoretical detail, relying heavily on empirical data and offering limited predictive capability. In particular, the complex surface chemistry governing SiC growth remains poorly understood. This study addresses these challenges by employing density functional theory (DFT) to investigate key heterogeneous reactions involving hydrocarbon intermediates on SiC surfaces, including dehydrogenation, hydrogenation, and carbon deposition. Transition state searches were conducted to identify reaction pathways and energy barriers. While first-principles calculations offer high accuracy, they are computationally intensive. To extend the utility of these first-principles results, vibrational analyses were performed using phonon-based statistical thermochemistry to compute temperature-dependent reaction rates which were used to develop Arrhenius-type surrogate kinetic models. Furthermore, the resulting framework provides a more rigorous, physically grounded basis for integrating atomistic insights into continuum-scale modeling, ultimately enabling improved prediction and optimization of SiC film growth in high-performance material systems.

Density Functional Theory↗

DOE OSTI2026

Fluid-kinetic modeling of a high power density radio frequency inductively coupled positive hydrogen ion source

High power density radio-frequency (RF) inductively coupled positive ion sources are attractive candidates for next-generation neutral beam injection (NBI) systems, where higher injected power and longer pulse lengths are desired without sacrificing source reliability. Operating at absorbed power densities of order $\gt 1~\mathrm{W\,cm}^{-3}$ places these sources in a regime with stronger gas heating, higher dissociation, and non-Maxwellian electron energy distributions. The Large Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma source designed to explore this high power density regime and to provide guidance for a positive ion source upgrade for the DIII-D NBI system. LUPIN is designed to operate at up to 20 kW of RF power at 2 MHz, coupling energy through a cylindrical quartz vessel to achieve target ion current densities of $2100\,\mathrm{A\,m}^{-2}$ . This paper presents fluid-kinetic modeling of LUPIN using the hybrid plasma equipment model where electrons are treated kinetically, and the simulations reveal that electron energy distribution function transitions from nearly Maxwellian in the core to bi-Maxwellian towards the edge. Parametric simulations investigate the effects of RF power, gas pressure, and frequency on plasma density, ion flux, and uniformity. Parametric sweeps reveal that increasing power shifts the primary ionization channel from molecular to atomic with diminishing flux gains due to skin-depth contraction and gas rarefaction. Higher frequency localizes heating and increases $\mathrm{H}_2^+$ and $\mathrm{H}_3^+$ delivery to the grid, while elevated pressure boosts ionization yet hinders ion transport due to increase in collisionality.

inductively coupled plasma↗

DOE OSTI2025

Fermi-Based Kinetic Model for the Sabatier Reaction: Sabatier Principle and Beyond It

The CO 2 methanation reaction (aka the Sabatier reaction) is one of the important carbon utilization technologies to reduce atmospheric greenhouse gases. In this study, kinetic modeling of the Sabatier reaction proceeding via the CO route has been presented by considering weak and strong chemisorption concepts rooted in electronic theory. Based on the single-route mechanism of the Sabatier reaction, the dependence of the reaction rate on the position of the catalyst Fermi level with respect to the antibonding states of the reagent (CO 2 ) and the product (CH 4 ), is demonstrated with the classic Sabatier volcano-shape. Furthermore, the model was applied to explain experimental results of reaction rates for Ni- and Ru-based catalytic systems and to understand origin of enhancement of turnover frequencies established for dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials

Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.

Kerisit, Sebastien N.↗