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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2018

A Drift-Kinetic Analytical Model for SOL Plasma Dynamics at Arbitrary Collisionality

A drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite-Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously-derived fluid models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

The influence of hydrothermal temperature and time toward crystallinity of zeolite X supported on glass wool for CO{sub 2} adsorption

In this study, the influence of hydrothermal temperature and time at zeolite X supported on glasswool were investigated. The results of characterization using XRD showed that a single phase zeolite X with highest crystallinity was obtained when hydrothermal temperature and time at 100°C during 24 hours (ZXF100-24H). The CO{sub 2} adsorption capacity of ZXF100-24H has reached up to 10.15 wt. %. Kinetics of CO{sub 2} adsorption onto zeolite X supported on glasswool was investigated using pseudo-first-order, pseudo-second-order and intra-particle diffusion kinetic models. After evaluating three kinetic models for CO{sub 2} adsorption at adsorption temperatures of 30°C, 40°C and 50°C, it was found that intra-particle diffusion kinetic model provided the best fitting for the adsorption data. Furthermore, the thermodynamic parameters of CO{sub 2} adsorption were obtained as follows, Gibbs free energy change (ΔG°) are −0.409 kJ/mol at 30°C, −0.274 kJ/mol at 40°C and −0.138 kJ/mol at 50 °C, whereas the enthalpy change (ΔH°) is −4.53 kJ/mol and the entropy change (ΔS°) is −0.0135 kJ/(mol K).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

A drift-kinetic analytical model for scrape-off layer plasma dynamics at arbitrary collisionality

We report a drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite–Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously derived fluid models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Isotherm, Kinetic, Process Modeling, and Techno-Economic Analysis of a Diamine-Appended Metal-Organic Framework for CO 2 Capture Using Fixed Bed Contactors

In this report, diamine-appended metal-organic frameworks exhibiting step-shaped CO 2 adsorption are exceptional candidates for energy-efficient carbon capture. However, there are few studies examining their performance in real-world capture scenarios, in part due to the challenge inherent in modeling their CO 2 uptake behavior. Here, we develop a dual-site Sips model to fit experimental CO 2 adsorption data for dmpn-Mg 2 (dobpdc) (dmpn = 2,2-dimethyl-1,3-diaminopropane; dobpdc 4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) and develop a linear driving force model for the adsorption kinetics based on available experimental data. These models are used to develop a dynamic, fixed bed, nonisothermal contactor model using shaped particles of the material, which is validated with experimental breakthrough data. We also examine the effects of the high heat of adsorption of the material on CO 2 uptake performance and find that heat removal is essential to maximize capture performance. We finally investigate "basic"(no bed cooling during adsorption) and "modified"(bed cooling during adsorption) temperature swing adsorption (TSA) processes using dmpn-Mg 2 (dobpdc), and their process economics are compared to a state-of-the-art monoethanolamine (MEA) capture system with and without heat recovery. In the absence of heat recovery, the adsorbent systems are more costly than established technology. However, with 85% heat recovery, both adsorbent-based TSA processes are projected to cost less than the MEA system. This work highlights that thermal management is vital for implementation of dmpn-Mg 2 (dobpdc) as a viable CO 2 capture technology. Investigation of other contactor technologies that can provide unique ways to manage system heat represent promising future areas of study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2015

Energy and momentum preserving Coulomb collision model for kinetic Monte Carlo simulations of plasma steady states in toroidal fusion devices

A kinetic Monte Carlo model suited for self-consistent transport studies is proposed and tested. The Monte Carlo collision operator is based on a widely used model of Coulomb scattering by a drifting Maxwellian and a new algorithm enforcing the momentum and energy conservation laws. The difference to other approaches consists in a specific procedure of calculating the background Maxwellian parameters, which does not require ensemble averaging and, therefore, allows for the use of single-particle algorithms. This possibility is useful in transport balance (steady state) problems with a phenomenological diffusive ansatz for the turbulent transport, because it allows a direct use of variance reduction methods well suited for single particle algorithms. In addition, a method for the self-consistent calculation of the electric field is discussed. Results of testing of the new collision operator using a set of 1D examples, and preliminary results of 2D modelling in realistic tokamak geometry, are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

A moment-based kinetic theory model for polydisperse gas–particle flows

Starting from a generalized population balance equation and the Boltzmann–Enskog collision model for hard spheres, a kinetic theory model for polydisperse gas–particle flows is introduced. Here, polydispersity results from spherical particles with the same material density but different diameters. The particle size distribution (PSD) of the particles is allowed to evolve in space and time due to physical processes such as mixing. In order to treat a continuous PSD, the particle-phase model is formulated in terms of the moments of the PSD, and velocity moments conditioned on the particle size. Velocity moments up to second order are included, resulting in transport equations for the mass, momentum and granular temperature, all conditioned on the particle size. In the numerical implementation, the PSD is represented using quadrature-based moment methods (QBMM). With QBMM, a continuous PSD can be treated using a relatively small number of moments as compared to class or sectional methods. In this work, a realizable numerical algorithm for solving the moment system in a finite-volume code is proposed, valid for dilute systems wherein frictional forces are negligible. The ability of the proposed model to describe polydisperse gas–particle systems is demonstrated using cluster-induced turbulence and riser flow.

