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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Thermodynamics and Kinetics of Phase Transformations in Energy Materials (Final Report)

The primary goal of this program is to develop computational tools and methodology for enabling a shift from empirical trial-and-error driven materials discovery towards data and information technology driven design of materials with targeted functional properties. Both thermodynamic and kinetic factors need to be taken into account when designing new materials. Thermodynamics sets the necessary-but-not-sufficient conditions for the suitability of a particular material in a given application, while the kinetics controls the rates of microscopic processes and hence the performance in energy applications (e.g., power in energy storage, conversion efficiencies in solar cells and thermoelectrics, chemical activity in catalysis). The complexity of these phenomena makes computational design difficult, not least due to the lack of suitable theoretical methods and computational tools. To address these needs, we develop systematic, quantitatively accurate first-principles methods to address two important problems: (i) predict and optimize the thermodynamics of materials over a large composition and structure space, and (ii) model kinetic processes such as mass and heat transport in crystalline solids under non-equilibrium conditions. These tasks are addressed by combining high-throughput computation with automatic construction of accurate lattice models using mathematically rigorous methods of information theory, such as compressive sensing and Bayesian inference.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental and modeling study of C2–C4 alcohol autoignition at intermediate temperature conditions

C2–C4 alcohols are advantageous blendstocks identified by many research groups, including the U.S. Department of Energy Co-Optima Initiative, towards enabling efficient, boosted Spark-Ignition (SI) engines. Their use in advanced engine applications requires a comprehensive understanding of their intermediate-temperature autoignition behavior. This work reports an experimental and modeling study covering their fundamental autoignition characteristics in a twin-piston rapid compression machine at pressures of 20 and 40 bar, intermediate temperatures from 750 to 980 K, and two fuel loading conditions representative of boosted SI engines. Direct comparison between these alcohols is made, where the order of reactivity is established across different thermodynamic and fuel loading conditions. Changes in preliminary exothermicity (or intermediate-temperature heat release) displayed in single-stage autoignition across different alcohols and conditions are also quantified. This provides insight into fuel-to-fuel differences, and how these could affect advanced combustion concepts such as spark-assisted compression ignition. Kinetic models are used to simulate the experiments, and reasonable agreement is obtained. Furthermore, the sensitivity analysis results demonstrate the importance of accurately capturing the autoignition kinetics, particularly H-abstraction reactions on the parent fuels by OH and HO 2 , and the branching ratio associated with these.

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

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗
NASA NTRS2016

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗
DOE OSTI2021

Fuel Properties and Kinetics

The Co-Optimization of Fuels and Engines (CoOptima) program is focused on a two-pronged approach to develop new, lower net-carbon fuels and advanced engine technologies in tandem. This effort brings nine national labs and several university and industry partners together to simultaneously identify new fuel candidates and to further engine technology development in the light-, medium-, and heavy-duty vehicle sectors with a focus on fuel economy and emission reduction. The fuel properties and kinetics project is aimed at achieving a better understanding of fuel property kinetics and how they impact engine performance such that more efficient and lower emitting engines can be made a reality. Several promising biofuel candidates were identified through tiered screening approaches and were evaluated utilizing unique bench-scale measurement techniques. Kinetic models were developed to describe autoignition and soot pre-cursor formation mechanisms to identify the most promising candidates that advanced various combustion strategies.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

Here, the growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O 2 ] n , and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

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

Coarse-Grained Modeling Elucidates Differential Metabolism of Saccharomyces cerevisiae under Varied Nutrient Limitations

Microorganisms such as Saccharomyces cerevisiae have a native ability to adapt their metabolism to varying nutrient conditions. Understanding their responses to nutrient limitations is critical for decoding cellular physiology and designing strategies for metabolic engineering. While the influence of carbon availability on yeast metabolism has been extensively studied, the role of nitrogen availability remains relatively underexplored. In this study, we utilized a coarse-grained kinetic model to systematically analyze and compare the effects of carbon and nitrogen limitations on yeast metabolism. Our model successfully revealed the differential metabolic characteristics of S. cerevisiae under carbon- and nitrogen-limited chemostat conditions. It also highlighted the significance of protein activity regulation at varying carbon-to-nitrogen ratios, and elucidated distinct strategies employed to maintain ATP homeostasis. Furthermore, this study provides a computational tool for investigating yeast physiology under nutrient limitations and offers quantitative and mechanistic insights into yeast metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

System of kinetic equations describing large-scale processes in collisionless space plasma

A system of kinetic equations describing relatively slow large-scale processes in collisionless magnetoplasma structures with a spatial resolution of about the characteristic gyroradius is derived. Plasma is assumed to be quasineutral, while the magnetic and electric fields are determined by the instantaneous distributions of the particle and current densities and the stress tensor of all plasma components in the longrange instantaneous interaction approximation. A special version of equations is derived for the case of magnetized electrons described by the Vlasov equation in the drift approximation. The obtained system of equations can be used to develop a global numerical kinetic model of the Earth’s magnetosphere with a spatial resolution of about 100 km, as well as local models of certain regions of the Earth’s magnetosphere with a higher resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Magnetogasdynamic Power Extraction and Flow Conditioning for a Gas Turbine

