An experimental realization of an unstrained, planar diffusion flame

A unique burner was constructed to experimentally realize a one-dimensional unstrained planar non-premixed flame, previously considered only in idealized theoretical models. One reactant, the fuel mixture in the current experiments, is supplied through a porous plug at the bottom of the combustion chamber and flows vertically up towards the horizontal flame. The crux of the design is the introduction of the oxidizer from above in such a way that its diffusion against the upward product flow is essentially one-dimensional, i.e., uniform over the burner cross-section. This feature was implemented by introducing the oxidizer into the burner chamber from the top through an array of 625 closely spaced hypodermic needles, and allowing the hot products to escape vertically up through the space between the needles. Due to the injection of oxidizer through discrete tubes, a three-dimensional “injection layer” exists below the exit plane of the oxidizer supply tubes. Experimental evidence suggests that this layer is thin and that oxidizer is supplied to the flame by 1-D counterdiffusion, producing a nearly unstrained flame. To characterize the burner, flame position measurements were conducted for different compositions and flowrates of H₂–CO₂ and O₂–CO₂ mixtures. The measured flame locations are compared to an idealized one-dimensional model in which only diffusion of oxidizer against the product flow is considered. The potential of the new burner is demonstrated by a study of cellular structures forming near the extinction limit. Consistent with previous investigations, cellular instabilities are shown to become more prevalent as the initial mixture strength and/or the Damköhler number are decreased. As the extinction limit is approached, the number of cells was observed to decrease progressively.     

Cell formation in non-premixed,axisymmetric jet flames near extinction

Systematic experiments with CO₂ diluted H₂–O₂ circular jet diffusion flames have been undertaken to study the formation of cellular flames, which occur for relatively low reactant Lewis numbers and near the extinction limit. The jet Reynolds number for all experiments was about 500, based on the centreline velocity, jet diameter and ambient fuel properties. The Lewis numbers, based on the initial mixture strength φ and ambient conditions of the investigated near-extinction mixtures, vary in the range 1.1–1.3 for oxygen and 0.25–0.29 for hydrogen (φ is defined here as the fuel-to-oxygen molar ratio normalized by the stoichiometric value). Various conditions near the extinction limit were investigated by fixing the fuel composition (H₂–CO₂ mixture), and systematically lowering the oxygen concentration in the co-flowing oxidizer stream past the point where cellular structures formed, until extinction occurred. The observed different instability states were correlated with the initial mixture strength and the proximity to the extinction limit.

The parameter space for cellularity was found to increase with decreasing initial mixture strength. For a given initial mixture strength, several cellular states were found to co-exist near the extinction limit, and the preferred number of cells (the azimuthal wave number) was observed to decrease with decreasing oxygen concentration (Damköhler number). These trends are consistent with previous theoretical work and our own stability analysis that will be reported elsewhere.     

Edge-flames and cellular structures in oscillatory premixed counterflows

We examine edge-flames in the context of symmetric counterflows of fresh mixture, and examine their dynamics when the rate of strain is varied periodically in time. For a Lewis number of 1, extinction boundaries and zero (in the mean) edge propagation speed boundaries are constructed in the forcing amplitude—Damköhler number plane. Forcing can turn advancing edges into retreating edges. For a Lewis number of 0.3, the fundamental distinction is between flames that display cellular structures and flames that do not. Forcing can convert a non-cellular flame into one with cells; and it can strongly affect the dynamics of cellular structures that exist in the absence of forcing.     

On intrinsic oscillation in radiation-affected diffusion flames

A linear stability analysis is conducted to study the onset of near-limit flame oscillation with radiative heat loss in 1-D chambered planar flames using multi-scale activation-energy asymptotics. The oscillatory instability near the radiation-induced extinction limit at large Damköhler numbers is identified, in additional to the one near the kinetic limit at small Damköhler numbers. It is shown that radiative loss assumes a similar role as varying the thermal diffusivity of the reactants. Thus, flame oscillation near the radiative limit is still thermal-diffusive in nature although it may develop under unity Lewis numbers. The unstable range of Damköhler numbers near the radiative limit shows quite similar parametric dependence on the Lewis numbers of reactants, Le_F and Le_O, the stoichiometry, ?, and the radiative loss as that near the kinetic limit. They both increase monotonically with Le_O and ? and increase then decrease with Le_F. Increasing radiative loss extends the parameter range under which flame oscillations may develop. However, they show different dependence on the temperature difference between the supplying reactants. Unless radiative loss approaches its maximum value the system can sustain, flame oscillation near the radiative limit is only possible within a limited range of ΔT, whereas it is promoted monotonically with decreasing ΔT near the kinetic limit. Furthermore, while radiative loss shows small effect on the nondimensional oscillation frequency, the dimensional frequency of flame oscillations near the radiative limit can be substantially smaller than that near the kinetic limit.     

