SEARCH · Engineering Papers

Results for “Kinetic modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Analysis of Nitrogen Shockwave Radiation at Velocity Around 7km/s

In the last decade, the Electric Arc Shock Tube (EAST) facility has been widely used to provide experimental data aiming to validate radiative heating models used for atmospheric entry, mostly at velocities above 9km/s in air. In more recent years, another experimental campaign aimed to look at velocity below 9 km/s and its analysis resulted in a significant improvement of the air kinetic model using new rates found in the literature but discrepancies were still observed, especially for the N2 second positive system. Despite these previous studies, the overall uncertainties related to surface heating for lunar return remain around 50% To reduce these uncertainties, the last experimental campaign at EAST, named test 62, was conducted in pure nitrogen at velocity ranging from 6 to 11 km/s at a freestream pressure of 0.2 Torr to provide more fundamental data on nitrogen. A previous study analyzed the data obtained at a velocity of 10.3 km/s, where atomic radiation is dominant, and provided spatially resolved density measurements of different states of nitrogen. This study will focus on experimental data obtained at shock velocity of 6.88 and 7.15 km/s where molecular nitrogen is expected to be the dominant source of radiation.

Augustin Tibere-Inglesse↗

DOE OSTI2020

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Effects of local order parameter dependent transport coefficient in diblock copolymers under applied electric fields

In this work, we present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics, and, in this work, it is applied to study the effects of electric fields on thin films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of the transport coefficient on the local order parameter is used to relate it with the diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. In addition, transient behavior is found to be affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric field induced alignment of lamellae domains due to the field are found to be in qualitative agreement with experiments. This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin films containing diblock copolymers in the presence of applied electric fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 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↗

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 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 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↗

DOE OSTI2022

Mechanism, thermochemistry, and kinetics of the reversible reactions: C 2 H 3 + H 2 ⇌ C 2 H 4 + H ⇌ C 2 H 5

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C 2 H 3 , C 2 H 4 , and C 2 H 5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol -1 , respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion – where experimental work is very challenging with a consequent dearth of accurate data –here we provide high-level theoretical results for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