Modeling CO With Flamelet-Generated Manifolds: Part 2 — Application

Author(s):  
Graham Goldin ◽  
Zhuyin Ren ◽  
Hendrik Forkel ◽  
Liuyan Lu ◽  
Venkat Tangirala ◽  
...  

Conventional Flamelet Generated Manifold (FGM) closure of the mean progress variable reaction rate assumes PDF shapes to account for turbulent fluctuations. The FGM parameters are commonly assumed to be statistically independent, and the marginal PDFs invariably require second moments, which are difficult to model accurately and have limited coefficients that can be adjusted to calibrate the simulation. A new model is presented which locates the flame brush with a turbulent flame speed model, and applies the FGM kinetic rate to model kinetically limited processes, such as CO quenching, behind the flame-front. The model is applied to 3D RANS simulations of an equivalence ratio sweep in the GE Entitlement Rig perfectly premixed combustor experiment. Calculating the mean FGM reaction progress source term with standard assumed shape PDFs leads to a narrow flame brush and equilibrium CO outlet emissions. By limiting the mean FGM reaction progress source term by the turbulent flame speed model, the flame brush is broadened and super-equilibrium CO is predicted at the outlet. Good agreement with measurement is obtained with default model coefficients. Since the majority of the mean reaction progress source term is limited by the turbulent flame speed reaction rate, it is demonstrated that the model is relatively insensitive to assumed shape PDFs for the FGM rate, as well as the parameter used to determine the turbulent flame leading edge.

Author(s):  
George Mallouppas ◽  
Graham Goldin ◽  
Yongzhe Zhang ◽  
Piyush Thakre ◽  
Jim Rogerson

Abstract Three Flamelet Generated Manifold reaction source term closure options and two different reactor types are examined with Large Eddy Simulation of an industrial gas turbine combustor operating at 3 bar. This work presents the results for the SGT-100 Dry Low Emission (DLE) gas turbine provided by Siemens Industrial Turbomachinery Ltd. The related experimental study was performed at the German Aerospace Centre, DLR, Stuttgart, Germany. The FGM model approximates the thermo-chemistry in a turbulent flame as that in a simple 0D constant pressure ignition reactors and 1D strained opposed-flow premixed reactors, parametrized by mixture fraction, progress variable, enthalpy and pressure. The first objective of this work is to compare the flame shape and position predicted by these two FGM reactor types. The Kinetic Rate (KR) model, studied in this work, uses the chemical rate from the FGM with assumed shapes, which are a Beta function for mixture fraction and delta functions for reaction progress variable and enthalpy. Another model investigated is the Turbulent Flame-Speed Closure (TFC) model with Zimont turbulent flame speed, which propagates premixed flame fronts at specified turbulent flame speeds. The Thickened Flame Model (TFM), which artificially thickens the flame to sufficiently resolve the internal flame structure on the computational grid, is also explored. Therefore, a second objective of this paper is to compare KR, TFC and TFM with the available experimental data.


Author(s):  
Graham Goldin ◽  
Federico Montanari ◽  
Sunil Patil

LES and RANS simulations of a Siemens scaled combustor are compared against comprehensive experimental data. The steady RANS simulation modeled one quarter of the geometry with 8M polyhedral cells using the SST-k-ω model. Unsteady LES simulations were performed on the quarter geometry (90°, 8M cells) as well as the full geometry (360°, 32M cells) using the WALE sub-grid model and dynamic evaluation of model coefficients. Aside from the turbulence model, all other models are identical for the RANS and LES. Combustion was modeled with the Flamelet Generated Manifold (FGM) model, which represents the thermo-chemistry by mixture fraction and reaction progress. RANS simulations are performed using Zimont and Peters turbulent flame speed (TFS) expressions with default model constants, as well as the kinetic rate from the FGM. The flame speed stalls near the wall with the TFS models, predicting a flame brush that extends to the combustor outlet, which is inconsistent with measurements. The FGM kinetic source model shows improved flame position predictions. The LES predictions of mean and rms axial velocity, mixture fraction and temperature do not show improvement over the RANS. All three simulations over-predict the turbulent mixing in the inner recirculation zone, causing flatter profiles than measurements. This over-mixing is exacerbated in the 900 case. The experiments show evidence of heat loss and the adiabatic simulations presented here might be improved by including wall heat-loss and radiation effects.


