Paquet d'onde et dispersion

2.a. Ondes périodiques

Considérons tout d'abord une onde plane sinusoïdale (en notation complexe) :

$$u(x,t)=Aexpi(kx-\omega t)\text{(1)}$$

On appelle relation de dispersion la fonction k(ω) donnant le nombre d'onde en fonction de la pulsation. La relation de dispersion dépend du milieu dans lequel l'onde de propage. On utilisera les relations de dispersion suivantes :

• Vide $$k\left(\omega \right)=\frac{\omega }{c}\text{(2)}$$
• Milieu à coupure (plasma) $$k\left(\omega \right)=\frac{1}{c}\sqrt{{\omega }^2-{\omega }_c^2}\text{(3)}$$
• Milieu transparent (indice optique) $$k\left(\omega \right)=\frac{\omega }{c}(1+b{\omega }^2)\text{(4)}$$
• Onde d'une particule non relativiste libre (onde de De Broglie) $$k\left(\omega \right)=\sqrt{\frac{2m\omega }{\hslash }}\text{(5)}$$

Soit une onde générée à partir d'une source périodique, de période T. La pulsation fondamentale est définie par :

$${\omega }_1=\frac{2\pi }{T}\text{(6)}$$

La série de Fourier permet d'écrire l'onde comme une somme d'ondes sinusoïdales :

$$u(x,t)=\sum _{n=1}^P{A}_{n}expi(k\left(n{\omega }_1\right)x-n{\omega }_{1}t)\text{(7)}$$

L'harmonique de rang n est une onde sinusoïdale de pulsation nω₁, dont le nombre d'onde est k(nω₁). L'amplitude A_n de cette harmonique est en général un nombre complexe.

La vitesse de phase est la vitesse de déplacement d'une surface de phase constante :

$$V_{\phi }=\frac{\omega }{k\left(\omega \right)}\text{(8)}$$

À l'exception du cas d'une onde électromagnétique dans le vide, la vitesse de phase dépend de la pulsation. Cela signifie que les différentes harmoniques ne se propagent pas exactement à la même vitesse. En conséquence, la forme de l'onde est en général modifiée au cours du temps. Ce phénomène est appelé dispersion de l'onde.

2.b. Paquet d'onde

En principe, un paquet d'onde est une onde localisée dans le temps (et dans l'espace), c'est-à-dire une impulsion plus ou moins brève. Il ne s'agit pas d'une onde périodique. On peut néanmoins construire une onde périodique constituée d'une succession périodique de paquets d'onde, de la manière suivante.

Pour simplifier la description, on pose T=1, ce qui revient à identifier l'indice n de l'harmonique à la fréquence. On considère un spectre constitué d'une fréquence centrale f (un nombre entier) assez élevée (au moins 100), et de P raies de part et d'autres de cette fréquence. Le spectre est donc constitué des 2P+1 raies f-P, f-P+1, ..f+P. Toutes les autres harmoniques sont nulles. On suppose de plus que P est petit devant f. Le nombre de raies est impair si P est entier, mais on peut choisir P demi-entier pour avoir un nombre de raies pair.

Comme l'indice n dans la série de Fourier reste relativement proche de f, on fait le développement limité suivant (à l'ordre 1) :

