Transported PDF simulation of turbulent CH4/H2 flames under MILD conditions with particle-level sensitivity analysis

Transported probability density function (TPDF) simulation with sensitivity analysis has been conducted for turbulent non-premixed CH₄/H₂ flames of the jet-into-hot-coflow (JHC) burner, which is a typical model to emulate moderate or intense low oxygen dilution combustion (MILD). Specifically, two cases with different levels of oxygen in the coflow stream, namely HM1 and HM3, are simulated to reveal the differences between MILD and hot-temperature combustion. The TPDF simulation well predicts the temperature and species distributions including those of OH, CO and NO for both cases with a 25-species mechanism. The reduced reaction activity in HM1 as reflected in the peak OH concentration is well correlated to the reduced oxygen in the coflow stream. The particle-level local sensitivities with respect to mixing and chemical reaction further show dramatic differences in the flame characteristics. HM1 is less sensitive to mixing and reaction parameters than HM3 due to the suppressed combustion process. Specifically, for HM1 the sensitivities to mixing and chemical reactions have comparable magnitude, indicating that the combustion progress is controlled by both mixing and reaction in MILD combustion. For HM3, there is however a change in the combustion mode: during the flame initialization, the combustion progress is more sensitive to chemical reactions, indicating that finite-rate chemistry is the controlling process during the autoignition process for flame stabilization; at further downstream where the flame has established, the combustion progress is controlled by mixing, which is characteristic of nonpremixed flames. An examination of the particles with the largest sensitivities reveals the difference in the controlling mixtures for flame stabilization, namely, the stoichiometric mixtures are important for HM1, whereas, fuel-lean mixtures are controlling for HM3. The study demonstrates the potential of TPDF simulations with sensitivity analysis to investigate the effects of finite-rate chemistry on the flame characteristics and emissions, and reveal the controlling physio-chemical processes in MILD combustion.     

Generalisation of the eddy-dissipation concept for jet flames with low turbulence and low Damköhler number

Moderate or intense low oxygen dilution (MILD) combustion has been the focus of a range of fundamental experimental and numerical studies. Reasonable agreement between experimental and numerical investigations, however, requires finite-rate chemistry models and, often, ad hoc model adjustment. To remedy this, an adaptive eddy dissipation concept (EDC) combustion model has previously been developed to target conditions encountered in MILD combustion; however, this model relies on a simplified, pre-defined assumption about the combustion chemistry. The present paper reports a generalised version of the modified EDC model without the need for an assumed, single-step chemical reaction or ad hoc coefficient tuning. The results show good agreement with experimental measurements of two CH₄/H₂ flames in hot coflows, showing improvements over the standard EDC model as well as the previously published modified EDC model. The updated version of the EDC model also demonstrates the capacity to reproduce the downstream transition in flame structure of a MILD jet flame seen experimentally, but which has previously proven challenging to capture computationally. Analyses of the previously identified dominant heat-release reactions provide insight into the structural differences between a conventional autoignitive flame and a flame in the MILD combustion regime, whilst highlighting the requirement for a generalised EDC combustion model.     

Strain rate and fuel composition dependence of chemiluminescent species profiles in non-premixed counterflow flames: comparison with model results

A detailed comparison has been conducted between chemiluminescence (CL) species profiles of OH^?, CH^?, and C₂ ^?, obtained experimentally and from detailed flame kinetics modeling, respectively, of atmospheric pressure non-premixed flames formed in the forward stagnation region of a fuel flow ejected from a porous cylinder and an air counterflow. Both pure methane and mixtures of methane with hydrogen (between 10 and 30 % by volume) were used as fuels. By varying the air-flow velocities methane flames were operated at strain rates between 100 and 350 s^?1, while for methane/hydrogen flames the strain rate was fixed at 200 s^?1. Spatial profiles perpendicular to the flame front were extracted from spectrograms recorded with a spectrometer/CCD camera system and evaluating each spectral band individually. Flame kinetics modeling was accomplished with an in-house chemical mechanism including C₁–C₄ chemistry, as well as elementary steps for the formation, removal, and electronic quenching of all measured active species. In the CH₄/air flames, experiments and model results agree with respect to trends in profile peak intensity and position. For the CH₄/H₂/air flames, with increasing H₂ content in the fuel the experimental CL peak intensities decrease slightly and their peak positions shift towards the fuel side, while for the model the drop in mole fraction is much stronger and the peak positions move closer to the fuel side. For both fuel compositions the modeled profiles peak closer to the fuel side than in the experiments. The discrepancies can only partly be attributed to the limited attainable spatial resolution but may also necessitate revised reaction mechanisms for predicting CL species in this type of flame.     

