Elsevier

Combustion and Flame

Volume 172, October 2016, Pages 153-161
Combustion and Flame

Joint experimental and numerical study of the influence of flame holder temperature on the stabilization of a laminar methane flame on a cylinder

https://doi.org/10.1016/j.combustflame.2016.06.025Get rights and content

Abstract

The mechanisms controlling laminar flame anchoring on a cylindrical bluff-body are investigated using DNS and experiments. Two configurations are examined: water-cooled and uncooled steel cylinders. Comparisons between experimental measurements and DNS show good agreement for the flame root locations in the two configurations. In the cooled case, the flame holder is maintained at about 300 K and the flame is stabilized in the wake of the cylinder, in the recirculation zone formed by the products of combustion. In the uncooled case, the bluff-body reaches a steady temperature of about 700 K in both experiment and DNS and the flame is stabilized closer to it. The fully coupled DNS of the flame and the temperature field in the bluff-body also shows that capturing the correct radiative heat transfer from the bluff-body is a key ingredient to reproduce experimental results.

Introduction

The burnt gas temperatures reached in combustion chambers usually exceeds the maximum temperatures which can be sustained by most materials, especially metals used in engines. Therefore, cooling these walls as well as all chamber elements in contact with the flame is mandatory for combustion chamber designers. While cooling is obviously needed to preserve walls, its effects on the flames themselves has received less attention and is usually neglected in many CFD approaches. Flame/wall interaction, for example, is a field of combustion which has not been investigated yet with sufficient care [1], [2], [3], [4], [5], [6]. In most cases, authors measure or compute the maximum wall heat fluxes induced by the flame but do not investigate the effects of the wall on the flame itself.

In the field of simulation, most models [7], [8], [9], [10], [11] assume adiabatic flows. For premixed flames, the famous BML (Bray Moss Libby) approach, for example, which is the workhorse of many theories for turbulent premixed flames [12], [13] assumes that a single variable (the progress variable c) is sufficient to describe the flow: this is true only when the flow is adiabatic. In the same way, many usual methods for chemistry tabulation such as FPV [14], FPI [15] or FGM [16] assume that chemistry can be described using only two variables, the mixture fraction z and the progress variable c, which implies that the flames must be adiabatic.1 Considering that wall heat fluxes in most chambers correspond to approximately 5–40% of the chamber total power, assuming adiabaticity is clearly not compatible with the high-precision methods which are sought today. Note that computing the interaction between the flame and the wall requires to compute both the flow and the temperature within the walls simultaneously: the LES code must be coupled with a heat transfer code within the combustor walls. This task is not simple [19], [20] because time scales are usually very different (a few milliseconds in the flow and a few minutes in the walls).

Among all walls present in a chamber, flame holders play a special role because they control the most sensitive zone of the chamber: the place where the flames are anchored. Any temperature change of the flame holder will induce a change of position for the flame roots and therefore a change in stability and efficiency. The coupling mechanisms between heat transfer within flameholder and flame stabilization have not been analyzed in detail yet. In a series of recent papers [21], [22], [23], the MIT group has numerically studied the stabilization of premixed flames on square flame holders and shown that the location of the flame roots but also the blow-off limits were strongly affected by the temperature of the flame holder.

The present study focuses on a similar question: which differences in flame anchoring are observed when the temperature of the flame holder varies from a low (typically 300 K) to a high value (700 K). To obtain such a large variation in temperature, a premixed laminar methane/air flame is stabilized on a cylindrical flame holder. Two flame holders are used, with exactly the same external shape. The first one has an internal water cooling system, leading to a surface temperature close to 300 K. The second one is a full, solid cylinder which is uncooled, leading to a surface temperature close to 700 K.

Both experiments and DNS are used to analyze the differences in flame structure near the flame holder. Simulations are performed in dual mode: the flow is computed with DNS using a 13 species kinetic scheme for CH4/air flames [24] while the temperature in the solid is computed with a heat transfer solver, coupled to the flow solver. The simulations, performed for cooled and uncooled flame holders, reveal drastic differences in flame root location and flow topologies. They also show that radiative heat transfer must be taken into account to predict the flame topology for the uncooled case.

Section 2 presents the experimental setup. The tools used for the coupled flow/solid simulation are described in Section 3. Results for the cooled flame holder are discussed in Section 4 before presenting the uncooled case in Section 5. Finally Section 6 discusses the influence of radiative heat fluxes on the flame stabilization when the flame holder is uncooled.

Section snippets

Experimental configuration

