A computational analysis of strained laminar flame propagation in a stratified CH4/H2/air mixture

Propagation of a H₂-added strained laminar CH₄/air flame in a rich-to-lean stratified mixture is numerically studied. The back-support effect, which is known to enhance the consumption speed of a flame propagating into a leaner mixture compared to that into a homogeneous mixture, is evaluated. A new method is devised to characterize unsteady reactant-to-reactant counterflow flames under transiently decreasing equivalence ratio, in order to elucidate the influence of flow strain on the back-support effect. In contrast to the conventional reactant-to-product configurations, the current configuration is more relevant to unsteady stratified flames back-supported by their own combustion products. Moreover, since H₂ distribution downstream of the flame is known to play a crucial role in back-supported CH₄/air flames, the influence of H₂ addition in the upstream mixture is examined. The results suggest that a larger strain rate leads to a larger equivalence ratio gradient at the reaction zone through increased flow divergence, which amplifies the back-support. Meanwhile, since H₂ addition in the upstream mixture does not affect the downstream H₂ content, the relative increase in the consumption speed, i.e. the back-support, is suppressed with larger H₂ addition. Especially, when the upstream H₂ content decreases with the equivalence ratio, the H₂ preferentially diffuses toward the unburned gas, which mitigates H₂ accumulation in the preheat zone and further weakens the back-support.     

Effect of Soret diffusion on lean hydrogen/air flames at normal and elevated pressure and temperature

The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H₂ while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H₂ is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H₂ on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H₂ caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.     

Effects of equivalence ratio variation on lean,stratified methane–air laminar counterflow flames

The effects of equivalence ratio variations on flame structure and propagation have been studied computationally. Equivalence ratio stratification is a key technology for advanced low emission combustors. Laminar counterflow simulations of lean methane–air combustion have been presented which show the effect of strain variations on flames stabilized in an equivalence ratio gradient, and the response of flames propagating into a mixture with a time-varying equivalence ratio. ‘Back supported’ lean flames, whose products are closer to stoichiometry than their reactants, display increased propagation velocities and reduced thickness compared with flames where the reactants are richer than the products. The radical concentrations in the vicinity of the flame are modified by the effect of an equivalence ratio gradient on the temperature profile and thermal dissociation. Analysis of steady flames stabilized in an equivalence ratio gradient demonstrates that the radical flux through the flame, and the modified radical concentrations in the reaction zone, contribute to the modified propagation speed and thickness of stratified flames. The modified concentrations of radical species in stratified flames mean that, in general, the reaction rate is not accurately parametrized by progress variable and equivalence ratio alone. A definition of stratified flame propagation based upon the displacement speed of a mixture fraction dependent progress variable was seen to be suitable for stratified combustion. The response times of the reaction, diffusion, and cross-dissipation components which contribute to this displacement speed have been used to explain flame response to stratification and unsteady fluid dynamic strain.     

Numerical investigation of flame–vortex interactions in laminar cross-flow non-premixed flames in the presence of bluff bodies

Flame stabilisation in a combustor having vortices generated by flame holding devices constitutes an interesting fundamental problem. The presence of vortices in many practical combustors ranging from industrial burners to high speed propulsion systems induces vortex–flame interactions and complex stabilisation conditions. The scenario becomes more complex if the flame sustains after separating itself from the flame holder. In a recent study [P.K. Shijin, S.S. Sundaram, V. Raghavan, and V. Babu, Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body, Combust. Theory Model. 18, 2014, pp. 692–710], the authors reported details of the regimes of flame stabilisation of non-premixed laminar flames established in a cross-flow combustor in the presence of a square cylinder. In that, the separated flame has been shown to be three dimensional and highly unsteady. Such separated flames are investigated further in the present study. Flame–vortex interactions in separated methane–air cross flow flames established behind three bluff bodies, namely a square cylinder, an isosceles triangular cylinder and a half V-gutter, have been analysed in detail. The mixing process in the reactive flow has been explained using streamlines of species velocities of CH₄ and O₂. The time histories of z-vorticity, net heat release rate and temperature are analysed to reveal the close relationship between z-vorticity and net heat release rate spectra. Two distinct fluctuating layers are visible in the proper orthogonal decomposition and discrete Fourier transform of OH mass fraction data. The upper fluctuating layer observed in the OH field correlates well with that of temperature. A detailed investigation of the characteristics of OH transport has also been carried out to show the interactions between factors affecting fluid dynamics and chemical kinetics that cause multiple fluctuating layers in the OH.     

