Large-eddy simulation of a bluff-body stabilised turbulent premixed flame using the transported flame surface density approach

A premixed propane–air flame stabilised on a triangular bluff body in a model jet-engine afterburner configuration is investigated using large-eddy simulation (LES). The reaction rate source term for turbulent premixed combustion is closed using the transported flame surface density (TFSD) model. In this approach, there is no need to assume local equilibrium between the generation and destruction of subgrid FSD, as commonly done in simple algebraic closure models. Instead, the key processes that create and destroy FSD are accounted for explicitly. This allows the model to capture large-scale unsteady flame propagation in the presence of combustion instabilities, or in situations where the flame encounters progressive wrinkling with time. In this study, comprehensive validation of the numerical method is carried out. For the non-reacting flow, good agreement for both the time-averaged and root-mean-square velocity fields are obtained, and the Karman type vortex shedding behaviour seen in the experiment is well represented. For the reacting flow, two mesh configurations are used to investigate the sensitivity of the LES results to the numerical resolution. Profiles for the velocity and temperature fields exhibit good agreement with the experimental data for both the coarse and dense mesh. This demonstrates the capability of LES coupled with the TFSD approach in representing the highly unsteady premixed combustion observed in this configuration. The instantaneous flow pattern and turbulent flame behaviour are discussed, and the differences between the non-reacting and reacting flow are described through visualisation of vortical structures and their interaction with the flame. Lastly, the generation and destruction of FSD are evaluated by examining the individual terms in the FSD transport equation. Localised regions where straining, curvature and propagation are each dominant are observed, highlighting the importance of non-equilibrium effects of FSD generation and destruction in the model afterburner.     

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.     

Dynamic behavior of a freely-propagating turbulent premixed flame under global stretch-rate oscillations

Dynamic features of a freely propagating turbulent premixed flame under global stretch rate oscillations were investigated by utilizing a jet-type low-swirl burner equipped with a high-speed valve on the swirl jet line. The bulk flow velocity, equivalence ratio and the nominal mean swirl number were 5 m/s, 0.80 and 1.23, respectively. Seven velocity forcing amplitudes, from 0.09 to 0.55, were examined with a single forcing frequency of 50 Hz. Three kinds of optical measurements, OH-PLIF, OH^* chemiluminescence and PIV, were conducted. All the data were measured or post-processed in a phase-locked manner to obtain phase-resolved information. The global transverse stretch rate showed in-phase oscillations centering around 60 (1/s). The oscillation amplitude of the stretch rate grew with the increment of the forcing amplitude. The turbulent flame structure in the core flow region varied largely in axial direction in response to the flowfield oscillations. The flame brush thickness and the flame surface area oscillated with a phase shift to the stretch rate oscillations. These two properties showed a maximum and minimum values in the increasing and decreasing stretch periods, respectively, for all the forcing amplitudes. Despite large variations in flame brush thickness at different phase angles, the normalized profiles collapse onto a consistent curve. This suggests that the self-similarity sustains in this dynamic flame. The global OH^* fluctuation response (i.e. response of global heat-release rate fluctuation) showed a linear dependency to the forcing velocity oscillation amplitudes. The flame surface area fluctuation response showed a linear tendency as well with a slope similar to that of the global OH^* fluctuation. This indicated that the flame surface area variations play a critical role in the global flame response.     

A theoretical study of premixed turbulent flame development

Flame development in a statistically stationary and uniform, planar, one-dimensional turbulent flow is theoretically studied. A generalized balance equation for the mean combustion progress variable, which includes turbulent diffusion and pressure-driven transport terms, as well as the mean rate of product creation, is introduced and analyzed by invoking the sole assumption of a self-similar flame structure, well-supported by numerous experiments. The assumption offers the opportunity to simplify the problem by splitting the aforementioned partial differential equation into two ordinary differential equations, which separately model spatial variations of the progress variable and time variations of flame speed and thickness. The self-similar profile of the progress variable, obtained in numerous experiments, is theoretically predicted. Closures of the normalized pressure-driven transport term and mean rate of product creation are obtained. The closed balance equation shows that turbulent diffusion dominates during the initial stage of flame development, followed by the transition to counter-gradient transport in a sufficiently developed flame. A criterion of the transition is derived. The transition is promoted by the heat release and pressure-driven transport. Fully developed mean flame brush thickness and speed are shown to decrease when either density ratio or pressure-driven transport increases. Solutions for the development of the thickness are obtained. The development is accelerated by the pressure-driven transport and heat release.     

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.     

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.     

