A Study of premixed propagating flame vortex interaction

Experimental data is presented for the interaction between a propagating flame and a simple vortex flow field structure generated in the wake of solid obstacles. The interaction between gas movement and obstacles creates vortex shedding forming a simple flow field recirculation. The presence of the simple turbulent structure within the gas mixture curls the flame front increasing curvature and enhancing burning rate. A novel twin camera Particle Image Velocimetry, PIV, was employed to characterise the flow field recirculation and the interaction with the flame front. The technique allowed the quantification of the flame/vortex interaction. The twin camera technique provides data to define the spatial variation of both the velocity of the flow field and flame front. Experimentally obtained values of local flame displacement speed and flame stretch rate are presented for simple flame/vortex interactions.     

Effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes

Experiments and theoretical analysis are presented to clarify the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames in tubes. The experiments reveal that the longer tubes have higher instability compared to shorter tubes and the lower diameter tubes have higher instability compared to higher diameter tubes. The secondary instability leading to turbulent burning is found to be more sensitive to change in geometrical parameters compared to primary instability (oscillating flat flame). The secondary instability is re-stabilized for some intermediate burning velocity conditions even though lower and higher burning velocity conditions show secondary instability. The appearance of such re-stabilization is only observed for some specific lengths of the tube. Present experimental observations pertaining to the effect of geometrical parameters is found to be contradicting the theoretical predictions based on pressure coupling mechanism. To clear the underlying mechanism, analytical growth rate is computed considering velocity coupling mechanism. The computed growth rates correctly predict the effect of geometrical parameters on thermo-acoustic instability of downward propagating flames. This work provides further evidence to believe that the flame -acoustic coupling in downward propagating flames is due to flame area modulation (leading to heat release modulation) through action of acoustic acceleration.     

Instability of a flame front propagating through a fuel-rich droplet–vapour–air cloud

A simple model of a flame front propagating through a fuel-rich droplet–vapour–air mixture is presented in which the fuel droplets are assumed to evaporate in a sharp front ahead of the reaction front. By performing a linear stability analysis neutral stability boundaries are determined. It is shown that the presence of the spray of droplets in the fresh mixture can have a profound effect by causing cellularization of the flame front. Specifically, we demonstrate that under certain circumstances a spray flame can be cellular when its equivalent non-spray flame is completely stable. Furthermore, it is shown that even when the non-spray flame is itself cellular the equivalent spray flame will have a finer cellular structure. These theoretical predictions verify qualitatively for the first time independent experimental observations from the literature. It is thus shown that the primary effect of the spray on the stability of these flames is due to heat loss from the absorption of heat by the droplets for vaporization. The influence of the initial liquid fuel loading and the latent heat of vaporization on the critical wavenumber associated with cellularity provide further evidence of the responsibility of the heat loss mechanism for these spray-related phenomena. Finally, the cellularity of the spray flames with their attendant increase in flame front area suggest a plausible rationale for the experimentally observed burning velocity enhancement induced by the use of a spray of fuel droplets.     

An assessment of large eddy simulations of premixed flames propagating past repeated obstacles

This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.     

On laminar premixed flame propagating into autoigniting mixtures under engine-relevant conditions

