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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Beyond thermodynamic defect models: A kinetic simulation of arsenic activation in CdTe

To overcome the limitations of thermodynamic defect models, we introduce a multilevel kinetic approach to simulate the activation of p doping in CdTe by group V elements using arsenic as a “model” dopant. On the lowest level, we calculate thermodynamic and kinetic parameters of point defects, complexes, and reactions from first principles. On the intermediate level, we use these parameters to calculate the kinetic rates of defect reactions. Finally, we simulate the time evolution of defects and free carriers. Our results show the importance of kinetic factors in defect chemistry models. We reveal the primary arsenic activation pathway to be a fast reaction in which the tellurium atoms get kicked out and replaced on the regular anion sites by interstitial arsenic species. We discover the important role of (As i As Te ) and (As i As i ) complexes that arise during activation anneal to form kinetically stabilized transient states that not only compensate the doping but also can produce deep recombination levels. Lastly, we expect that our modeling approach and the gained insight into the atomic processes behind the doping formation will advance the defect chemistry modeling of electronic properties of materials.

14 SOLAR ENERGY↗
DOE OSTI2022

Studying Cloud and Radiative Impacts Through an Improved Physically Based Representation of Organic Aerosol in a Large-Scale Model (WRF-Chem) (Final Report)

The work performed under this award resulted in the development of a computationally-efficient kinetic model to simulate the chemistry, thermodynamics, and microphysics of organic aerosol. Versions of the kinetic model were used to study a suite of aerosol processes including, multigenerational gas-phase chemistry, phase-state influenced gas/particle partitioning, heterogeneous oxidation, oligomerization reactions, and experimental artifacts. The kinetic model has now been implemented into a regional chemistry-climate model (WRF-Chem) to further study the formation, evolution, properties, and spatiotemporal distribution of organic aerosol at urban, regional, and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of an Optimized Gasoline Surrogate Formulation for PACE Experiments and Simulations

New powertrain solutions are needed to address societal challenges stemming in part from a transportation sector which relies heavily on combustion of conventional hydrocarbon fuels like gasoline. One way to advance new solutions is through predictive simulations. Predictive simulations can potentially limit the extent of experimental validations, reduce time spent in engineering design cycles, unlock new strategies for high efficiency combustion with power density, and minimize tailpipe emissions. However, the simulation tools currently available are either too computationally expensive, inadequate in their accuracy, or a combination of the two. A DOE-funded consortium of national laboratories called the Partnership to Advance Combustion Engines (PACE) seeks to address this gap by rapidly delivering new knowledge and tools which enable market-competitive powertrain solutions for light-duty vehicles. Approaches to modeling combustion in powertrain systems typically incorporate computational fluid dynamics (CFD) simulations with chemical kinetic models. These CFD simulations generally use surrogate fuels featuring a limited number of existing components in a reduced kinetic model to limit computational costs. Such surrogate fuels often sacrifice matching several or more combustion and physical properties of the target gasoline fuel. Validation of surrogate fuels is also often only pursued for a small subset of standard metrics such as research or motor octane number (RON, MON). The purpose of this project is to develop an optimal E10 gasoline surrogate fuel suitable for tasks across PACE to facilitate rapid common analysis and progression toward consortium goals. This project received contributions from tasks led by eight PACE principal investigators. The project funding here reflects all tasks under Pitz, whose tasks also contributed to additional PACE projects. Based on the optimal gasoline surrogate fuel composition, a reduced chemical kinetic model will be shared with PACE researchers, industry, and the broader combustion community.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1992

Multiple lesion track structure model

A multilesion cell kinetic model is derived, and radiation kinetic coefficients are related to the Katz track structure model. The repair-related coefficients are determined from the delayed plating experiments of Yang et al. for the C3H10T1/2 cell system. The model agrees well with the x ray and heavy ion experiments of Yang et al. for the immediate plating, delaying plating, and fractionated exposure protocols employed by Yang. A study is made of the effects of target fragments in energetic proton exposures and of the repair-deficient target-fragment-induced lesions.

Wilson, John W.↗
DOE OSTI2019

Autoignition behavior of iso-olefins

