Large Eddy Simulation of Turbulent Premixed Swirling Flames Using Dynamic Thickened Flame with Tabulated Detailed Chemistry

A sub-grid scale (SGS) combustion model by combining dynamic thickened flame (DTF) with flamelet generated manifolds (FGM) tabulation approach (i.e. DTF-FGM) is developed for investigating turbulent premixed combustion. In contrast to the thickened flame model, the dynamic thickening factor of the DTF model is determined from the flame sensor, which is obtained from the normalized gradient of the reaction progress variable from the one-dimensional freely propagating premixed flame simulations. Therewith the DTF model can ensure that the thickening of the flame is limited to the regions where it is numerically necessary. To describe the thermo-chemistry states, large eddy simulation (LES) transport equations for two characteristic scalars (the mixture fraction and the reaction progress variable) and relevant sub-grid variances in the DTF-FGM model are presented. As to the evaluation of different SGS combustion models, another model by utilizing the combination of presumed probability density function (PPDF) and FGM (i.e. PPDF-FGM) is also described. LES of two cases with or without swirl in premixed regime of the Cambridge swirl burner flames are performed to evaluate the developed SGS combustion model. The predicted results are compared with the experimental data in terms of the influence of different LES grids, model sensitivities to the thickening factor, the wrinkling factor, and the PPDF of characteristic scalars, the evaluation of different modelling approaches for the sub-grid variances of characteristic scalars, and the predictive capability of different SGS combustion models. It is shown that the LES results with the DTF-FGM model are in reasonable agreement with the experimental data, and better than the results with the PPDF-FGM approach due to its ability to predict better in regions where flame is not resolved.     

当量比对预混燃烧影响的大涡模拟研究

本文对在突扩燃烧室内甲烷和空气的预混燃烧进行了大涡模拟(LES)研究,考虑预混燃料的当量比对燃烧室提供的动力及产生的污染物的影响．利用LES计算了不同当量比条件下燃烧室内湍流预混燃烧反应流场的温度、浓度、涡量和压力分布,最后对当量比0.5时B点和C点的温度和速度进行EMD分解,得到了温度场和速度场的各阶模态的平均周期．结果表明：随着当量比从0.5增加至0.7,燃烧反应趋于剧烈,燃烧室的最高温度提高了350K,平均压力从32.876 Pa增大到34.833Pa,燃烧产生的瞬态径向最高浓度从0.5%增加到0.95%．     

Assessment of Finite Rate Chemistry Large Eddy Simulation Combustion Models

Large Eddy Simulations (LES) of a swirl-stabilized natural gas-air flame in a laboratory gas turbine combustor is performed using six different LES combustion models to provide a head-to-head comparative study. More specifically, six finite rate chemistry models, including the thickened flame model, the partially stirred reactor model, the approximate deconvolution model and the stochastic fields model have been studied. The LES predictions are compared against experimental data including velocity, temperature and major species concentrations measured using Particle Image Velocimetry (PIV), OH Planar Laser-Induced Fluorescence (OH-PLIF), OH chemiluminescence imaging and one-dimensional laser Raman scattering. Based on previous results a skeletal methane-air reaction mechanism based on the well-known Smooke and Giovangigli mechanism was used in this work. Two computational grids of about 7 and 56 million cells, respectively, are used to quantify the influence of grid resolution. The overall flow and flame structures appear similar for all LES combustion models studied and agree well with experimental still and video images. Takeno flame index and chemical explosives mode analysis suggest that the flame is premixed and resides within the thin reaction zone. The LES results show good agreement with the experimental data for the axial velocity, temperature and major species, but differences due to the choice of LES combustion model are observed and discussed. Furthermore, the intrinsic flame structure and the flame dynamics are similarly predicted by all LES combustion models examined. Within this range of models, there is no strong case for deciding which model performs the best.     

Subgrid Scale Modeling in Large-Eddy Simulation of Turbulent Combustion Using Premixed Flamelet Chemistry

