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 ...
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
- Page 1 and 2: THÈSE Pour l’obtention du Grade
- Page 3 and 4: Acknowledgements Acknowledgements T
- Page 5 and 6: Résumé __________________________
- Page 7 and 8: Nomenclature Nomenclature Roman tel
- Page 9 and 10: Nomenclature Subscripts tel-0062309
- Page 11 and 12: Contents tel-00623090, version 1 -
- Page 13 and 14: Contents 6.4. SYNGAS FUELLED-ENGINE
- Page 15 and 16: Introduction CHAPTER 1 INTRODUCTION
- Page 17 and 18: Introduction proves to have higher
- Page 19 and 20: Introduction Chapter 3 - Experiment
- Page 21 and 22: Bibliographic revision CHAPTER 2 BI
- Page 23 and 24: Bibliographic revision point today
- Page 25 and 26: Bibliographic revision - Boudouard
- Page 27 and 28: Bibliographic revision Table 2.1 -
- Page 29 and 30: Bibliographic revision Biomass Dryi
- Page 31 and 32: Bibliographic revision Circulating
- Page 33 and 34: Bibliographic revision or eliminate
- Page 35 and 36: Bibliographic revision established
- Page 37 and 38: Bibliographic revision Hydrogen Hyd
- Page 39 and 40: Bibliographic revision of low moist
- Page 41 and 42: Bibliographic revision scrubbing an
- Page 43 and 44: Bibliographic revision suggests tha
- Page 45: Bibliographic revision 1 d( δ A) 1
- Page 49 and 50: Bibliographic revision 2 ( rsr ) 2
- Page 51 and 52: Bibliographic revision This evoluti
- Page 53 and 54: Bibliographic revision The burning
- Page 55 and 56: Bibliographic revision δVG = − a
- Page 57 and 58: Bibliographic revision 2 1 − −
- Page 59 and 60: Bibliographic revision where the su
- Page 61 and 62: Bibliographic revision the stretche
- Page 63 and 64: Bibliographic revision burning velo
- Page 65 and 66: Experimental set ups and diagnostic
- Page 67 and 68: Experimental set ups and diagnostic
- Page 69 and 70: Experimental set ups and diagnostic
- Page 71 and 72: Experimental set ups and diagnostic
- Page 73 and 74: Experimental set ups and diagnostic
- Page 75 and 76: Experimental set ups and diagnostic
- Page 77 and 78: Experimental set ups and diagnostic
- Page 79 and 80: Experimental set ups and diagnostic
- Page 81 and 82: Experimental set ups and diagnostic
- Page 83 and 84: Experimental set ups and diagnostic
- Page 85 and 86: Experimental set ups and diagnostic
- Page 87 and 88: Experimental set ups and diagnostic
- Page 89 and 90: Chapter 4 CHAPTER 4 EXPERIMENTAL AN
- Page 91 and 92: Chapter 4 4.1 Laminar burning veloc
- Page 93 and 94: Chapter 4 4.1.1.1 Flame morphology
- Page 95 and 96: Chapter 4 P i = 1.0 bar, Ti = 293 K
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