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.     

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.     

多孔障碍物对预混火焰传播的影响

以甲烷为代表性气体，研究了半封闭管道中设置多孔障碍物对可燃气体爆炸火焰传播的影响，基于大涡模拟对实验进行了重现，对比了实验与模拟中火焰传播过程的形状、位置及速度，分析了模拟结果中火焰穿过障碍物前后的流场和表面积变化，给出了衡量火焰褶皱程度的指标及算法。结果表明:大涡模拟结果与实验结果有较好的一致性；火焰在存在障碍物的管道内传播，经历层流快速膨胀、受阻回流、湍流快速发展和脉动减速4个阶段，各阶段火焰依次分别呈现加速、减速、二次加速、二次减速的波动变化；当可燃气体在开口与点火位置同端的管道内爆炸，火焰在接近障碍物时，受管道封闭端和障碍物约束显著，而出现脉动回流现象；火焰穿过多孔障碍物后，传播速度骤升至峰值，较未穿过障碍物前的最大速度可增加58.7%；障碍物是导致火焰面破碎以及面积褶皱率增大的直接原因，火焰褶皱率最大可达44.8%，比未穿过障碍物前的最大褶皱率增大39.27%。     

Experimental investigation of flame/solid interactions in turbulent premixed combustion

An experimental study has been carried out to investigate the interaction between propagating turbulent premixed flames and solid obstacles. The experimental rig was configured specifically to allow detailed measurements with laser-based optical diagnostics. A wall-type solid obstacle was mounted inside a laboratory-scale combustion chamber with rectangular cross-section. The flame was initiated, by igniting a combustible mixture of methane in air at the center of the closed end of the combustion chamber. The flame front development was visualized by a high-speed (9000 frame/s) digital video camera and flame images were synchronized with ignition timing and chamber pressure data. The tests were carried out with lean, stoichiometric and rich mixtures of methane in air. The images were used to calculate highly resolved temporal and spatial data for the changes in flame shape, speed, and the length of the flame front. The results are discussed in terms of the influence of mixture equivalence ratio on the flame structure and resulting overpressure. The reported data revealed significant changes in flame structure as a result of the interaction between the propagating flame front and the transient recirculating flow formed behind the solid obstacle. Combustion images show that the flame accelerates and decelerates as it impinges on the obstacle wall boundaries. It is also found that the mixture concentrations have a significant influence on the nature of the flame/solid interactions and the resulting overpressure. The highest flame speed of 40 m/s was obtained with the unity fuel–air equivalence ratio. Burning of trapped mixture behind the solid obstruction was found to be highly correlated with the flame front length and the rate of pressure rise.     

连续圆孔障碍物对油气泄压爆炸火焰特性影响大涡模拟

采用WALE模型和Zimont预混火焰模型对内置圆孔障碍物油气泄压爆炸火焰特性进行了大涡模拟,并将大涡模拟计算结果和RNG k-ε湍流模型计算结果以及实验结果进行对比分析,验证了大涡模拟的精确性。结果表明:（1）大涡模拟在预测油气爆炸超压、火焰传播速度以及火焰形态变化等方面比RNG k-ε湍流模型精确度更高,且能表现出更多流场的精细化结构;（2）障碍物诱导管道内形成湍流度较高的流场区域,导致火焰产生褶皱弯曲变形,增大火焰面积,加速火焰传播;（3）爆炸超压、火焰传播速度和火焰面积内在联系密切,具有显著的耦合性,且随时间的变化趋势存在高度的一致性。     

方孔障碍物对瓦斯火焰传播影响的实验与大涡模拟

为揭示置障管道内甲烷/空气预混火焰传播特性,运用高速摄影技术对甲烷/空气预混火焰的形状变化和火焰前锋的速度特性进行实验,并利用大涡模拟对管道内的流场结构进行数值分析。结果表明:置障管道内依次出现了球形火焰、指尖形火焰及“蘑菇”状火焰,且“蘑菇”状火焰出现之后,火焰开始反向传播;“蘑菇”状火焰是双涡旋结构与火焰前锋面相互作用的结果,而火焰的反向传播是由流场中出现逆流结构引起的;障碍物对火焰前锋有明显的加速作用;大涡模拟成功再现了实验中观察到的火焰形状、火焰前锋速度及流场结构,说明大涡模拟适用于置障管道内预混火焰传播特性的研究。     

Large Eddy Simulation of Diesel Fuel Injection and Mixing in a HSDI Engine

