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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Investigation of Iso-octane Ignition and Validation of a Multizone Modeling Method in an Ignition Quality Tester

An ignition quality tester was used to characterize the autoignition delay times of iso-octane. The experimental data were characterized between temperatures of 653 and 996 K, pressures of 1.0 and 1.5 MPa, and global equivalence ratios of 0.7 and 1.05. A clear negative temperature coefficient behavior was seen at both pressures in the experimental data. These data were used to characterize the effectiveness of three modeling methods: a single-zone homogeneous batch reactor, a multizone engine model, and a three-dimensional computational fluid dynamics (CFD) model. A detailed 874 species iso-octane ignition mechanism (Mehl, M.; Curran, H. J.; Pitz, W. J.; Westbrook, C. K.Chemical kinetic modeling of component mixtures relevant to gasoline. Proceedings of the European Combustion Meeting; Vienna, Austria, April 14-17, 2009) was reduced to 89 species for use in these models, and the predictions of the reduced mechanism were consistent with ignition delay times predicted by the detailed chemical mechanism across a broad range of temperatures, pressures, and equivalence ratios. The CFD model was also run without chemistry to characterize the extent of mixing of fuel and air in the chamber. The calculations predicted that the main part of the combustion chamber was fairly well-mixed at longer times (> ~30 ms), suggesting that the simpler models might be applicable in this quasi-homogeneous region. The multizone predictions, where the combustion chamber was divided into 20 zones of temperature and equivalence ratio, were quite close to the coupled CFD-kinetics results, but the calculation time was ~11 times faster than the coupled CFD-kinetics model. Although the coupled CFD-kinetics model captured the observed negative temperature coefficient behavior and pressure dependence, discrepancies remain between the predictions and the observed ignition time delays, suggesting improvements are still needed in the kinetic mechanism and/or the CFD model. This approach suggests a combined modeling approach, wherein the CFD calculations (without chemistry) can be used to examine the sensitivity of various model inputs to in-cylinder temperature and equivalence ratios. In conclusion, these values can be used as inputs to the multizone model to examine the impact on ignition delay. Additionally, the speed of the multizone model also makes it feasible to quickly test more detailed kinetic mechanisms for comparison to experimental data and sensitivity analysis.

33 ADVANCED PROPULSION SYSTEMS↗
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 OSTI2014

Comparison of two-fluid and gyrokinetic models for kinetic Alfvén waves in solar and space plasmas

An analytical comparative study of a two-fluid and a gyrokinetic model of kinetic Alfvén waves (KAWs) is presented for various solar and space plasma environments. Based on the linear KAW dispersion relation for gyrokinetics (Howes et al. 2006), the wave group velocity and electromagnetic polarizations are obtained analytically. Then the gyrokinetic wave properties are compared with those of the two-fluid model. The results show that both models agree well with each other not only in the long wavelength regime (>> the ion gyroradius ρ {sub i}) for all cases considered, but also in wavelengths ∼ρ {sub i} and <<ρ {sub i} (still much larger than the electron gyroscale) for a moderate or low (≲ 1) and a high (>>1) ion/electron temperature ratio T {sub 0i}/T {sub 0e}, respectively. However, the fluid model calculations deviate strongly from the gyrokinetic model at scales <ρ {sub i} for a relatively low T {sub 0i}/T {sub 0e} due to the electron gyroradius effect. Meanwhile, the plasma β {sub i} can make the gyrokinetic dispersion relation of KAWs become complex and sometimes have an oscillation-like structure. With the inherent simplicity of the fluid theory, these results may improve our understanding of the applicability of the two-fluid model, and may have important implications for computer simulation studies of KAWs in the solar and space plasma surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2013

A primitive kinetic-fluid model for quasi-parallel propagating magnetohydrodynamic waves

