Bifurcation and extinction limit of stretched premixed flames with chain-branching intermediate kinetics and radiative loss

Premixed counterflow flames with thermally sensitive intermediate kinetics and radiation heat loss are analysed within the framework of large activation energy. Unlike previous studies considering one-step global reaction, two-step chemistry consisting of a chain branching reaction and a recombination reaction is considered here. The correlation between the flame front location and stretch rate is derived. Based on this correlation, the extinction limit and bifurcation characteristics of the strained premixed flame are studied, and the effects of fuel and radical Lewis numbers as well as radiation heat loss are examined. Different flame regimes and their extinction characteristics can be predicted by the present theory. It is found that fuel Lewis number affects the flame bifurcation qualitatively and quantitatively, whereas radical Lewis number only has a quantitative influence. Stretch rates at the stretch and radiation extinction limits respectively decrease and increase with fuel Lewis number before the flammability limit is reached, while the radical Lewis number shows the opposite tendency. In addition, the relation between the standard flammability limit and the limit derived from the strained near stagnation flame is affected by the fuel Lewis number, but not by the radical Lewis number. Meanwhile, the flammability limit increases with decreased fuel Lewis number, but with increased radical Lewis number. Radical behaviours at flame front corresponding to flame bifurcation and extinction are also analysed in this work. It is shown that radical concentration at the flame front, under extinction stretch rate condition, increases with radical Lewis number but decreases with fuel Lewis number. It decreases with increased radiation loss.     

Influences of heat flux on extinction characteristics of steady/unsteady premixed stagnation flames

This paper investigates the extinction characteristics of premixed stagnation flames (PSFs) with controlled heat losses and flow disturbances. The low-frequency air flow pulsations that imitate the operational transients in practical combustors were specially introduced. The tunable diode-laser absorption spectroscopy (TDLAS) measurement was applied to obtain the temperature profile and wall heat flux. It is found that, for steady flame with a fixed equivalence ratio, the extinction stretch rate dramatically increases as the wall heat flux decreases. The extinction criterion is summarized as a global Karlovitz number of 0.57 by establishing a relationship between the global and local stretch rates. Numerical simulations reveal that the local extinction Karlovitz number of steady PSFs is approximately 1.0 regardless of the conditions such as heat flux and equivalence ratio. Further experiments present that the air pulsations with a repetition of ～5 Hz significantly deteriorate the flame stability. Particularly, for unsteady perturbed flames, the extinction stretch rate exhibits a nonlinear trend, yielding two regimes with discrepant sensitivities to wall heat flux. The unsteady simulation then highlights a local stretch rate overshoot in the presence of pulsation. It is caused by the time delay between the inlet velocity and flame front movement that eventually leads to poor flame stability. Moreover, in the high heat-flux regime, a smaller local stretch rate overshoot results in the weak dependence of extinction limits on heat fluxes.     

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.     

Stationary combustion regimes and extinction limits of one-dimensional stretched premixed flames in a gap between two heat conducting plates

Stationary combustion regimes, their linear stability and extinction limits of stretched premixed flames in a narrow gap between two heat conducting plates are studied by means of numerical simulations in the framework of one-dimensional thermal-diffusion model with overall one-step reaction. Various stationary combustion modes including normal flame (NF), near-stagnation plane flame (NSF), weak flame (WF) and distant flame (DF) are detected and found to be analogous to the same-named regimes of conventional counterflow flames. For the flames stabilized in the vicinity of stagnation plane at moderate and large stretch rates (which are NF, NSF and WF) the effect of channel walls is basically reduced to additional heat loss. For distant flame characterized by large flame separation distance and small stretch rates intensive interphase heat transfer and heat recirculation are typical. It is shown that in mixture content / stretch rate plane the extinction limit curve has ε-shape, while for conventional counterflow flames it is known to be C-shaped. This result is quite in line with recent experimental findings and is explained by extension of extinction limits at small stretch rates at the expense of heat recirculation. Analysis of the numerical results makes possible to reveal prime mechanisms of flame quenching on different branches of ε-shaped extinction limit curve. Namely, two upper limits are caused by stretch and heat loss. These limits are direct analogs of the upper and lower limits on conventional C-shaped curve. Two other limits are related with weakening of heat recirculation and heat dissipation to the burner. Thus, the present study provides a satisfactory explanation for the recent experimental observations of stretched flames in narrow channel.     

燃料与自由基的Lewis数对预混气体点火的影响

通过渐近理论分析研究了燃料与自由基的Lewis数对预混气体点火的影响。采用包含自由基的两步化学反应,基于火焰球模型,推导出了描述火焰球半径随点火能以及燃料与自由基的Lewis数而变化的关系式。并在此基础上发现不同参数条件下成功点火的三种情况,研究了燃料与自由基的Lewis数对最小点火能的影响。研究结果表明:随着燃料Lewis数的增大,最小点火能增大;随着自由基Lewis数的增大,最小点火能减小。     

The role of duct thickness on the quenching process of premixed flame propagation

