B.6.

Problème complémentaire

B.6. Problème complémentaire

Remarque préliminaire : Ce problème sur les fonctions de Green prolonge le TD sur la relation
entre distributions et fonctions de Green. L’important est le lien que l’on peut établir entre un
énoncé de physique et le modèle mathématique correspondant.

Problème d’acoustique sous-marine

On considère une source de très petite taille (assimilée à une source ponctuelle omnidirectionnelle
unitaire) et monochromatique en temps. On note f sa fréquence et ω sa pulsation. On suppose en
premier lieu que cette source est placée dans un milieu infini, isotrope et homogène. On note c la
célérité du son dans l’eau.

1.) Écrire l’EDP et les CF correspondant à cette problématique.

2.) Montrer que l’expression suivante est bien solution du problème :

ei(ωt−kr)
p (r, t) =
4πr

Remarque B.1 On admettra qu’il s’agit d’une égalité entre distributions spatiales (la variable
temps n’introduisant pas de difficulté particulière). De même, suite à l’exercice en électrostatique,
on ne considérera que des fonctions test à symétrie sphérique (indépendantes de θ et de ψ).

3.) Donner la fonction de Green du problème.

4.) On s’intéresse maintenant au champ de pression produit par cette même source mais cette
fois-ci, placée dans un domaine semi-infini en immersion (z allant de 0 pour la surface, à l’infini
pour le fond). Les conditions aux limites sont ainsi donnée en z = 0 par les caractéristiques de
réflexion/transmission du dioptre air/eau. On admettra que cette interface peut être considérée
comme molle pour une source située dans l’eau.

La source est immergée une profondeur zs . L’étude de ce type de scénario est d’une utilité relati-
vement pratique puisqu’il se trouve dans les problèmes de propagation d’ondes sonores par grand
fond.

Donner la nouvelle liste des contraintes imposées à la nouvelle solution du problème physique.

Donner l’expression du champ acoustique produit par la source dans ces conditions. Déterminer la
nouvelle fonction de Green du problème. Que se passe-t-il si la source n’est plus ponctuelle ?

131
B. Énoncés des TDs

Corrigé

1)

La pression devra suivre l’équation des onde, la ponctualité spatiale se traduit par un dirac et la
fréquence pure par une exponentielle complexe. Pour un domaine infini, lorsque les sources sont
dans un domaine borné, la condition au frontière est la condition de Sommerfeld. Au final, on a :
� � �
∂2
EDP : Δ − c12 ∂t 2 P = −δo eiωt
(B.12)
CF : condition de Sommerfeld (voir poly)

2.)
ej(ωt−kr)
P (r, θ, φ, t) =
4πr

On vérifie avant tout que la solution proposée vérifie l’équation des ondes inhomogènes

ΔP + k 2 P = −δ (r) · ejωt

Mathématiquement, cette équation n’a de sens que dans l’espace des distributions.

Le système est à symétrie centrale donc indépendant de θ et ϕ :

� �
1 ∂ ∂P
ΔP = 2 r2
r ∂r ∂r

Concrètement, il faut donc vérifier au sens des distributions que :

� �
1 ∂ 2 ∂P
r + k 2 P = −δ (r) · ejωt
r2 ∂r ∂r

Dans premier temps, nous allons traiter la partie �ΔP, ϕ� :

� +∞ � �
�ΔP, ϕ� = 4π P 2rϕ� (r) + r2 ϕ�� (r) dr (B.13)
0
� +∞
ei(ωt−kr) � �
= 4π 2rϕ� (r) + r2 ϕ�� (r) dr (B.14)
0 4πr
� +∞ � +∞
−ikr �
= 2e iωt
e ϕ (r) dr + e iωt
e−ikr rϕ�� (r) dr (B.15)
0 0
� �� � � �� �
A1 A2

