Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics”

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

Plasma-based global pathway analysis to understand the chemical kinetics of plasma-assisted combustion and fuel reforming

The Global Pathway Analysis (GPA) algorithm helps analyze the chemical kinetics of complex combustion systems by identifying important global reaction pathways connecting a source species to a sink species through various important intermediate species (i.e., hub species). Here, the present work aims to extend GPA algorithm to plasma-assisted combustion and fuel reforming systems to identify the dominant global pathways in such systems at various conditions. In addition, the present study extends the ability of GPA algorithm to identify reaction cycles involving the excitation of high-concentration species (e.g., O 2 , N 2 , and fuel) to their vibrational and electronic states and the subsequent de excitation to their ground state, based on their significance on the reactivity of plasma-assisted systems in terms of gas heating and radical production. Provisions are made in the GPA algorithm to evaluate the reactivity of identified re action pathways and cycles based on the element-flux transfer (i.e., dominance), heat release, and radical production rate. The newly developed Plasma-based Global Pathway Analysis (PGPA) algorithm is then used to analyze the plasma assisted combustion of ammonia and reforming of methane. The PGPA analyses elucidated the significance of vibrational-translational cycles on the reactivity of NH 3 /air mixtures. Further, analyses on the production of NO ascribed the early reforming of NH 3 to N 2 and H 2 in impeding the production of NO during plasma-assisted NH 3 ignition. Lastly, the enhanced reforming of CH 4 /N 2 mixtures using plasma has been attributed to electron impact dissociation of CH 4 when compared to thermal reforming. In contrast, conventional path-Flux analysis (PFA) was found to require significant manual effort and pre-analysis intuitions from expert knowledge, making it arduous to provide valuable in sights into plasma chemistry. The user-friendly and automated nature of PGPA thus provides a valuable tool for assessing the kinetics of plasma-assisted systems helpful in analyzing and, further, a foundation in reducing plasma-assisted chemistry, without the needs of expert knowledge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Genetic algorithm optimization of a chemical kinetic mechanism for propane at engine relevant conditions

Propane has demonstrated significant potential for reductions in greenhouse gas and pollutant emissions in medium- and heavy-duty engine applications, but further improvements require accurate, compact, and scalable chemical kinetic mechanisms to design the next generation of propane fueled engines, particularly at the boosted operating conditions necessary to meet the power density demand of medium- and heavy-duty applications. In this work, six key chemical reactions were identified in a reduced mechanism with 70 species and 352 reactions through a sensitivity analysis performed at conditions typical of thermodynamic trajectories observed in a high compression ratio, long stroke engine operated on propane from throttled to boosted operating conditions. While the original mechanism was validated against rapid compression machine (RCM) data, it was found to overpredict experimental autoignition tendencies in 2-zone, 0-D SI engine simulations performed in Chemkin Pro. Subsequently, a genetic algorithm approach was used to optimize the six reaction rate parameters within established uncertainty bounds by performing RCM simulations and comparing to two independent sets of literature ignition delay times for propane, thus generating two new kinetic mechanisms. The first optimization achieved a mean absolute percent error (MPE) reduction in 2nd stage ignition delay of 61.4% in seven generations, while the second optimization utilized a newer experimental RCM dataset, and achieved MPE reduction of 56.7% in seven generations, and further marginal improvement to 57.8% reduction in 34 generations. Finally, the two mechanisms were then evaluated again in the 2-zone 0-D SI engine model in Chemkin Pro comparing typical mean and knocking cycle trajectories, and it was found that the second optimized mechanism provided better prediction of knock onset at the representative conditions evaluated in this work, particularly for higher load operating conditions.

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

CHEMNODE: CHEMICAL KINETICS SOLVER APPROACH BASED ON NEURAL ORDINARY DIFFERENTIAL EQUATIONS

This software contains an algorithm to model and predict the time evolution of chemical kinetics in reacting flow simulations using neural ordinary differential equations (NODEs). The approach works by using artificial neural networks to predict the chemical source terms. It achieves this by using forward-mode automatic differentiation and the Levenberg-Marquardt algorithm to adjust the neural network parameters, such that the discrepancies between the actual and predicted species profiles are minimized.

OWOYELE, OPEOLUWA↗
DOE OSTI2021

Burning velocities of R-32/O 2 /N 2 mixtures: Experimental measurements and development of a validated detailed chemical kinetic model

