CFD Simulations AC2-09
SANDIA Flame D
Application Challenge AC2-09 © copyright ERCOFTAC 2024
Overview of CFD Simulations
All the calculations presented below were obtained within the MOLECULES FP5 project Contract N° G4RD-CT-2000-00402 by the team of the Institute of Thermal Machinery, Częstochowa University of Technology. The computations were performed with BOFFIN-LES code developed at Imperial College by the group of Prof. W.P. Jones. BOFFIN-LES computer code utilizes a boundary conforming general curvilinear coordinate system with a co-located storage arrangement. It incorporates a fully implicit formulation and is second order accurate in space and time. For the convection terms an energy conserving discretization scheme is used and matrix preconditioned conjugate gradient methods are used to solve the equations for pressure and velocity etc. The CFD simulations are all LES predictions with various subgrid scale models and turbulence/combustion interaction approaches and neither RANS nor URANS methods are studied in this document.
In the LES calculations two models of turbulence/combustion interaction were applied: steady flamelet model and simplified Conditional Moment Closure (CMC) neglecting the convection term in physical space (The CMC module was developed by Prof. E. Mastorakos from Cambridge University). In both cases the standard subgrid-scale (SGS) Smagorinsky model was used. Then in order to evaluate the importance of the subgrid-scale models the LES calculations were also performed using steady flamelet and dynamic (Germano) SGS model.
SIMULATION CASE CFD1
Solution Strategy
In the CFD1 the steady flamelet concept was applied with the standard Smagorinsky SGS model for turbulence.
Computational Domain
The CFD1 results were obtained with computational meshes 80×80×160 nodes. The computational domain at the inlet and outlet plane extended to 5.5D and 18.3D respectively in both horizontal directions. The length of the domain was equal to 50D. The mesh was stretched in axial direction by exponential function and in radial directions by hyperbolic tangent function. The grid refinement studies for the LES calculations showed that the grid resolution with 80×80×160 nodes in the proposed computational domain is sufficient and further grid refinement leads to minor changes of the statistically converged parameters. The computational domain is shown in Fig.4.
![]() |
![]() |
![]() | |
Fig. 4. Computational domain for Sandia flame D (left); mesh resolution in the inlet plane (right) |
Boundary Conditions
The boundary conditions in the inlet plane were assumed to be as follows:
- the mean and RMS profiles of the axial velocity component were interpolated from experimental data - Fig.3 presents comparison of the experimental data with boundary profiles applied in computations; the random disturbances, introduced as a white noise, were scaled by RMS profile and next they were superimposed on the mean profile;
- the mean and RMS values of the radial velocity components were assumed equal to zero;
- the mixture fraction was assumed equal to 1.0 in the main jet; 0.27 in the pilot jet and zero in the co-flowing air.
At the lateral boundaries the axial velocity was assumed equal to the velocity of coflowing air (0.9 m/s) while the remaining components were equal to zero. At the outlet the convective type boundary conditions were assumed which do not require specification of any variables.
Application of Physical Models
In the most general case modeling of the combustion processes is very expensive computationally since together with the solution of the flow field it requires solution of additional transport equations for particular N species (e.g. CO, CO2, H2O, H2, etc.) produced in chemical reactions. The transport equations for species have the following form:
where:
is the density | |
is the velocity component | |
is the mass fraction of species | |
is the reaction rate (speed of creation/destruction of a given species) | |
is the diffusion coefficient usually taken the same (denoted by ) for each species and defined as , where is the molecular viscosity and is the Prandtl number. |
Reaction rate of a given species k is a sum of the reaction rates in all M reactions in which species k occurs. It is defined as:
where:
is the reaction rate of species in reaction | |
are the molar stochiometric coefficients after and before reaction respectively | |
is the atomic weight of species | |
is the rate of progress of reaction ; it is a function of temperature, density and species mass fraction and may be obtained from chemical kinetics or experiment; methods of determination of the rate of progress of reaction are beyond the scope of this report. |
With the assumption of low Mach number flow the equation of energy may be expressed as the transport equation for the temperature which is given as [1]:
where:
is the temperature | |
is the heat diffusion coefficient | |
is the specific heat | |
is the heat release defined as a sum of product of reaction rates and formation enthalpies |
The equations given above together with the Navier-Stokes equations, the equation of state and the continuity equation form a closed system which allows computing the flow field together with combustion process. However, their direct implementation in a computer code with regard for tens of species and tens (or even several hundred) of chemical reactions is still impossible from the point of view of capability of available computers, and for this reason significant simplifications have to be made.
The turbulence/combustion interaction steady flamelet concept implemented in the BOFFIN code was introduced by Peters [2], which stated that the flame can be seen as an ensemble of laminar flamelets. It used the equation for the conserved scalar referred to as the mixture fraction. The equation for the mixture fraction (denoted as ) has the simple form of convection-diffusion equation and is given as:
The mixture fraction is a normalized quantity
() and represents a
local fuel to oxidizer ratio ([pic]- pure oxidizer; [pic]- pure fuel).
The assumption that one conserved scalar is sufficient to describe
thermochemical state of the flow decouples the modelling of reactive
phenomena from that of flow modelling. Assuming that particular species
and temperature are functions of the mixture fraction ([pic], [pic]) the
equations (1) and (3) may be transformed [1] into the mixture fraction
space resulting in:
Numerical Accuracy
CFD Results
References
3. Branley N., Jones W.P.: "Large eddy simulations of a turbulent non-premixed flame", Combustion and Flame, vol. 127, 2001
4. Perry R.H. and Green D.W.: Chemical Engineer's Handbook, McGraw-Hill, 1993
SIMULATION CASE CFD2
(as per CFD 1)
Contributed by: Andrzej Boguslawski — Technical University of Częstochowa
© copyright ERCOFTAC 2024