Experimental and numerical study of product gas characteristics of ammonia/air premixed laminar flames stabilized in a stagnation flow

Ammonia has widely attracted interest as a potential candidate not only as a hydrogen energy carrier but also as a carbon free fuel for internal combustion engines, such as gas turbines. Because ammonia contains a nitrogen atom in its molecule, nitrogen oxides (NO_x) and other pollutants may be formed when it burns. Therefore, understanding the fundamental product gas characteristics of ammonia/air laminar flames is important for the design of ammonia-fueled combustors to meet stringent emission regulations. In this study, the product gas characteristics of ammonia/air premixed laminar flames for various equivalence ratios were experimentally and numerically investigated up to elevated pressure conditions. In the experiments, a stagnation flame configuration was employed because an ammonia flame can be stabilized by using such a configuration without a pilot flame. The experimental results showed that the maximum NO mole fraction was about 3,500 ppmv, at an equivalence ratio of 0.9 at 0.1 MPa. The NO mole fraction decreased as the equivalence ratio increased. In addition, the maximum value of the NO mole fraction decreased with an increase in mixture pressure. Furthermore, it was experimentally clarified that the simultaneous reduction of NO and unburnt ammonia can be achieved at an equivalence ratio of about 1.06, which is the target equivalence ratio for emission control in rich-lean two-stage ammonia combustors. Comparison of experimental and numerical results showed that even though the reaction mechanisms employed have been optimized for predicting the laminar burning velocity of ammonia/air flames, they failed to satisfactorily predict the measured species in this study. Sensitivity analysis was used to identify elementary reactions that control the species profiles but have negligible effects on the burning velocity. It is considered that these reaction models need to be updated for accurate prediction of product gas characteristics of ammonia/air flames.     

On the combined effects of compositional inhomogeneity and ammonia addition to turbulent flames of ethylene

This paper is part of a broader program aimed at investigating the effects of co-firing clean fuels such as ammonia or hydrogen with hydrocarbons. The focus is on soot formation as well as flame stability in turbulent mixed-mode combustion, which is highly relevant in practical combustors. Ammonia substitution for nitrogen results in reduced flame stability, and this is correlated to differences in flame speed and extinction strain rate. While it is known that the addition of ammonia suppresses soot, visual inspection of compositionally inhomogeneous flames of ethylene-ammonia indicates a reduction in ammonia's ability to suppress soot formation. Measurements of soot volume fraction and laser-induced fluorescence in selected UV and visible bands are made along the centreline in selected flames to test this hypothesis. Experimental results are then compared to simulations in laminar diffusion flames, stratified counterflow flames, and partially premixed flames. All results confirm the soot-inhibiting ability of ammonia. Increasing inhomogeneity, leading to higher centreline mixture fractions, enhances soot formation, and the level of enhancement is greater for flames with ammonia than without. Moreover, it is found that partial premixing is ultimately responsible for determining the amount of soot formed as opposed to stratification of fuel mixtures near the pilot.     

Quantitative measurement of temperature in oxygen enriched CH4/O2/N2 premixed flames using Laser Induced Thermal Grating Spectroscopy (LITGS) up to 1.0 MPa

The application of laser diagnostics to high pressure combustion phenomena is particularly challenging, especially in practical combustors such as rocket motors. In this study, temperature measurements using Laser Induced Thermal Grating Spectroscopy (LITGS) are demonstrated in oxygen enriched CH₄/O₂/N₂ premixed laminar flames at pressures up to 1.0?MPa. We use a previously developed OH absorption LITGS technique to determine product gas temperatures from 0.3 to 1.0?MPa, for both high temperature oxygen-enriched and pure-oxygen flames, for measurements up to 3000?K. Further, we demonstrate how it is necessary to correct the measurements for the local absorption of laser light to obtain accurate temperatures, and offer a technique for producing the correction by using different laser energies. Once the correction is applied, we demonstrate that the measurements at 0.5?MPa are within 1.6% of the adiabatic non-strained flame temperatures, with a standard deviation of about 160?K, thus offering a competitive method for the challenging conditions at high pressures and temperatures. The values obtained at derived temperatures at 1.0?MPa were lower than the adiabatic unstrained flame temperatures, which could possibly be attributed to loss mechanisms.     