Strain rates,flow patterns and flame surface densities in hydrodynamically unstable,weakly turbulent premixed flames

Recent numerical and experimental studies have unveiled a potentially marked difference between the laminar as well as turbulent propagation of premixed flames exhibiting Darrieus–Landau (DL) (or hydrodynamic) instabilities from flames for which instabilities are inhibited. In this study we utilize two-dimensional numerical simulations of slot burner flames as well as experimental Propane–Air Bunsen flames to analyse differences in turbulent propagation, strain rate and induced flow patterns of hydrodynamically stable and unstable flames. We also investigate the effects of hydrodynamic instability on quantities which are directly related to reaction rate closure models, such as flame surface density and stretch factor. A clear enhancement of turbulent flame speed can be observed for unstable flames, generally mitigated at higher turbulence intensity, which is attributed to a flame area increase induced by the characteristic cusp-like DL-induced corrugation, absent in stable flames, which occurs concurrently and in synergy with turbulent wrinkling. Unstable flames also exhibit, both numerically and experimentally, a different correlation between strain rate and flame curvature and are observed to give rise to a channeling of the induced flow in the fresh mixture. Conditionally averaged flame surface density is also observed to attain smaller values in unstable flames, as a result of the thicker turbulent flame brush, indicating that closure models should incorporate instability-related parameters in addition to turbulence-related parameters.     

A DNS study of the physical mechanisms associated with density ratio influence on turbulent burning velocity in premixed flames

Data obtained in 3D direct numerical simulations of statistically planar, 1D weakly turbulent flames characterised by different density ratios σ are analysed to study the influence of thermal expansion on flame surface area and burning rate. Results show that, on the one hand, the pressure gradient induced within a flame brush owing to heat release in flamelets significantly accelerates the unburned gas that deeply intrudes into the combustion products in the form of an unburned mixture finger, thus causing large-scale oscillations of the burning rate and flame brush thickness. Under the conditions of the present simulations, the contribution of this mechanism to the creation of the flame surface area is substantial and is increased by σ, thus implying an increase in the burning rate by σ. On the other hand, the total flame surface areas simulated at σ = 7.53 and 2.5 are approximately equal. The apparent inconsistency between these results implies the existence of another thermal expansion effect that reduces the influence of σ on the flame surface area and burning rate. Investigation of the issue shows that the flow acceleration by the combustion-induced pressure gradient not only creates the flame surface area by pushing the finger tip into the products, but also mitigates wrinkling of the flame surface (the side surface of the finger) by turbulent eddies. The latter effect is attributed to the high-speed (at σ = 7.53) axial flow of the unburned gas, which is induced by the axial pressure gradient within the flame brush (and the finger). This axial flow acceleration reduces the residence time of a turbulent eddy in an unburned zone of the flame brush (e.g. within the finger). Therefore, the capability of the eddy for wrinkling the flamelet surface (e.g. the side finger surface) is weakened owing to a shorter residence time.     

The impact of dilatation,scrambling, and pressure transport in turbulent premixed flames

Premixed turbulent flames feature strong interactions between chemical reactions and turbulence that affect scalar and turbulence statistics. The focus of the present work is on clarifying the impact of pressure dilatation/flamelet scrambling effects with a comprehensive second-moment closure used for evaluation purposes. Model extensions that take into account flamelet orientation and molecular diffusion are derived. Isothermal pressure transport is included with an additional variable density contribution derived for the flamelet regime of combustion. Full closure is assessed by comparisons with Direct Numerical Simulations (DNSs) of statistically ‘steady’ fully developed premixed turbulent planar flames at different expansion ratios. Subsequently, the prediction of lean premixed turbulent methane–air flames featuring fractal grid generated turbulence in an opposed jet geometry is considered. The overall agreement shows that ‘dilatation’ effects contribute to counter-gradient transport and can also increase the turbulent kinetic energy significantly. Levels of anisotropy are broadly consistent with the DNS data and key aspects of opposed jet flames are well predicted. However, it is also shown that complications arise due to interactions between the imposed pressure gradient and combustion and that redistribution is affected along with the scalar flux at the leading edge. The latter is strongly affected by the reaction rate closure and, potentially, by pressure transport. Overall, the derived models offer significant improvements and can readily be applied to the modelling of premixed turbulent flames at practical rates of heat release.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