Large-eddy simulation (LES) of turbulent combustion with premixed flamelets is investigated in this paper. The approach solves the filtered Navier–Stokes equations supplemented with two transport equations, one for the mixture fraction and another for a progress variable. The LES premixed flamelet approach is tested for two flows: a premixed preheated Bunsen flame and a partially premixed diffusion flame (Sandia Flame D). In the first case, we compare the LES with a direct numerical simulation (DNS). Four non-trivial models for the chemical source term are considered for the Bunsen flame: the standard presumed beta-pdf model, and three new propositions (simpler than the beta-pdf model): the filtered flamelet model, the shift-filter model and the shift-inversion model. A priori and a posteriori tests are performed for these subgrid reaction models. In the present preheated Bunsen flame, the filtered flamelet model gives the best results in a priori tests. The LES tests for the Bunsen flame are limited to a case in which the filter width is only slightly larger than the flame thickness. According to the a posteriori tests the three models (beta-pdf, filtered flamelet and shift-inversion) show more or less the same results as the trivial model, in which subgrid reaction effects are ignored, while the shift-filter model leads to worse results. Since LES needs to resolve the large turbulent eddies, the LES filter width is bounded by a maximum. For the present Bunsen flame this means that the filter width should be of the order of the flame thickness or smaller. In this regime, the effects of subgrid reaction and subgrid flame wrinkling turn out to be quite modest. The LES-results of the second case (Sandia Flame D) are compared to experimental data. Satisfactory agreement is obtained for the main species. Comparison is made between different eddy-viscosity models for the subgrid turbulence, and the Smagorinsky eddy-viscosity is found to give worse results than eddy-viscosities that are not dominated by the mean shear. Paper presented on the Eccomas Thematic Conference Computational Combustion 2007, submitted for a special issue of Flow, Turbulence and Combustion.     

Large Eddy Simulation of a Propagating Turbulent Premixed Flame

A large eddy simulation of a turbulent premixed flame propagatingthrough a chamber containing a square obstruction is presented anddiscussed. The governing equations for compressible, reacting flowsare Favre filtered and turbulence closure is achieved using thedynamic Smagorinsky subgrid model. A simple flame surface densitymodel based on the flamelet concept is employed for the subgrid scalereaction rate. The simulation gives very good agreement with experimentalresults for the speed and the shape of the flame as it propagates throughthe chamber. The peak pressures, however, are underpredicted by20–30%. Reasons for this are discussed and it is concluded that amore sophisticated combustion model is required for complex flowssuch as this one, if a more accurate prediction of the pressureis to be achieved.     

Large Eddy Simulation of Premixed Turbulent Flames Using the Probability Density Function Approach

In the present study a Large Eddy Simulation and Filtered Density Function model is applied to three premixed piloted turbulent methane flames at different Reynolds Numbers using the Eulerian stochastic fields approach. The model is able to reproduce the flame structure and flow characteristics with a low number of fields (between 4 and 16 fields). The results show a good agreement with experimental data with the same closures employed in non-premixed combustion without any adjustment for combustion regime. The effect of heat release on the flow field is captured correctly. A wide range of sensitivity studies is carried out, including the number of fields, the chemical mechanism, differential diffusion effects and micro-mixing closures. The present work shows that premixed combustion (at least in the conditions under study) can be modelled using LES-PDF methods.. Finally, the ability of the model to predict flame quenching is studied. The model can accurate capture the conditions at which combustion is not sustainable and large pockets of extinction appear.     

Large Eddy Simulation of a Bluff Body Stabilised Premixed Flame Using Flamelets

Large Eddy Simulations of an unconfined turbulent lean premixed flame, which is stabilised behind a bluff body, are conducted using unstrained flamelets as the sub-grid scale combustion closure. The statistics from the simulations are compared with the corresponding data obtained from the experiment and it is demonstrated that the experimental observations are well captured. The relative positioning of the shear layers and the flame brush are analysed to understand the radial variations of the turbulent kinetic energy at various streamwise locations. These results are also compared to confined bluff body stabilised flames, to shed light on the relative role of incoming and shear driven turbulence on the behaviour of the flame brush and the turbulent kinetic energy variation across it.     

Scalar Dissipation Rate Transport in the Context of Large Eddy Simulations for Turbulent Premixed Flames with Non-Unity Lewis Number

