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 ...
Numerical simulation of a syngas-fuelled engine The rate of internal energy change for a mixture is given by: dh dm dT = u + mc dθ dθ dθ ∑ ∑ (6.15) i i i i vi i where m i is the mass of species i (O 2 , N 2 , CO 2 ,H 2 O, N, NO, OH, H, O, etc.) and c v is the specific heat under constant volume (a function of temperature only c v =du/dT), with ⎡⎛ n 1 ⎞ ⎤ cvi( T) = Rsi ⎢⎜ ∑ ai, nT − ⎟−1⎥ ⎣⎝ n ⎠ ⎦ (6.16) with the values of m i , dm i , T, dT found from the corresponding first-law analysis of the cylinder contents. The rate of entropy change is: tel-00623090, version 1 - 13 Sep 2011 dS dm m dT V dp = + − dθ dθ T dθ T dθ i i ∑ Si ( T, xip) ∑ cpi (6.17) i i With ' ⎛xi p ⎞ Si ( T, xip) = Si ( T, pi) −Rsi ln⎜ ⎟ (6.18) ⎝ pi ⎠ ' and S ( T, p ) i i the standard state entropy of species i, which is a function of temperature only, with x i the molar fraction of species i in the mixture (Ferguson, 1986; Heywood, 1988), given by the following property relation: 5 n−1 ⎡ ' ⎛ T ⎞ ⎤ si( T, pi) = Rsi ⎢ai1lnT ⎜ ∑ ai, n ⎟+ ai7⎥ ⎣ ⎝ n= 2 n −1⎠ ⎦ (6.19) For the Gibbs free enthalpy or energy: dG dm = ∑ (6.20) i i dθ i dθ μ where μ i =g i (T,p i ) is the chemical potential of species i in the mixture, with gi( T, pi) = gi( T, xip) = hi( T) −Tsi( T, xip) ⎡ ⎛xi p ⎞⎤ = hi( T) −T ⎢si( T, pi) −Rsi ln⎜ ⎟⎥ ⎢⎣ ⎝ pi ⎠⎥⎦ (6.21) 172
Chapter 6 For all the above expressions, it is assumed that the unburned mixture is frozen in composition and the burned mixture is always in equilibrium. Finally, the well-known ideal gas relation is given by: pV = mR T (6.22) s 6.1.3 Heat Transfer The heat transfer from gas to the walls is formulated as: ( ) Q = h A T −T (6.23) g g w tel-00623090, version 1 - 13 Sep 2011 where h g is the heat transfer coefficient, A is the area in contact with the gas, T g is the gas temperature and T w is the wall temperature. In the single zone analysis, the heat transfer coefficient is the same for all surfaces in the cylinder. In general, a classic global heat transfer model is applied to calculate the heat transfer coefficient and an area-averaged heat transfer rate. Several correlations for calculating the heat transfer coefficient in SI and CI engines have been published in the literature. These studies have generally relied on dimensional analysis for turbulent flow that correlates the Nusselt, Reynolds, and Prandtl numbers. Using experiments in spherical vessels or engines and applying the assumption of quasi-steady conditions has led to empirical correlations for both SI and CI engine heat transfer. These correlations provide a heat transfer coefficient representing a spatially-averaged value for the cylinder. Hence, they are commonly referred to as global heat transfer models, e.g. Woschni (1967), Annand (1963), or Hohenberg (1979). In this code one applies the classical Woschni’s correlation. The Woschni heat transfer correlation is given as: h () t a B P() t T() t v() t g −02 . 08 . −055 . 08 . = s (6.24) where a s is a scaling factor used for tuning of the coefficient to match a specific engine geometry calculated and used by Hohenberg, (1979) as 130, B is the bore (m), P and T are the instantaneous cylinder pressure (bar) and gas temperature (K), respectively. The instantaneous characteristic velocity, v is defined as: 173
- Page 125 and 126: Chapter 4 Notice the similar behavi
- Page 127 and 128: Chapter 4 A very good agreement bet
- Page 129 and 130: Chapter 4 ( ) Q = h T − T (4.21)
- Page 131 and 132: Chapter 4 tel-00623090, version 1 -
- Page 133 and 134: Chapter 4 are tested and discussed.
- Page 135 and 136: Chapter 4 7 500 Pressure (bar) 6 5
- Page 137 and 138: Chapter 4 Pressure (bar) 7 6 5 4 3
- Page 139 and 140: Chapter 4 4.2.3.4 Quenching distanc
- Page 141 and 142: Chapter 4 10000 Quenching distance
- Page 143 and 144: Chapter 5 CHAPTER 5 EXPERIMENTAL ST
- Page 145 and 146: Chapter 5 30 10 25 8 Pressure (bar)
- Page 147 and 148: Chapter 5 30 Piston position (mm) 2
- Page 149 and 150: Chapter 5 5.1.1.4 In-cylinder press
- Page 151 and 152: Chapter 5 estimation of various par
- Page 153 and 154: Chapter 5 TDC 1.25 ms 2.5 ms 3.75 m
- Page 155 and 156: Chapter 5 Piston position (mm) 500
- Page 157 and 158: Chapter 5 tel-00623090, version 1 -
- Page 159 and 160: Chapter 5 From figure 5.15 is possi
- Page 161 and 162: Chapter 5 From figure 5.16 is obser
- Page 163 and 164: Chapter 5 80 Pressure (bar) 70 60 5
- Page 165 and 166: Chapter 5 80 10 Pmax (bar) 70 60 50
- Page 167 and 168: Chapter 5 -5.0 ms -3.75 ms -2.5 ms
- Page 169 and 170: Chapter 5 observation emphasis the
- Page 171 and 172: Chapter 6 CHAPTER 6 NUMERICAL SIMUL
- Page 173 and 174: Chapter 6 centered at the spark plu
- Page 175: Chapter 6 H 2 O, (3) N 2 , (4) O 2
- Page 179 and 180: Chapter 6 motions within the cylind
- Page 181 and 182: Chapter 6 tel-00623090, version 1 -
- Page 183 and 184: Chapter 6 Heat transfer Wei et al.,
- Page 185 and 186: Chapter 6 The calibration coefficie
- Page 187 and 188: Chapter 6 6.3.2.2 In-cylinder volum
- Page 189 and 190: Chapter 6 40 Experimental 30 Numeri
- Page 191 and 192: Chapter 6 80 70 60 Numerical Experi
- Page 193 and 194: Chapter 6 downdraft syngas than for
- Page 195 and 196: Chapter 6 80 23º BTDC 60 29.5º BT
- Page 197 and 198: Chapter 6 with experimental results
- Page 199 and 200: Conclusions CHAPTER 7 CONCLUSIONS 7
- Page 201 and 202: Conclusions radius and time for syn
- Page 203 and 204: Conclusions conditions, therefore s
- Page 205 and 206: Conclusions tel-00623090, version 1
- Page 207 and 208: References References tel-00623090,
- Page 209 and 210: References tel-00623090, version 1
- Page 211 and 212: References tel-00623090, version 1
- Page 213 and 214: References tel-00623090, version 1
- Page 215 and 216: References tel-00623090, version 1
- Page 217 and 218: References tel-00623090, version 1
- Page 219 and 220: References tel-00623090, version 1
- Page 221 and 222: Appendix A - Overdetermined linear
- Page 223 and 224: Appendix A - Overdetermined linear
- Page 225 and 226: Appendix B- Syngas-air mixtures pro
Chapter 6<br />
For all the above expressions, it is assumed that the unburned mixture is frozen in<br />
composition and the burned mixture is always in equilibrium. Finally, the well-known<br />
i<strong>de</strong>al gas re<strong>la</strong>tion is given by:<br />
pV<br />
= mR T<br />
(6.22)<br />
s<br />
6.1.3 Heat Transfer<br />
The heat transfer from gas to the walls is formu<strong>la</strong>ted as:<br />
( )<br />
Q<br />
= h A T −T (6.23)<br />
g g w<br />
tel-00623090, version 1 - 13 Sep 2011<br />
where h g is the heat transfer coefficient, A is the area in contact with the gas, T g is the<br />
gas temperature and T w is the wall temperature. In the single zone analysis, the heat<br />
transfer coefficient is the same for all surfaces in the cylin<strong>de</strong>r. In general, a c<strong>la</strong>ssic<br />
global heat transfer mo<strong>de</strong>l is applied to calcu<strong>la</strong>te the heat transfer coefficient and an<br />
area-averaged heat transfer rate.<br />
Several corre<strong>la</strong>tions for calcu<strong>la</strong>ting the heat transfer coefficient in SI and CI engines<br />
have been published in the literature. These studies have generally relied on<br />
dimensional analysis for turbulent flow that corre<strong>la</strong>tes the Nusselt, Reynolds, and<br />
Prandtl numbers. Using experiments in spherical vessels or engines and applying the<br />
assumption of quasi-steady conditions has led to empirical corre<strong>la</strong>tions for both SI and<br />
CI engine heat transfer. These corre<strong>la</strong>tions provi<strong>de</strong> a heat transfer coefficient<br />
representing a spatially-averaged value for the cylin<strong>de</strong>r. Hence, they are commonly<br />
referred to as global heat transfer mo<strong>de</strong>ls, e.g. Woschni (1967), Annand (1963), or<br />
Hohenberg (1979). In this co<strong>de</strong> one applies the c<strong>la</strong>ssical Woschni’s corre<strong>la</strong>tion.<br />
The Woschni heat transfer corre<strong>la</strong>tion is given as:<br />
h () t a B P() t T() t v()<br />
t<br />
g<br />
−02 . 08 . −055 . 08 .<br />
=<br />
s<br />
(6.24)<br />
where a s is a scaling factor used for tuning of the coefficient to match a specific engine<br />
geometry calcu<strong>la</strong>ted and used by Hohenberg, (1979) as 130, B is the bore (m), P and T<br />
are the instantaneous cylin<strong>de</strong>r pressure (bar) and gas temperature (K), respectively.<br />
The instantaneous characteristic velocity, v is <strong>de</strong>fined as:<br />
173