25.02.2013 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!