This work entails characterizing the flammability of the refrigerant R-32 (CH 2 F 2 ) by both experimental measurements and modeling. Burning velocities S u were measured using a constant-volume spherical-flame method for R-32/O2/N 2 mixtures with O 2 /N 2 ratios ranging from 21% (synthetic air) to 40%, pressures of (1 to 3) bar, and equivalence ratios $\phi$ of (0.8 to 1.3). Based on a critical assessment of available data, and extended by our own calculations, a detailed chemical kinetic model was developed and key reactions determined using reaction path and sensitivity analyses. Initiation and combustion were identified as distinct kinetic regimes and burning velocities were found to be controlled by two primary reactions: unimolecular decomposition of CH 2 F 2 → CHF + HF and the subsequent reaction, CHF + O 2 → CHFO + O, the latter reaction initiating the radical chain propagating and branching by producing O atoms. Sensitive rate constants in the kinetic model were critically adjusted within their uncertainties and current knowledge bounds to best fit the experimental burning velocities. We found that rate constants in the model could be adjusted to match a given experimental S u for specific conditions (O 2 loading, P, T, $\phi$). This, however, then fixes predicted burning velocities for other all conditions within (3 to 4)% if physically realistic rate parameters are maintained. Thus, the entire set of experimental data is fit, not just to particular conditions. Relative random uncertainties in the experimental Su measurements were (4 to 6)%, but assumptions made for thermal radiation lost by the burned gas in the spherical-flame experiments add an additional systematic uncertainty. Systematic differences between the limiting cases of adiabatic (no thermal radiation lost) and optically-thin (all thermal radiation lost) varied significantly with conditions and ranged from (4 to 30)% at high to low velocities, respectively, translating into uncertainties of (2 to 15)% considering the average of two limiting cases. Comparison of experimental and kinetically modeled Su values suggests that the burned gas tends towards the optically-thin limit at the lowest pressures and fuel loadings and toward the adiabatic limit at the highest pressures and loadings. We tested and found support for this conclusion with a detailed analysis as a function of all the conditions (T, P, % O 2 , $\phi$). This behavior appears to transition from optically-thin to adiabatic as the density of the initial fuel increases, which results in increased CO 2 in the burned gas and thus increased absorption of the thermal radiation (consistent with the Beer-Lambert Law). The validated detailed model based on evaluated kinetics is shown to accurately predict burning velocities for R-32 O 2 /N 2 mixtures over a wide range of conditions and provides a reliable basis for extrapolation to other conditions.

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

ReactionMechanismSimulator.jl: A modern approach to chemical kinetic mechanism simulation and analysis

Abstract We present ReactionMechanismSimulator.jl (RMS), a modern differentiable software for the simulation and analysis of chemical kinetic mechanisms, including multiphase systems. RMS has already been applied to problems in combustion, pyrolysis, polymers, pharmaceuticals, catalysis, and electrocatalysis. RMS is written in Julia, making it easy to develop and allowing it to take advantage of Julia's extensive numerical computing ecosystem. In addition to its extensive library of optimized analytic Jacobians, RMS can generate and use Jacobians computed using automatic differentiation and symbolically generated analytic Jacobians. RMS is demonstrated to be faster than Cantera and Chemkin in several benchmarks. RMS also implements an extensive set of features for analyzing chemical mechanisms, including a library of easy‐to‐call plotting functions, molecular structure resolved flux diagram generation, crash analysis, traditional sensitivity analysis, transitory sensitivity analysis, and an automatic mechanism analysis toolkit. RMS implements efficient adjoint and parallel forward sensitivity analyses. We also demonstrate the ease of adding new features to RMS.

Johnson, Matthew S.↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2024

A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

02 PETROLEUM↗
DOE OSTI2022

Autoignition Delay Time Measurements and Chemical Kinetic Modeling of Hydrogen/Ammonia/Natural Gas Mixtures

Abstract In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels’ characteristics for use in high-pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study CO formation profiles as well as the ignition delay times behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida’s shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles. Experimental results were then compared with the GRI 3.0 mechanism, as well as the proprietary UCF 2022 mechanism utilizing CHEMKIN-Pro software. In general, both models were able to capture the trend in autoignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia-hydrogen mixtures, GRI 3.0 failed to predict ignition delay times whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Finally, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.

Baker, Jessica B.↗
DOE OSTI2022

Experimental Ignition Delay Time Measurements and Chemical Kinetics Modeling of Hydrogen/Ammonia/Natural Gas Fuels

