Effects of heat conduction and radical quenching on premixed stagnation flame stabilised by a wall

The premixed stagnation flame stabilised by a wall is analysed theoretically considering thermally sensitive intermediate kinetics. We consider the limit case of infinitely large activation energy of the chain-branching reaction, in which the radical is produced infinitely fast once the cross-over temperature is reached. Under the assumptions of potential flow field and constant density, the correlation for flame position and stretch rate of the premixed stagnation flame is derived. Based on this correlation, the effects of heat conduction and radical quenching on the wall surface are examined. The wall temperature is shown to have great impact on flame bifurcation and extinction, especially when the flame is close to the wall. Different flame structures are observed for near-wall normal flame, weak flame, and critically quenched flame. The fuel and radical Lewis numbers are found to have opposite effects on the extinction stretch rate. Moreover, it is also demonstrated that only when the flame is close to the wall does the radical quenching strongly influence the flame bifurcation and extinction. The extinction stretch rate is shown to decrease with the amount of radical quenching for different fuel and radical Lewis numbers. Besides, the coupling between the wall heat conduction and radical quenching is found to greatly influence the bifurcation and extinction of the premixed stagnation flame.     

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.     

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.     

Opposed flow flame spread over an array of thin solid fuel sheets in a microgravity environment

In this work a numerical study has been carried out to gain physical insight into the phenomena of opposed flow flame spread over an array of thin solid fuel sheets in a microgravity environment. The two-dimensional (2D) simulations show that the flame spread rates for the multiple-fuel configuration are higher than those for the flame spreading over a single fuel sheet. This is due to reduced radiation losses from the flame and increased heat feedback to the solid fuel. The flame spread rate exhibits a non-monotonic variation with decrease in the interspace distance between the fuel sheets. Higher radiation heat feedback primarily as gas/flame radiation was found to be responsible for the increase in the flame spread rate with the reduction of the interspace distance. It was noted that as the interspace distance between the fuel sheets was reduced below a certain value, no steady solution could be obtained. However, at very small interspace distances, steady state spread rates were obtained. Here, due to oxygen starvation the flame spread rate decreased and eventually at some interspace distance the flame extinguished. With fuel emittance (equal to absorptance) reduced to ‘0’ the flame spread rate was nearly independent of the interspace distance, except at very small distances where the flame spread rate dropped due to oxygen starvation. A flame extinction plot with the extinction oxygen level was constructed for the multiple-fuel configuration at various interspace distances. The default fuel with an emittance of 0.92 was found to be more flammable in the multiple-fuel configuration than in a single fuel sheet configuration. For a fuel emittance equal to zero, the extinction oxygen limit decreases for both the single and the multiple fuel sheet configurations. However, the two flammability curves cross over at a certain fuel separation distance. The multiple-fuel configurations become less flammable compared to the single fuel sheet configuration below a certain separation distance.     

Effects of endothermic chain-branching reaction on spherical flame initiation and propagation

A theoretical model is developed to describe the spherical flame initiation and propagation. It considers endothermic chain-branching reaction and exothermic recombination reaction. Based on this model, the effects of endothermic chain-branching reaction on spherical flame initiation and propagation are assessed. First, the analytical solutions for the distributions of fuel and radical mass fraction as well as temperature are obtained within the framework of large activation energy and quasi-steady assumption. Then, a correlation describing spherical flame initiation and propagation is derived. Based on this correlation, different factors affecting spherical flame propagation and initiation are examined. It is found that endothermicity of the chain-branching reaction suppresses radical accumulation at the flame front and thus reduces flame intensity. With the increase of endothermicity, the unstretched flame speed decreases while both flame ball radius and Markstein length increases. Endothermicity has a stronger effect on the stretched flame speed with larger fuel Lewis number. The Markstein length is found to increase monotonically with endothermicity. Furthermore, the endothermicity of the chain-branching reaction is shown to affect the transition among different flame regimes including ignition kernel, flame ball, propagating spherical flame, and planar flame. The critical ignition power radius increases with endothermicity, indicating that endothermicity inhibits the ignition process. The influence of endothermicity on ignition becomes relatively stronger at higher crossover temperature or higher fuel Lewis number. Moreover, one-dimensional transient simulations are conducted to validate the theoretical results. It is shown that the quasi-steady-state assumption used in theoretical analysis is reasonable and that the same conclusion on the effects of endothermic chain-branching reaction can be drawn from simulation and theoretical analysis.     

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.     

