27.12.2013 Views

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 ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 6<br />

centered at the spark plug (Verhelst and Sheppard, 2009). The multi-zone approach is<br />

retained throughout the expansion phase, i.e. from the end of <strong>combustion</strong> until the<br />

EVO.<br />

6.1.1 Conservation and state equations<br />

The governing equations are presented for the <strong>combustion</strong> process. During<br />

compression and expansion the same equations are retained, with adaptations<br />

regarding the number of zones inclu<strong>de</strong>d.<br />

Assuming that ‘n’ is the number of the current burned zones, with the <strong>la</strong>test-generated<br />

burned zone indicated by subscript ‘n’, the energy conservation equation is applied for<br />

the unburned zone, the n-1 already burned zones and the total cylin<strong>de</strong>r charge,<br />

yielding, respectively.<br />

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

dQu dUu dVu dmb<br />

= + p + hu<br />

(6.1)<br />

dθ dθ dθ dθ<br />

dQbi ,<br />

dUbi ,<br />

dVbi<br />

,<br />

= + p ( i = 1,..., n−1)<br />

(6.2)<br />

dθ dθ dθ<br />

dQ dU dV<br />

= + p<br />

(6.3)<br />

dθ dθ dθ<br />

Zero blow-by rate is assumed during the whole closed cycle period. The instantaneous<br />

cylin<strong>de</strong>r volume is equal to the sum of the volumes of the unburned and burned zones,<br />

resulting in the following re<strong>la</strong>tion after differentiating with respect to crank angle:<br />

dV dV dV<br />

dθ dθ dθ<br />

n<br />

= ∑<br />

bi ,<br />

+<br />

u<br />

(6.4)<br />

i = 1<br />

where V is the instantaneous cylin<strong>de</strong>r volume, calcu<strong>la</strong>ted from the engine geometric<br />

characteristics as a function of crank angle:<br />

2<br />

( sin )<br />

⎧ rc<br />

−1 ⎡ 1<br />

⎤⎫<br />

V = VTDC<br />

⎨1+ 1− cosθ + 1− 1−ϕ θ<br />

2<br />

⎢<br />

⎬<br />

ϕ<br />

⎥<br />

⎩ ⎣<br />

⎦⎭<br />

(6.5)<br />

with ϕ being the crank radius to piston rod length ratio. The rate of change of the<br />

instantaneous cylin<strong>de</strong>r volume against crank angle is<br />

169

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

Saved successfully!

Ooh no, something went wrong!