Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗
NASA NTRS2024

Hybrid Kinetic Modeling of the Magnetosheath Impulsive Plasma Cloud Penetration Through the Magnetopause and Comparison with MMS and Other Spacecraft Observations

This research examines the plasma processes under penetration of the plasma clouds (plasmoids) across the magnetopause which is modeled as a tangential discontinuity (TD). Cases with the parallel magnetic field in both sides out of the TD are under investigation. Plasma parameters and magnetic field were chosen from the MMS mission and other spacecraft observations. The results are important for understanding the following basic space plasma physics problems: (a) plasma cloud deformation and strong phase mixing with magnetospheric plasma; (b) the transfer of mass, momentum and energy of magnetosheath and magnetic cloud plasma into magnetospheric plasmas; (c) necessary conditions for plasma cloud penetration via the magnetopause; (d) wave generation by plasma clouds inside the magnetopause.

A S Lipatov↗
NASA NTRS1986

A model for the kinetics of a solar-pumped long path laser experiment

A kinetic model for a solar-simulator pumped iodine laser system is developed and compared to an experiment in which the solar simulator output is dispersed over a large active volume (150 cu cm) with low simulator light intensity (approx. 200 solar constants). A trace foreign gas which quenches the upper level is introduced into the model. Furthermore, a constant representing optical absorption of the stimulated emission is introduced, in addition to a constant representing the scattering at each of the mirrors, via the optical cavity time constant. The non-uniform heating of the gas is treated as well as the pressure change as a function of time within the cavity. With these new phenomena introduced into the kinetic model, a best reasonable fit to the experimental data is found by adjusting the reaction rate coefficients within the range of known uncertainty by numerical methods giving a new bound within this range of uncertainty. The experimental parameters modeled are the lasing time, laser pulse energy, and time to laser threshold.

