26.07.2013 Views

Analyse de la propagation acoustique à bassealtitude par équation ...

Analyse de la propagation acoustique à bassealtitude par équation ...

Analyse de la propagation acoustique à bassealtitude par équation ...

SHOW MORE
SHOW LESS

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

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

Saved successfully!

Ooh no, something went wrong!