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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluating Possible Formation Mechanisms of Criegee Intermediates during the Heterogeneous Autoxidation of Squalene

Organic molecules in the environment oxidatively degrade by a variety of free radical, microbial, and biogeochemical pathways. A significant pathway is heterogeneous autoxidation, in which degradation occurs via a network of carbon and oxygen centered free radicals. Recently, we found evidence for a new heterogeneous autoxidation mechanism of squalene that is initiated by hydroxyl (OH) radical addition to a carbon–carbon double bond and apparently propagated through pathways involving Criegee Intermediates (CI) produced from β-hydroxy peroxy radicals (β-OH-RO 2 •). It remains unclear, however, exactly how CI are formed from β-OH-RO 2 •, which could occur by a unimolecular or bimolecular pathway. Combining kinetic models and multiphase OH oxidation measurements of squalene, we evaluate the kinetic viability of three mechanistic scenarios. Scenario 1 assumes that CI are formed by the unimolecular bond scission of β-OH-RO 2 •, whereas Scenarios 2 and 3 test bimolecular pathways of β-OH-RO 2 • to yield CI. Scenario 1 best replicates the entire experimental data set, which includes effective uptake coefficients vs [OH] as well as the formation kinetics of the major products (i.e., aldehydes and secondary ozonides). Finally, although the unimolecular pathway appears to be kinetically viable, future high-level theory is needed to fully explain the mechanistic relationship between CI and β-OH-RO 2 • in the condensed phase.

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

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

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

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Integrative mapping reveals molecular features underlying the mechanism of nucleocytoplasmic transport

Nuclear pore complexes (NPCs) enable rapid, selective, and robust nucleocytoplasmic transport. To explain how transport emerges from the system components and their interactions, we used experimental data and theoretical information to construct an integrative Brownian dynamics model of transport through an NPC, coupled to a kinetic model of transport in the cell. The model recapitulates key aspects of transport for a wide range of molecular cargoes, including preribosomes and viral capsids. Our model quantifies how flexible phenylalanine-glycine (FG) repeat proteins create an entropic barrier to passive diffusion and how this barrier is selectively lowered in facilitated diffusion by the many transient interactions of nuclear transport receptors with the FG repeats. Selective transport is enhanced by “fuzzy” multivalent interactions, redundant FG repeat mass, coupling to the energy-dependent RanGTP concentration gradient, and exponential dependence of transport kinetics on the transport barrier. Our model will facilitate rational modulation of the NPC and its artificial mimics.

