22.10.2014 Views

Rapport de stage de DEA - référence CERFACS : WN/CFD/03/75 ...

Rapport de stage de DEA - référence CERFACS : WN/CFD/03/75 ...

Rapport de stage de DEA - référence CERFACS : WN/CFD/03/75 ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3.2 Analyse <strong>de</strong>s champs moyens<br />

3.2.1 Champ moyen <strong>de</strong> richesse en sortie d’injecteur<br />

A partir <strong>de</strong>s concentrations massiques d’air et <strong>de</strong> méthane, il est aisé <strong>de</strong> reconstruire la<br />

richesse du mélange. Comme l’on s’intéresse ici à la qualité du mélange en sortie d’injecteur,<br />

il est naturel <strong>de</strong> commencer par regar<strong>de</strong>r la variable richesse φ moyennée dans le plan <strong>de</strong> sortie<br />

<strong>de</strong> l’injecteur. Cette variable est définie <strong>de</strong> la manière suivante :<br />

φ = s Y F<br />

Y O<br />

où s = ν Ow O<br />

ν F w F<br />

= 4 (3.1)<br />

avec Y k les fractions massiques, ν k les coefficients stoechiométriques et w k les masses molaires<br />

<strong>de</strong>s espèce indicées O pour l’oxygène et F pour le méthane. Cependant, il est important <strong>de</strong> se<br />

rappeler que le maillage utilisé est périodique sur un angle <strong>de</strong> 30 <strong>de</strong>grés. Ce détail fait <strong>de</strong> l’axe<br />

central un axe qu’aucune vitesse ne peut couper. Autrement dit, on interdit aux particules <strong>de</strong><br />

flui<strong>de</strong> <strong>de</strong> traverser l’axe central <strong>de</strong> la géométrie. La conséquence <strong>de</strong> cet écart à la réalité est une<br />

surestimation <strong>de</strong> la taille <strong>de</strong> la zone <strong>de</strong> recirculation à la pointe du cône <strong>de</strong> l’injecteur. De plus,<br />

le maillage ne traduit pas l’élargissement <strong>de</strong> la géométrie réelle lors du passage <strong>de</strong> l’injecteur à<br />

la chambre <strong>de</strong> combustion, ce qui change complètement la dynamique du "vortex breakdown"<br />

(dégénérescence <strong>de</strong> tourbillon, voir [?]) dans la chambre. Enfin, il n’y a pas <strong>de</strong> réaction dans le<br />

calcul contrairement au cas réel où la zone <strong>de</strong> réaction se trouve dans la chambre <strong>de</strong> combustion<br />

mais également à quelques millimètres en amont <strong>de</strong> la pointe du cône <strong>de</strong> l’injecteur. Notons<br />

tout <strong>de</strong> même que notre calcul n’a jamais eu la prétention <strong>de</strong> simuler le "vortex breakdown" ni<br />

la chimie <strong>de</strong> la réaction <strong>de</strong> combustion. Quoiqu’il en soit, pour toutes ces raisons, notre calcul<br />

est certainement erroné à l’extrémité du cône <strong>de</strong> l’injecteur, ce qui nous engage à nous focaliser<br />

sur le plan P exhibé par la figure ?? lors <strong>de</strong> la visualisation <strong>de</strong> nos résultats en sortie d’injecteur.<br />

Ce plan orthogonal à l’axe central se situe à 5,5mm en amont <strong>de</strong> la pointe <strong>de</strong> l’injecteur.<br />

Portons ainsi notre attention sur le champ moyen <strong>de</strong> richesse dans le plan P en nous rapportant<br />

à la figure ??. Le tracé <strong>de</strong> la richesse en fonction du rayon pour différents rayons montre<br />

déjà que l’on peut considérer la répartition <strong>de</strong> la richesse comme quasiment axisymétrique en<br />

moyenne. Cependant, un gradient <strong>de</strong> richesse apparait dans la direction radiale puisque l’on<br />

passe <strong>de</strong> phi=0,5 contre le cône à phi=0,8 à la périphérie <strong>de</strong> la sortie <strong>de</strong> l’injecteur avec une<br />

moyenne globale pour le dispositif <strong>de</strong> phi=0,70. Force est donc <strong>de</strong> constater que cela ne ressemble<br />

pas à un profil <strong>de</strong> prémélange parfait, même en moyenne.<br />

3.2.2 Champs moyens <strong>de</strong> vitesse dans le plan P<br />

Une autre variable à regar<strong>de</strong>r dans le plan <strong>de</strong> sortie <strong>de</strong> l’injecteur est la vitesse. La figure<br />

?? présente <strong>de</strong> gauche à droite les composantes axiales, radiales et orthoradiales du champ <strong>de</strong><br />

vitesse moyen dans le plan P. On retrouve ici aussi une propriété d’axisymétrie du champ <strong>de</strong><br />

vitesse comme on l’a déjà remarqué pour le champ moyen <strong>de</strong> richesse. Mais là encore, le champ<br />

<strong>de</strong> vitesse moyen n’a rien d’uniforme et présente sans équivoque <strong>de</strong>s gradients. La condition aux<br />

limites <strong>de</strong> murs adhérents est également bien visible.<br />

16

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

Saved successfully!

Ooh no, something went wrong!