Hydrodynamic and chemical scaling for blow-off dynamics of lean premixed flames stabilized on a meso-scale bluff-body

Direct numerical simulations were conducted to investigate the effect of two parameters, density ratio and laminar flame speed, on the conditions of the onset of local extinction and blow-off of lean premixed flames, stabilized on a meso-scale bluff-body in hydrogen-air and syngas–air mixtures. A total of six simulation cases were considered as isolated comparison of the two parametric effects of the fluid dynamic instability and flame time scale. For all cases under study, the general flame development towards the blow-off limit showed a sequence of five distinct modes, with possible cyclic patterns among the different modes for a range of velocity conditions. The onset of local extinction was observed during the asymmetric vortex shedding and vortex street mode. As the density ratio is decreased, the flow inunder reviewstability is promoted through the increased sinuous mode, and such behavior was properly scaled by the Strouhal number. Although the blow-off velocity is altered by the fluid dynamic effects, the condition for the onset of local extinction and blow-off was mainly dictated by the competition between flow residence time associated with the lateral flame motion and ignition delay of the local mixtures. Time scale analysis supported the validity of the findings across all the cases investigated.     

Computational study of lean premixed turbulent flames using RANSPDF and LESPDF methods

A computational study is performed on a series of four piloted, lean, premixed turbulent jet flames. These flames use the Sydney Piloted Premixed Jet Burner (PPJB), and with jet velocities of 50, 100, 150 and 200 m/s are denoted PM150, PM1100, PM1150 and PM1200, respectively. Calculations are performed using the RANSPDF and LESPDF methodologies, with different treatments of molecular diffusion, with detailed chemistry and flamelet-based chemistry modelling, and using different imposed boundary conditions. The sensitivities of the calculations to these different aspects of the modelling are compared and discussed. Comparisons are made to experimental data and to previously-performed calculations. It is found that, given suitable boundary conditions and treatment of molecular diffusion, excellent agreement between the calculations and experimental measurements of the mean and variance fields can be achieved for PM150 and PM1100. The application of a recently developed implementation of molecular diffusion results in a large improvement in the computed variance fields in the LESPDF calculations. The inclusion of differential diffusion in the LESPDF calculations provides insight on the behaviour in the near-field region of the jet, but its effects are found to be confined to this region and to the species CO, OH and H₂. A major discrepancy observed in many previous calculations of these flames is an overprediction of reaction progress in PM1150 and PM1200, and this discrepancy is also observed in the LESPDF calculations; however, a parametric study of the LESPDF mixing model reveals that, with a sufficiently large mixing frequency, calculations of these two flames are capable of yielding improved reaction progress in good qualitative agreement with the mean and RMS scalar measurements up to an x/D of 30. Lastly, the merits of each computational methodology are discussed in light of their computational costs.     

An experimental study of partially premixed flame sound

The occurrence of oscillating combustion and combustion instability has led to resurgence of interest in the causes, mechanisms, suppression, and control of combustion noise. Noise generated by enclosed flames is of greater practical interest but is more complicated than that by open flames, which itself is not clearly understood. Studies have shown that different modes of combustion, premixed and non-premixed, differ in their sound generation characteristics. However, there is lack of understanding of the region bridging these two combustion modes. This study investigates sound generation by partially premixed flames. Starting from a non-premixed flame, air was gradually added to achieve partial premixing while maintaining the fuel flow rate constant. Methane, ethylene, and ethane partially premixed flames were studied with hydrogen added for flame stabilization. The sound pressure generated by methane partially premixed flames scales with M⁵ compared to M³ for turbulent non-premixed methane flames. Also, the sound pressure generated by partially premixed flames of ethane and ethylene scales as M^4.5. With progressive partial premixing, spectra level increases at all frequencies with a greater increase in the high-frequency region compared to the low-frequency region; flames develop a peak and later a constant level plateau in the low frequency region. The partially premixed flames of methane, ethylene, and ethane generate a similar SPL as a function of equivalence ratio when the fuel volume flow rate is matched. However, when fuel mass flow rate is matched, the ethane and ethylene flames produce a similar SPL, which is lower than that produced by the methane flame.     

