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

An experimental and kinetic modeling study of the ignition of 2-methyl decane

Ignition delay times (IDTs) of 2-methyl decane (C 11 H 24 -2) are measured in a high-pressure shock tube and in a rapid compression machine at equivalence ratios in the range 0.5–2.0 at 90 % dilution, at temperatures in the range 600–1430 K and at pressures of 15 and 30 bar. To clarify the effect of the branched methyl group on fuel oxidation, IDTs of n-undecane (nC 11 H 24 ) are also measured at similar conditions to those measured for C 11 H 24 -2. A new chemical kinetic mechanism, using C3MechV4.0.1 as the core chemistry, is developed and validated against the new experimental data. The thermodynamic properties of the fuel (RH), alkyl (Ṙ), alkyl peroxy (RȮ 2 ), hydroperoxy-alkyl (Q̇OOH), and peroxy hydroperoxy alkyl (Ȯ 2 QOOH) radicals are updated in both the C 11 H 24 -2 and nC 11 H24 models using THERM25. A reaction path flux analysis for C 11 H 24 -2 at different temperatures was conducted. Compared to nC 11 H 24 , C 11 H 24 -2 shows slower reactivity. At low and intermediate temperatures, the chain propagation pathway Q̇OOH ↔ C 11 cyclic ether + ȮH is favored for C 11 H 24 -2, while the chain branching pathway Q̇OOH ↔ Ȯ 2 QOOH ↔ C 11 carbonyl hydroperoxide (KHP) + ȮH is suppressed, leading to lower reactivity compared to nC 11 H 24 . At high temperatures, the presence of the branched methyl group inhibits the direct decomposition of the fuel, resulting in reduced C 2 H 4 formation, which in turn suppresses the reactivity of the fuel.

2-methyl decane↗
DOE OSTI2026

A Parallel Kinetic Model for Surface and Bulk Charge Storage in ε -MnO 2 Pseudocapacitors

Pseudocapacitive materials such as manganese dioxide (MnO 2 ) are attractive for energy storage applications due to their ability to combine the fast kinetics of capacitors with the higher energy density of battery-type systems. However, the electrochemical behavior of MnO 2 remains difficult to interpret mechanistically, in part because existing models often fail to distinguish between surface-based redox processes and bulk intercalation mechanisms. In this work, we develop a physics-based model that represents MnO 2 pseudocapacitance as a linear combination of two independent, parallel electrochemical processes: (i) the surface or near-surface redox storage and (ii) lithium ion intercalation into the bulk material. These two processes are treated with distinct kinetic and thermodynamic parameters and are assumed to proceed independently. The total measured current is assumed to be the sum of these two partial currents. We validate the model using rate-dependent cyclic voltammetry experiments, demonstrating that it captures key trends and provides physically interpretable parameters reflecting the relative contributions of surface and bulk processes. By enabling a clear separation between these mechanisms, the model offers a useful framework for analyzing pseudocapacitive materials and can guide the rational design of high-performance energy storage electrodes.

Energy - Storage↗
DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗
DOE OSTI2026

Operando neutron radiography validates a parameter-free transport–kinetics model for thick solid-state battery cathodes

Tortuosity-weighted interfacial flux for lithium (TWIF-Li) predicts through-thickness Li gradients in thick composite all-solid-state cathodes without fitted parameters. Image-derived microstructures, GITT-derived concentration-dependent solid diffusion, and tortuosity-weighted interfacial kinetics reproduce operando neutron radiography across practical rates, delivering transferable design rules to suppress transport-limited reaction fronts.

Adam, Andre [ORNL] (ORCID:0000000245023033)↗
DOE OSTI2025

X-ray crystallographic and hydrogen deuterium exchange studies confirm alternate kinetic models for homolog insulin monomers