Olefins, and particularly iso-olefins constitute a significant portion of petroleum-derived gasoline, while they also represent a production route for bio-derived transportation-relevant fuels. Although linear structures have been studied for some time, the combustion behavior of iso-olefins is not as well understood. To appropriately mimic gasoline in fundamental tests and in computational models, this lack of understanding must be addressed. Towards this end, rapid compression machine (RCM) experiments were performed for 2-methyl-2-butene (2M2B) and 2-methyl-2-pentene (2M2P) as two representative structures. A detailed chemical kinetic model for these iso-olefins is under development to gain further insight into their autoignition chemistry. In the kinetic model, reaction types and associated rate constants for olefins are less clear than those for n-alkanes. The structural similarity of the two studied iso-olefins allows the development of consistent reaction rate rules to accurately simulate their behavior. The RCM tests cover temperatures from 700 to 1000 K, pressures of 25 and 45 bar, with lean, stoichiometric and rich conditions, and indicate that 2M2P is more reactive than 2M2B due to its longer hydrocarbon backbone. The difference is more extreme at rich conditions. These trends are investigated by the kinetic model where the presence of allylic carbons in the molecule limits the contribution to the low temperature reactivity.

2-methyl-2-butene↗
NASA NTRS2023

The Gamma-Ray Pulsar Phenomenology in View of 3D Kinetic Global Magnetosphere Models

We develop kinetic plasma models of pulsar magnetospheres with magnetic-field-line-dependent plasma injection that reveal the importance of various magnetosphere regions in regulating the gamma-ray emission. We set different particle injection rates for the so-called open, closed, and separatrix zones. Moderate particle injection rates in open and closed zones ensure a global field structure close to the force-free one, while the dissipation occurs mainly in and around the equatorial current sheet. The particles injected in the separatrix zone affect the particle populations that enter the equatorial current sheet region and, therefore, the corresponding accelerating electric fields, particle energies, the spectral cutoff energy, and gamma-ray efficiency. The separatrix zone models reproduce the recently discovered fundamental plane of gamma-ray pulsars consistent with curvature radiation emission, the gamma-ray light-curve shapes, and the radio-lag vs. peak-separation correlation reported in the Fermi second pulsar catalog. The model beaming factors indicate that the pulsar total gamma-ray luminosities listed in the Fermi catalogs are overestimations of the actual ones. We find that the radiation reaction limited regime starts ceasing to govern the high-energy emission for spin-down powers less than 1034 erg/s. Our results also indicate that toward high magnetic inclination angles, the "Y point" around the rotational equator migrates well inside the light cylinder sparking additional peaks in the gamma-ray pulse profiles. We find that an equivalent enhanced particle injection beyond the Y point strengthens these features making the model gamma-ray light curves inconsistent with those observed.

Neutron stars↗
DOE OSTI2020

A numerical study on the sooting tendencies of bio-derived fuels for spark-ignition engines

While more efficient vehicles have been developed over the decades, their performance is limited by the properties of existing fuels. In response, the Co-Optimization of Fuels and Engines Initiative (Co-Optima) under the U.S. Department of Energy (DOE) has developed a rigorous screening approach to identify the most promising biomass-derived blendstocks that are suitable for advanced spark-ignition (SI) engines. A detailed kinetic model has also been developed to predict the combustion properties of the selected blendstocks. This kinetic model has been designed mainly targeting ignition and flame propagation properties, and it has not previously been validated for soot formation. In this work, we numerically predicted the sooting tendency of 20 Co-Optima SI blendstocks using the Co-Optima kinetic model. The sooting tendencies are determined quantitatively using the Yield Sooting Index (YSI) methodology. As shown in Fig. 1, the predicted YSIs show good agreement with measurements, except for 2,5-dimethylfuran. We also quantify the sensitivity of the YSI predictions to aromatic growth reactions and fuel decomposition reactions of the test fuels. It is found that PAH growth reactions have little impact on YSIs for all 20 fuels under investigation. On the other hand, fuel decomposition reactions have a significant influence on the YSI of 2,5-dimethylfuran. Perturbing their reaction rates by a factor of two is shown to lead to 11 % YSI prediction uncertainty.

09 BIOMASS FUELS↗
DOE OSTI2023

