Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Finite-dimensional collisionless kinetic theory

A collisionless kinetic plasma model may often be cast as an infinite-dimensional noncanonical Hamiltonian system. Here I show that, when this is the case, the model can be discretized in space and particles while preserving its Hamiltonian structure, thereby producing a finite-dimensional Hamiltonian system that approximates the original kinetic model. I apply the general theory to two example systems: the relativistic Vlasov-Maxwell system with spin and a gyrokinetic Vlasov-Maxwell system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Analysis of a general circulation model. II - Distribution of kinetic energy in the South Atlantic and Kuroshio/Oyashio systems

It was found (Garzoli et al., 1992) that the general circulation model of Semtner and Chervin (1992) provides accurate descriptions of the Brazil-Malvinas and the Kuroshio/Oyashio confluence systems, except for the fact that the model prediction shows less variability than that present in observations. This paper investigates the problem of model variability by analyzing the mean and the eddy kinetic energy from the model and comparing the values with the Geosat altimeter observations for the South Atlantic Ocean and for the Kuroshio system. It is found that, while the model shows transient eddy activity in the areas that overlap the Geosat observations, the energy level of the model transient motions is considerably smaller following an arch along the bottom topography. The same was found from the comparisons made with values obtained from FGGE and surface drifters. It is suggested that the model is poorly resolving instabilities in the confluence front, and is not resolving other transients appearing in regions of marked topography.

Garraffo, Zulema↗
DOE OSTI2018

Phase-field model of pitting corrosion kinetics in metallic materials

This study presents a thermodynamically consistent phase field model for the quantitative prediction of the pitting corrosion kinetics in metallic materials. An order parameter is introduced to represent the metal’s physical state at each point in the entire domain. The free energy of the metal-electrolyte system is defined in terms of its metal ion concentration and the order parameter. The mass transport in the electrolyte and the electrochemical reactions in the electrolyte/metal interface are taken into consideration in the model to simulate metal corrosion in a corrosive environment. The governing equations for mass transport and the order parameter are solved in such a manner that the system’s free energy is reduced as a result of diffusion and migration processes, and the distribution of electrostatic potential is governed by Poisson’s equation. A calibration study is performed to couple the kinetic interface parameter with the corrosion current density, which results in a direct relationship between overpotential and the kinetic interface parameter. A comparative study is performed to validate the phase field model against the experimental results. Several case studies are also examined to understand the corrosion behavior of closely located pits, stressed material, ceramic particles-reinforced steel, and crystallographic orientation dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The dissolution behavior of borosilicate glasses in far-from equilibrium conditions

In far-from-equilibrium conditions, the dissolution of borosilicate glasses used to immobilize nuclear waste is known to be a function of both temperature and pH. The aim of this paper is to study effects of these variables on three model waste glasses (SON68, ISG, AFCI). To do this, experiments were conducted at temperatures of 23, 40, 70, and 90 °C and pH(RT) values of 9, 10, 11, and 12 with the single-pass flow-through (SPFT) test method. The results from these tests were then used to parameterize a kinetic rate model based on transition state theory. Both the absolute dissolution rates and the rate model parameters are compared with previous results. Discrepancies in the absolute dissolution rates as compared to those obtained using other test methods are discussed. Rate model parameters for the three glasses studied here are nearly equivalent within error and in relative agreement with previous studies. The results were analyzed with a linear multivariate regression (LMR) and a nonlinear multivariate regression performed with the use of the Glass Corrosion Modeling Tool (GCMT), which is capable of providing a robust uncertainty analysis. This robust analysis highlights the high degree of correlation of various parameters in the kinetic rate model. As more data are obtained on borosilicate glasses with varying compositions, the effect of glass composition on the rate parameter values could possibly be obtained. This would allow for the possibility of predicting the forward dissolution rate of glass based solely on composition

Neeway, James J.↗
DOE OSTI2025

Theoretically Informed Kinetics (ThInK): Establishing a modern C 0 -C 3 mechanism for combustion modeling

In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.

Energy transfer↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2022

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

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

Geologic hydrogen as an emerging fuel: experimental insights, thermodynamics, kinetics, and reactive transport modeling

Geologic hydrogen (GeoH 2 ) is emerging as a viable clean energy source. It is largely produced through serpentinization, a geological process in which ultramafic rocks react with water under suitable temperature and pressure. Here, this review synthesizes the current understanding of H 2 generation by serpentinization, with an emphasis on reaction mechanisms, kinetics, and thermodynamics, as well as on modeling flow and transport of reacting fluids in geological formations. We describe the role of mineral assemblages, such as olivine and pyroxene, fluid-rock interactions, and catalytic surfaces, in influencing GeoH 2 yield and reaction rates. By integrating models of reaction kinetics, subsurface reactive flow and transport, and the serpentinization process, and by accounting for the thermodynamic state of the system, this review aims to guide future GeoH 2 research and to evaluate the potential of natural hydrogen as a sustainable clean energy source.

