Large Eddy Simulation of ALSTOM’s Reheat Combustor Using Tabulated Chemistry and Stochastic Fields-Combustion Model

Author(s):  
Rohit Kulkarni ◽  
Birute Bunkute ◽  
Fernando Biagioli ◽  
Michael Duesing ◽  
Wolfgang Polifke

Large Eddy Simulations (LES) of natural gas ignition and combustion in turbulent flows are performed using a novel combustion model based on a composite progress variable, a tabulated chemistry ansatz and the stochastic-fields turbulence-chemistry interaction model. It is a significant advantage of this approach that it can be applied to industrial configurations with multi-stream mixing at relatively low computational cost and modeling complexity. The computational cost is independent of the chemical mechanism or the type of fuel, but increases linearly with the number of streams. The model is validated successfully against the Cabra methane flame and Delft Jet in Hot Coflow (DJFC) flame. Both cases constitute fuel jets in a vitiated coflow. The DJFC flame coflow has a non-uniform mixture of air and hot gases. The model considers this non-uniformity by an additional mixture fraction dimension, emulating a ternary mixing case. The model not only predicts flame location, but also the temperature distribution quantitatively. The LES combustion model is further extended to consider four stream mixing. It has been successfully validated for ALSTOM’s reheat combustor at atmospheric conditions. Compared to the past steady-state RANS (Reynolds Averaged Navier-Stokes) simulations [1], the LES simulations provide an even better understanding of the turbulent flame characteristics, which helps in the burner optimization.

2019 ◽  
Vol 141 (11) ◽  
Author(s):  
Alexander Avdonin ◽  
Alireza Javareshkian ◽  
Wolfgang Polifke

Abstract This paper demonstrates that a large Eddy simulation (LES) combustion model based on tabulated chemistry and Eulerian stochastic fields can successfully describe the flame dynamics of a premixed turbulent swirl flame. The combustion chemistry is tabulated from one-dimensional burner-stabilized flamelet computations in dependence on progress variable and enthalpy. The progress variable allows to efficiently include a detailed reaction scheme, while the dependence on enthalpy describes the effect of heat losses on the reaction rate. The turbulence-chemistry interaction is modeled by eight Eulerian stochastic fields. An LES of a premixed swirl burner with a broadband velocity excitation is performed to investigate the flame dynamics, i.e., the response of heat release rate to upstream velocity perturbations. In particular, the flame impulse response and the flame transfer function (FTF) are identified from LES time series data. Simulation results for a range of power ratings are in good agreement with the experimental data.


2012 ◽  
Vol 2012 ◽  
pp. 1-11 ◽  
Author(s):  
Rohit Kulkarni ◽  
Wolfgang Polifke

The potential of a progress variable formulation for predicting autoignition and subsequent kernel development in a nonpremixed jet flame is explored in the LES (Large Eddy Simulation) context. The chemistry is tabulated as a function of mixture fraction and a composite progress variable, which is defined as a combination of an intermediate and a product species. Transport equations are solved for mixture fraction and progress variable. The filtered mean source term for the progress variable is closed using a probability density function of presumed shape for the mixture fraction. Subgrid fluctuations of the progress variable conditioned on the mixture fraction are neglected. A diluted hydrogen jet issuing into a turbulent coflow of preheated air is chosen as a test case. The model predicts ignition lengths and subsequent kernel growth in good agreement with experiment without any adjustment of model parameters. The autoignition length predicted by the model depends noticeably on the chemical mechanism which the tabulated chemistry is based on. Compared to models using detailed chemistry, significant reduction in computational costs can be realized with the progress variable formulation.


Fluids ◽  
2020 ◽  
Vol 5 (3) ◽  
pp. 126 ◽  
Author(s):  
Kai Zhang ◽  
Ali Ghobadian ◽  
Jamshid M. Nouri

The scale-resolving simulation of a practical gas turbine combustor is performed using a partially premixed finite-rate chemistry combustion model. The combustion model assumes finite-rate chemistry by limiting the chemical reaction rate with flame speed. A comparison of the numerical results with the experimental temperature and species mole fraction clearly showed the superiority of the shear stress transport, K-omega, scale adaptive turbulence model (SSTKWSAS). The model outperforms large eddy simulation (LES) in the primary region of the combustor, probably for two reasons. First, the lower amount of mesh employed in the simulation for the industrial-size combustor does not fit the LES’s explicit mesh size dependency requirement, while it is sufficient for the SSTKWSAS simulation. Second, coupling the finite-rate chemistry method with the SSTKWSAS model provides a more reasonable rate of chemical reaction than that predicted by the fast chemistry method used in LES simulation. Other than comparing with the LES data available in the literature, the SSTKWSAS-predicted result is also compared comprehensively with that obtained from the model based on the unsteady Reynolds-averaged Navier–Stokes (URANS) simulation approach. The superiority of the SSTKWSAS model in resolving large eddies is highlighted. Overall, the present study emphasizes the effectiveness and efficiency of coupling a partially premixed combustion model with a scale-resolving simulation method in predicting a swirl-stabilized, multi-jets turbulent flame in a practical, complex gas turbine combustor configuration.


