Effect of flame retardants on side-wall quenching of partially premixed laminar flames

A combined experimental and numerical investigation of partially premixed laminar methane-air flames undergoing side-wall quenching (SWQ) is performed. A well-established SWQ burner is adapted to allow the seeding of the main flow with additional gaseous products issued from a (secondary) wall inlet close to the flame’s quenching point. First, the characteristics of the partially premixed flame that quenches at the wall are assessed using planar laser-induced fluorescence measurements of the OH radical, and a corresponding numerical simulation with fully-resolved transport and chemistry is conducted. A boundary layer of enriched mixture is formed at the wall, leading to a reaction zone parallel to the wall for high injection rates from the wall inlet. Subsequently, in a numerical study, the wall inflow is mixed with dimethylmethylphosphonat (DMMP), a phosphor-based flame retardant. The DMMP addition allows the assessment of the combined effects of heat loss and flame retardants on the flame structure during flame-wall interaction. With an increasing amount of DMMP in the injected mixture, the flame stabilizes further away from the wall and shows a decrease in the local heat-release rate. Thereby, the maximum wall heat flux is significantly reduced. That results in a lower thermal load on the quenching wall. The flame structure analysis shows an accumulation of the intermediate species HOPO at the wall similar to the CO accumulation during the quenching of premixed flames without flame retardant addition. The study shows how the structure of a partially premixed flame is influenced by a wall that releases either additional fuel or a mixture of fuel and flame retardant. The insights gained from the canonical configuration can lead to a better understanding of the combined effects of flame retardants and heat losses in near-wall flames.     

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.     

Turbulent burning velocity and its related statistics of ammonia‐hydrogen‐air jet flames at high Karlovitz number: Effect of differential diffusion

To clarify the role of differential diffusion in highly turbulent premixed flames, a series of turbulent premixed ammonia/hydrogen/air flames were investigated using the NH-PLIF diagnostics. The investigated flames have almost the same laminar burning velocity, S_L, but are characterized by different Lewis number, Le, from 0.56 to 1.77. The Karlovitz number, Ka, of these flames ranges from 11 to 1052, and the turbulence intensity, u’/S_L, covers from 10 to 156. It is observed that the global consumption speed, S_T,GC/S_L, of sub-unity Le flames is much larger than that of super-unity Le flames at high Ka, indicating that the differential diffusion plays a significant role in highly turbulent flames. The flame surface density and the area ratio of turbulent flames with different Le are, however, similar under wide turbulent conditions. The stretch factor of sub-unity Le flames is estimated to be significantly larger than that of super-unity Le cases. The enhanced S_T,GC of sub-unity Le flames is suggested to be attributed to the promotion of local burning rates by the couple effect of differential diffusion and turbulent flame stretch within the flame brush, rather than the enlargement of flame surface area at high Ka. Furthermore, three correlations for the S_T,GC were developed based on Damkohler's second hypothesis with consideration of the Le effect. The correlation of S_T,GC/S_L ∼ (Re_T·Le^-2)^0.5 is further validated by using small-scale methane/air and large-scale ammonia/air flames at high Ka, where Re_T is turbulent Reynolds number. It suggests that the S_T,GC is roughly inversely proportional to the Le, and the differential diffusion effect should be included in the theoretical analysis and numerical simulation of highly turbulent flames.     

Fundamental burning velocities of meso-scale propagating spherical flames with H2, CH4 and C3H8 mixtures

This study is performed to experimentally examine the fundamental burning velocity characteristics of meso-scale outwardly propagating spherical laminar flames in the range of flame radius r_f approximately from 1 to 5 mm for hydrogen, methane and propane mixtures, in order to make clear a method for improving combustion of micro–meso scale flames. Macro-scale laminar flames with r_f > 7 mm are also examined for comparison. The mixtures have nearly the same laminar burning velocity (S_L0 = 25 cm/s) for unstretched flames and different equivalence ratios ?. The radius r_f and the burning velocity S_Ll of meso-scale flames are estimated by using sequential schlieren images recorded under appropriate ignition conditions. It is found that S_Ll of hydrogen and methane premixed meso-scale flames at the same r_f or the Karlovitz number Ka shows a tendency to increase with decreasing ?, whereas S_Ll of propane flames increases with ?. However, S_Ll tends to decrease with the Lewis number Le and the Markstein number Ma, irrespective of the type of fuel and ?. It also becomes clear that the optimum flame size and Ka to improve the burning velocity exist for some mixtures depending on Le and fuel types.     