The propagation of premixed laminar flame in ducts of circular cross-section considering a thermal-diffusive model is investigated numerically. Heat losses by conduction to the channels walls are taken into account using the thermally thin wall regime. The effects and the relationship between thickness and diameter of the tube with the flame speed propagation are studied and the quenching condition is obtained as a function of the heat-loss parameter. The mathematical model employs the axisymmetric energy and species equations. The calculations are based on a two-step chemistry, with an Arrhenius, energetically neutral, radical production reaction followed by an exothermic radical recombination reaction. For large values of the heat-loss parameter, the wall temperature is close to the free stream temperature and all the heat losses through the wall are convected away. No heat feedback occurs. On the other hand, for small values of the heat-loss parameter, a feedback mechanism occurs by transferring heat from the burned gas to the fresh mixture along the tube wall. For values of the heat-loss parameter of order unity, the heat feedback mechanism is able to sustain the flame propagation and the quenching condition disappears, producing an almost planar flame front as the propagation velocity reduces. For this two-step reaction mechanism, the radical species behaviour at the duct walls seems to have negligible effect on the quenching process.     

Oscillatory thermal-diffusive instability of combustion waves in a model with chain-branching reaction and heat loss

In this paper we investigate the properties and the linear stability of premixed combustion waves in a non-adiabatic thermal-diffusive model with a two-step chain-branching reaction mechanism. Here we focus only on the emergence of the pulsating instabilities, and the stability analysis is carried out for Lewis numbers for fuel greater than one, and various values of Lewis number for radicals. We consider the problem in two spatial dimensions to allow perturbations of a multidimensional nature. It is demonstrated that the flame speed as a function of the parameters is a double-valued C-shaped function, i.e. for a given set of parameter values there are either two solutions, fast and slow solution branches, propagating with different speed, or the combustion wave does not exist. The extinction of combustion waves occurs at finite values of the parameters and non-zero flame speed. The flame structure demonstrates a slow recombination regime behaviour with negligible fuel leakage for the fast solution branch away from the extinction condition. For parameter values close to the extinction condition and on the slow solution branch, the fuel leakage is significant and a fast recombination regime is observed. It is demonstrated that two types of instabilities emerge in the model: the uniform planar and the travelling instability. The slow solution branch is always unstable due to the uniform perturbations. The fast solution branch is either stable or loses stability due to the travelling or uniform perturbations. The switching between the onset of various regimes of instability is due to the bifurcation of co-dimension two. In the adiabatic limit this bifurcation is found for Lewis number for fuel equal to one, whereas in the non-adiabatic case it moves towards values above unity. The properties of the travelling instability are studied in detail.     

The role of radical wall quenching in flame stability and wall heat flux: hydrogen-air mixtures

The role of wall quenching of radicals in ignition, extinction and autothermal behaviour of premixed H₂–air flames impinging on a flat surface was studied using numerical bifurcation techniques, with detailed gas-phase chemistry and surface radical recombination reactions. Quenching out of radicals was found to retard the system at ignition due solely to the kinetics of the surface reactions. While kinetically extinction is also retarded, the thermal feedback from the wall recombination of radicals can render the flame more stable and lead to a higher wall heat flux as a function of wall temperature compared to an inert surface under some conditions. It is also shown that the combined kinetic and thermal effects of wall radical quenching can expand the autothermal regime. Implications for estimating flammability limits near reactive surfaces of tubes are finally discussed.

M This article features multimedia enhancements available from the abstract page in the online journal; see http://www.iop.org.     

Concept and combustion characteristics of ultra-micro combustors with premixed flame

Under micro-scale combustion influenced by quenching distance, high heat loss, shortened diffusion characteristic time, and flow laminarization, we clarified the most important issues for the combustor of ultra-micro gas turbines (UMGT), such as high space heating rate, low pressure loss, and premixed combustion. The stability behavior of single flames stabilized on top of micro tubes was examined using premixtures of air with hydrogen, methane, and propane to understand the basic combustion behavior of micro premixed flames. When micro tube inner diameters were smaller than 0.4 mm, all of the fuels exhibited critical equivalence ratios in fuel-rich regions, below which no flame formed, and above which the two stability limits of blow-off and extinction appeared at a certain equivalence ratio. The extinction limit for very fuel-rich premixtures was due to heat loss to the surrounding air and the tube. The extinction limit for more diluted fuel-rich premixtures was due to leakage of unburned fuel under the flame base. This clarification and the results of micro flame analysis led to a flat-flame burning method. For hydrogen, a prototype of a flat-flame ultra-micro combustor with a volume of 0.067 cm³ was made and tested. The flame stability region satisfied the optimum operation region of the UMGT with a 16 W output. The temperatures in the combustion chamber were sufficiently high, and the combustion efficiency achieved was more than 99.2%. For methane, the effects on flame stability of an upper wall in the combustion chamber were examined. The results can be explained by the heat loss and flame stretch.     

Structure and dynamic response of radiative diffusion flames