A DNS assessment of linear relations between filtered reaction rate,flame surface density,and scalar dissipation rate in a weakly turbulent premixed flame

Linear relations between (i) filtered reaction rate and filtered flame surface density (FSD) and (ii) filtered reaction rate and filtered scalar dissipation rate (SDR), which are widely used in Large Eddy Simulation (LES) research into premixed turbulent combustion, are examined by processing DNS data obtained from a statistically 1D planar flame under weakly turbulent conditions that are most favourable for the two approaches (flamelet combustion regime, single-step chemistry, equidiffusive mixture, adiabatic burner, and low Mach number). The analysis well supports the former approach provided that the filtered reaction rate is combined with filtered molecular transport term. In such a case, both the RANS and LES FSD approaches are based on local relations valid within weakly perturbed flamelets. Accordingly, simply recasting RANS expressions to a filtered form works well. On the contrary, while the FSD and SDR approaches appear to be basically similar at first glance, the analysis does not support the latter one, but shows that a ratio of the filtered reaction rate to the filtered SDR is strongly scattered within the studied flame brush, with its conditionally mean value varying significantly with Favre-filtered combustion progress variable. As argued in the paper, these limitations of the LES SDR approach stem from the fact that it is based on a relation valid after integration over weakly perturbed flamelets, but this relation does not hold locally within such flamelets. Consequently, when a sufficiently small filter is applied to instantaneous fields, the filter may contain only a part of the local flamelet, whereas the linear relation holds solely for the entire flamelet and may not hold within the filtered flamelet volume. Thus, the present study implies that straightforwardly recasting well-established RANS equations to a filtered form is a flawed approach if the equations are based on integral features of local burning.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

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.     

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.     

Effects of heat conduction and radical quenching on premixed stagnation flame stabilised by a wall

The premixed stagnation flame stabilised by a wall is analysed theoretically considering thermally sensitive intermediate kinetics. We consider the limit case of infinitely large activation energy of the chain-branching reaction, in which the radical is produced infinitely fast once the cross-over temperature is reached. Under the assumptions of potential flow field and constant density, the correlation for flame position and stretch rate of the premixed stagnation flame is derived. Based on this correlation, the effects of heat conduction and radical quenching on the wall surface are examined. The wall temperature is shown to have great impact on flame bifurcation and extinction, especially when the flame is close to the wall. Different flame structures are observed for near-wall normal flame, weak flame, and critically quenched flame. The fuel and radical Lewis numbers are found to have opposite effects on the extinction stretch rate. Moreover, it is also demonstrated that only when the flame is close to the wall does the radical quenching strongly influence the flame bifurcation and extinction. The extinction stretch rate is shown to decrease with the amount of radical quenching for different fuel and radical Lewis numbers. Besides, the coupling between the wall heat conduction and radical quenching is found to greatly influence the bifurcation and extinction of the premixed stagnation flame.     

Strain rate and flame orientation statistics in the near-wall region for turbulent flame-wall interaction

Flame–wall interaction (FWI) in premixed turbulent combustion has been analysed based on a counter-flow like configuration at the statistically stationary state. For the present configuration, the two FWI sub-zones, i.e the influence zone and the quenching zone, can be identified from the DNS results. Detailed analysis of the important quantities, such as the flame temperature, flame–wall distance, wall heat flux, flame curvature and dilatation (including the flame normal and tangential strain rates), and some orientation relations between the flame normal and the principal strain rate directions, have been reported, together with the physical explanations. All these statistical results are determined by the relative strengths of the wall heat flux, thermal expansion and the flame–turbulence interaction.     

Effects of subgrid scale and combustion modelling on flame structure of a turbulent premixed flame within LES and tabulated chemistry framework

The effects of combustion and SubGrid Scale (SGS) modelling on the overall flame characteristics of a turbulent premixed flame are investigated. This is achieved in terms of mean flow statistics, variances and flame surfaces. In particular, the chemical flame structure is analysed and compared. The Artificially Thickened Flame (ATF) approach coupled with the Flamelet Generated Manifolds (FGMs) and Filtered TAbulated Chemistry for LES (F-TACLES) approaches are used for this investigation. A Germano like procedure for dynamical calculation of SGS wrinkling is used which ensures the conservation of the total flame surface for both models. It turns out that using the dynamic SGS wrinkling model improves the results. Although the results of both combustion models in terms of statistics, mean and variances show very good agreement, the resolved flame surfaces hide different dynamic behaviour.