On the stability of one-dimensional diffusion flames established between plane,parallel, porous walls

A one-dimensional, non-premixed flame stability analysis is undertaken.Oscillatory and cellular flame instabilities are identified by a careful studyof the numerically calculated eigenvalues of the linearized system of equations. The numerical investigation details the critical locations for changes in flame behaviour, as well as the critical values of variousparameters that affect flame stability. A critical Lewis number, greaterthan unity, is identified as the value where unstable oscillations maybegin to appear (Le?>?Le _c) and for which cellular flames can exist(Le?<?Le _c). Some prior discussions are clarified regarding theaforementioned critical values, as well as the role of convection inproducing flame instabilities. The methodology of the stability analysis isdiscussed in detail.     

A review of chemical diffusion: Criticism and limits of simplified methods for diffusion coefficient calculation

A review of the physics and modelling of mass diffusion involving different gaseous chemical species is firstly proposed. Both accurate and simplified models for mass diffusion involve the calculation of individual species diffusion coefficients. Since these are computationally expensive, in CFD they are commonly estimated by assuming constant Lewis or Schmidt numbers for each chemical species. The constant Lewis number assumption is particularly used. As a matter of fact, these assumptions have never been theoretically justified nor verified in practical flames. The only published information are the first observations by Smooke and Giovangigli about the Lewis number against temperature distributions in methane–air premixed and counterflow diffusion one-dimensional flames. The aim of this work is to verify these assumptions. Functional dependences of molecular properties appearing in these numbers are made explicit to show that while Sc _i depends only on composition, Le _i depends also on temperature and therefore it certainly cannot be assumed constant in a flame. Then, accurately calculating molecular properties, distributions of these characteristic numbers against temperature are obtained a posteriori from numerical simulations of different flames, premixed and non-premixed, and burning different fuels. For non-premixed flames, individual species Lewis number distributions are broad for most of the species considered in this article, whilst they are tight for premixed flames. Some attention is focused on the particular shape of Lewis distributions in non-premixed flames: they are characterized by four or five (when extinction is experienced) branches associated to precise regions in the flame (basically, lean, rich and stoichiometric combusting zones). Instead, the Schmidt distributions are always tighter, also when extinctions take place: for many species they can be approximatively assumed constant. Finally, a simplified procedure to estimate individual species diffusion coefficients is suggested, assuming the median of non-premixed flame Schmidt distributions has a constant value for each chemical species.     

Measurements of three-dimensional mean flame surface area ratio in turbulent premixed Bunsen flames

A data processing scheme with particular emphasis on proper flame contour smoothing is developed and applied to measure the three-dimensional mean flame surface area ratio in turbulent premixed flames. The scheme is based on the two-sheet imaging technique such that the mean flame surface area ratio is an average within a window covering a finite section of the turbulent flame brush. This is in contrast to the crossed-plane tomograph technique which applies only to a line. Two sets of Bunsen flames have been investigated in this work with the turbulent Reynolds number up to 4000 and the Damköhler number ranging from less than unity to close to 10. The results show that three-dimensional effects are substantial. The measured three-dimensional mean flame surface area ratio correlates well with a formula similar to the Zimont model for turbulent burning velocity but with different model constants. Also, the mean flame surface area ratio displays a weak dependency on turbulence intensity but a strong positive dependency on the turbulence integral length scale.     

Further studies of the reaction kernel structure and stabilization of jet diffusion flames

