Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “premixed flames”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of non-thermal termolecular reactions on detonation development in hydrogen (H 2 )/methane (CH 4 ) - air mixtures

The binary fuel blend of H 2 /CH 4 is one of the most promising hydrogen-enriched hydrocarbon fuels in spark-ignition (SI) engines. Yet, the undesirable phenomenon of super-knock, which can severely and in-stantaneously damage an SI engine, limits its widespread adoption. Moreover, there is still a lack of con-sensus on the precise mechanism by which this phenomenon occurs i.e. via flame acceleration or spon-taneous ignition, despite numerous previous investigations. At the same time, recent studies [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078 -1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522] have demonstrated a high probability of occurrence of non-thermal reactions in premixed flames of such H 2 /CH 4 fuel blends with air due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 . The present study focuses on the evolution of an initial deflagration front to a detonation wave in H 2 /CH 4 - air mixtures under SI engine relevant conditions through fully resolved, constant volume 1D simulations with and without non-thermal reactivity. Non-thermal reactions were included in the macroscopic kinetics model as chemically termolecular reactions facilitated by the H + CH 3 and H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Further, the nonthermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions. Therefore, an additional set of simulations were performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515-522]. Compared to the baseline case, the onset of spontaneous ignition in the end-gas region was observed to be delayed in the presence of non-thermal termolecular reactions. Concurrently, the developing detonation was observed to be significantly stronger. In contrast, applying corrections to the recombi-nation/association rate constants resulted in a completely different behavior. Specifically, detonation was observed to occur due to self acceleration of the primary flame in the absence of spontaneous ignition in the end-gas region. Sensitivity analysis was performed to quantify the effects of non-thermal reactions on the duration of heat release rate and thereby the mechanism of detonation formation. In addition, chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions re-sponsible for the observed results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Heat and mass transfer in combustion - Fundamental concepts and analytical techniques

Fundamental combustion phenomena and the associated flame structures in laminar gaseous flows are discussed on physical bases within the framework of the three nondimensional parameters of interest to heat and mass transfer in chemically-reacting flows, namely the Damkoehler number, the Lewis number, and the Arrhenius number which is the ratio of the reaction activation energy to the characteristic thermal energy. The model problems selected for illustration are droplet combustion, boundary layer combustion, and the propagation, flammability, and stability of premixed flames. Fundamental concepts discussed include the flame structures for large activation energy reactions, S-curve interpretation of the ignition and extinctin states, reaction-induced local-similarity and non-similarity in boundary layer flows, the origin and removal of the cold boundary difficulty in modeling flame propagation, and effects of flame stretch and preferential diffusion on flame extinction and stability. Analytical techniques introduced include the Shvab-Zeldovich formulation, the local Shvab-Zeldovich formulation, flame-sheet approximation and the associated jump formulation, and large activation energy matched asymptotic analysis. Potentially promising research areas are suggested.

Law, C. K.↗
DOE OSTI2022

The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations

Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.

42 ENGINEERING↗
DOE OSTI2021

MACHINE LEARNING BASED CHEMICAL EXPLOSIVE MODE ANALYSIS (ML-CEMA)

The software consists of a Machine Learning based Chemical Explosive Mode Analysis (ML-CEMA) tool for advanced computational flame diagnostics. CEMA, originally based on eigen-analysis of the local thermochemical system, is capable of identifying reaction fronts and limit phenomena such as auto-ignition and extinction in practical combustion systems, such as internal combustion engines and gas turbine combustors. However, the original CEMA is computationally expensive for large reaction mechanisms that are typically needed to describe fuel chemistry of practical large-hydrocarbon fuels. This novel ML-CEMA tool employs a ML technique to accelerate the eigen-analysis of the basic CEMA approach by orders of magnitude, thus making it suitable for practical fuels. In ML-CEMA,zero-dimensional (0D) reactors and one-dimensional (1D) premixed flames are first used to generate a large number of data points for neural network based ML training. The trained ML model is then used to perform CEMA prediction. This ML-CEMA tool has been demonstrated in canonical 0D and 1D configurations as well as highly-transient three-dimensional spray flames exhibiting multi-mode turbulent combustion, showing promising results. ML-CEMA, as a standalone tool, can be used for computationally-efficient diagnostics of massive datasets generated from both experiments and simulations. For example, based on spatially resolved measurements of a small set of reactive scalars(such as temperature, hydroxyl radical and formaldehyde), ML-CEMA can effectively identify flame fronts and rare events. ML-CEMA also provides a robust online or offline flame feature detection tool. When used for on-the-fly simulations, ML-CEMA further enables zone-adaptive combustion modeling, in which the predicted eigenvalue is used as a robust mode indicator for judicious assignment of locally-valid combustion models. This ML-CEMA based zone-adaptive model can lead to substantial computational cost savings when used for large-scale simulations of multi-mode combustion systems. Third Party Code Web Page to Download Code Web Page Location of Third Party License

Xu, Chao↗
NASA NTRS1992

Two-dimensional stability of laminar flames