Author(s):  
Murase Kagenobu ◽  
Oshima Nobuyuki ◽  
Takahashi Yusuke

This paper focuses on the numerical simulation of Sandia National Laboratories “the piloted methane/air burner flame D.” Large Eddy Simulation and 2-scalar flamelet approach are applied for the turbulent and partially premixed combustion field, which is expressed by the LES filtered equations of scalar G for tracking the flame surfaces and mixture fraction of a fuel and an oxidizer. The flamelet data consists of temperature, specific volume and laminar flame speed are calculated by the detail chemical reaction with GRI-Mech 3.0. Two kinds of flamelet data are validated; one is “equilibrium flamelet data” calculated by 0-dimensional equilibrium solution based on equilibrium model; the other is “diffusion flamelet data” calculated by 1-dimensional counter flow solution based on laminar flamelet model. Consequently, the “diffusion flamelet data” gives better result in this type of combustion field.


2000 ◽  
Vol 1 (2) ◽  
pp. 209-227 ◽  
Author(s):  
S Menon

Next-generation gas turbine and internal combustion engines are required to reduce pollutant emissions significantly and also to be fuel efficient. Accurate prediction of pollutant formation requires proper resolution of the spatio-temporal evolution of the unsteady mixing and combustion processes. Since conventional steady state methods are not able to deal with these features, methodology based on large-eddy simulations (LESs) is becoming a viable choice to study unsteady reacting flows. This paper describes a new LES methodology developed recently that has demonstrated a capability to simulate reacting turbulent flows accurately. A key feature of this new approach is the manner in which small-scale turbulent mixing and combustion processes are simulated. This feature allows proper characterization of the effects of both large-scale convection and small-scale mixing on the scalar processes, thereby providing a more accurate prediction of chemical reaction effects. LESs of high Reynolds number premixed flames in the flamelet regime and in the distributed reaction regime are used to describe the ability of the new subgrid combustion model.


Author(s):  
Savvas Gkantonas ◽  
Jenna M. Foale ◽  
Andrea Giusti ◽  
Epaminondas Mastorakos

Abstract The simulation of soot evolution is a problem of relevance for the development of low-emission aero-engine combustors. Apart from detailed CFD approaches, it is important to also develop models with modest computational cost so a large number of geometries can be explored, especially in view of the need to predict engine-out soot particle size distributions (PSDs) to meet future regulations. This paper presents an approach based on Incompletely Stirred Reactor Network (ISRN) modeling that simplifies calculations, allowing for the use of very complex chemistry and soot models. The method relies on a network of Incompletely Stirred Reactors (ISRs), which are inhomogeneous in terms of mixture fraction but characterized by homogeneous conditional averages, with the conditioning performed on the mixture fraction. The ISRN approach is demonstrated here for a single sector lean-burn model combustor operating on Jet-A1 fuel in pilot-only mode, for which detailed CFD and experimental data are available. Results show that reasonable accuracy is obtained at a significantly reduced computational cost. Real fuel chemistry and a detailed physicochemical sectional soot model are consequently employed to investigate the sensitivity of ISRN predictions to the chemical mechanism chosen and to provide an estimate of the soot particle size distribution at the combustor exit.


2020 ◽  
Author(s):  
Rui Zhang ◽  
Chengwen Zhong ◽  
Sha Liu ◽  
Congshan Zhuo

Abstract In this paper, we introduce the incompressible discrete Maxwellian equilibrium distribution function and external forces into the two-stage third-order Discrete Unified Gas-Kinetic Scheme (DUGKS) for simulating low-speed incompressible turbulent flows with forcing term. The Wall-Adapting Local Eddy-viscosity (WALE) and Vreman sub-grid models for Large-Eddy Simulations (LES) of wall-bounded turbulent flows are coupled within the present framework. In order to simulate the three-dimensional turbulent flows associated with great computational cost, a parallel implementation strategy for the present framework is developed, and is validated by three canonical wall-bounded turbulent flows, viz., the fully developed turbulent channel flow at a friction Reynolds number (Re) about 180, the turbulent plane Couette flow at a friction Re number about 93 and three-dimensional lid-driven cubical cavity flow at a Re number of 12000. The turbulence statistics are computed by the present approach with both WALE and Vreman models, and their predictions match precisely with each other. Especially, the predicted flow physics of three-dimensional lid-driven cavity are consistent with the description from abundant literatures. While, they have small discrepancies in comparison to the Direct Numerical Simulation (DNS) due to the relatively low grid resolution. The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.