In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels' characteristics for use in high pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study carbon monoxide (CO) formation profiles as well as the ignition delay times (IDTs) behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida's shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles over a temperature range between 1200 K and 1800 K with an average pressure of 2.2 atm. Experimental mixtures included both neat and combination natural gas/hydrogen with ammonia addition, with all mixtures except one having an equivalence ratio of 1. Results were then compared with the GRI 3.0 mechanism, as well as the newly developed UCF 2022 mechanism utilizing chemkin-pro software. In general, both models were able to capture the trend in auto-ignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia–hydrogen mixtures, GRI 3.0 failed to predict ignition delay times, whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Lastly, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Evaluating Chemical Kinetics Predictions for Propane Using 3-D and 0-D Models in a Boosted Spark-Ignited Engine

Propane has been shown to be a promising alternative fuel to reduce emissions while simultaneously achieving high efficiencies in medium- and heavy-duty engines. These high-power density applications require boosted engines which, combined with high compression ratio, can lead to auto-ignition and knock. While three-dimensional (3-D) computational fluid dynamics (CFD) models are often used for resolving the complex fluid flow in engines, these models can become computationally expensive when simulating detailed chemical kinetics. Likewise, zero-dimensional (0-D) models are computationally concise enough for kinetics development, but lack any flow-field information which governs the flame propagation processes in spark ignition (SI) engines. This work presents a comprehensive comparison between 3-D and 0-D closed cycle simulations at knocking conditions in a high compression ratio high stroke-to-bore ratio propane engine. In order to initialize the flow-field for the 3-D closed cycle (intake valve closing, (IVC) to exhaust valve opening, (EVO)) simulation, a motored multi-cycle 3-D model was run using Converge to create a map at IVC, reducing the computational time. The map allowed a non-homogeneous 3-D closed cycle simulation to be satisfactorily validated against experiments, while a homogeneous case using only the turbulence field mapping was also simulated, mimicking 0-D modeling. The 3-D simulations were used to prescribe the initial conditions (e.g., IVC thermodynamics, speciation, burn-rate profile) for a 2-zone 0-D SI engine model in Chemkin Pro for both cases. It was found that 2-zone 0-D modeling underpredicted the knock onset timing, likely due to the lack of thermal stratification in the unburned gas region. Future work will carry multi-zone 0-D modeling to capture the fuel auto-ignition in the unburned region.