The stability of laminar flames was studied numerically and the dependence of stability on finite rate chemistry with low activation energy and variable thermodynamic and transport properties was addressed. The calculations show that activation energy and details of chemistry play a minor role in altering the linear neutral stability results from asymptotic analysis. Variable specific heat makes a marginal change to the stability; variable transport properties, on the other hand, tend to substantially enhance the stability from a critical wave number of about 0.50 to 0.20. Also, the effects of variable properties tend to nullify the effects of nonunity Lewis number. When the Lewis number of a single species is different from unity, as is true in a hydrogen-air premixed flame, the stability results remain close to that of unity Lewis number.

Mukunda, H. S.↗
NASA NTRS1992

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗
NASA NTRS2004

Suppression of Low Strain Rate Nonpremixed Flames by an Agent

The extinction and structure of non-premixed methane/air flames were investigated in normal gravity and microgravity through the comparison of experiments and calculations using a counterflow configuration. From a fire safety perspective, low strain rate conditions are important for several reasons. In normal gravity, many fires start from small ignition sources where the convective flow and strain rates are weak. Fires in microgravity conditions, such as a manned spacecraft, may also occur in near quiescent conditions where strain rates are very low. When designing a fire suppression system, worst-case conditions should be considered. Most diffusion flames become more robust as the strain rate is decreased. The goal of this project is to investigate the extinction limits of non-premixed flames using various agents and to compare reduced gravity and normal gravity conditions. Experiments at the NASA Glenn Research Center's 2.2-second drop tower were conducted to attain extinction and temperature measurements in low-strain non-premixed flames. Extinction measurements using nitrogen added to the fuel stream were performed for global strain rates from 7/s to 50/s. The results confirmed the "turning point" behavior observed previously by Maruta et al. in a 10 s drop tower. The maximum nitrogen volume fraction in the fuel stream needed to assure extinction for all strain rates was measured to be 0.855+/-0.016, associated with the turning point determined to occur at a strain rate of 15/s. The critical nitrogen volume fraction in the fuel stream needed for extinction of 0-g flames was measured to be higher than that of 1-g flames.

Olson, Sandra L.↗
NASA NTRS1993

Studies of Premixed Laminar and Turbulent Flames at Microgravity

The work of the Principal Investigator (PI) has encompassed four topics related to the experimental and theoretical study of combustion limits in premixed flames at microgravity, as discussed in the following sections. These topics include: (1) radiation effects on premixed gas flames; (2) flame structure and stability at low Lewis number; (3) flame propagation and extinction is cylindrical tubes; and (4) experimental simulation of combustion processes using autocatalytic chemical reactions.

Ronney, Paul D.↗
DOE OSTI2020

Oxidation of methylamine

A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. In this work, the rate coefficients for the reactions of CH 3 NH 2 + O 2 → CH 2 NH 2 / CH 3 NH + HO 2 , CH 3 NH 2 + H → CH 3 + NH 3 , CH 3 NH → CH 2 NH 2 , and CH 3 NH + O 2 → CH 2 NH + HO 2 were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CH 3 NH 2 and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH 2 group in CH 3 NH 2 compared to CH 4 . More work is desirable to confirm the products of the CH 3 NH + NO reaction and to look into possible pathways to NH 3 in methylamine oxidation.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Model assessment of synthetic jets for turbulent combustion experiments

Abstract Understanding turbulent premixed flames is essential to predict and optimize advanced combustion strategies, but critical capability gaps exist for collecting and validating measurements such as turbulent flame speed. Here, we evaluate synthetic jets as a new, promising turbulence generation device for constant-volume combustion chambers, quantitatively assessing turbulence intensity and spatial uniformity in a hypothetical 4,189-cm 3 vessel for various premixture conditions. Graphic Abstract

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Machine learning assisted modeling of mixing timescale for LES/PDF of high-Karlovitz turbulent premixed combustion

Accurate modeling of mixing in the transported probability density function (PDF) method remains a great challenge, especially for turbulent premixed combustion under extreme conditions such as high Karlovitz number Ka. Recently, a power-law based mixing timescale model was developed for the large-eddy simulations (LES)/PDF modeling of high-Ka number turbulent premixed flames. It is found in this work that the power-law mixing timescale model is highly sensitive to the model parameters. It is thus critically needed to develop accurate calibration of these model parameters. The empirical specification of the model parameters developed in Zhang et. al. is found to be inadequate for accurate modeling of the mixing timescale. Here, machine learning is introduced as an attractive alternative in this work for the specification of the model parameters. A high-Ka number DNS jet flame is used as the training and validation of the machine learning models. The choices of the input parameters are discussed and compared for the machine learning models. The effect of differential molecular diffusion on mixing is examined by including the effect of the Lewis number in the training of the machine learning models. The performance of different machine learning algorithms is compared for the specification of the mixing model parameters. Overall, excellent performance of the machine learning models is observed for assisting the mixing modeling. The feasibility, interpretability, applicability, generality, and portability of using machine learning are discussed in general to provide a perspective on applying data-driven machine learning for turbulent combustion modeling studies.

42 ENGINEERING↗
DOE OSTI2022

Turbulence-induced bias in time-averaged laser absorption tomography of correlated concentration and temperature fields with a first-order correction

