CFD Simulations AC2-09: Difference between revisions

From KBwiki
Jump to navigation Jump to search
Line 37: Line 37:
The CFD1 results were  obtained  with  computational  meshes  80×80×160
The CFD1 results were  obtained  with  computational  meshes  80×80×160
nodes. The computational domain at the inlet and outlet plane  extended
nodes. The computational domain at the inlet and outlet plane  extended
to 5.5D and 18.3D respectively  in  both  horizontal  directions.  The
to 5.5''D'' and 18.3''D'' respectively  in  both  horizontal  directions.  The
length of the domain was equal to 50D. The mesh was stretched in  axial
length of the domain was equal to 50''D''. The mesh was stretched in  axial
direction  by  exponential  function  and  in  radial  directions  by
direction  by  exponential  function  and  in  radial  directions  by
hyperbolic tangent function. The grid refinement studies  for  the  LES
hyperbolic tangent function. The grid refinement studies  for  the  LES
Line 45: Line 45:
refinement leads  to  minor  changes  of  the  statistically  converged
refinement leads  to  minor  changes  of  the  statistically  converged
parameters. The computational domain is shown in Fig.4.
parameters. The computational domain is shown in Fig.4.
===Boundary Conditions===
===Boundary Conditions===
<!--{{Demo_AC_CFD_BC}}-->
<!--{{Demo_AC_CFD_BC}}-->

Revision as of 10:41, 29 April 2011

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

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.

Boundary Conditions

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.

Numerical Accuracy

CFD Results

References

SIMULATION CASE CFD2

(as per CFD 1)



Contributed by: Andrzej Boguslawski — Technical University of Częstochowa

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice


© copyright ERCOFTAC 2024