An EDC-based turbulent premixed combustion model

In the present study, a new turbulent premixed combustion model is proposed by integrating the Coherent Flame Model with the modified eddy dissipation concept, and relating the fine structure mass fraction to the flame surface density. First, experimental results of turbulent flame speed available from literature are compared with the predicted results at different turbulence intensities to validate the flame surface density model. It is observed that the model is able to predict the turbulent burning speeds accurately. Then, a comprehensive validation is carried out utilizing data on a turbulent lifted methane flame issuing into a vitiated co-flow. Detailed comparison of temperature and species concentrations between experiment and simulation is performed at different heights of the flame. Overall, the model is found to predict both the spatial variation and peak values of the scalars at various heights satisfactorily.     

Effects of confinement on premixed turbulent swirling flame using large Eddy simulation

This paper utilises large eddy simulation (LES) to study swirling reacting flows by comparison with experimental observations. The purpose is to provide further insights in engineering designs, as well as to improve modelling. A reduced-scale swirl burner has been developed for the experiments. Comparison of particle image velocimetry (PIV) measurements with LES results using finite rate chemistry shows that LES captures all the salient features of an unconfined flame including velocity and temperature distributions. However, when the flame is confined within a cylindrical combustor, the simulated flame shape is initially not consistent with experimental observation. Investigations show that the discrepancy is caused by the often practised assumption of adiabatic wall temperature. With the use of an assumed wall temperature distribution guided by laboratory observation, results of LES are consistent with experiments. Although the latter LES approach requires more computational resources, the improvement is found to be justified.     

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.     

The role of duct thickness on the quenching process of premixed flame propagation

The propagation of premixed laminar flame in ducts of circular cross-section considering a thermal-diffusive model is investigated numerically. Heat losses by conduction to the channels walls are taken into account using the thermally thin wall regime. The effects and the relationship between thickness and diameter of the tube with the flame speed propagation are studied and the quenching condition is obtained as a function of the heat-loss parameter. The mathematical model employs the axisymmetric energy and species equations. The calculations are based on a two-step chemistry, with an Arrhenius, energetically neutral, radical production reaction followed by an exothermic radical recombination reaction. For large values of the heat-loss parameter, the wall temperature is close to the free stream temperature and all the heat losses through the wall are convected away. No heat feedback occurs. On the other hand, for small values of the heat-loss parameter, a feedback mechanism occurs by transferring heat from the burned gas to the fresh mixture along the tube wall. For values of the heat-loss parameter of order unity, the heat feedback mechanism is able to sustain the flame propagation and the quenching condition disappears, producing an almost planar flame front as the propagation velocity reduces. For this two-step reaction mechanism, the radical species behaviour at the duct walls seems to have negligible effect on the quenching process.     

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.     

Characterization of premixed swirling methane/air diffusion flame through filtered Rayleigh scattering

Characteristics of a premixed, swirl methane/air diffusion flame at atmospheric pressure are measured by filtered Rayleigh scattering (FRS). Three operating conditions are investigated with the equivalence ratios of the methane/air flame covering a range of 0.67—0.83. Under each condition, single-shot and averaged FRS images over a region measured 39.3×65.6 mm² at seven cross sections of the flame are collected to demonstrate the flame behavior. A gradient calculation algorithm is applied to identify reaction zone locations and structures in the instantaneous FRS measurements. Statistical analysis for the mean FRS measurements is performed by means of joint probability density functions. The experimental results indicate that thermochemical state of the swirl flame is strongly influenced by equivalence ratio, leading to varieties of flame structures and temperature distributions. The gradient of the instantaneous FRS images clearly illustrates the characteristics of the reaction zone. The results also demonstrate that FRS can provide detailed insights into the behavior of turbulent flames.     

Effect of stoichiometric mixture fraction on soot fraction and emission spectra with application to oxy-combustion

Many proposed oxy-combustion concepts for carbon capture incorporate the recycling of flue gas which is used as a dilution gas to aid in the control of temperature and heat flux. Improvements in efficiency may be realized by significantly reducing the recycle flue gas (RFG), however, in application, care must be taken to avoid excessive radiant heat flux and gas temperature. One of the features oxy-combustion, unlike air-fired combustion, is that the oxygen and dilution gases are initially separated. RFG can, for example, be strategically blended with either the fuel stream, or oxidizer stream, or both, which affects the stoichiometric mixture fraction, Z_st. In this work, the effects of the amount of dilution, or RFG, and Z_st on soot fraction are experimentally investigated in a laminar coflow flame. Carbon dioxide is employed as the dilution gas to simulate the recycling of dry flue gas. Soot fraction and temperature are quantitatively determined by a flame image processing technique. In addition, the visible and near-IR emission spectra are given. When dilution, or RFG, is reduced, while holding Z_st constant, soot formation and thermal radiation increase due to higher temperature. However, high temperature flames with reduced or zero soot are achieved by increasing Z_st via the combination of fuel dilution and oxygen enrichment. This study highlights the inherent flexibility of oxy-fuel combustion, which offers the opportunity to control flame temperature and gas volume while independently controlling soot formation and radiant heat transfer.     

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.     

Nonlinear hydrodynamic and thermoacoustic oscillations of a bluff-body stabilised turbulent premixed flame