The structure of axisymmetric laminar jet diffusion flames of ethane, ethylene, acetylene, and propane in quasi-quiescent air has been studied numerically in normal earth gravity (1g) and zero gravity (0g). The time-dependent full Navier–Stokes equations with buoyancy were solved using an implicit, third-order accurate numerical scheme, including a C₃-chemistry model and an optically thin-media radiation model for heat losses. Observations of the flames were also made at the NASA Glenn 2.2-Second Drop Tower. For all cases of the fuels and gravity levels investigated, a peak reactivity spot, i.e., reaction kernel, was formed in the flame base, thereby holding a trailing diffusion flame. The location of the reaction kernel with respect to the burner rim depended inversely on the reaction-kernel reactivity or velocity. In the C₂ and C₃ hydrocarbon flames, the H₂–O₂ chain reactions were important at the reaction kernel, yet the CH₃ + O → CH₂O + H reaction, a dominant contributor to the heat-release rate in methane flames studied previously, did not outweigh other exothermic reactions. Instead of the C₁-route oxidation pathway in methane flames, the C₂ and C₃ hydrocarbon fuels dehydrogenated on the fuel side and acetylene was a major hydrocarbon fragment burning at the reaction kernel. The reaction-kernel correlations between the reactivity (the heat-release or oxygen-consumption rate) and the velocity, obtained previously for methane, were developed further for various fuels in more universal forms using variables related to local Damköhler numbers and Peclet numbers.     

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.     

Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body

A knowledge of flame stability regimes in the presence of cylindrical bluff-bodies of various dimensions is essential to design non-premixed burners. The reacting flow field in such cases is reported to be three-dimensional and unsteady. In the literature, only a few experimental investigations with limited measurements are available. Therefore, in this work, a detailed numerical study of laminar cross-flow non-premixed methane–air flames in the presence of a square cylinder is presented. The flow, temperature, species and reaction fields have been predicted using a comprehensive transient three-dimensional reacting flow model with detailed chemical kinetics and variable thermo-physical properties, in order to get a good insight into the flame stabilisation phenomena. Further, analyses of quantities such as local equivalence ratio, cell Damköhler number, species velocity, net consumption rate of methane, which are not easily obtained through experiments even with detailed diagnostics, have been carried out. The influence of the flow field due to varying inlet velocity of the oxidiser, in the presence of the bluff-body, on flame anchoring location has been analysed in detail. Local equivalence ratio contours obtained from non-reacting flow calculations are seen to be quite useful in analysing the mixing process and in the prediction of flame anchoring locations when the flames are not separated. Cell Damköhler number has been calculated using cell size, species velocity of the fuel, which is a derived quantity, and the net reaction rate of the fuel. The flame zone, which is customarily inferred from the contours of temperature, CO and OH, is also shown to be predicted well by the contour line corresponding to a Damköhler number equal to unity. The net reaction rate of CH₄ and the net rates of two dominant reactions, which consume methane, show clearly the variation in the flame anchoring locations in these three cases. Further, the three-dimensionality of these flames are analysed by plotting the mean temperature contours in y–z planes. Finally, the unsteadiness in the separated flame case is analysed.     

Structure and dynamics of edge flames in the near wake of unequal merging shear flows

We examine in this study the structure and dynamic properties of an edge flame formed in the near-wake of two initially separated shear flows, one containing fuel and the other oxidiser. A comprehensive study is carried out within the diffusive-thermal framework where the flow field, computed a-priori, is used for the determination of the combustion field. Our focus is on the effects of three controlling parameters: the Damköhler number controlling the overall flow rate, the oxidiser-to-fuel strain rate ratio of the supply streams that determines the extent of oxidiser entrainment towards the mixing zone, and the Lewis number, assumed equal for the fuel and oxidiser, that depends on the mixture composition. Response curves, representing the edge flame standoff distance as a function of Damköhler number, exhibit two distinct shapes: C-shaped and U-shaped curves characterising the response of low and high Lewis number flames, respectively. Stability considerations show that the upper solution branch of the C-shaped response curve is unstable and hence corresponds to physically unrealistic states, but due to heat conduction toward the cold plate the lower solution branch is always stable. The states forming this solution branch correspond to flame attachment, where the edge flame remains practically attached to the tip of the plate until it is blown off by the flow when the velocity exceeds a critical value. The U-shaped response, on the other hand, consists of equilibrium states that are globally stable. Thus, high Lewis number flames can be always stabilised near the splitter plate, with the edge held stationary or undergoing a back and forth motion, or lifted and stabilised downstream by the flow. Insight into the distinct stabilisation characteristics, exhibited by the different Lewis number cases, is given by examining the relationship between the local flow velocity and the edge propagation speed.     

