Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A-Priori Tuning of Modified Magnussen Combustion Model

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for that particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor simulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for simulations.

Norris, A. T.↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
NASA NTRS1982

Ion distribution effects of turbulence on a kinetic auroral arc model

An inverted-V auroral arc structure plasma-kinetic model is extended to phenomenologically include the effects of electrostatic turbulence, with k-parallel/k-perpendicular being much less than unity. It is shown that, unless plasma sheet ions are very much more energetic than the electrons, anomalous resistivity is not a large contributor to parallel electrostatic potential drops, since the support of the observed potential drop requires a greater dissipation of energy than can be provided by the plasma sheet. Wave turbulence can, however, be present, with the ion cyclotron turbulence levels suggested by the ion resonance broadening saturation mechanism of Dum and Dupree (1970) being comparable to those observed on auroral field lines. The diffusion coefficient and net growth rate are much smaller than estimates based solely on local plasma properties.

Cornwall, J. M.↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Elucidating the Mechanism for the Reaction of $o$-Phthalaldehyde with Primary Amines in the Presence of Thiols

The use of o-phthalaldehyde (OPA) in combination with a thiol reagent is a common method for detecting primary amines in amino acids, peptides, and proteins. Despite its widespread use, the exact reaction mechanism has been debated since the 1980s. Here, we measure the kinetics of the reaction between OPA, alanine, and a dithiol (1,4-dithiolthreitol, DTT) as a function of pH and reagent concentration. Using these new measurements and accompanying kinetic models, we find evidence that the pH dependence of the kinetics arises from both the protonation states of alanine and DTT, the hydration state of OPA, and the unproductive equilibrium with DTT, all of which are pH-dependent. Importantly, these results support the mechanism originally proposed by Sternson [Rational design and evaluation of improved o-phthalaldehyde-like fluorogenic reagents. Anal. Biochem. 1985, 144, 233–246] and Wong [Reaction of o-phthalaldehyde with alanine and thiols: kinetics and mechanism. J. Am. Chem. Soc. 1985, 107, 6421–6422], in which the primary amine first reacts with OPA, followed by a reaction with the thiol to form the fluorescent isoindole product.

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

A cure-rate model for the Shuttle filament-wound case

An epoxy and carbon fiber composite has been used to produce a light-weight rocket case for the Space Shuttle. A kinetic model is developed which can predict the extent of epoxy conversion during the winding and curing of the case. The model accounts for both chemical and physical kinetics. In the model, chemical kinetics occur exclusively up to the time the transition temperature equals the reaction temperature. At this point the resin begins to solidify and the rate of this process limits the rate of epoxy conversion. A comparison of predicted and actual epoxy conversion is presented for isothermal and temperature programmed cure schedules.

Cagliostro, D. E.↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Mixing and non-equilibrium chemical reaction in a compressible mixing layer

The effects of compressibility, chemical reaction exothermicity, and non-equilibrium chemical modeling in a reacting plane mixing layer were investigated by means of two dimensional direct numerical simulations. The chemical reaction was irreversible and second order of the type A + B yields Products + Heat. The general governing fluid equations of a compressible reacting flow field were solved by means of high order finite difference methods. Physical effects were then determined by examining the response of the mixing layer to variation of the relevant non-dimensionalized parameters. The simulations show that increased compressibility generally results in a suppressed mixing, and consequently a reduced chemical reaction conversion rate. Reaction heat release was found to enhance mixing at the initial stages of the layer growth, but had a stabilizing effect at later times. The increased stability manifested itself in the suppression or delay of the formation of large coherent structures within the flow. Calculations were performed for a constant rate chemical kinetics model and an Arrhenius type kinetic prototype. The choice of the model was shown to have an effect on the development of the flow. The Arrhenius model caused a greater temperature increase due to reaction than the constant kinetic model. This had the same effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗
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 OSTI2021

Plasma-assisted catalysis for ammonia synthesis in a dielectric barrier discharge reactor: key surface reaction steps and potential causes of low energy yield

