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.

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

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

Saved successfully!

Ooh no, something went wrong!