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

158 Barreau-profil LES
- un schéma centré du second ordre pour les flux diffusifs;
- un avancement en temps selon une procédure de Runge-Kutta à 5 pas du second
ordre.
Par rapport au calcul RANS, le maillage a été densifié au niveau des faces du
profil (figure 5.1: domaine 2), passant de 253 à 353 points pour la circonférence, et
de 81 à 91 points dans la direction normale. En outre il a été développé dans la
troisième direction sur 3 diamètres selon 31 points. Ces évolutions, en accord avec
l’expérience acquise sur le calcul du barreau isolé, sont destinées à une représentation
des structures turbulentes ainsi que de leur dynamique qui est naturellement tridi-
mensionnelle. Comme pour le calcul sur le barreau isolé, l’envergure est limitée, en
comparaison à l’expérience qui a été menée sur une longueur de 30 diamètres. Ceci
permet de reduire les coûts de calcul. On notera que l’envergure simulée correspond
à environ une longueur de corrélation pour le barreau isolé. Le maillage est ainsi
constitué d’environ 2.2 10 6 points (∼ 72 000 × 31), et les mailles en paroi du profil
se caractérisent par les longueurs suivantes:
∆y + < 1.25 dans la direction normale à la paroi;
∆x + < 300 (< 85 pour x/c < 0.13, i.e. au voisinage du bord d’attaque) parallèlement
à la paroi;
∆z + < 350 dans la direction de l’envergure.
Le maillage était divisé en 7 domaines dans un premier temps, qui ont ensuite
été redécoupés pour obtenir 32 domaines et ainsi accroître la parallélisation. Le
découpage en domaines s’est fait en gardant la direction z intacte, afin de produire
des blocs ayant environ le même nombre de points dans chaque direction, ce qui
permet une parallélisation plus efficace par réduction des échanges entre domaines.
On rappelle que l’axe du barreau est décalé de ∆y = −0.1 d par rapport à l’ordonnée
du bord d’attaque du profil, afin de reproduire le montage expérimental. Néanmoins,
les résultats RANS ont indiqué que les effets aérodynamiques résultant sont faibles.
Aux parois des deux corps, on impose la condition d’adhérence adiabatique.
Par ailleurs, compte tenu de l’expérience acquise sur le calcul du barreau isolé, la
condition de glissement est appliquée aux limites en envergure afin de permettre le
développement tridimensionnel.
Enfin, la condition limite de non réflexion est utilisée en périphérie du domaine
de calcul pour limiter les réflexions parasites des ondes acoustiques. En outre, une
zone éponge similaire à celle utilisée pour le calcul ’barreau isolé’ est constituée par