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 δr δr dA t d d δν δη () = ( ν )( η ) (2.6) * * Now let another surface that is close to the flame surface be represented by r ( νη , , n ) such that: , * * ∂r ∂r ∂a ∂r ∂r ∂a = + , = + , = + ∂ν ∂ν ∂ν ∂η ∂η ∂η * r r a (2.7) Where a is a small magnitude displacement vector; then: tel-00623090, version 1 - 13 Sep 2011 ∂ Where ∇ t = eν + e ∂ν * * * * * ⎛δr δr ⎞ * δr δr dA t ⎜ ⎟ n d d d d ⎝ δν δη ⎠ δν δη δr δr ≅ ⎡1+∇ t. a⎤ d d = ⎡1+∇t. a⎤dA t δν δη ⎣ ⎦ ⎣ ⎦ () = × ⋅ ( ν)( η) ( ν )( η ) η ( ν)( η) ( ) (2.8) ∂ , which represents the gradient operator along the tangential ∂η plane of the flame surface (Chung and Law, 1988). It is useful to note that we can always decompose any arbitrary velocity into two components: a tangential to the flame surface and another normal to the flame surface as follows: V = V + V = V n n+ V ( . ) n t t (2.9) Now considering a curved flame surface A(t) moving in the space with local velocity W (note that each spatial point has its own velocity). Then, the rate of change of the flux of a vector G across the flame surface is given by the following expression based on the Reynolds’ transport theorem. d dt ∫ . ⎡∂G GndA W . G G . W G . ⎤ = ∫ ⎢ + ∇ − ∇ + ∇W⎥ . ndA ⎢⎣ ∂t ⎥⎦ At ( ) At ( ) (2.10) By specifying G = n , Eq. (2.10) yields: d dt ∫ ⎡∂n ⎤ dA = ∫ ⎢ + W. ∇n −n. ∇ W + n∇. W . n dA ∂t ⎥ ⎣ ⎦ At ( ) At ( ) (2.11) 43

Bibliographic revision Since n is a unit normal vector, we have: ( ) 2 ∂n 1 ∂n . n = = 0 ∂t 2 ∂t ∇ = ∇ ⎜ ⎟ = 0 ⎝ 2 ⎠ 2 n W. n . n W. ⎛ ⎞ (2.12) (2.13) Therefore, the transport equation (2.11) can be written as follows: d dt ∫ dA = ⎡−n n : ∇ W +∇. W ⎤ ∫ ⎣ ⎦ dA At ( ) At ( ) (2.14) Note that in tensor notation, the first term of the right side of the Eq. (2.14) can be written as: tel-00623090, version 1 - 13 Sep 2011 −nn: ∇ W = −n n i j ∂w ∂x j i (2.15) If we consider a surface element, then A(t) can be substituted by δΑ. The flame stretch rate κ can be expressed as: 1 d( δ A) κ = = −nn: ∇ W +∇. W (2.16) δ A dt The general velocity W can be considered to have two components: one in the normal direction at a speed S u and the other is the local fluid velocity. Thus: W = S + S n u (2.17) Several vector operations can be applied to Eq. (2.17) to give: ∇ W = ∇ S + Su ∇ n+ n∇Su nn : ∇ W = nn : ∇ S + n∇Su ∇ . W = ∇ . S + S ∇ . n+ n. ∇S u u (2.18) (2.19) (2.20) Substituting these expressions in the Eq. (2.16), we have: 1 d( δA) ∂Si 1 ∂ρ S κ = = −nn: ∇ S+∇ . S+ Su∇ . n = − ni nj +− + δA dt ∂x ρ ∂t r j u c (2.21) 44

Bibliographic revision<br />

Since n is a unit normal vector, we have:<br />

<br />

( )<br />

2<br />

∂n<br />

1 ∂n<br />

. n = = 0<br />

∂t<br />

2 ∂t<br />

<br />

∇ = ∇ ⎜ ⎟ = 0<br />

⎝ 2 ⎠<br />

2<br />

n<br />

W. n . n W.<br />

⎛ ⎞<br />

(2.12)<br />

(2.13)<br />

Therefore, the transport equation (2.11) can be written as follows:<br />

d<br />

dt<br />

∫<br />

<br />

dA = ⎡−n n : ∇ W +∇.<br />

W ⎤<br />

∫ ⎣<br />

⎦<br />

dA<br />

At ( ) At ( )<br />

(2.14)<br />

Note that in tensor notation, the first term of the right si<strong>de</strong> of the Eq. (2.14) can be<br />

written as:<br />

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

<br />

−nn:<br />

∇ W = −n n<br />

i<br />

j<br />

∂w<br />

∂x<br />

j<br />

i<br />

(2.15)<br />

If we consi<strong>de</strong>r a surface element, then A(t) can be substituted by δΑ. The f<strong>la</strong>me stretch<br />

rate κ can be expressed as:<br />

1 d( δ A)<br />

<br />

κ = = −nn: ∇ W +∇.<br />

W<br />

(2.16)<br />

δ A dt<br />

The general velocity W can be consi<strong>de</strong>red to have two components: one in the normal<br />

direction at a speed S u and the other is the local fluid velocity. Thus:<br />

<br />

W = S + S n<br />

u<br />

(2.17)<br />

Several vector operations can be applied to Eq. (2.17) to give:<br />

<br />

∇ W = ∇ S + Su<br />

∇ n+ n∇Su<br />

<br />

nn : ∇ W = nn : ∇ S + n∇Su<br />

<br />

∇ . W = ∇ . S + S ∇ . n+ n.<br />

∇S<br />

u<br />

u<br />

(2.18)<br />

(2.19)<br />

(2.20)<br />

Substituting these expressions in the Eq. (2.16), we have:<br />

1 d( δA)<br />

∂Si<br />

1 ∂ρ<br />

S<br />

κ = = −nn: ∇ S+∇ . S+ Su∇ . n = − ni nj<br />

+− +<br />

δA dt ∂x ρ ∂t r<br />

j<br />

u<br />

c<br />

(2.21)<br />

44

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

Saved successfully!

Ooh no, something went wrong!