Experimental analysis of local flame extinction in a turbulent jet diffusion flame by high repetition 2-D laser techniques and multi-scalar measurements

In this paper, we present a detailed experimental study of turbulence chemistry interactions in the “DLR_B” turbulent jet diffusion flame. The flame operates on mixtures of CH₄, H₂, and N₂ in the fuel stream at Re = 22,800 and is a target flame within the TNF workshop. Extinction and re-ignition events can be tracked in real time and related to the underlying flow field phenomena and temperature fields. Time resolved measurements of OH radical concentration fields are performed in combination with temperature and velocity field measurements. For this purpose, we combined high repetition rate (33 kHz) PLIF imaging with stereoscopic PIV and double pulse Rayleigh imaging techniques. Comparisons are made with results from multi-scalar Raman/Rayleigh/LIF point measurements that reveal the thermochemical state of the flame. The large deviations from equilibrium observed on resulting OH/temperature joint pdfs could be related to strain rate and Damköhler number variations caused by turbulent flow structures leading to frequent extinctions. The 2D measurement series uniquely reveal the underlying mechanism that can lead to such events. Finally, comparisons are made to strained laminar flame calculations, which are generally found to be in good agreement with the measured data.     

Diffusion flame instabilities in a 0.75 mm non-premixed microburner

We examine the cellular instabilities of laminar non-premixed diffusion flames that arise in a polycrystalline alumina microburner with a channel wall gap of dimension 0.75 mm. Changes in the flame structure are observed as a function of the fuel type (H₂, CH₄, and C₃H₈) and diluent. The oxidizer is O₂/inert. In contrast to previous observations on laminar diffusion flame instabilities, the current instabilities occur in the direction of flow above the splitter plate, and only occur for the heavier fuel types. They are not observed in a H₂–O₂ mixture, which will only support a continuous laminar flame inside our burner, regardless of the initial mixture strength and whether or not the flame is in near-quenching conditions. The only exception is when helium is added to the H₂–O₂ mixture, raising the effective Lewis numbers of both components.     

Gravity effects on partially premixed flames: an experimental-numerical investigation

While premixed and nonpremixed microgravity flames have been extensively investigated, the corresponding literature regarding partially premixed flames (PPFs) is sparse. We report the first experimental investigation of burner-stabilized microgravity PPFs. Partially premixed flames with multiple reaction zones are established in microgravity on a Wolfhard–Parker slot burner in the 2.2 s drop tower at the NASA Glenn Research Center. Microgravity measurements include flame imaging, and thermocouple and radiometer data. Detailed simulations are also used to provide further insight into the steady and transient response of these flames to variations in g. The flame topology and interactions between the various reaction zones are strongly influenced by gravity. The flames widen substantially in microgravity. During the transition from normal to microgravity, the flame structure experiences a fast change and another relatively slower transient change. The fast response is due to the altered advection as the value of g is reduced, while the slow response is due to the changes in the diffusive fluxes. The radiative heat loss from the flames increases in microgravity. A scaling analysis based on a radiation Damköhler number is able to characterize the radiation heat loss.     

Experimental investigation of combustion instabilities of a mesoscale multinozzle array in a lean-premixed combustor

Self-excited combustion instabilities in a mesoscale multinozzle array, also referred to as a micromixer-type injector, have been experimentally investigated in a lean-premixed tunable combustor operating with preheated methane and air. The injector assembly consists of sixty identical swirl injectors of 6.5 mm inner diameter, which are evenly distributed across the combustor dump plane. Their flow paths are divided into two groups – inner and outer stages – to form radially stratified reactant stoichiometry for the control of self-excited instabilities. OH PLIF measurements of stable flames reveal that the presence of radial staging has a remarkable influence on stabilization mechanisms, reactant jet penetration/merging, and interactions between adjacent flame fronts. In an inner enrichment case, two outer (leaner) streams merge into a single jet structure, whereas the inner (richer) reactant jets penetrate far downstream without noticeable interactions between neighboring flames. The constructed stability map in the 〈?_i, ?_o〉 domain indicates that strong self-excited instabilities occur under even split and outer enrichment conditions at relatively high global equivalence ratios. This is attributed to large-scale flame surface deformation in the streamwise direction, as manifested by vigorous detachment/attachment movements. The use of the inner fuel staging method was found, however, to limit the growth of large-amplitude heat release rate fluctuations, because the center flames are securely anchored during the whole period of oscillation, giving rise to a moderate lateral motion. We demonstrate that the collective motion of sixty flames – rather than the individual local flame dynamics – play a central role in the development of limit cycle oscillations. This suggests that the distribution pattern of the injector array, in combination with the radial fuel staging scheme, is the key to the control of the instabilities.     