Hyperacceleration effects on turbulent combustion in premixed step-stabilized flames

Experimental results are presented from an investigation of the effects of large transverse accelerations on flame propagation and blowout limits in premixed step-stabilized flames. The accelerations, which exceed ±10,000 g in the present study, induce large body forces on the high-density reactants and low-density products. These body forces can substantially alter the flame propagation mechanisms and dramatically increase the flame blowout limits. Sustained centripetal accelerations a_c ≡ U²/R are created by flowing a premixed propane–air reactant stream with equivalence ratios 0.7  Φ  1.9 at various speeds U through a semicircular channel with radius R. A backward-facing step of height h on the radially outer (a_c > 0) or inner (a_c < 0) wall stabilizes the flame. For a_c > 0 the acceleration acts to force high-density reactants into the recirculation zone and low-density products into the reactant stream, while a_c < 0 forces hot products into the recirculation zone and impedes cold reactants from entering this zone. An otherwise identical straight channel provides corresponding baseline (a_c = 0) results for comparison. The flow speed U, equivalence ratio Φ, and step height h are systematically varied for a_c = 0, a_c > 0, and a_c < 0. Shadowgraph and chemiluminescence imaging show that as a_c→ +∞ the propagation of the flame across the channel becomes independent of the flame burning velocity and instead is primarily due to large-scale “centrifugal pumping” driven by the induced body forces. For a_c → −∞ the body forces effectively segregate reactants and products to produce a nearly flat flame. In both cases, for large |a_c| values the resulting blowout limits can be substantially higher than those at a_c = 0.     

Prediction of electron and ion concentrations in low-pressure premixed acetylene and ethylene flames

Flame stabilisation and extinction in a number of different flows can be affected by application of electric fields. Electrons and ions are present in flames, and because of charge separation, weak electric fields can also be generated even when there is no externally applied electric field. In this work, a numerical model incorporating ambipolar diffusion and plasma kinetics has been developed to predict gas temperature, species, and ion and electron concentrations in laminar premixed flames without applied electric fields. This goal has been achieved by combining the existing CHEMKIN-based PREMIX code with a recently developed methodology for the solution of electron temperature and transport properties that uses a plasma kinetics model and a Boltzmann equation solver. A chemical reaction set has been compiled from seven sources and includes chemiionisation, ion-molecule, and dissociative–recombination reactions. The numerical results from the modified PREMIX code (such as peak number densities of positive ions) display good agreement with previously published experimental data for fuel-rich, non-sooting, low-pressure acetylene and ethylene flames without applied electric fields.     

Flame sheet dynamics of bluff-body stabilized flames during longitudinal acoustic forcing

Bluff-body stabilized flames are susceptible to combustion instabilities due to interactions between acoustics, vortical disturbances, and the flame. In order to elucidate these flow-flame interactions during an instability, an experimental and computational investigation of the flame-sheet dynamics of a harmonically excited flame was performed. It is shown that the flame dynamics are controlled by three key processes: excitation of shear layer instabilities by the axially oscillating flow, anchoring of the flame at the bluff body, and the kinematic response of the flame to this forcing. The near-field flame features are controlled by flame anchoring and the far-field by kinematic restoration. In the near-field, the flame response grows with downstream distance due to flame anchoring, which prevents significant flame movement near the attachment point. Theory predicts that this results in linear flame response characteristics as a function of perturbation amplitude, and a monotonic growth in magnitude of the flame-sheet fluctuations near the stabilization point, consistent with the experimental data. Farther downstream, the flame response reaches a maximum and then decays due to the dissipation of the vortical disturbances and action of flame propagation normal to itself, which acts to smooth out the wrinkles generated by the harmonic flow forcing. This behavior is strongly non-linear, resulting in significant variation in far-field flame-sheet response with perturbation amplitude.     