The effect of fuel inlet turbulence intensity on H2/CH4 flame structure of MILD combustion using the LES method

Large eddy simulations (LES) are employed to investigate the effect of the inlet turbulence intensity on the H₂/CH₄ flame structure in a hot and diluted co-flow stream which emulates the (Moderate or Intense Low-oxygen Dilution) MILD combustion regime. In this regard, three fuel inlet turbulence intensity profiles with the values of 4%, 7% and 10% are superimposed on the annular mixing layer. The effects of these changes on the flame structure under the MILD condition are studied for two oxygen concentrations of 3% and 9% (by mass) in the oxidiser stream and three hot co-flow temperatures 1300, 1500 and 1750 K. The turbulence-chemistry interaction of the numerically unresolved scales is modelled using the (Partially Stirred Reactor) PaSR method, where the full mechanism of GRI-2.11 represents the chemical reactions. The influences of the turbulence intensity on the flame structure under the MILD condition are studied by using the profile of temperature, CO and OH mass fractions in both physical and mixture fraction spaces at two downstream locations. Also, the effects of this parameter are investigated by contours of OH, HCO and CH₂O radicals in an area near the nozzle exit zone. Results show that increasing the fuel inlet turbulence intensity has a profound effect on the flame structure particularly at low oxygen mass fraction. This increment weakens the combustion zone and results in a decrease in the peak values of the flame temperature and OH and CO mass fractions. Furthermore, increasing the inlet turbulence intensity decreases the flame thickness, and increases the MILD flame instability and diffusion of un-burnt fuel through the flame front. These effects are reduced by increasing the hot co-flow temperature which reinforces the reaction zone.     

Evaluation of heat release indicators in lean premixed H2/Air cellular tubular flames

Heat release rate in combustion systems must be understood in order to control thermoacoustic instabilities, flame extinction, and heat losses. Traditionally OH chemiluminescence (OH*) is used to trace heat release rate (HRR) in H₂/air flames, but its accuracy as a tracer has not been assessed. Lean premixed H₂/air cellular tubular flames are a good test case to evaluate HRR tracers due to the presence of highly reactive flame cells surrounded by regions of near extinction. Comparing the calculated heat release rate to OH* concentration, one finds that [OH*] profiles correlate with the regions of high reactivity (flame cells) but the correlation fails in the low reactivity regions where the HRR is much higher than the [OH*] value indicates. Alternate HRR tracers including [H] and pixel-by-pixel products of [O₂]x[H], [OH]x[H₂], and [O]x[H₂] are analyzed with detailed numerical simulations. The chosen products derive from the main chain reaction steps that contribute to overall HRR in lean, premixed H₂/air flames. Findings suggest that [H] is an accurate yet simple way of tracking HRR. Planar measurements of HRR are possible if LIF measurements of [H] are improved.     

Flame structure analysis and flamelet progress variable modelling of strained coal flames

Strained two-phase pulverised coal flames in a counterflow configuration are investigated numerically. Three operating conditions with different coal-to-primary-air ratios and inlet velocities were evaluated in order to establish different flame regimes. At first, the two-phase flow of the fully resolved reference cases is calculated solving the transport equation for the species and directly evaluating the reaction rates. Different flame structures are identified using the heat release rate and the chemical explosive mode as markers, showing that complex structures with a combination of lean premixed and non-premixed flames can be observed in strained counterflow coal flames. In addition to the fully resolved simulation, the suitability of the Flamelet-Progress Variable (FPV) model is investigated. Both premixed and non-premixed tables are employed. At first, the suitability of the look-up tables is evaluated by means of an a priori analysis, using the fully resolved simulations as reference solutions, showing that the non-premixed flamelet table correctly predicts the structure of the strained coal flames, while the premixed table shows sensible deviations in terms of temperature and species, especially at rich conditions. Finally, the a posteriori analysis shows that the fully coupled FPV model with a non-premixed flamelet look-up table can accurately predict strained coal flames.     