Problems of predicting turbulent burning rates

Different approaches to the modelling of turbulent combustion first are reviewed briefly. A unified, stretched flamelet approach then is presented. With Reynolds stress modelling and a generalized probability density function (PDF) of strain rate, it enables a source term, in the form of a probability of burning function, P_b, to be expressed as a function of Markstein numbers and the Karlovitz stretch factor. When P_b is combined with some turbulent flame fractal considerations, an expression is obtained for the turbulent burning velocity. When it is combined with the profile of the unstretched laminar flame volumetric heat release rate plotted against the reaction progress variable and the PDF of the latter, an expression is obtained for the mean volumetric turbulent heat release rate. Through these relationships experimental values of turbulent burning velocity might be used to evaluate P_b and hence the CFD source term, the mean volumetric heat release rate.

Different theoretical expressions for the turbulent burning velocity, including the present one, are compared with experimental measurements. The differences between these are discussed and this is followed by a review of CFD applications of these flamelet concepts to premixed and non-premixed combustion. The various assumptions made in the course of the analyses are scrutinized in the light of recent direct numerical simulations of turbulent flames and the applications to the flames of laser diagnostics. Remaining problem areas include a sufficiently general combination of strain rate and flame curvature PDFs to give a single PDF of flame stretch rate, the nature of flame quenching under positive and negative stretch rates, flame responses to changing stretch rates and the effects of flame instabilities.     

Framework for simulating stationary spherical flames

Understanding and quantifying the effects of flame stretch rate on the laminar flame speed and flame structure plays an important role from interpreting experimentally-measured laminar burning velocities to characterizing the impact of turbulence on premixed flames. Unfortunately, accounting for these effects often requires an unsteady reacting flow solver and may be computationally expensive. In this work, we propose a mathematical framework to perform simulations of stationary spherical flames. The objective is to maintain the flame at a constant radius (and hence a constant stretch rate) by performing a coordinate change. The governing equations in the new flame-attached frame of reference resemble the original equations for freely-propagating spherical flames. The only difference is the presence of additional source terms whose purpose is to drive the numerical solution to a steady state. These source terms involve one free parameter: the flame stretch rate, which may either be computed in real time or imposed by the user. This parameter controls ultimately the steady state flame radius and the steady state flame speed. That is why, at a given stretch rate, the results of the stationary spherical flame simulations match those of a freely-expanding spherical flame. As an illustration, the dependence of the laminar flame speed on the stretch rate is leveraged to extract Markstein lengths for hydrogen/air mixtures at different equivalence ratios, as well as for hydrocarbon/air mixtures (CH₄ and C₇H₁₆). Numerical predictions are in good agreement with experimental measurements (within experimental uncertainties). Finally, the proposed methodology is implemented in the chemical kinetic software FlameMaster. The use of a dedicated steady-state solver with a non-uniform optimized mesh leads to significant reductions in the computational cost, highlighting that the proposed methodology is ideally suited for other chemical kinetic software such as Chemkin/Premix and Cantera.     

Effect of Lewis number on premixed laminar lean-limit flames stabilized on a bluff body

