SEARCH · Engineering Papers

Results for “combustion 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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Probing the Effects of NO x and SO x Impurities on Oxy-Fuel Combustion in Supercritical CO 2 : Shock Tube Experiments and Chemical Kinetic Modeling

The direct-fired supercritical carbon dioxide cycles are one of the most promising power generation methods in terms of their efficiency and environmental friendliness. Two important challenges in implementing these cycles are the high pressure (300 bar) and high CO 2 dilution (>80%) in the combustor. The design and development of supercritical oxy-combustors for natural gas require accurate reaction kinetic models to predict the combustion outcomes. The presence of a small amount of impurities in natural gas and other feed streams to oxy-combustors makes these predictions even more complex. During oxy-combustion, trace amounts of nitrogen present in the oxidizer is converted to NO x and gets into the combustion chamber along with the recirculated CO 2 . Similarly, natural gas can contain a trace amount of ammonia and sulfurous impurities that get converted to NO x and SO x and get back into the combustion chamber with recirculated CO 2 . In this work, a reaction model is developed for predicting the effect of impurities such as NO x and SO x on supercritical methane combustion. The base mechanism used in this work is GRI Mech 3.0. H 2 S combustion chemistry is obtained from Bongartz et al. while NOx chemistry is from Konnov. The reaction model is then optimized for a pressure range of 30–300 bar using high-pressure shock tube data from the literature. It is then validated with data obtained from the literature for methane combustion, H 2 S oxidation, and NO x effects on ignition delay. The effect of impurities on CH 4 combustion up to 16 atm is validated using NO x -doped methane studies obtained from the literature. In order to validate the model for high-pressure conditions, experiments are conducted at the UCF shock tube facility using natural gas identical mixtures with N 2 O as an impurity at ~100 bar. Current results show that there is a significant change in ignition delay with the presence of impurities. A comparison is made with experimental data using the developed model and predictions are found to be in good agreement. The model developed was used to study the effect of impurities on CO formation from sCO 2 combustors. It was found that NO x helps in reducing CO formation while the presence of H 2 S results in the formation of more CO. The reaction mechanism developed herein can also be used as a base mechanism to develop reduced mechanisms for use in CFD simulations.

Energy & Fuels↗

DOE OSTI2022

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

An Assessment of Kinetic Models for Ammonia Flame Extinction

The growing interest in using ammonia as a carbon-neutral replacement fuel has prompted numerous recent research efforts towards understanding fundamental combustion characteristics and chemical kinetics of ammonia oxidation. Despite these recent efforts, there remains a large deficiency in extinction limit data for ammonia flames, particularly under elevated pressure and temperature conditions relevant to gas turbines. Even the literature reaction mechanisms developed specifically for ammonia have been untested under gas turbine relevant conditions, and existing kinetic models have shown large variations in laminar flame speed, extinction strain rate, and speciation predictions for NH3–H2 mixtures at ambient conditions. Recent work by Thomas and co-workers, for example, compared different ammonia kinetic models in the literature for predicting their extinction strain rate measurements for non-premixed NH3-H2 counterflow flames, and again have highlighted the discrepancies in model predictions. Here, our work has undertaken to further highlight important reaction pathways for ammonia oxidation that may be responsible for these discrepancies for extinction strain rate predictions, as well as elucidate sensitive reaction steps with large uncertainties that require further attention. This work points to a need for an improved understanding of ammonia combustion chemistry as well as new datasets at elevated pressure and temperature conditions to better constrain rate expressions for key reaction steps.

ammonia, flame, combustion, kinetic, kinetics↗

NASA NTRS1996

Propellant Chemistry for CFD Applications