Author(s):  
Sunil Patil ◽  
Federico Montanari

Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulations (LES) of a Siemens scaled combustor are compared against comprehensive experimental data. The steady RANS simulation modeled one quarter of the geometry with 8 M polyhedral cells using the shear stress transport (SST) k-ω model. Unsteady LES were performed on the quarter geometry (90 deg, 8 M cells) as well as the full geometry (360 deg, 32 M cells) using the wall-adapting local eddy-viscosity (WALE) subgrid model and dynamic evaluation of model coefficients. Aside from the turbulence model, all other models are identical for the RANS and LES. Combustion was modeled with the flamelet generated manifold (FGM) model, which represents the thermochemistry by mixture fraction and reaction progress. RANS simulations are performed using Zimont and Peters turbulent flame-speed (TFS) expressions with default model constants, as well as the kinetic rate from the FGM. The flame-speed stalls near the wall with the TFS models, predicting a flame brush that extends to the combustor outlet, which is inconsistent with measurements. The FGM kinetic source model shows improved flame position predictions. The LES predictions of mean and rms axial velocity, mixture fraction, and temperature do not show improvement over the RANS. All three simulations overpredict the turbulent mixing in the inner recirculation zone, causing flatter profiles than measurements. This overmixing is exacerbated in the 90 deg case. The experiments show evidence of heat loss, and the adiabatic simulations presented here might be improved by including wall heat-loss and radiation effects.


2004 ◽  
Vol 126 (4) ◽  
pp. 701-707 ◽  
Author(s):  
Ulf Engdar ◽  
Per Nilsson ◽  
Jens Klingmann

Most of the common modeling approaches to premixed combustion in engineering applications are either based on the assumption of infinitely fast chemistry or the flamelet assumption with simple chemistry. The level-set flamelet library approach (FLA) has shown great potential in predicting major species and heat release, as well as intermediate and minor species, where more simple models often fail. In this approach, the mean flame surface is tracked by a level-set equation. The flamelet libraries are generated by an external code, which employs a detailed chemical mechanism. However, a model for the turbulent flame speed is required, which, among other considerations, depends on the turbulence intensity, i.e., these models may show sensitivity to turbulence modeling. In this paper, the FLA model was implemented in the commercial CFD program Star-Cd, and applied to a lean premixed flame stabilized by a triangular prism (bluff body). The objective of this paper has been to investigate the impact on the mean flame position, and hence on the temperature and species distribution, using three different turbulent flame speed models in combination with four different turbulence models. The turbulence models investigated are: the standard k-ε model, a cubic nonlinear k-ε model, the standard k-ω model and the shear stress transport (SST) k-ω model. In general, the computed results agree well with experimental data for all computed cases, although the turbulence intensity is strongly underestimated at the downstream position. The use of the nonlinear k-ε model offers no advantage over the standard model, regardless of flame speed model. The k-ω based turbulence models predict the highest turbulence intensity with the shortest flame lengths as a consequence. The Mu¨ller flame speed model shows the least sensitivity to the choice of turbulence model.  


Author(s):  
Ulf Engdar ◽  
Per Nilsson ◽  
Jens Klingmann

Most of the common modeling approaches to premixed combustion in engineering applications are either based on the assumption of infinitely fast chemistry or the flamelet assumption with simple chemistry. The level-set flamelet library approach (FLA) has shown great potential in predicting major species and heat release, as well as intermediate and minor species, where more simple models often fail. In this approach, the mean flame surface is tracked by a level-set equation. The flamelet libraries are generated by an external code, which employs a detailed chemical mechanism. However, a model for the turbulent flame speed is required, which, amongst other considerations, depends on the turbulence intensity, i.e. these models may show sensitivity to turbulence modeling. In this paper, the FLA model was implemented in the commercial CFD program Star-CD, and applied to a lean premixed flame stabilized by a triangular prism (bluff body). The objective of this paper has been to investigate the impact on the mean flame position, and hence on the temperature and species distribution, using three different turbulent flame speed models in combination with four different turbulence models. The turbulence models investigated are: the standard k-ε model, a cubic non-linear k-ε model, the standard k-ω model and the Shear Stress Transport (SST) k-ω model. In general, the computed results agree well with experimental data for all computed cases, although the turbulence intensity is strongly underestimated at the downstream position. The use of the non-linear k-ε model offers no advantage over the standard model, regardless of flame speed model. The k-ω based turbulence models predict the highest turbulence intensity with the shortest flame lengths as a consequence. The Mu¨ller flame speed model shows the least sensitivity to the choice of turbulence model.


