these approches numeriques pour la simulation du bruit a large ...

130 Ecoulements avec parois en LES
résolution des plus grosses structures de paroi dont l’espacement caractéristique est
de l’ordre de ∆z + = 100.
Ce maillage a ensuite été découpé selon y en 4 domaines de dimensions topologiques
similaires afin de mener des calculs parallèles sur 4 processeurs et réduire ainsi les
temps d’obtention des résultats (voir figure 4.1).
Les conditions aux limites sont périodiques selon x et z, et une condition d’adhérence
adiabatique est imposée aux parois (y = ±δ).
1er calcul: amortissement selon Van Driest en paroi
Cette modification du modèle de Smagorinsky consiste seulement à appliquer,
au terme de dissipation de sous-maille, une fonction d’amortissement suivant une loi
de Van Driest. Ainsi on compense près de la paroi la sur-estimation de la viscosité
turbulente engendrée par l’influence du gradient moyen de vitesse dans les Sij.
avec:
et:
νsgs = � fd(y + ) � 2 C 2 s ∆ 2 c | � S|
| � �
S| = 2� Sij � Sij
�Sij = 1
�
∂ �ui
+
2 ∂xj
∂ �uj
∂xi
Cs = 0.18
fd(y + ) = 1 − e −y+ /25
Ce modèle est uniquement appliqué au calcul de canal plan, où il sert de référence.
En effet, il est conçu pour un écoulement proche d’une couche limite bidimension-
nelle et se trouve inadapté dans le cas de décollements, d’effets tridimensionnels ou
de parois non planes. En outre, il est difficilement utilisable en géométrie complexe
car sa formulation nécessite l’estimation de la distance à la paroi.
Ce calcul a été initialisé en interpolant sur notre maillage un résultat instan-
tané de calcul effectué sur la même configuration mais avec un maillage légèrement
différent (calcul LES par E. Sergent [98]). Il aurait été possible d’initialiser le calcul
par un profil laminaire et d’attendre la transition puis l’établissement de la turbu-
lence, mais cette procédure aurait été coûteuse en temps et inutile, le calcul de la
transition n’étant pas au centre de notre étude. Ici le champ initial dispose déja de
