Flame fingers and interactions of hydrodynamic and thermodiffusive instabilities in laminar lean hydrogen flames

Lean hydrogen/air flames are prone to hydrodynamic and thermodiffusive instabilities. In this work, the contribution of each instability mechanism is quantified separately by performing detailed simulations of laminar planar lean hydrogen/air flames with different diffusivity models and equations of state to selectively suppress the hydrodynamic or thermodiffusive instability mechanism.From the analysis of the initial phase of the simulations, the thermodiffusive instability is shown to dominate the flame dynamics. If differential diffusion and, hence, the thermodiffusive instability is suppressed, the flame features a strong reduction of the instability growth rates, whereas if present, a wide range of unstable wave numbers is observed due to the strong destabilizing nature of differential diffusion. When instabilities are fully developed, lean hydrogen/air flames feature the formation of small-scale cellular structures and large-scale flame fingers. While the size of the former is known to be close to the most unstable wave length of a linear stability analysis, this work shows that flame fingers also originate from the thermodiffusive instability and most noteworthy, are not linked to an interaction of the two instability mechanisms. They are stable with respect to external perturbations and feature an enhanced flame propagation as the formation of a central cusp at their tip enables the co-existence of two strongly curved leading edges with high reactivity. The thermodiffusive instability is shown to significantly affect the flames’ consumption speed, while the consumption speed enhancement caused by the hydrodynamic instability is significantly smaller. Further, the surface area increase due to wrinkling is strongly diminished if one of the two instability mechanisms is missing. This is linked to a synergistic interaction between the two mechanisms, as the propagation of flame fingers is enhanced by the presence of the hydrodynamic instability due to a widening of the streamlines ahead of the flame fingers.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

Nonlinear equation for curved nonstationary flames and flame stability

A time-dependent nonlinear equation for a nonstationary curved flame front of an arbitrary expansion coefficient is derived under the assumptions of a small but finite flame thickness and weak nonlinearity. On the basis of the derived equation, stability of two-dimensional curved stationary flames propagating in tubes with ideally adiabatic and slip walls is studied. The stability analysis shows that curved stationary flames become unstable for sufficiently wide tubes. The obtained stability limits are in a good agreement with the results of numerical simulations of flame dynamics and with semiqualitative stability analysis of curved stationary flames. Possible outcomes of the obtained instability at the nonlinear stage are discussed. The instability may result in extra wrinkles at a flame front close to the stability limits and in self-turbulization of the flame far from the limits. The self-turbulization can also be interpreted as a fractal structure. The fractal dimension of a flame front and velocity of a self-turbulized flame are evaluated.     

Computational and experimental study of oxygen-enhanced axisymmetric laminar methane flames

Three axisymmetric laminar coflow diffusion flames, one of which is a nitrogen-diluted methane/air flame (the ‘base case’) and the other two of which consist of nitrogen-diluted methane vs. pure oxygen, are examined both computationally and experimentally. Computationally, the local rectangular refinement method is used to solve the fully coupled nonlinear conservation equations on solution-adaptive grids. The model includes C2 chemistry (GRI 2.11 and GRI 3.0 chemical mechanisms), detailed transport, and optically thin radiation. Because two of the flames are attached to the burner, thermal boundary conditions at the burner surface are constructed from smoothed functional fits to temperature measurements. Experimentally, Raman scattering is used to measure temperature and major species concentrations as functions of the radial coordinate at various axial positions. As compared to the base case flame, which is lifted, the two oxygen-enhanced flames are shorter, hotter, and attached to the burner. Computational and experimental flame lengths show excellent agreement, as do the maximum centreline temperatures. For each flame, radial profiles of temperature and major species also show excellent agreement between computations and experiments, when plotted at fixed values of a dimensionless axial coordinate. Computational results indicate peak NO levels in the oxygen-enhanced flames to be very high. The majority of the NO in these flames is shown to be produced via the thermal route, whereas prompt NO dominates for the base case flame.     

Critical conditions at liftoff limit of a laminar n-butane-air jet diffusion flame