Despite the crucial role of various insulin analogs in achieving satisfactory glycemic control, a comprehensive understanding of their in-solution dynamic mechanisms still holds the potential to further optimize rapid insulin analogs, thus significantly improving the well-being of individuals with Type 1 Diabetes. Here, we employed hydrogen-deuterium exchange mass spectrometry to decipher the molecular dynamics of newly modified and functional insulin analog. A comparative analysis of H/D dynamics demonstrated that the modified insulin exchanges deuterium atoms faster and more extensively than the intact insulin aspart. Additionally, we present new insights derived from our 2.5 Å resolution X-ray crystal structure of modified hexamer insulin analog at ambient temperature. Furthermore, we obtained a distinctive side-chain conformation of the Asn3 residue on the B chain (AsnB3) by operating a comparative analysis with a previously available cryogenic rapid-acting insulin structure (PDB_ID: 4GBN). The experimental conclusions have demonstrated compatibility with modified insulin’s distinct cellular activity, comparably to aspart. Additionally, the hybrid structural approach combined with computational analysis employed in this study provides novel insight into the structural dynamics of newly modified and functional insulin vs insulin aspart monomeric entities. It allows further molecular understanding of intermolecular interrelations driving dissociation kinetics and, therefore, a fast action mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Microwave-Assisted Plastic Upcycling: Dynamic Data Reconciliation, Parameter Estimation, and Kinetic Modeling

Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.