The paper aims to assess the performance of large eddy simulation (LES) in predicting the unsteady reacting flows in internal combustion engines. The incompatibility due to the turbulence dissipation was avoided in the k-equation LES formulation. Two versions of the LES have been tested with different filtering. The cell-specific filtering was found to give realistic prediction of the instantaneous temperature and pressure field during the combustion process. The coupling of combustion heat release, temperature field and turbulent flow field was found to be strong in the LES predicted flow and combustion fields which showed wrinkled flame structures. The formulation gives improved agreement with available experimental data.     

Large eddy simulation of turbulent premixed piloted flame using artificial thickened flame model coupled with tabulated chemistry

A sub-grid scale (SGS) combustion model, which combines the artificial thickened flame (ATF) model with the flamelet generated manifold (FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations (LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3. More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.     

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.     

Evaluations of SGS Combustion,Scalar Flux and Stress Models in a Turbulent Jet Premixed Flame

A newly developed fractal dynamic SGS (FDSGS) combustion model and a scale self-recognition mixed (SSRM) SGS stress model are evaluated along with other SGS combustion, scalar flux and stress models in a priori and a posteriori manners using DNS data of a hydrogen-air turbulent plane jet premixed flame. A posteriori tests reveal that the LES using the FDSGS combustion model can predict the combustion field well in terms of mean temperature distributions and peak positions in the transverse distributions of filtered reaction progress variable fluctuations. A priori and a posteriori tests of the scalar flux models show that a model proposed by Clark et al. accurately predicts the counter-gradient transport as well as the gradient diffusion, and introduction of the model of Clark et al. into the LES yields slightly better predictions of the filtered progress variable fluctuations than that of a gradient diffusion model. Evaluations of the stress models reveal that the LES with the SSRM model predicts the velocity fluctuations well compared to that with the Smagorinsky model.     

方形管内楔形障碍物对火焰结构与传播的影响

通过实验与数值模拟方法对CH4/空气预混火焰在有楔形障碍物的卧式燃烧方管内的传播进行了研究。采用多镜头Cranz Schardin高速摄像机和压力传感器等实验设备获得了高清晰度的障碍物诱导火焰失稳的分幅时序照片以及障碍物背风表面压力变化曲线。数值模拟则基于RANS方法与EDU-Arrhenius燃烧模型,计算结果与实验结果基本相符,反映了火焰在管内传播与变形的详细过程。通过综合分析实验与计算结果,得到了由楔形障碍物导致的火焰加速与变形的内在机理,揭示了火焰传播过程中由层流燃烧向湍流燃烧转捩的本质。     

LES Modeling of the Impact of Heat Losses and Differential Diffusion on Turbulent Stratified Flame Propagation: Application to the TU Darmstadt Stratified Flame

Common combustion chambers often exhibit turbulent flames propagating in partially-premixed mixtures. This propagation is generally governed by aerodynamics, unsteady mixing and chemical processes and may also be affected by conductive heat losses when the reactive zone develops close to the burner lips. The Filtered TAbulated Chemistry for Large Eddy Simulation (F-TACLES) model has been recently developed to include tabulated chemistry in Large Eddy Simulation (LES) of adiabatic stratified flames in flamelet regimes. The present article proposes a modeling approach to account for both differential diffusion and non-adiabatic effects on flame consumption speed following the F-TACLES formalism. The adiabatic F-TACLES model is first detailed using a generalized formalism for diffusive fluxes allowing either to account for differential diffusion or not. The F-TACLES model is then extended to non-adiabatic situations. A correction factor based on the non-adiabatic consumption rate is introduced to recover a realistic filtered flame consumption speed. The objective is here to tackle flame stabilization mechanisms when heat losses affect the reaction zone. The proposed approach is validated through the simulation of the unconfined stratified turbulent jet flame TSF-A for which stabilization process is affected by heat losses. Five simulations are performed for both adiabatic and non-adiabatic flow conditions comparing unity Lewis number and complex diffusion assumptions. The adiabatic F-TACLES model predicts a flame anchored at the burner lip disagreeing with experimental data. The non-adiabatic simulation exhibits local extinction due to heat losses near the burner exit. The flame is then lifted improving the comparison with experiments. Results also show a significant impact of molecular diffusion model on both mean flame consumption rate and angle.     

LES of Low to High Turbulent Combustion in an Elevated Pressure Environment

A subgrid scale flame surface density combustion model for the Large Eddy Simulation (LES) of premixed combustion is derived and validated. The model is based on fractal characteristics of the flame surface, assuming a self similar wrinkling of the flame between smallest and largest wrinkling length scales. Experimental and direct numerical simulation databases as well as theoretical models are used to derive a model for the fractal parameters, namely the cut-off lengths and the fractal dimension suitable in the LES context. The combustion model is designed with the intent to simulate low as well as high Reynolds number premixed turbulent flame propagation and with a focus on correct scaling with pressure. The combustion model is validated by simulations of turbulent Bunsen flames with methane and propane fuel at pressure levels between 0.1 MPa and 2 MPa and at turbulence levels of $0 < u^{\prime }/s_{L}^{0} < 11$ , conditions typical for spark ignition engines. The predicted turbulent flame speed is in a very good agreement with the experimental data and a smooth transition from resolved flame wrinkling to fully modelled, nearly subgrid-only wrinkling is realized. Evaluating the influence of mesh resolution shows a predicted mean flame surface and turbulent flame speed independent of mesh resolution for cases with 9–86 % resolved flame surface. Additional simulations of a highly turbulent jet flame at 0.1 MPa and 0.5 MPa and the comparison with experimental data in terms of flame shape, velocity field and turbulent fluctuations validates the model also at conditions typical for gas turbines.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

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.     

Experimental investigation on micro-dynamic behavior of gas explosion suppression with SiO2 fine powders

To study the effect of inert dust on gas explosion suppression mechanism, SiO₂ fine powders were sprayed to suppress premixed CH₄-Air gas explosion in a 20 L spherical experimental system. In the experiment, high speed schlieren image system was adopted to record explosion flame propagation behaviors, meanwhile, pressure transducers and ion current probes were used to clearly record the explosion flame dynamic characteristics. The experimental results show that the SiO₂ fine powders suppressed evidently the gas explosion flame, and reduced the peak value of pressure and flame speed by more than 40 %. The ion current result shows that the SiO₂ super fine powders were easy to contact with and absorb free radicals near the combustion reaction region, which greatly reduced the combustion reaction intensity, and in turn influenced the flame propagation and pressure rising.     

The propagation mechanism of high speed turbulent deflagrations

The propagation regimes of combustion waves in a 30 cm by 30 cm square cross–sectioned tube with an obstacle array of staggered vertical cylindrical rods (with BR=0.41 and BR=0.19) are investigated. Mixtures of hydrogen, ethylene, propane, and methane with air at ambient conditions over a range of equivalence ratios are used. In contrast to the previous results obtained in circular cross–sectioned tubes, it is found that only the quasi–detonation regime and the slow turbulent deflagration regimes are observed for ethylene–air and for propane–air. The transition from the quasi–detonation regime to the slow turbulent deflagration regime occurs at (where D is the tube “diameter” and is the detonation cell size). When , the quasi–detonation velocities that are observed are similar to those in unobstructed smooth tubes. For hydrogen–air mixtures, it is found that there is a gradual transition from the quasi–detonation regime to the high speed turbulent deflagration regime. The high speed turbulent deflagration regime is also observed for methane–air mixtures near stoichiometric composition. This regime was previously interpreted as the “choking” regime in circular tubes with orifice plate obstacles. Presently, it is proposed that the propagation mechanism of these high speed turbulent deflagrations is similar to that of Chapman–Jouguet detonations and quasi-detonations. As well, it is observed that there exists unstable flame propagation at the lean limit where . The local velocity fluctuates significantly about an averaged velocity for hydrogen–air, ethylene–air, and propane–air mixtures. Unstable flame propagation is also observed for the entire range of high speed turbulent deflagrations in methane–air mixtures. It is proposed that these fluctuations are due to quenching of the combustion front due to turbulent mixing. Quenched pockets of unburned reactants are swept downstream, and the subsequent explosion serves to overdrive the combustion front. The present study indicates that the dependence on the propagation mechanisms on obstacle geometry can be exploited to elucidate the different complex mechanisms of supersonic combustion waves. Received 5 November 2001 / Accepted 12 June 2002 / Published online 4 November 2002 Correspondence to: J. Chao (e-mail: jenny.chao@mail.mcgill.ca) An abridged version of this paper was presented at the 18th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Seattle, USA, from July 29 to August 3, 2001.     

Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method

Large eddy simulations (LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density (FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds (FGM) method is employed to determine major scalars. This coupled sub-grid scale (SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence. The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.