Moradi, Rasoul [Univ. of Southern California, Los ↗
DOE OSTI2018

A new chemical kinetic method of determining RON and MON values for single component and multicomponent mixtures of engine fuels

A new method of using chemical kinetic reaction modeling to predict the Research Octane Number (RON) and Motor Octane Number (MON) of single component fuels and fuel mixtures is described and illustrated via comparisons between computed and experimental values obtained using the well-established ASTM test procedures in a Cooperative Fuels Research (CFR) engine. Comparisons include predictions of RON and MON for a large variety of neat fuels, studies determining the RON and MON of mixtures of primary reference fuels (PRF) and toluene, and studies of RON and MON for mixtures of single-component and multiple-component gasoline surrogate mixtures with ethanol. Advantages in costs, time, and experimental complexity of the kinetic modeling approach compared to the existing engine test procedures are discussed.

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

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the role of reactive nitrogen intermediates in hetero-cyclization during hydrothermal liquefaction of food waste

Hydrothermal liquefaction (HTL) has emerged as a promising strategy for converting abundant, water-rich organic streams into an energy-dense, sustainable, biocrude. Despite major strides in improving oil yields and process intensification, a key underlying challenge persists in relating the composition of the feedstock to the ultimate fate of nitrogen, which is especially important for biocrude quality. To elucidate how nitrogen partitions into the respective gas, aqueous, oil and solid hydrothermal products, food waste was treated under HTL conditions (15 wt% solids, 575 K reaction temperature, and 1 hour reaction time) and elemental, chromatographic and mass spectroscopy analyses were performed on the products to establish dominant product classes and to close mass, carbon, and nitrogen balances. Liquid products (aqueous and biocrude) were found to contain nearly 60% of the nitrogen, with a majority of the biocrude-phase nitrogen in the form of amides. Ab initio thermodynamic and kinetic simulations on model reactants were used to evaluate potential reaction pathways involving reactive nitrogen intermediates. A subsequent kinetic model was evaluated at reaction temperatures ranging 300–1000 K and for a range of feedstock compositions to identify primary reaction pathways. The Maillard reaction is revealed to be the dominant pathway for converting reactive nitrogen intermediates into observed products, especially those in which nitrogen is present as a heterocyclic aromatic. Subsequent product analysis by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) confirms strong agreement among the experimental data, kinetic model and proposed Maillard pathways.

09 BIOMASS FUELS↗
DOE OSTI2018

A detailed reaction mechanism for oxidative coupling of methane over Mn/Na 2 WO 4 /SiO 2 catalyst for non-isothermal conditions

Direct production of ethylene from natural gas or other methane sources has proven to be challenging from catalytic, process design and economic perspectives, with no reports of a commercially viable process to date. Understanding of the underlying non-ideal homogeneous and heterogeneous reactions involved could lead to success, enabling utilization of vast natural gas resources as a replacement for petrochemicals. In this regards, an improved reaction mechanism and computational model, validated experimentally, for oxidative coupling of methane (OCM) over a Mn/Na 2 WO 4 /SiO 2 catalyst is developed, taking into account the non-isothermal behavior exhibited by the gas and surface reactions. The reaction mechanism involves detailed gas-phase and surface elementary steps to describe the OCM chemistry and employing a reduced model for computational efficiency. Our model uses laboratory-scale packed-bed reactor experiments for kinetic mechanism development. The modeling and simulation is particularly focused on a one-dimensional packed-bed reactor. The kinetic model is validated over a wide range of experimental conditions. The model-predicted results and their comparisons with the experimental data are presented wide temperature ranges and methane-to-oxygen feed ratios under adiabatic reactor conditions. The kinetic model is first of its kind in which non-isothermal OCM kinetic behavior is captured using detailed gas-phase and heterogeneous reaction kinetics.

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

Elementary solutions of coupled model equations in the kinetic theory of gases

The method of elementary solutions is employed to solve two coupled integrodifferential equations sufficient for determining temperature-density effects in a linearized BGK model in the kinetic theory of gases. Full-range completeness and orthogonality theorems are proved for the developed normal modes and the infinite-medium Green's function is constructed as an illustration of the full-range formalism. The appropriate homogeneous matrix Riemann problem is discussed, and half-range completeness and orthogonality theorems are proved for a certain subset of the normal modes. The required existence and uniqueness theorems relevant to the H matrix, basic to the half-range analysis, are proved, and an accurate and efficient computational method is discussed. The half-space temperature-slip problem is solved analytically, and a highly accurate value of the temperature-slip coefficient is reported.

Kriese, J. T.↗
NASA NTRS2011

Acceleration of the KINETICS Integrated Dynamical/Chemical Computational Model Using MPI

Understanding the evolution of a planet's atmosphere not only provides a better theoretical understanding of planetary physics and the formation of planets, but also grants useful insight into Earth's own atmosphere. One of the tools used at JPL for the modeling of planetary atmospheres and protostellar disks is KINETICS. KINETICS can simulate years of complex dynamics and chemistry.

planetary physics↗
DOE OSTI2018

The Homogeneous Cooling State as a Verification Test for Kinetic Theory-Based Continuum Models of Gas–Solid Flows

Granular and multiphase (gas–solids) kinetic theory-based models have emerged a leading modeling strategy for the simulation of particle flows. Similar to the Navier–Stokes equations of single-phase flow, although substantially more complex, kinetic theory-based continuum models are typically solved with computational fluid dynamic (CFD) codes. Under the assumptions of the so-called homogeneous cooling state (HCS), the governing equations simplify to an analytical solution describing the “cooling” of fluctuating particle velocity, or granular temperature. The HCS is used here to verify the implementation of a recent multiphase kinetic theory-based model in the open source mfix code. In this work, results from the partial verification test show that the available implicit (backward) Euler time integration scheme converges to the analytical solution with the expected first-order rate. A second-order accurate backward differentiation formula (BDF) is also implemented and observed to converge at a rate consistent with its formal accuracy.

Computational fluid dynamics↗
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