The extension and limitation of the existing one-dimensional kinetic-fluid model (Vlasov-MHD (magnetohydrodynamic) model), which has been used to analyze parametric instabilities of parallel propagating Alfvén waves, are discussed. The inconsistency among the given velocity distribution functions in the past studies is resolved through the systematic derivation of the multi-dimensional Vlasov-MHD model. The linear dispersion analysis of the present model indicates that the collisionless damping of the slow modes is adequately evaluated in low beta plasmas, although the deviation between the present model and the full-Vlasov theory increases with increasing plasma beta and increasing propagation angle. This is because the transit-time damping is not correctly evaluated in the present model. It is also shown that the ponderomotive density fluctuations associated with the envelope-modulated quasi-parallel propagating Alfvén waves derived from the present model is not consistent with those derived from the other models such as the Landau-fluid model, except for low beta plasmas. The result indicates the present model would be useful to understand the linear and nonlinear development of the Alfvénic turbulence in the inner heliosphere, whose condition is relatively low beta, while the existing model and the present model are insufficient to discuss the parametric instabilities of Alfvén waves in high beta plasmas and the obliquely propagating waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Status of Chronic Oxidation Studies of Graphite

Graphite will undergo extremely slow, but continuous oxidation by traces of moisture that will be present, albeit at very low levels, in the helium coolant of HTGR. This chronic oxidation may cause degradation of mechanical strength and thermal properties of graphite components if a porous oxidation layer penetrates deep enough in the bulk of graphite components during the lifetime of the reactor. The current research on graphite chronic oxidation is motivated by the acute need to understand the behavior of each graphite grade during prolonged exposure to high temperature chemical attack by moisture. The goal is to provide the elements needed to develop predictive models for long-time oxidation behavior of graphite components in the cooling helium of HTGR. The tasks derived from this goal are: (1) Oxidation rate measurements in order to determine and validate a comprehensive kinetic model suitable for prediction of intrinsic oxidation rates as a function of temperature and oxidant gas composition; (2) Characterization of effective diffusivity of water vapor in the graphite pore system in order to account for the in-pore transport of moisture; and (3) Development and validation of a predictive model for the penetration depth of the oxidized layer, in order to assess the risk of oxidation caused damage of particular graphite grades after prolonged exposure to the environment of helium coolant in HTGR. The most important and most time consuming of these tasks is the measurement of oxidation rates in accelerated oxidation tests (but still under kinetic control) and the development of a reliable kinetic model. This report summarizes the status of chronic oxidation studies on graphite, and then focuses on model development activities, progress of kinetic measurements, validation of results, and improvement of the kinetic models. Analysis of current and past results obtained with three grades of showed that the classical Langmuir-Hinshelwood model cannot reproduce all data collected so far. Starting from here we propose a modification of the LH model to include temperature activation of graphite surface as a Boltzmann activation function. The enhanced Boltzmann-Langmuir-Hinshelwood model (BLH) was tested successfully on three grades of graphite. The model is a robust, comprehensive mathematical function that allows better fitting of experimental results spanning a wide range of temperature and partial pressures of water vapor and hydrogen. However, the model did not fit satisfactorily the data extracted from the old report on graphite H-451 oxidation by water.

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

Mesoscale Modeling of Kinetic Phase Behaviors in Mg-B-H (Subcontract Report)

Storage of hydrogen on board vehicles is one of the critical enabling technologies for creating hydrogenfueled transportation systems that can reduce oil dependency and mitigate the long-term effects of fossil fuels on climate change. Stakeholders in developing hydrogen infrastructure are currently focused on highpressure storage at 350 bar and 700 bar, in part because no viable solid-phase storage material has emerged. Nevertheless, solid-state materials, including high-density hydrides, remain of interest because of their unique potential to meet all DOE targets and deliver hydrogen at lower pressures and higher on-board densities. A successful solution would significantly reduce costs and ensure the economic viability of a U.S. hydrogen infrastructure. The Mg(BH 4 ) 2 -MgB 2 system represents a highly promising solution because of its reasonable reaction enthalpy, high intrinsic capacity, and demonstrated reversibility, yet suffers from poor reaction kinetics. This subcontract aims to deliver a phase-field model for the kinetics of the evolution of the relevant phases within the Mg-B-H system during hydrogenation and dehydrogenation. This model will be used within a broader theory, synthesis, and characterization framework to study the properties of geometry-selected nanoparticles of pristine and doped MgB 2 /Mg(BH 4 ) 2 with two aims: (1) understand the intrinsic limitations in (de)hydrogenation; (2) devise strategies for improving thermodynamics and kinetics through nanostructuring.