The experimental rig is shown in Fig. 1: a lean premixed methane-air V-flame is stabilized over steel cylindrical bluff body (radius of r=4 mm). The burner has a constant cross section of h=34 by l=94 mm so that the flame remains two-dimensional. Individual mass flow meters are used to control air and methane flow rates. Fuel and oxidizer are premixed before entering the injection chamber though six holes. Glass wool, small glass balls and two honeycombs panels are used to laminarize the flow.

Numerical strategy

To capture the effects of flame holder cooling on the flame, a coupled DNS of the flow and of the temperature field within the flame holder is performed.

CBB configuration

Since the bluff-body temperature is controlled by a water flow which is not computed, a boundary condition at the inner diameter of the flame holder is required. The convective inner flux ϕs → w (Fig. 2) is modeled through a Newton law at the solid/water boundary: ϕsw=hturb(TsTw)where Ts is the local inside skin temperature of the cylinder, Tw is the mean temperature of the cooling water in the outer passage and hturb is the heat transfer coefficient. The water flow is turbulent with a

UBB configuration

For the uncooled flame holder, a steady symmetrical flame is also observed. The comparison between DNS and experiments is very good (Fig. 9, left). Compared to the CBB case (Fig. 5, left), the flame is much closer to the flame holder. The radial heat flux in the fluid region can be used to determine its angle. The angle θUBB corresponds to the azimuthal point where the heqt flux changes sign: θUBB=θ/ϕfluid.n=0 where n is the normal unit vector pointing inside of the cylinder. As shown in

Influence of the cylinder emissivity

Section 5 showed that radiative heat transfer represents 45% of the bluff-body heat losses for the UBB case. This suggests that changes in stabilization mechanisms may be induced by changing the flame holder emissivity and its temperature. Different computations have been carried out with emissivity ranging between 0.02 and 1.

DNS results show that the flame root position (Fig. 14) is roughly independent of the emissivity between ϵ=0.8 (weakly oxidized bluff-body) and ϵ=1.0 (perfect black body).

Conclusions

A comparison between DNS and experiments has been carried out to study the anchoring mechanism of a flame attached on a cylindrical bluff-body. Two distinct configurations have been scrutinized. In the cooled bluff-body case (CBB) the flame holder temperature is about 300 K and the flame is stabilized approximately one radius downstream of the cylinder at an angle θCBB ≈ 0.15π and attached at the center of the recirculation zone in which products and reactants are mixed. In the uncooled

Aknowledgments

This work was granted access to the high-performance computing resources of CINES under the allocation x20152b7036 made by Grand Equipement National de Calcul Intensif. The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) / ERC Grant agreement ERC-AdG 319067-INTECOCIS. The authors also thank A. Felden, who implemented the LU13 scheme in the flow solver and Dr. A. Ghani for his helpful

References (47)

  • R.B. Price et al.

    Optical studies of the generation of noise in turbulent flames

    Symp. (Int.) Combust.

    (1969)
  • O. Colin et al.

    Development of high-order Taylor Galerkin schemes for LES

    J. Comput. Phys.

    (2000)
  • V. Moureau et al.

    Numerical methods for unsteady compressible multi-component reacting flows on fixed and moving grids

    J. Comput. Phys.

    (2005)
  • T. Poinsot et al.

    Boundary conditions for direct simulations of compressible viscous flows

    J. Comput. Phys.

    (1992)
  • J.O. Hirschfelder et al.

    Molecular theory of gases and liquids

    (1954)
  • I. Dyakov et al.

    Measurement of adiabatic burning velocity in methane–oxygen–nitrogen mixtures

    Combust. Sci. Technol.

    (2001)
  • V.A.L. Erie et al.

    A set of GMRES routines for real and complex arithmetics on high performance computers

    ACM Trans.

    (2005)
  • DouH.S. et al.

    Critical condition for flow transition in a full- developed annulus flow

    J. Pet. Sci. Eng.

    (2000)
  • G.M. Makhviladze et al.

    Flame propagation in a closed channel with cold side walls

    Chem. Eng. Sci.

    (1991)
  • A. Gruber et al.

    Direct numerical simulation of premixed flame boundary layer flashback in turbulent channel flow

    J. Fluid Mech.

    (2012)
  • LuJ.H. et al.

    Unsteady heat transfer during side wall quenching of a laminar flame

    Symp. (Int.) Combust.

    (1990)
  • G. Bruneaux et al.

    Premixed flame–wall interaction in a turbulent channel flow: budget for the flame surface density evolution equation and modelling

    J. Fluid Mech.

    (1997)
  • A. Gruber et al.

    Turbulent flame wall interaction: a direct numerical simulation study

    J. Fluid Mech.

    (2010)
  • Cited by (0)

    View full text