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.
Ecoulements avec parois en LES 139<br />
Uoo<br />
d<br />
Φ<br />
y<br />
d<br />
ο<br />
x<br />
Ac 1<br />
y =<br />
185d<br />
x = 15d<br />
Fig. 4.10 – Configuration et notations <strong>pour</strong> le cas ’barreau isolé’ [ : surface 2<br />
(paroi <strong>du</strong> barreau), · : surface 3 (à 1 d de <strong>la</strong> paroi) ]. Remarque: on définit r<br />
comme le rayon associé aux coordonnées cylindriques.<br />
4.2.1 Configuration numérique<br />
Le code de calcul utilise:<br />
- un schéma centré <strong>pour</strong> les flux convectifs (quatrième ordre en monodimensionnel),<br />
sans viscosité numérique;<br />
- un schéma centré <strong>du</strong> second ordre <strong>pour</strong> les flux diffusifs;<br />
- un avancement en temps selon une procé<strong>du</strong>re de Runge-Kutta à 5 pas <strong>du</strong> second<br />
ordre.<br />
Le mail<strong>la</strong>ge s’étend sur un rayon de 20 diamètres autour de l’axe <strong>du</strong> barreau,<br />
et sur une envergure de 3 diamètres. Cette dernière longueur est limitée afin de<br />
ré<strong>du</strong>ire les coûts de calcul, elle permet toutefois <strong>la</strong> représentation d’une longueur<br />
de corré<strong>la</strong>tion (Lc ∼ 3 d selon Szepessy [107] <strong>pour</strong> Red = 4.3 10 4 ). La figure 4.11<br />
présente <strong>la</strong> vue d’une section <strong>du</strong> mail<strong>la</strong>ge à envergure fixée. Cette dernière, constituée<br />
de 197 points selon <strong>la</strong> circonférence et 193 selon le rayon, est répétée à intervalles<br />
réguliers le long de l’envergure afin de <strong>la</strong> discrétiser en 31 points. Le nombre total<br />
de points est d’environ 1.2 millions. Ainsi, les mailles en paroi se caractérisent par<br />
les longueurs suivantes:<br />
∆y + < 3 (< 2 dans <strong>la</strong> zone décollée turbulente) dans <strong>la</strong> direction normale à <strong>la</strong> paroi;<br />
∆x + < 200 (< 90 dans <strong>la</strong> zone décollée turbulente) parallèlement à <strong>la</strong> paroi;<br />
∆z + < 550 (< 350 dans <strong>la</strong> zone décollée turbulente) dans <strong>la</strong> direction de l’envergure.<br />
Ces valeurs sont sensiblement plus élevées que celles utilisées dans le canal p<strong>la</strong>n, ceci