2020 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Alain Fossi ◽  
Alain DeChamplain ◽  
Benjamin Akih-Kumgeh ◽  
Jeffrey Bergthorson

Purpose This study aims to deal with the large eddy simulation (LES) of an ignition sequence and the resulting steady combustion in a swirl-stabilized liquid-fueled combustor. Particular attention is paid to the ease of handling the numerical tool, the accuracy of the results and the reasonable computational cost involved. The primary aim of the study is to appraise the ability of the newly developed computational fluid dynamics (CFD) methodology to retrieve the spark-based flame kernel initiation, its propagation until the full ignition of the combustion chamber, the flame stabilization and the combustion processes governing the steady combustion regime. Design/methodology/approach The CFD model consists of an LES-based spray module coupled to a subgrid-scale ignition model to capture the flame kernel initiation and the early stage of the flame kernel growth, and a combustion model based on the mixture fraction-progress variable formulation in the line of the flamelet generated manifold (FGM) method to retrieve the subsequent flame propagation and combustion properties. The LES-spray module is based on an Eulerian-Lagrangian approach and includes a fully two-way coupling at each time step to account for the interactions between the liquid and the gaseous phases. The Wall-Adapting Local Eddy-viscosity (WALE) model is used for the flow field while the eddy diffusivity model is used for the scalar fluxes. The fuel is liquid kerosene, injected in the form of a polydisperse spray of droplets. The spray dynamics are tracked using the Lagrangian procedure, and the phase transition of droplets is calculated using a non-equilibrium evaporation model. The oxidation mechanism of the Jet A-1 surrogate is described through a reduced reaction mechanism derived from a detailed mechanism using a species sensitivity method. Findings By comparing the numerical results with a set of published data for a swirl-stabilized spray flame, the proposed CFD methodology is found capable of capturing the whole spark-based ignition sequence in a liquid-fueled combustion chamber and the main flame characteristics in the steady combustion regime with reasonable computing costs. Research limitations/implications The proposed CFD methodology simulates the whole ignition sequence, namely, the flame kernel initiation, its propagation to fully ignite the combustion chamber, and the global flame stabilization. Due to the lack of experimental ignition data on this liquid-fueled configuration, the ability of the proposed CFD methodology to accurately predict ignition timing was not quantitatively assessed. It would, therefore, be interesting to apply this CFD methodology to other configurations that have experimental ignition data, to quantitatively assess its ability to predict the ignition timing and the flame characteristics during the ignition sequence. Such further investigations will not only provide further validation of the proposed methodology but also will potentially identify its shortfalls for better improvement. Practical implications This CFD methodology is developed by customizing a commercial CFD code widely used in the industry. It is, therefore, directly applicable to practical configurations, and provides not only a relatively straightforward approach to predict an ignition sequence in liquid-fueled combustion chambers but also a robust way to predict the flame characteristics in the steady combustion regime as significant improvements are noticed on the prediction of slow species. Originality/value The incorporation of the subgrid ignition model paired with a combustion model based on tabulated chemistry allows reducing computational costs involved in the simulation of the ignition phase. The incorporation of the FGM-based tabulated chemistry provides a drastic reduction of computing resources with reasonable accuracy. The CFD methodology is developed using the platform of a commercial CFD code widely used in the industry for relatively straightforward applicability.


Author(s):  
D. Mira ◽  
M. Vázquez ◽  
G. Houzeaux ◽  
S. Gövert ◽  
J. W. B. Kok ◽  
...  

The primary purpose of this study is to evaluate the ability of LES, with a turbulent combustion model based on steady flamelets, to predict the flame stabilization mechanisms in an industrial can combustor at full load conditions. The test case corresponds to the downscaled Siemens can combustor tested in the high pressure rig at the DLR. The effects of the wall temperature on the prediction capabilities of the codes is investigated by imposing several heat transfer conditions at the pilot and chamber walls. The codes used for this work are Alya and OpenFOAM, which are well established CFD codes in the fluid mechanics community. Prior to the simulation, results for 1-D laminar flames at the operating conditions of the combustor are compared with the detailed solutions. Subsequently, results from both codes at the mid-plane are compared against the experimental data available. Acceptable results are obtained for the axial velocity, while discrepancies are more evident for the mixture fraction and the temperature, particularly with Alya. However, both codes showed that the heat losses influence the size and length of the pilot and main flame.


Sign in / Sign up

Export Citation Format

Share Document