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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts

Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.

42 - ENGINEERING↗
DOE OSTI2024

Thermodynamic control on the decomposition of organic matter across different electron acceptors

The increasing availability of high-resolution characterization of natural organic matter (OM) data has shifted the paradigm of lumped descriptions of OM components and potential microbial activities. Our recent development of a substrate-explicit thermodynamic model uniquely enables incorporating complex OM pools to formulate biogeochemical reaction models based on their elemental compositions. While this previous work facilitates prediction of aerobic respiration of complex OM, it is equally imperative to consider the role of non-oxygenic electron acceptors in regulating OM turnover and the fate of carbon. In this study, we significantly expand our previous model by flexibly incorporating both detailed OM chemistry and electron acceptors other than oxygen. Here, our modeling analysis has revealed substantial variations in the energy status of OM molecules across different soils, which drive the co-occurrence of different electron-accepting processes. We demonstrated the effectiveness of the proposed model using a consistency check with experimental data. Through systematic evaluation of the impact of diverse chemical inputs (both electron donors and acceptors) on OM decomposition, the new model also revealed how key microbial growth parameters such as carbon use efficiency (CUE) and reaction rates vary across different electron-accepting processes. Our model provides a unified framework integrating thermodynamic and kinetic constraints on microbial metabolic activities. It complements traditional kinetic models, which are often designed solely to capture mass fluxes. We conclude that thermodynamic modeling emerges as a powerful tool for describing the mechanisms underlying the interplay between microbial growth and OM chemistry and cycling across different electron acceptors, enhancing our ability to project complex ecosystem behaviors in dynamic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

High-pressure CO 2 dissociation with nanosecond pulsed discharges

Herein we investigate the conversion of CO 2 into CO and O 2 with nanosecond repetitively pulsed (NRP) discharges in a high-pressure batch reactor. Stable discharges are obtained at up to 12 bar. By-products are measured with gas chromatography. The energy efficiency is determined for a range of processing times, pulse energy, and fill pressures. It is only weakly sensitive to the plasma operating parameters, i.e the extent of CO 2 conversion is almost linearly-dependent on the specific energy invested. A conversion rate as high as 14% is achieved with an energy efficiency of 23%. For long processing times, saturation in the yield and a drop in efficiency are observed, due to the increasing role of three-body recombination reactions, as described by zero-dimensional detailed kinetic modeling. The modeling reveals the presence of three-stage kinetics between NRP pulses, controlled by electron-impact CO 2 dissociation, vibrational relaxation, and neutral elementary kinetics. Transport effects are shown to be important for CO 2 conversion at high pressures. For fill pressures beyond 10 bar, CO 2 may locally transit into supercritical states. The supercritical plasma kinetics may bypass atomic oxygen pathways and directly convert CO 2 into O 2 . This work provides a detailed analysis of plasma-based high-pressure CO 2 conversion, which is of great relevance to future large-scale sustainable carbon capture, utilization, and storage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

LOKI

LOKI is a simulation code that models weakly-coupled plasma evolution using Eulerian kinetic models in four phase space dimensions. The code is useful to understand basic plasma nonlinear plasma physics phenomena - in particular, wave-particle and wave-wave resonances that occur in laser plasma interactions. LOKI can simulate fundamental physical process in low-density plasmas on a scale of hundreds of light wavelengths. It does so by solving the relativistic Vlasov and Valsov-Fokker-Planck differential transport equations describing the collision less and low-collisionality evolution of charged particles in electromagnetic fields in two configuration space and two velocity space dimensions. Two principle electromagnetic models are provided: electrostatic (Poisson's Equation) and electromagnetic (Maxwell's equations).

DIMITS, ANDRIS↗
DOE OSTI2024

SHOCK TUBE IGNITION STUDIES OF RENEWABLE DIESEL FUELS FOR MEDIUM AND HEAVY-DUTY TRANSPORTATION