Damahe, Harish [West Virginia Univ., Morgantown, W↗
DOE OSTI2025

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Observation-constrained kinetic modeling of isoprene SOA formation in the atmosphere

Isoprene has the largest global non-methane hydrocarbon emission, and the oxidation of isoprene plays a crucial role in the formation of secondary organic aerosol (SOA). Two primary processes are known to contribute to SOA formation from isoprene oxidation: (1) the reactive uptake of isoprene-derived epoxides on acidic or aqueous particle surfaces and (2) the absorptive gas–particle partitioning of low-volatility oxidation products. In this study, we developed a new multiphase condensed isoprene oxidation mechanism that includes these processes with key molecular intermediates and products. The new mechanism was applied to simulate isoprene gas-phase oxidation products and SOA formation from previously published chamber experiments under a variety of conditions and atmospheric observations during the Southern Oxidant and Aerosol Studies (SOAS) field campaign. Our results show that SOA formation from most of the chamber experiments is reasonably reproduced using our mechanism, except when the concentration ratios of initial nitric oxide to isoprene exceed ~ 2, the formed SOA is significantly underpredicted. The SOAS simulations also reasonably agree with the measurements regarding the diurnal pattern and concentrations of different product categories, while the total isoprene SOA remains underestimated. The molecular compositions of the modeled SOA indicate that multifunctional low-volatility products contribute to isoprene SOA more significantly than previously thought, with a median mass contribution of ~ 57 % to the total modeled isoprene SOA. However, this contribution is intricately intertwined with IEPOX-derived SOA (IEPOX: isoprene-derived epoxydiols), posing challenges for their differentiation using bulk aerosol composition analysis (e.g., the aerosol mass spectrometer with positive matrix factorization). Furthermore, the SOA from these pathways may vary greatly, mainly dependent on the volatility estimation and treatment of particle-phase processes (i.e., photolysis and hydrolysis). Our findings emphasize that the various pathways to produce these low-volatility species should be considered in models to more accurately predict isoprene SOA formation. The new condensed isoprene chemical mechanism can be further incorporated into regional-scale air quality models, such as the Community Multiscale Air Quality Modelling System (CMAQ), to assess isoprene SOA formation on a larger scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A standardized workflow for kinetic metabolic model curation and dissemination

Kinetic metabolic models provide invaluable insights into cellular metabolism, supporting applications in synthetic biology, metabolic engineering, and systems biology. However, reproducibility and utility of these models hinge on clear and rigorous documentation, standardized annotation, and accessible visualization. This paper presents a workflow for building, annotating, visualizing, and sharing kinetic metabolic models. Our method integrates community standards and open-source tools to ensure reproducibility, interoperability, and user accessibility. This procedure enables researchers to produce reusable and well-documented kinetic models, advancing their role as powerful tools in metabolic research.

Cook, Margaret [Univ. of Washington, Seattle, WA (↗
DOE OSTI2025

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗
DOE OSTI2025

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗
DOE OSTI2024

Towards UV-Models of Kinetic Mixing and Portal Matter VII: A Light Dark Photon in the 3c3L1A1B Model

The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new U(1)DU(1)D​ gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter (DM) with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter (PM), also exist. The presence of such types of states and the strong suggestion of a need to embed U(1)DU(1)D​ into a non-abelian gauge structure not too far above the TeV scale based on the RGE running of the U(1)DU(1)D​ gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group G=SU(3)c×SU(3)L×U(1)A×U(1)B=3c3L1A1BG=SU(3)c​×SU(3)L​×U(1)A​×U(1)B​=3c​3L​1A​1B​ is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-abelian fashion that is not a simple product group of the form G=GSM×GDG=GSM​×GD​ encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.

Rizzo, Thomas G. [SLAC National Accelerator Labora↗
DOE OSTI2025

Toward UV models of kinetic mixing and portal matter. VII. A light dark photon in the 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 model

The kinetic mixing (KM) portal, by which the Standard Model (SM) photon mixes with a light dark photon arising from a new 𝑈⁢(1) 𝐷 gauge group, allows for the possibility of viable scenarios of sub-GeV thermal dark matter with appropriately suppressed couplings to the SM. This KM can only occur if particles having both SM and dark quantum numbers, here termed portal matter, also exist. The presence of such types of states and the strong suggestion of a need to embed 𝑈⁢(1) 𝐷 into a non-Abelian gauge structure not too far above the TeV scale based on the renormalization-group equation running of the 𝑈⁢(1)𝐷 gauge coupling is potentially indicative of an enlarged group linking together the visible and dark sectors. The gauge group 𝐺 = 𝑆⁢𝑈⁢(3) 𝑐 × 𝑆⁢𝑈⁢(3) 𝐿 × 𝑈⁢(1) 𝐴 × 𝑈⁢(1) 𝐵 = 3 𝑐 ⁢3 𝐿 ⁢1 𝐴 ⁢1 𝐵 is perhaps the simplest setup wherein the SM and dark interactions are partially unified in a non-Abelian fashion that is not a simple product group of the form 𝐺 = 𝐺 SM × 𝐺 𝐷 encountered frequently in earlier work. The present paper describes the implications and phenomenology of this type of setup.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A Kinetic Model-Driven Techno-Economic Analysis of Plastic Pyrolysis: Linking Process Dynamics to Economic Viability

This study employs a kinetic model integrated into Aspen Plus to predict pyrolysis product distribution under various conditions. A techno-economic assessment calculated the minimum selling price (MSP) of pyrolysis oil under different operating conditions for the baseline capacity of 100 kta, and across eight processing capacities ranging from 30 to 150 kta. The lowest MSP under the baseline capacity is estimated at $\$$420/ton, which is 33% lower than the 2023 average US crude oil price ($\$$74.6/bbl, equivalent to $\$$634/ton based on the density of pyrolysis oil). Under Monte Carlo simulation, accounting for variability in key economic and technical parameters, the mean MSP is estimated at $\$$1137/ton. The economic viability depends on feedstock price remaining below $\$$320/ton, defining the break-even feedstock price threshold. Sensitivity analysis further identifies capital investment and transportation cost as key economic drivers. Capacities beyond 90 kta show limited economies of scale benefits. Reducing product storage time cuts capital costs by 7% but raises operational risk. Uncertainty analysis suggests the economic feasibility of pyrolysis oil is unlikely to compete with crude oil without policy incentives.

petrochemicals↗