Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High-Temperature Fe-Based Fischer–Tropsch Synthesis: Experimentally Validated Kinetic Models Implemented at Pellet and Reactor Scales

Pellet- and reactor-scale models for Fischer–Tropsch synthesis (FTS) with a Fe–K/silica catalyst were developed to investigate the sensitivity of the hydrocarbon products and carbon dioxide selectivity to process conditions and feed composition at high temperature (350–400 °C), moderate pressure (1–10 bar), and a range of H 2 /CO ratios (3–1). The major objective of this paper is to develop, validate, and evaluate a high-temperature FTS model that is then used to assess the feasibility of process integration with syngas production. Since there is limited kinetic data available, in literature at these conditions, bench-scale reactor tests were conducted to obtain operational data for parameter fitting of kinetic expressions used in the model. This resulting kinetic model demonstrated agreement with the experimental data with an R 2 of 0.97 to the testing data set and, thus, was feasible to apply at pellet and reactor scales. Here, multiple pellet sizes were modeled to detail the role of transport limitations as the sphere’s diameter approached and exceeded 1 mm. Application of the reactor model indicated that hydrocarbon selectivity depended strongly on temperature, whereas the ratio of olefin to paraffin products decreased with increasing temperature, pressure, and H 2 /CO ratio. Product selectivity was not sensitive to the conversion of carbon monoxide. Furthermore, the roles of the pressure and H 2 /CO ratio were closely coupled. At a H 2 /CO ratio of 3, only slight variations in selectivity occurred over a pressure range of 1–20 bar, whereas at a ratio of 1, selectivity could vary by as much as 30% over the same pressure range. At pressures below 5 bar and temperatures above 350 °C, minimal selectivity to heavy hydrocarbons (C 12+ ) is obtained, and selectivity to midrange products (C 5–11 ) rapidly declined as pressure dropped below 5 bar, which indicated that an operational pressure of at least 5 bar is needed to achieve reasonable yields in this temperature range. These results, while tentative, provide guidelines for further experimentation and evaluation of integrated FTS processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2024

Chemical insights into ethyl acetate flames from experiment and kinetic modeling: Laminar burning velocity, speciation and NO$_x$ emission

Oxygenated fuels, such as alcohols, ethers, and esters, are promising alternatives to conventional fuels. These fuels can help reduce detrimental emissions like carbon monoxide and unburned hydrocarbons and enhance octane ratings. Among these oxygenates, ethyl acetate (EA), a small alkyl ester sourced from biomass, emerges as a clean, promising energy carrier. It serves as a surrogate fuel to facilitate investigations into the combustion behaviours of biodiesel. Despite its importance, the literature knowledge of EA combustion characteristics is limited. Therefore, this study aims to broaden the knowledge of the combustion behaviour of this type of oxygenated fuel compound. In this study, we measured the laminar burning velocities of EA by employing a heat flux burner and a closed combustion vessel over the equivalence ratios of 0.7 – 1.7, pressures of 1 – 10 bar and temperatures ranging from 353 – 423 K. Further, we also measured the NOx emissions in exhaust gas of the premixed flames fueled by EA/air for the first time over the equivalence ratio of 0.8 – 1.2. Additionally, we employed a non-premixed counterflow flame setup for extensive characterisation of species and their concentration under diverse conditions encompassing various strain rates and oxygen concentrations. Finally, we utilized these newly measured data to construct and validate a detailed kinetic model developed as part of this work. The newly developed model will help characterize the combustion properties of EA.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

An experimental, theoretical and kinetic modeling study of the N 2 O-H 2 system: Implications for N 2 O + H