Usually premixed flame propagation and laminar burning velocity are studied for mixtures at normal or elevated temperatures and pressures, under which the ignition delay time of the premixture is much larger than the flame resistance time. However, in spark-ignition engines and spark-assisted compression ignition engines, the end-gas in the front of premixed flame is at the state that autoignition might happen before the mixture is consumed by the premixed flame. In this study, laminar premixed flames propagating into an autoigniting dimethyl ether/air mixture are simulated considering detailed chemistry and transport. The emphasis is on the laminar burning velocity of autoigniting mixtures under engine-relevant conditions. Two types of premixed flames are considered: one is the premixed planar flame propagating into an autoigniting DME/air without confinement; and the other is premixed spherical flame propagating inside a closed chamber, for which four stages are identified. Due to the confinement, the unburned mixture is compressed to high temperature and pressure close to or under engine-relevant conditions. The laminar burning velocity is determined from the constant-volume propagating spherical flame method as well as PREMIX. The laminar burning velocities of autoigniting DME/air mixture at different temperatures, pressures, and autoignition progresses are obtained. It is shown that the first-stage and second-stage autoignition can significantly accelerate the flame propagation and thereby greatly increase the laminar burning velocity. When the first-stage autoignition occurs in the unburned mixture, the isentropic compression assumption does not hold and thereby the traditional method cannot be used to calculate the laminar burning velocity. A modified method without using the isentropic compression assumption is proposed. It is shown to work well for autoigniting mixtures. Besides, a power law correlation is obtained based on all the laminar burning velocity data. It works well for mixtures before autoignition while improvement is still needed for mixtures after autoignition.     

Resonance of a turbulent flame in a high-frequency acoustic wave

Resonance of a weakly turbulent flame in a high-frequency acoustic wave is obtained. Because of the resonance, an acoustic wave may increase noticeably the amplitude of flame wrinkles, and the respective increase in propagation velocity of the turbulent flame front becomes larger by a factor of 10-20. The effect of resonance is especially important for turbulent flames with realistic thermal expansion propagating in a closed burning chamber, which may account for considerable scattering of experimental results on turbulent flame velocity.     

Spatiotemporal measurements of flame stretch and propagation rates for lean and rich CH4/air premixed flames interacting with a turbulent-wake

A 1.5 m long turbulent-wake combustion vessel with a 0.15 m × 0.15 m cross-sectional area is proposed for spatiotemporal measurements of curvature, strain, dilatation and burning rates along a freely downward-propagating premixed flame interacting with a parallel row of staggered vortex pairs having both compression (negative) and extension (positive) strains simultaneously. The wanted wake is generated by rapidly withdrawing an electrically-controlled, horizontally-oriented sliding plate of 5 mm thickness for flame–wake interactions. Both rich and lean CH₄/air flames at the equivalence ratios  = 1.4 and  = 0.7 with nearly the same laminar burning velocity are studied, where flame–wake interactions and their time-dependent velocity fields are obtained by high-speed, high-resolution DPIV and laser-tomography. Correlations among curvature, strain, stretch, and dilatation rates along wrinkled flame fronts at different times are measured and thus their influences on front propagation rates can be analyzed. It is found that strain-related effects have significant influence on front propagation rates of rich CH₄/air (diffusionally stable) flames even when the curvature weights more in the total stretch than the strain rate does. The local propagation rates along the wrinkled flame front are more intense at negative strain rates corresponding to positive peak dilatation rates but the global propagation rate averaged along the rich flame front remains constant during all period of flame–wake interaction. For lean CH₄/air (diffusionally unstable) flames, the curvature becomes a dominant parameter influencing the structure and propagation of the wrinkled flame front, where both local and global propagation rates increase significantly with time, showing unsteady flame propagation. These experimental results suggest that the theory of laminar flame stretch can be applicable to a more complex flame–wake interaction involving unsteadiness and multitudinous interactions between vortices.     

Stabilized flames in hybrid aluminum-methane-air mixtures

A premixed methane–air bunsen-type flame is seeded with micron-sized (d₃₂ = 5.6 μm) atomized aluminum powder over a wide range of solid fuel concentrations. The burning velocities of the resulting two-phase hybrid flame are determined using the total surface area of the inner flame cone and the known volumetric flow rate, and spatially resolved flame spectra are obtained with a spectral scanning system. Flame temperatures are derived through polychromatic fitting of Planck’s law to the continuous part of the spectrum. It is found that an increase in the solid fuel concentration changes the aluminum combustion regime from low temperature oxidation to full-fledged flame front propagation. For stoichiometric methane–air mixtures, the transition occurs in the aluminum concentration range of 140–220 g/m³ and is manifested by the appearance of AlO sub-oxide bands and an increase in the flame temperature to 2500 K. The flame burning velocity is found to decrease only slightly with an increase in aluminum concentration, in contrast to the rapid decrease in flame speed, followed by quenching, that is observed for flames seeded with inert SiC particles. The observed behavior of the burning velocity and flame temperature leads to the conclusion that intense aluminum combustion in a hybrid flame only occurs when the flame front propagating through the aluminum suspension is coupled to the methane–air flame.     