The stability mechanism of laminar coflow jet diffusion flames in normal gravity has been studied computationally and experimentally. N-butane, the heaviest alkane in a gaseous state at ambient temperature and pressure, is used as the fuel since the reaction mechanism is similar to that of higher (liquid) hydrocarbons. The critical mean n-butane jet and coflowing air velocities at flame stability limits are measured using a small fuel tube burner (0.8 mm inner diameter). The time-dependent, axisymmetric numerical code with a detailed reaction mechanism (58 species and 540 reactions), molecular diffusive transport, and a radiation model, reveals a flame structure. A fuel-lean peak reactivity spot (i.e., reaction kernel), possessing the hybrid nature of diffusion-premixed flame structure at a constant temperature of ≈1560 K, is formed at the flame base and controls the flame stability. In a near-quiescent environment, the flame base resides below the fuel tube exit plane and thereby premixing is limited. As the coflowing air velocity is increased incrementally under a fixed fuel jet velocity, the flame base moves slightly above (≈1 mm) the burner exit and vigorous premixed combustion becomes prevailing. The local heat-release rate at the reaction kernel nearly doubles due to the increased convective oxygen flux (i.e., a blowing effect). The local Damköhler number, newly defined as a ratio of the square root of the local heat-release rate and the local velocity, decreases gradually first and drops abruptly at a critical threshold value and the flame base lifts off from the burner rim. The calculated coflow air velocity at liftoff is ≈0.38 m/s at the fuel jet velocity of 2 m/s, which is consistent with an extrapolated measured value of 0.41 m/s. This work has determined the critical Damköhler number at the stability limit quantitatively, for the first time, for laminar jet diffusion flames.     

Large eddy simulation of piloting effects on turbulent swirling flames

Pilot flames, created by additional injectors of pure fuel, are often used in turbulent burners to enhance flame stabilization and reduce combustion instabilities. The exact mechanisms through which these additional rich zones modify the flame anchoring location and the combustion dynamics are often difficult to identify, especially when they include unsteady hydrodynamic motion. This study presents Large Eddy Simulations (LES) of the reacting flow within a large-scale gas turbine burner for two different cases of piloting, where either 2 or 6% of the total methane used in the burner is injected through additional pilot flame lines. For each case, LES shows how the pilot fuel injection affects both flame stabilization and flame stability. The 6% case leads to a stable flame and limited hydrodynamic perturbations in the initial flame zone. The 2% case is less stable, with a small-lift-off of the flame and a Precessing Vortex Core (PVC) in the cold stabilization zone. This PVC traps some of the lean cold gases issuing from the pilot passage stream, changes the flame stabilization point and induces instability.     

Numerical studies of curved stationary flames in wide tubes

The nonlinear problem of the propagation of curved stationary flames in tubes of different widths is studied by means of direct numerical simulation of the complete system of hydrodynamic equations including thermal conduction, viscosity, fuel diffusion and chemical kinetics. While only a planar flame can propagate in a narrow tube of width smaller than half of the cut–off wavelength determined by the linear theory of the hydrodynamic instability of a flame front, in wider tubes stationary curved flames propagate with velocities considerably larger than the corresponding velocity of a planar flame. It is shown that only simple ‘single-hump’ slanted stationary flames are possible in wide tubes, and ‘multi–hump’ flames are possible in wide tubes only as a nonstationary mode of flame propagation. The stability limits of curved stationary flames in wider tubes and the secondary Landau–Darrieus instability are investigated. The dependence of the velocity of the stationary flame on the tube width is studied. The analytical theory describes the flame reasonably well when the tube width does not exceed some critical value. The dynamics of the flame in wider tubes is shown to be governed by a large–scale stability mechanism resulting in a highly slanted flame front. In wide tubes, the skirt of the slanted flame remains smooth with the length of the skirt and the flame velocity increasing progressively with the increase of the tube width above the second critical value. Results of the analytical theory and numerical simulations are discussed and compared with the experimental data for laminar flames in wide tubes.     

Turbulent propagation of premixed flames in the presence of Darrieus–Landau instability

