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

motions within the cylin<strong>de</strong>r, which in turn <strong>de</strong>termines the extent of mixing between the<br />

fresh charge and the residuals, as well as internal and external heat transfer rates. The<br />

key to the premixed <strong>combustion</strong> mo<strong>de</strong>ling is the prediction of S te, the turbulent f<strong>la</strong>me<br />

speed normal to the surface of the f<strong>la</strong>me. In turbulent f<strong>la</strong>mes, the f<strong>la</strong>me speed <strong>de</strong>pends<br />

on both chemical kinetics and the local turbulence characteristics.<br />

Many methods for <strong>de</strong>scribing and calcu<strong>la</strong>ting the turbulent f<strong>la</strong>me speed have been<br />

<strong>de</strong>veloped (see for instance the excellent review of Lipatnikov and Chomiak, (2002)).<br />

The goal of this work is to <strong>de</strong>velop a fast simu<strong>la</strong>tion program for the <strong>combustion</strong> of<br />

syngas in spark ignition engines. The main interest is the pressure <strong>de</strong>velopment in the<br />

engine cylin<strong>de</strong>rs, which is directly re<strong>la</strong>ted to the power output and the efficiency.<br />

Therefore, in this work the so-called DamkÖhler method is used and according to this<br />

mo<strong>de</strong>l turbulent f<strong>la</strong>me speed is as follows (Blizard and Keck, 1974):<br />

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

S = C u + S<br />

(6.26)<br />

te<br />

2<br />

'<br />

Where, u’ is the root mean square (rms) turbulent velocity, C 2 a calibration constant<br />

<strong>de</strong>pen<strong>de</strong>nt of the engine geometry and S u the <strong>la</strong>minar burning velocity.<br />

Obviously, proper in-cylin<strong>de</strong>r turbulence mo<strong>de</strong>ling needs to be estimated. For this<br />

propose, a simple turbulence mo<strong>de</strong>l, firstly proposed by Hall and Bracco, (1987) and<br />

used by several authors [Verhelst and Sierens, (2007); Farhad et al., (2009); Fe<strong>de</strong>rico<br />

et al., (2010)] has been consi<strong>de</strong>red:<br />

⎛ θ − 360 ⎞<br />

u' TDC<br />

= 0.75up = 0.75(2 sn), u' = u' TDC ⎜1−0.5<br />

45<br />

⎟<br />

(6.27)<br />

⎝<br />

⎠<br />

u<br />

where u’ TDC is the rms turbulent velocity at TDC, taken to be 0.75 times the mean<br />

piston speed; θ is the crank angle and, s, is the stoke. A linear <strong>de</strong>cay of the rms<br />

turbulent velocity u’ from top <strong>de</strong>ad center is imposed.<br />

6.2. Numerical solution procedure<br />

The basic concept of the mo<strong>de</strong>l is the division of the burned gas region into several<br />

distinct zones for taking into account the temperature stratification of the burned gas.<br />

The multi-zone simu<strong>la</strong>tion mo<strong>de</strong>l is applied throughout the closed part of the engine<br />

cycle, between IVC and EVO, i.e. compression, <strong>combustion</strong> and expansion. Admission<br />

phase is also inclu<strong>de</strong>d in the co<strong>de</strong> in or<strong>de</strong>r to take into account the heating of the<br />

175

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

Saved successfully!

Ooh no, something went wrong!