2000 ◽  
Vol 123 (1) ◽  
pp. 175-181 ◽  
Author(s):  
D. S.-K. Ting ◽  
M. D. Checkel

This paper presents a flame growth model based on experimental measurements of flame speed and mean turbulent strain rate. Methane/air mixtures of 0.7 and 0.9 equivalence ratios were centrally spark-ignited in a 125 mm cubical chamber. Based on schlieren images and combustion pressure traces, a linear correlation was found between the turbulent flame speed and the turbulent strain rate. For these unity-Lewis-number and near-zero-Markstein-number flames, the effectiveness of turbulent strain in enhancing the flame speed was found to increase linearly with the mean flame radius over the range of conditions tested.


2012 ◽  
Vol 709 ◽  
pp. 191-222 ◽  
Author(s):  
T. D. Dunstan ◽  
N. Swaminathan ◽  
K. N. C. Bray

AbstractThe sensitivity of the turbulent flame speed to the geometry of the flame is investigated using direct numerical simulations of turbulent premixed flames in three canonical configurations: freely propagating statistically planar flames, planar flames stabilized in stagnating flows, and rod-stabilized V-flames. We consider both the consumption speed, which measures the integrated rate of burning, and the propagation speed, which measures the speed of an isosurface within the flame brush. An algebraic model for the propagation speed of the leading edge of the flame brush, which is blind to flame geometry, is also applied to the data for the purposes of establishing its range of validity and the causes of its failure. The turbulent consumption speed is found to be strongly geometry dependent, primarily due to the continuous growth of the flame brush thickness. Changes in the structure and consumption speed of instantaneous flame fronts are found to be only weakly sensitive to flame geometry. The turbulent propagation speed is analysed in terms of its reactive, diffusive and turbulent flux components. All three terms are shown to be significant, both through the flame brush and along the leading edge. The leading-edge propagation speed is found to be sensitive to flame geometry only in the V-flames under certain conditions. It is suggested that this apparent geometry dependence, which the model cannot capture, results from the relation between the turbulence and mean flow time scales in these particular cases, and is not intrinsic to the flame geometry itself.


1997 ◽  
Vol 353 ◽  
pp. 83-114 ◽  
Author(s):  
DENIS VEYNANTE ◽  
THIERRY POINSOT

In most practical situations, turbulent premixed flames are ducted and, accordingly, subjected to externally imposed pressure gradients. These pressure gradients may induce strong modifications of the turbulent flame structure because of buoyancy effects between heavy cold fresh and light hot burnt gases. In the present work, the influence of a constant acceleration, inducing large pressure gradients, on a premixed turbulent flame is studied using direct numerical simulations.A favourable pressure gradient, i.e. a pressure decrease from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. It also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. As proposed by Libby (1989), the turbulent flame speed is modified by a buoyancy term linearly dependent on both the imposed pressure gradient and the integral length scale lt.A simple model for the turbulent flux u″c″ is also proposed, validated from simulation data and compared to existing models. It is shown that turbulent premixed flames can exhibit both gradient and counter-gradient transport and a criterion integrating the effects of pressure gradients is derived to differentiate between these regimes. In fact, counter-gradient diffusion may occur in most practical ducted flames.


Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 146 ◽  
Author(s):  
Aaron Endres ◽  
Thomas Sattelmayer

Boundary layer flashback from the combustion chamber into the premixing section is a threat associated with the premixed combustion of hydrogen-containing fuels in gas turbines. In this study, the effect of pressure on the confined flashback behaviour of hydrogen-air flames was investigated numerically. This was done by means of large eddy simulations with finite rate chemistry as well as detailed chemical kinetics and diffusion models at pressures between 0 . 5 and 3 . It was found that the flashback propensity increases with increasing pressure. The separation zone size and the turbulent flame speed at flashback conditions decrease with increasing pressure, which decreases flashback propensity. At the same time the quenching distance decreases with increasing pressure, which increases flashback propensity. It is not possible to predict the occurrence of boundary layer flashback based on the turbulent flame speed or the ratio of separation zone size to quenching distance alone. Instead the interaction of all effects has to be accounted for when modelling boundary layer flashback. It was further found that the pressure rise ahead of the flame cannot be approximated by one-dimensional analyses and that the assumptions of the boundary layer theory are not satisfied during confined boundary layer flashback.


Sign in / Sign up

Export Citation Format

Share Document