On the global structure and asymptotic stability of low-stretch diffusion flame: Forced convection

The present work analyzes cylindrical diffusion flames (Tsuji burner) under low stretch condition, considering fuel injection also from the backward region of the burner. To highlight the fundamental aspects of this flame, some assumptions are imposed, like constant thermodynamic and transport coefficients, unitary Lewis number and no radiative heat loss. It is also considered potential flow model and incompressible Navier–Stokes model. Despite the simplicity of the former model, results from both models show good agreement. Also, an asymptotic analysis describing the problem far from the burner is able to capture the most important mechanisms controlling the flame, then the flame shape is determined and the dependence of the characteristic length scales on Peclet number (based on the burner properties), free stream velocity and stoichiometry is revealed. The results show that the flame width is proportional to the mass stoichiometric coefficient and reciprocal to the Peclet number the 1/4 power and free stream velocity the 3/4 power, and that the flame height is proportional to the square of the mass stoichiometric coefficient and to the square root of the ratio of Peclet number to free stream velocity. In addition, an asymptotic stability analysis reveals low-stretch flame extinction to be caused by reduction in fuel and oxidizer concentrations, which provides the range of the stoichiometric coefficient for stable regime, and at the same time the range of heat released.     

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.     

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.     

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.     

Edge flames stabilized in a non-premixed microcombustor

The dynamics of an edge flame confined in a non-premixed microcombustor model is studied numerically within the context of a diffusive-thermal model. Fuel and oxidizer, separated upstream by a thin plate, flow through a channel with a prescribed velocity. At the tip of the plate, the fuel and oxidizer mix and, when ignited, an edge flame is sustained at some distance from the plate. The objective in this work is to consider the effects of confinement, differential diffusion, and heat loss on the dynamics of an edge flame in a narrow channel. We consider a wide range of channel widths and allow for changing Lewis numbers, and both adiabatic conditions and heat losses along the channel walls. The results illustrate how the flame shape and standoff distance are affected by the channel width, by mixture composition through variations in Lewis numbers and by heat losses. Conditions for flame stabilization, flame oscillations and flame extinction or blowoff are predicted.     

Dynamics and burning limits of near-limit hot,warm, and cool diffusion flames of dimethyl ether at elevated pressures

The near-limit diffusion flame regimes and extinction limits of dimethyl ether at elevated pressures and temperatures are examined numerically in the counterflow geometry with and without radiation at different oxygen concentrations. It is found that there are three different flame regimes—hot flame, warm flame, and cool flame—which exist, respectively, at high, intermediate, and low temperatures. Furthermore, they are governed by three distinct chain-branching reaction pathways. The results demonstrate that the warm flame has a double reaction zone structure and plays a critical role in the transition between cool and hot flames. It is also shown that the cool flame can be formed in several different ways: by either radiative extinction or stretch extinction of a hot flame or by stretch extinction of a warm flame. A warm flame can also be formed by radiative extinction of a hot flame or ignition of a cool flame. A general €-shaped flammability diagram showing the burning limits of all three flame regimes at different oxygen mole fractions is obtained. The results show that thermal radiation, reactant concentration, temperature, and pressure all have significant impacts on the flammable regions of the three flame regimes. Increases in oxidizer temperature, oxygen concentration, and pressure shift the cool flame regime to higher stretch rates and cause the warm flame to have two extinction limits. At elevated temperatures, it is found that there is a direct transition between the hot flame and warm flame at low stretch rates. The results also show that, unlike the hot flame, the cool flame structure cannot be scaled by using pressure-weighted stretch rates due to the its significant reactant leakage and strong dependence of reactivity on pressure. The present results advance the understanding of near-limit flame dynamics and provide guidance for experimental observation of different flame regimes.     

Propagation and extinction of an unsteady spherical spray flame front