Effects of OH concentration and temperature on NO emission characteristics of turbulent non-premixed CH4/NH3/air flames in a two-stage gas turbine like combustor at high pressure

This study examined the effects of OH concentration and temperature on the NO emission characteristics of turbulent, non-premixed methane (CH₄)/ammonia (NH₃)/air swirl flames in two-stage combustors at high pressure. Emission data were obtained using large-eddy simulations with a finite-rate chemistry method from model flames based on the energy fraction of NH₃ (E_NH3) in CH₄/NH₃ mixtures. Although NO emissions at the combustor exit were found to be significantly higher than those generated by CH₄/air and NH₃/air flames under both lean and stoichiometric primary zone conditions, these emissions could be lowered to approximately 300 ppm by employing far-rich equivalence ratios (?) of 1.3 to 1.4 in the primary zone. This effect was possibly due to the lower OH concentrations under far-rich conditions. An analysis of local flame characteristics using a newly developed mixture fraction equation for CH₄/NH₃/air flames indicated that the local temperature and NO and OH concentration distributions with local ? were qualitatively similar to those in NH₃/air flames. That is, the maximum local NO and OH concentrations appeared at local ? of 0.9, although the maximum temperature was observed at local ? of 1.0. Both the temperature and OH concentration were found to gradually decrease with the partial replacement of CH₄ with NH₃. Consequently, NO emissions from CH₄/NH₃ flames were maximized at E_NH3 in the range of 20% to 30%, after which the emissions decreased. Above 2100 K, the NO emissions from CH₄/NH₃ flames increased exponentially with temperature, which was not observed in NH₃/air flames because of the lower flame temperatures in the latter. But, the maximum NO concentration in CH₄/NH₃ flames was occurred at a temperature slightly below the maximum temperature, just as in NH₃/air flames. The apparent exponential increase in NO emissions from CH₄/NH₃ flames is attributed to a similar trend in the OH concentration at high temperatures.     

Heat release rate as represented by [OH] × [CH2O] and its role in autoignition

Data from a recent instantaneous, simultaneous, high-resolution imaging experiment of Rayleigh temperature and laser induced fluorescence (LIF) of OH and CH₂O at the base of a turbulent lifted methane flame issuing into a hot vitiated coflow are analysed and contrasted to reference flames to further investigate the stabilization mechanisms involved. The use of the product of the quantified OH and semi-quantified CH₂O images as a marker for heat release rate is validated for transient autoigniting laminar flames. This is combined with temperature gradient information to investigate the flame structure. Super-equilibrium OH, the nature of the profiles of heat release rate with respect to OH mole fraction, and comparatively high peak heat release rates at low temperature gradients is found in the kernel structures at the flame base, and found to be indicative of autoignition stabilization.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

Chemical structure of atmospheric pressure premixed laminar formic acid/hydrogen flames

The work presents an experimental and kinetic modeling study of laminar premixed formic acid [HC(O)OH]/H₂/O₂/Ar flames at different equivalence ratios (φ=0.85, 1.1 and 1.3) stabilized on a flat burner at atmospheric pressure, as well as laminar flame speed of HC(O)OH/O₂/Ar flames (φ=0.5–1.5) at 1 atm. Flame structure as well as laminar flame speed were simulated using three different detailed chemical kinetic mechanisms proposed for formic acid oxidation. The components in the fuel blends show different consumption profiles, namely, hydrogen is consumed slower than formic acid. According to kinetic analysis, the reason of the observed phenomenon is that the studied flames have hydrogen as a fuel but also as an intermediate product formed from HC(O)OH decomposition. Comparison of the measured and simulated flame structure shows that all the mechanisms satisfactorily predict the mole fraction profiles of the reactants, main products, and intermediates. It is noteworthy that the mechanisms proposed by Glarborg et al., Konnov et al. and the updated AramcoMech2.0 adequately predict the spatial variations in the mole fractions of free radicals, such as H, OH O and HO₂. However, some drawbacks of the mechanisms used were identified; in particular, they predict different concentrations of CH₂O. As for laminar flame speed simulations, the Konnov et al. mechanism predicts around two times higher values than in experiment, while the Glarborg et al. and updated AramcoMech2.0 show good agreement with the experimental data.     

Laser-induced fluorescence determination of flame temperatures in comparison with CARS measurements