Combined experimental and computational study of laminar, axisymmetric hydrogen–air diffusion flames

We investigate the structure of two-dimensional, axisymmetric, laminar hydrogen–air flames in which a cylindrical fuel stream is surrounded by coflowing air, using laser-diagnostic and computational methods. Spontaneous Raman scattering and coherent anti-Stokes Raman scattering (CARS) are used to measure the distributions of major species and temperature. Computationally, we solve the governing conservation equations for mass, momentum, energy, and species, using detailed chemistry and transport. The fuel is diluted with nitrogen (1:1) to reduce heat transfer to the burner, to match the zero temperature gradient at the fuel exit. Three average fuel exit velocities are studied: 18, 27, and 50 cm/s. Comparisons of the measured and computed results are performed for radial profiles at a number of axial positions, and along the axial centerline. Peak major species mole fractions and temperatures are quantitatively predicted by the computations, and the axial species profiles are predicted to within the experimental uncertainty. In the radial profiles studied, base-case computations excluding thermal diffusion of light species were in excellent agreement with the measurements. While the addition of thermal diffusion led to some discrepancy with the measured results, the magnitude of the differences was no more than 25%. The computations predicted the axial centerline profiles from the burner exit to the maximum temperature well, though the experimental temperatures in the downstream mixing region decreased somewhat faster than the computed profiles. Radiative losses are seen to be negligible in these flames, and changes in transport properties and variations in initial flow velocities generally led to only modest changes in the axial profiles. The results also show that the detailed axial profiles of major species and temperature at different fuel jet velocities scale quantitatively with the jet velocity.     

Experimental measurements of geometric properties of turbulent stratified flames

Combustion under stratified conditions is common in many systems. However, relatively little is known about the structure and dynamics of turbulent stratified flames. Two-dimensional imaging diagnostics are applied to premixed and stratified V-flames at a mean equivalence ratio of 0.77, and low turbulent intensity, within the corrugated flame range. The present results show that stratification affects the mean turbulent flame speed, structure and geometric properties. Stratification increases the flame surface density above the premixed flame levels in all cases, with a maximum reached at intermediate levels of stratification. The flame surface density (FSD) of stratified flames is higher than that of premixed flames at the same mean equivalence ratio. Under the present conditions, the FSD peaks at a stratification ratio around 3.0. The FSD curves for stratified flames are further skewed towards the product side. The distribution of flame curvature in stratified flames is broader and more symmetric relative to premixed flames, indicating an additional mechanism of curvature generation, which is not necessarily due to cusping. These experiments indicate that flame stratification affects the intrinsic behaviour of turbulent flames and suggest that models may need to be revised in the light of the current evidence.     

Nonlinear kinematic response of premixed flames to harmonic velocity disturbances

This paper analyzes the nonlinear dynamics of premixed flames responding to harmonic velocity disturbances. These nonlinear dynamics were studied by solving a constant flame speed front tracking equation for the flame’s response to harmonically oscillating velocity disturbances. The solution to these equations is used to quantify the transfer function relating the ratio of the normalized flame area to velocity fluctuations, G = (A′/A_o)/(u′/u_o), upon the amplitude of velocity oscillations, ε = u′/u_o. Due to nonlinearities, the amplitude of this transfer function relative to its linear value decreases with increasing amplitude of velocity oscillation, u′/u_o. In contrast, the transfer function phase exhibits almost no amplitude dependence. The velocity amplitude where transfer function nonlinearities become significant depends strongly upon three parameters: a Strouhal number, St = ωL_f/u_o (where L_f is the flame length), the ratio of the flame length to width, β = L_f/R, and the flame shape in the absence of perturbations (i.e., conical, inverted wedge, etc.). In the linear case, the transfer function, G, depends only upon an algebraic combination of the first two parameters, given by St₂ = St (1 + β²)/β². In general, however, G exhibits a distinct dependence upon both parameters St and β. In particular, we show that the nonlinear response of G is an intrinsically dynamic phenomenon; i.e., its quasi-steady response (St 1) is purely linear. As such, nonlinearity is enhanced with increasing Strouhal numbers. In contrast, nonlinearity is suppressed at large β values; as such, the response of a long flame remains quite similar to its linear value, even at large ε values where the flame front exhibits substantial corrugation and cusping. Finally, we show that the response of conical flames remains much more linear at comparable disturbance amplitudes than for “V” or wedge-shaped flames. These predictions are shown to be consistent with available experimental data.     