Nitrogen-diluted hydrogen burning in air is modeled numerically using a constant density and one-step reaction model in a plane two-dimensional counterflow configuration. An optically thin assumption is used to investigate the effects of radiation on the dynamics, structure, and extinction of diffusion flames. While there exist dual steady-state extinction limits for the 1D radiative flame response, it is found that as the 1D radiative extinction point is approached the 1D low-stretch diffusion flame exhibits oscillatory response, even with sub-unity Lewis number fuel. These radiation-induced limit cycle oscillations are found to have increasing amplitude and decreasing frequency as the stretch rate is reduced. Flame oscillation eventually leads to permanent extinction at the stretch rate which is larger than the steady-state radiative extinction value. Along the 1D radiative response curve, the transition from 1D flame to 2D structure and the differences in the resulting 2D flame patterns are also examined using a variety of initial profiles, with special emphasis on the comparison of using the initial profiles with and without a flame edge. Similar to the previous studies on the high-stretch adiabatic edge flames using the same configuration, the high-stretch radiative flames are found to resist 1D blow-off quenching through various 2D structures, including propagating front and steady cellular flames for initial profiles with and without flame edges. For all initial profiles studied, the low-stretch radiative flames are also found to exhibit different 2D flame phenomena near the 1D radiative extinction limit, such as transient cellular structures, steady cellular structures, and pulsating ignition fronts. Although the results demonstrate the presence of low-stretch and high-stretch 2D bifurcation branches close to the corresponding 1D extinction limits irrespective of the initial profile used, particular 2D flame structures in certain stretch rate range are initial profile dependent. The existence of two-dimensional flame structures beyond the 1D steady-state radiative extinction limit suggests that the flammable range is expanded as compared to that predicted by the 1D model. Hence, multi-dimensional flame patterns need to be accounted for when determining the flammability limits for a given system.     

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.     

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.     

Interaction of a premixed flame with a liquid fuel film on a wall

In piston engines and in gas turbines, the injection of liquid fuel often leads to the formation of a liquid film on the combustor wall. If a flame reaches this zone, undesired phenomena such as coking may occur and diminish the lifetime of the engine. Moreover, the effect of such an interaction on maximum wall heat fluxes, flame quenching, and pollutant formation is largely unknown. This paper presents a numerical study of the interaction of a premixed flame with a cold wall covered with a film of liquid fuel. Simulations show that the presence of the film leads to a very rich zone at the wall in which the flame cannot propagate. As a result, the flame wall distance remains larger with liquid fuel than it is for a dry wall, and maximum heat fluxes are smaller. The nature of the interaction of flame wall interaction with a liquid fuel is also different from the classical flame/dry wall interaction: it is controlled mainly by chemical mechanisms and not by the thermal quenching effect observed for flames interacting with dry walls: the existence of a very rich zone created above the liquid film is the main mechanism controlling quenching.     

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.     

Oscillating edge-flames

It has been known for some years that when a near-limit flame spreads over a liquid pool of fuel, the edge of the flame can oscillate. It is also known that when a near-asphyxiated candle-flame burns in zero gravity, the edge of the (hemispherical) flame can oscillate violently prior to extinction. We propose that these oscillations are nothing more than a manifestation of the large Lewis number instability well known in chemical reactor studies and in combustion studies, one that is exacerbated by heat losses. As evidence of this we examine an edge-flame confined within a fuel-supply boundary and an oxygen-supply boundary, anchored by a discontinuity in data at the fuel-supply boundary. We show that when the Lewis number of the fuel is 2, and the Lewis number of the oxidizer is 1, oscillations of the edge occur when the Damköhler number is reduced below a critical value. During a single oscillation period there is a short premixed propagation stage and a long diffusion stage, behaviour that has been observed in flame spread experiments. Oscillations do not occur when both Lewis numbers are equal to 1.     

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.     

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.     

Flame structure and flame spread rate over a solid fuel in partially premixed atmospheres

We have investigated the downward flame spread over a thin solid fuel. Hydrogen, methane, or propane, included in the gaseous product of pyrolysis reaction, is added in the ambient air. The fuel concentration is kept below the lean flammability limit to observe the partially premixing effect. Both experimental and numerical studies have been conducted. Results show that, in partially premixed atmospheres, both blue flame and luminous flame regions are enlarged, and the flame spread rate is increased. Based on the flame index, a so-called triple flame is observed. The heat release rate ahead of the original diffusion flame is increased by adding the fuel, and its profile is moved upstream. Here, we focus on the heat input by adding the fuel in the opposed air, which could be a direct factor to intensify the combustion reaction. The dependence of the flame spread rate on the heat input is almost the same for methane and propane/air mixtures, but larger effect is observed for hydrogen/air mixture. Since the deficient reactant in lean mixture is fuel, the larger effect of hydrogen could be explained based on the Lewis number consideration. That is, the combustion is surely intensified for all cases, but this effect is larger for lean hydrogen/air mixture (Le < 1), because more fuel diffuses toward the lean premixed flame ahead of the original diffusion flame. Resultantly, the pyrolysis reaction is promoted to support the higher flame spread rate.