these approches numeriques pour la simulation du bruit a large ...
these approches numeriques pour la simulation du bruit a large ...
these approches numeriques pour la simulation du bruit a large ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
46 Méthodes numériques<br />
et � Pk = τtij<br />
∂ eui<br />
∂xj<br />
Modèles de turbulence<br />
désigne <strong>la</strong> pro<strong>du</strong>ction d’énergie turbulente.<br />
Le modèle de turbulence choisi détermine <strong>la</strong> construction de τtij ainsi que le calcul<br />
des coefficients intervenant dans les équations de transport (k et ω) et dans les<br />
termes modélisés. Trois modèles ont été mis en œuvre dans le cadre de cette étude<br />
et sont présentés ici par ordre croissant de complexité.<br />
• Modèle linéaire<br />
Il s’agit <strong>du</strong> modèle k − ω de base, développé par D.C. Wilcox [115]. L’utilisation<br />
de <strong>la</strong> variable ω comme seconde variable turbulente est moins répen<strong>du</strong>e que ɛ (<strong>la</strong><br />
dissipation turbulente) car son interprétation physique est plus complexe. Néanmoins<br />
cette formu<strong>la</strong>tion permet un traitement plus aisé des parois: <strong>la</strong> valeur de ω est<br />
connue sur <strong>la</strong> surface (→ ∞) et il n’apparaît pas nécessaire d’utiliser de loi empirique<br />
en proche paroi [114]. Ce comportement autorise alors généralement une densité<br />
<strong>du</strong> mail<strong>la</strong>ge de l’ordre de ∆y + = 5 perpendicu<strong>la</strong>irement à <strong>la</strong> surface. En outre <strong>la</strong><br />
variable ω permet de limiter les erreurs par rapport à ɛ en écoulement incompressible<br />
par une meilleure estimation des termes dominants [103]. De manière générale, les<br />
comparaisons avec les modèles k − ɛ con<strong>du</strong>isent à l’équivalence des modèles, voire à<br />
un meilleur comportement des modèles k − ω en écoulement décollé [79].<br />
Le modèle se définit de <strong>la</strong> manière suivante:<br />
avec:<br />
et les coefficients:<br />
τtij<br />
�<br />
∂ �ui<br />
= µt +<br />
∂xj<br />
∂ �uj<br />
−<br />
∂xi<br />
2<br />
3 δij<br />
�<br />
∂ �uk<br />
∂xk<br />
µt = Cµρ � k<br />
�ω<br />
− 2<br />
3 ρ� kδij<br />
Cµ = 1 σk = 2.0 σω = 2.0 Ck = 0.09 Cω1 = 5<br />
9<br />
Cω2 = 3<br />
40<br />
calibrés en fonction d’écoulements de référence (turbulence homogène isotrope,<br />
couche limite turbulente sur p<strong>la</strong>que p<strong>la</strong>ne...)[115].