Measurements of flame orientation and scalar dissipation in turbulent partially premixed methane flames

Three turbulent flames were studied using a new experimental facility developed at Sandia National Laboratories. Line imaging of Raman and Rayleigh scattering and CO laser-induced fluorescence (LIF) yielded information on all major species, temperature, mixture fraction, and a 1D surrogate measure of scalar dissipation. Simultaneously, crossed planar OH LIF imaging provided information on the instantaneous flame orientation, allowing estimation of the full 3D (flame-normal) scalar dissipation rate. The three flames studied were methane–air piloted jet flames (Sandia flames C, D, and E), which cover a range in Reynolds number from 13,400 to 33,600. The statistics of the instantaneous flame orientation are examined in the different flames, with the purpose of studying the prevailing kinematics of isoscalar contours. The 1D and 3D results for scalar dissipation rate are examined in detail, both in the form of conditional averages and in the form of probability density functions. The effect of overall strain and Reynolds number on flame suppression and eventual extinction is also investigated, by examining the doubly conditional statistics of temperature in the form of S-shaped curves. This latter analysis reveals that double conditioning of temperature on both mixture fraction and scalar dissipation does not collapse the data from these flames onto the same curve at low scalar dissipation rates, as might be expected from simple flamelet concepts.     

Large eddy simulation of turbulent premixed flames using level-set G-equation

Level-set G-equation and stationary flamelet chemistry are used in large eddy simulation of a propane/air premixed turbulent flame stabilized by a bluff body. The aim was to study the interaction between the flame front and turbulent eddies, and in particular to examine the effect of sub-grid scale (SGS) eddies on the wrinkling of the flame surface. The results indicated that the two types of turbulence eddies—the resolved large scale eddies and the unresolved SGS eddies—have different effects on the flame. The fluctuation of the flame surface, which is responsible for the broadening of the time averaged mean flame brush by turbulence, depends on the large resolved turbulence eddies. Time averaged mean flow velocity, temperature, and major species concentrations mainly depend on the large scale resolved eddies. The unresolved SGS eddies contribute to the wrinkling at the SGS level and play an important role in the enhancement of the propagation speed of the resolved flame front. In addition, the spatially filtered intermediate species, such as radicals, and the spatially filtered reaction rates strongly depend on the small SGS eddies. The asymptotic behavior of flame wrinkling by the SGS eddies, with respect to the decrease in filter size and grid size, is investigated further using a simplified level-set equation in a model shear flow. It is shown that to minimize the influence of the SGS eddies, fine grid and filter size may have to be used.     

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.     

Wall-impinging laminar premixed n-dodecane flames under autoignitive conditions

