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.
82 Méthodes numériques<br />
est supposée indépendante de z (i.e. l’espérance mathématique est supposée égale à<br />
<strong>la</strong> moyenne selon l’envergure).<br />
Alors, si on note Spp(f) <strong>la</strong> densité spectrale de puissance en champ lointain, et si<br />
on marque par exp et sim respectivement les conditions expérimentales et simulées,<br />
on peut appliquer <strong>la</strong> correction suivante aux résultats de <strong>la</strong> simu<strong>la</strong>tion:<br />
(Spp(f)) exp = (Spp(f)) sim + 10.0 log<br />
� � Lexp<br />
� Lexp<br />
0 0<br />
� � Lsim Lsim<br />
0 0<br />
Γ(∆z,f)dz1dz2<br />
Γ(∆z,f)dz1dz2<br />
�<br />
(2.27)<br />
Il suffira donc d’estimer Γ(∆z,f) sur <strong>la</strong> base des mesures ou <strong>du</strong> calcul. Nous<br />
utiliserons le calcul, afin de ne pas intro<strong>du</strong>ire d’information venant de l’expérience,<br />
puisque <strong>la</strong> LES est supposée prendre en compte l’aspect aléatoire. Toutefois, dans<br />
les simu<strong>la</strong>tions présentées ici, Γ n’est connue que sur envergure Lsim = 3 d, alors qu’il<br />
faudrait <strong>la</strong> connaître sur toute l’envergure expérimentale. Ce problème sera résolu en<br />
approchant ∆z → Γ(∆z,f) par une courbe exponentielle: ∆z → exp(−|∆z|/Le(f))<br />
sur l’envergure simulée. La connaissance de Le(f) ainsi obtenue <strong>pour</strong> chaque fréquence<br />
f permettra d’extrapoler <strong>la</strong> Γ sur l’envergure voulue.<br />
Remarque: <strong>pour</strong> un fréquence f donnée, si on considère Γ(∆z,f) sous <strong>la</strong> forme<br />
d’une fonction échelon (Γ(∆z,f) = 1 <strong>pour</strong> ∆z ≤ Lc, et Γ(∆z,f) = 0 <strong>pour</strong> ∆z > Lc),<br />
on retrouve <strong>la</strong> correction proposée par Kato et al. [58, 57]:<br />
Si Lc < Lsim : Sppexp(f) = Sppsim(f) + 10 log(Lexp/Lsim)<br />
Si Lsim < Lc < Lexp : Sppexp(f) = Sppsim(f) + 20 log(Lc/Lsim) + 10 log(Lexp/Lc)<br />
Si Lexp < Lc : Sppexp(f) = Sppsim(f) + 20 log(Lexp/Lsim)<br />
2.4 Intégration <strong>du</strong> processus de calcul<br />
Pour conclure ce chapitre de présentation des méthodes numériques, <strong>la</strong> figure 2.11<br />
propose un schéma rappe<strong>la</strong>nt <strong>la</strong> combinaison des procé<strong>du</strong>res de calcul aérodynamiques<br />
et acoustiques.<br />
Comme il a déja été suggéré, <strong>la</strong> première phase <strong>du</strong> calcul acoustique, notée Ad-<br />
vantia (1) et correspondant au calcul des intégrandes, peut s’effectuer en parallèle<br />
<strong>du</strong> calcul aérodynamique. La phase Advantia (2) regroupe ensuite les contributions<br />
acoustiques calculées afin de créer le signal intégral. Toutefois, dans le cadre de<br />
cette étude <strong>la</strong> parallélisation n’a été que simulée: les champs aérodynamiques uti-<br />
lisés <strong>pour</strong> le calcul des intégrandes ont été stockés jusqu’à <strong>la</strong> fin <strong>du</strong> calcul. Ensuite