Etude de la combustion de gaz de synthèse issus d'un processus de ...

Etude de la combustion de gaz de synthèse issus d'un processus de ... Etude de la combustion de gaz de synthèse issus d'un processus de ...

tel.archives.ouvertes.fr
from tel.archives.ouvertes.fr More from this publisher
27.12.2013 Views

Chapter 2 fundamental processes such as ignition, NO, and soot formation, and flame quenching. Moreover, some turbulent flame models prescribe the turbulent burning velocity as a function of laminar burning velocity. Thus, detailed information describing the dependence of the laminar burning velocity, flame thickness, ignition temperature, heat release rate and flame quenching on various system parameters can be a valuable diagnostic and design aid. tel-00623090, version 1 - 13 Sep 2011 There is a significant discrepancy in measuring burning velocities, which gives an indication of the difficulties and uncertainties associated with experimental determination of flame properties. In the light of the earlier experimental studies of Markstein, (1964), the asymptotic analysis of Klimov, (1963), and the computations of laminar flame structure with detailed chemical kinetics by several researchers, all of which show the importance of the flame stretch rate (Dixon-Lewis, 1991). There can be little doubt that this is often neglected key variable (Law, 1989). It follows that any experimental or computed value of laminar burning velocity should be associated with a value of the flame stretch rate. Ideally, the stretch-free value of the burning velocity should be quoted and the influence of stretch rate upon this value should be indicated by the value of the appropriated Markstein length. For these reasons this chapter begins with the definition of stretch rate and the corresponding Karlovitz and Markstein numbers. Following, the theory evolved with the burning velocity determination are described with emphasis for the constant volume and constant pressure methods due to be extensively used. This part of the chapter ends with the flammability limits description as is another important parameter of premixed laminar flames. 2.5.1 Flame stretch A flame surface propagating in a uniform flow field is submitted to strain and curvature effects leading to changes in the frontal area. Karlovitz et al., (1953) and Markstein, (1964) initiated the study of stretched premixed flames and demonstrated the importance of the aerodynamic stretching and the preferential diffusion on the flame response in terms of flame front instability. The flame stretch factor (κ) is defined as the relative rate of change of flame surface area (A) (Williams, 1985): 41

Bibliographic revision 1 d( δ A) 1 dA κ = = (2.2) δ A dt A dt The effect of stretch on the flame is to reduce the thickness of the flame front and hence the flame speed and influence the flame structure through its coupled effect with mass and heat diffusion. The concept of flame stretch can be applied to laminar flame speed; flame stabilization; flammability limits; and modeling of turbulent flames. Let’s first derive the basic relationship between the stretch rate and the strain rate, dilatation of the fluid element, and curvature of the flame surface. The three perpendicular coordinates on a curved flame surface are shown in Fig. 2.7, which has two unit vectors (ν and η ) tangent to the flame surface and an outward normal unit vector n , at the spatial point r( νη , , n) as a function of the three independent coordinates. tel-00623090, version 1 - 13 Sep 2011 Figure 2.7 - Curved laminar flame front with three perpendicular curvilinear coordinates. The elemental arc ( ds ) ν are given by: in the directionν and the elemental arc ( ) ⎛δr ⎞ ⎛δr ⎞ ds = d ds = dη ν ⎜ η ⎜ ⎟ δν ⎟ ⎝ ⎠ ⎝δη ⎠ ( ) ν , ( ) n η ν ds η in the direction η (2.3) The elemental flame surface can be calculated by: ⎛δr δr ⎞ dA t ⎜ ⎟ n d d ⎝δν δη ⎠ () = × ⋅ ( ν )( η ) (2.4) In the orthogonal curvilinear coordinates system, the two unit vectors e , e can be given by: δr δν δr δη e = , e and e e n ν = × = δr δν η δr δη ν η (2.5) ν η Thus: 42

Bibliographic revision<br />

1 d( δ A)<br />

1 dA<br />

κ = = (2.2)<br />

δ A dt A dt<br />

The effect of stretch on the f<strong>la</strong>me is to reduce the thickness of the f<strong>la</strong>me front and<br />

hence the f<strong>la</strong>me speed and influence the f<strong>la</strong>me structure through its coupled effect with<br />

mass and heat diffusion. The concept of f<strong>la</strong>me stretch can be applied to <strong>la</strong>minar f<strong>la</strong>me<br />

speed; f<strong>la</strong>me stabilization; f<strong>la</strong>mmability limits; and mo<strong>de</strong>ling of turbulent f<strong>la</strong>mes.<br />

Let’s first <strong>de</strong>rive the basic re<strong>la</strong>tionship between the stretch rate and the strain rate,<br />

di<strong>la</strong>tation of the fluid element, and curvature of the f<strong>la</strong>me surface. The three<br />

perpendicu<strong>la</strong>r coordinates on a curved f<strong>la</strong>me surface are shown in Fig. 2.7, which has<br />

two unit vectors (ν and η ) tangent to the f<strong>la</strong>me surface and an outward normal unit<br />

<br />

vector n , at the spatial point r( νη , , n)<br />

as a function of the three in<strong>de</strong>pen<strong>de</strong>nt coordinates.<br />

tel-00623090, version 1 - 13 Sep 2011<br />

Figure 2.7 - Curved <strong>la</strong>minar f<strong>la</strong>me front with three perpendicu<strong>la</strong>r curvilinear coordinates.<br />

<br />

The elemental arc ( ds ) ν<br />

are given by:<br />

in the directionν and the elemental arc ( )<br />

<br />

<br />

⎛δr<br />

⎞<br />

⎛δr<br />

⎞<br />

ds = d ds = dη<br />

ν ⎜<br />

η ⎜ ⎟<br />

δν ⎟<br />

⎝ ⎠ ⎝δη<br />

⎠<br />

( ) ν , ( )<br />

n <br />

η ν <br />

<br />

ds η<br />

in the direction η<br />

(2.3)<br />

The elemental f<strong>la</strong>me surface can be calcu<strong>la</strong>ted by:<br />

<br />

⎛δr<br />

δr<br />

⎞ <br />

dA t ⎜ ⎟ n d d<br />

⎝δν<br />

δη ⎠<br />

() = × ⋅ ( ν )( η )<br />

(2.4)<br />

<br />

In the orthogonal curvilinear coordinates system, the two unit vectors e , e can be<br />

given by:<br />

δr<br />

<br />

<br />

δν δr<br />

δη <br />

e = , e and e e n<br />

ν = × =<br />

δr<br />

δν η <br />

δr<br />

δη<br />

ν η<br />

(2.5)<br />

ν<br />

η<br />

Thus:<br />

42

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!