Currently extensive research on alternative fuels is being conducted due to their increasing demand to reduce greenhouse emissions. One renewable fuel studied in this work is dimethyl ether (DME) blended with propane(C3H8) as a potential mixture for heavy-duty engines used in semi-trucks. The blend has the potential to drastically reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. To develop the use of mixture, one must conduct detailed conceptual and simulation studies before progressing to detail studies in CFD, engine modifications, and live testing. For simulations, accurate high-fidelity chemical kinetic models are necessary. However, the validity of the chemical kinetic mechanism for operating conditions of a heavy-duty mixing-controlled compression (MCCI) engine was widely unknown until recent work presented here and published. In this work, we studied the ignition of DME and propane blends in a shock tube under MCCI engine conditions. Ignition delay time (IDT) gathered behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Testing was conducted for undiluted varieties spanning from temperatures of 700 to 1100 K at pressures ranging from 55 to 84 bar for various blends (100% CH3OCH3, 100% C3H8, 60% CH3OCH3/ 40% C3H8) of DME and propane were combusted in synthetic air (21% O2/ 79% N2). Several experiments were conducted at higher pressures (90-120 bar) to improve the model performance and accuracy. The ignition delay times (IDTs) were compared to recent mechanisms, including Aramco3.0, NUIG, and Dames et al. A common trend among the mechanisms was overpredicted experimental IDTs. Further studies were conducted by a sensitivity analysis using the Dames et al. model, and critical reactions sensitive to IDTs of DME-propane mixture near 60 bar are outlined. Chemical analysis was conducted on the NTC region to explain chemical kinetics which is critical for developing MCCI heavy duty engines.

Mohammed, Zuhayr Pasha↗
DOE OSTI2022

2D Reconstruction of Magnetotail Electron Diffusion Region Measured by MMS

Models for collisionless magnetic reconnection in near-Earth space are distinctly characterized as 2D or 3D. In 2D kinetic models, the frozen-in law for the electron fluid is usually broken by laminar dynamics involving structures set by the electron orbit size, while in 3D models the width of the electron diffusion region is broadened by turbulent effects. We present an analysis of in situ spacecraft observations from the Earth's magnetotail of a fortuitous encounter with an active reconnection region, mapping the observations onto a 2D spatial domain. While the event likely was perturbed by low-frequency 3D dynamics, the structure of the electron diffusion region remains consistent with results from a 2D kinetic simulation. As such, the event represents a unique validation of 2D kinetic, and laminar reconnection models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Investigating the Adsorption–Desorption Kinetics of a Molecular Water Oxidation Catalyst at an Electrode Interface

Probing the dynamics of molecular catalysts at electrode–electrolyte interfaces is essential for understanding catalytic mechanisms. Structure-specific spectroscopic methods are particularly powerful for examining electrocatalytic interfaces but are mostly used under steady-state conditions. Herein, we combined surface-enhanced infrared absorption spectroscopy (SEIRAS) with phase-sensitive detection (PSD) to investigate the dynamics of a molecular Ir-based water oxidation catalyst at the Au–electrolyte interface. We found that the amplitude of the absorbance of the catalyst is anticorrelated to that of interfacial water. This anticorrelation can be understood by the adsorption of the electrooxidized catalyst on the electrode and concurrent displacement of interfacial water. The infrared signals from the interface exhibit an increasing phase lag with respect to the electrode potential with an increasing scan rate of the potential. Kinetic modeling suggests that the potential-dependent adsorption–desorption kinetics of the molecular catalyst on the electrode gives rise to this phase lag. Furthermore, this study shows that PSD-SEIRAS is a powerful tool for investigating the interfacial dynamics of electrocatalytic systems.

Catalysts↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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

Heat Integration Optimization and Dynamic Modeling Investigation for Advancing the Coal-Direct Chemical Looping Process