Numerical study of the effect of turbulence on rate of reactions in the MILD combustion regime

In this paper, the importance of fluctuations in flow field parameters is studied under MILD combustion conditions. In this way, a turbulent non-premixed CH₄+H₂ jet flame issuing into a hot and deficient co-flow air is modeled using the RANS Axisymmetric equations. The modeling is carried out using the EDC model to describe the turbulence-chemistry interaction. The DRM-22 reduced mechanism and the GRI2.11 full mechanism are used to represent the chemical reactions of H₂/methane jet flame. Results illustrate that although the fluctuations in temperature field are small and the reaction zone volume are large in the MILD regime, the fluctuations in temperature and species concentrations are still effective on the flow field. Also, inappropriate dealing with the turbulence effect on chemistry leads to errors in prediction of temperature up to 15% in the present flame. By decreasing of O₂ concentration of hot co-flow air, the effect of fluctuations in flow field parameters on flame characteristics are still significant and its effect on species reaction rates does not decrease. On the other hand, although decreasing of jet inlet Reynolds number at constant inlet turbulence intensity addresses to smaller fluctuations in flow filed, it does not lead to lower the effect of turbulence on species distribution and temperature field under MILD combustion conditions.     

Simultaneous measurement of localized heat release with OH/CH2O-LIF imaging and spatially integrated OH∗ chemiluminescence in turbulent swirl flames

The in-situ and localized observation of heat release in turbulent flames is important for the validation of computational modeling of turbulent flows with combustion. In the present work we obtain localized information on heat release rate (HRR) by the commonly accepted technique of the simultaneous and single-shot planar imaging of OH and CH₂O concentrations by laser-induced fluorescence (LIF). Additionally, we combine this with the simultaneous line-of-sight and temporally resolved chemiluminescence detection of OH^?, spatially integrated within the flame volume, interrogated by the laser sheets used for the HRR imaging technique. The combined diagnostic methods are demonstrated for a swirl-stabilized, premixed turbulent methane/air flame of 30-kW thermal power, and they show the existence of correlations between both HRR-sensitive diagnostic techniques.     

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

The growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O₂]ⁿ, and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.     

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.     

Simultaneous stereo‐PIV and OH×CH2O PLIF measurements in turbulent ultra lean CH4/H2 swirling wall‐impinging flames

We present experimental results from turbulent low-swirl lean H₂/CH₄ flames impinging on an inclined, cooled iso-thermal wall, based on simultaneous stereo-PIV and OH×CH₂O PLIF measurements. By increasing the H₂ fraction in the fuel while keeping Karlovitz number (Ka) fixed in a first series of flames, a fuel dependent near-wall flame structure is identified. Although Ka is constant, flames with high H₂ fraction exhibit significantly more broken reaction zones. In addition, these high H₂ fraction flames interact significantly more with the wall, stabilizing through the inner shear layer and well inside the near-wall swirling flow due to a higher resistance to mean strain rate. This flame-wall interaction is argued to increase the effective local Ka due to heat loss to the wall, as similar flames with a (near adiabatic) ceramic wall instead of a cooled wall exhibit significantly less flame brokenness. A second series of leaner flames were investigated near blow-off limit and showed complete quenching in the inner shear layer, where the mean strain rate matches the extinction strain rate extracted from 1D flames. For pure CH₄ flames (Ka ≈ 30), the reaction zone remains thin up to the quenching point, while conversely for the 70% H₂ flames (Ka ≈ 1100), the reaction zone is highly fragmented. Remarkably, in all near blow-off cases with CH₄ in the fuel, a large cloud of CH₂O persists downstream the quenching point, suggesting incomplete combustion. Finally, ultra lean pure hydrogen flames were also studied for equivalence ratios as low as 0.22, and through OH imaging, exhibit a clear transition from a cellular flame structure to a highly fragmented flame structure near blow-off.     

Flame studies of C2 hydrocarbons