Stock, L. V.↗
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 (↗
NASA NTRS2012

Partial Overhaul and Initial Parallel Optimization of KINETICS, a Coupled Dynamics and Chemistry Atmosphere Model

KINETICS is a coupled dynamics and chemistry atmosphere model that is data intensive and computationally demanding. The potential performance gain from using a supercomputer motivates the adaptation from a serial version to a parallelized one. Although the initial parallelization had been done, bottlenecks caused by an abundance of communication calls between processors led to an unfavorable drop in performance. Before starting on the parallel optimization process, a partial overhaul was required because a large emphasis was placed on streamlining the code for user convenience and revising the program to accommodate the new supercomputers at Caltech and JPL. After the first round of optimizations, the partial runtime was reduced by a factor of 23; however, performance gains are dependent on the size of the data, the number of processors requested, and the computer used.

KINETICS↗
NASA NTRS1989

Turbulent hydrocarbon combustions kinetics - Stochastic modeling and verification

Idealized reactors, that are designed to ensure perfect mixing and are used to generate the combustion kinetics for complex hydrocarbon fuels, may depart from the ideal and influence the kinetics model performance. A complex hydrocarbon kinetics model that was established by modeling a jet-stirred combustor (JSC) as a perfectly stirred reactor (PSR), is reevaluated with a simple stochastic process in order to introduce the unmixedness effect quantitatively into the reactor system. It is shown that the comparisons of the predictions and experimental data have improved dramatically with the inclusion of the unmixedness effect in the rich combustion region. The complex hydrocarbon kinetics is therefore verified to be mixing effect free and be applicable to general reacting flow calculations.

Wang, T. S.↗
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↗
NASA NTRS1981

Effects of turbulence on a kinetic auroral arc model

A plasma kinetic model of an inverted-V auroral arc structure which includes the effects of electrostatic turbulence is proposed. In the absence of turbulence, a parallel potential drop is supported by magnetic mirror forces and charge quasi neutrality, with energetic auroral ions penetrating to low altitudes; relative to the electrons, the ions' pitch angle distribution is skewed toward smaller pitch angles. The electrons energized by the potential drop form a current which excites electrostatic turbulence. In equilibrium the plasma is marginally stable. The conventional anomalous resistivity contribution to the potential drop is very small. Anomalous resistivity processes are far too dissipative to be powered by auroral particles. It is concluded that under certain circumstances equilibrium may be impossible and relaxation oscillations set in.

Cornwall, J. M.↗
DOE OSTI2024

Reaction Pathways and Energy Consumption in NH 3 Decomposition for H 2 Production by Low Temperature, Atmospheric Pressure Plasma

Pathways for NH 3 decomposition to N 2 and N 2 H 4 by atmospheric pressure nonthermal plasma are analyzed using a combination of molecular beam mass spectrometry measurements and zero-dimensional kinetic modeling. Experimental measurements show that NH 3 conversion and selectivity towards N 2 formation scale monotonically with the specific energy input into the plasma with ~ 100% selectivity to N 2 formation achieved at specific energy inputs above 0.12 J cm −3 (3.1 eV (molecule NH 3 ) −1 ). The kinetic model recovers these trends, although it underpredicts N 2 selectivity at low specific energy input. These discrepancies can be explained by the underestimation of reaction rate coefficients for reactions that consume N 2 H x species in collisions with H radicals and/or radial nonuniformities in power deposition, gas temperature, and species concentrations that are not represented by the plug flow approximation used in the model. The kinetic model shows that N 2 formation proceeds through N 2 H x decomposition pathways rather than NH x decomposition pathways in low temperature, atmospheric pressure plasma. Higher selectivity toward N 2 production can be achieved by operating at higher NH 3 conversion and with a higher gas temperature. Furthermore, the high energy cost of NH 3 decomposition by atmospheric pressure nonthermal plasma found in this work (25–50 eV (molecule NH 3 converted) −1 ; 17–33 eV (molecule H 2 formed) −1 ) is a result of the energy requirement for electron-impact dissociation of NH 3 and the significant re-formation of NH 3 by three-body recombination reactions between NH 2 and H.

Nonthermal plasma↗
DOE OSTI2025

Real-fluid behavior in rapid compression machines: Does it matter?

Rapid compression machines (RCMs) have been extensively used to quantify fuel autoignition chemistry and validate chemical kinetic models at high-pressure conditions. Historically, the analyses of experimental and modeling RCM autoignition data have been conducted based on the adiabatic core hypothesis with ideal gas assumption, where real-fluid behavior has been completely overlooked, though this might be significant at common RCM test conditions. Here, this work presents a first-of-its-kind study that addresses two significant but overlooked questions for autoignition studies within RCMs in the fundamental combustion community: (i) experiment-wise, can unaccounted-for real-fluid behavior in RCMs affect the interpretation and analysis of RCM experimental data? and (ii) simulation-wise, can unaccounted-for real-fluid behavior in RCMs affect RCM autoignition modeling and the validation of chemical kinetic models? To this end, theories for real-fluid isentropic change are newly proposed and derived based on high-order Virial EoS, and are further incorporated into an effective-volume real-fluid autoignition modeling framework newly developed for RCMs. With detailed analyses, the strong real-fluid behavior in representative RCM tests is confirmed, which can greatly influence the interpretation of RCM autoignition experiments, particularly the determination of end-of-compression temperature and evolution of the adiabatic core in the reaction chamber. Furthermore, real-fluid RCM modeling results reveal that considerable error can be introduced into simulating RCM autoignition experiments when following the community-wide accepted effective-volume approach by assuming ideal-gas behavior, which can be as high as 64% in the simulated ignition delay time at compressed pressure of 125 bar and lead to contradictory validation results of chemical kinetic models. Therefore, we recommend the community to adopt frameworks with real-fluid behavior fully accounted for (e.g., the one developed in this study) to analyze and simulate past and future RCM experiments, so as to avoid misinterpretation of RCM autoignition experiments and eliminate the potential errors that can be introduced into the simulation results with the existing RCM modeling frameworks.

High-order Virial equation of state↗
DOE OSTI2026

Kinetic and modeling studies of the mechanism of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2

Since its discovery over 15 years ago, the reversible dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 has remained one of the more intriguing hydrogen-cycling systems. While the mechanism of this reaction has been the subject of a good deal of speculation and computational studies, prior to this work it had not been probed through kinetic studies. Previous reports of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 have not included kinetic studies. The present studies have shown that the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 is suppressed by hydrogen pressure indicating that the rate-limiting step in this process involves hydrogen elimination. Computational modeling of kinetic data obtained from monitoring the hydrogen elimination from Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 under static vacuum over a range of temperatures supports that the dehydrogenation occurs through a reversible three-step process in which the elimination of hydrogen from the [B 3 H 10 ] − intermediate is rate limiting. A mechanism involving the low energy transfer of neighboring BH3 groups is proposed to account for the formation of [B 3 H 8 ] − at relatively low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Modeling of the HiPco process for carbon nanotube production. II. Reactor-scale analysis

The high-pressure carbon monoxide (HiPco) process, developed at Rice University, has been reported to produce single-walled carbon nanotubes from gas-phase reactions of iron carbonyl in carbon monoxide at high pressures (10-100 atm). Computational modeling is used here to develop an understanding of the HiPco process. A detailed kinetic model of the HiPco process that includes of the precursor, decomposition metal cluster formation and growth, and carbon nanotube growth was developed in the previous article (Part I). Decomposition of precursor molecules is necessary to initiate metal cluster formation. The metal clusters serve as catalysts for carbon nanotube growth. The diameter of metal clusters and number of atoms in these clusters are some of the essential information for predicting carbon nanotube formation and growth, which is then modeled by the Boudouard reaction with metal catalysts. Based on the detailed model simulations, a reduced kinetic model was also developed in Part I for use in reactor-scale flowfield calculations. Here this reduced kinetic model is integrated with a two-dimensional axisymmetric reactor flow model to predict reactor performance. Carbon nanotube growth is examined with respect to several process variables (peripheral jet temperature, reactor pressure, and Fe(CO)5 concentration) with the use of the axisymmetric model, and the computed results are compared with existing experimental data. The model yields most of the qualitative trends observed in the experiments and helps to understanding the fundamental processes in HiPco carbon nanotube production.

Validation Studies↗