Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Summary of Simplified Two Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydrogen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (greater than l x 10(exp -20)) moles per cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T(sub 4)). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/Air fuel and for H2/O2. A similar correlation is also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T(sub 4)) as a function of overall fuel/air ratio, pressure and initial temperature (T(sub 3)). High values of the regression coefficient R squared are obtained.

Marek, C. John↗
NASA NTRS2005

Simplified Two-Time Step Method for Calculating Combustion Rates and Nitrogen Oxide Emissions for Hydrogen/Air and Hydorgen/Oxygen

A simplified single rate expression for hydrogen combustion and nitrogen oxide production was developed. Detailed kinetics are predicted for the chemical kinetic times using the complete chemical mechanism over the entire operating space. These times are then correlated to the reactor conditions using an exponential fit. Simple first order reaction expressions are then used to find the conversion in the reactor. The method uses a two-time step kinetic scheme. The first time averaged step is used at the initial times with smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, temperature, and pressure. The second instantaneous step is used at higher water concentrations (> 1 x 10(exp -20) moles/cc) in the mixture which gives the chemical kinetic time as a function of the instantaneous fuel and water mole concentrations, pressure and temperature (T4). The simple correlations are then compared to the turbulent mixing times to determine the limiting properties of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. This time is regressed over the complete initial conditions using the Excel regression routine. Chemical kinetic time equations for H2 and NOx are obtained for H2/air fuel and for the H2/O2. A similar correlation is also developed using data from NASA s Chemical Equilibrium Applications (CEA) code to determine the equilibrium temperature (T4) as a function of overall fuel/air ratio, pressure and initial temperature (T3). High values of the regression coefficient R2 are obtained.

Molnar, Melissa↗
DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗
DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

60 - APPLIED LIFE SCIENCES↗
NASA NTRS1995

An innovative approach to the development of a portable unit for analytical flame characterization in a microgravity environment

The availability of manned laboratory facilities in space offers wonderful opportunities and challenges in microgravity combustion science and technology. In turn, the fundamentals of microgravity combustion science can be studied via spectroscopic characterization of free radicals generated in flames. The laser-induced fluorescence (LIF) technique is a noninvasive method of considerable utility in combustion physics and chemistry suitable for monitoring not only specific species and their kinetics, but it is also important for imaging of flames. This makes LIF one of the most important tools for microgravity combustion science. Flame characterization under microgravity conditions using LIF is expected to be more informative than other methods aimed at searching for effects like pumping phenomenon that can be modeled via ground level experiments. A primary goal of our work consisted in working out an innovative approach to devising an LIF-based analytical unit suitable for in-space flame characterization. It was decided to follow two approaches in tandem: (1) use the existing laboratory (non-portable) equipment and determine the optimal set of parameters for flames that can be used as analytical criteria for flame characterization under microgravity conditions; and (2) use state-of-the-art developments in laser technology and concentrate some effort in devising a layout for the portable analytical equipment. This paper presents an up-to-date summary of the results of our experiments aimed at the creation of the portable device for combustion studies in a microgravity environment, which is based on a portable UV tunable solid-state laser for excitation of free radicals normally present in flames in detectable amounts. A systematic approach has allowed us to make a convenient choice of species under investigation, as well as the proper tunable laser system, and also enabled us to carry out LIF experiments on free radicals using a solid-state laser tunable in the UV.

Dubinskiy, Mark A.↗
DOE OSTI2022

Investigating the kinetic effect of prenol on iso -octane auto-ignition by means of an experimental and modeling study

A detailed experimental and kinetic modeling study was dedicated to understand the reported octane hyperboosting effect of prenol, by means of the measurement of the ignition delay times of its blends with iso -octane, and measurement of the mole fraction profiles of the fuels and intermediates inside the ULille rapid compression machine. Here, these results show that prenol addition leads to a reduction of the first-stage ignition phenomena and negative temperature coefficient behavior, which is only qualitatively captured by the model and is consistent with knock resistance improvement. It is suggested that this behavior is caused by two different factors. The first originates from gas-phase reactivity of prenol, and spans from the formation of unreactive unsaturated species through resonance-stabilized radicals, thereby constituting a competitive pathway for the radical pool generated by iso -octane. The second is of catalytic nature and cannot be captured by means of gas-phase kinetic modeling, but could also play an important role in the behavior of prenol in internal combustion engines.

