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

Experimental and numerical laminar syngas combustion Heat radiation of the burned gases to the chamber walls are modeled neglecting the radiation of the particles. This assumption is reasonable due to the high purity of the tested mixtures. During combustion, heat radiation can have origin on unburned gases, which contains: H 2 , CO, CH 4 , CO 2 , N 2 and O 2 , and also on the burned gases, which mainly contains: CO 2 , H 2 O and N 2 , being irrelevant the remaining combustion products. Only molecules that have a non null dipolar moment are susceptible to emit thermal radiation (Boust, 2006). Therefore only CO 2 , H 2 O and CH 4 are considered. In practical terms, the radiation of the unburned gases heated by compression is insignificant comparatively with the burned gases, which temperature is 6-7 times higher. Thus, it is assumed that only CO 2 and H 2 O radiates significantly. tel-00623090, version 1 - 13 Sep 2011 The radiation heat transfer is modeled by the Stefan’s law considering the burned gases as a grey body with uniform temperature T g and ε the apparent grey-body emissivity calculated from the contributions of H 2 O and CO 2 . As the spectral emissivities of these species are similar, the emissivity variation term Δε is included. ε = ε + ε −Δ ε (4.25) CO2 H2O At the end of combustion, only burned gases are inside the chamber. Then, the net superficial radiative flow Q ,r received by the wall, with absorption factor α, from the burned gases is given by the Stefan’s constant. 4 4 ( ) Q = αεσ T − T (4.26) r g w When the sphere of burned gases (radius r) does not occupy the entire chamber (radius R), the sphere surface ratio gives the radiative flow. 2 4π r Qr( r) = Q 2 r (4.27) 4π R The emissivity of H 2 O and CO 2 as well as the variation term are calculated using the correlation of Leckner, (1972). This correlation reproduces the gases temperature influence, the partial pressure of each species and the length of the average radius. 4.2.2 Calculation procedure In the multi-zone model, flame propagation is seen as the consecutive consumption of unburned mixture within the zones with an equal mass distribution between the zones 126

Chapter 4 tel-00623090, version 1 - 13 Sep 2011 in the spherical vessel (Fig. 4.37). Before ignition, the mass in the spherical vessel is divided into n zones. At the time when combustion has just begun in the bomb, the flame front will consume zone 1 first. As a result, the temperature and hence pressure of zone 1 will increase, thereby compressing the rest of the unburned gas (considered as single entity) and increasing the pressure inside the vessel to a higher value. After the consumption of the first zone, combustion of the second and subsequent zones will take place at a higher pressure than the initial pressure. At any instant of time when the flame front is passing through the n th zone, the combustion of this zone takes place at a temperature of T u , n−1 (> T i ) and a constant pressure P u,n−1 (> P i ). The combustion within a given zone takes place progressively. After the flame has consumed the n th zone, it is then assumed to be adiabatic. Subsequent combustion further compresses the burned gas and the unburned gas. As a result, temperature and density gradients are established in the burned gas region. At the end of combustion, the burned gas cooling is computed. Figure 4.37 – Radial distribution of the multiple zones inside a spherical vessel. Hatched portion indicates the position of the flame front at an instant of time. Figure 4.38 shows the flowchart of the adapted Fortran code. The output data of the code are: the burned gas temperature, flame radius and flame speed, pressure, as well as the wall thermal flux. At the end of combustion, an energy balance is made. 127

Chapter 4<br />

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

in the spherical vessel (Fig. 4.37). Before ignition, the mass in the spherical vessel is<br />

divi<strong>de</strong>d into n zones. At the time when <strong>combustion</strong> has just begun in the bomb, the<br />

f<strong>la</strong>me front will consume zone 1 first. As a result, the temperature and hence pressure<br />

of zone 1 will increase, thereby compressing the rest of the unburned gas (consi<strong>de</strong>red<br />

as single entity) and increasing the pressure insi<strong>de</strong> the vessel to a higher value. After<br />

the consumption of the first zone, <strong>combustion</strong> of the second and subsequent zones will<br />

take p<strong>la</strong>ce at a higher pressure than the initial pressure. At any instant of time when the<br />

f<strong>la</strong>me front is passing through the n th zone, the <strong>combustion</strong> of this zone takes p<strong>la</strong>ce at a<br />

temperature of T u , n−1 (> T i ) and a constant pressure P u,n−1 (> P i ). The <strong>combustion</strong> within<br />

a given zone takes p<strong>la</strong>ce progressively. After the f<strong>la</strong>me has consumed the n th zone, it is<br />

then assumed to be adiabatic. Subsequent <strong>combustion</strong> further compresses the burned<br />

gas and the unburned gas. As a result, temperature and <strong>de</strong>nsity gradients are<br />

established in the burned gas region. At the end of <strong>combustion</strong>, the burned gas cooling<br />

is computed.<br />

Figure 4.37 – Radial distribution of the multiple zones insi<strong>de</strong> a spherical vessel. Hatched portion<br />

indicates the position of the f<strong>la</strong>me front at an instant of time.<br />

Figure 4.38 shows the flowchart of the adapted Fortran co<strong>de</strong>. The output data of the<br />

co<strong>de</strong> are: the burned gas temperature, f<strong>la</strong>me radius and f<strong>la</strong>me speed, pressure, as well<br />

as the wall thermal flux. At the end of <strong>combustion</strong>, an energy ba<strong>la</strong>nce is ma<strong>de</strong>.<br />

127

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

Saved successfully!

Ooh no, something went wrong!