Experimental evidence seems to indicate that the life of a laminar spherical flame front propagating through a fresh mixture of air and liquid fuel droplets can be roughly split into three stages: (1) ignition, (2) radial propagation with a smooth flame front and (3) propagation with flame front cellularization and/or pulsation. In this work, the second stage is analysed using the slowly varying flame approach, for a fuel rich flame. The droplets are presumed to vaporize in a sharp front ahead of the reaction front. Evolution equations for the flame and evaporation fronts are derived. For the former the combined effect of heat loss due to droplet vaporization and radiation plays a dominant explicit role. In addition, the structure of the evaporation front is deduced using asymptotics based on a large parameter associated with spray vaporization. Numerical calculations based on the analysis point to the way in which the spray modifies conditions for flame front extinction. Within the framework of the present simplified model the main relevant parameters turn out to be the initial liquid fuel load in the fresh mixture and/or the latent heat of vaporization of the fuel.     

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.     

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.     

Experimental low-stretch gaseous diffusion flames in buoyancy-induced flowfields

The present study experimentally investigates the structure and instabilities associated with extremely low-stretch (1 s⁻¹) gaseous diffusion flames. Ultra-low-stretch flames are established in normal gravity by bottom burning of a methane/nitrogen mixture discharged from a porous spherically symmetric burner of large radius of curvature. OH-PLIF and IR imaging techniques are used to characterize the reaction zone and the burner surface temperature, respectively. A flame stability diagram mapping the response of the ultra-low-stretch diffusion flame to varying fuel injection rate and nitrogen dilution is explored. In this diagram, two main boundaries are identified. These boundaries separate the stability diagram into three regions: sooting flame, non-sooting flame, and extinction. Two distinct extinction mechanisms are noted. For low fuel injection rates, flame extinction is caused by heat loss to the burner surface. For relatively high injection rates, at which the heat loss to burner surface is negligible, flame radiative heat loss is the dominant extinction mechanism. There also exists a critical inert dilution level beyond which the flame cannot be sustained. The existence of multi-dimensional flame phenomena near the extinction limits is also identified. Various multi-dimensional flame patterns are observed, and their evolutions are studied using direct chemiluminescence and OH-PLIF imaging. The results demonstrate the usefulness of the present burner configuration for the study of low-stretch gaseous diffusion flames.     

Cellular and sporadic flame regimes of low-Lewis-number stretched premixed flames

3D structure and dynamical behavior of low-Lewis-number stretched premixed flames are numerically simulated within the framework of a thermo-diffusive model with one-step chemical reaction. The results are compared with microgravity experiments at qualitative level. The influence of Lewis number, equivalence ratio, and heat loss intensity on flame structure is investigated. It is experimentally and numerically found that lean counterflow flames can appear as a set of separate ball-like flames in a state of chaotic motion. It is shown that the time averaged flame balls coordinate may be considered as important characteristic similar to coordinate of continuous flame front. Numerical simulations reveal essential incompleteness of combustion at high level of heat losses. This incompleteness occurs in the process of lean mixtures combustion and is caused by fuel leakage through the gaps among ball-like flames.     

Triple flame: Inherent asymmetries and pentasectional character

A two-dimensional triple-flame numerical model of a laminar combustion process reflects flame asymmetric structural features that other analytical models do not generate. It reveals the pentasectional character of the triple flame, composed of the central pure diffusion-flame branch and the fuel-rich and fuel-lean branches, each of which is divided into two sections: a near-stoichiometric section and a previously unreported near-flammability-limits section with combined diffusion and premixed character. Results include propagation velocity, fuel and oxidiser mass fractions, temperature and reaction rates. Realistic stoichiometric ratios and reaction orders match experimental planar flame characteristics. Constant density, a one-step reaction, and a mixture fraction gradient at the inlet as the simulation parameter are imposed. The upstream equivalence ratio or the upstream reactant mass fractions are linear or hyperbolic functions of the transverse coordinate. The use here of experimental kinetics data differs from previous analytical works and results in flame asymmetry and different flammability limits. Upstream mixture composition gradient affects propagation velocity, flame curvature, diffusion flame reaction rate, and flammability limits. Flammability limits extend beyond those of a planar flame due to transverse heat and mass diffusion causing the pentasectional character.     

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.