Laminar one-dimensional (1D) flames in a stagnation flow stabilised at a wall are used to study flame–wall interaction under diesel engine conditions. The thermochemical conditions correspond to that of the Engine Combustion Network (ECN) Spray A reference case. A range of inflow velocities is considered, where the lowest inflow velocity is chosen such that the flame is detached from the inlet. The presence of a wall is shown to have a significant impact on the flame structure and emission formation. The 1D flame and homogeneous reactor results exhibit two distinct reaction zones due to low- and high-temperature chemistry (LTC and HTC, respectively). The burner-stabilised flames are overall dominated by autoignition for all inflow velocities. For the impinging jet flames, the response of the LTC reaction zone follows closely that of the burner-stabilised flames up to relatively high inflow velocities. The HTC reaction zone, however, deviates strongly from the burner-stabilised flames, already at low inflow velocities and quenches at high inflow velocities. A budget analysis revealed a strong contribution from diffusion in the HTC reaction zone, resulting in an increasing importance of deflagrative combustion as opposed to autoignition. This trend was attributed to enhanced strain rates at higher inlet velocity leading to higher gradients. Wall heat transfer was also investigated. The highest wall heat transfer rates were observed for mixtures between $\Phi =1.0$ and $\Phi =1.5$ and for inlet velocities just below the quenching limit. This was attributed jointly to the higher peak product temperatures for these mixtures and to their enhanced resilience to quenching under strain which leads to higher temperature gradients at the wall just before quenching. NO formation was studied. The highest NO formation was observed near $\Phi =1.0,$ though the response to strain rate was different for stoichiometric and rich mixtures, which was attributed to differing NO formation pathways.     

Validation of a mixture-averaged thermal diffusion model for premixed lean hydrogen flames

The mixture-averaged thermal diffusion model originally proposed by Chapman and Cowling is validated using multiple flame configurations. Simulations using detailed hydrogen chemistry are done on one-, two-, and three-dimensional flames. The analysis spans flat and stretched, steady and unsteady, and laminar and turbulent flames. Quantitative and qualitative results using the thermal diffusion model compare very well with the more complex multicomponent diffusion model. Comparisons are made using flame speeds, surface areas, species profiles, and chemical source terms. Once validated, this model is applied to three-dimensional laminar and turbulent flames. For these cases, thermal diffusion causes an increase in the propagation speed of the flames as well as increased product chemical source terms in regions of high positive curvature. The results illustrate the necessity for including thermal diffusion, and the accuracy and computational efficiency of the mixture-averaged thermal diffusion model.     

Unsteady propagation of premixed methane/propane flames in a mesoscale disk burner of variable-gaps

Unsteady flame propagation of air-premixed methane and propane flames was investigated in a new mesoscale disk burner, of which disk-gap could be precisely varied. To begin with, the quenching disk-gaps on the flammability limits were measured. In most cases, with the slight increase of the disk-gap, cellular flame structures could be generated. The initiation of such cellular structures could be explained by the thermally induced hydrodynamic instability, and it could be enhanced if the Lewis number was sufficiently small. When the disk-gap was sufficiently larger than a critical value (approximately 1.5 times the quenching distance), the cellular flame structure changed into a smooth one in the azimuthal direction. With a further increase of the disk-gap, the flame propagation velocities approached to constant values. These values were comparable to the laminar burning velocities except for the propane-rich conditions, in which much larger propagation velocities were observed. The flame stretch effects (coupled with Le-effects) within a narrow space were suspected as the reason. The structural transition of the premixed flame could be investigated successfully through various disk-gaps, from the smallest quenching-scale to the ordinary large scales via whole mesoscales including Hele–Shaw scales.     

Particle formation in opposed-flow diffusion flames of ethylene: An experimental and numerical study

An experimental and numerical study on particles inception and growth is performed in opposed-flow diffusion flames of ethylene and air characterized by different sooting tendencies. Spectrally resolved UV-visible laser induced fluorescence, laser induced incandescence and laser light scattering measurements are used to characterize different classes of combustion-generated compounds. A detailed kinetic model accounting for both gas-phase and particle formation is used. Comparison between experimental results and numerical predictions gives a qualitative view of the mechanism of particle formation in opposed-flow flames.Particle inception is the result of both chemical growth and coagulation of aromatic compounds. In the region close to the flame front where the temperature is relatively high and radicals are abundant, the particle inception is due to a chemical growth mechanism by which aromatic molecules add aromatic radicals leading to the formation of biphenyl-like structures. The growth process continues as high-molecular mass aromatics are moved away from the flame zone towards the stagnation plane by the addition of acetylene and other aromatics forming particles of increasing sizes. Graphitization of these particles and thermal annealing lead to the formation of soot particles. At relatively lower temperatures, found across the stagnation plane, particles growth still occurs and it is mainly due to a process of physical coagulation of PAHs.The experimental and numerical results obtained in this work demonstrate and explain the sensibility of inception and growth of particles to radical concentration and temperature in opposed-flow flame configurations.     