We investigate the role played by hydrodynamic instability in the wrinkled flamelet regime of turbulent combustion, where the intensity of turbulence is small compared to the laminar flame speed and the scale large compared to the flame thickness. To this end the Michelson–Sivashinsky (MS) equation for flame front propagation in one and two spatial dimensions is studied in the presence of uncorrelated and correlated noise representing a turbulent flow field. The combined effect of turbulence intensity, integral scale, and an instability parameter related to the Markstein length are examined and turbulent propagation speed monitored for both stable planar flames and corrugated flames for which the planar conformation is unstable. For planar flames a particularly simple scaling law emerges, involving quadratic dependence on intensity and a linear dependence on the degree of instability. For corrugated flames we find the dependence on intensity to be substantially weaker than quadratic, revealing that corrugated flames are more resilient to turbulence than planar flames. The existence of a threshold turbulence intensity is also observed, below which the corrugated flame in the presence of turbulence behaves like a laminar flame. We also analyze the conformation of the flame surface in the presence of turbulence, revealing primary, large-scale wrinkles of a size comparable to the main corrugation. When the integral scale is much smaller than the characteristic corrugation length we observe, in addition to primary wrinkles, secondary small-scale wrinkles contaminating the surface. The flame then acquires a multi-scale, self-similar conformation, with a fractal dimension, for one-dimensional flames, plateauing at 1.23 for large intensities. The existence of an intermediate integral scale is also found at which the turbulent speed is maximized. When two-dimensional flames are subject to turbulence, the primary wrinkling patterns give rise to polyhedral-cellular structures which bear a very close resemblance to those observed in experiments on hydrodynamically unstable expanding spherical flames.     

Effect of thermal expansion on the linear stability of planar premixed flames for a simple chain-branching model: The high activation energy asymptotic limit

The linear stability of freely propagating, adiabatic, planar premixed flames is investigated in the context of a simple chain-branching chemistry model consisting of a chain-branching reaction step and a completion reaction step. The role of chain-branching is governed by a crossover temperature. Hydrodynamic effects, induced by thermal expansion, are taken into account and the results compared and contrasted with those from a previous purely thermal-diffusive constant density linear stability study. It is shown that when thermal expansion is properly accounted for, a region of stable flames predicted by the constant density model disappears, and instead the flame is unstable to a long-wavelength cellular instability. For a pulsating mode, however, thermal expansion is shown to have only a weak effect on the critical fuel Lewis number required for instability. These effects of thermal expansion on the two-step chain-branching flame are shown to be qualitatively similar to those on the standard one-step reaction model. Indeed, as found by constant density studies, in the limit that the chain-branching crossover temperature tends to the adiabatic flame temperature, the two-step model can be described to leading order by the one-step model with a suitably defined effective activation energy.     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Stabilization, propagation and instability of tribrachial triple flames

A tribrachial (or triple) flame is one kind of edge flame that can be encountered in nonpremixed mixing layers, consisting of a lean and a rich premixed flame wing together with a trailing diffusion flame all extending from a single point. The flame could play an important role on the characteristics of various flame behaviors including lifted flames in jets, flame propagation in two-dimensional mixing layers, and autoignition fronts. The structure of tribrachial flame suggests that the edge is located along the stoichiometric contour in a mixing layer due to the coexistence of all three different types of flames. Since the edge has a premixed nature, it has unique propagation characteristics. In this review, the propagation speed of tribrachial flames will be discussed for flames propagating in mixing layers, including the effects of concentration gradient, velocity gradient, and burnt gas expansion. Based on the tribrachial edge structure observed experimentally in laminar lifted flames in jets, the flame stabilization characteristics including liftoff height, reattachment, and blowout behaviors and their buoyancy-induced instability will be explained. Various effects on liftoff heights in both free and coflow jets including jet velocity, the Schmidt number of fuel, nozzle diameter, partial premixing of air to fuel, and inert dilution to fuel are discussed. Implications of edge flames in the modeling of turbulent nonpremixed flames and the stabilization of turbulent lifted flames in jets are covered.     

The unstable behavior of cellular premixed flames induced by intrinsic instability