Current concepts for reusable launch vehicle design have created renewed interest in the use of RP-1 fuels for high pressure and tri-propellant propulsion systems. Such designs require the use of an analytical technology that accurately accounts for the effects of real fluid properties, combustion of large hydrocarbon fuel modules, and the possibility of soot formation. These effects are inadequately treated in current computational fluid dynamic (CFD) codes used for propulsion system analyses. The objective of this investigation is to provide an accurate analytical description of hydrocarbon combustion thermodynamics and kinetics that is sufficiently computationally efficient to be a practical design tool when used with CFD codes such as the FDNS code. A rigorous description of real fluid properties for RP-1 and its combustion products will be derived from the literature and from experiments conducted in this investigation. Upon the establishment of such a description, the fluid description will be simplified by using the minimum of empiricism necessary to maintain accurate combustion analyses and including such empirical models into an appropriate CFD code. An additional benefit of this approach is that the real fluid properties analysis simplifies the introduction of the effects of droplet sprays into the combustion model. Typical species compositions of RP-1 have been identified, surrogate fuels have been established for analyses, and combustion and sooting reaction kinetics models have been developed. Methods for predicting the necessary real fluid properties have been developed and essential experiments have been designed. Verification studies are in progress, and preliminary results from these studies will be presented. The approach has been determined to be feasible, and upon its completion the required methodology for accurate performance and heat transfer CFD analyses for high pressure, tri-propellant propulsion systems will be available.

Farmer, R. C.↗

NASA NTRS2008

Gaseous Surrogate Hydrocarbons for a Hifire Scramjet that Mimic Opposed Jet Extinction Limits for Cracked JP Fuels

This paper describes, first, the top-down methodology used to define simple gaseous surrogate hydrocarbon (HC) fuel mixtures for a hypersonic scramjet combustion subtask of the HiFIRE program. It then presents new and updated Opposed Jet Burner (OJB) extinction-limit Flame Strength (FS) data obtained from laminar non-premixed HC vs. air counterflow diffusion flames at 1-atm, which follow from earlier investigations. FS represents a strain-induced extinction limit based on cross-section-average air jet velocity, U(sub air), that sustains combustion of a counter jet of gaseous fuel just before extinction. FS uniquely characterizes a kinetically limited fuel combustion rate. More generally, Applied Stress Rates (ASRs) at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or t) can directly be compared with extinction limits determined numerically using either a 1-D or (preferably) a 2-D Navier Stokes simulation with detailed transport and finite rate chemistry. The FS results help to characterize and define three candidate surrogate HC fuel mixtures that exhibit a common FS 70% greater than for vaporized JP-7 fuel. These include a binary fuel mixture of 64% ethylene + 36% methane, which is our primary recommendation. It is intended to mimic the critical flameholding limit of a thermally- or catalytically-cracked JP-7 like fuel in HiFIRE scramjet combustion tests. Our supporting experimental results include: (1) An idealized kinetically-limited ASR reactivity scale, which represents maximum strength non-premixed flames for several gaseous and vaporized liquid HCs; (2) FS characterizations of Colket and Spadaccini s suggested ternary surrogate, of 60% ethylene + 30% methane + 10% n-heptane, which matches the ignition delay of a typical cracked JP fuel; (3) Data showing how our recommended binary surrogate, of 64% ethylene + 36% methane, has an identical FS; (4) Data that characterize an alternate surrogate of 44% ethylene + 56% ethane with identical FS and nearly equal molecular weights; this could be useful when systematically varying the fuel composition. However, the mixture liquefies at much lower pressure, which limits on-board storage of gaseous fuel; (5) Dynamic Flame Weakening results that show how oscillations in OJB input flow (and composition) can weaken (extinguish) surrogate flames up to 200 Hz, but the weakening is 2.5x smaller compared to pure methane; and finally, (6) FS limits at 1-atm that compare with three published 1-D numerical OJB extinction results using four chemical kinetic models. The methane kinetics generally agree closely at 1-atm, whereas, the various ethylene models predict extinction limits that average ~ 45% high, which represents a significant problem for numerical simulation of surrogate-based flameholding in a scramjet cavity. Finally, we continue advocating the FS approach as more direct and fundamental for assessing idealized scramjet flameholding potentials than measurements of "unstrained" premixed laminar burning velocity or blowout in a Perfectly Stirred Reactor.

Pellett, Gerald L.↗

DOE OSTI2025

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗

DOE OSTI2020

Probing the low-temperature chemistry of methyl hexanoate: Insights from oxygenate intermediates

