Propane/air diffusion flame

Authors: Steward F.R., Tennankore K.N.

Type: Experimental

Status: tip



Flow Parameters

Reference Publications



The main objective of the present research is to study a turbulent burning propane/air diffusion flame in a cylindrical combustion chamber. The results obtained are compared with experimental data [Stew79]. The flame under investigation is a turbulent jet diffusion flame showing an axisymmetric behaviour. The fuel enters the computational domain (see figure 1) via a centred inlet (tube), while the combustion air is injected as an axial co-flow. Due to the axisymmetric character of the flame, the mesh size can be significantly reduced rotating the 2-D plane of figure 1 by only some degrees. The lateral faces are considered to be rotational boundary conditions.

Figure 1: Schematic of the computational domain.

Flow Parameters

In this study, the simulations performed are based on either the 2-equation k-ε model or the two-equation k-ω model.
The radiative heat transfer is calculated by means of the discrete ordinate method. Nowadays, this method represents one of the most used in CFD simulations. In order to optimize the correlation between the method’s correctness and the numerical effort, the number of considered propagation directions has been limited. In the present sutdy, 32 directions are considered to describe the spatial propagation of radiation. The absorption efficiency of the gas phase is described by Modak’s approximation [Mod79]. Conductive heat transfer is not considered.
Physical parameters:
  • dynamic viscosity : %$\mu = 1.6 \times10^{-5} kg.m^{-1}.s^{-1}$%
  • density : %$\rho = 1.17861 kg.m^{-3}$%
  • specific heat capacity : %$ c_{p} = 1017.24 J.K^{-1}.kg^{-1}$%
  • reference pressure : %$ p = 101300Pa$%
  • reference temperature : %$ T = 298.15K$%

Boundary conditions

Fuel inlet:
* Total flowrate: %$2.62609\times10^{-4}kg.s^{-1}$%
* Temperature : %$436 K$%
* Mixing Rate : %$1$%

Air inlet:
* Total flowrate: %$4.282\times10^{-3}kg.s^{-1}$%
* Temperature : %$353 K$%
* Mixing Rate : %$0$%

The temperature profile as well as the wall specific radiation emissivities are given as function of the y-coordinate in the following table:

Distance y (m) Temperature (K) Emissivity
from 0 to 0.09125 415 0.550
from 0.09125 to 0.18250 415 0.535
from 0.18250 to 0.27375 426 0.560
from 0.27375 to 0.36500 460 0.565
from 0.36500 to 0.45625 497 0.570
from 0.45625 to 0.54750 541 0.575
from 0.54750 to 0.63875 586 0.580
from 0.63875 to 2.50000 611 0.585

Lateral faces:
The lateral faces are defined as symmetric boundary conditions.

Reference Publications

[Stew79] Steward F.R., Tennankore K.N. : "Towards a finite difference solution coupled with the zone method for radiative transfer for a cylindrical combustion chamber." Journal of the Institute of Energey, pp.107-114, september 1979- [Mod79] Modak, A. T. : “Exponantial wide band parameters for the pure rotational band of water vapor“ Journal of quantitative Spectroscopy and radiative transfer, vol 21, pp.131 - 142,1979.


Simulation results available for this case:
Code Version Author Restrictions
Code_Saturne 2.0-beta2 M. Hassanaly AccessEDFGroup
Number of topics: 1


Add New Results Here.

To add results please fill in the form by clicking on the CREATE button below. You will need to provide all the information on the form.

Please select which group or groups can access the information from the list. If you don't want to restrict the access, please select None.

  • New Test Case Result topic:

Current Tags:
Saturne_Validation_2.01Add my vote for this tag create new tag
, view all tags
Topic attachments
I Attachment Action Size Date Who Comment
pngpng geom.png manage 4.0 K 2010-09-10 - 13:10 MarcusCharwath  
Topic revision: r60 - 2017-05-19 - 03:29:30 - AllenZhang

Computational Fluid Dynamics and Turbulence Mechanics
@ the University of Manchester
Copyright &© by the contributing authors. Unless noted otherwise, all material on this web site is the property of the contributing authors.