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.
100 Barreau-profil RANS<br />
La concentration en paroi des régions de sous- et sur-vitesse, par les simu<strong>la</strong>tions, peut<br />
être reliée au processus de moyenne d’ensemble utilisé par les <strong>approches</strong> RANS ins-<br />
tationnaires. Les structures turbulentes ne sont pas représentées, seuls ressortent les<br />
tourbillons de von Kármán. Ainsi, des calculs de jets ont montré une sous-estimation<br />
de <strong>la</strong> diffusion turbulente. Par ailleurs <strong>la</strong> convection des tourbillons de von Kármán<br />
se fait selon <strong>la</strong> trajectoire moyenne: il n’y a pas de dispersion des trajectoires par<br />
<strong>la</strong> turbulence contrairement aux visualisation expérimentales de Cantwell et Coles<br />
[13], ce qui ré<strong>du</strong>it <strong>la</strong> zone d’influence des tourbillons.<br />
Enfin, <strong>la</strong> section [C] fait apparaître des résultats à nouveau consistants <strong>pour</strong> les cal-<br />
culs et représentatifs <strong>du</strong> sil<strong>la</strong>ge mesuré. On retrouve un é<strong>la</strong>rgissement supérieur <strong>du</strong><br />
sil<strong>la</strong>ge dans les mesures, les tourbillons étant encore présents après le profil dans les<br />
mesures [80] et les calculs (figures 3.3). On observe toutefois une dissymétrie dans<br />
les résultats expérimentaux qui ne se retrouve pas dans les calculs. Le mail<strong>la</strong>ge <strong>du</strong><br />
calcul prend en compte le déca<strong>la</strong>ge de ∆y = 1 mm de l’axe <strong>du</strong> barreau par rap-<br />
port au bord d’attaque <strong>du</strong> profil, mais ce<strong>la</strong> n’a que peu d’effet sur <strong>la</strong> symétrie de<br />
l’écoulement calculé. La dissymétrie observée dans les mesures semble plutôt reliée<br />
à une incidence parasite <strong>du</strong> profil expérimental (∼ 2 deg) d’après l’étude faite sur le<br />
Cp.<br />
Vitesse fluctuante<br />
L’anémométrie fil chaud à un fil permet de connaître <strong>la</strong> valeur efficace u ′ de <strong>la</strong><br />
vitesse fluctuante dans <strong>la</strong> direction de <strong>la</strong> vitesse moyenne. On se basera donc sur<br />
cette grandeur <strong>pour</strong> l’étude des calculs.<br />
Il faut noter que dans l’approche RANS une grande partie de <strong>la</strong> vitesse fluctuante<br />
n’est pas résolue mais se trouve modélisée au travers <strong>du</strong> tenseur de Reynolds. On<br />
reconstitue donc ici u ′ de manière à prendre en compte à <strong>la</strong> fois <strong>la</strong> valeur fluc-<br />
tuante résolue, mais également l’énergie cinétique turbulente estimée au travers des<br />
équations de transport:<br />
u ′ =<br />
�<br />
< (�u.�n − | < �u >t |) 2 >t + 2<br />
< k > (3.1)<br />
3<br />
où �u représente <strong>la</strong> vitesse résolue et �n =< �u >t /| < �u >t | <strong>la</strong> direction de <strong>la</strong> vitesse<br />
moyenne.<br />
Le terme < (�u.�n − | < �u >t |) 2 >t désigne donc <strong>la</strong> fluctuation de vitesse résolue<br />
dans <strong>la</strong> direction de l’écoulement moyen, tandis que 2<br />
3 < k >t représente <strong>la</strong> contribution<br />
de <strong>la</strong> turbulence, supposée isotrope. Cette reconstruction n’est donc qu’ap-<br />
proximative.