The purpose of the project is to address the optimization and startup operation of a modular coal direct chemical looping (CDCL) combustion system integrated with a steam cycle for power generation to reduce the risks involved in further scale-up of the technology. The modular reactor design of the CDCL process provides flexibility in the fabrication of the reactor and in its operating capacity (i.e. turndown ratio) at the cost of a more complex heat exchange network (HEN) design and integration. To address the technology gaps and advance the efficiency and economic feasibility of the CDCL technology, the project will perform a detailed and comprehensive analysis of the integration of a modular CDCL reactor system and a steam cycle system under both static and transient conditions via HEN process performance simulations and system dynamic modeling, respectively. The scope of work consists of 1) Experimental and computational studies of the CDCL combustor reactor 2) Comprehensive static (i.e. steady-state) system HEN design analysis in CDCL 550 MWe commercial unit for power generation and 3) Dynamic modeling of site specific design of 10MWe CDCL large pilot plant. The project team has successfully developed and validated a kinetic model for the oxidation of oxygen carriers in the combustor using the unreacted shrinking core model (UCSM). The model is capable of capturing the oxidation kinetics of fully or partially reduced oxygen carrier particles. A computational fluid dynamics (CFD) model is developed to simulate the hydrodynamics, heat transfer, and chemical reaction occurring in the CDCL combustor. The model is developed in MFIX and ANSYS Fluent. Key aspects of CDCL combustor operation, including heat transfer, oxygen carrier oxidation, and the transport of oxygen carrier particles, are simulated using this CFD model. The HEN for a commercial scale 550 MWe CDCL power plant is simulated and optimized using ASPEN Plus. Practical design considerations are incorporated based on industrial experiences. The performance and cost for the commercial CDCL plant is updated based on these analyses. A dynamic model for the 10 MWe CDCL pilot plant is developed in ProTRAX simulation software. The model is based on the pilot plant design developed in project DE-FE0027654 “10 MWe CDCL Large Pilot Plang – Pre-FEED Study” and the steam cycle data obtained from Dover Light & Power plant. The transient behaviors during pilot plant load variation are simulated using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗
DOE OSTI2025

Active Site Dynamics in Molybdenum-Based Silica-Supported Olefin Metathesis Catalysts: Site Renewal and Decay Beyond the Chauvin Cycle

Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.

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

New physics in multi-lepton tau decays

Dark particles with lepton-flavor-violating couplings to the tau lepton can induce rare neutrinoless $τ$ decays with large final state multiplicities. We study models where transitions of the type $τ^\pm\to \ell^\pm\,ϕ$, with $ϕ$ a light new particle, initiate a chain of decays in the dark sector that terminate with decays into electrons, muons, or pions. These decay cascades appear as rare five or even seven-body $τ$ decays with multiple reconstructable resonances. We survey several representative models: kinetically mixed dark photon, gauged $L_i-L_j$ models, and other more exotic charge assignments such as chiral $U(1)'$ extensions of the Standard Model. The main new ingredient is the possibility of flavor violation at very high scales. In these models, a number of channels that have not yet been searched for experimentally, such as $τ\to 5μ$, $τ\to 3μ\,2e$, $τ\to μ\,4e$, and hadronic channels like $τ\to μ\,4π$, typically dominate over the previously-considered signatures such as $τ\to 3μ$. While some of the models, such as the gauged $L_i-L_j$ ones, also contain more challenging channels with missing energy due to decays to neutrinos, they can still be searched for via fully visible channels.

Ema, Yohei [Florida U.]↗
DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Efficient kinetic thermal inverse modeling for organic material decomposition

The prevalent use of organic materials in manufacturing is a fire safety concern, and motivates the need for predictive thermal decomposition models. A critical component of predictive modeling is numerical inference of kinetic parameters from bench scale data. Currently, an active area of computational pyrolysis research focuses on identifying efficient, robust methods for optimization. This paper demonstrates that kinetic parameter calibration problems can successfully be solved using classical gradient-based optimization. We explore calibration examples that exhibit characteristics of concern: high nonlinearity, high dimensionality, complicated schemes, overlapping reactions, noisy data, and poor initial guesses. The examples demonstrate that a simple, non-invasive change to the problem formulation can simultaneously avoid local minima, avoid computation of derivative matrices, achieve a computational efficiency speedup of 10x, and make optimization robust to perturbations of parameter components. Techniques from the mathematical optimization and inverse problem communities are employed. By re-examining gradient-based algorithms, we highlight opportunities to develop kinetic parameter calibration methods that should outperform current methods.

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

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