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 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

NASA NTRS1990

Introducing the VRT gas turbine combustor

Melconian, Jerry O.↗

DOE OSTI2018

Experimental and modeling study of the pyrolysis and oxidation of an iso-paraffinic alcohol-to-jet fuel

Single pulse shock tube experiments were performed to determine the pyrolytic and oxidative decomposition products of an alcohol-to-jet fuel (ATJ) composed primarily of about 90% iso-dodecane. Experiments were conducted at 4 bar nominal pressure, 2 ms nominal reaction time, and temperatures ranging from 900 K to 1550 K. For oxidation, the equivalence ratio was 0.27. Gas chromatography (GC) was used to analyze the products of the chemical reactions and to determine their mole fractions. The experimental results were compared against mole fractions computed using a detailed iso-alkane kinetic model, developed for the present study. The model was able to accurately reproduce the experimental results of both ATJ pyrolysis and oxidation. The experimental data and the kinetic modelintroduced in the current study shall aid in progress in understanding the oxidation of the alternative jet fuels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Pulse current auxiliary forging of sintered TiAl alloys based on the investigation of dynamic recrystallization and phase transformation

Highlights: • DRX kinetics model of sintered near γ TiAl alloy was established. • DRX process was accompanied by twinning when compressed at 1125 °C/0.1 s{sup − 1}. • Hot deformation has a promotion to the phase transformation of γ → α. - Abstract: Pulse current auxiliary forging is a new and efficient way to fabricate high performance TiAl based alloys. Dynamic recrystallization (DRX) and phase transformation have a pronounced effect on the subsequent deformability of TiAl based alloys. Hence, the dynamic recrystallization kinetics were determined at the temperature range of 1125–1275 °C and strain rate range of 0.1–0.001 s{sup − 1}. In addition, the DRX mechanism and phase transformation of TiAl based alloys were investigated to model the microstructure evolution in order to provide guidance for the pulse current auxiliary forging process. When the specimens were compressed at 1125 °C/0.1 s{sup − 1}, twinning played an important role in the DRX process. With increasing deformation temperature, the microstructure transformed from near gamma into duplex microstructure. In addition, according to the kinetics of dynamic recrystallization and phase transformation, fully DRXed TiAl based alloys were fabricated through pulse current auxiliary forging. The as-forged TiAl-based alloys show good superplastic performance with an elongation over 410%.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Pyrolysis of Cyclohexane and 1-Hexene at High Temperatures and Pressures—A Photoionization Mass Spectrometry Study

Cycloalkanes are important components of a wide range of fuels. However, there are few experimental data at simultaneously high temperatures and pressures similar to those found in practical systems. Such data are necessary for developing and testing chemical kinetic models. In this study, data relevant to cycloalkane pyrolysis were obtained from high repetition rate shock tube experiments coupled with synchrotron-based photoionization mass spectrometry diagnostics. The pyrolysis of cyclohexane was studied over 1270–1550 K and ~9 bar, while the more reactive primary decomposition product, 1-hexene, was studied at 1160–1470 K and ~5 bar. Insights into the decomposition of the parent molecules, the formation of primary products and the production of aromatic species were gained. Simulations were performed with models for cyclohexane and 1-hexene that were based on literature models. The results indicate that over several hundred microseconds reaction time at high pressures and temperatures the pyrolysis of cyclohexane is largely dominated by reactions initiated by cyclohexyl radicals. Furthermore, good agreement between the simulations and the experiments were observed for cyclohexane and 1-hexene with a modified version of the cyclohexane model. Conversely, the 1-hexene model did not reproduce the experimental observations.

1-alkenes↗

DOE OSTI2019

Validation of Wichner Predictive Model for Chronic Oxidation by Moisture of Nuclear Graphite

