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 ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
58 Méthodes numériques<br />
On suppose que <strong>la</strong> turbulence est homogène isotrope, que le spectre de l’énergie<br />
turbulente peut être approché par le spectre de Kolmogorov: E(k) = Ckɛ 2/3 k −5/3<br />
où Ck ∼ 1.4, et on fait certaines hypothèse de modélisation [94] dont les principales<br />
sont:<br />
- les petites échelles ont seulement un rôle dissipatif, on néglige le phénomène de<br />
backscatter [5, 94] c’est à dire de remontée d’énergie vers les grandes échelles;<br />
- l’écoulement est en équilibre local permanent: <strong>la</strong> dissipation, effectuée en totalité<br />
par les petites échelles, est supposée égale au flux à <strong>la</strong> coupure descendant des grandes<br />
échelles.<br />
On trouve alors que νsgs, <strong>pour</strong> assurer <strong>la</strong> dissipation de l’énergie turbulente venant<br />
des grosses échelles, doit respecter <strong>la</strong> loi suivante:<br />
avec:<br />
νsgs = C 2 s ∆ 2 c| � S| (2.12)<br />
| � �<br />
S| = 2� Sij � Sij<br />
�Sij = 1<br />
�<br />
∂ �ui<br />
+<br />
2 ∂xj<br />
∂ �uj<br />
∂xi<br />
Cs = 0.18<br />
Toutefois ce modèle, conçu en turbulence homogène isotrope, in<strong>du</strong>it une sures-<br />
timation de <strong>la</strong> dissipation en présence de gradients moyens. Ces derniers, in<strong>du</strong>its<br />
au voisinage d’une paroi par exemple, interviennent de manière impromptue dans le<br />
terme | � S| et augmentent sensiblement sa valeur. Une solution consiste alors à ré<strong>du</strong>ire<br />
<strong>la</strong> valeur de Cs afin de maintenir le bon niveau de dissipation. Cette technique n’as-<br />
sure toutefois l’équilibre énergétique qu’au niveau global, en outre <strong>la</strong> valeur de Cs est<br />
alors généralement estimée de manière empirique sans construction physique faite<br />
a priori. L’idée suivante consiste à utiliser une valeur de Cs variable en espace, sui-<br />
vant par exemple une loi de Van Driest en paroi. Toutefois, <strong>la</strong> loi de Van Driest est<br />
va<strong>la</strong>ble en couche limite bidimensionnelle, mais est discutable en situation tridimen-<br />
sionnelle, ainsi qu’en présence des phénomènes complexes tels que le décollement.<br />
Elle nécessite en outre de disposer de <strong>la</strong> distance à <strong>la</strong> paroi, ce qui est ambigu en<br />
géométrie complexe. Une autre solution consiste à évaluer l’énergie à dissiper au<br />
moyen d’un filtrage supérieur, tel que proposé par Germano et al. [43] par exemple.<br />
Ceci est efficace mais coûteux en temps de calcul vis-à-vis des applications en tur-<br />
bomachines.<br />
