Physics:Premixed flame

From HandWiki
Different flame types of a Bunsen burner depend on oxygen supply. On the left a rich fuel mixture with no premixed oxygen produces a yellow sooty diffusion flame, and on the right a lean fully oxygen premixed flame produces no soot and the flame color is produced by molecular radical band emission.

A premixed flame is a flame formed under certain conditions during the combustion of a premixed charge (also called pre-mixture) of fuel and oxidiser. Since the fuel and oxidiser—the key chemical reactants of combustion—are available throughout a homogeneous stoichiometric premixed charge, the combustion process once initiated sustains itself by way of its own heat release. The majority of the chemical transformation in such a combustion process occurs primarily in a thin interfacial region which separates the unburned and the burned gases. The premixed flame interface propagates through the mixture until the entire charge is depleted.[1] The propagation speed of a premixed flame is known as the flame speed (or burning velocity) which depends on the convection-diffusion-reaction balance within the flame, i.e. on its inner chemical structure. The premixed flame is characterised as laminar or turbulent depending on the velocity distribution in the unburned pre-mixture (which provides the medium of propagation for the flame).

Premixed flame propagation

Laminar

Under controlled conditions (typically in a laboratory) a laminar flame may be formed in one of several possible flame configurations. The inner structure of a laminar premixed flame is composed of layers over which the decomposition, reaction and complete oxidation of fuel occurs. These chemical processes are much faster than the physical processes such as vortex motion in the flow and, hence, the inner structure of a laminar flame remains intact in most circumstances. The constitutive layers of the inner structure correspond to specified intervals over which the temperature increases from the specified unburned mixture up to as high as the adiabatic flame temperature (AFT). In the presence of volumetric heat transfer and/or aerodynamic stretch, or under the development intrinsic flame instabilities, the extent of reaction and, hence, the temperature attained across the flame may be different from the AFT.

Laminar burning velocity

For a one-step irreversible chemistry, i.e., [math]\displaystyle{ \nu_F \rm{F} + \nu_O \rm{O}_2 \rightarrow \rm{Products} }[/math], the planar, adiabatic flame has explicit expression for the burning velocity derived from activation energy asymptotics when the Zel'dovich number [math]\displaystyle{ \beta\gg 1. }[/math] The reaction rate [math]\displaystyle{ \omega }[/math] (number of moles of fuel consumed per unit volume per unit time) is taken to be Arrhenius form,

[math]\displaystyle{ \omega = B \left(\frac{\rho Y_F}{W_{F}}\right)^m \left(\frac{\rho Y_{O_2}}{W_{O_2}}\right)^n e^{-E_a/RT}, }[/math]

where [math]\displaystyle{ B }[/math] is the pre-exponential factor, [math]\displaystyle{ \rho }[/math] is the density, [math]\displaystyle{ Y_F }[/math] is the fuel mass fraction, [math]\displaystyle{ Y_{O_2} }[/math] is the oxidizer mass fraction, [math]\displaystyle{ E_a }[/math] is the activation energy, [math]\displaystyle{ R }[/math] is the universal gas constant, [math]\displaystyle{ T }[/math] is the temperature, [math]\displaystyle{ W_F\ \& \ W_{O_2} }[/math] are the molecular weights of fuel and oxidizer, respectively and [math]\displaystyle{ m\ \& \ n }[/math] are the reaction orders. Let the unburnt conditions far ahead of the flame be denoted with subscript [math]\displaystyle{ u }[/math] and similarly, the burnt gas conditions by [math]\displaystyle{ b }[/math], then we can define an equivalence ratio [math]\displaystyle{ \phi }[/math] for the unburnt mixture as

[math]\displaystyle{ \phi = \frac{\nu_{O_2} W_{O_2}}{\nu_F W_F}\frac{Y_{F,u}}{Y_{O_2,u}} }[/math].

Then the planar laminar burning velocity for fuel-rich mixture ([math]\displaystyle{ \phi\gt 1 }[/math]) is given by[2][3]

[math]\displaystyle{ S_L=\left\{\frac{2 B\lambda_b \rho_b^{m+n} \nu_F^m Y_{O_2,u}^{m+n-1} G(n,m,a)}{ c_{p,b} \rho_u^2 \nu_{O_2} W_{O_2}^{m+n-1} \beta^{m+n+1} \mathrm{Le}_{O_2}^{-n} \mathrm{Le}_F^{-m}}\right\}^{1/2} e^{-E_a/2RT_b} + O(\beta^{-1}), }[/math]

where

[math]\displaystyle{ G(n,m,a) = \int_0^\infty y^n (y+a)^m\ dy }[/math]

and [math]\displaystyle{ a=\beta(\phi-1)/\mathrm{Le}_F }[/math]. Here [math]\displaystyle{ \lambda }[/math] is the thermal conductivity, [math]\displaystyle{ c_p }[/math] is the specific heat at constant pressure and [math]\displaystyle{ \mathrm{Le} }[/math] is the Lewis number. Similarly one can write the formula for lean [math]\displaystyle{ \phi\lt 1 }[/math] mixtures. This result is first obtained by T. Mitani in 1980.[4] Second order correction to this formula with more complicated transport properties were derived by Forman A. Williams and co-workers in the 80s.[5][6][7]

Variations in local propagation speed of a laminar flame arise due to what is called flame stretch. Flame stretch can happen due to the straining by outer flow velocity field or the curvature of flame; the difference in the propagation speed from the corresponding laminar speed is a function of these effects and may be written as: [8][9]

[math]\displaystyle{ U_L = S_L - S_L \mathcal{M}_c\delta_L \kappa + \mathcal{M}_a \delta_L\mathbf{n}\cdot\nabla\mathbf{n}\cdot\mathbf{n} }[/math]