Structural response of different Lewis number premixed flames interacting with a toroidal vortex

Simultaneous measurements of temperature, CH* and OH* chemiluminescent species are carried out to explore the impact of stretch rate and curvature on the structure of premixed flames. The configuration of an initially flat premixed flame interacting with a toroidal vortex is selected for the present study and reasons for this choice are discussed. Lewis number effects are assessed by comparing methane and propane flames. It is emphasized that the flame structure experiences very strong variations. In particular, the flame is shrunk (broadened) in the initial (final) period of the interaction with the vortex where strain rate (curvature) contribution of the stretch rate is predominant. By further analysing independently the thickness of the preheat and reaction zones, it is shown that for propane flames, not only the former but also the latter is significantly altered in zones where the flame curvature is negative. Changes in the reaction zone properties are further emphasized using CH* and OH* radicals. It is demonstrated that higher thermal diffusivity plays a significant role around curved regions, in which the enhanced diffusion of heat leads to a strong increase of CH* compared to OH* intensity. As an overall conclusion, this study suggests that it would be interesting to reassess the internal flame structure at lower and moderate Karlovitz numbers since changes might appear for a moderate vortex intensity with typical size much larger than the flame thickness.     

A numerical study on the effect of CO addition on extinction limits and NO x formation in lean counterflow CH4/air premixed flames

The effect of CO addition on extinction and NO _x formation in lean premixed counterflow CH₄/air flames was investigated by numerical simulation. Detailed chemistry and complex thermal and transport properties were employed. A method that gradually switched off the initial reactions of NO formation from different routes was used to analyse the variation of NO formation mechanism. The results indicate that the addition of certain amount of CO increases the strain extinction limits and reduces the radiation extinction limits. As a result, the lean flammability limit of CH₄/air premixed flame is extended to leaner side by the addition of CO. The formation of NO in a flame is increased with the addition of CO at a constant equivalence ratio. For an ultra-lean flame, the increase in the formation of NO is mainly because of the increase in the contribution from the NNH intermediate route, while for a near stoichiometric flame, this increase is mainly attributable to the rise in the contribution from the thermal route. With the fraction of added CO being gradually increased, the formation of NO₂ in a flame first decreases and then increases at a given equivalence ratio. The addition of CO reduces the formation N₂O in an ultra-lean flame, while affects little on the formation of N₂O in a near stoichiometric flame.     

A detailed chemical kinetic reaction mechanism for oxidation of four small alkyl esters in laminar premixed flames

A detailed chemical kinetic reaction mechanism has been developed for a group of four small alkyl ester fuels, consisting of methyl formate, methyl acetate, ethyl formate, and ethyl acetate. This mechanism is validated by comparisons between computed results and recently measured intermediate species mole fractions in fuel-rich, low-pressure, premixed laminar flames. The model development employs a principle of similarity of functional groups in constraining the H atom abstraction and unimolecular decomposition reactions for each of these fuels. As a result, the reaction mechanism and formalism for mechanism development are suitable for extension to larger oxygenated hydrocarbon fuels, together with an improved kinetic understanding of the structure and chemical kinetics of alkyl ester fuels that can be extended to biodiesel fuels. Variations in concentrations of intermediate species levels in these flames are traced to differences in the molecular structure of the fuel molecules.