An extension of the Russian AJAX concept to a turbojet is being explored. This magnetohydrodynamic (MHD) energy bypass engine cycle incorporating conventional gas turbine technology has MHD flow conditioning at the inlet to electromagnetically extract part of the inlet air kinetic energy. The electrical power generated can be used for various on-board vehicle requirements including plasma flow control around the vehicle or it may be used for augmenting the expanding flow in the high speed nozzle by MHD forces to generate more thrust. In order to achieve this interaction, the air needs to be ionized by an external means even up to fairly high flight speeds, and the leading candidates may be classified as electrical discharge devices. The present kinetic modeling calculations suggest that the use of electron beams with characteristics close to the commercially available e-beam systems (electron energy approx. 60 keV, beam current approx. 0.2 mA/sq cm) to sustain ionization in intermediate pressure, low-temperature (P = 0.1 atm, T = 300 K) supersonic air flows allows considerable reduction of the flow kinetic energy (up to 10 to 20 percent in M = 3 flows). The calculations also suggest that this can be achieved at a reasonable electron beam efficiency (eta approx. 5), even if the e-beam window losses are taken into account. At these conditions, the exit NO and O atom concentrations due to e-beam initiated chemical reactions do not exceed 30 ppm. Increasing the beam current up to approx. 2 mA/sq cm, which corresponds to a maximum electrical conductivity of sigma(sub max) approx. 0.8 mho/m at the loading parameter of K = 0.5, would result in a much greater reduction of the flow kinetic energy (up to 30 to 40 percent). The MHD channel efficiency at these conditions would be greatly reduced (to eta approx. 1) due to increased electron recombination losses in the channel. At these conditions, partial energy conversion from kinetic energy to heat would result in a significant total pressure loss (P(sub 0)/P(sub 0i) approx. 0.3). The total pressure loss can be reduced operating at the loading parameter closer to unity, at the expense of the reduced electrical power output. Raising the beam current would also result in the increase of the exit O atom concentrations (up to 600 ppm) and NO (up to 150 ppm).

Adamovich, Igor V.↗
NASA NTRS2010

Numerical Simulation of Vitiation Effects on a Hydrogen-Fueled Dual-Mode Scramjet

The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.

Vyas, Manan A.↗
DOE OSTI2024

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. Here, we develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

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

NMR spin-rotation relaxation and diffusion of methane

The translational diffusion-coefficient and the spin-rotation contribution to the 1 H NMR relaxation rate for methane (CH 4 ) are investigated using MD (molecular dynamics) simulations, over a wide range of densities and temperatures, spanning the liquid, supercritical, and gas phases. Here, the simulated diffusion-coefficients agree well with measurements, without any adjustable parameters in the interpretation of the simulations. A minimization technique is developed to compute the angular velocity for non-rigid spherical molecules, which is used to simulate the autocorrelation function for spin-rotation interactions. With increasing diffusivity, the autocorrelation function shows increasing deviations from the single-exponential decay predicted by the Langevin theory for rigid spheres, and the deviations are quantified using inverse Laplace transforms. The 1 H spin-rotation relaxation rate derived from the autocorrelation function using the “kinetic model” agrees well with measurements in the supercritical/gas phase, while the relaxation rate derived using the “diffusion model” agrees well with measurements in the liquid phase. 1 H spin-rotation relaxation is shown to dominate over the MD-simulated 1 H- 1 H dipole-dipole relaxation at high diffusivity, while the opposite is found at low diffusivity. At high diffusivity, the simulated spin-rotation correlation time agrees with the kinetic collision time for gases, which is used to derive a new expression for 1 H spin-rotation relaxation, without any adjustable parameters.

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

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

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

Understanding fundamental effects of biofuel structure on ignition and physical fuel properties

The development of biochemicals and biofuels with advantaged properties of higher combustion efficiency at lower emission levels can be aided by a priori prediction of critical global characteristics based on molecular structure. Existing predictive techniques are incomplete for comprehensive ignition and physical fuel property prediction. Here this investigation focuses on filling this knowledge gap by developing a priori prediction techniques for predicting important ignition as well as physical fuel properties of promising advantaged biofuels. We demonstrate the utility of this a priori prediction technique on saturated (≤C 5 ) alcohols, both linear and branched in nature by correlating their molecular structure with both important ignition and physical properties of the fuels. This work utilizes a semi-automated method of predicting ignition characteristics via accurate fundamental-based calculations of the fuel radical cascades that govern low-temperature combustion. This numerical methodology is supported by the relevant fuel ignition delay time (IDT) and research octane number (RON) data experimentally measured by the Advanced Fuel Ignition Delay Analyzer (AFIDA). Phyiscal fuel property trends are also analyzed by correlating the molecular structure of the fuel to its physical characteristics. Based on their molecular structure, recommendations are made on an a priori method for selecting biofuels that have advantaged properties for the selected fuel application. Lastly, the utility of this investigation is extended beyond the subset of alcohols considered by analyzing the IDT trends of i-Pentanol and n-Hexanol, revealing the fidelity of the novel kinetic modelling approach and the accuracy of the predicted trends in this study.

09 BIOMASS FUELS↗