Dispersion : plasma dilué, effet de peau, paquet d'ondes
Cours · ondes électromagnétiques (MP), chapitre 2 · CPGE MP (2e année) · le début du chapitre, la suite sur Adloun
Dans le vide, l'équation de d'Alembert décrit une propagation non dispersive et non absorbante : toutes les composantes harmoniques d'un signal se déplacent à la même célérité , et un paquet d'ondes garde sa forme. Les milieux matériels s'écartent de cet idéal. Dans un plasma dilué, la vitesse de propagation dépend de la fréquence (dispersion) et une onde de trop basse fréquence ne se propage plus (coupure). Dans un conducteur ohmique, le champ ne pénètre que sur une faible épaisseur, qu'il atténue exponentiellement (effet de peau). Ces deux situations, au programme de MP, illustrent la démarche générale : chercher des ondes planes harmoniques , dont le vecteur d'onde est a priori complexe, et lire dans la relation de dispersion tout le comportement du milieu. On y adjoint la distinction essentielle entre vitesse de phase et vitesse de groupe, qui gouverne la propagation d'un paquet d'ondes réel.
2.1 Relation de dispersion d'un milieu linéaire
2.1.1 Milieu linéaire et décomposition de Fourier
Une équation aux dérivées partielles portant sur un champ est linéaire si toute combinaison linéaire de solutions est encore solution. Elle ne fait intervenir et ses dérivées qu'au premier degré, sans terme constant ni produit du type .
La linéarité autorise la décomposition de Fourier : un signal quelconque étant une superposition d'ondes harmoniques, il suffit d'étudier la propagation de chaque composante indépendamment, puis de sommer.
Méthode : Établir une relation de dispersion
Pour une équation linéaire à coefficients constants :
- injecter l'onde plane progressive harmonique complexe ;
- utiliser les correspondances et ;
- factoriser par : il reste une relation algébrique entre et , la relation de dispersion.
La relation de dispersion est la relation imposée par l'équation de propagation aux ondes harmoniques. Le vecteur d'onde est a priori complexe :
La partie réelle fixe la longueur d'onde et la propagation de la phase ; la partie imaginaire fixe l'atténuation en amplitude (avec la convention et une onde vers les , ).
Sur :
- si dépend de (soit non proportionnel à ) : le milieu est dispersif ;
- si alors que : l'onde se propage en s'atténuant exponentiellement, le milieu est absorbant ;
- si est imaginaire pur (, ) : l'onde est évanescente, son amplitude décroît sans propagation de phase, et sans dissipation (voir le plasma sous la coupure) : décroissance n'est pas absorption.
Un milieu peut être l'un, l'autre, les deux, ou aucun (d'Alembert : réel).
2.1.2 Vitesse de phase
Pour une onde harmonique de pulsation et de vecteur d'onde réel , la vitesse de phase est la vitesse de déplacement d'un plan de phase constante :
C'est la vitesse d'une crête d'une onde monochromatique idéale.
La vitesse de phase peut dépasser (c'est le cas dans le plasma) sans violer la relativité : une onde strictement monochromatique est infinie dans l'espace et le temps, elle ne transporte aucune information. Dans le plasma, l'énergie et l'information voyagent à la vitesse de groupe, inférieure à ; ce n'est pas une règle générale (au voisinage d'une raie d'absorption, peut dépasser ), mais aucun signal ne va jamais plus vite que .
2.2 Onde électromagnétique dans un plasma dilué
2.2.1 Modèle du plasma dilué
Un plasma dilué est un gaz d'électrons libres (densité volumique ) sur un fond d'ions positifs immobiles assurant la neutralité électrique globale. Le plasma est dilué : on néglige les collisions (donc l'absorption dissipative) et l'action de la force magnétique de l'onde sur les électrons (). On étudie une onde plane progressive harmonique transverse électrique .
Démonstration (Conductivité complexe du plasma)
Le principe fondamental de la dynamique appliqué à un électron (charge , masse ), soumis au seul champ électrique de l'onde, s'écrit en notation complexe :
La densité de courant volumique fait apparaître la conductivité complexe
imaginaire pure : le courant est en quadrature avec le champ, donc la puissance volumique moyenne cédée aux électrons . Le plasma dilué ne dissipe pas d'énergie.
2.2.2 Relation de dispersion et pulsation plasma
Démonstration (Relation de dispersion du plasma)
Dans le plasma la charge volumique est nulle (neutralité), . En combinant les équations de Maxwell, l'équation de propagation du champ transverse est
On injecte et , avec :
où l'on a utilisé .
La pulsation plasma est la pulsation propre d'oscillation collective des électrons du plasma. Elle joue le rôle d'une pulsation de coupure : elle sépare un domaine de transparence d'un domaine réactif.
n'existe : l'onde est évanescente. Pour l'onde se propage ; la courbe,
convexe, tend vers l'asymptote non dispersive .2.2.3 Transparence, coupure et ondes évanescentes
La relation donne deux comportements :
- (transparence) : est réel, l'onde se propage. Alors
- (domaine réactif) : , donc avec ; l'onde est évanescente (), sans propagation de phase. Elle ne transporte aucune puissance moyenne : le plasma réfléchit totalement l'onde.
À : , la coupure.
dans le domaine de transparence, et avec ; les deux
tendent vers à haute fréquence.Dans le domaine réactif, . La relation de structure avec imaginaire pur montre que est en quadrature temporelle avec . Le vecteur de Poynting moyen : aucune énergie n'est transportée en moyenne. C'est la signature du transport nul par une onde évanescente, à distinguer de l'onde progressive absorbée (effet de peau, ci-après), qui transporte de l'énergie en la dissipant.
La haute atmosphère est ionisée par le rayonnement solaire : à , soit à et à . Les ondes de fréquence (grandes ondes AM, ondes courtes) sont réfléchies : cela permet les liaisons radio par-delà l'horizon (rebonds successifs sol--ionosphère). Les ondes FM (), la télévision et les liaisons satellites () la traversent.
Un métal est un plasma dense : pour le cuivre donne , soit , dans l'ultraviolet. Comme le visible vérifie , il est réfléchi : c'est l'origine de l'éclat métallique. Au-delà de , le métal devient transparent : c'est ce qu'on observe pour les métaux alcalins (sodium, ) ou l'aluminium (). Le modèle des électrons libres explique l'éclat, pas la couleur : celle du cuivre vient de l'absorption par des électrons liés, qu'il ignore.
2.3 Effet de peau dans un conducteur ohmique
2.3.1 Équation de diffusion du champ
Dans un conducteur ohmique de conductivité réelle (loi d'Ohm locale ), en régime lentement variable, le courant de conduction domine largement le courant de déplacement.
Démonstration (Approximation des régimes quasi stationnaires)
Ordre de grandeur des deux courants : conduction , déplacement . Leur rapport vaut
Pour le cuivre (), ce rapport vaut pour (rayons X). Jusqu'à des fréquences très élevées, on néglige le courant de déplacement devant le courant de conduction. En pratique, c'est la loi d'Ohm à réel elle-même qui cesse de valoir bien avant, dans l'infrarouge (vers pour le cuivre), quand la période de l'onde devient comparable à la durée moyenne entre deux collisions d'un électron : l'effet de peau tel qu'il est décrit ici vaut donc de l'électrotechnique aux micro-ondes.
En négligeant le courant de déplacement, les équations de Maxwell donnent pour le champ dans le conducteur une équation de diffusion :
Ce n'est pas une équation de d'Alembert : il n'y a pas de propagation à célérité définie, mais une diffusion amortie du champ dans le métal.
Démonstration (Établissement de l'équation)
On prend le rotationnel de Maxwell--Faraday :
Le conducteur est localement neutre (), et Maxwell--Ampère sans courant de déplacement s'écrit . Il vient , d'où le résultat.
2.3.2 Épaisseur de peau
Démonstration (Relation de dispersion et épaisseur de peau)
On cherche . L'équation de diffusion donne
On extrait la racine de partie imaginaire négative (onde amortie vers ). Comme , on a , soit
L'épaisseur de peau (ou profondeur de pénétration) d'un conducteur ohmique est
Le champ s'écrit : il est à la fois amorti (facteur ) et propagatif (phase ). Sur une distance l'amplitude est divisée par . Cette atténuation s'accompagne d'une dissipation par effet Joule (contrairement à l'onde évanescente du plasma).
, , : à , à , à . À haute fréquence le courant est confiné dans une pellicule superficielle : la résistance effective augmente (on utilise des câbles multibrins, fil de Litz), et un blindage métallique mince suffit à confiner les ondes.
Quand , : le champ ne pénètre pas. On retrouve le conducteur parfait, dans lequel ; les courants et charges se répartissent en surface. C'est l'idéalisation utilisée pour les cavités et les guides d'ondes.
2.4 Paquet d'ondes et vitesse de groupe
Les deux programmes demandent de calculer la vitesse de groupe à partir de la relation de dispersion. Le programme de MP demande en outre de l'associer à la propagation de l'enveloppe d'un paquet d'ondes dans un milieu faiblement dispersif ; celui de MPI se limite à en citer une interprétation sur un exemple qualitatif (battement de deux ondes, vitesse de l'énergie). La notion de paquet d'ondes continu et son étalement sont donc des compléments en MPI.
2.4.1 Notion de paquet d'ondes
Un signal physique réel est limité dans le temps et dans l'espace : c'est un paquet d'ondes, superposition d'ondes harmoniques dont les pulsations se répartissent autour d'une valeur centrale sur une largeur .
La durée d'un paquet d'ondes et la largeur de son spectre sont liées par une inégalité de type Fourier :
Un paquet bref est spectralement large ; un signal quasi monochromatique ( petit) est nécessairement long. C'est l'analogue classique des inégalités de Heisenberg (relation retrouvée en physique quantique).
2.4.2 Vitesse de groupe
Considérons la superposition de deux ondes de pulsations voisines et de vecteurs d'onde , d'égale amplitude :
La porteuse se déplace à ; l'enveloppe se déplace à . En passant à la limite continue :
La vitesse de groupe d'un paquet d'ondes centré sur est
C'est la vitesse de propagation de l'enveloppe du paquet, donc de l'énergie et de l'information dans un milieu transparent.
Méthode : Calculer et
À partir de la relation de dispersion :
- exprimer (ou ) dans le domaine propagatif ( réel) ;
- vitesse de phase : ;
- vitesse de groupe : , soit .
Milieu non dispersif .
2.4.3 Étalement d'un paquet dans un milieu dispersif
Dans un milieu dispersif, dépend de . Les différentes composantes spectrales du paquet se propageant à des vitesses différentes, le paquet s'étale : sa largeur croît, son amplitude crête décroît (à énergie conservée). L'étalement est d'autant plus rapide que la courbure (dispersion de vitesse de groupe) est grande. Le programme se limite au régime faiblement dispersif, où l'on assimile le paquet à une enveloppe rigide se déplaçant à sur des distances modérées.
Capacité numérique : Propagation d'un paquet d'ondes dispersif
Complément : cette capacité numérique figure aux programmes de PC et de PSI ; elle n'est pas exigible en MP ni en MPI. On simule la propagation d'un paquet gaussien dans un plasma (), en unités réduites (). Le paquet initial, écrit avec la convention du cours , est ; il est complexe, ce qui ne garde que les ondes qui avancent vers les (un signal réel contiendrait aussi la moitié qui part vers les ). On le décompose par transformée de Fourier discrète, on multiplie chaque composante par , on reconstruit le signal et l'on trace sa partie réelle (le signal physique) et son module (l'enveloppe). On visualise l'étalement et l'on compare la position du centre du paquet à .
import numpy as np
import matplotlib.pyplot as plt
c, omega_p = 1.0, 1.0 # unites reduites
k0, sigma = 2.0, 2.0 # vecteur d'onde central, largeur initiale
def omega(k): # relation de dispersion du plasma (paire en k)
return np.sqrt(omega_p**2 + (c*k)**2)
N = 8192
z = np.linspace(-200.0, 800.0, N, endpoint=False) # boite longue : pas de repliement
dz = z[1] - z[0]
# Convention du cours exp(j(omega t - k z)) : a t = 0, le paquet vaut
# exp(-z^2/2 sigma^2) exp(-j k0 z).
s0 = np.exp(-z**2/(2*sigma**2)) * np.exp(-1j*k0*z)
S0 = np.fft.fft(s0)
q = 2*np.pi*np.fft.fftfreq(N, d=dz) # la FFT decompose sur exp(+j q z), soit q = -k
def paquet(t): # chaque composante recoit le facteur exp(+j omega t)
return np.fft.ifft(S0 * np.exp(1j*omega(q)*t))
def centre_largeur(s): # centre et largeur quadratique de |s|^2
p = np.abs(s)**2
zm = np.sum(z*p)/np.sum(p)
return zm, np.sqrt(np.sum((z - zm)**2*p)/np.sum(p))
v_phi, v_g = omega(k0)/k0, c**2*k0/omega(k0)
print("v_phi = %.3f v_g = %.3f v_phi*v_g = %.3f" % (v_phi, v_g, v_phi*v_g))
_, larg0 = centre_largeur(paquet(0.0))
fig, ax = plt.subplots(3, 1, figsize=(8, 6), sharex=True, sharey=True)
for a, t in zip(ax, [0.0, 90.0, 180.0]):
s = paquet(t)
zm, larg = centre_largeur(s)
print("t = %3.0f : centre %5.1f (v_g*t = %5.1f), largeur x %.2f, amplitude %.2f"
% (t, zm, v_g*t, larg/larg0, np.abs(s).max()))
a.plot(z, s.real, color="#0f7387", lw=0.6) # signal physique
a.plot(z, np.abs(s), color="crimson", lw=1.2) # enveloppe
a.axvline(v_g*t, color="crimson", ls="--", lw=0.8)
a.set_ylabel("t = %.0f" % t)
ax[-1].set_xlim(-20, 200); ax[-1].set_xlabel("z")
ax[0].set_title("Etalement d'un paquet d'ondes dans un plasma")
plt.tight_layout(); plt.savefig("paquet_dispersif.png", dpi=130)
Le programme affiche :
v_phi = 1.118 v_g = 0.894 v_phi*v_g = 1.000
t = 0 : centre -0.0 (v_g*t = 0.0), largeur x 1.00, amplitude 1.00
t = 90 : centre 79.8 (v_g*t = 80.5), largeur x 2.58, amplitude 0.68
t = 180 : centre 159.7 (v_g*t = 161.0), largeur x 4.87, amplitude 0.51
Observations. Un unique paquet avance vers la droite ; son centre suit à près (trait tireté de la figure qui suit), alors que les crêtes de la porteuse vont plus vite () et traversent l'enveloppe d'arrière en avant ; on vérifie . Surtout, le paquet s'étale : sa largeur est multipliée par en et par en , et son amplitude crête est divisée par deux. La porteuse n'est plus régulière : les composantes de grand , plus rapides ( croît avec ), sont passées en tête. Le choix des paramètres compte : avec un paquet plus large () centré sur , le même code ne montre aucun étalement visible (largeur multipliée par seulement en ). L'étalement est d'autant plus marqué que le spectre du paquet est large ( petit) et que la courbure est grande (ici, près de la coupure).
et : signal (trait fin) et enveloppe (trait épais). Le trait tireté marque la position . Le paquet s'élargit et s'affaisse :
c'est l'étalement dispersif.2.5 Synthèse du chapitre
- Méthode : (, ) relation de dispersion . phase et dispersion ; atténuation (absorption si , évanescence si ).
- Plasma dilué : (imaginaire pure, non dissipatif) ; , .
- : transparence, , , ; : onde évanescente, réflexion totale, Poynting moyen nul. Ionosphère : quelques MHz.
- Effet de peau (conducteur ohmique, courant de déplacement négligé) : équation de diffusion , , ; onde amortie et propagative, dissipation Joule ().
- Vitesses : (crêtes, peut dépasser ) ; (enveloppe, énergie, information ; dans le plasma).
- Paquet d'ondes : ; en milieu dispersif il s'étale ().
Exercices du chapitre
- Lire une relation de dispersion : dispersion, absorption, coupure corrigé
- Un électron dans l'onde : conductivité complexe, approximations justifiées, énergie cinétique corrigé
- Une formule, deux mondes : l'ionosphère et la transparence ultraviolette des métaux alcalins corrigé
- Où est l'énergie d'une onde dans un plasma ? Vitesse de l'énergie et onde évanescente corrigé
- Effet de peau : sous-marins, papier d'aluminium et cave à vin corrigé
- Deux ondes voisines : la porteuse traverse l'enveloppe corrigé