27.06.2013 Views

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

Simulation numérique de l'essorage et du refroidissement d'un film ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

3.3 Validations 131<br />

En sortie <strong>de</strong> buse, on constate la présence d’un secteur angulaire d’intensité Iw<br />

très faible. Il correspond au noyau potentiel. L’épaississement <strong>de</strong>s <strong>de</strong>ux couches <strong>de</strong><br />

mélange <strong>de</strong> part <strong>et</strong> d’autre <strong>du</strong> j<strong>et</strong> sont le lieu <strong>de</strong> l’intensité maximale <strong>et</strong> tra<strong>du</strong>isent<br />

son épanouissement. On note aussi la présence <strong>de</strong> <strong>de</strong>ux grosses recirculations.<br />

On observe le maximum <strong>de</strong> la composante Iu <strong>de</strong> l’intensité turbulente au niveau <strong>du</strong><br />

point <strong>de</strong> stagnation <strong>du</strong> j<strong>et</strong>. Sur les côtés apparaissent les <strong>de</strong>ux grosses recirculations.<br />

L’intensité turbulente Iu augmente dans la zone <strong>de</strong> développement <strong>du</strong> j<strong>et</strong> puis di-<br />

minue à mesure qu’on approche <strong>de</strong> la paroi, pour atteindre sa valeur maximale au<br />

point <strong>de</strong> stagnation <strong>du</strong> j<strong>et</strong>, autour <strong>de</strong> 0,165. C<strong>et</strong>te information est corroborée par<br />

les simulations LES réalisées par Beaubert [Beaubert 02] pour différents Reynolds,<br />

<strong>et</strong> son évolution longitudinale est tracée à la figure 3.39 <strong>et</strong> comparée à <strong>de</strong>s valeurs<br />

expérimentales.<br />

Iu<br />

0.3<br />

0.25<br />

0.2<br />

0.15<br />

0.1<br />

0.05<br />

0<br />

0 2 4 6 8 10<br />

z/d<br />

Figure 3.39 Répartition <strong>de</strong> la composante selon x <strong>de</strong> l’intensité turbulente Iu, : données<br />

expérimentales [Maurel 01b], trait : simulation<br />

L’évolution globale est respectée même si la chute d’intensité turbulente autour <strong>de</strong><br />

z/d = 2 est surestimée par le simulation. En outre, la valeur <strong>du</strong> pic d’intensité à<br />

la paroi n’est pas correctement ren<strong>du</strong>, les simulations LES menées par Beaubert<br />

[Beaubert 02] soulignent ce problème, invoquant une lacune <strong>de</strong> maillage dans c<strong>et</strong>te<br />

zone pour y capter les gradients élevés.<br />

Concernant la composante Iw, on note là encore la présence <strong>du</strong> noyau potentiel. Le<br />

caractère fortement tridimensionnel est mis en exergue notamment au niveau <strong>de</strong>s<br />

recirculations ainsi qu’à la paroi.

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

Saved successfully!

Ooh no, something went wrong!