This document reports the results of multi-year extensive research on chronic oxidation of nuclear graphite by traces of moisture that might be present in the helium coolant of high-temperature gas-cooled reactors (HTGRs). The research focuses on four grades of nuclear graphite that are considered potential candidates for component fabrication in advanced nuclear reactors. The four grades are (1) PCEA (medium-fine grain, extruded, manufactured by GrafTech Inc., USA); (2) NBG-17 (medium-fine grain, vibrationally molded, manufactured by SGL Carbon, Germany/France); (3) IG-110 (microfine grain, isostatically molded, manufactured by Toyo Tanso, Japan); and (4) 2114 (microfine grain, isostatically pressed, manufactured by Mersen, USA). The experimental part of this research comprised three objectives that were consistently pursued for each graphite material: (1) accelerated oxidation tests in well-controlled conditions (temperature, gas composition) and analysis of chemical kinetic models and of their parameters for each graphite grade; (2) measurements of water vapor transport rates from H 2 O/He mixtures through thin graphite slabs for estimation of the water effective diffusivity specific to each material; (3) controlled oxidation of graphite specimens in H 2 O/He mixtures and characterization of the effects of oxidation temperatures on the density profiles of oxidized layers. These three groups of experimental results were interpreted using the model developed by Wichner et al. for predicting the penetration depth of oxidants (moisture or oxygen) as a function of oxidation conditions. Basically, starting from the kinetic and transport properties of each graphite grade, the ultimate goal was predicting the stable density profile in the oxidized layer and its dependence on the oxidation temperature. The profiles obtained this way were then compared against direct observations of oxidized layers in specimens obtained in well-controlled oxidation conditions. The general agreement between predictions and observations validated the general applicability of the Wichner model. The results obtained in this project were presented at annual Department of Energy–Advanced Reactor Technologies program reviews and at specialist meetings and conferences and were published in Oak Ridge National Laboratory (ORNL) reports and in professional journals. The results published in 2018 had an immediate impact in 2019, when a group from three US universities used ORNL kinetic data as input information in a complex simulation of chronic oxidation by moisture using the same graphite grades during long-time operation in a modular high-temperature gas-cooled reactor. Their conclusion, in line with the conclusions obtained at the end of this project, was that chronic graphite oxidation by moisture during normal operation conditions of an HTGR will be limited to a thin surface layer of graphite blocks at the reactor bottom, where the highest temperatures are. That will not affect the integrity of graphite components and reactor safety, provided the impurity concentration in the helium coolant is rigorously controlled and maintained at the design specifications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1993

Explicit least squares system parameter identification for exact differential input/output models

The equation error for a class of systems modeled by input/output differential operator equations has the potential to be integrated exactly, given the input/output data on a finite time interval, thereby opening up the possibility of using an explicit least squares estimation technique for system parameter identification. The paper delineates the class of models for which this is possible and shows how the explicit least squares cost function can be obtained in a way that obviates dealing with unknown initial and boundary conditions. The approach is illustrated by two examples: a second order chemical kinetics model and a third order system of Lorenz equations.

Pearson, A. E.↗

DOE OSTI2022

Combustion chemistry of alkenes and alkadienes

Alkenes formed, during the refining of crude oil, by cracking the heavier fractions are present in transportation fuels in significant amounts, up to as much as 15–20% in gasoline. Moreover, alkenes are also the major intermediate products of the oxidation of alkanes, which play a significant role in autoignition chemistry. This review has assessed the recent progress in gas-phase detailed kinetic model development for species with C=C double bond, mostly C2–C10 alkenes and 1,3-butadiene. The compiled knowledge on alkene combustion chemistry enabled a better understanding of the influence of the number and the position of the C=C double bond on the chemical kinetics and hence combustion behavior of alkenes in engines. At first, the article gives an extensive overview of fundamental combustion experiments by considering studies of C2–C10 alkenes and 1,3-butadiene in shock tubes, rapid compression machines, laminar flames, and jet-stirred and flow reactors. The value of the data from such experiments is critically discussed. Secondly, this article highlights the important reaction classes involved in alkene oxidation over low-, intermediate- and high-temperature ranges. Combustion chemistry covering C2 to C10 alkenes, with a special emphasis on C2 to C7 isomers is discussed by presenting a large body of experimental and modeling investigations. Detailed chemistry differences between alkene isomers and also between alkenes and alkanes are also addressed. Thirdly, the article presents important reaction pathways for PAH precursor formation in different alkenes. Finally, a summary of the distinguishing features of alkene combustion chemistry and an outlook towards future research in this area are presented. In conclusion, this review is focused on linear and branched chain alkenes, and the chemistry of cyclo-alkenes is not included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗

NASA NTRS1994

Induced nucleation of carbon dust in red giant stars

This study quantitatively tests the proposed model of induced nucleation of carbonaceous grains in carbon-rich red giant stars. Induced nucleation is the process of grain growth initiated by the presence of reactive surfaces provided by seed particles. The numerical study was performed using a deailed chemical kinetic model of carbon deposition, grain coagulation, and homogeneous nucleation of polycyclic aromatic hydrocarbons (PAHs). The model uses a method of moments to keep track of developing grain population in the forming dust shell. We test the efficiency of grain formation for large ranges of dust shell parameters typical for carbon stars. Our model is capable of producing a range of optically thick and thin dust shells in carbon stars. Results are in accord with (IRAS) spectral classes of carbon stars. The resulting composite grains produced are consistent with those recently found in ancient meteorites. This model also provides a realistic explanation for high abundances of (PAHs) in the interstellar medium and some planetary nebulae.