Effect of the ..beta..-Hydroxy Group on Ester Reactivity: Combustion Kinetics of Methyl Hexanoate and Methyl 3-Hydroxyhexanoate

One of the major biofuels used today is biodiesel, composed of fatty acid methyl esters. The fats and oils feedstocks used to make biodiesel are in relatively short supply. To meet increasing global demand, engineered microorganisms have been developed to catalyze conversion of sugar to fatty acid esters that contain unique ..beta..-hydroxy-esters. This study investigated the effect of a hydroxyl group on the combustion characteristics of methyl esters by comparing the chemical behavior of methyl hexanoate (MHx) and methyl 3-hydroxyhexanoate (M3OHHx) - used as surrogates for diesel-boiling-range esters with longer fatty acid chains. The oxidation of these esters was studied experimentally in a flow reactor at 0.84 and 10 bar, 600 to 1,100 K, and stoichiometric conditions; and in a constant volume combustion chamber at 5 and 10 bar, 600 to 900 K, and equivalence ratios of 0.3 and 0.6. MHx was more reactive in the constant volume chamber at temperatures below 800 K at 10 bar and equivalence ratio of 0.6. MHx also exhibited a higher indicated cetane number (16.4) than M3OHHx (8.1). We investigated this reactivity trend using an updated MHx kinetic model and M3OHHx model developed as part of this work. The kinetic models predicted that radicals formed from MHx were consumed by low-temperature chemistry reactions, whereas the M3OHHx reactivity was significantly governed by the ..beta..-radical (on the same carbon as the OH group) chemistry which mainly terminated via the less reactive chain propagation pathway to methyl-3-oxohexanoate + HO2. This study extends our understanding of structural effects on ester reactivity, which will allow for accurate surrogate formulation and performance simulations.

AFIDA↗
DOE OSTI2023

Impact of Radiolysis on Iodine Speciation in a Variety of Matrices

This report details the work conducted in fiscal year 2023 with the aim of providing a better understanding of the speciation of iodine in support of domestic 99 Mo production Shine who uses a sulfate matrix in their process. Laboratory experiments to study the behavior of iodine speciation were conducted using spectroscopic methods including UV-vis and Raman as well as electrochemical studies. Radiolysis occurring in solution presents a unique challenge that is largely unaddressed with thermodynamic modeling, therefore a kinetic model was developed to address these issues. The focus of the studies conducted this year used surrogate radiolytic conditions induced by hydrogen peroxide for comparison to radiolytic conditions from direct gamma irradiation and neat acid solutions. For each case a model incorporating the kinetics was used to predict iodine speciation in the defined conditions where the output was directly compared to the experimental conditions. Further modeling was conducted incorporating the process conditions that are specific to Shine with incorporation of radiolysis and gas sparging.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1986

Distribution of nitrogen and defects in SiO(x)N(y)/Si structures formed by the thermal nitridation of SiO2/Si

Previously reported nitrogen distributions in SiO2 films on Si which have been thermally nitrided at 1000 C have been explained by a kinetic model of the nitridation process which rests upon the effects of interfacial strain. A critical test of this kinetic model is the validity of the predictions regarding nitrogen distributions obtained at other nitridation temperatures. In this work, nitrogen distributions determined via X-ray photoelectron spectroscopy are reported for samples nitrided at 800 and 1150 C, and are shown to be consistent with the kinetic model. In addition, the intensity of a fluorine marker is found to correlate with the nitrogen distribution, and is postulated to be related to kinetically generated defects in the dielectric film, consistent with the strain-dependent energy of formation of defects proposed recently to explain electrical data.

Vasquez, R. P.↗
DOE OSTI2020

Revealing the critical role of radical-involved pathways in high temperature cyclopentanone pyrolysis

