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

tel.archives.ouvertes.fr
from tel.archives.ouvertes.fr More from this publisher
27.12.2013 Views

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

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

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

Saved successfully!

Ooh no, something went wrong!