Cadwell, Brian J.↗

NASA NTRS1986

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗

DOE OSTI2019

Cell free biosynthesis of isoprenoids from isopentenol

Cell-free systems are growing in importance for the biosynthesis of complex molecules. These systems combine the precision of traditional chemistry with the versatility of biology in creating superior overall processes. Recently, a new synthetic pathway for the biosynthesis of isoprenoids using the substrate isopentenol, dubbed the isopentenol utilization pathway (IUP), was demonstrated to be a promising alternative to the native 2C-methyl-d-erythritol-4-phosphate (MEP) and mevalonate (MVA) pathways. This simplified pathway, which contains a minimum of four enzymes to produce basic monoterpenes and only depends on ATP and isopentenol as substrates, allows for a highly flexible approach to the commercial synthesis of isoprenoid products. In this work, we use metabolic reconstitution to characterize this new pathway in vitro and demonstrate its use for the cell-free synthesis of mono-, sesquit-, and diterpenoids. Kinetic modeling and sensitivity analysis were also used to identify the most significant parameters for taxadiene productivity, and metabolic control analysis was employed to elucidate protein-level interactions within this pathway, which demonstrated that the IUP enzymatic system is primarily controlled by the concentration and kinetics of choline kinase (CK) and not regulated by any pathway intermediates. We report this is a significant advantage over the natural MEP or MVA pathways as it greatly simplifies future metabolic engineering efforts, both in vitro and in vivo, aiming at improving the kinetics of CK. Finally, we used the insights gathered to demonstrate an in vitro IUP system that can produce 220 mg/L of the diterpene taxadiene, in 9 hr, almost 3-fold faster than any system reported thus far.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

A quantitative and spatially resolved analysis of the performance-bottleneck in high efficiency, planar hybrid perovskite solar cells

Hybrid perovskites represent a potential paradigm shift for the creation of low-cost solar cells. Current power conversion efficiencies (PCEs) exceed 22%. However, despite this, record PCEs are still far from their theoretical Shockley–Queisser limit of 31%. To increase these PCE values, there is a pressing need to understand, quantify and microscopically model charge recombination processes in full working devices. Here, we present a complete microscopic account of charge recombination processes in high efficiency (18–19% PCE) hybrid perovskite (mixed cation and methylammonium lead iodide) solar cells. We employ diffraction-limited optical measurements along with relevant kinetic modeling to establish, for the first time, local photoluminescence quantum yields, trap densities, trapping efficiencies, charge extraction efficiencies, quasi-Fermi-level splitting, and effective PCE estimates. Correlations between these spatially resolved parameters, in turn, allow us to conclude that intrinsic electron traps in the perovskite active layers limit the performance of these state-of-the-art hybrid perovskite solar cells.

14 SOLAR ENERGY↗

DOE OSTI2020

Molecular recognition of a host protein by NS1 of pandemic and seasonal influenza A viruses

The 1918 influenza A virus (IAV) caused the most severe flu pandemic in recorded human history. Nonstructural protein 1 (NS1) is an important virulence factor of the 1918 IAV. NS1 antagonizes host defense mechanisms through interactions with multiple host factors. One pathway by which NS1 increases virulence is through the activation of phosphoinositide 3-kinase (PI3K) by binding to its p85β subunit. Here we present the mechanism underlying the molecular recognition of the p85β subunit by 1918 NS1. Using X-ray crystallography, we determine the structure of 1918 NS1 complexed with p85β of human PI3K. We find that the 1918 NS1 effector domain (1918 NS1 ED ) undergoes a conformational change to bind p85β. Using NMR relaxation dispersion and molecular dynamics simulation, we identify that free 1918 NS1 ED exists in a dynamic equilibrium between p85β-binding–competent and –incompetent conformations in the submillisecond timescale. Moreover, we discover that NS1 ED proteins of 1918 (H1N1) and Udorn (H3N2) strains exhibit drastically different conformational dynamics and binding kinetics to p85β. These results provide evidence of strain-dependent conformational dynamics of NS1. Using kinetic modeling based on the experimental data, we demonstrate that 1918 NS1 ED can result in the faster hijacking of p85β compared to Ud NS1 ED , although the former has a lower affinity to p85β than the latter. Our results suggest that the difference in binding kinetics may impact the competition with cellular antiviral responses for the activation of PI3K. We anticipate that our findings will increase the understanding of the strain-dependent behaviors of influenza NS1 proteins.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