The effects of global Lewis number Le on the statistical behaviour of the unclosed terms in the transport equation of the Favre-filtered scalar dissipation rate (SDR) Ñ _c have been analysed using a Direct Numerical Simulation (DNS) database of freely propagating statistically planer turbulent premixed flames with Le ranging from 0.34 to 1.2. The DNS data has been explicitly filtered to analyse the statistical behaviour of the unclosed terms in the SDR transport equation arising from turbulent transport T ₁, density variation due to heat release T ₂, scalar-turbulence interaction T ₃, reaction rate gradient T ₄, molecular dissipation (?D ₂) and diffusivity gradients f(D) in the context of Large Eddy Simulations (LES). It Le has significant effects on the magnitudes of T ₁, T ₂, T ₃, T ₄, (?D ₂) and f(D). Moreover, both qualitative and quantitative behaviours of the unclosed terms T ₁, T ₂, T ₃, T ₄, (?D ₂) and f(D) are found to be significantly affected by the LES filter width Δ, which have been explained based on a detailed scaling analysis. Both scaling analysis and DNS data suggest that T ₂, T ₃, T ₄, (?D ₂) and f(D) remain leading order contributors to the SDR \(\tilde {{N}}_{c} \) transport for LES. The scaling estimates of leading order contributors to the SDR \(\tilde {{N}}_{c} \) transport has been utilised to discuss the possibility of extending an existing SDR model for Reynolds Averaged Navier Stokes (RANS) simulation for SDR \(\tilde {{N}}_{c} \) closure in the context of LES of turbulent premixed combustion.     

Assessment of the Dynamic SGS Wrinkling Combustion Modeling Using the Thickened Flame Approach Coupled with FGM Tabulated Detailed Chemistry

The dynamic power-law wrinkling model proposed by Charlette et al. is coupled with Flamelet Generated Manifolds (FGM) tabulated chemistry combined with an artificially thickened flame model (ATF) for large eddy simulation. The dynamic formulation is similar to the “Germano” procedure and uses Taylor series based Gaussian filter. Thereby, the power-law wrinkling model parameter is considered to have both temporal and spatial dependency. Series of simulations are conducted for a lean premixed turbulent flame, using both dynamic and non-dynamic versions of the wrinkling model under different grid levels. The simulation results applying the non-dynamic wrinkling model show different behavior for each particular flame resolution, where none of the simulations could deliver the correct flame statistics, such as flame height. The dynamic version of the power-law wrinkling model improves the results independently of the flame resolution, as a consequence of the conservation of the total flame surface.     

Large Eddy Simulation of a Novel Gas-Assisted Coal Combustion Chamber

In this work a recently presented combustion chamber that is specifically designed for the investigation of gas-assisted coal combustion and the validation of models is simulated under reactive conditions for the first time. In the configuration coal combustion is assisted and stabilized by a methane flame. In the course of the investigation, the configuration’s complexity is increased successively. Results of the isothermal single-phase flow are discussed first. Subsequently, reproducibility of the single-phase methane flame by means of the applied modeling approach is evaluated. In a further step, coal particles having the same thermal power as the methane flame are injected into the configuration. Particle histories, the conversion of the coal particles as well as its retroactive effect on the gas phase are investigated. Experimental results based on laser diagnostics are provided for all operating points and used for comparison with numerical results. Gas phase velocity fields for all operating points are available. In order to identify the reaction in the reactive single-phase case planar laser induced fluorescence of the OH-radical (OH-PLIF) was applied. Overall good agreement between numerical and experimental results could be obtained. In the Large Eddy Simulation (LES) a Flamelet Generated Manifold (FGM) based model is utilized. The four-dimensional manifold is spanned by two mixture fractions, a reaction progress variable and the enthalpy on which the gas phase chemistry gets mapped onto. Thereby, the model accounts for both, volatiles reaction and char conversion. Furthermore, finite rate chemistry effects as well as non-adiabatic physics are considered.     

Large Eddy Simulation and PIV Measurements of Unsteady Premixed Flames Accelerated by Obstacles

