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.     

Compressible turbulent flame speeds of highly turbulent standing flames

In this paper we present the first measurement of turbulent burning velocities of a highly turbulent compressible standing flame induced by shock-driven turbulence in a Turbulent Shock Tube. High-speed schlieren, chemiluminescence, PIV, and dynamic pressure measurements are made to quantify flame–turbulence interaction for high levels of turbulence at elevated temperatures and pressure. Distributions of turbulent velocities, vorticity and turbulent strain are provided for regions ahead and behind the standing flame. The turbulent flame speed is directly measured for the high-Mach standing turbulent flame. From measurements of the flame turbulent speed and turbulent Mach number, transition into a non-linear compressibility regime at turbulent Mach numbers above 0.4 is confirmed, and a possible mechanism for flame generated turbulence and deflagration-to-detonation transition is established.     

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.     

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.     

Using self-similar properties of turbulent premixed flames to downsize chemical tables in high-performance numerical simulations

Detailed chemical mechanisms have to be incorporated in turbulent combustion modelling to predict flame propagation, ignition, extinction or pollutant formation. Unfortunately, hundreds of species and thousands of elementary reactions are involved in hydrocarbon chemical schemes and cannot be handled in practical simulations, because of the related computational costs and the need to model the complexity of their interaction with turbulent motions. Detailed chemistry may be handled using look-up tables, where chemical parameters such as reaction rates and/or species mass fractions are determined from a reduced set of coordinates, progress variables or mixture fractions, as proposed in ILDM, FPI or FGM methods. Nevertheless, these tables may require large computer memory spaces and non-negligible access times. This issue becomes of crucial importance when running on massively parallel computers: to implement these databases in shared memories would induce a large number of data exchanges, reducing the overall code performance; on the other hand duplicating databases in every local processor memory may become impossible either for large databases or small local memories. This work proposes to take advantage of the self-similar behaviour of turbulent premixed flames to reduce the size of these chemical databases, specifically when running on massively parallel machines, under the FPI (Flame Prolongation of ILDM) framework. Several approaches to reduce the database are investigated and discussed both in terms of memory requirements and access times. A very good compromise is obtained for methane–air turbulent premixed flames, where the size of the database is decreased by a factor of 1000, while the access time is reduced by about 60%.     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

Large eddy simulations of turbulent flames using the filtered density function model

Large eddy simulations (LES) of the Sandia/Sydney swirl burners (SM1 and SMA1) and the Sandia/Darmstadt piloted jet diffusion flame (Flame D) are performed. These flames are part of the database of turbulent reacting flows widely considered as benchmark test cases for validating turbulent-combustion models. In the simulations presented in this paper, the subgrid scale (SGS) closure model adopted for turbulence-chemistry interactions is based on the transport filtered density function (FDF) model. In the FDF model, the transport equation for the joint probability density function (PDF) of scalars is solved. The main advantage of this model is that the filtered reaction rates can be exactly computed. However, the density field, computed directly from the FDF solver and needed in the hydrodynamic equations, is noisy and causes numerical instability. Two numerical approaches that yield a smooth density field are examined. The two methods are based on transport equations for specific sensible enthalpy (h_s) and RT, where R is the gas constant and T is the temperature. Consistency of the two methods is assessed in a bluff-body configuration using Reynolds averaged Navier-Stokes (RANS) methodology in conjunction with the transported PDF method. It is observed that the h_s method is superior to the RT method. Both methods are used in LES of the SM1 burner. In the near-field region, the h_s method produces better predictions of temperature. However, in the far-field region, both methods show deviation from data. Simulations of the SMA1 burner and Flame D are also presented using the h_s method. Some deficiencies are seen in the predictions of the SMA1 burner that may be related to the simple chemical kinetics model and mixing model used in the simulations. Simulations of Flame D show good agreement with data. These results indicate that, while further improvements to the methodology are needed, the LES/FDF method has the potential to accurately predict complex turbulent flames.     

Large scale effects in weakly turbulent premixed flames

