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

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

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

36 Chapitre 2 Modèles physiques<br />

Trois versions <strong>du</strong> modèle ont été successivement développées.<br />

– Le modèle original <strong>de</strong> Durbin [Durbin 95] comprend une condition limite <strong>de</strong> paroi<br />

pour f, le terme <strong>de</strong> contrainte <strong>de</strong> pression obtenu à partir d’une équation <strong>de</strong> relaxa-<br />

tion elliptique. Il nécessite aussi la distance aux parois pour déterminer le coefficient<br />

Cε1 qui régule la pro<strong>du</strong>ction <strong>de</strong> dissipation <strong>de</strong> l’énergie cinétique turbulente dans<br />

l’équation sur ε.<br />

– La version développée par Parneix & Durbin [S. Parneix 97] intro<strong>du</strong>it une nouvelle<br />

expression pour le coefficient Cε1 permutant la distance à la paroi <strong>et</strong> le rapport v2/k.<br />

– Finalement, Lien & Durbin donnent une nouvelle condition <strong>de</strong> paroi pour f <strong>de</strong> façon<br />

à augmenter la robustesse <strong>numérique</strong> <strong>du</strong> modèle [Lien 96].<br />

Une comparaison systématique <strong>de</strong> ces trois versions a été réalisée [Kalitzin 99] sur<br />

un écoulement autour d’une aile d’avion <strong>et</strong> sur un obstacle. L’auteur montre qu’elles<br />

prédisent <strong>de</strong>s résultats très similaires dans la région potentielle <strong>de</strong> l’écoulement <strong>et</strong> dans<br />

la majeure partie <strong>de</strong> la couche limite. Cependant, la définition <strong>de</strong> Cε1 v 2 /k provoque<br />

une surestimation <strong>du</strong> taux <strong>de</strong> croissance <strong>de</strong> la couche limite pour les écoulements à haut<br />

Reynolds avec forts gradients <strong>de</strong> pression. Imposer f nul aux parois r<strong>et</strong>ar<strong>de</strong> la transition,<br />

ce qui est souvent intéressant dans les prédictions d’écoulements <strong>de</strong> transition.<br />

Le modèle v 2 − f est composé <strong>de</strong> trois équations <strong>de</strong> transport, une sur k (l’énergie<br />

cinétique turbulente), une sur ε (la dissipation <strong>de</strong> l’énergie cinétique turbulente), <strong>et</strong> une<br />

sur v 2 .<br />

Le terme <strong>de</strong> pro<strong>du</strong>ction dans l’équation sur v 2 est un terme <strong>de</strong> contrainte <strong>de</strong> pression.<br />

Sa dépendance non locale dans l’écoulement, plus particulièrement en présence <strong>de</strong> parois<br />

soli<strong>de</strong>s, est modélisée avec une équation <strong>de</strong> relaxation elliptique pour la quantité f.<br />

La viscosité turbulente cinématique s’exprime comme le carré d’une vitesse v <strong>et</strong> d’un<br />

temps Tt caractéristiques <strong>de</strong> l’écoulement :<br />

µt = ρCµv 2 Tt<br />

(2.79)<br />

On tient ainsi compte <strong>de</strong> l’anisotropie turbulente près <strong>de</strong>s parois à la différence <strong>du</strong><br />

modèle k − ε.<br />

Le temps caractéristique Tt <strong>de</strong> la turbulence est défini par :<br />

<br />

k ν<br />

Tt = max ; CTt(<br />

ε ε )1/2<br />

<br />

(2.80)

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

Saved successfully!

Ooh no, something went wrong!