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 Where, r c represents the local radius of the flame curvature. The first term of Eq. (2.21) represents the effect of the strain; the second term represents the effect of the fluid expansion of fluid (dilatation) and, the third term, the effect of the flame curvature. It is evident that the flame can be stretched by the combined effect of strain, volume expansion of the fluid and the curvature of the flame, which arises from the nonuniformities of the flow and the normal propagation of the flame front. For the present purposes, an appropriate unified tensor expression, in terms of strain rate, κ s , and the stretch rate due to flame curvature, κ c , is that of Candel and Poinsot, (1990). κ = κ + κ (2.22) s c tel-00623090, version 1 - 13 Sep 2011 With κ s =−nn: ∇ S+∇. S (2.23) κ = S ∇. n (2.24) c u In spherical expanding flames it is convenient to use spherical coordinates (r, θ, φ), the components of n and S are written as (n r , n θ , n φ ) and (s r , s θ , s φ ), respectively. Then: ⎡ 2⎛∂sr ⎞ 2⎛1 ∂s 2 1 s θ sr ⎞ ⎛ ∂ φ sr sθ cotθ ⎞ κs =− ⎢nr ⎜ n n r ⎟+ θ ⎜ + φ r θ r ⎟+ ⎜ + + ⎟ ⎣ ⎝ ∂ ⎠ ⎝ ∂ ⎠ ⎝r sinφ ∂φ r r ⎠ ⎛∂s 1 s s r s φ 1 s s θ ∂ θ ⎞ ⎛∂ ∂ r φ ⎞ + nn r θ ⎜ + − nn r φ r r θ r ⎟+ ⎜ + − ⎟ ⎝ ∂ ∂ ⎠ ⎝ ∂r r sinθ ∂φ r ⎠ ⎛1 ∂sφ 1 ∂sθ cotθ ⎞⎤ + nn θ φ ⎜ + − sφ ⎟⎥ ⎝r ∂θ r sinθ ∂φ r ⎠⎦⎥ 2 1 ∂ ( r sr ) 1 ∂ ( sθ sinθ ) 1 ∂sφ + + + 2 r ∂r r sinθ ∂θ r sinθ ∂ φ (2.25) 2 ( r nr ) ∂ ( n sinθ ) ⎡ 1 ∂ 1 1 ∂n ⎤ θ φ κc = un + ⎢ + + ⎥ 2 ⎢r ∂r r sinθ ∂θ r sinθ ∂φ ⎥ ⎣ ⎦ (2.26) For an outward spherically propagating flame, the flame surface is identified by the cold front of radius r u , n r =1, n θ =n φ =0, and s r =s g , s θ =s φ =0 (Bradley et al., 1996). The burning velocity, s u , is associated with this surface and the gas velocity ahead of it is s g . The flame speed, dr u /dt, is equal to s g + n r s u , and is indicated by S n . Applying this conditions to Eqs. 2.25 and 2.26 gives 45

Bibliographic revision 2 ( rsr ) 2 ⎛∂s ⎞ 1 ∂ κs = nr ⎜ + = ⎝ ∂ s r g 2 2 r ⎟ ⎠ r ∂r ru 2 ( r nr ) 1 ∂ s κc = un = 2n 2 r r ∂r r u u (2.27) (2.28) and the total stretch rate is: s κ = (2.29) 2 n ru The outwardly propagating flame, ignited from a central ignition point, is the most common spherical flame and because the stretch rates to which it is subjected are well defined, it is well suited to burning velocity measurements. tel-00623090, version 1 - 13 Sep 2011 2.5.1.1 Karlovitz number The Karlovitz number (Ka) is defined as the nondimensional stretch factor, using the thickness of the unstretched flame (δ u0 ) and the normal unstretched laminar flame velocity (S u0 ) to form a reference time (Kuo, 2005): Ka δ u0 = = (2.30) S u0 κ Residence time for crossing an unstretched flame Characteristic time for flame stretching For convenience, the Karlovitz number can be written as a sum of two parts due to the contribution from strain and curvature as follows: Ka = Ka + Ka (2.31) s c Where δ δ δ Ka = SS + : ∇ S ; Ka = S ∇ n = ( ηη) u 0 u 0 u 0 u s c u Su0 Su0 R Su0 S (2.32) 2.5.1.2 Markstein number Flame stretch can change the burning velocity of premixed flames significantly. It is important to know the relation between these two parameters, because the burning velocity is a crucial parameter in premixed combustion processes. Markstein, (1964) was the first to propose a phenomenological relation between the laminar burning velocity and the curvature of the flame front. Notice that the straining of the flow is not 46