08 HYDROGEN↗
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↗
DOE OSTI2013

An Experimental and Modeling Study of NO x -CO Formation in High Hydrogen Content (HHC) Fuels Combustion in Gas Turbine Applications

Species concentration measurements specifically those associated with NO x can act as important validation targets for developing kinetic models to predict NO x emissions under syngas as well as natural gas combustion accurately. In this collaborative research effort that included both experimental measurements, model development and simulations a comprehensive kinetic model and a multiphysics computational fluid dynamics platform has been developed and validated against the experimental data available in the literature as well as those acquired under this project. The experimental data provide the necessary NO x and speciation data for conditions relevant to gas turbine operations but are not readily available in the literature. The comprehensive chemical kinetic model consists of CO/H 2 /NO x oxidation with the full implementation of NO x evolution pathways, including thermal, prompt, N 2 O and NNH paths. The experiments conducted included NO x perturbed oxidation of natural gas at elevated pressure in laminar flow reactor and syngas/air combustion in a McKenna Burner – Flow Tube setup. A wide range of equivalence ratio, operating pressure as well as H 2 /CO ratio (for syngas only) was investigated. Temperature and NO x concentrations were measured in the flame and post-combustion zone. Experiments were also conducted for seeded syngas where trace hydrocarbon was introduced. The proposed model has been extensively tested. Predictions from the model are compared against multiple experimental datasets over a wide range of venues and operating conditions. The experimental venues include shock tube, plug flow reactor, and stirred reactor experiments that cover pressures from 1 to 100 bar and equivalence ratios from 0.5 to 1.5. In general, the overall model predictions are in good agreement with global combustion targets, such as ignition delay time, as well as with more detailed measurements from flow reactors and stirred reactors. Simulations are conducted for a wide range of reacting mixtures (H 2 /O 2 /N 2 , CO/H 2 /O 2 and CO/H 2 O/O 2 /N 2 ) with initial NO and NO 2 perturbations to consider exhaust gas recirculation (EGR) conditions.

08 HYDROGEN↗
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↗
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 OSTI2014

Finite-density effects in the Fredrickson-Andersen and Kob-Andersen kinetically-constrained models

We calculate the corrections to the thermodynamic limit of the critical density for jamming in the Kob-Andersen and Fredrickson-Andersen kinetically-constrained models, and find them to be finite-density corrections, and not finite-size corrections. We do this by introducing a new numerical algorithm, which requires negligible computer memory since contrary to alternative approaches, it generates at each point only the necessary data. The algorithm starts from a single unfrozen site and at each step randomly generates the neighbors of the unfrozen region and checks whether they are frozen or not. Our results correspond to systems of size greater than 10{sup 7} × 10{sup 7}, much larger than any simulated before, and are consistent with the rigorous bounds on the asymptotic corrections. We also find that the average number of sites that seed a critical droplet is greater than 1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
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 OSTI2015

Development of Kinetic Mechanisms for Next-Generation Fuels and CFD Simulation of Advanced Combustion Engines

Predictive chemical kinetic models are needed to represent next-generation fuel components and their mixtures with conventional gasoline and diesel fuels. These kinetic models will allow the prediction of the effect of alternative fuel blends in CFD simulations of advanced spark-ignition and compression-ignition engines. Enabled by kinetic models, CFD simulations can be used to optimize fuel formulations for advanced combustion engines so that maximum engine efficiency, fossil fuel displacement goals, and low pollutant emission goals can be achieved.

30 DIRECT ENERGY CONVERSION↗
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↗