42 ENGINEERING↗
DOE OSTI2026

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2020

Development of a high temperature pyrolysis mechanism for cyclopentanone, a potential biofuel derived from biomass

Cyclopentanone (CPO) is a promising biofuel for spark-ignition engines due to its ring strain and high auto-ignition resistance. The Department of Energy has identified CPO as a potential biofuel candidate because of its unique chemical properties. Understanding CPO decomposition is crucial to building a high-temperature combustion model. In this study, 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 (C2H4), and CPO absorbances over the temperature range of 1156-1416 K and 8.53-10.06 atm were measured. A corresponding kinetic model was generated using the Reaction Mechanism Generator (RMG) with the parameters of dominant pathways from either previous studies or quantum calculation within current work. The final model containing 821 species and 79,859 reactions exhibited a good agreement with the experimental results. In this study, the C2H4/CO ratio was used as an essential factor to validate the model and to prove that radical-involved bimolecular pathways were as significant as unimolecular decomposition of CPO. The rate of production (ROP) analysis showed that H radicals play major roles in both primary reactions and product pyrolysis. The insights from this work can be used to generate a better CPO combustion model and help evaluate CPO as a biofuel. Such understanding can also be helpful when proposing pyrolysis or combustion models for other biofuels.

cyclopentanone, reaction mechanism, pyrolysis, sho↗
DOE OSTI2025

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Combustion characteristics and detailed simulations of surrogates for a Tier II gasoline certification fuel

An experimental and numerical study of combustion of a gasoline certification fuel (‘indolene’), and four (S4) and five (S5) component surrogates for it, is reported for the configurations of an isolated droplet burning with near spherical symmetry in the standard atmosphere, and a single cylinder engine designed for advanced compression ignition of pre-vaporized fuel. The intent was to compare performance of the surrogate for these different combustion configurations and to assess the broader applicability of the kinetic mechanism and property database for the simulations. A kinetic mechanism comprised of 297 species and 16797 reactions was used in the simulations that included soot chemistry and accounted for unsteady transport, liquid diffusion inside the droplet, radiative heat transfer, and variable properties. The droplet data showed a clear preference for the S5 surrogate in terms of burning rate. The simulations showed generally very good agreement with measured droplet, flame, and soot shell diameters. Measurements of combustion timing, in-cylinder pressure, and mass-averaged gas temperature were also well predicted with a slight preference for the S5 surrogate. Preferential vaporization was not evidenced from the evolution of droplet diameter but was clearly revealed in simulations of the evolution of mixture fractions inside the droplets. The influence of initial droplet diameter (Do) on droplet burning was strong, with S5 burning rates decreasing with increasing Do due to increasing radiation losses from the flame. Flame extinction was predicted for Do =3.0 mm as a radiative loss mechanism but not predicted for smaller Do for the conditions of the simulations.

Guo, Songtao↗
DOE OSTI2022

Combustion characteristics and detailed simulations of surrogates for a Tier II gasoline certification fuel

An experimental and numerical study of combustion of a gasoline certification fuel (‘indolene’), and four (S4) and five (S5) component surrogates for it, is reported for the configurations of an isolated droplet burning with near spherical symmetry in the standard atmosphere, and a single cylinder engine designed for advanced compression ignition of pre-vaporized fuel. The intent was to compare performance of the surrogate for these different combustion configurations and to assess the broader applicability of the kinetic mechanism and property database for the simulations. A kinetic mechanism comprised of 297 species and 16,797 reactions was used in the simulations that included soot formation and evolution, and accounted for unsteady transport, liquid diffusion inside the droplet, radiative heat transfer, and variable properties. The droplet data showed a clear preference for the S5 surrogate in terms of burning rate. The simulations showed generally very good agreement with measured droplet, flame, and soot shell diameters. Measurements of combustion timing, in-cylinder pressure, and mass-averaged gas temperature were also well predicted with a slight preference for the S5 surrogate. Preferential vaporization was not evidenced from the evolution of droplet diameter but was clearly revealed in simulations of the evolution of mixture fractions inside the droplets. Here, the influence of initial droplet diameter (D o ) on droplet burning was strong, with S5 burning rates decreasing with increasing D o due to increasing radiation losses from the flame. Flame extinction was predicted for D o =3.0 mm as a radiative loss mechanism but not predicted for smaller D o for the conditions of the simulations.