where [math]\displaystyle{ \delta_L }[/math] is the laminar flame thickness, [math]\displaystyle{ \kappa }[/math] is the flame curvature, [math]\displaystyle{ \mathbf{n} }[/math] is the unit normal on the flame surface pointing towards the unburnt gas side, [math]\displaystyle{ \mathbf{v} }[/math] is the flow velocity and [math]\displaystyle{ \mathcal{M}_c\ \& \ \mathcal{M}_a }[/math] are the respective Markstein numbers of curvature and strain.

Turbulent

In practical scenarios, turbulence is inevitable and, under moderate conditions, turbulence aids the premixed burning process as it enhances the mixing process of fuel and oxidiser. If the premixed charge of gases is not homogeneously mixed, the variations on equivalence ratio may affect the propagation speed of the flame. In some cases, this is desirable as in stratified combustion of blended fuels.

A turbulent premixed flame can be assumed to propagate as a surface composed of an ensemble of laminar flames so long as the processes that determine the inner structure of the flame are not affected.[10] Under such conditions, the flame surface is wrinkled by virtue of turbulent motion in the premixed gases increasing the surface area of the flame. The wrinkling process increases the burning velocity of the turbulent premixed flame in comparison to its laminar counterpart.

The propagation of such a premixed flame may be analysed using the field equation called as G equation[11][12] for a scalar [math]\displaystyle{ G }[/math] as:

[math]\displaystyle{ \frac{\partial G}{\partial t} + \mathbf{v} \cdot \nabla G = U_L |\nabla G| }[/math],

which is defined such that the level-sets of G represent the various interfaces within the premixed flame propagating with a local velocity [math]\displaystyle{ U_L }[/math]. This, however, is typically not the case as the propagation speed of the interface (with resect to unburned mixture) varies from point to point due to the aerodynamic stretch induced due to gradients in the velocity field.

Under contrasting conditions, however, the inner structure of the premixed flame may be entirely disrupted causing the flame to extinguish either locally (known as local extinction) or globally (known as global extinction or blow-off). Such opposing cases govern the operation of practical combustion devices such as SI engines as well as aero-engine afterburners. The prediction of the extent to which the inner structure of flame is affected in turbulent flow is a topic of extensive research.

Premixed flame configuration

The flow configuration of premixed gases affects the stabilization and burning characteristics of the

Bunsen flame

In a Bunsen flame, a steady flow rate is provided which matches the flame speed so as to stabilize the flame. If the flow rate is below the flame speed, the flame will move upstream until the fuel is consumed or until it encounters a flame holder. If the flow rate is equal to the flame speed, we would expect a stationary flat flame front normal to the flow direction. If the flow rate is above the flame speed, the flame front will become conical such that the component of the velocity vector normal to the flame front is equal to the flame speed.

Stagnation flame

Here, the pre-mixed gases flow in such a way so as to form a region of stagnation (zero velocity) where the flame may be stabilized.

Spherical flame

In this configuration, the flame is typically initiated by way of a spark within a homogeneous pre-mixture. The subsequent propagation of the developed premixed flame occurs as a spherical front until the mixture is transformed entirely or the walls of the combustion vessel are reached.

Applications

Since the equivalence ratio of the premixed gases may be controlled, premixed combustion offers a means to attain low temperatures and, thereby, reduce NOx emissions. Due to improved mixing in comparison with diffusion flames, soot formation is mitigated as well. Premixed combustion has therefore gained significance in recent times. The uses involve lean-premixed-prevaporized (LPP) gas turbines and SI engines.

See also

References

  1. Lewis, Bernard; Elbe, Guenther von (2012). Combustion, Flames and Explosions of Gases. Elsevier. ISBN 9780323138024. https://books.google.com/books?id=xZL3V0w45aQC&q=Academic+Press+combustion+flames+and+explosions+of+gases&pg=PP1. 
  2. Williams, F. A. (2018). Combustion theory. CRC Press.
  3. Linan, A., & Williams, F. A. (1993). Fundamental aspects of combustion.
  4. MITANI, T. (1980). Propagation velocities of two-reactant flames. Combustion Science and Technology, 21(3-4), 175-177.
  5. Rogg, B., & Williams, F. A. (1985). Asymptotic analysis of laminar flame propagation with variable transport coefficients. Combustion science and technology, 42(5-6), 301-316.
  6. Chelliah, H. K., & Williams, F. A. (1987). Asymptotic analysis of two-reactant flames with variable properties and Stefan-Maxwell transport. Combustion science and technology, 51(4-6), 129-144.
  7. Rogg, B. (1986). On the accuracy of asymptotic flame speed predictions for two-reactant flames. Combustion science and technology, 45(5-6), 317-329.
  8. Clavin, P., & Graña-Otero, J. C. (2011). Curved and stretched flames: the two Markstein numbers. Journal of Fluid Mechanics, 686, 187-217.
  9. Clavin, Paul, and Geoff Searby. Combustion Waves and Fronts in Flows: Flames, Shocks, Detonations, Ablation Fronts and Explosion of Stars. Cambridge University Press, 2016.
  10. Peters, Norbert (2000). Turbulent combustion. Cambridge University Press. ISBN 9780511612701. OCLC 56066895. 
  11. Williams, F. A. (1985). Turbulent combustion. In The mathematics of combustion (pp. 97-131). Society for Industrial and Applied Mathematics.
  12. Kerstein, Alan R. (1988-01-01). "Field equation for interface propagation in an unsteady homogeneous flow field". Physical Review A 37 (7): 2728–2731. doi:10.1103/PhysRevA.37.2728. PMID 9899999.