Propagation d’un paquet d’ondes dans un milieu dispersif
Buts : Simuler la propagation d’un paquet d’ondes dans un milieu dispersif et visualiser le phénomène
d’étalement.
1 La relation de dispersion
On considère le système vu en cours d’une corde tendue selon l’axe x, effectuant de petites oscillations suivant
l’axe y et rappelée élastiquement vers la position x = 0. L’équation qui gouverne ces oscillations, dite équation
de Klein-Gordon, est
∂2 y 2
2∂ y
= c − ω20 y,
∂t 2 ∂x 2
avec c = T0 /µ et ω0 = κ/µ. T0 = 500 N est la norme de la tension dans la corde, µ = 10−3 kg/m la masse
p p
linéique de la corde et κ = 0, 1 N/m2 (ou Pa) la constante de rappel par unité de longueur. L’équation ci-dessus
conduit à la relation de dispersion
1q 2
k(ω) = ω − ω20
c
pour ω > ω0 . Nous ne considérerons que des pulsations supérieures à ω0 pour simplifier.
1. Sur Pyzo, importer les packages numpy et matplotlib.pyplot. Rentrer les valeurs numériques et définir les
paramètres utiles. Créer une fonction qui renvoie k(ω).
2 Le spectre du signal
Nous voulons visualiser la propagation d’une onde localisée dans le temps. Nous allons considérer une onde
avec un spectre gaussien, centré sur une pulsation ωm et de largeur ∆ω. Idéalement, il faudrait considérer un
spectre continu, avec une amplitude associée à la pulsation ω valant
(ω − ωm )2
−
a(ω) = e 2∆ω2 .
En pratique, nous allons considérer une somme sur un grand nombre de fréquences dans l’intervalle [ωm −
3∆ω; ωm + 3∆ω] réparties régulièrement, ce qui donnera une bonne approximation du spectre continu.
1
Lycée Laetitia Bonaparte, PSI 2023-2024
2. Créer la fonction a(ω). Créer également un array w qui contient N pulsations, de ωm − 3∆ω à ωm + 3∆ω.
On pourra prendre comme paramètres ωm = 25 rad/s, ∆ω = 5 rad/s et N = 200. Chaque pulsation ωn dans
w a une amplitude associée a(ωn ). Tracer le spectre correspondant.
3. Définir ensuite la fonction
N
X
y(t , x) = a(ωn ) cos [ωn t − k(ωn )x]
n=1
qui décrit la propagation de l’onde avec le spectre précédent. On pourra utiliser une boucle for.
3 Visualisation de l’évolution du paquet d’ondes
Nous allons ensuite visualiser le paquet d’onde à différents instants.
4. Créer un array x pour la coordonnée x. On pourra prendre 2000 valeurs allant de −1000 à 12000 m. Tracer
ensuite y(t , x) pour différents instants t , par exemple t = 0, 7, 15 s.
5. Commenter l’évolution du paquet d’onde : largeur, hauteur, symétrie...
Il est intéressant de comparer l’évolution dans un milieu dispersif avec celle dans un milieu non-dispersif. Pour
cela, on peut définir
ω
k ′ (ω) = ,
c
relation non-dispersive qui correspond au cas où on n’a pas de force de rappel.
6. Créer une fonction y ′ (t , x) avec le même spectre que précédemment, mais basée sur la relation de
dispersion k ′ (ω). Visualiser son évolution au cours du temps et comparer au cas dispersif.
L’équation de Klein-Gordon permet également, et de manière plus intéressante, de décrire la propagation d’une
onde électromagnétique dans un plasma (gaz d’électrons et d’atomes ionisés). En physique des particules, les
bosons de masse non-nulle (bosons W, Z et de Higgs) sont des ondes décrites par l’équation de Klein-Gordon.
Dans ce dernier cas, cela signifie que même le vide peut être un milieu dispersif !