Understanding the combustion of methyl esters is crucial to elucidate kinetic pathways and predict combustion parameters, soot yields, and fuel performance of biodiesel, however most kinetic studies of methyl esters have focused on smaller, surrogate model esters. Methyl hexanoate is a larger methyl ester approaching the chain length of methyl esters found in biodiesel and has not received as much research attention as other smaller esters. The purpose of this work is to present the first atmospheric pressure combustion data of methyl hexanoate, CH 3 CH 2 CH 2 CH 2 CH 2 COOCH 3 . Mixtures of 2% methyl hexanoate in O 2 and N 2 are studied using a plug flow reactor at atmospheric pressure, wall temperatures from 573 to 973 K, residence times from roughly 1-2 s., and fuel equivalence ratios of 1, 1.5, and 2. Exhaust gases are analyzed by a gas chromatograph-mass spectrometer system and species mole fractions are presented. Here, the literature model shows satisfactory agreement with the experimental species profiles and improvements for future mechanistic studies are suggested. In particular, this work proposes new unimolecular decomposition pathways of methyl hexanoate to form methanol or methyl acetate. Furthermore, the experiment detected three unsaturated esters that are direct products of the low temperature oxidation chemistry and it provides more insight into branching ratios for the formation of methyl hexanoate radicals and for the decomposition of hydroperoxyalkyl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Impact of gasoline composition on the effects of nitric oxide on autoignition and knock in a DISI engine

Modern spark-ignition engines use exhaust gas recirculation (EGR) to dilute the charge and suppress knock, enabling the use of higher compression ratios and/or more optimum combustion phasing for higher efficiency. The effectiveness of EGR is affected by the composition of the fuel and its chemical-kinetic interactions with combustion products. Among those, nitric oxide (NO) has been shown to strongly affect autoignition reactivity. However, the impact of fuel composition of the effect of NO on reactivity is not well-understood. Here, in this study, engine experiments were conducted to assess the impact of NO seeded to the intake on knock-limited operation of two gasoline fuels (high cycloalkane content, or HCA, and high olefin content, or HO). Results showed that compositionally-different fuels responded differently to NO. HCA, which was less knock-limited than HO for NO < 200 ppm, became more knock-limited for NO > 200 ppm. Moreover, it was found that differences in knock between fuels were caused by differences in autoignition chemistry and not in the sequential autoignition process of the end gas that occurs due to thermal stratification. Chemical kinetic simulations were performed to better understand the experimental results. For HCA, intermediate-temperature heat release had a greater impact on autoignition reactivity than low-temperature heat release, while the opposite was observed for HO. For both fuels, NO enhances the magnitude of low-temperature heat release via NO + HO 2 → NO 2 + OH. The effect of NO on reactivity was stronger for HCA because OH produced from NO helped to overcome the OH quenching effect of cyclopentane, a main species in HCA. In contrast, HO had relatively strong inherent low-temperature chemistry arising from iso-octane, which reduced the impact of NO on reactivity. For the range of NO mole fractions tested in this study, in-cylinder NO increased fuel’s knock propensity, especially for fuels with mild low-temperature chemistry.

Autoignition↗

NASA NTRS1997

Heterogeneous Combustion of Porous Graphite Particles in Microgravity