Douvry-Rabjeau, Julien [Oakland University, Roches↗
DOE OSTI2026

Impact of Nd ions on the chemical kinetic behavior of radiolytic transients in molten LiCl-KCl-NdCl3 salt mixtures

Pyrochemical reprocessing technologies can be used to recover valuable materials from used nuclear fuel (UNF), such as uranium. However, other critical materials, such as fission product neodymium, are challenging to recover from molten chloride eutectic mixtures (LiCl-KCl) due to the presence of multivalent states (Nd2+/Nd3+), disproportionation reactions, and re-dissolution in the molten chloride salt. Understanding, predicting, and controlling these processes is further complicated by the presence of ionizing radiation fields, for which little is known on the interaction of neodymium ions with the radiolytic transients (es– and Cl2.–) in molten LiCl-KCl salt mixtures. Here, chemical kinetics, from integrated electron pulse irradiation and transient absorption spectroscopy, are presented for the reaction of Nd3+ ions with es– and Cl2.– as a function of temperature (400–600 ?).

36 - MATERIALS SCIENCE↗
DOE OSTI2020

The Sensitivity of Chemical Kinetic Models on Flame Transfer Functions in Acoustic Fluctuation Environments with Varying Equivalence Ratios

Various factors affect the heat release response of flames – thermodynamics, hydrodynamics, acoustics, etc. The coupling of heat release rate and acoustics, in particular, is essential to the design of combustion devices so as to avoid adverse effects during their operation. Modeling flames in acoustic environments requires an accurate understanding of the fundamental physics that occurs when flames interact with a fluctuating pressure field. Over the past couple of decades, researchers have been quantifying the impact of this fluctuating pressure field on the unsteady heat release rate. The flame transfer functions computed have also been shown to depend on the level of detail in the chemical kinetic model. This paper attempts to bridge the gap between these studies to understand how sensitive the flame transfer function is to various kinetic models. The problem is studied by performing multiple one-dimensional simulations using the low-Mach solver NGA. The detailed GRI-Mech 3.0 will be used for comparison. Also, the pressure amplitude and reactant equivalence ratio are varied to understand how the environment affects this sensitivity.

Girdhar, Anant↗
DOE OSTI2021

A chemical kinetics perspective on thermoelectric transport

There is a growing community of chemists and materials scientists exploring thermoelectric properties of organic, polymeric, hybrid, and composite materials. Many of these materials are nondegenerate semiconductors, meaning that the Fermi and transport charge transport energy levels are significantly offset. Others are ionic conductors. While the meaning of the essential parameters electrical conductivity (σ) and Seebeck coefficient (S, the voltage difference per degree of temperature difference) are accessible, the origins of S are not readily apparent to one trained mainly in chemistry. The purpose of this manuscript is to illustrate S as being the result of a steady state difference between concentrations of two designations of charge carriers, those on hotter and those on colder sides of a material sample, analogous to the equilibrium among interconverting chemical species that differ in enthalpy and entropy. The ion-based Seebeck coefficient, also known as the Soret effect, can be explained using principles similar to those applicable to electrons and holes. We hope that this analysis leads to wider understanding of the origins of S through an explanation using the language of chemical kinetics and thermodynamics, and appreciation of ways that thermoelectric efficiency can be enhanced in emerging materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Chemical Kinetics Underlying the Sooting Tendency and Auto-Ignition Characteristics of Linear, Branched, and Cyclic Ether Compounds

Biofuels present opportunities for improving the performance and reducing emissions from internal combustion engines by incorporating oxygenated functional groups to the fuels. Among various oxygenates, ethers have been recognized as promising candidates for an alternative to conventional diesel fuel owing to their higher reactivity and lower sooting tendency. The detailed guidelines for designing ethers, however, have not been fully discussed, even though their combustion characteristics are sensitive to the molecular structure. This study was devoted to exploring the structure-property relationships, particularly focusing on the cetane number and yield sooting index, using five linear, branched, and cyclic ethers: di-amyl-ether, 4-butoxy-heptane, 3,3-dimethyl-oxetane, 2-ethyl-4-methyl-1,3-dioxolane, and 2-isopropyl-4-methyl-1,3-dioxolane. First, we examined the chemical kinetics underlying the sooting tendency of the test fuels. The combustion product distribution was measured from flow reactor experiments at 750-1100 K, F=3, at atmospheric pressure. As a result, it was revealed that the sooting tendency is closely related to the size of hydrocarbon intermediates in the high-temperature regime (>1000 K); that is, larger hydrocarbons lead to more soot precursor formation. The underlying chemistry determining the size of the hydrocarbons from the tested fuels was analyzed using reaction pathway analysis and quantum mechanics calculations, which showed that the branched and cyclic ether structures form abundant C3-C4 compounds. Moreover, the auto-ignition characteristics of the test fuels were studied using the flow reactor at low-temperature (400-700 K) and F=1. We found a clear difference in the combustion-product distribution from high and low reactivity fuels, which was then correlated to the systematic analysis of the key reaction energy barriers with the varying molecular structure.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
DOE OSTI2023

Infrared compatible rapid mixer to probe millisecond chemical kinetics

Fast microfluidic mixers are a valuable tool for studying solution-phase chemical reaction kinetics and molecular processes with spectroscopy. However, microfluidic mixers that are compatible with infrared vibrational spectroscopy have seen only limited development due to the poor infrared transparency of the current microfabrication material. We describe the design, fabrication, and characterization of CaF2-based continuous flow turbulent mixers, which are capable of measuring kinetics in the millisecond time window with infrared spectroscopy, when integrated into an infrared microscope. Kinetics measurements demonstrate the ability to resolve relaxation processes with 1 millisecond time resolution, and straightforward improvements are described that should result in sub-100 µs time-resolution.

Instruments & Instrumentation↗
DOE OSTI2021

A chemical kinetic mechanism for combustion and flame propagation of CH 2 F 2 /O 2 /N 2 mixtures

Here in this work, we compiled and evaluated rate expressions for reactions relevant to the decomposition and combustion of CH 2 F 2 (difluoromethane, refrigerant R-32) in CH 2 F 2 /O 2 /N 2 flames. The recommended values have been used in premixed flame calculations, reported elsewhere, to model experimentally derived burning velocities determined using a constant volume spherical flame method. In this work, we also provide a detailed description of the reaction pathways for decomposition and combustion of CH 2 F 2 . This work is part of a larger effort at NIST to characterize and predict the flammability of new refrigerant working fluids and their blends for consideration as replacements of current refrigerants with high global warming potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