42 ENGINEERING↗
NASA NTRS2015

Development of a Premixed Combustion Capability for Scramjet Combustion Experiments

Hypersonic air-breathing engines rely on scramjet combustion processes, which involve high speed, compressible, and highly turbulent flows. The combustion environment and the turbulent flames at the heart of these engines are difficult to simulate and study in the laboratory under well controlled conditions. Typically, wind-tunnel testing is performed that more closely approximates engine testing rather than a careful investigation of the underlying physics that drives the combustion process. The experiments described in this paper, along with companion data sets being developed separately, aim to isolate the chemical kinetic effects from the fuel-air mixing process in a dual-mode scramjet combustion environment. A unique fuel injection approach is taken that produces a nearly uniform fuel-air mixture at the entrance to the combustor. This approach relies on the precombustion shock train upstream of the dual-mode scramjet combustor. A stable ethylene flame anchored on a cavity flameholder with a uniformly mixed combustor inflow has been achieved in these experiments allowing numerous companion studies involving coherent anti-Stokes Raman scattering (CARS), particle image velocimetry (PIV), and planar laser induced fluorescence (PLIF) to be performed.

Rockwell, Robert D.↗
DOE OSTI2021

Multidimensional Numerical Modeling of Combustion Dynamics in a Non-Premixed Rotating Detonation Engine With Adaptive Mesh Refinement

In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

42 ENGINEERING↗
DOE OSTI2020

Effects of Zr substitution on soot combustion over cubic fluorite-structured nanoceria: Soot-ceria contact and interfacial oxygen evolution

Ceria is widely used as a catalyst for soot combustion, but effects of Zr substitution on the reaction mechanism is ambiguous. The present work elucidates effects of Zr substitution on soot combustion over cubic fluorite-structured nanoceria. The nanostructured CeO 2 , Ce 0.92 Zr 0.08 O 2 , and Ce0.84Zr0.16O2 composed of 5–6 nm crystallites display Tm-CO 2 (the temperature at maximum CO 2 yield) at 383, 355, and 375°C under 10 vol.% O 2 /N 2 , respectively. Additionally, the size of agglomerate decreases from 165.5 to 51.9–57.3 nm, which is beneficial for the soot-ceria contact. Moreover, Zr increases the amount of surface oxygen vacancies, generating more active oxygen (O2- and O-) for soot oxidation. Thus, the activities of Ce 0.92 Zr 0.08 O 2 and Ce 0.84 Zr 0.16 O 2 in soot combustion are better than that of CeO2. Although oxygen vacancies promote the migration of lattice O 2- , the enriched surface Zr also inhibits the mobility of lattice O 2- . Therefore, the Tm-CO 2 of Ce 0.84 Zr 0.16 O 2 is higher than that of Ce 0.92 Zr 0.08 O 2 . Based on reaction kinetic study, soot in direct contact with ceria preferentially decomposes with low activation energy, while the oxidation of isolated soot occurs through diffusion with high activation energy. The obtained findings provide new understanding on the soot combustion over nanoceria.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane

The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature.

Harstad, Kenneth G.↗
DOE OSTI2020

Experimental and modeling study of C2–C4 alcohol autoignition at intermediate temperature conditions

C2–C4 alcohols are advantageous blendstocks identified by many research groups, including the U.S. Department of Energy Co-Optima Initiative, towards enabling efficient, boosted Spark-Ignition (SI) engines. Their use in advanced engine applications requires a comprehensive understanding of their intermediate-temperature autoignition behavior. This work reports an experimental and modeling study covering their fundamental autoignition characteristics in a twin-piston rapid compression machine at pressures of 20 and 40 bar, intermediate temperatures from 750 to 980 K, and two fuel loading conditions representative of boosted SI engines. Direct comparison between these alcohols is made, where the order of reactivity is established across different thermodynamic and fuel loading conditions. Changes in preliminary exothermicity (or intermediate-temperature heat release) displayed in single-stage autoignition across different alcohols and conditions are also quantified. This provides insight into fuel-to-fuel differences, and how these could affect advanced combustion concepts such as spark-assisted compression ignition. Kinetic models are used to simulate the experiments, and reasonable agreement is obtained. Furthermore, the sensitivity analysis results demonstrate the importance of accurately capturing the autoignition kinetics, particularly H-abstraction reactions on the parent fuels by OH and HO 2 , and the branching ratio associated with these.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