In this study we numerically investigate large scale premixed flames in weakly turbulent flow fields. A large scale flame is classified as such based on a reference hydrodynamic lengthscale being larger than a neutral (cutoff) lengthscale for which the hydrodynamic or Darrieus–Landau (DL) instability is balanced by stabilizing diffusive effects. As a result, DL instability can develop for large scale flames and is inhibited otherwise. Direct numerical simulations of both large scale and small scale three-dimensional, weakly turbulent flames are performed at constant Karlovitz and turbulent Reynolds number, using two paradigmatic configurations, namely a statistically planar flame and a slot Bunsen flame. As expected from linear stability analysis, DL instability induces its characteristic cusp-like corrugation only on large scale flames. We therefore observe significant morphological and topological differences as well as DL-enhanced turbulent flame speeds in large scale flames. Furthermore, we investigate issues related to reaction rate modeling in the context of flame surface density closure. Thicker flame brushes are observed for large scale flames resulting in smaller flame surface densities and overall larger wrinkling factors.     

Modelling of turbulent scalar flux in turbulent premixed flames based on DNS databases

A transport equation for scalar flux in turbulent premixed flames was modelled on the basis of DNS databases. Fully developed turbulent premixed flames were obtained for three different density ratios of flames with a single-step irreversible reaction, while the turbulent intensity was comparable to the laminar burning velocity. These DNS databases showed that the countergradient diffusion was dominant in the flame region. Analyses of the Favre-averaged transport equation for turbulent scalar flux proved that the pressure related terms and the velocity–reaction rate correlation term played important roles on the countergradient diffusion, while the mean velocity gradient term, the mean progress variable gradient term and dissipation terms suppressed it. Based on these analyses, modelling of the combustion-related terms was discussed. The mean pressure gradient term and the fluctuating pressure term were modelled by scaling, and these models were in good agreement with DNS databases. The dissipation terms and the velocity–reaction rate correlation term were also modelled, and these models mimicked DNS well.     

Three-dimensional direct simulations and structure of expanding turbulent methane flames

Direct numerical simulations (DNS) are ideally suited to investigate in detail turbulent reacting flows in simple geometries. For an increasing number of applications, detailed models must be employed to describe the chemical processes with sufficient accuracy. Despite the huge cost of such simulations, recent progress has allowed the direct numerical simulation of turbulent premixed flames while employing complete reaction schemes. We briefly describe our own developments in this field and use the resulting DNS code to investigate more extensively the structure of premixed methane flames expanding in a three-dimensional turbulent velocity field, initially homogeneous and isotropic. This situation typifies, for example, the initial flame development after spark ignition in a gas turbine or an internal combustion engine. First investigation steps have been carried out at low turbulence levels on this same configuration in the past Symposium, and we build on top of these former results. Here, a considerably higher Reynolds number is considered, the simulation has been repeated twice in to limit the possibility of spurious, very specific results, and several complementary post-processing steps are carried out. Characteristic features concerning the observed combustion regime are presented. We then investigate in a quantitative manner the evolution of flame surface area, global stretch-rate, flame front curvature, flame thickness, and correlation between thickness and curvature. The possibility of obtaining reliable information on flame front curvature from two-dimensional slices is checked by comparison with the exact procedure.     

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 numerical study of auto-ignition in turbulent lifted flames issuing into a vitiated co-flow

This paper presents a numerical study of auto-ignition in simple jets of a hydrogen–nitrogen mixture issuing into a vitiated co-flowing stream. The stabilization region of these flames is complex and, depending on the flow conditions, may undergo a transition from auto-ignition to premixed flame propagation. The objective of this paper is to develop numerical indicators for identifying such behavior, first in well-known simple test cases and then in the lifted turbulent flames. The calculations employ a composition probability density function (PDF) approach coupled to the commercial CFD code, FLUENT. The in-situ-adaptive tabulation (ISAT) method is used to implement detailed chemical kinetics. A simple k–ε turbulence model is used for turbulence along with a low Reynolds number model close to the solid walls of the fuel pipe.