Recent theoretical investigations on graphite particle combustion have employed several levels of heterogeneous reaction models, ranging from global to elementary models, to describe the oxidation of carbon to gaseous products. Unlike the counterpart homogeneous reaction models, these heterogeneous reaction models are not well developed because of the difficulties associated with decoupling the physical characteristics of the solid (e.g. surface area taking part in combustion) from the chemical kinetic data. This is certainly true for porous graphite particle combustion, where heterogeneous and homogeneous reactions occur within the pores and play an important role in the overall oxidation process. As a result, there are considerable uncertainties of physical phenomena predicted using different heterogeneous kinetic models available in the literature. A good example, discussed later in this paper, is the predicted critical particle size below which the mass burning rate becomes exponentially small. The main goal of this study is to understand the basic mechanism controlling such rapid changes in burning rates, by developing a model where physical contributions are decoupled from chemical rate constants in a consistent manner. Another important goal of the proposed study is to develop a truly intrinsic, detailed heterogeneous reaction model for porous graphite combustion at high-temperatures, and to derive a systematically reduced heterogeneous reaction model in terms of the elementary reaction rate constants of the detailed model. The validation of chemical kinetic models describing the heterogeneous and homogeneous combustion in and around a spherically symmetric porous graphite particle can be considerably simplified by experimental measurements obtained under microgravity conditions. A vital component of this study is to conduct such supporting experiments on particle burning rate and surface temperature using NASA microgravity facilities, in close coordination with the theoretical effort. The basic understanding obtained and models developed as part of this project will be useful for optimal design of coal combustion devices. These models can also be extended to investigate the role of heterogeneous chemistry on pollutant formation pathways in combustion devices. The theoretical approach developed here, with pore diffusion effects decoupled from the chemical effects, can also be extended to understand the heterogeneous combustion of other porous fuels, for example, combustion of magnesium in a CO2 environment for propulsion in the Martian atmosphere.

Chelliah, Harsha K.↗

DOE OSTI2021

Design and Experimental Demonstration of an Additive Manufactured Smart Oxy-Methane Burner for High Pressure and Supercritical Combustion

Pressurized oxy combustion-based systems can improve efficiency by recovering latent heat of the steam in the Flue Gas and achieving 90% CO2 capture. In addition, the novel Directly Heated Supercritical Carbon Dioxide (DH-SCO2) power cycles can achieve high thermal efficiencies and provide nearly full carbon capture. Additionally, due to the reduction of flue gas at higher pressure, smaller system size and capital cost reductions are also possible. Recent thermodynamic analysis of the DH-SCO2 cycle performed by the UTEP research team shows that combustion conditions in the vicinity of 300 bar pressure and 1000-1400 K temperature allow for relatively high system efficiencies while operating within the limit of available combustor materials. However, the realization of a directly heated supercritical power cycle requires combustion systems to operate in supercritical conditions and at temperature far below the blowout limit of conventional flames (above 1500 K). The thermodynamic properties along with the combustion properties and kinetics are unexplored at such conditions. Additionally, the interaction of the supercritical environment with energy components is also unknown. High-pressure combustion tests are performed using a smart burner at some intermediate pressure ranges (<20 bar) to help minimize these knowledge gaps. The knowledge obtained from the high-pressure test will assist in understanding the combustion chamber pressurization mechanism, ignition and flame behavior at the elevated pressure. The obtained data will act as a systematic first step in testing at higher pressures of 100 and 300 bar pressures.The primary purpose of this Dissertation is to demonstrate the operability of a low NOx smart burner for oxy-methane combustion at high pressure (< 20 bar) and scalable up to supercritical conditions. A shear co-axial smart burner is designed with a real-time temperature monitoring capability to understand the burner face interaction at high-pressure conditions. In addition, the burner has four independent injection ports to allow the independent injection of fuel and dilution gases in the combustor. The maximum operating capability of the burner is 575 kWth. The burner was fabricated using a Laser Powder Bed Fusion process with Nickel Alloy 718. A powder removal technology was invented comprising ultrasonic vibration, liquid nitrogen exposure and media blasting to remove powders from internal channels of the burner. The burner operability tests are performed in the high-pressure combustor and the swirl combustor. The high-pressure combustor was used to investigate the burner operability and thermal soak back at different pressurized conditions. The experimental tests in high-pressure combustor up to 275 kWth input resulted in 16.5 bar chamber pressure and 198°C thermal soaks back to the burner. The burner was capable of providing the required thermal input within a 3% deviation range. In addition, soot formation occurred at high-pressure tests. The swirl combustor was used to observe the flame stability of the burner. Flame lift-off was observed for jet velocities above 450 m/s. Additionally, lift-off decreased for low co-flow velocities. CO2 dilution experiments showed increased flame instabilities for all conditions above 50% dilution ratios. At lower thermal inputs, partial flame blow-offs occurred for dilution ratios above 50%. All conditions significantly reduced the flame temperature and increased the flame lift-off height. Finally, a 2nd generation AM smart burner was designed using the knowledge from the 1st generation burner experiments. The 2nd generation burner incorporated two sets of swirlers with 0.9 swirl no. A cooling system was also designed for long-duration tests at higher pressures. The thermal input and division of the burner's power are kept the same as the 1st generation burner. The burner is to be fabricated using nickel-alloy 718 for high-pressure handling capability. The design can sustain at high-pressure conditions up to 100 bar.