In this paper, we present a study on the effect of Lewis number, Le, on the stabilization and blow-off of laminar lean limit premixed flames stabilized on a cylindrical bluff body. Numerical simulations and experiments are conducted for propane, methane and two blends of hydrogen with methane as fuel gases, containing 20% and 40% of hydrogen by volume, respectively. It is found that the Le?>?1 flame blows-off via convection from the base of the flame (without formation of a neck) when the conditions for flame anchoring are not fulfilled. Le?≤?1 flames exhibit a necking phenomenon just before lean blow-off. This necking of the flame front is a result of the local reduction in mass burning rates causing flame merging and quenching of the thin flame tube formed. The structure of these flames at the necking location is found to be similar to tubular flames. It is found that extinction stretch rates for tubular flames closely match values at the neck location of bluff-body flames of corresponding mixtures, suggesting that excessive flame stretch is directly responsible for blow-off of the studied Le?≤?1 flames. After quenching of the neck, the upstream part forms a steady and stable residual flame in the wake of the bluff body while the downstream part is convected away.     

Characterization of the effects of pressure and hydrogen concentration on laminar burning velocities of methane–hydrogen–air mixtures

The aim of the present work was to characterize both the effects of pressure and of hydrogen addition on methane/air premixed laminar flames. The experimental setup consists of a spherical combustion chamber coupled to a classical shadowgraphy system. Flame pictures are recorded by a high speed camera. Global equivalence ratios were varied from 0.7 to 1.2 for the initial pressure range from 0.1 to 0.5 MPa. The mole fraction of hydrogen in the methane + hydrogen mixture was varied from 0 to 0.2. Experimental results were compared to calculations using a detailed chemical kinetic scheme (GRIMECH 3.0). First, the results for atmospheric laminar CH₄/air flames were compared to the literature. Very good agreements were obtained both for laminar burning velocities and for burned gas Markstein length. Then, increasing the hydrogen content in the mixture was found to be responsible for an increase in the laminar burning velocity and for a reduction of the flame dependence on stretch. Transport effects, through the reduction of the fuel Lewis number, play a role in reducing the sensitivity of the fundamental flame velocity to the stretch. Finally, when the pressure was increased, the laminar burning velocity decreased for all mixtures. The pressure domain was limited to 0.5 MPa due to the onset of instabilities at pressures above this value.     

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.     

The stability of laminar explosion flames

The paper presents the results of a fundamental experimental and theoretical study of Darrieus–Landau, thermo-diffusive, instabilities in atmospheric explosions, and, on a smaller scale, in laboratory explosions in closed vessels. Pressure dependencies were sought to exploit the leading role of the Peclet number in the phenomena, so that similar Peclet numbers were achieved in both instances. However, in large atmospheric explosions large Peclet numbers were achieved by the size of the fireball, whereas in the closed vessel explosion it was achieved at a higher pressure by a much smaller flame, but because of the higher pressure, one endowed with a small laminar flame thickness. This study covers a much wider range of fuels and of pressures and the dependencies of the phenomena on both of these were carefully studied, although, for the atmospheric explosions, the data only covered propane and methane. The roles of both Markstein and Peclet numbers become clear and give rise to a more fundamental correlating parameter, a critical Karlovitz number, K_cl, for flame stability. This is based on the flame stretch rate, normalised by its multiplication by the chemical reaction time in a laminar flame. The experimentally measured dependencies of this key parameter on pressure and Markstein number are reported for the first time for so many different fuels. The critical Karlovitz number for flame stability decreases with increase in the strain rate Markstein number. As a result, it is possible to predict the extent of the unstable regime for laminar flames as a function of Ma_sr and pressure. Such data can be used to estimate the severity of large scale atmospheric explosions. As Ma_sr becomes highly negative, the regime of stability is markedly reduced.     

Laser-based investigation of flame surface density and mean reaction rate during flame-wall interaction at elevated pressure

Velocities and flame front locations are measured simultaneously in a turbulent, side-wall quenching (SWQ) V-shaped flame during flame-wall interaction (FWI) at 1 and 3 bar by means of particle image velocimetry (PIV) and planar laser-induced fluorescence of the OH radical (OH-PLIF). The turbulent flame brush is characterized based on the spatial distribution of the mean reaction progress variable and a common direct method is used to derive the flame surface density (FSD) from the two-dimensional data by image processing. As the near-wall reaction zone is limited to a smaller region closer to the wall at higher pressure, higher peak values are observed in the FSD at 3 bar. A second definition of the FSD adapted for flames exposed to quenching is utilized similar to previous studies emphasizing the impact of FWI. The influence of the wall on the flame front topology is investigated based on a flame front-conditioned FSD and its variability within the data set. In a last step, an estimate of the mean reaction rate is deduced using an FSD model and evaluated in terms of integral and space-averaged values. A decreasing trend of integral mean reaction rate in regions with increasing flame quenching is observed for both operating conditions, but more pronounced at 3 bar. Space-averaged mean reaction rates, however, increase in the quenching region, as the size of the reaction zone decreases.     