The first indicator is based on an analysis of the species transport with respect to the budget of convection, diffusion and chemical reaction terms. This is a powerful tool for investigating aspects of turbulent combustion that would otherwise be prohibitive or impossible to examine experimentally. Reaction balanced by convection with minimal axial diffusion is taken as an indicator of auto-ignition while a diffusive–reactive balance, preceded by a convective–diffusive balanced pre-heat zone, is representative of a premixed flame. The second indicator is the relative location of the onset of creation of certain radical species such as HO₂ ahead of the flame zone. The buildup of HO₂ prior to the creation of H, O and OH is taken as another indicator of autoignition.

The paper first confirms the relevance of these indicators with respect to two simple test cases representing clear auto-ignition and premixed flame propagation. Three turbulent lifted flames are then investigated and the presence of auto-ignition is identified. These numerical tools are essential in providing valuable insights into the stabilization behaviour of these flames, and the demarcation between processes of auto-ignition and premixed flame propagation.     

Zone conditional analysis of a freely propagating one-dimensional turbulent premixed flame

The zone conditional conservation equations are derived and validated against the DNS data of a freely propagating one-dimensional turbulent premixed flame. Conditional flow velocities are calculated by the conditional continuity and momentum equations, and a modeled transport equation for the Reynolds average reaction progress variable. An asymptotic formula for turbulent burning velocity is obtained with the effects of a finite Damköhler number accounted for as an additional factor. It is shown that flame generated turbulence is primarily due to correlations between fluctuating gas velocities and fluctuating unit normal vector on a flame surface. More investigation is required to validate general predictive capability of the derived conditional conservation equations and the relationships modeled for closure.     

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.     

The compositional structure of highly turbulent piloted premixed flames issuing into a hot coflow

Simultaneous line measurements of major species and temperature by the Raman–Rayleigh technique, combined with CO two-photon laser-induced fluorescence and crossed-plane OH planar laser-induced fluorescence have been applied to a series of flames in the Piloted Premixed Jet Burner (PPJB). The PPJB is capable of stabilizing highly turbulent premixed jet flames through the use of a stoichiometric pilot and a large coflow of hot combustion products. Four flames with increasing jet velocities and constant jet equivalence ratios are examined in this paper. The characteristics of these four flames range from stable flame brushes with reaction zones that can be described as thin and “flamelet-like” to flames that have thickened reaction zones and exhibit extinction re-ignition behaviour. Radial profiles of the mean temperature are reported, indicating the mean thermal extent of the pilot and spatial location of the mean flame brush. Measurements of carbon monoxide (CO) and the hydroxyl radical (OH) reveal a gradual decrease in the conditional mean as the jet velocity is increased and the flame approaches extinction. Experimental results for the conditional mean temperature gradient show a progressive trend of reaction zone thickening with increasing jet velocities, indicating the increased interaction of turbulence with the reaction zone at higher turbulence levels. For the compositions examined, the product of CO and OH mole fractions ([CO][OH]) is shown to be a good qualitative indicator for the net rate of production of carbon dioxide. The axial variation of [CO][OH] is shown to correlate well with the mean chemi-luminescence of the flames including the extinction re-ignition regions. The experimental findings reported in this paper further support the hypothesis of an initial ignition region followed by extinction and re-ignition regions for certain PPJB flames.     

The stabilization mechanism and structure of turbulent hydrocarbon lifted flames

The structure and stabilization mechanism of turbulent lifted non-premixed hydrocarbon flames have been investigated using combined laser imaging techniques. The techniques include Rayleigh scattering, laser induced predissociation fluorescence of OH, LIF of PAH, LIF of CH₂O, and planar imaging velocimetry. The geometrical structure of multi-reaction zones and flow field at the stabilization region have been simultaneously measured in 16 hydrocarbon flames. The data reveal the existence of triple flame structure at the stabilization region of turbulent lifted flames. Increasing the jet velocity leads to an increase of the lift-off height and to a broadening of the lift-off region. Further analysis of the stabilization criterion at the lift-off height based on the premixed nature of triple-flame propagation and flow field data has been presented and discussed.     

Axis switching in turbulent buoyant diffusion flames