Temperature profiles in several premixed low pressure H₂/O₂/N₂ flames and in an atmospheric pressure CH₄/air flame were determined by laser-induced fluorescence (LIF) and by CARS experiments. In the LIF study, temperatures were derived from OH excitation spectra, CARS temperatures were deduced from N₂ Q-branch spectra. The present study is the first quantitative comparison of these two methods for temperature determination in flames burning at pressures up to 1 bar. The resulting temperatures showed good agreement.     

Numerical simulation of premixed H2–air cellular tubular flames

The detailed flame structure of laminar premixed cellular flames in the tubular domain is simulated in 2D using a fully-implicit primitive variable finite difference formulation that includes multicomponent transport and detailed chemical kinetics. Numerical results for H₂/air flames are presented and compared against spatially resolved experimental measurements of temperature and chemical species including atomic H and OH. The experimental results compare well for flame structure and cell number, despite the numerical model under-predicting the peak temperature by 200 K. Numerical experiments were performed to assess the ability for cellular tubular flames to impact experimental and numerical investigations of practical flames. The cellular flame structure is found to provide a highly sensitive geometry that is useful for validating diffusive transport modelling approximations. This capability is exemplified through the development of a simple and accurate approximation for thermal diffusion (i.e. the Soret effect) that is suitable for practical combustion codes.     

Effects of Soret diffusion on the laminar flame speed and Markstein length of syngas/air mixtures

The effects of Soret diffusion on premixed syngas/air flames at normal and elevated temperatures and pressures are investigated numerically including detailed chemistry and transport. The emphasis is placed on assessing and interpreting the influence of Soret diffusion on the unstretched and stretched laminar flame speed and Markstein length of syngas/air mixtures. The laminar flame speed and Markstein length are obtained by simulating the unstretched planar flame and positively-stretched spherical flame, respectively. The results indicate that at atmospheric pressure the laminar flame speed of syngas/air is mainly reduced by Soret diffusion of H radical while the influence of H₂ Soret diffusion is negligible. This is due to the facts that the main reaction zone and the Soret diffusion for H radical (H₂) are strongly (weakly) coupled, and that Soret diffusion reduces the H concentration in the reaction zone. Because of the enhancement in the Soret diffusion flux of H radical, the influence of Soret diffusion on the laminar burning flux increases with the initial temperature and pressure. Unlike the results at atmospheric pressure, at elevated pressures the laminar flame speed is shown to be affected by the Soret diffusion of H₂ as well as H radical. For stretched spherical flame, it is shown that the Soret diffusion of both H and H₂ should be included so that the stretched flame speed can be accurately predicted. Similar to the laminar flame speed, the Markstein length is also reduced by Soret diffusion. However, the reduction is found to be mainly caused by Soret diffusion of H₂ rather than that of H radical. Moreover, the influence of Soret diffusion on the Markstein length is demonstrated to decrease with the initial temperature and pressure.     

An investigation of flame zones and burning velocities of laminar unconfined methane-oxygen premixed flames

Numerical and experimental investigations of unconfined methane-oxygen laminar premixed flames are presented. In a lab-scale burner, premixed flame experiments have been conducted using pure methane and pure oxygen mixtures having different equivalence ratios. Digital photographs of the flames have been captured and the radial temperature profiles at different axial locations have been measured using a thermocouple. Numerical simulations have been carried out with a C2 chemical mechanism having 25 species and 121 reactions and with an optically thin radiation sub-model. The numerical results are validated against the experimental and numerical results for methane-air premixed flames reported in literature. Further, the numerical results are validated against the results from the present methane-oxygen flame experiments. Visible regions in digital flame photographs have been compared with OH isopleths predicted by the numerical model. Parametric studies have been carried out for a range of equivalence ratios, varying from 0.24 to 1.55. The contours of OH, temperature and mass fractions of product species such as CO, CO₂ and H₂O, are presented and discussed for various cases. By using the net methane consumption rate, an estimate of the laminar flame speed has been obtained as a function of equivalence ratio.     

Studies of autoignition-assisted nonpremixed cool flames