Turbulent premixed flames often experience thermoacoustic instabilities when the combustion heat release rate is in phase with acoustic pressure fluctuations. Linear methods often assume a priori that oscillations are periodic and occur at a dominant frequency with a fixed amplitude. Such assumptions are not made when using nonlinear analysis. When an oscillation is fully saturated, nonlinear analysis can serve as a useful avenue to reveal flame behaviour far more elaborate than period-one limit cycles, including quasi-periodicity and chaos in hydrodynamically or thermoacoustically self-excited system. In this paper, the behaviour of a bluff-body stabilised turbulent premixed propane/air flame in a model jet-engine afterburner configuration is investigated using computational fluid dynamics. For the frequencies of interest in this investigation, an unsteady Reynolds-averaged Navier–Stokes approach is found to be appropriate. Combustion is represented using a modified laminar flamelet approach with an algebraic closure for the flame surface density. The results are validated by comparison with existing experimental data and with large eddy simulation, and the observed self-excited oscillations in pressure and heat release are studied using methods derived from dynamical systems theory. A systematic analysis is carried out by increasing the equivalence ratio of the reactant stream supplied to the premixed flame. A strong variation in the global flame structure is observed. The flame exhibits a self-excited hydrodynamic oscillation at low equivalence ratios, becomes steady as the equivalence ratio is increased to intermediate values, and again exhibits a self-excited thermoacoustic oscillation at higher equivalence ratios. Rich nonlinear behaviour is observed and the investigation demonstrates that turbulent premixed flames can exhibit complex dynamical behaviour including quasiperiodicity, limit cycles and period-two limit cycles due to the interactions of various physical mechanisms. This has implications in selecting the operating conditions for such flames and for devising proper control strategies for the avoidance of thermoacoustic instability.     

Stabilization mechanism of turbulent premixed flames in strongly swirled flows

The stabilization of turbulent premixed flames in strongly swirled flows undergoing vortex breakdown is studied in the case of the ALSTOM En-Vironmental (EV) double cone burner using a simple one-dimensional boundary layer type model and computational fluid dynamics, mainly at the level of large-eddy simulation. The analysis shows that, due to flame curvature effects, the flame speed on the combustor axis is 2 D _t/R _F lower than the turbulent burning rate, where D _t is a characteristic turbulent diffusion coefficient and R _F the flame radius of curvature. Flame propagation with negative speed observed in the experiments, i.e. the flame completely embedded in the central recirculation zone on the symmetry axis, is explained with the one-dimensional model as caused by the factor 2 D _t/R _F being larger than the characteristic turbulent burning rate. A peculiar sudden displacement of the flame anchoring location deep into the burner, which takes place experimentally at a critical value of the equivalence ratio, cannot however be explained with the present one-dimensional approach due to the modelling assumptions. The mathematical analysis is supported in this case with large-eddy simulation which can accurately reproduce the flame behaviour across the full operating range. It is finally shown that steady RANS methods cannot cope with the problem due to their inability to correctly predict the velocity flowfield in this burner.     

Simulations and experiments on the ignition probability in turbulent premixed bluff-body flames

The ignition characteristics of a premixed bluff-body burner under lean conditions were investigated experimentally and numerically with a physical model focusing on ignition probability. Visualisation of the flame with a 5 kHz OH* chemiluminescence camera confirmed that successful ignitions were those associated with the movement of the kernel upstream, consistent with previous work on non-premixed systems. Performing many separate ignition trials at the same spark position and flow conditions resulted in a quantification of the ignition probability P_ign, which was found to decrease with increasing distance downstream of the bluff body and a decrease in equivalence ratio. Flows corresponding to flames close to the blow-off limit could not be ignited, although such flames were stable if reached from a richer already ignited condition. A detailed comparison with the local Karlovitz number and the mean velocity showed that regions of high P_ign are associated with low Ka and negative bulk velocity (i.e. towards the bluff body), although a direct correlation was not possible. A modelling effort that takes convection and localised flame quenching into account by tracking stochastic virtual flame particles, previously validated for non-premixed and spray ignition, was used to estimate the ignition probability. The applicability of this approach to premixed flows was first evaluated by investigating the model's flame propagation mechanism in a uniform turbulence field, which showed that the model reproduces the bending behaviour of the S_T-versus-u′ curve. Then ignition simulations of the bluff-body burner were carried out. The ignition probability map was computed and it was found that the model reproduces all main trends found in the experimental study.     

Enhancement of convective heat transfer by electrically-induced swirling effect in laminar and fully-developed internal flows

A computational and experimental approach is undertaken to study the enhancement of convective heat transfer in fully-developed internal flows by an electrically-induced secondary flow field. Using longitudinal flat electrodes along a parallel-plate configuration, the corona discharge is employed to generate an electrically induced secondary flow on the cross section of the flow passage. The electrically-induced secondary flow forms a swirling flow field in the fully-developed condition and enhances the heat transfer significantly. The flow field was solved computationally and the results were verified and validated by grid refinement study and computational error analysis.     

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.     

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.     

Diffusive–thermal oscillations of rich premixed hydrogen–air flames in a microflow reactor

In this paper the dynamics of rich hydrogen–air flames in a microflow reactor with controlled temperature of the walls is investigated numerically using the thermal-diffusion model with two-step kinetics in one spatial dimension. It is found that as the parameters of the system are varied the sequence of bifurcation occurs leading to the formation of complex spatio-temporal patterns. These include pulsating, chaotic, mixed-mode and FREI (Flames with Repetitive Extinction and Ignition) oscillations. The critical parameter values for the existence of different dynamical regimes are found in terms of equivalence ratio and flow velocity.     

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.     

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.