Dynamics of buoyant diffusion flames from rectangular, square, and round fuel sources were investigated using direct numerical simulation (DNS). Fully three-dimensional simulations were performed employing high-order numerical methods and boundary conditions to solve governing equations for variable-density flow and finite-rate Arrhenius chemistry. Significant differences among the different cases were revealed in the vortex dynamics, entrainment rate, small-scale mixing, and consequently flame structures. Mixing and entrainment enhancement in non-circular flames in comparison with circular ones was explained using the Biot–Savart instability theory, which relates vortex dynamics to the local azimuthal curvature. An extension of the theory elucidated why rectangular flames entrain more efficiently and spread wider than square ones, although both configurations have corners. It also provided an explanation for the aspect ratio effects in the near field. In the far field, nonlinear effects were dominant and the general transport equations for vorticity were analyzed in detail. The corner effects and aspect ratio effects were shown to be augmented by the intricate interactions among vortex dynamics, combustion, and buoyancy through the various terms in the equations. The presence of corners in non-circular flames led to concentrated regions of fine-scale mixing and intense reactions centered around the corners. Moreover, the rectangular flames exhibited a different dynamic behavior from even the square one, by creating discrepancies in entrainment, mixing, and combustion between the minor and major axis directions. Increasing the aspect ratio exacerbated such directional discrepancies, and ultimately led to axis switching. It was the first time that axis switching was observed by DNS in a rectangular flame of aspect ratio 3, which raised further questions in combustion prediction and control. Finally, a unified explanation for corner and aspect ratio effects was given on the basis of the Biot–Savart instability theory and the vorticity transport equations.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

Premixed flames modelled with thermally sensitive intermediate branching kinetics

The foundations of a relatively simple two-step kinetic scheme for flame chemistry are outlined, involving a model chain branching process that should adopt the activation temperature of a rate-limiting branching reaction in order to offer a broad approximation for hydrocarbon flames. A model energetic intermediate reactant then acts as a buffer between fuel consumption and the release of heat, as the intermediate is converted into products through a completion reaction step. By taking the rate of the latter reaction to be linear in the concentration of the intermediate, which is consistent with the final state being an equilibrium in a broader chemical system, a form of the model is arrived at which admits asymptotic solutions in a thermodiffusive context with constant coefficients. These are developed to second order for large values of the activation energy of the branching reaction and are found to involve the same trends that are seen for lean methane and hydrogen flames calculated using detailed chemical and transport models. Linear stability analysis identifies the ranges of Lewis numbers in which cellular or oscillatory instability can arise, with the latter form of instability disappearing above a threshold heat of reaction. These and the underlying flame solutions themselves depend on the heat of reaction and the degree of heat loss but not on the activation temperature of the branching reaction, to leading order. Near the limit of flammability a direct parallel arises with one-step kinetic models for premixed flames.     

Nano organic carbon and soot in turbulent non-premixed ethylene flames

Spectral optical techniques are combined to characterise the distribution of large-molecule soot precursors, nanoparticles of organic carbon, and soot in two turbulent non-premixed ethylene flames with differing residence times. Laser-induced fluorescence, laser-induced incandescence and light scattering are used to define distributions across the particle size distribution. From the scattering and laser-induced emission measurements it appears that two classes of particles are formed. The first ones are preferentially formed in the fuel-rich region of the flame closer to the nozzle, have sizes of the order of few nanometers but are not fully solid particles, because the constituent molecules still maintain their individual identity exhibiting strong broadband fluorescence in the UV. The second class of particles constituted by solid particles, with sizes of the order of tens of nanometers are able to absorb a sufficient number of photons to be heated to incandescent temperatures. These larger particles are formed at larger residence times in the flame since they are the result of slow growth processes such as coagulation or carbonization. The flames are also modeled in order to produce mixture fraction maps. A new discovery is that nanoparticles of organic carbon concentration, unlike soot, does correlate well with mixture fraction, independent of position in the flame. This is likely to be a significant benefit to future modelling of soot inception processes in turbulent non-premixed flames.