22.10.2014 Views

simulation numerique de la combustion diphasique - cerfacs

simulation numerique de la combustion diphasique - cerfacs

simulation numerique de la combustion diphasique - cerfacs

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

– f<strong>la</strong>mme <strong>de</strong> prémé<strong>la</strong>nge : carburant et air sont mé<strong>la</strong>ngés en amont <strong>de</strong> <strong>la</strong> f<strong>la</strong>mme ;<br />

– <strong>la</strong> réaction chimique est décrite par un schéma cinétique à une seule étape irréversible<br />

(cf. section 1.5).<br />

4.3.2 Caractéristiques <strong>de</strong> <strong>la</strong> f<strong>la</strong>mme monophasique<br />

La première étape <strong>de</strong> <strong>la</strong> mise en oeuvre du calcul <strong>diphasique</strong> est l’étu<strong>de</strong> <strong>de</strong> <strong>la</strong> f<strong>la</strong>mme <strong>de</strong><br />

prémé<strong>la</strong>nge monophasique.<br />

On fait le choix d’un carburant qui est liqui<strong>de</strong> dans les conditions standards : l’éthanol<br />

(C 2 H 5 OH).<br />

La réaction chimique <strong>de</strong> <strong>combustion</strong> <strong>de</strong> l’éthanol dans l’air est <strong>la</strong> suivante :<br />

C 2 H 5 OH + 3 (O 2 + 3.76 N 2 ) −→ 2 CO 2 + 3 H 2 O + 11.28 N 2 (4.4)<br />

L’écoulement étudié est schématisé sur <strong>la</strong> figure 4.2.<br />

U 1<br />

S L<br />

gaz frais<br />

gaz brules<br />

x<br />

entree<br />

front <strong>de</strong> f<strong>la</strong>mme<br />

sortie<br />

Fig. 4.2 – Schéma <strong>de</strong>scriptif <strong>de</strong> <strong>la</strong> f<strong>la</strong>mme <strong>de</strong> prémé<strong>la</strong>nge.<br />

Dans une f<strong>la</strong>mme <strong>la</strong>minaire <strong>de</strong> prémé<strong>la</strong>nge, <strong>la</strong> vitesse <strong>de</strong> propagation du front <strong>de</strong> f<strong>la</strong>mme (vitesse<br />

<strong>de</strong> f<strong>la</strong>mme <strong>la</strong>minaire) est principalement fonction <strong>de</strong> <strong>la</strong> cinétique <strong>de</strong> réaction et <strong>de</strong>s conditions<br />

<strong>de</strong> température et <strong>de</strong> pression <strong>de</strong>s gaz frais.<br />

C’est une caractéristique majeure <strong>de</strong> <strong>la</strong> f<strong>la</strong>mme, que l’on définit par :<br />

∫ +∞<br />

−∞<br />

s L =<br />

˙ω F dx<br />

ρ 1 (YF 1 − Y F 2) (4.5)<br />

où l’indice 1 correspond aux gaz frais et l’indice 2 aux gaz brûlés.<br />

Le schéma cinétique réduit est donc déterminé <strong>de</strong> façon à reproduire <strong>la</strong> vitesse <strong>de</strong> f<strong>la</strong>mme <strong>la</strong>minaire<br />

mesurée expérimentalement.<br />

La procédure d’ajustement du schéma opérée ici est <strong>la</strong> suivante :<br />

1. On dispose <strong>de</strong> données expérimentales (Hayashi et Kumagai [6]).<br />

2. On utilise un outil <strong>de</strong> référence pour le calcul <strong>de</strong> f<strong>la</strong>mmes p<strong>la</strong>nes <strong>la</strong>minaires prémé<strong>la</strong>ngées :<br />

PREMIX. Il s’agit d’un co<strong>de</strong> <strong>de</strong> calcul développé par SANDIA National Laboratories<br />

résolvant les écoulements gazeux réactifs stationnaires et incompressibles en une dimension<br />

avec <strong>de</strong>s propriétés thermophysiques détaillées.<br />

Les coefficients n F , n O , β, A et E a <strong>de</strong> l’équation (1.32) sont ajustés <strong>de</strong> manière à faire<br />

coïnci<strong>de</strong>r les valeurs <strong>de</strong> s L calculées par PREMIX et les données exprérimentales pour une<br />

gamme <strong>de</strong> richesse al<strong>la</strong>nt <strong>de</strong> 0.6 à 1.1.<br />

45

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

Saved successfully!

Ooh no, something went wrong!