Influence of radiation losses on the stability of premixed flames on a porous-plug burner

Analysis of the planar premixed flames on a porous plug was performed numerically for finite activation energy within the diffusive-thermal model. The paper is focused on the influence of radiation heat loses on the flame standoff distance and its linear stability. We show that the presence of volumetric heat losses limits the range of the mass flow range as well as it can promote the flame instabilities of different kinds, both oscillatory and cellular. The oscillatory instability, which for freely propagating flames can be usually observed for the Lewis number larger than one, in the porous-plug case occurs also for flames with unity and lower than unity Lewis number. For flames with Le < 1 both cellular and oscillatory instabilities can be observed simultaneously in a certain range of the mass flow rate.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

The role of unequal diffusivities in ignition and extinction fronts in strained mixing layers

We have studied flame propagation in a strained mixing layer formed between a fuel stream and an oxidizer stream, which can have different initial temperatures. Allowing the Lewis numbers to deviate from unity, the problem is first formulated within the framework of a thermo-diffusive model and a single irreversible reaction. A compact formulation is then derived in the limit of large activation energy, and solved analytically for high values of the Damköhler number. Simple expressions describing the flame shape and its propagation velocity are obtained. In particular, it is found that the Lewis numbers affect the propagation of the triple flame in a way similar to that obtained in the studies of stretched premixed flames. For example, the flame curvature determined by the transverse enthalpy gradients in the frozen mixing layer leads to flame-front velocities which grow with decreasing values of the Lewis numbers.

The analytical results are complemented by a numerical study which focuses on preferential-diffusion effects on triple flames. The results cover, for different values of the fuel Lewis number, a wide range of values of the Damköhler number leading to propagation speeds which vary from positive values down to large negative values     

Numerical simulation of premixed H2–air cellular tubular flames

The detailed flame structure of laminar premixed cellular flames in the tubular domain is simulated in 2D using a fully-implicit primitive variable finite difference formulation that includes multicomponent transport and detailed chemical kinetics. Numerical results for H₂/air flames are presented and compared against spatially resolved experimental measurements of temperature and chemical species including atomic H and OH. The experimental results compare well for flame structure and cell number, despite the numerical model under-predicting the peak temperature by 200 K. Numerical experiments were performed to assess the ability for cellular tubular flames to impact experimental and numerical investigations of practical flames. The cellular flame structure is found to provide a highly sensitive geometry that is useful for validating diffusive transport modelling approximations. This capability is exemplified through the development of a simple and accurate approximation for thermal diffusion (i.e. the Soret effect) that is suitable for practical combustion codes.     

Unsteady effects on flame extinction limits during gaseous and two-phase flame/vortex interactions

In highly fluctuating flows, it happens that high values of the strain-rate do not induce extinction of the flame front. Unsteady effects minimize the flame response to rapidly varying strain fields. In the present study, the effects of time-dependent flows on non-premixed flames are investigated during flame/vortex interactions. Gaseous flames and spray flames in the external sheath combustion regime are considered. To analyse the flame/vortex interaction process, the velocity field and the flame geometry are simultaneously determined using particle imaging velocimetry and laser-induced fluorescence of the CH radical. The influence of vortex flows on the extinction limits for different vortex parameters and for different gaseous and two-phase flames is examined. If the external perturbation is applied over an extended period of time, the extinction strain-rate is that corresponding to the steady-state flame, and this critical value mainly depends on the fuel and oxidizer compositions and the injection temperature. If the external perturbation is applied during a short period of time, extinction occurs at strain-rates above the steady-state extinction strain-rate. This deviation appears for flow fluctuation timescales below steady flame diffusion timescales. This behaviour is induced by diffusive processes, limiting the ability of the flame to respond to highly fluctuating flows. With respect to unsteady effects, the spray flames investigated in this article behave essentially like gaseous flames, because evaporation takes place in a thin layer before the flame front. Extinction limits are only slightly modified by the spray, controlling process being the competition between aerodynamic and diffusive timescales.