Stratified turbulent Bunsen flames: flame surface analysis and flame surface density modelling

In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification.

All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1?φ?1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate.

In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1?φ?1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4?φ?1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1?φ?1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data, thereby rendering this approach promising for application in LES.     

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.     

Effects of equivalence ratio variation on lean,stratified methane–air laminar counterflow flames

The effects of equivalence ratio variations on flame structure and propagation have been studied computationally. Equivalence ratio stratification is a key technology for advanced low emission combustors. Laminar counterflow simulations of lean methane–air combustion have been presented which show the effect of strain variations on flames stabilized in an equivalence ratio gradient, and the response of flames propagating into a mixture with a time-varying equivalence ratio. ‘Back supported’ lean flames, whose products are closer to stoichiometry than their reactants, display increased propagation velocities and reduced thickness compared with flames where the reactants are richer than the products. The radical concentrations in the vicinity of the flame are modified by the effect of an equivalence ratio gradient on the temperature profile and thermal dissociation. Analysis of steady flames stabilized in an equivalence ratio gradient demonstrates that the radical flux through the flame, and the modified radical concentrations in the reaction zone, contribute to the modified propagation speed and thickness of stratified flames. The modified concentrations of radical species in stratified flames mean that, in general, the reaction rate is not accurately parametrized by progress variable and equivalence ratio alone. A definition of stratified flame propagation based upon the displacement speed of a mixture fraction dependent progress variable was seen to be suitable for stratified combustion. The response times of the reaction, diffusion, and cross-dissipation components which contribute to this displacement speed have been used to explain flame response to stratification and unsteady fluid dynamic strain.     

Premixed turbulent flame front structure investigation by Rayleigh scattering in the thin reaction zone regime

Premixed turbulent flames of methane–air and propane–air stabilized on a bunsen type burner were studied using planar Rayleigh scattering and particle image velocimetry. The fuel–air equivalence ratio range was from lean 0.6 to stoichiometric for methane flames, and from 0.7 to stoichiometric for propane flames. The non-dimensional turbulence rms velocity, u′/S_L, covered a range from 3 to 24, corresponding to conditions of corrugated flamelets and thin reaction zones regimes. Flame front thickness increased slightly with increasing non-dimensional turbulence rms velocity in both methane and propane flames, although the flame thickening was more prominent in propane flames. The probability density function of curvature showed a Gaussian-like distribution at all turbulence intensities in both methane and propane flames, at all sections of the flame.The value of the term Dκ, the product of molecular diffusivity evaluated at reaction zone conditions and the flame front curvature, has been shown to be smaller than the magnitude of the laminar burning velocity. This finding questions the validity of extending the level set formulation, developed for corrugated flames region, into the thin reaction zone regime by increasing the local flame propagation by adding the term Dκ to laminar burning velocity.     

A flamelet analysis of the burning velocity of premixed turbulent expanding flames

Direct numerical simulation is a very powerful tool to evaluate the validity of new models and theories for turbulent combustion. In this paper, direct numerical simulations of spherically expanding premixed turbulent flames in the corrugated flamelet regime are performed. The flamelet-generated manifold method is used to deal with detailed reaction kinetics. The numerical method is validated for both laminar and turbulent expanding flames. The computational results are analyzed by using an extended flame stretch theory. It is investigated whether this theory is able to describe the influence of flame stretch and curvature on the local burning velocity of the flame. If the full profiles of flame stretch and curvature through the flame front are included in the theory, the local mass burning rate is predicted accurately. The influence of several approximations, which are used in other existing theories, is studied. When flame stretch is assumed to be constant through the flame front or when curvature of the flame front is neglected, the theory fails to predict the local mass burning rate.     

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.     

