Analyse de la propagation acoustique à bassealtitude par équation ...
Analyse de la propagation acoustique à bassealtitude par équation ...
Analyse de la propagation acoustique à bassealtitude par équation ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6.1.1 - Introduction <strong>de</strong>s conditions aux limites dans l'algorithme "Split-Step Fourier"<br />
La métho<strong>de</strong> "Split-Step" est bien adaptée <strong>à</strong> <strong>la</strong> <strong>propagation</strong> sous-marine où les profils <strong>de</strong> célérité<br />
du son créent un chenal <strong>acoustique</strong>. L'énergie <strong>acoustique</strong> est canalisée <strong>à</strong> faible profon<strong>de</strong>ur vis-<strong>à</strong>-vis <strong>de</strong>s<br />
distances <strong>de</strong> <strong>propagation</strong>. Ce phénomène est dû au fait que <strong>la</strong> célérité du son croît régulièrement <strong>à</strong> <strong>par</strong>tir<br />
d'une certaine profon<strong>de</strong>ur.<br />
La métho<strong>de</strong> "Split-Step" impose que <strong>la</strong> fonction qi et ses dérivées successives <strong>par</strong> rapport <strong>à</strong> z<br />
ten<strong>de</strong>nt vers zéro avec <strong>la</strong> profon<strong>de</strong>ur. En effet, le domaine d'intégration étant limité en profon<strong>de</strong>ur, il est<br />
impératif <strong>de</strong> ne pas créer d'effets <strong>de</strong> repliement lors <strong>de</strong> l'application <strong>de</strong> transformées <strong>de</strong> Fourier discrètes.<br />
Par ailleurs, dans le cadre <strong>de</strong> <strong>la</strong> <strong>propagation</strong> sous-marine, <strong>la</strong> condition <strong>de</strong> réflexion totale <strong>à</strong> <strong>la</strong><br />
surface est prise en compte en ajoutant symétriquement <strong>à</strong> <strong>la</strong> surface un <strong>de</strong>mi-espace fictif. Le profil <strong>de</strong><br />
célérité est choisi symétrique et on impose l'im<strong>par</strong>ité <strong>de</strong> <strong>la</strong> fonction qi ou z. On crée ainsi une source<br />
virtuelle image <strong>de</strong> <strong>la</strong> source réelle et déphasée <strong>de</strong> n.<br />
6.1.2 - Pas d'échantillonnage vertical et pas d'intégration<br />
Les résultats que nous présentons ici ont été obtenus <strong>par</strong> Cristol [271.<br />
Le théorème <strong>de</strong> Shannon impose un pas d'intégration vertical:<br />
2.h rn<br />
= N-1<br />
si hrn est <strong>la</strong> profon<strong>de</strong>ur du domaine et N le nombre <strong>de</strong> points <strong>de</strong> l'échantillonnage vertical. Ceci afin <strong>de</strong><br />
calculer correctement <strong>la</strong> transformée <strong>de</strong> Fourier:<br />
(r,$)F1w(r,z)}<br />
Pour le calcul <strong>de</strong> <strong>la</strong> transformée <strong>de</strong> Fourier inverse, le théorème <strong>de</strong> Shannon impose:<br />
m<br />
Pour le pas r, il n'existe pas <strong>de</strong> critère général, <strong>la</strong> pratique montre que pour les profils d'indice<br />
les plus fréquemment observés, un pas <strong>de</strong> l'ordre <strong>de</strong> plusieurs centaines <strong>de</strong> mètres suffit, pour le domaine<br />
<strong>de</strong>s basses fréquences (f 200 Hz).<br />
IN \2 2 F (N<br />
=<br />
\2<br />
I - - i<br />
1N2h<br />
sur le domae - - 1<br />
t 2<br />
\2 (N \2<br />
iNh -<br />
rn<br />
De même, <strong>la</strong> pratique montre que le théorème <strong>de</strong> Shannon est légèrement insuffisant ; le pas<br />
d'échantillonnage vertical doit être <strong>de</strong> quelques fois inférieur <strong>à</strong> <strong>la</strong> longueur d'on<strong>de</strong>.<br />
S, - 1 ¡<br />
2 iNhi<br />
rn