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.
Méthodes numériques 71<br />
solide (ui.ni = VSi.ni sur <strong>la</strong> surface), on retrouve bien l’expression (2.21) obtenue<br />
par Ffowcs Williams & Hawkings.<br />
Cette formu<strong>la</strong>tion permet d’appliquer <strong>la</strong> physique issue de l’équation de Ffowcs<br />
Williams & Hawkings (prise en compte explicite des sources, pas de calcul de gra-<br />
dient normal à <strong>la</strong> surface) à une surface d’intégration indépendante de <strong>la</strong> surface<br />
solide. Elle permet en outre de prendre en compte d’éventuelles sources volumiques<br />
par les termes calculés sur une surface englobant ces sources. Ceci tend donc à se<br />
rapprocher de l’équation de Kirchhoff, qui propage en champ lointain une onde es-<br />
timée depuis une surface initiale quelquonque, <strong>la</strong>quelle doit se situer en dehors de <strong>la</strong><br />
région perturbée de l’écoulement.<br />
L’équation 2.22 peut ensuite être mise sous forme intégrale par convolution<br />
avec une fonction de Green. On note: �x et t <strong>la</strong> position et le temps de l’auditeur<br />
(réception), �y et τ <strong>la</strong> position et le temps d’intégration des sources (émission). Le<br />
vecteur radial est alors: �r = �x − �y, et les composantes <strong>du</strong> vecteur radial unitaire:<br />
ˆri = ri/r = (xi − yi)/r. On peut ainsi obtenir une expression de ρ ′ (�x,t) avec <strong>la</strong><br />
fonction de Green va<strong>la</strong>ble en champ libre G = δ(g)/r où g = t − τ − r/a0, et le<br />
changement de variable:<br />
�<br />
Q(τ)δ(g(τ)) =<br />
<strong>la</strong> somme étant prise sur tous les zéros τ n ret<br />
N�<br />
n=1<br />
Q<br />
|∂g/∂τ| (τ n ret )<br />
de l’équation au temps retardé: g = 0. Si<br />
on considère que Mr = Miˆri < 1 au niveau des sources, comme ce sera toujours le<br />
cas dans <strong>la</strong> présente étude, l’équation au temps retardé ne possède qu’une solution<br />
τ = t − ||�x − �y(τ)||/a0 et on a:<br />
4πa 2 0ρ ′ = ∂2<br />
∂xi∂xj<br />
�<br />
f>0<br />
�<br />
Tij<br />
r(1 − Mr)<br />
�<br />
+ ∂<br />
� � �<br />
Q<br />
∂t f=0 r(1 − Mr)<br />
ret<br />
dS<br />
ret<br />
dV − ∂<br />
∂xi<br />
�<br />
f=0<br />
�<br />
Li<br />
r(1 − Mr)<br />
�<br />
où le symbole [.]ret signifie que les grandeurs sont évaluées à: τ = t − r/a0.<br />
En champ lointain on peut ensuite écrire: a 2 0 ρ′ = p ′ .<br />
ret<br />
dS<br />
(2.23)<br />
La solution est ici obtenue <strong>pour</strong> un milieu extérieur libre au repos. Dans le cas plus<br />
général (écoulement uniforme, écoulement extérieur non homogène, présence de pa-<br />
rois...), une solution de Green adaptée doit être utilisée à moins de justifier que les