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

140 Ecoulements avec parois en LES
afin de limiter le temps de calcul. Toutefois, il faut noter que dans le cas du barreau
présenté ici, la transition ne survient qu’après le décollement et l’écoulement est
alors dominé par le lâcher tourbillonnaire, dont la longueur caractéristique est égale
au diamètre du barreau.
y
z
x
Fig. 4.11 – Vue du maillage au niveau de la section (z = 0).
L’adhérence adiabatique du fluide est imposée sur la paroi du barreau.
Afin limiter le confinement des ondes acoustiques à l’intérieur du domaine de calcul,
la condition de non réflexion est appliquée à la frontière circonférentielle. En outre
une zone éponge est constituée en étirant les mailles en périphérie (figure 4.11) et
en ajoutant une viscosité numérique dans le schéma en espace 2 selon la loi:
ɛze = 0.1 + 0.22
π
× arctan (50 × (r − 0.13)) (4.2)
où r est le rayon au point considéré dans le domaine de calcul (en m). La figure
4.12 présente un tracé de ɛze(r). On remarque que: ɛze(0.5d) = 10 −3 , la visco-
sité numérique est donc très faible au niveau de la paroi du barreau. Ceci permet
d’espérer que l’influence sur les résultats LES reste faible, ce qui sera vérifié a poste-
riori en compararant des grandeurs aérodynamiques calculées avec l’expérience (en
particulier les fluctuations, mais aussi la position du décollement qui est sensible à
Red et donc à la viscosité). Par contre, dans la région périphérique du domaine de cal-
cul (15 ≥ r/d ≤ 20), ɛze devient relativement forte et atteint le niveau généralement
utilisé pour stabiliser des calculs Euleriens à M∞ ∼ 0.2 [103, 104].
2. cf. chapitre 2, paragraphe sur les schémas pour les flux convectifs.