How fast can we burn, 2.0

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

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.     

Experimental measurements of two-dimensional planar propagating edge flames

The study of edge flames has received increased attention in recent years. This work reports the results of a recent study into two-dimensional, planar, propagating edge flames that are remote from solid surfaces (called here, “free-layer” flames, as opposed to layered flames along floors or ceilings). They represent an ideal case of a flame propagating down a flammable plume, or through a flammable layer in microgravity. The results were generated using a new apparatus in which a thin stream of gaseous fuel is injected into a low-speed laminar wind tunnel thereby forming a flammable layer along the centerline. An airfoil-shaped fuel dispenser downstream of the duct inlet issues ethane from a slot in the trailing edge. The air and ethane mix due to mass diffusion while flowing up towards the duct exit, forming a flammable layer with a steep lateral fuel concentration gradient and smaller axial fuel concentration gradient. We characterized the flow and fuel concentration fields in the duct using hot wire anemometer scans, flow visualization using smoke traces, and non-reacting, numerical modeling using COSMOSFloWorks. In the experiment, a hot wire near the exit ignites the ethane-air layer, with the flame propagating downwards towards the fuel source. Reported here are tests with the air inlet velocity of 25 cm/s and ethane flows of 967-1299 sccm, which gave conditions ranging from lean to rich along the centerline. In these conditions the flame spreads at a constant rate faster than the laminar burning rate for a premixed ethane-air mixture. The flame spread rate increases with increasing transverse fuel gradient (obtained by increasing the fuel flow rate), but appears to reach a maximum. The flow field shows little effect due to the flame approach near the igniter, but shows significant effect, including flow reversal, well ahead of the flame as it approaches the airfoil fuel source.     

Analysis of the flame thickness of turbulent flamelets in the thin reaction zones regime

The thickness of the instantaneous flamelets in a turbulent flame brush on a weak-swirl burner burning in the thin reaction zones regime has been analysed experimentally, theoretically, and numerically. The experimental flame thickness has been measured correlating two simultaneous Rayleigh images and one OH-image from two closely spaced cross sections in the flame. It appears that the low temperature edge of the flame is thickened by turbulent eddies but that these structures cannot penetrate far enough into the flame front to distort the inner layer for the moderate Karlovitz numbers used. The flame front based on the temperature gradient at the inner layer becomes thinner for lean flames and thicker for rich methane–air flames. This has been explained theoretically and numerically by studying the influence of flame stretch and preferential diffusion on the flame thickness. It appears that the flame front thickness at the inner layer (and mass burning rate) is not influenced by turbulent mixing processes, and it seems that eddies of the size of the inner layer have to be used to change this picture. Experiments closer to the boundary of the broken reaction zones regime have to confirm this in the future.     

Partially premixed flamelet in LES of acetone spray flames

An effective partially premixed flamelet model for large eddy simulation (LES) of turbulent spray combustion is formulated. Different flame regimes are identified with a flame index defined by budget terms in a 2-D multi-phase flamelet formulation, and the application in LES of partially pre-vaporized spray flames shows a favorable agreement with experiments. Simulations demonstrate that, compared to the conventional single-regime flamelets, the present partially premixed flamelet formulation shows its ability in capturing the subgrid regime transitions, yielding a well prediction of peak gas temperature and the downstream flame spreading. A propagating premixed flame front is found coupled with a trailing diffusion burning through the spray evaporation, and the spray effect on regime discrimination is manifested with transport budget analysis. A two-phase regime indicator is then proposed, by which the evaporation-dictated regime is properly described. Its intended use will rely on both gas and spray flamelet structures.