Cyclopentanone (CPO) is a promising biofuel for spark-ignition engines due to its ring strain and high auto-ignition resistance. Understanding CPO decomposition is crucial for building a high-temperature combustion model. Here we present a comprehensive kinetic model for high-temperature pyrolysis of CPO with verified results from high-pressure shock tube (HPST) measurements. The time- histories of carbon monoxide (CO), ethylene (C 2 H 4 ), and CPO absorbances over the temperature range of 1156-1416 K and pressure range of 8.53-10.06 atm were measured during current experiments. A corresponding detailed kinetic model was generated using the Reaction Mechanism Generator (RMG) with dominant unimolecular/radical-involved decomposition pathways from either previous studies or quantum calculations within the current work. The obtained model containing 821 species and 79,859 reactions exhibited a good agreement with the experimental results. In this study, the absorbance ratio between C 2 H 4 and CO was used as an important factor to validate models and to prove that radical-involved bimolecular pathways were as significant as unimolecular decomposition of CPO. The rate of production (ROP) analysis showed H radicals play a major role in the decomposition, and the whole decomposition process could be divided into three stages based on the H radical concentration. Finally, the insights from present work can be used to generate a better CPO combustion model and help evaluate CPO as an advanced biofuel.

09 BIOMASS FUELS↗
DOE OSTI2020

General H2O Outgassing Model Applicable to Silica-Filled Silicones

Moisture-induced agining and degradation are continuing concerns in long-term storage and use of silicone materials. The silica fillers frequently used to provide mechanical reinforcement in silicone formulations have a high affinity for H2O and act as a source of moisture outgassing later when placed in dry or vacuum applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Approaching PetaVolts per Meter Plasmonics Using Structured Semiconductors

A newly uncovered class of plasmons in the strongly excited limit opens access to unprecedented Petavolts per meter electromagnetic fields with wide-ranging, transformative impact. Unlike conventional plasmons, such plasmons are constituted by non-perturbative, large-amplitude oscillations of the ultradense, delocalized free electron Fermi gas inherent in conductive media. Here structured semiconductors doped to have an appropriate conduction electron density are introduced to tune the properties of the Fermi gas for matched excitation of large-amplitude plasmons using readily available electron beams which enables immediate experimental validation. Specifically, an electrostatic, surface “crunch-in” plasmon is collisionlessly excited by the beam launched inside a tube. Strong excitation due to matching results in relativistic oscillations of the electron gas and unravels unique phenomena. Relativistically induced ballistic electron transport comes about due to relativistic multifold increase in the mean free path and also leads to unconventional heat deposition beyond Ohm’s law. This explains the absence of observed damage or solid-plasma formation in past experiments on conductive samples interacting with electron bunches shorter than 10-13 seconds. Furthermore, relativistic momentum leads to copious tunneling of electron gas across the surface, which then crunches inside the tube. Relativistic effects along with large, localized electron density variations underlying these modes necessitate kinetic approach to theoretical and computational modeling. Kinetic model presented here demonstrates experimental viability of observing tens of gigavolts per meter plasmonic fields excited by matching readily available electron beams to plasmons in semiconductors with 1018cm-3 free electron density, and paves the way for Petavolts per meter plasmonics.

42 ENGINEERING↗
DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

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

Bimolecular Peroxy Radical (RO 2 ) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO 2 ) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH 3 O 2 + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO 2 + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO 2 + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO 2 + OH reaction is a significant sink of RO 2 radicals at the early stage of autoignition, affecting fuel oxidation through RO 2 ↔ QOOH, RO 2 ↔ alkene + HO 2 , or RO 2 + HO 2 ↔ ROOH + O 2 . Additionally, the kinetic stability of the trioxide formed from RO 2 + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). In conclusion, the present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

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

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗
DOE OSTI2018

Nuclear yield reduction in inertial confinement fusion exploding-pusher targets explained by fuel-pusher mixing through hybrid kinetic-fluid modeling

In inertial confinement fusion implosion experiments involving fast laser heating of glass capsules containing gaseous nuclear fuel, some nuclear measurements cannot be rendered by standard hydrodynamics numerical simulations. Here, the calculated values of the nuclear yield are found to be too high by more than a decade in some cases. The first kinetic simulations of the enclosed reacting gas have bridged only half that gap. In this Rapid Communication, using a hybrid ion kinetic-fluid numerical model which takes into account the kinetic interaction between the fuel and the pusher, we reach full agreement with the experimental yield data. We thus unambiguously demonstrate that the kinetic behavior of the system, resulting in both nonthermodynamic equilibrium of the ions comprising the fuel and strong pusher-fuel mixing, is the correct interpretation of the data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