The oxidation characteristics of C₂ hydrocarbons were revisited in flames established in the counterflow configuration. Laminar flame speeds of ethane/air, ethylene/air, and acetylene/oxygen/nitrogen mixtures as well as extinction strain rates of non-premixed ethane/air flames were measured using digital particle image velocimetry. The experiments were modeled using three different kinetic models. While the experimental and computed laminar flame speeds agreed closely for all C₂ hydrocarbons under fuel-lean conditions, notable discrepancies were identified under fuel-rich conditions. Using the computed flame structures, insight was provided into the controlling mechanisms that could be responsible for the observed discrepancies. More specifically, the uncertainties associated with the kinetics of the thermal decomposition of the ethyl radical were found to be a potential source of the observed discrepancies for ethane flames. It was shown also by using alternative rate constants for the ethyl radical decomposition, the rate of flame propagation and the extinction propensity are affected notably. Furthermore, the values of the branching ratio of acetylene consumption reactions involving atomic oxygen were found to have a significant effect on the propagation of rich acetylene flames.     

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.     

Role of radicals on MILD combustion inception

The criterion used to define MILD combustion in non-premixed condition is analysed using Direct Numerical Simulation (DNS) of MILD combustion of methane-diluted air established with internal exhaust gas recirculation. The simulations reveal multiple interacting reaction zones in MILD combustion which are extremely different from conventional combustion. Furthermore, DNS deduced S-curves highlight the role of chemically active species. Specifically, the temperature rise is accompanied with an increase in the scalar dissipation rate of mixture fraction, which is quite contrasting to the classical S-curve from the classical flame theories. This observation is explained on a physical basis.     

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.     

Near blowout dynamics of a premixed,swirl stabilized flame

This paper analyzes the flame dynamics of near LBO (lean blowoff) swirl stabilized flames, using simultaneous OH and CH₂O (formaldehyde) PLIF (planar laser induced fluorescence) measurements. Prior studies have shown that recirculation stabilized flames approach blowoff through two distinct stages – “stage 1” characterized by local extinction, where the overall flame and flow field remain largely unchanged, and “stage 2”, characterized by a fundamental change in the flow field, accompanied by violent flame flapping and wake disruption. This paper quantifies extinction spots along the flame edge, and entrained reactants within the combustion product region to analyze these stages in greater detail. Extinction spots were quantified by the overlapping regions of OH and CH₂O – numerous such spots were found near blowoff. The entraining of unburnt reactants into the recirculation zone was quantified by detecting low intensity OH pockets that were not surrounded by CH₂O. As expected, the flame near blowoff displayed significantly more entrained reactant pockets relative to a stable flame. Unexpectedly however, the volume of these pockets is tiny compared to the products, even on the edge of blowoff. Once they enter the wake, they are short-lived, suggesting that they are diluted and/or quickly react. This was surprising given the non-trivial baseflow and flame position disruption at these conditions, suggesting a striking similarity between the average composition of the wake, to that of a stable flame.     

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.     

Computational study on lean and rich combustion of flame ball,counterflow flame and planar flame: Their limits and stoichiometries

To investigate (fuel-)lean/rich limits and essential stoichiometries, i.e., the borders of lean/rich combustion, one-dimensional steady computations with detailed chemistry for flame balls, counterflow flames, and stretch-free planar flames were conducted using a CH₄/O₂/Xe mixture that has been used in microgravity experiments. As continuous converged solutions were obtained under lean/rich conditions, it was suggested that the existence of flame ball not only under lean but also under rich condition. Flame radii and temperatures of flame balls decreased and increased toward the lean/rich limits from their maximum and minimum values, respectively. The lean limits were wider in the order of the flame ball, counterflow flame, and stretch-free planar flame. Therefore, the lean flammability limit corresponded to the lean limit of the flame ball in the mixture. Conversely, the rich limits were wider in the order of the counterflow flame, stretch-free planar flame, and flame ball. Thus, the rich flammability limit corresponded to the rich limit of the counterflow flame in the mixture. Essential stoichiometry, which represents the actual stoichiometry depending on the dominant transport in near-flame front, was not uniquely determined as conventional stoichiometry (ϕ = 1); it was located between the equivalence ratio of ϕ = 1 and ϕ_c, where ϕ _c denotes the critical equivalence ratio is evaluated using the fuel and oxidizer Lewis number of a target mixture. The results indicated that the essential stoichiometry of the stretch-free planar flame corresponded to ϕ = 1, that of the flame ball corresponded to ϕ = ϕ _c, and that of the stretched flame was located between ϕ = 1 and ϕ _c depending on the stretch rate.