$$k\left(n{\omega }_1\right)=k\left(f{\omega }_1\right)+(n-f){\omega }_{1}k\text{'}\left(f{\omega }_1\right)\text{(9)}$$

où l'on a introduit la fonction k'(ω), dérivée de la relation de dispersion.

On pose q=n-f, un indice qui varie de -P à +P dans la somme. La phase de l'harmonique de rang n s'écrit alors :

$$\phi (x,t)=k\left(f{\omega }_1\right)x-f{\omega }_{1}t+q{\omega }_{1}k\text{'}\left(f{\omega }_1\right)x-q{\omega }_{1}t\text{(10)}$$

Après factorisation, on obtient donc l'expression suivante de l'onde :

$$u(x,t)=expi(k\left(f{\omega }_1\right)x-f{\omega }_{1}t)\sum _{q=-P}^{P}expiq{\omega }_1(k\text{'}\left(f{\omega }_1\right)x-t)\text{(11)}$$

L'exponentielle en facteur est une onde sinusoïdale se propageant à la vitesse de phase de la fréquence f. La somme constitue une enveloppe, dont la vitesse de propagation est :

$$V_g=\frac{1}{k\text{'}\left(f{\omega }_1\right)}=\frac{d\omega }{dk}\left(f{\omega }_1\right)\text{(12)}$$

Cette vitesse est la vitesse de groupe. C'est la vitesse de déplacement des paquets d'onde. Dans un milieu dispersif, la vitesse de groupe est différente de la vitesse de phase.

Nous avons défini le paquet d'onde à partir de son spectre temporel : les fréquences temporelles ω sont multiples d'une fréquence fondamentale et les fréquences spatiales k sont calculées avec la relation de dispersion k(ω). On peut aussi définir le paquet d'onde à partir de son spectre spatial : les fréquences spatiales k sont alors multiples d'une fréquence fondamentale et les fréquences temporelles sont calculées avec la relation ω(k). Dans le vide, les deux descriptions sont équivalentes. En radio-transmission, on utilise plutôt un spectre en fréquence correspondant au spectre du signal émis.

import numpy
import math
from matplotlib.pyplot import *
import matplotlib.animation as animation
from onde import Onde
            

4.a. Dispersion avec coupure

On choisit une fréquence de coupure de 50 :

onde = Onde(suivi=True,dispersion='coupure',coupure=50)
                

On définit un paquet d'onde avec une fréquence centrale de 100 et 10 fréquences de part et d'autre, avec un fenêtrage de Hamming :

onde.paquet(100,10,window='hamming') 
                

Voyons tout d'abord le spectre :

figure(figsize=(12,5))
stem(onde.freq,onde.amp)
xlabel("f")
ylabel("A")
grid()
axis([0,150,0,1.0])
                

figA.pdf

Voici l'onde à l'instant initial :

temps = 0.0
N = 5000
xmin = -1.0
xmax = 1.0
(x,y) = onde.echantillons(xmin,xmax,temps,N)
figure(figsize=(12,5))
plot(x,y)
xlabel("x")
ylabel("u")
axis([xmin,xmax,-15,15])
grid()
                

figB.pdf

Les paquets se déplacent à la vitesse de groupe. Comme le suivi à la vitesse de groupe est activé, on peut voir la forme des paquets un instant plus tard :

temps=3.0
(x,y) = onde.echantillons(xmin,xmax,temps,N)
figure(figsize=(12,5))
plot(x,y)
xlabel("x")
ylabel("u")
axis([xmin,xmax,-15,15])
grid()
                 

figC.pdf

On observe un étalement des paquets.

Il est possible d'obtenir une animation montrant l'évolution dans le temps des paquets. Voici comment procéder :

temps = 0.0
dt = 0.05
N = 5000
xmin = -2.0
xmax = 2.0

(x,y) = onde.echantillons(xmin,xmax,temps,N)

fig, ax = subplots()
line, = ax.plot(x,y)
ax.grid()

def animate(i):
    global temps,xmin,xmax,N
    temps += dt
    (x,y) = onde.echantillons(xmin,xmax,temps,N)
    line.set_xdata(x)
    line.set_ydata(y)
    return line,

ani = animation.FuncAnimation(fig,animate,1000,interval=40)
    
show()
                 

La fonction animate est appelée 1000 fois, toutes les 40 millisecondes. À chaque appel, on incrémente le temps de dt. En modififant cette valeur, on peut changer la vitesse de l'animation.

L'animation permet de bien voir la différence entre la vitesse de groupe et la vitesse de phase, et permet de suivre l'étalement progressif des paquets d'onde.

4.b. Onde de De Droglie

L'onde de De Broglie associée à une particule libre définit une amplitude de probabilité. Son module au carré donne la probabilité de trouver la particule à l'abscisse x.

onde = Onde(suivi=True,dispersion='debroglie')
onde.paquet(100,10,window='hamming')
                 

Voici la densité de probabilité à l'instant initial :

temps = 0.0
N = 5000
xmin = -1.0
xmax = 1
(x,y) = onde.echantillons_proba(xmin,xmax,temps,N)
figure(figsize=(12,5))
plot(x,y)
xlabel("x")
ylabel("u")
axis([xmin,xmax,-15,15])
grid()
                 

figD.pdf

Voyons la densité de probabilité un instant plus tard :

temps=3.0
(x,y) = onde.echantillons_proba(xmin,xmax,temps,N)
figure(figsize=(12,5))
plot(x,y)
xlabel("x")
ylabel("u")
axis([xmin,xmax,-15,15])
grid()
                 

figE.pdf