Autoignition-assisted nonpremixed cool flames of diethyl ether (DEE) are investigated in both laminar counterflow and turbulent jet flame configurations. First, the ignition and extinction limits of laminar nonpremixed cool flames of diluted DEE are measured and simulated using detailed kinetic models. The laminar flame measurements are used to validate the kinetic models and guide the turbulent flame measurements. The results show that, below a critical mixture condition, for elevated temperature and dilute mixtures, the cool flame extinction limit and the low-temperature ignition limit merge, leading to autoignition-assisted cool flame stabilization without hysteresis. Based on the findings from the laminar flame experiments, autoignition-assisted turbulent lifted cool flames are established using a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner. The lift-off heights of the turbulent cool flames are quantified using formaldehyde planar laser-induced fluorescence. Based on an analogy with autoignition-assisted lifted hot flames, a correlation is proposed such that the autoignition-assisted cool flame lift-off height scales with the product of the flow velocity and the square of the first-stage ignition delay time. Using this scaling, we demonstrate that the kinetic mechanism that most accurately predicts the laminar flame ignition and extinction limits also best predicts the turbulent cool flame lift-off height.     

CO2 addition and pressure effects on laminar and turbulent lean premixed CH4 air flames

An experimental study on CH₄–CO₂–air flames at various pressures is conducted by using both laminar and turbulent Bunsen flame configurations. The aim of this research is to contribute to the characterization of fuel lean methane/carbon dioxide/air premixed laminar and turbulent flames at different pressures, by studying laminar and turbulent flame propagation velocities, the flame surface density and the instantaneous flame front wrinkling parameters. PREMIX computations and experimental results indicate a decrease of the laminar flame propagation velocities with increasing CO₂ dilution rate. Instantaneous flame images are obtained by Mie scattering tomography. The image analysis shows that although the height of the turbulent flame increases with the CO₂ addition rate, the flame structure is quite similar. This implies that the flame wrinkling parameters and flame surface density are indifferent to the CO₂ addition. However, the pressure increase has a drastic effect on both parameters. This is also confirmed by a fractal analysis of instantaneous images. It is also observed that the combustion intensity S_T/S_L increases both with pressure and the CO₂ rate. Finally, the mean fuel consumption rate decreases with the CO₂ addition rate but increases with the pressure.     

Spectroscopic observations of methane-pulverized coal flames

A premixed laminar flame burner was used to study physical and chemical effects of burning methane and pulverized coal simultaneously. Spectral emissions obtained as a function of height in the flame, for OH, CH, C₂, and CO were used to calculate excited state species number densities. Premixed methane-air flames with a number of equivalence ratios and methane-air-coal dust flames with three equivalence ratios were studied and differences were noted. Spatial temperatures were measured in all flames by means of a thermocouple probe.The introduction of a pulverized coal additive (30mg/min) leads to changes in OH, Ch, C₂, and CO number densities and temperature at various heights, which depend on the stoichiometry of the methane-air flame. Possible chemical and physical mechanisms responsible for these changes are discussed qualitatively. The data pose a severe test for any quantitative model which includes the fluid dynamics, gas-phase chemistry, heterogeneous processes, and radiation transport.     

On laminar premixed flame propagating into autoigniting mixtures under engine-relevant conditions

Usually premixed flame propagation and laminar burning velocity are studied for mixtures at normal or elevated temperatures and pressures, under which the ignition delay time of the premixture is much larger than the flame resistance time. However, in spark-ignition engines and spark-assisted compression ignition engines, the end-gas in the front of premixed flame is at the state that autoignition might happen before the mixture is consumed by the premixed flame. In this study, laminar premixed flames propagating into an autoigniting dimethyl ether/air mixture are simulated considering detailed chemistry and transport. The emphasis is on the laminar burning velocity of autoigniting mixtures under engine-relevant conditions. Two types of premixed flames are considered: one is the premixed planar flame propagating into an autoigniting DME/air without confinement; and the other is premixed spherical flame propagating inside a closed chamber, for which four stages are identified. Due to the confinement, the unburned mixture is compressed to high temperature and pressure close to or under engine-relevant conditions. The laminar burning velocity is determined from the constant-volume propagating spherical flame method as well as PREMIX. The laminar burning velocities of autoigniting DME/air mixture at different temperatures, pressures, and autoignition progresses are obtained. It is shown that the first-stage and second-stage autoignition can significantly accelerate the flame propagation and thereby greatly increase the laminar burning velocity. When the first-stage autoignition occurs in the unburned mixture, the isentropic compression assumption does not hold and thereby the traditional method cannot be used to calculate the laminar burning velocity. A modified method without using the isentropic compression assumption is proposed. It is shown to work well for autoigniting mixtures. Besides, a power law correlation is obtained based on all the laminar burning velocity data. It works well for mixtures before autoignition while improvement is still needed for mixtures after autoignition.