Par IPP, on a :
� +∞
� −ikr
�+∞
A1 = −2e iωt
e ϕ (r) 0
+ 2e iωt
ik e−ikr ϕ (r) dr (B.16)
0
� +∞
= −2eiωt ϕ (0) + 2eiωt ik e−ikr ϕ (r) dr (B.17)
0

132
 B.6. Problème complémentaire

Par IPP, on a aussi :

� +∞
� �+∞
A2 = eiωt re−ikr ϕ� (r) 0 − eiωt (1 − ikr) e−ikr ϕ� (r) dr (B.18)
0
� +∞
A1 iωt −ikr �
=0− +e ikre ϕ (r) dr (B.19)
2 0

Par une seconde IPP, on a :

� +∞
A1 � −ikr �+∞
A2 = − + ike iωt
re ϕ (r) 0 − ike iωt
(1 − ikr) e−ikr ϕ (r) dr (B.20)
2 0
� +∞ � +∞
A1
=− + 0 − ikeiωt e−ikr ϕ (r) dr − k 2 eiωt e−ikr rϕ (r) dr (B.21)
2 0 0

Donc :

�ΔP, ϕ� = A1 + A2 (B.22)
� +∞
= −eiωt ϕ (0) − k 2 eiωt e−ikr rϕ (r) dr (B.23)
0

Par ailleurs, on a :
� +∞ � +∞
1 ∂2P ω 2 iωt
− 2 ϕ (r) 4πr 2
dr = e e−ikr ϕ (r) r dr (B.24)
c 0 ∂t2 c2 0

Donc :
�� � � � � � +∞
1 ∂2 ω2
Δ− P, ϕ = −eiωt ϕ (0) + − k 2
e iωt
e−ikr rϕ (r) dr (B.25)
c2 ∂t2 c2 0

Dans le cas général, on ne peut pas aller plus loin. Par contre si on impose k 2 = ω 2 /c2 (ce qui est
le cas classique), alors on a :
�� � �
1 ∂2
Δ − 2 2 P, ϕ = −eiωt ϕ (0) (B.26)
c ∂t
� �
= −δo eiωt , ϕ (B.27)

Pour que notre solution soit la bonne, il faut en plus vérifier les conditions aux frontières. Dans
notre cas (domaine infini), ces conditions ne sont que la condition de Sommerfeld dans toutes les
directions : � �
∂P 1 ∂P
lim r + =0
r→∞ ∂r c ∂t

Cette condition est parfaitement vérifiée car :

� �
∂P 1 jk
= − − ej(ωt−kr)
∂r 4πr2 4πr
et donc :
∂P 1 ∂P 1 j(ωt−kr)
+ =− e
∂r c ∂t 4πr2
La limite quand r devient grand est bien 0.

133
B. Énoncés des TDs

3.) La fonction de Green vérifie :

� � �
1 ∂2
EDP : Δ− c2 ∂t2 P = −δr0 eiωt
(B.28)
CF : condition de Sommerfeld (voir poly)

Par changement de repère, on a :

ei(ωt−k|r−r0 |)
G (r, r0 ) = (B.29)
4π |r − r0 |

4.) On vérifie rapidement que dans le cas du demi-espace, la fonction de Green du problème est :

ej (ωt−k|r−r0 |)
�
ej(ωt−k|r−r0 |)
G (r, r0 ) = −
4π |r − r0 | 4π |r − r0� |

où r0� est l’image par la symétrie du plan (interface) de la position de la source r0 . Physiquement
cela s’interprète comme une image virtuelle de la source.

Dans le cas présent et en pointant la direction z vers le fond, la source est située à la profondeur
zs donnera lieu au champ de pression acoustique :
� √ � � √ �
j ωt−k x2 +y 2 +(z−zs )2 j ωt−k x2 +y 2 +(z+zs )2
e e
P (r) = � − �
2 2
4π x2 + y 2 + (z − zs ) 4π x2 + y 2 + (z + zs )

134