Ammonia synthesis experiments were carried out in a coaxial dielectric barrier discharge (DBD) reactor packed with several different supports and metal catalysts. There was a marked increase in the reaction rate, over that obtained in an empty DBD plasma reactor, upon introduction of a packed bed of γ-Al 2 O 3 , Ru/γ-Al 2 O 3 , or SiO 2 particles. The difference in the reaction rates over γ-Al 2 O 3 and Ru/γ-Al 2 O 3 was minimal. Complementary zero-dimensional plasma kinetic model analysis was also performed using inputs from experimental data. This kinetic analysis allowed for gas phase reactions, Eley–Rideal (E–R) reactions, and direct adsorption of radical species on the γ-Al 2 O 3 surface. On the metal surface, dissociative adsorption of N 2 and H 2 , and Langmuir–Hinshelwood reactions were also included. This analysis revealed that, under the conditions of our experiments, ammonia synthesis proceeds principally by the formation of reactive radicals in the gas phase, which then adsorb and participate in E–R reactions on both the metal and support material surfaces. Furthermore, this finding illustrates a challenge for substantially increasing the energy yield for plasma-assisted ammonia synthesis in typical DBD reactors containing packed catalyst beads.

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

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

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

Observation and Modeling of the Thermal Runaway of High-Capacity Pouch Cells Due to an Internal Short Induced by an Indenter

Owing to the complexity of coupling mechanical and electrical solvers for finite element analysis (FEA) of mechanical abuse of a battery, the simulation should be performed by incorporating a local short circuit model activated based on the status of battery elements. However, the approach should also consider various other conditions of the battery simulated, such as the initial state of charge (SOC), for precise simulations of the internal short. This study proposes an approach developed with a mechanical model and single-particle battery model, along with a reaction kinetics model for thermal abuse using LS-DYNA. Several key parameters were selected, such as the activation criteria, kinetic abuse and element integration points during the internal short circuit induced by the indenter. The impact of these key parameters in the FEA on thermal runaway response for various initial SOCs from 25 to 100 % were analyzed. The FEA results obtained were compared against experimental results with various parameters such as activation criterion of internal short, kinetic abuse based on accelerating rate calorimetry (ARC) data. Finally, the root causes of the discrepancy between the FEA and experimental results are discussed regarding integration point what solver can make difference during calculation.

ENERGY STORAGE↗
NASA NTRS1995

Composition/Structure/Dynamics of comet and planetary satellite atmospheres

This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.

Combi, Michael R.↗
DOE OSTI2020

Inorganic Sulfur Species Formed upon Heterogeneous OH Oxidation of Organosulfates: A Case Study of Methyl Sulfate

In this report recent studies reveal that organosulfates at the particle surface can be oxidized by gas-phase OH radicals with significant rates. Inorganic sulfur species, such as the bisulfate ion (HSO 4 - ) and sulfate ion (SO 4 2- ), can be formed upon these heterogeneous oxidation processes through the formation and subsequent reactions of sulfate radical anions (SO 4 •- ) in the particle phase. However, the amount of inorganic sulfur species produced in these heterogeneous oxidation reactions is not known. We investigate the heterogeneous OH oxidation of sodium methyl sulfate (CH 3 SO 4 Na), the smallest organosulfate detected in atmospheric particles, using an oxidation flow reactor at a relative humidity of 75%. We quantify the kinetics by measuring the decay of CH 3 SO 4 Na and the amount of HSO 4 - and SO 4 2- formed upon oxidation using ion chromatography. Kinetic measurements determine the heterogeneous OH reaction rate to be (5.72 ± 0.14) x 10 -13 cm 3 molecule -1 s -1 , with an effective OH uptake coefficient, γ eff , of 0.31 ± 0.06. The molar yield of inorganic sulfur species, defined as the total number of moles of HSO 4 - and SO 4 2- formed per mole of CH 3 SO 4 Na consumed upon oxidation, is found to be significant and has an average value of 0.62 ± 0.18 upon oxidation. A kinetic model is developed to describe the kinetics and inorganic sulfur species formation upon oxidation. Model simulations suggest that CH 3 SO 4 Na tends to decompose rapidly into formaldehyde and SO 4 •- , and the reaction of SO 4 •- with CH 3 SO 4 Na plays a significant role in both governing the kinetics and the formation of inorganic sulfur species.

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

Partial equilibrium model for predicting concentration of CO in combustion

A simple kinetic model, based upon the concept of partial equilibrium, is developed for predicting carbon monoxide concentrations in such steady-flow hydrocarbon-air combustion systems as gas-turbine combustors. The only two kinetic constraints used in the model are on the rate of change of the total number of gaseous particles in the system, and on the rate of change of the CO concentration. The accuracy of the model is verified by comparison with experimental results obtained by burning kerosene with air in a 130,000 Btu/hr atmospheric-pressure steady-flow burner. Burnt-gas cooling rates of 1,000,000 K/sec, about the magnitude encountered in the dilution zone of gas-turbine combustors, were obtained with a compact water-cooled heat exchanger. Good agreement between measured CO concentrations and the values predicted by the partial-equilibrium model was obtained.

Morr, A. R.↗
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↗