Raveh, Barak [School of Computer Science and Engin↗
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 OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2024

Catalytic production of δ-valerolactone (DVL) from biobased 2-hydroxytetrahydropyran (HTHP) – Combined experimental and modeling study

δ-Valerolactone (DVL) is a five-carbon (C5) cyclic ester that can undergo ring-opening polymerization to yield high-performance, biocompatible polyesters. But current market prices of C5 chemicals like DVL are very high due to poor availability of C5 feedstock in petroleum. Herein, we demonstrate a novel route to DVL synthesis via dehydrogenation of biomass-derived 2-hydroxytetrahydropyran (HTHP) over Cu/SiO 2 without the use of toxic reagents. Since HTHP exists in thermal equilibrium with 3,4-dihydropyran (DHP) via dehydration, and with 2,2’-oxybis(tetrahydropyran) and 5-(tetrahydropyran-2-yloxy)pentanal via acetalization, we have also determined the thermochemistry (ΔH rxn and ΔG rxn ) of each competing reaction using density functional theory (DFT) calculations at the M06-2X/cc-pVTZ level. Finally, by developing a kinetic model of all 8 reactions involved, we have achieved 84% selectivity to DVL at 150°C in a packed bed reactor for over 72 hours of time on stream.

42 ENGINEERING↗
DOE OSTI2025

Emergence and disruption of cooperativity in a denitrifying microbial community

Anthropogenic perturbations to the nitrogen cycle, primarily through use of synthetic fertilizers, is driving an unprecedented increase in the emission of nitrous oxide (N 2 O), a potent greenhouse gas and an ozone depleting substance, causing urgency in identifying the sources and sinks of N 2 O. Microbial denitrification is a primary contributor to biotic production of N 2 O in anoxic regions of soil, marine systems, and wastewater treatment facilities. Here, through comprehensive genome analysis, we show that pathway partitioning is a ubiquitous mechanism of complete denitrification within microbial communities. We have investigated mechanisms and consequences of process partitioning of denitrification through detailed physiological characterization and kinetic modeling of a synthetic community of Rhodanobacter thiooxydans FW510-R12 and Acidovorax sp. GW101-3H11. We have discovered that these two bacterial isolates, from a heavily nitrate (NO 3 − ) contaminated superfund site, complete denitrification through the exchange of nitrite (NO 2 − ) and nitric oxide (NO). The process partitioning of denitrification and other processes, including amino acid metabolism, contribute to increased cooperativity within this denitrifying community. We demonstrate that certain contexts, such as high NO 3 − , cause unbalanced growth of community members, due to differences in their substrate utilization kinetics. The altered growth characteristics of community members drives accumulation of toxic NO 2 − , which disrupts denitrification causing N 2 O off gassing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Catalytic Reduction of Esters over Zirconia-Supported Metal Catalysts

Esters are often produced as unwanted byproducts during the catalytic upgrading of ethanol to diesel fuel precursors through Guerbet coupling. Removal of esters from the product stream is important to prevent the loss of downstream catalyst activity from ester-derived carboxylic acids. In this work, we studied ester hydrogenolysis to the parent alcohols as a viable route for enhanced diesel fuel production. Specifically, we investigated the reduction of hexyl acetate in butanol over ZrO 2 -supported Ni, Co, Cu, Rh, Pd, and Pt catalysts, where Cu/ZrO 2 was the most selective catalyst for the hydrogenolysis of hexyl acetate into hexanol and ethanol. Thermodynamic analysis reveals that a 90% alcohol yield can be obtained at 200 °C, 30 bar, and a relatively high H 2 :hexyl acetate molar ratio of 480:1. Experimentally, an alcohol yield of 88% yield was obtained with a 10 wt % Cu/ZrO 2 catalyst at these conditions with a residence time of 5.4 h kg cat kmol gas –1 . Catalytic tests on the support revealed that ZrO 2 catalyzes the transesterification reaction between hexyl acetate and butanol. However, only the Cu sites can catalyze the hydrogenolysis of the esters into the final alcohols. We developed a kinetic model for our experimental results, which shows that the transesterification and hydrogenolysis reactions run at two different timescales, the former being 10 times faster than the latter. Data regression has been used to develop a model to predict the mole fraction distribution of ester hydrogenolysis products over a wide range of contact times. Cu/ZrO 2 loses half its catalytic activity after 80 h of time on stream. Modeling of deactivation data reveals that the ZrO 2 support conserves a residual activity due to external active sites, while active sites over the Cu surface deactivate at different rates. Furthermore, the catalytic conversion of esters into their parent alcohols is relevant to the production of surrogate liquid fuels since alcohols can be bimolecularly dehydrated to produce a blend of ethers with diesel fuel-like properties.

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

Axisymmetric hybrid Vlasov equilibria with applications to tokamak plasmas

We derive axisymmetric equilibrium equations in the context of the hybrid Vlasov model with kinetic ions and massless fluid electrons, assuming isothermal electrons and deformed Maxwellian distribution functions for the kinetic ions. The equilibrium system comprises a Grad–Shafranov partial differential equation and an integral equation. These equations can be utilized to calculate the equilibrium magnetic field and ion distribution function, respectively, for given particle density or given ion and electron toroidal current density profiles. The resulting solutions describe states characterized by toroidal plasma rotation and toroidal electric current density. Additionally, due to the presence of fluid electrons, these equilibria also exhibit a poloidal current density component. This is in contrast to the fully kinetic Vlasov model, where axisymmetric Jeans equilibria can only accommodate toroidal currents and flows, given the absence of a third integral of the microscopic motion.

Physics↗
DOE OSTI2024

Kinetic Study of Polyvinyl Chloride Pyrolysis with Characterization of Dehydrochlorinated PVC

In this paper, we study the kinetics of polyvinyl chloride (PVC) decomposition using a combination of experimental and computational approaches. Here, we develop a simplified kinetic model that contains only two steps: dehydrochlorination of PVC and further decomposition of the PVC residue. The model is consistent with density functional theory (DFT) calculations and experimental data. Dehydrochlorination is an autocatalytic reaction that starts with a tertiary chloride (Cl) and generates hydrogen chloride (HCl) and benzene as the main products. Benzene and HCl formation rates showed similar trends, indicating that HCl likely catalyzes a homolytic carbon-carbon (C-C) bond cleavage, which gives rise to benzene and an aliphatic fragment. We characterized the structure of dehydrochlorinated PVC (PVC residue) by using thermal gravimetric analysis (TGA), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance spectroscopy (NMR). FTIR and NMR results indicate that the PVC residue contains 20% quaternary carbon content, indicating a high concentration of cross-linked molecules. We predict that the most probable structure in the cross-linked centers of the PVC residue is cyclohexadiene, which is supported by DFT calculations, FTIR, and NMR.

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

Bisphenol A degradation by manganese oxides at circumneutral pH: Quantitative evaluation of dissolved Mn(III) species with pyrophosphate

Although trivalent manganese (Mn(III)) species have been recognized as crucial intermediates in the degradation of organic contaminants by Mn oxides, quantitative research on their specific roles remains scarce. Here, our study investigated the degradation processes of an organic pollutant, Bisphenol A (BPA), by dissolved Mn(III) and Mn(III)-bearing oxides, and elucidated the differences of the underlying mechanisms and reaction pathways between several Mn oxides and dissolved Mn(III). Our results indicated that BPA degradation rates with Mn(III)-bearing oxides alone follow the order: δ-MnO 2 >> γ-MnOOH > Mn 3 O 4 . Adding pyrophosphate (PP) significantly enhanced BPA degradation by promoting the formation of Mn(III)-PP complexes and exposing more reactive sites, achieved through destabilizing the crystal structure and mitigating of Mn(II) readsorption, particularly in γ-MnOOH and Mn 3 O 4 . Our kinetic model revealed that heterogeneous degradation by Mn oxides is the predominant reaction pathway, accounting for 61.4 %, 87.8 %, and 73.8 % of the total degraded BPA for δ-MnO 2 , γ-MnOOH, and Mn 3 O 4 , respectively, even in the presence of significant amount of dissolved Mn(III) intermediates due to high PP concentrations. These results offer mechanistic details on BPA degradation by Mn oxides and the influence of ligand concentration, providing helpful insights for optimizing degradation strategies of organic pollutants.

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

External perturbation-driven Sabatier breakthrough

The Sabatier reaction (CO 2 + 4H 2 → CH 4 + 2H 2 O) is gaining renewed interest due to its potential to reduce energy carrier storage costs, serve as a feedstock for various organic chemicals, and supply in-situ propellant and life-support resources for long-duration Mars missions. This study demonstrates that combining a modest 2 mA electric field with H 2 feed modulation markedly elevates the CO 2 hydrogenation activity of 2 wt% Ru/CeO 2 catalyst. CO 2 conversion reaches 88 % and 93 % with a CH 4 yield of 83 % and 89 % at 350 °C and 450 °C, respectively. A simple lumped kinetic model reveals that the combined external perturbations not only shift the reaction mechanism but also redistribute key surface-adsorbed intermediates such as hydrogen adatoms and hydrogen-activated CO 2 among the Ru clusters, Ru/CeO 2 interface, and ceria surface. The electric field accelerates the conversion of adsorbed CO 2 to the hydrogenated CO 2 species on Ru and boosts CH 4 formation rate constant, while simultaneously suppresses the formation of undesired, non-reactive surface intermediates. Degree-of-rate-control analysis pinpoints proton migration across the metal-support interface as the decisive lever under these coupled perturbations. In conclusion, these findings establish that rational pairing of metal-support design with well-tuned electric fields and feed oscillations can unlock unprecedented Sabatier rates, guiding the development of next-generation reactors for efficient CO 2 to CH 4 conversion.

10 - SYNTHETIC FUELS↗