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 ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
70 Méthodes numériques<br />
2.3.2 Imp<strong>la</strong>ntation numérique<br />
di Francescantonio [29] et Brentner et Farassat [10] ont repris l’équation de<br />
Ffowcs Williams & Hawkings en supprimant <strong>la</strong> condition d’imperméabilité de <strong>la</strong><br />
surface.<br />
Tij<br />
n<br />
U = 0<br />
0<br />
V s<br />
Pij<br />
Vs i<br />
propagation<br />
Fig. 2.10 – Formu<strong>la</strong>tion de l’équation de Ffowcs Williams & Hawkings par di Francescantonio,<br />
utilisant une surface arbitraire et perméable à l’écoulement.<br />
Ainsi, en utilisant une surface englobant les solides et une partie de l’air, et en<br />
remp<strong>la</strong>çant l’intérieur par de l’air au repos, on obtient <strong>pour</strong> l’équation de Ffowcs<br />
Williams & Hawkings:<br />
∂2ρ ′<br />
∂t2 − a2 ∂<br />
0<br />
2ρ ′<br />
∂x2 j<br />
(1)<br />
� �� �<br />
∂<br />
=<br />
2<br />
(2)<br />
(3)<br />
� �� � � �� �<br />
∂<br />
∂<br />
{Tij H(f)} − {Li δ(f)} + {Q δ(f)} (2.22)<br />
∂xi∂xj<br />
∂xi<br />
∂t<br />
où: Li = Pijnj + ρui(uj − VSj)nj, Q = ρ0{(1 − ρ/ρ0)VSi + ρui/ρ0}ni<br />
Dans ce cas, les termes surfaciques 2 et 3 prenent également en compte les qua-<br />
drupôles inclus dans le domaine intérieur à <strong>la</strong> surface perméable choisie. Le terme<br />
volumique 1 ne représente plus que les quadrupôles extérieurs à <strong>la</strong> surface.<br />
La surface choisie peut être en mouvement re<strong>la</strong>tif par rapport au solide. Par<br />
exemple, <strong>pour</strong> le calcul <strong>du</strong> <strong>bruit</strong> d’un rotor, on <strong>pour</strong>ra utiliser une surface fixe en-<br />
globant <strong>la</strong> totalité <strong>du</strong> rotor. D’un autre côté on vérifie que si on utilise <strong>la</strong> surface<br />
x