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 6.3.2 RCM As shown in chapter 5, RCM simulates a single engine cycle and some adaptations to the model presented above are needed in order to take into account the dynamic and geometrical specificness of the RCM. Three main aspects should be considered, the instantaneous volume, heat transfer and burning rate. Two distinctive cases are considered in the code: single compression, and compression and expansion strokes. 6.3.2.1 Flame propagation tel-00623090, version 1 - 13 Sep 2011 The flame development in the combustion chamber of a RCM was shown above for the compression and expansion strokes. In order two verify the assumption frequently used for engine combustion models, where the flame area is a spherical flame front truncated by the cylinder walls and the piston, centered at the spark plug adopted in this work, the flame images are used in an image treatment Matlab code developed by Strozzi, (2008), which allows visualizing the flame contours. Figures 6.5 show the flame contours in a RCM combustion chamber for downdraft syngas-air stoichiometric mixture. The window has 613×337 pixels with resolution of 2.835 pixels/mm. The frequency of the each contour is kept to 1.4 KHz. Figure 6.5 - RCM flame contours of stoichiometric downdraft syngas-air mixture with ignition at 12.5 ms BTDC. The earlier stage of flame kernel with influence of the spark plug was removed from the contours due to lack of definition. After around 5.0 ms (the duration of the spark), the flame contours are well defined and shows that the propagation is spherical. At this moment the piston is reaching TDC and the flame shows a remarkable slow down due to piston deadening, which is felt on the pressure signal by a decreasing gradient. After TDC, the flame contour shows to near double the displacement due to expansion. 182
Chapter 6 6.3.2.2 In-cylinder volume The instantaneous cylinder volume is calculated based on the piston position as a function of crank angle and fitted by polynomial functions in order to be implemented in the code. Thus, first one evaluates the error in the polynomial approach. Figure 6.6 shows the comparison between the experimental volume inside the cylinder and the corresponding six degree polynomial expression applied up to top dead center to better fitting. The remaining part is constant and equal to the clearance volume. The maximum error is about 0.5% (5 cm 3 ). 1200 1000 Experimental Polynomial tel-00623090, version 1 - 13 Sep 2011 Volume (cm 3 ) 800 600 400 200 0 180 210 240 270 300 330 360 390 420 450 480 510 540 Crank Angle (degrees) Figure 6.6 – In-cylinder volume polynomial fitting: single compression downdraft syngas case, ignition timing at 12.5 ms BTDC. Figure 6.7 shows in-cylinder volume variation during compression and expansion and the corresponding six degree polynomial. The maximum error is 1.5% (15 cm 3 ) obtained at TDC. 183
- 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 and 176: Chapter 6 H 2 O, (3) N 2 , (4) O 2
- Page 177 and 178: Chapter 6 For all the above express
- 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: Chapter 6 The calibration coefficie
- 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
- Page 227 and 228: Appendix C -Rivère model Heat flux
- Page 229 and 230: Appendix C -Rivère model The work
- Page 231 and 232: tel-00623090, version 1 - 13 Sep 20
Chapter 6<br />
6.3.2.2 In-cylin<strong>de</strong>r volume<br />
The instantaneous cylin<strong>de</strong>r volume is calcu<strong>la</strong>ted based on the piston position as a<br />
function of crank angle and fitted by polynomial functions in or<strong>de</strong>r to be implemented in<br />
the co<strong>de</strong>. Thus, first one evaluates the error in the polynomial approach.<br />
Figure 6.6 shows the comparison between the experimental volume insi<strong>de</strong> the cylin<strong>de</strong>r<br />
and the corresponding six <strong>de</strong>gree polynomial expression applied up to top <strong>de</strong>ad center<br />
to better fitting. The remaining part is constant and equal to the clearance volume. The<br />
maximum error is about 0.5% (5 cm 3 ).<br />
1200<br />
1000<br />
Experimental<br />
Polynomial<br />
tel-00623090, version 1 - 13 Sep 2011<br />
Volume (cm 3 )<br />
800<br />
600<br />
400<br />
200<br />
0<br />
180 210 240 270 300 330 360 390 420 450 480 510 540<br />
Crank Angle (<strong>de</strong>grees)<br />
Figure 6.6 – In-cylin<strong>de</strong>r volume polynomial fitting: single compression downdraft syngas case,<br />
ignition timing at 12.5 ms BTDC.<br />
Figure 6.7 shows in-cylin<strong>de</strong>r volume variation during compression and expansion and<br />
the corresponding six <strong>de</strong>gree polynomial. The maximum error is 1.5% (15 cm 3 )<br />
obtained at TDC.<br />
183