Islam, Md Nawshad Arslan↗

DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1992

An analysis of combustion studies in shock expansion tunnels and reflected shock tunnels

The effect of initial nonequilibrium dissociated air constituents on the combustion of hydrogen in high-speed flows for a simulated Mach 17 flight condition was investigated by analyzing the results of comparative combustion experiments performed in a reflected shock tunnel test gas and in a shock expansion tunnel test gas. The results were analyzed and interpreted with a one-dimensional quasi-three-stream combustor code that includes finite rate combustion chemistry. The results of this study indicate that the combustion process is kinetically controlled in the experiments in both tunnels and the presence of the nonequilibrium partially dissociated oxygen in the reflected shock tunnel enhances the combustion. Methods of compensating for the effect of dissociated oxygen are discussed.

Jachimowski, Casimir J.↗

NASA NTRS2003

A Design Methodology for Rapid Implementation of Active Control Systems Across Lean Direct Injection Combustor Platforms

The VACCG team is comprised of engineers at Virginia Tech who specialize in the subject areas of combustion physics, chemical kinetics, dynamics and controls, and signal processing. Currently, the team's work on this NRA research grant is designed to determine key factors that influence combustion control performance through a blend of theoretical and experimental investigations targeting design and demonstration of active control for three different combustors. To validiate the accuracy of conclusions about control effectiveness, a sequence of experimental verifications on increasingly complex lean, direct injection combustors is underway. During the work period January 1, 2002 through October 15, 2002, work has focused on two different laboratory-scale combustors that allow access for a wide variety of measurements. As the grant work proceeds, one key goal will be to obtain certain knowledge about a particular combustor process using a minimum of sophisticated measurements, due to the practical limitations of measurements on full-scale combustors. In the second year, results obtained in the first year will be validated on test combustors to be identified in the first quarter of that year. In the third year, it is proposed to validate the results at more realistic pressure and power levels by utilizing the facilities at the Glenn Research Center.

Baumann, William T.↗

NASA NTRS1973

Experimental and theoretical studies of NO/x/ formation in a jet-stirred combustor.

As an approach to testing kinetic models which couple NO formation to combustion reactions, a jet-stirred combustor was used to study NO(x) formation under kinetically controlled combustion conditions. The experimental data were compared with kinetic calculations for a well-stirred reactor. The agreement between theory and experiment was quite good for the combustion of hydrogen/air and carbon monoxide/air, for which detailed mechanisms were used. A quasi-global mechanism for the combustion of propane/air resulted in substantial underprediction of NO(x) formation.

Engleman, V. S.↗

NASA NTRS2012

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗

DOE OSTI2021

CFD modeling of pre-spark heat release in a boosted direct-injection spark-ignition engine

Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.

42 ENGINEERING↗

DOE OSTI2024

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1992

Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method - Combustion of methane

A method of systematic optimization, solution mapping, as applied to a large-scale dynamic model is presented. The basis of the technique is parameterization of model responses in terms of model parameters by simple algebraic expressions. These expressions are obtained by computer experiments arranged in a factorial design. The developed parameterized responses are then used in a joint multiparameter multidata-set optimization. A brief review of the mathematical background of the technique is given. The concept of active parameters is discussed. The technique is applied to determine an optimum set of parameters for a methane combustion mechanism. Five independent responses - comprising ignition delay times, pre-ignition methyl radical concentration profiles, and laminar premixed flame velocities - were optimized with respect to thirteen reaction rate parameters. The numerical predictions of the optimized model are compared to those computed with several recent literature mechanisms. The utility of the solution mapping technique in situations where the optimum is not unique is also demonstrated.

Frenklach, Michael↗