Chapter 2<br />

Where, r c represents the local radius of the f<strong>la</strong>me curvature. The first term of Eq. (2.21)<br />

represents the effect of the strain; the second term represents the effect of the fluid<br />

expansion of fluid (di<strong>la</strong>tation) and, the third term, the effect of the f<strong>la</strong>me curvature. It is<br />

evi<strong>de</strong>nt that the f<strong>la</strong>me can be stretched by the combined effect of strain, volume<br />

expansion of the fluid and the curvature of the f<strong>la</strong>me, which arises from the<br />

nonuniformities of the flow and the normal propagation of the f<strong>la</strong>me front.<br />

For the present purposes, an appropriate unified tensor expression, in terms of strain<br />

rate, κ s , and the stretch rate due to f<strong>la</strong>me curvature, κ c , is that of Can<strong>de</strong>l and Poinsot,<br />

(1990).<br />

κ = κ + κ<br />

(2.22)<br />

s<br />

c<br />

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

With<br />

<br />

κ<br />

s<br />

=−nn: ∇ S+∇.<br />

S<br />

(2.23)<br />

κ = S ∇.<br />

n <br />

(2.24)<br />

c<br />

u<br />

In spherical expanding f<strong>la</strong>mes it is convenient to use spherical coordinates (r, θ, φ), the<br />

components of n and S are written as (n r , n θ , n φ ) and (s r , s θ , s φ ), respectively. Then:<br />

⎡<br />

2⎛∂sr ⎞ 2⎛1 ∂s<br />

2 1 s<br />

θ<br />

sr ⎞ ⎛ ∂<br />

φ sr<br />

sθ<br />

cotθ<br />

⎞<br />

κs<br />

=− ⎢nr<br />

⎜ n n<br />

r<br />

⎟+ θ ⎜ +<br />

φ<br />

r θ r<br />

⎟+ ⎜ + + ⎟<br />

⎣ ⎝ ∂ ⎠ ⎝ ∂ ⎠ ⎝r sinφ ∂φ<br />

r r ⎠<br />

⎛∂s<br />

1 s<br />

s<br />

r<br />

s<br />

φ 1 s s<br />

θ<br />

∂<br />

θ ⎞ ⎛∂<br />

∂<br />

r φ ⎞<br />

+ nn<br />

r θ ⎜ + − nn<br />

r φ<br />

r r θ r<br />

⎟+ ⎜ + − ⎟<br />

⎝ ∂ ∂ ⎠ ⎝ ∂r r sinθ ∂φ<br />

r ⎠<br />

⎛1 ∂sφ<br />

1 ∂sθ<br />

cotθ<br />

⎞⎤<br />

+ nn<br />

θ φ ⎜ + − sφ<br />

⎟⎥<br />

⎝r ∂θ r sinθ ∂φ<br />

r ⎠⎦⎥<br />

2<br />

1 ∂ ( r sr<br />

) 1 ∂ ( sθ<br />

sinθ<br />

) 1 ∂sφ<br />

+ +<br />

+<br />

2<br />

r ∂r<br />

r sinθ<br />

∂θ r sinθ ∂ φ<br />

(2.25)<br />

2<br />

( r nr<br />

) ∂ ( n sinθ<br />

)<br />

⎡ 1 ∂ 1 1 ∂n<br />

⎤<br />

θ<br />

φ<br />

κc<br />

= un<br />

+ ⎢ + + ⎥<br />

2<br />

⎢r<br />

∂r r sinθ ∂θ r sinθ ∂φ<br />

⎥<br />

⎣<br />

⎦<br />

(2.26)<br />

For an outward spherically propagating f<strong>la</strong>me, the f<strong>la</strong>me surface is i<strong>de</strong>ntified by the cold<br />

front of radius r u , n r =1, n θ =n φ =0, and s r =s g , s θ =s φ =0 (Bradley et al., 1996). The burning<br />

velocity, s u , is associated with this surface and the gas velocity ahead of it is s g . The<br />

f<strong>la</strong>me speed, dr u /dt, is equal to s g + n r s u , and is indicated by S n . Applying this conditions<br />

to Eqs. 2.25 and 2.26 gives<br />

45

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

Saved successfully!

Ooh no, something went wrong!