The influence of correlated scalar fluctuations on time-averaged laser absorption tomography measurements of temperature and species in a piloted turbulent premixed flame was examined using a coupled spectroscopic and fluid-dynamic analysis. To understand bias associated with turbulence, spatio-temporally resolved temperature and species mole fraction profiles predicted by large eddy simulations (LES) were used to synthetically generate time-resolved line-of-sight absorption measurements at short time scales (microsecond) to reflect the unsteady nature of a canonical jet burner across various transverse measurement planes. Inversion methods were employed on the time-averaged line-of-sight data to produce radially-resolved temperature and mole fraction profiles, analogous to those produced by laser absorption tomography performed on a time-averaged axisymmetric flowfield. It is shown that bias in the measurements compared to true time-averaged scalar fields is a function primarily of temperature dependence in absorptivity and non-zero correlation between temperature and species concentration scalars. Finally, a first-order correction to tomography measurements is proposed to account for the bias based on estimated correlations and the known spectroscopic parameters of the probed absorption transitions.

42 ENGINEERING↗
DOE OSTI2022

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Identity and Chemistry of C 7 H 7 Radicals Observed during Soot Formation

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C 7 H 7 radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H 7 , the results clearly demonstrate that the isomeric composition of C 7 H 7 strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C 7 H 7 isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C 7 H 7 reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C 7 H 7 isomers and predict negligible contributions from other C 7 H 7 isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

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

Cleavage of an aromatic ring and radical migration

The present study undertakes a theoretical evaluation of thermal decomposition of aromatic-ring radicals. Potential energy surfaces and associated reaction rate coefficients were calculated for 1- and 2-naphthalenyl, acetanaphthylenyl, and pyrenyl radicals. Kinetic Monte-Carlo simulations were performed to examine the rates of ring cleavage in two sooting laminar premixed flames of ethylene. The simulations showed that the thermal aromatic-ring cleavage is comparable in rate to oxyradical decomposition in a heavier-sooting flame. The simulation also revealed, unexpectedly, fast internal ring radical migration, comparable in frequency to reaction events of aromatic growth.

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

Optimum Injector Parameters for Thermoacoustic Stability in a Multi-Nozzle Can Combustion System

Abstract High-frequency transverse instabilities are an important concern in can combustor configurations. In these configurations which are typically operated with multiple injectors around a central injector, each injector is subjected to different parts of the acoustic mode shape and thus respond differently for the same instability mode. Recent work by the author has modeled the response of premixed flames to excitation by natural high-frequency transverse modes in a can combustor both in the center and outer nozzles. The stability of these acoustically non-compact flames was assessed using the Rayleigh criterion (Rayleigh Integral denoted as RI) and not the overall unsteady heat release as is the case for compact flames. Several key control parameters were studied, namely — flame angle, swirling strength, nozzle location. For non-axisymmetric modes such as the commonly occurring 1-T mode, both radial and azimuthal offsets of the nozzle location affected stability. The framework was applied to an optimization study to identify the optimal combination of parameters that minimizes RI for the different nozzles in the multi-nozzle system. In this study, a N-around-1 configuration was studied, and the results indicated that the different nozzles needed to be operated at different flame angles and swirl numbers to result in an overall minimum RI. However, the specific response of the different injectors was not considered. The helical mode distribution at each injector varies as we azimuthally go around the combustor’s injector distribution and thus the most amplified mode and the resulting flame response would be different. To minimize RI, it is important to determine the injector configurations that result in a hydrodynamic profile that minimizes the individual RI for each nozzle. The resulting relationship between the injector’s flow and local hydrodynamics can then be used in a hydrodynamics study of an individual injector so that the most optimal injector is chosen depending on its location in the combustor dump plane.

Acharya, Vishal↗
DOE OSTI2020

Development of a Physics-Based Combustion Model for Engine Knock Prediction

The objective of this project is to improve the prediction of engine knock by developing a new combustion modeling framework. Engine knock is a limiting factor to constrain the increase of fuel efficiency for spark ignition (SI) engines in most passenger cars. Efforts to increase fuel efficiency, increasing the compression ratio or downsizing, lead to the increase in the tendency of the knock occurrence. The knock is an undesired ignition of the end-gas, unburned fuel/air mixture ahead of the spark-ignited premixed flame, resulting in rapid in-cylinder pressure rises and engine damages. The combustion modeling framework developed in this project can consider turbulence-chemistry interactions during end-gas ignition, while using a reasonably detailed chemical mechanism developed for ignition and combustion reactions under engine relevant conditions, and the subtle characteristics of spark-ignited flame propagation. It is developed in the context of large eddy simulation (LES), which can capture stochastic in-cylinder processes. The developed model is incorporated into a commercial software for engine simulation, CONVERGE CFD, as a user defined function, and validated. Engine knock and knock-free experiments as well as direct numerical simulation (DNS) of end-gas ignition in homogeneous turbulence are performed to help model development and provide data sets for model validation. With further validation, the developed model is expected to advance the predictive capability for engine knock simulations and thus contribute to improving the fuel efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