The unstable behavior of cellular premixed flames induced by intrinsic instability is studied by two-dimensional unsteady calculations of reactive flows. In the present numerical simulation, the compressible Navier–Stokes equation including a one-step irreversible chemical reaction is employed. We consider two basic types of phenomena to account for the intrinsic instability of premixed flames, i.e., hydrodynamic and diffusive-thermal effects. The hydrodynamic effect is caused by the thermal expansion through the flame front; the diffusive-thermal effect is caused by the preferential diffusion of mass versus heat. A disturbance with several wavelength components is superimposed on a planar flame, and the formation of a cellular flame induced by hydrodynamic and diffusive-thermal effects is numerically simulated. After the cellular-flame formation, the combination and division of cells are observed. The behavior of cellular-flame fronts becomes more unstable when the Lewis number is lower than unity, since the diffusive-thermal effect has a great influence on the unstable behavior. The cell size changes with time, and its average is greater than the critical wavelength and becomes smaller by decreasing the Lewis number. The flame velocity of cellular flames depends strongly on the length of computational domain in the direction tangential to the flame front. As the length of computational domain increases, the flame velocity becomes larger. This is because the long-wavelength components of disturbances play an important role in the shape of cellular flames, i.e., in the flame-surface area.     

Low pressure premixed CH4/air flames with forced periodic mixture fraction oscillations: experimental approach

An experimental setup for the generation and investigation of periodic equivalence ratio oscillations in laminar premixed flames is presented. A special low-pressure burner was developed which generates stable flames in a wide pressure range down to 20 mbar and provides the possibility of rapid mixture fraction variations. The technical realization of the mixture fraction variations and the characteristics of the burner are described. 1D laser Raman scattering was applied to determine the temperature and concentration profiles of the major species through the flame front in correlation to the phase-angle of the periodic oscillation. OH* chemiluminescence was detected to qualitatively analyze the response of the flame to mixture fraction variations by changing shape and position. Exemplary results from a flame at p=69 mbar, forced at a frequency of 10 Hz, are shown and discussed. The experiments are part of a cooperative research project including the development of kinetic models and numerical simulation tools with the aim of a better understanding and prediction of periodic combustion instabilities in gas turbines. The focus of the current paper lies on the presentation of the experimental realization and the measuring techniques.     

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.     

Numerical studies of flames in wide tubes: stability limits of curved stationary flames

Flame dynamics in wide tubes with ideally adiabatical and slip walls is studied by means of direct numerical simulations of the complete set of hydrodynamical equations including thermal conduction, fuel diffusion, viscosity, and chemical kinetics. Stability limits of curved stationary flames in wide tubes and the hydrodynamic instability of these flames (the secondary Darrieus-Landau instability) are investigated. The stability limits found in the present numerical simulations are in a very good agreement with the previous theoretical predictions. It is obtained that close to the stability limits the secondary Darrieus-Landau instability results in an extra cusp at the flame front. It is shown that the curved flames subject to the secondary Darrieus-Landau instability propagate with velocity considerably larger than the velocity of the stationary flames.     

Stabilization mechanism of lifted flame edge in the near field of coflow jets for diluted methane

The stabilization mechanism of lifted flames in the near field of coflow jets has been investigated experimentally and numerically for methane fuel diluted with nitrogen. The lifted flames were observed only in the near field of coflow jets until blowout occurred in the normal gravity condition. To elucidate the stabilization mechanism for the stationary lifted flames of methane having the Schmidt number smaller than unity, the behavior of the flame in the buoyancy-free condition, and unsteady propagation characteristics after ignition were investigated numerically at various conditions of jet velocity. It has been found that buoyancy plays an important role for flame stabilization of lifted flames under normal gravity, such that the flame becomes attached to the nozzle in microgravity. The stabilization mechanism is found to be due to the variation of the propagation speed of the lifted flame edge with axial distance from the nozzle in the near field of the coflow as compared to the local flow velocity variation at the edge.     

Characteristics of microjet methane diffusion flames