In gas explosions, the unsteady coupling of the propagating flame and the flow field induced by the presence of blockages along the flame path produces vortices of different scales ahead of the flame front. The resulting flame–vortex interaction intensifies the rate of flame propagation and the pressure rise. In this paper, a joint numerical and experimental study of unsteady premixed flame propagation around three sequential obstacles in a small-scale vented explosion chamber is presented. The modeling work is carried out utilizing large eddy simulation (LES). In the experimental work, previous results (Patel et al., Proc Combust Inst 29:1849–1854, 2002) are extended to include simultaneous flame and particle image velocimetry (PIV) measurements of the flow field within the wake of each obstacle. Comparisons between LES predictions and experimental data show a satisfactory agreement in terms of shape of the propagating flame, flame arrival times, spatial profile of the flame speed, pressure time history, and velocity vector fields. Computations through the validated model are also performed to evaluate the effects of both large-scale and sub-grid scale (SGS) vortices on the flame propagation. The results obtained demonstrate that the large vortical structures dictate the evolution of the flame in qualitative terms (shape and structure of the flame, succession of the combustion regimes along the path, acceleration-deceleration step around each obstacle, and pressure time trend). Conversely, the SGS vortices do not affect the qualitative trends. However, it is essential to model their effects on the combustion rate to achieve quantitative predictions for the flame speed and the pressure peak.     

Large Eddy Simulation of Stratified and Sheared Flames of a Premixed Turbulent Stratified Flame Burner Using a Flamelet Model with Heat Loss

This paper presents large eddy simulations (LES) of the Darmstadt turbulent stratified flame burner (TSF) at different operating conditions including detailed heat loss modeling. The target cases are a non-reacting and two reacting cases. Both reacting cases are characterized by stratification, while one flame additionally features shear. In the regime diagram for premixed combustion, the studied flames are found at the border separating the thin reaction zones regime and the broken reaction zones regime. A coupled level set/progress variable model is utilized to describe the combustion process. To account for heat loss, an enthalpy defect approach is adopted and reformulated to include differential diffusion effects. A novel power-law rescaling methodology is proposed to integrate the enthalpy defect approach into the level set/progress variable model which is extensively validated in two validation scenarios. It is demonstrated that the LES with the newly developed model captures the influence of heat loss well and that the incorporation of heat loss effects improves the predictions of the TSF-burner over adiabatic simulations, while reproducing the experimentally observed flame lift-off from the pilot nozzle.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

Investigation of Two-Fluid Methods for Large Eddy Simulation of Spray Combustion in Gas Turbines

An extension of the large eddy simulation (LES) technique to two-phase reacting flows, required to capture and predict the behavior of industrial burners, is presented. While most efforts reported in the literature to construct LES solvers for two-phase flow focus on Euler–Lagrange formulation, the present work explores a different solution (‘two-fluid’ approach) where an Eulerian formulation is used for the liquid phase and coupled with the LES solver of the gas phase. The equations used for each phase and the coupling terms are presented before describing validation in two simple cases which gather some of the specificities of real combustion chamber: (1) a one-dimensional laminar JP10/air flame and (2) a non-reacting swirled flow where solid particles disperse (Sommerfeld and Qiu, Int. J. Multiphase Flow 19(6):1093–1127, 1993). After these validations, the LES tool is applied to a realistic aircraft combustion chamber to study both a steady flame regime and an ignition sequence by a spark. Results bring new insights into the physics of these complex flames and demonstrate the capabilities of two-fluid LES.     

3D Numerical Simulation of a Laminar Experimental SWQ Burner with Tabulated Chemistry

Flame-wall interaction (FWI) plays an important role in enclosed combustion systems. For avoiding the complexity of close to reality combustors, in this study an atmospheric premixed V-shaped flame interacting with an isothermal cold wall in a side wall quenching (SWQ) configuration is investigated. A stoichiometric methane/air mixture is used as fuel. A three-dimensional (3D) numerical simulation, which resolves all flow structures is combined with a tabulated chemistry approach (flamelet generated manifold, FGM). Results are compared with experimental data and two-dimensional simulations. The FGM approach is a suitable trade-off between computationally expensive detailed chemistry simulations and over simplified single step mechanisms. 2D simulations are used to investigate the influence of the uncertainty of the wall temperature, to show that the resolution in 3D is sufficient and that the influence of the flame thickening on the wall heat fluxes can be determined. Our results show that the 3D FGM approach is in close agreement to experimentally obtained flow and temperature fields. The dimensionless wall heat flux and Péclet number matches the expected values of 0.16 and 7, respectively. However, during FWI the measured CO mole fractions are not reproduced accurately showing that the transported variables in the present approach of tabulated chemistry do not recover premixed flame structures near walls.