The reaction of N 2 O with His the key step in consumption of nitrous oxide in thermal processes. The major product channel is N 2 + OH, while NH + NO constitute minor products. In addition, a pathway involving HNNO, initiated by N 2 O + H (+M)$\rightleftarrows$HNNO (+M) (R3, R4), has been inferred from experiment and theory by Burke and coworkers. At longer reaction times, the reaction may reach partial equilibration, and in addition to k 3 and k 4 the importance of this channel depends on the thermodynamic properties of HNNO and its consumption reactions, mainly HNNO + H. In the present work, we re-examined the thermochemistry of HNNO and calculated rate constants and branching fractions for the HNNO + H reaction. Experiments on the N 2 O-H 2 system were conducted in a high-pressure flow reactor at 100 atm as a function of temperature (600-925 K) and stoichiometry and explained in terms of an updated chemical kinetic model. The results support the importance of the HNNO pathway, which results in inhibition of N 2 O consumption and formation of NH 3 . In addition, selected literature results on the N 2 O-H 2 system are re-examined and the implications for the other product channels of N 2 O + H, in particular NH + NO, are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 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 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.↗
DOE OSTI2025

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Kinetics Modeling and Reactor Design Study of Glucose-to-Terpenes Cell-Free Conversion

Cell-free systems offer many advantages over traditional biological conversion by eliminating biological growth constraints. It also offers easy manipulation and finetuning of the reaction conditions for each individual enzyme. The conversion of cellulosic glucose to Limonene, a terpene, is a promising pathway for producing fuels and chemicals. Recent advances in developing cell-free systems focuses on bench scale optimization of terpene yield and to demonstrate its feasibility towards commercialization [1,2]. There is significant knowledge gap regarding reaction kinetics of these cell-free systems to further study how it will perform at larger scale. We present here, our studies on reaction kinetics and reactor design implications of cell-free glucose to Limonene conversion to facilitate the further development and commercialization of this process. We developed a novel kinetic model based on the metabolic-network structure of the cell-free system with multi-substrate reversible Michaelis-Menten rate law. To estimate kinetic parameters for this system of rate equations, we employed Bayesian optimization to perform global search with the assistance of gaussian processes to balance exploration and exploitation. The model parameters estimated showed good results compared with experimental data. The estimated parameters were used to perform sensitivity analysis. We found that Hexokinase is one of the most critical enzymes that affect the conversion of the glucose. We also observed that abundance of co-factors is also critical to the conversion of glucose to limonene. We investigated packed bed reactors with enzymes immobilized on the surface of particles to convert glucose stream into Limonene for larger scale production. The reactor design such as particle size, enzyme loading, and flow rate are found to be critical for improving yields. [1] Dudley, Q.M., Nash, C.J. and Jewett, M.C., 2019. Synthetic Biology, 4(1), p.ysz003. [2] Korman, T.P., Opgenorth, P.H. and Bowie, J.U., 2017. Nature communications, 8(1), p.15526.

09 BIOMASS FUELS↗
DOE OSTI2025

VERITAS : A density-functional theory-based multiband kinetic model for understanding x-ray spectroscopy of dense plasmas

X-ray spectroscopy has long been a powerful diagnostic tool for hot, dilute plasmas, providing insights into plasma conditions by measuring line shifts and broadenings of atomic transitions. The technique critically depends on the accuracy of atomic physics models used to interpret spectroscopic measurements for inferring plasma properties such as free-electron density and temperature. Over the past decades, the atomic and plasma physics communities have developed robust atomic physics models to account for various processes in hot, dilute classical plasmas. While these models have been successful in that regime, their applicability becomes uncertain when interpreting x-ray spectroscopy experiments of above-solid-density plasmas. Given that finite-temperature density-functional theory (DFT) offers a more accurate description of dense plasma environments, we present the development of a DFT-based multi-band kinetic model, VERITAS, designed to improve the interpretation of x-ray spectroscopic measurements in high-density plasmas produced by laser-driven spherical implosions. This work details the VERITAS model and its application to both time-integrated and time-resolved x-ray spectra from implosion experiments on OMEGA. The advantages and limitations of the VERITAS model will also be discussed, along with potential directions for advancing x-ray spectroscopy of dense and superdense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling kinetic effects of charged vacancies on electromechanical responses of ferroelectrics: Rayleighian approach

Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025

Kumar, Rajeev (ORCID:0000000194943488)↗