42 ENGINEERING↗

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↗

DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Thermodynamics and kinetics of precipitation and austenite reversion during aging of Ti-free grade 300 maraging steel manufactured by laser powder bed fusion (LPBF)

Austenite reversion and intermetallic precipitation are two key phase transformations in maraging steels which determine the strength and ductility. Usually, an optimum aging temperature/time is determined using trial and error approach with inputs from thermodynamic databases and a-posteriori microscopy/hardness evaluations. In this study, a thermo-kinetic model is developed to predict the precipitation and austenite reversion kinetics during various post-fabrication heat treatments in a Ti-free variant of grade 300 maraging steel manufactured using laser powder bed fusion. The thermo-kinetic model was based on the principles of Scheil solidification and the classical nucleation theory. The predictions from the thermo-kinetic model were verified using suitable microscopy techniques, and the results from the thermo-kinetic model agree well with the experimental measurements. Overall, results highlight the fact that an Integrated Computational Materials Engineering approach can be used to predict the precipitation and austenite reversion kinetics, and be used as an input for an efficient design of heat treatments for achieving optimum mechanical properties.

36 MATERIALS SCIENCE↗

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 OSTI2017

Subsequent R + O 2 Chemistry of Intermediates Formed in Low-Temperature R + O 2 Reactions: Potential Importance in Modeling Autoignition Behavior

Comprehensive chemical kinetics models used in the simulation of hydrocarbon and biofuel oxidation rely on accurate prescription of the underlying reaction mechanisms and rate parameters of associated elementary reactions. For practical transportation fuels, such models contain thousands of elementary reactions, which collectively define chain-initiation, -propagation, -branching, and -inhibition pathways. In the low-temperature regime, below approximately 1000 K where R + O 2 reactions dominate, primary oxidation intermediates including cyclic ethers, carbonyls, and conjugate alkenes are formed in abundance via unimolecular decomposition of either chemically activated or thermalized radicals, specifically organic peroxy (ROO) or hydroperoxyalkyl species (QOOH). Experimental results from multiplexed photoionization mass spectrometry (MPIMS) experiments are detailed herein for several intermediates, derived initially from R + O 2 reactions of hydrocarbons and biofuels, and show that intermediate species formed in the initial steps of oxidation undergo similar reactions to those of the parent molecule, including through QOOH-mediated pathways. Products from QOOH decomposition via chain-inhibition and chain-propagation pathways, namely conjugate alkenes, carbonyls, and cyclic ethers, are detected directly. Despite such rich chemistry involving QOOH radicals, most comprehensive chemical kinetics models neglect the complete description of primary oxidation intermediates, and rather consider a restricted number of reaction pathways. It is suggested that exclusion of the details of the oxidation of these intermediate products may affect the interpretation of combustion simulations using such models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2024

Forward and inverse predictive transient models of TREAT using surrogate reactivity models

In this work, we present two novel approaches for predicting the power evolution and control rod height of the Transient Reactor Test Facility (TREAT) to support experiment modeling; specifically, transient analysis of the NASA-sponsored Sirius series of experiments. These approaches utilize steady-state Monte Carlo model, point kinetics model, and surrogate models to predict power evolution and control rod axial position during transient experiments. Both approaches were tested and validated against several Sirius experiments that were performed in the TREAT facility at different power levels. The validation test results show very good agreement with the experimental data, and the models were able to accurately predict the power evolution and the axial control rod position with an average error within 3.0%. Finally, this indicates that these approaches will help the reactor engineering team of the TREAT facility in preparing and predicting the power and temperature of the experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

maranasgroup/yeastsMFA

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

Dinh, Hoang↗

DOE OSTI2023

maranasgroup/k-sacce306

Dinh, Hoang↗

DOE OSTI2016

Eigenvalue Solver for Fluid and Kinetic Plasma Models in Arbitrary Magnetic Topology

ArbiTER (Arbitrary Topology Equation Reader) is a new code for solving linear eigenvalue problems arising from a broad range of physics and geometry models. The primary application area envisioned is boundary plasma physics in magnetic confinement devices; however ArbiTER should be applicable to other science and engineering fields as well. The code permits a variable numbers of dimensions, making possible application to both fluid and kinetic models. Additionally, the use of specialized equation and topology parsers permits a high degree of flexibility in specifying the physics and geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