An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

A physical relationship between consumption and displacement speed for premixed flames with finite thickness

The characterization of premixed flames by a flame speed has been a subject that has occupied much interest in the literature in many systematic studies on combustion phenomena. Consumption and displacement speeds are two such flame speeds that are understood to describe the flame dynamics under the effect of flame curvature, flow non-uniformities, Lewis number and turbulence effects along with heat transfer with flame holders and cold walls. As such, much work has been done in the past where either one of these two speeds has been employed along with a linear sensitivity coefficient (Markstein length) for describing different sensitivities to stretch effects. However, despite recent attempts using the asymptotic theory, the relationship between these two quantities has only been clarified in a limited manner for flames of finite thickness. In this study, we use flame stretch theory that takes into account changes of stretch, curvature, heat transfer and Lewis number effects throughout the pre-heat zone and its integral effect on the flame reaction zone. A sound mathematical and physical basis is provided for understanding the two speeds that is valid for weak as well as strong stretch effects. Understanding from theory is further demonstrated by analysing several example 1D stretched flames along with a 2D bluff body flame near extinction.     

How fast can we burn, 2.0

We review the state of the art in measurements and simulations of the behavior of premixed laminar and turbulent flames, subject to differential diffusion, stretch and curvature. The first part of the paper reviews the behavior of premixed laminar flames subject to flow stretch, and how it affects the accuracy of measurements of unstrained laminar flame speeds in stretched and spherically propagating flames. We then examine how flow field stretch and differential diffusion interact with flame propagation, promoting or suppressing the onset of thermodiffusive instabilities. Secondly, we survey the methodology for and results of measurements of turbulent flame speeds in the light of theory, and identify issues of consistency in the definition of mean flame speeds, and their corresponding mean areas. Data for methane at a single operating condition are compared for a range of turbulent conditions, showing that fundamental issues that have yet to be resolved for Bunsen and spherically propagating flames. Finally, we consider how the laminar flame scale response of flames to flow perturbations interacting with differential diffusion leads to very different outcomes to the overall sensitivity of the burning rate to turbulence, according to numerical simulations (DNS). The paper concludes with opportunities for future measurements and model development, including the perennial recommendation for robust archival databases of experimental and DNS results for future testing of models.     

Extinction and flame bifurcations of stretched dimethyl ether premixed flames

Extinction limits and flame bifurcation of lean premixed dimethyl ether–air flames are numerically investigated using the counterflow flame with a reduced chemistry. Emphasis is paid to the combined effect of radiation and flame stretch on the extinction and flammability limits. A method based on the reaction front is presented to predict the Markstein length. The predicted positive Markstein length agrees well with the experimental data. The results show that flow stretch significantly reduces the flame speed and narrows the flammability limit of the stretched dimethyl ether–air flame. It is found that the combined effect of radiation and flow stretch results in a new flame bifurcation and multiple flame regimes. At an equivalence ratio slightly higher than the flammability limit of the planar flame, the distant flame regime appears at low stretch rates. With an increase in the equivalence ratio, in addition to the distant flame, a weak flame isola emerges at moderate stretch rates. With a further increase in the equivalence ratio, the distant flame and the weak flame branches merge together, resulting in the splitting of the weak flame branch into two weak flame branches, one at low stretch and the other at high stretch. Flame stability analysis demonstrates that the high stretch weak flame is also stable. Furthermore, a K-shaped flammability limit diagram showing various flame regimes and their extinction limits is obtained.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

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.