Characteristics of microjet methane diffusion flames stabilized on top of the vertically oriented, stainless-steel tubes with an inner diameter ranging from 186 to 778 μ m are investigated experimentally, theoretically and numerically. Of particular interest are the flame shape, flame length and quenching limit, as they may be related to the minimum size and power of the devices in which such flames would be used for future micro-power generation. Experimental measurements of the flame shape, flame length and quenching velocity are compared with theoretical predictions as well as detailed numerical simulations. Comparisons of the theoretical predictions with measured results show that only Roper's model can satisfactorily predict the flame height and quenching velocity of microjet methane flames. Detailed numerical simulations, using skeletal chemical kinetic mechanism, of the flames stabilized at the tip of d = 186, 324 and 529 μ m tubes are performed to investigate the flame structures and the effects of burner materials on the standoff distance near extinction limit. The computed flame shape and flame length for the d = 186 μm flame are in excellent agreement with experimental results. Numerical predictions of the flame structures strongly suggest that the flame burns in a diffusion mode near the extinction limit. The calculated OH mass fraction isopleths indicate that different tube materials have a minor effect on the standoff distance, but influence the quenching gap between the flame and the tube.     

Stability characteristics and flowfields of turbulent non-premixed swirling flames

A simple, yet representative, burner geometry is used for the investigation of highly swirling turbulent unconfined, non-premixed, flames of natural gas. The burner configuration comprises a ceramic faced bluff-body with a central fuel jet. The bluff-body is surrounded by an annulus that delivers a swirling primary flow of air. The entire burner assembly is housed in a wind tunnel providing a secondary co-flowing stream of air. This hybrid bluff-body/swirl burner configuration stabilizes complex turbulent flames not unlike those found in practical combustors, yet is amenable to modelling because of its well-defined boundary conditions. Full stability characteristics including blow-off limits and comprehensive maps of flame shapes are presented for swirling flames of three different fuel mixtures: compressed natural gas (CNG), CNG–air (1:2 by volume) and CNG–H₂ (1:1 by volume).

It is found that with increased fuel flow, flame blow-off mode may change with swirl number, S_g. At low swirl, the flame remains stable at the base but blows off in the neck region further downstream. At higher swirl numbers, the flames peel off completely from the burner's base. Swirling CNG–air flames are distinct in that they only undergo base blow-off. In the low range of swirl number, increasing S_g causes limited improvement in the blow-off limits of the flames investigated and (for a few cases) can even lead to some deterioration over a small intermediate range of S_g. It is only above a certain threshold of swirl that significant improvements in blow-off limits appear. Six flames are selected for further detailed flowfield and composition measurements and these differ in the combination of swirl number, primary axial velocity through the annulus, U_s, and bulk fuel jet velocity, U_j. Only velocity field measurements are presented in this paper. A number of flow features are resolved in these flames, which resemble those already associated with non-reacting swirling flows of equivalent swirl obtained with the present burner configuration. Additionally, asymmetric flowfields inherent to some flames are revealed where the fluidic centreline of the flow (defined in the two-dimensional (U–W velocity pair) velocity field by the ?ω? = 0 tangential velocity contour), meanders strongly on either side of the geometric centreline downstream by about one bluff-body diameter. Flow structures revealed by the velocity data are correlated to flame shapes to yield a better understanding of how the velocity field influences the flames physical characteristics.     

Numerical simulation of flames as gas-dynamic discontinuities

The dynamics of thin premixed flames is computationally studied within the context of a hydrodynamic theory. A level-set method is used to track down the flame, which is treated as a free-boundary interface. The flow field is described by the incompressible Navier–Stokes equations, with different densities for the burnt and unburnt gases, supplemented by singular source terms that properly account for thermal expansion effects. The numerical scheme has been tested on several benchmark problems and was shown to be stable and accurate. In particular, the propagation of a planar flame front and the dynamics of hydrodynamically unstable flames were successfully simulated. This includes recovering the planar front in narrow domains, the Darrieus–Landau linear growth rate for long waves of small amplitude, and the nonlinear development of cusp-like structures predicted by the Michelson–Sivashinsky equation for a small density change. The stationary flame of a Bunsen burner with uniform and parabolic outlet flows were also simulated, showing in particular a careful mapping of the flow field. Finally, the evolution of a hydrodynamically unstable flame was studied for finite amplitude disturbances and realistic values of thermal expansion. These results, which constitute one of the main objectives of this study, elucidate the effect of thermal expansion on flame dynamics.