Electromagnetisme Dans Le Vide

Électromagnétisme dans le vide
Préparation à l'agrégation de l'ENS de Lyon 2020-2021
Jérémy Ferrand
2
Bibliographie
Bibliographie classique
• Cap prépa - Physique - 1ère et 2nd année
• H. Gié - Électromagnétisme - 1ère et 2nd année
• P. Roux - Électromagnétisme
• Bertin Faroux Renault - Électromagnétisme Tomes 1, 3 et 4
• R. Feynman - Électromagnétisme 1
• H Prépa - Électromagnétisme - 1ère et 2nd année
• H Prépa - Ondes
Bibliographie pour aller plus loin

• J. D. Jackson - Électrodynamique classique
• Landeau - Électromagnétisme des milieux continus
• Purcell - Electromagnétisme - Cours de Berkley
Bibliographie pour méthodes mathématiques

• W. Appel - Mathématiques pour la physique
• A Alastuey - Physique et outils mathématiques
Épreuves d'agrégation lié à l'électromagnétisme

• Épreuve A 2018 : Interaction de Wan der Waals. Dipôles et interactions entre eux.
• Épreuve A 2013 : Équations de Maxwell, Force de Lorentz, Densité d'énergie électromagnétique,
Vecteur de Poynting, Condensateur, OPPM, Induction, Conducteur parfait, Pression de Radiation,
Rayonnement Dipolaire.
• Épreuve A 2011 : Modèle de l'électron élastiquement lié.
• Épreuve A 2008 : Contribution des noms de l'électromagnétisme, Équations de Maxwell, Conden-
sateur plan, Solénoïde inni, ARQS, Énergie électromagnétiques, Potentiels retardés, Propagation
dans un câble coaxial.
• Épreuve A 2007 : Dipôle magnétique, propagation dans un plasma, Loi d'Ohm, Magnétohydrodyna-
mique.
• Épreuve A 2003 : Équations de Maxwell, Potentiel et potentiel vecteur, Équations sur le potentiel,
Jauge de Lorentz, Force de Lorentz, Fil rectiligne chargé, Rayonnement dipolaire et champs à partir
de A, Énergie, Antenne, Modèle de l'électron élastiquement lié, propagation guidée et bre optique.
• Épreuve A 2001 : Force de Lorentz, Équations de Maxwell, Conservation de la charge, Énergie
volumique, Conservation de l'énergie et vecteur de Poynting, Interaction lumière matière, Magnétisme
classique atomique, RMN classique.
3
4
Table des matières
1 Rappels mathématiques 9
1.1 Analyses vectorielles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.1.1 Champ scalaire, champ vectoriel . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.1.2 Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.3 Divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.1.4 Rotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.1.5 Laplacien scalaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.1.6 Laplacien vectoriel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.2 Relations vectorielles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.1 Théorème de Stokes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.2 Théorème de Green-Ostrogradsky . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.3 Conséquences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.2.4 Égalités vectorielles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3 Diérents systèmes de coordonnées . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3.1 Cartésien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.3.2 Cylindrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.3.3 Sphérique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4 Fonction de Green . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2 Électrostatique 19
2.1 Introduction et historique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2 Cadre de l'électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.3 Approche intégrale et loi de Coulomb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3.1 Loi de Coulomb (1780) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3.2 Champ de Coulomb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3.3 Répartition continue de charges . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.4 Liens avec les équations de Maxwell . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.4.1 Compatibilité avec la force de Coulomb . . . . . . . . . . . . . . . . . . . . . . . . 23
2.4.2 Équivalence des deux approches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.5 Théorème de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.6 Structure du champ électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.6.1 Invariance et principe de Curie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.6.2 Exemples d'invariants . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.6.3 Propriétés de symétrie des champs électriques . . . . . . . . . . . . . . . . . . . . . 26
2.6.4 Calculs de champ électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.7 Potentiel électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.7.1 Dénition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.7.2 Exemples de calculs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.7.3 Signication physique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
5
6 TABLE DES MATIÈRES
2.7.4 Développement multipolaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

2.8 Équation de Poisson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.8.1 Mise en équation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.8.2 Application à l'écrantage de Debye [Exercice] . . . . . . . . . . . . . . . . . . . . . 33
2.8.3 Retour à la forme intégrale : Fonction de Green . . . . . . . . . . . . . . . . . . . . 35
2.9 Énergie des distributions de charges discrètes . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.9.1 Théorème de l'extrémum du potentiel . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.9.2 Énergie d'une distribution discrète . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.10 Énergie des distributions continues . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.11 Stockage de l'énergie électrostatique : les condensateurs . . . . . . . . . . . . . . . . . . . . 39
3 Magnétostatique 41
3.1 Introduction et rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.2 Lois locales et lois intégrales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.2.1 Lois locales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.2.2 Conséquences : lois globales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.2.3 Théorème d'Ampère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.3 Symétries et calculs de champs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.3.1 Symétries du champ magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.3.2 Exemples de calculs du champ B . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.4 Dipôle magnétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4 Relations de passage 51
4.1 Relation de passage pour le champ électrique . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.2 Relation de passage pour le champ magnétique . . . . . . . . . . . . . . . . . . . . . . . . 52
5 Approximation des régimes quasi-stationnaires 55

5.1 Équations de Maxwell et potentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.1.1 Pourquoi le courant de déplacement ? . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.1.2 D'où proviennent 0 et µ0 ? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.1.3 Équations de propagation des potentiels et des champs . . . . . . . . . . . . . . . . 57
5.1.4 Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.2 Démarche naïve de l'ARQS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
5.2.1 Méthode 1 : en utilisant le calcul de magnétostatique . . . . . . . . . . . . . . . . . 61
5.2.2 Méthode 2 : avec des ordres de grandeur . . . . . . . . . . . . . . . . . . . . . . . . 61
5.2.3 Contre exemple et erreurs commises . . . . . . . . . . . . . . . . . . . . . . . . . . 62
5.3 Démarche rigoureuse de l'ARQS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
5.3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
5.3.2 Régime quasi-stationnaire électrique ( 1, ξ 1) . . . . . . . . . . . . . . . . . . 63
5.3.3 Régime quasi-stationnaire magnétique ( 1, ξ 1) . . . . . . . . . . . . . . . . . 64
5.3.4 Remarques supplémentaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
6 Conducteurs 67
6.1 Régime statique - Conducteurs à l'équilibre . . . . . . . . . . . . . . . . . . . . . . . . . . 67
6.2 Régime dynamique - Conduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
6.3 Électroneutralité d'un conducteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
6.4 Retour sur l'ARQS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
6.5 Ondes et métaux conducteurs : eet de peau . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6.5.1 Métal imparfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6.5.2 Limite du métal parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
6.5.3 Force exercée sur le conducteur par le champ . . . . . . . . . . . . . . . . . . . . . 74
TABLE DES MATIÈRES 7
6.6 Magnétorésistance et eet Hall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75

6.6.1 Magnétorésistance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.6.2 Eet Hall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
6.7 Plasma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
7 Guidage d'ondes 81
7.1 Réexion sur un métal parfait : nécessité d'une onde rééchie . . . . . . . . . . . . . . . . 81
7.2 Réexion sous incidence oblique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
7.3 Guide d'ondes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
7.3.1 Plan-plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
7.3.2 Extension au guide rectangulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
8 Induction électromagnétique 89
8.1 Énoncé historique de la loi de Faraday . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
8.2 Cas de l'induction de Neumann . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
8.3 Induction de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
8.4 Transformations des champs E et B suivant les référentiels . . . . . . . . . . . . . . . . . . 91
8.4.1 Cadre de l'étude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
8.4.2 Transformation des champs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
8.5 Cadre général de l'induction de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
8.6 Courants de Foucault . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
9 Rayonnement dipolaire 103

9.1 Cadre de l'étude et hypothèses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
9.2 Expression du champ électromagnétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
9.3 Comportement énergétique du rayonnement . . . . . . . . . . . . . . . . . . . . . . . . . . 106
9.3.1 Le vecteur de Poynting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
9.3.2 Cas du rayonnement dipolaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
9.4 Rayonnement d'antenne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
9.4.1 Antenne électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
9.4.2 Antenne réelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
9.5 Rayonnement de la matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
9.5.1 Expérience introductive . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
9.5.2 Un modèle simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
9.5.3 Auto-consistance du modèle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
8 TABLE DES MATIÈRES
Chapitre 1
Rappels mathématiques
Ce premier chapitre consiste simplement à donner tous les rappels mathématiques nécessaires pour ce
cours d'électromagnétisme. En eet le formalisme de l'électromagnétisme nécessite des notions sur les
champs scalaires et vectoriels ainsi que toute l'analyse vectorielle de base et leurs signications physiques.
Nous eectuons également un rappel sur les fonctions de Green qui servira dans certains problèmes et
calculs.
Les vecteurs seront notés systématiquement en gras an d'alléger la lecture du document.
1.1 Analyses vectorielles

1.1.1 Champ scalaire, champ vectoriel
Tout d'abord, donnons quelques dénitions et caractéristiques des champs.
Un champ d'une certaine grandeur U dans un domaine de l'espace D à un instant donné t est
l'ensemble des grandeurs U aux divers points M de D. Ainsi un champ est une fonction à la fois
de la position r(x, y, z) et du temps t. On note alors ce champ U (x, y, z, t) ou U (r, t). On pourra
prendre comme exemple le champ de pression P (r, t) de température T (r, t) ou de vitesse v(r, t)
d'un uide.
On distinguera les champs scalaires dont la quantité physique peut être représenté par un scalaire
(pression, température) et les champs vectoriels ou la donnée d'un vecteur à 2 ou 3 coordonnées est
nécessaire pour décrire la quantité physique (vitesse).
On parle de champ uniforme lorsqu'il ne dépend pas de l'espace et permanent ou stationnaire

lorsqu'il ne dépend pas du temps.
On appelle ligne de champ, à un instant donné t, la courbe tangente en chaque point M de

l'espace au champ vectoriel V(r, t). Un ensemble de lignes de champ est alors appelé carte de
champ.
Enn on appelle tube de champ d'un champ vectoriel V(r, t), à un instant donné t, la surface
formée par l'ensemble des lignes de champ V qui s'appuient sur une courbe fermée.
9
10 CHAPITRE 1. RAPPELS MATHÉMATIQUES
1.1.2 Gradient
Une quantité importante pour un champ vectoriel V(r, t) est la circulation d'un point A vers
un point B le long de la courbe Γ qui est l'intégrale calculée le long de cette courbe :
ˆ B
C= V.ds (1.1)
A(Γ)
où ds est l'abscisse curviligne qui paramètre Γ. Cette circulation dépend a priori de la courbe Γ.
Un champ est dit à circulation conservative si sa circulation entre deux points ne dépend pas
du chemin suivi.
Soit un champ scalaire U (x, y, z, t), on dénit le gradient de U qu'on note ∇U le champ vectoriel dont les
coordonnées cartésiennes sont les dérivées partielles de U par rapport aux coordonnées correspondantes :
∂U ∂U ∂U
∇U = ex + ey + ez (1.2)
∂x ∂y ∂z
Notons deux choses :

• Le gradient dénit en tout point de l'espace un vecteur ∇U à partir de U . Ainsi le gradient n'est pas
juste un vecteur mais un opérateur.
• La dénition n'est pas la même dans d'autres systèmes de coordonnées (se reporter à la n du
chapitre) mais les résultats sur le gradient ne changeront pas.
Le gradient représente la variation d'un champ scalaire pour un déplacement élémentaire :
dU = dr.∇U (1.3)
Ainsi on a que la variation d'un champ scalaire est égale à la circulation de son gradient.
La conséquence directe de ce résultat est que le gradient est toujours perpendiculaire aux surfaces
où le champ scalaire est constant.
Il est immédiat qu'un champ uniforme à un gradient nul et que réciproquement un champ de gradient nul
dans tout un domaine est uniforme dans ce domaine. Le gradient caractérise ainsi le caractère non-uniforme
d'une champ scalaire U . Le gradient est toujours orienté dans le sens croissant de U .
1.1.3 Divergence
Comme la circulation pour un champ scalaire on peut dénir une grandeur pour un champ vectoriel
important en électromagnétisme. Si on considère une surface Σ s'appuyant sur une courbe fermée orientée
Γ.
Le ux d'un champ vectoriel V à travers Σ comme la somme de tous les ux élémentaires à travers
une surface élémentaire dS orienté :
¨ ¨
Φ= dΦ = V.dS (1.4)
Σ Σ
Le sens du vecteur de surface élémentaire dS se déduit de l'orientation de la courbe par la règle de la

main droite ou du tire-bouchon.
1.1. ANALYSES VECTORIELLES 11
Il existe des champs qui sont dits à ux conservatif dans un domaine D de l'espace si leur ux à travers
toute surface fermée Σ appartenant à D est nul :
‹
∀Σ ∈ D Φ= V.dS = 0 (1.5)
Σ
À présent, calculons le ux d'un champ vectoriel V(x, y, z) sortant d'un cube innitésimal de volume
dτ situé en un point M (x, y, z) représenté sur la Figure 1.1.
Figure 1.1 Calcul du ux sortant d'un cube innitésimal.

Pour cela, calculons le ux sortant à travers chaque face du cube. Si on prend la face marqué 1 et qu'on
calcule le ux à travers cette surface. Comme la taille du cube tend vers 0, on peut prendre en première
approximation la valeur du champ de vecteur au centre de la face an de calculer le ux. On a ainsi :

dy
φ1 = −n.V x, y − , z dxdz (1.6)
2
De même pour la face 2 :

dy
φ2 = n.V x, y + , z dxdz (1.7)
2
À présent si on eectue la somme des deux ux et qu'on écrit la diérence au premier ordre :

dy dy
Φ1 + Φ2 = n. V x, y + , z − V x, y − , z dxdz (1.8)
2 2

∂V(x, y, z)
= n. dy dxdz (1.9)
∂y
∂Vy (x, y, z)
= dτ (1.10)
∂y
où Vy = V.n est la composante du champ selon ey .

On peut alors répéter ce calcul pour les faces 3 et 4 ainsi que 5 et 6 (avec des vecteurs normaux n
diérents) et on peut sommer tous les ux :

∂Vx ∂Vy ∂Vz
Φ = φ1 + φ2 + φ3 + φ4 + φ5 + φ6 = + + dτ = ∇.Vdτ (1.11)
∂x ∂y ∂z
où ∇.V est appelé divergence du champ vectoriel V.
L'opérateur divergence au point M peut être interprété comme "l'écoulement" vers l'extérieur
par unité de volume du champ vectoriel V au voisinage de M .
De plus comme le ux sortant total sortant d'un volume ni est la somme des ux sortants de chacune
de ses parties on peut intégrer la divergence sur le volume entier. Ainsi pour un volume V délimité par
une surface Σ on peut écrire :
¨ ˚
V.n dΣ = ∇.V dτ (1.12)
Σ τ
1.1.4 Rotationnel
Revenons sur la circulation d'un champ de vecteurs sur une courbe fermé Γ. Si on découpe cette courbe
comme sur la Figure 1.2, il est assez évident que la circulation le long de la courbe Γ est alors la somme
Γ1 + Γ2 . En eet, la circulation sur Γ1 vaut la circulation sur Γa plus celle sur Γab . Nous avons de même
pour Γ2 mais avec un signe opposé pour celle sur Γab . Ainsi en sommant la circulation sur Γ1 et Γ2 on
obtient que le contour Γ car le découpage Γab disparait.
Γb
Γ2
Γab
Γa
Γ1 ds2
ds1
Figure 1.2 Calcul de la circulation en découpant la courbe en deux.

Prenons à présent une surface quelconque S qui s'appuie sur la courbe Γ et découpons cette surface en
petits carrés (assez petit pour qu'on puisse considérer qu'il s'agit bien de carrés plats). On obtient alors
un quadrillage de la surface et il nous reste à calculer la circulation sur un carré.
Prenons le schéma de la Figure 1.3. On suppose tout d'abord que le carré est assez petit pour que le
champ soit constant sur tout le carré. La courbe est orienté dans le sens trigonométrique et on calcule
facilement la circulation :
˛
dy dx dy dx
V.ds = Vx x, y − dx + Vy x + , y dy − Vx x, y + dx − Vy x − , y dy
2 2 2 2

dy dy dx dx
= − Vx x, y + − Vx x, y − dx + Vy x + , y − Vy x − , y dy
2 2 2 2

∂Vy ∂Vx
= − dxdy
∂x ∂y
où dxdy n'est que l'aire du carré considéré. On dénit alors un autre opérateur pour rendre compte de
cette circulation, le rotationnel dont la composante selon z (qui est la coordonnée normale au carré) s'écrit
justement
∂Vy ∂Vx
(∇ ∧ V)z = −
∂x ∂y
1.1. ANALYSES VECTORIELLES 13
y
y + dy/2
4 2
V(x,y)
1
y - dy/2
ds
x - dx/2 x + dx/2
x
Figure 1.3 Calcul de la circulation sur un carré.

Si on écrit le rotationnel pour toutes les composantes on a :

∂Vz ∂Vy ∂Vx ∂Vz ∂Vy ∂Vx
∇∧V = − ; − ; − (1.13)
∂y ∂z ∂z ∂x ∂x ∂y
L'expression du rotationnel dans d'autres systèmes de coordonnées n'est pas la même il faut faire atten-
tion.
On a donc trouvé que la circulation de tout vecteur V le long d'un carré innitésimal est le produit de
la composante normale à la surface du rotationnel de V par l'aire du carré. La circulation de la courbe
Γ peut être maintenant relié au rotationnel du champ de vecteur. En eet il nous reste plus qu'à somme
toutes les boucles et on trouve le théorème (de Stokes) :
˛ ¨
V.ds = (∇ ∧ V)n dS (1.14)
Γ S
où n est le vecteur unitaire normal en tout point à la surface S .
1.1.5 Laplacien scalaire

Avec les opérateurs qu'on vient de voir il est possible de construire un nouvel opérateur rencontré en
électromagnétisme. Si on prend la divergence du gradient d'un champ scalaire on montre que :
∂2U ∂2U ∂2U

∇.(∇U ) = ∇2 U = ∆U = + + (1.15)
∂x2 ∂y 2 ∂z 2
L'opérateur ∇2 est appelé Laplacien (scalaire dans ce cas car s'applique à un champ scalaire).
1.1.6 Laplacien vectoriel

On peut dénir la même chose pour un champ vectoriel mais il est pas aisé de savoir quel que le gradient
d'un champ de vecteur... Par analogie, mais uniquement en coordonnées cartésiennes, on peut le
dénir ainsi :
∇.(∇V) = ∆V = (∆Vx )ex + (∆Vy )ey + (∆Vz )ez (1.16)
Mais ceci n'est pas une bonne dénition car elle dépend d'un système de coordonnées. Pour trouver une
dénition indépendante du système de coordonnées il faut passer par le développement du double produit
vectoriel :
A ∧ (B ∧ C) = B(A.C) − (A.B)C (1.17)
Et donc :
∇ ∧ (∇ ∧ V) = ∇(∇.V) − (∇.∇)V = ∇(∇.V) − ∆V (1.18)
Ainsi on a une dénition intrinsèque :
∆V = ∇(∇.V) − ∇ ∧ (∇ ∧ V) (1.19)
1.2 Relations vectorielles

Certaines relations vectorielles sont utiles en électromagnétisme. Même si nous en avons déjà vu une
partie, il est bon de lister tout ce qui peut être utile.
1.2.1 Théorème de Stokes

Comme nous l'avons vu, si nous avons une courbe fermée Γ, pour toute surface Σ s'appuyant sur la
courbe Γ et pour tout champ vectoriel V :
˛ ¨
V.ds = (∇ ∧ V).dS (1.20)
Γ Σ
où dS est le vecteur normal à la surface élémentaire et proportionnel à la dite surface.
1.2.2 Théorème de Green-Ostrogradsky

Pour tout volume V délimité par une surface Σ, on peut écrire pour un champ vectoriel V :
¨ ˚
V.n dΣ = ∇.V dτ (1.21)
Σ τ
1.2.3 Conséquences
On peut dériver de ces deux théorèmes d'autres relations :
• Formule de Kelvin pour tout courbe fermée Γ et toute surface Σ s'appuyant dessus :
˛ ¨
U ds = − (∇U ) ∧ dΣ (1.22)
Γ Σ
• Formule du gradient pour une surface Σ enfermant le volume τ , pour tout champ scalaire U on a :
¨ ˚
U dS = (∇U )dτ (1.23)
Σ τ
• Formule du rotationnel pour une surface Σ enfermant le volume τ , pour tout champ vectoriel V on
a: ¨ ˚
V ∧ dS = − (∇ ∧ V)dτ (1.24)
Σ τ
1.3. DIFFÉRENTS SYSTÈMES DE COORDONNÉES 15
1.2.4 Égalités vectorielles

Tous ces opérateurs ont certaines propriétés et ainsi des égalités vectorielles sont utiles dans le cours
d'électromagnétisme. Pour tout champ vectoriel A, B, C et tout champ scalaire U :
A∧A=0 (1.25)
A.(A ∧ B) = 0 (1.26)
A ∧ (B ∧ C) = B(A.C) − C(A.B) (1.27)
∇ ∧ (∇U ) = 0 ⇒ si ∇ ∧ A = 0, ∃φ, A = ∇φ (1.28)
∇.(∇ ∧ A) = 0 ⇒ si ∇.A = 0, ∃V, A = ∇ ∧ V (1.29)
1.3 Diérents systèmes de coordonnées

Les diérents systèmes de coordonnées sont tous équivalents dans le sens où la donnée de trois coor-
données est nécessaire pour repérer un point dans l'espace mais les diérents systèmes sont plus ou moins
pratiques pour traiter des problèmes de physique (notamment d'électromagnétisme) en fonction de la
géométrie ou des symétries du système.
1.3.1 Cartésien
Le système de coordonnées cartésien est le plus naturel et le plus simple. Un point de l'espace M est
repéré par les coordonnées (x, y, z ) et où les vecteurs du repère sont xes.
z
M(x,y,z)
ez ey y
ex
x
Figure 1.4 Système de coordonnées cartésien.
∂U ∂U ∂U
∇U = ex + ey + ez (1.30)
∂x ∂y ∂z
∂Vx ∂Vy ∂Vz
∇.V = + + (1.31)
∂x ∂y ∂z

∂Vz ∂Vy ∂Vx ∂Vz ∂Vy ∂Vx
∇∧V = − ex + − ey + − ez (1.32)
∂y ∂z ∂z ∂x ∂x ∂y
∂2U ∂2U ∂2U
∆U = + + (1.33)
∂x2 ∂y 2 ∂z 2
1.3.2 Cylindrique
Le système de coordonnées cylindrique est utile lorsqu'un système a une symétrie de rotation. Un point
de l'espace M est repéré par les coordonnées (ρ, θ, z )
z
M(ρ,θ,z)
ρ ez
eθ
θ eρ
Figure 1.5 Système de coordonnées cylindrique.
∂U 1 ∂U ∂U
∇U = eρ + eθ + ez (1.34)
∂ρ ρ ∂θ ∂z
1 ∂(rVr ) 1 ∂Vθ ∂Vz
∇.V = + + (1.35)
r ∂r r ∂θ ∂z
1 ∂Vz ∂Vθ ∂Vr ∂Vz 1 ∂ ∂Vr
∇∧V = − er + − eθ + (rVθ ) − ez (1.36)
r ∂θ ∂z ∂z ∂r r ∂r ∂θ
1 ∂2U ∂2U ∂2U 1 ∂2U ∂2U

1 ∂ ∂U 1 ∂U
∆U = r + 2 2 + = + + + (1.37)
r ∂r ∂r r ∂θ ∂z 2 ∂r2 r ∂r r2 ∂θ2 ∂z 2
1.3.3 Sphérique
Le système de coordonnées sphérique est utile lorsqu'un système a deux symétrie de rotation. Un point
de l'espace M est repéré par les coordonnées (r, θ, ϕ)
∂U 1 ∂U 1 ∂U
∇U = er + eθ + eϕ (1.38)
∂r r ∂θ r sin θ ∂ϕ
1 ∂(r2 Vr ) 1 ∂ 1 ∂Vϕ
∇.V = 2 + (sin θVθ ) + (1.39)
r ∂r r sin θ ∂θ r sin θ ∂ϕ

1 ∂ ∂Aθ 1 ∂Vr 1 ∂
∇∧V = (sin θVϕ ) − er + − (rAϕ ) eθ (1.40)
r sin θ ∂θ ∂ϕ r sin θ ∂ϕ r ∂r

1 ∂ ∂Vr
+ (rVθ ) − eϕ (1.41)
r ∂r ∂θ
∂2U 2 ∂U 1 ∂2U 1 ∂U 1 ∂2U
∆U = + + + + (1.42)
∂r2 r ∂r r2 ∂ϕ2 r2 tan ϕ ∂ϕ r2 sin2 ϕ ∂θ2
1.4. FONCTION DE GREEN 17
M(r,θ,φ) er
eφ
θ
r
eθ
φ
Figure 1.6 Système de coordonnées sphérique.
1.4 Fonction de Green

Considérons un opérateur diérentiel linéaire 1 Dr,t . On considère alors une équation diérentielle de la
forme :
Dr,t [ψr,t ] = Sr,t (1.43)
où Sr,t est le terme de source. Une solution de cette équation s'obtient par la fonction de Green G(r, r0 , t, t0 )
qui vérie :
Dr,t [G(r, r0 , t, t0 )] = δ(r − r0 )δ(t − t0 ) (1.44)
Alors la fonction réponse vaut :
ˆ ˆ
ψG (r, t) = d3 r0 dt0 G(r, r0 , t, t0 )S(r0 , t0 ) (1.45)
Démonstration :
ˆ ˆ
3 0
Dr,t [ψG (r, t)] = d r dt0 Dr,t [G(r, r0 , t, t0 )]S(r0 , t0 ) (1.46)
ˆ ˆ
= d3 r0 dt0 δ(r − r0 )δ(t − t0 )S(r0 , t0 ) (1.47)
= S(r, t) (1.48)
1. Un opérateur diérentiel une combinaison linéaire des opérateurs qu'on vient de détailler.
Chapitre 2
Électrostatique
2.1 Introduction et historique
L'étude des interactions entre objets mais aussi entre constituants de la matière est un sujet récurrent
en physique. Le problème se pose à diérentes échelles, depuis les corps célestes avec des échelles typiques
gigantesques (taille système solaire 109 km), jusqu'aux constituants de la matière (∼ 10−15 m pour la taille
du noyau). La théorie actuelle de la physique arrive à rendre compte de l'ensemble des phénomènes actuel-
lement interprétés en faisant intervenir seulement quatre interactions, dites fondamentales : l'interaction
gravitationnelle, l'interaction électromagnétique, l'interaction faible et l'interaction forte.
L'électromagnétisme se préoccupe des interactions entre corps chargés ainsi que de la dynamique des
champs créés par ces charges, qu'elles soient immobiles (électrostatique) ou en mouvement (courants,
rayonnements, induction, ondes...).
Coté histoire, l'étude de l'électricité et de l'électromagnétisme est extrêmement ancienne car certains
phénomènes électromagnétiques (la foudre, les aurores boréales ou les propriétés magnétiques de certaines
substances) sont connus depuis la haute antiquité. Cependant leur étude détaillée dut attendre la possibilité
de produire à volonté de l'électricité, c'est-à-dire les machines électrostatiques au 17ème siècle. À partir
de ce moment l'évolution des connaissances se précipite et au début du 18ème siècle on découvre qu'il est
possible de transporter l'électricité par des ls. Mais c'est durant le 18ème siècle que les premières théories
scientiques (telles qu'on les connait) sont développées notamment par Cavendish (Potentiel) et Coulomb
(force électrostatique). Au 19ème, tout s'enchaine avec l'invention de la pile électrique par Volta (1800),
la découverte de l'induction par Faraday (1831) et l'unication de la théorie par Maxwell (1864) conrmé
par Hertz (1885).
La théorie unié de Maxwell est connue sous la forme des équations de Maxwell dans le vide a
qui sont au nombre de 4 :
ρ(r, t)
∇.E(r, t) = Équation de Maxwell-Gauss (2.1)
0
∇.B(r, t) = 0 Équation de Maxwell-Flux (Maxwell-Thomson) (2.2)
∂B(r, t)
∇ ∧ E(r, t) = − Équation de Maxwell-Faraday (2.3)
∂t
1 ∂E(r, t)
∇ ∧ B(r, t) = µ0 j(r, t) + Équation de Maxwell-Ampère (2.4)
c2 ∂t
a. En réalité ces équations restent valable tant qu'on se place pas dans un milieu matériel ou un milieu dilué
chargé comme un plasma.
19
20 CHAPITRE 2. ÉLECTROSTATIQUE
La contribution de Maxwell consiste essentiellement à l'ajout du terme c12 ∂E(r,t)

∂t dans la dernière équa-
tion.
Le couple champ électrique E/champ magnétique B constitue le champ électromagnétique. La densité
de charge ρ et le courant électrique volumique j constitue les sources du champ. La permittivité diélectrique
du vide 0 , la perméabilité magnétique du vide µ0 et la vitesse de la lumière dans le vide c constituent les
constantes des équations avec la relation
c2 0 µ0 = 1 (2.5)
À ces 4 équations on ajoute la loi d'interaction des charges avec le champ électromagnétique, qui
est la force de Lorentz :
F = q(E + v ∧ B) (2.6)
qui est la force exercée par un champ électromagnétique (E, B) sur une charge q munie d'une
vitesse v.
L'ensemble de ces 5 équations combiné avec le principe fondamental de la dynamique de Newton permet
de rendre compte de tous les mouvements des particules chargées et de tous les phénomènes électroma-
gnétiques associés.
Le but de ce cours est de partir des équations de Maxwell et monter qu'elles permettent de prédire tous
les phénomènes électromagnétiques importants observés. Ce cours suppose déjà une connaissance de base
sur l'électromagnétisme de niveau L1/L2 ce qui permettra simplement de faire des rappels sur certaines
points et permettre un approfondissement sur d'autres points importants à la bonne compréhension de
l'électromagnétisme et de sa théorie.
2.2 Cadre de l'électrostatique

Pour commencer nous allons nous placer dans le cas le plus simple de l'électromagnétisme à savoir
l'électrostatique dans le "vide". Cela consiste à considérer le cas statique, ou plutôt stationnaire, à savoir
que toutes les charges sont xes de façon permanente dans l'espace, ou bien si elles se déplacent, elle
se déplacent comme un courant continu dans un circuit (de sorte que ρ et j soient constants). Dans ces
conditions les 5 équations précédentes s'écrivent :
ρ(r)
∇.E(r) = (2.7)
0
∇.B(r) = 0 (2.8)
∇ ∧ E(r) = 0 (2.9)
∇ ∧ B(r) = µ0 j(r) (2.10)
F(r) = q(E(r) + v ∧ B(r)) (2.11)
On voit dans ce cas là que les équations sur le champ électrique et sur le champ magnétique sont totalement
découplés et qu'il est possible d'étudier les deux ainsi que leurs eets séparément. Dans ce premier chapitre
d'électrostatique, nous supposerons qu'il n'existe pas de champ magnétique. Les équations à disposition
sont :
ρ(r)
∇.E(r) =
0
(2.12)
∇ ∧ E(r) = 0
F(r) = qE(r)
2.3. APPROCHE INTÉGRALE ET LOI DE COULOMB 21
2.3 Approche intégrale et loi de Coulomb

2.3.1 Loi de Coulomb (1780)
Coulomb en 1780 propose à partir d'expérience que si deux particules chargées ponctuelles comme
représenté sur la Figure 2.1 alors la force F1/2 qu'exerce la particule 1 sur la particule 2 est :
• Proportionnel à q1 q2 ;
• Attractive ou répulsive en fonction du signe de q1 q2 ;
1
• Proportionnelle à : expérimentalement à l'époque, il n'était pas possible de montrer que
||r1 − r2 ||2
la puissance était -2, il l'a supposé par analogie à la loi de la gravitation.
(M2,q2)
r2-r1
r2 (M1,q1)
r1
Figure 2.1 Deux charges ponctuelles en interaction.

Ainsi on peut écrire que :
q1 q2
F1/2 = K (r2 − r1 ) (2.13)
||r1 − r2 ||3
où le coecient de proportionnalité K vaut :
1
K= (2.14)
4π0
dans le système d'unités actuel. Cette constante est complètement xé. En eet nous avons la relation
0 µ0 c2 = 1 et la vitesse de la lumière est xé à c = 299 792 458 m.s−1 et µ0 également à µ0 = 4π.10−7 H.m−1 .
On a alors une valeur approchée de 0 ' 8, 844.10−12 F.m−1 .
La force électrostatique (loi de Coulomb) s'écrit alors :
1 q1 q2
F1/2 = (r2 − r1 ) (2.15)
4π0 ||r1 − r2 ||3
Avec ces valeurs numériques on peut alors eectuer un ordre de grandeur pour comparer des forces.
Prenons par exemple l'atome d'hydrogène composé d'un proton avec un électron qui gravite autour. On
montre alors que la force électrostatique est 1040 supérieure à la force de gravitation entre les deux corps.
Ainsi l'attraction ou la répulsion entre charges est dominante à l'échelle microscopique et est responsable
des propriétés physiques et chimiques de la matière. Cependant à grand échelle la force de gravitation
reprend le dessus notamment car la "matière macroscopique" est neutre.
Remarques :
• Cette force respecte bien la réciprocité des forces ;
• Dans ce modèle de charges ponctuelles, une charge ne peut agir sur elle-même ;
• L'expression des forces entre atomes/électrons dans la matière est bien plus compliqué notamment
à cause des eets d'écrantage qui apparaissent.
Si au lieu d'avoir 2 charges, nous avons un ensemble D de N charges ponctuelles autour d'une autre
charge Q alors la force exerce cet ensemble sur la charge q est égale à la sommes des forces qu'exercent
ces charges séparément :
N
X 1 qqi
FD/q = (r − ri ) (2.16)
4π0 ||r − ri ||3
i=1
où r est la position de la charge q et ri de la ième charge de l'ensemble D de charge qi . On parle alors de
superposition ou de linéarité des forces électrostatiques. Ce principe de superposition provient directement
de la linéarité des équations de Maxwell.
2.3.2 Champ de Coulomb

On a vu qu'on ne pouvait pas avoir de force d'une charge sur elle-même (car Coulomb est vrai que si
la charge est strictement poncuelle). On peut alors par commodité introduire le champ électrique d'une
charge comme la force par unité de charge en un point de l'espace. Si on reprend la force précédente avec
les deux charges alors on a :
F1/2 1 q1
E1 (r) = = r r = r2 − r1 (2.17)
q2 4π0 ||r||3
et si on prend la distribution de charges :
N
FD/q X 1 qi
ED (r) = = (r − ri ) (2.18)
q 4π0 ||r − ri ||3
i=1
On retrouve ainsi que la force de Lorentz F = qE avec cette dénition du champ électrique.
2.3.3 Répartition continue de charges

Si au lieu d'avoir une répartition ponctuelle on suppose qu'on a une répartition continue dans un domaine
D. En chaque point du domaine on peut dénir une densité de charge volumique ρ(r) qui est moyenné sur
un volume mésoscopique 1 . Dans ce cas là, le champ de Coulomb s'écrit en un point M (r) quelconque de
l'espace : ˆ
1 ρ(r0 )
ED (r) = d3 r0 (r − r0 ) (2.19)
4π0 D ||r − r0 ||3
Dans la suite on utilisera toujours l'expression continue car on peut passer à l'expression ponctuelle
simplement en écrivant :
XN
ρ(r) = qi δ(r − ri ) (2.20)
i=1
On peut écrire les mêmes expressions pour des densités surfaciques ou linéiques.
Si on veut calculer quelque chose avec ces intégrales, il faut les écrire en détail. Par exemple si on se
place en cartésien et qu'on calcule la composante x du champ électrique, on a :
ˆ
(x1 − x2 )ρ(x2 , y2 , z2 )dx2 dy2 dz2
Ex (x1 , y1 , z1 ) = 2 2 2 3/2
(2.21)
R3 4π0 [(x1 − x2 ) + (y1 − y2 ) + (z1 − z2 ) ]
1. On peut dénir la taille de ce volume mésoscopique comme en mécanique des uides. Si on trace l'allure de la densité de
charges en fonction du volume sur lequel on moyenne alors : pour des petits volumes, les uctuations seront très importantes
car on verra chaque charge au fur et à mesure que le volume grandit ; puis la densité va atteindre un plateau sur lequel on
se placera pour dénir la taille du volume mésoscopique ; enn pour des tailles trop grandes, on verra des variations lentes
correspondant aux éventuelles variations de densité de charges macroscopique qu'on cherche à quantier avec la fonction
ρ(r).
2.4. LIENS AVEC LES ÉQUATIONS DE MAXWELL 23
Nous n'allons pas beaucoup utiliser cette formule. Nous l'écrivons ici que pour souligner le fait que nous
avons complètement résolu tous les problèmes d'électrostatique dans lesquels on connaît la position de
toutes les charges (ou la répartition). On se donne la répartition de charges, quel est le champ en tout
point de l'espace ? La réponse est simplement : calculez cette intégrale ! Il n'y a plus rien à dire ce le sujet
car il ne s'agit plus qu'une question de calcul d'intégrale faisable par un ordinateur. Il est important d'avoir
à l'esprit que tout ce qu'on pourra faire plus tard dans ce chapitre est déjà entièrement contenu dans ce
que nous venons de faire. Cependant est généralement plus simple de deviner la valeur de l'intégrale plutôt
que de la calculer mais pour deviner il faut des nouveaux outils. C'est ce que nous allons essayer de faire.
2.4 Liens avec les équations de Maxwell

2.4.1 Compatibilité avec la force de Coulomb
Quel est le lien entre la force de Coulomb décrite précédemment et les équations de Maxwell (version élec-
trostatique) de l'introduction ? On voit dans les deux cas on relie le champ électrique créé à la distribution
de charge. Essayons de voir si l'un n'est pas juste la réécriture de l'autre. Il faut tout d'abord remarquer
une petite égalité vectorielle (qu'on peut facilement retrouver avec le développement de Taylor) :
(r − r0 )

1 1
= −∇r = ∇r0 (2.22)
||r − r0 ||3 ||r − r0 || ||r − r0 ||
On a alors :
ˆ
1 ρ(r0 )
ED (r) = d3 r0 (r − r0 ) (2.23)
4π0 D ||r − r0 ||3
ˆ
1 3 0 0 1
= d r ρ(r ) −∇r (2.24)
4π0 D ||r − r0 ||
ˆ
ρ(r0 )

1
= −∇r d3 r0 (2.25)
4π0 D ||r − r0 ||
Or comme ∇r ∧ (∇r f ) = 0 pour tout f cela signie que :
∇r ∧ ED (r) = 0 (2.26)
Le champ de Coulomb vérie déjà une équation de Maxwell de l'électrostatique. Pour la seconde équation,
prenons la divergence de l'expression du champ de Coulomb 2 :
ˆ
1 3 0 0 1
∇r .ED (r) = ∇r . − d r ρ(r )∇r
4π0 D ||r − r0 ||
ˆ
1 3 0 0 2 1
= − d r ρ(r )∇r
4π0 D ||r − r0 ||
ˆ
1
= d3 r0 ρ(r0 )4πδ(r − r0 )
4π0 D
ρ(r)
=
0
Et donc ED (r) vérie bien la seconde équation de l'électrostatique.
Tout ce qu'on vient faire ne dépend pas de la nature du champ mais simplement de sa forme en 1/r2 .
Ainsi tous les champs dits newtoniens vérient ce type d'équation locale type Maxwell.

2. En utilisant le résultat suivant sur le laplacien de 1/r : −4πδ(r − r0 ) = ∇2r 1
||r−r0 ||
2.4.2 Équivalence des deux approches

Il est possible de prouver l'équivalence des approches intégrale et local grâce à un théorème.
Théorème d'Helmoltz-Hodge :
Si V(r) est un champ de vecteur qui est nul à l'inni alors V(r) est déni par la donnée de sa
divergence et de son rotationnel. La solution est unique.
Démonstration :
Si lim V(r) = 0
r→+∞
alors il existe V(k) tel que :
ˆ ˆ
1 −ik.r 1
V(k) = 3
dr V(r)e V(r) = dk3 V(k)eik.r (2.27)
(2π)3/2 R3 (2π)3/2 R3
Si on a :
∇.V = a(r) (2.28)

∇ ∧ V = b(r) (2.29)
avec a et b connus en tout point de l'espace. Si on prend la transformée de Fourier de ces deux équations
alors :
ik.V = A(k) (2.30)

ik ∧ V = B(k) (2.31)
La première équation nous renseigne sur V dans la direction de k (1 équation scalaire) et la seconde équation
nous renseigne sur la direction perpendiculaire à k (2 équations scalaire seulement pour une direction dans
un plan). Ainsi nous avons 3 équations pour 3 inconnues (les 3 directions de V ) donc on a une solution
unique si elle existe. On peut alors par transformée de Fourier inverse trouver une unique solution pour
V si elle existe.
Ainsi par ce théorème montre que les approches local et intégrale sont parfaitement équivalentes en
termes mathématiques et donc physiquement aussi.
2.5 Théorème de Gauss

Nous avons dit que les deux approches étaient équivalentes et que le calcul des champs de Coulomb
consistait en un calcul compliqué d'intégrales. Est-ce que l'approche locale des équations de Maxwell
peuvent permettre de faciliter le calcul ? La réponse est oui mais pas toujours. Si on reprend l'équation de
Maxwell-Gauss pour une distribution de charge ρ(r) dans un certain domaine D de l'espace. On intègre
alors sur ce domaine D de frontière Σ :
ρ(r)
∇.E =

ˆ ˆ 0
ρ(r)
∇.Edτ = dτ
0
Dˆ D ˆ
1
E.dS = ρ(r)
Σ 0 D
2.6. STRUCTURE DU CHAMP ÉLECTROSTATIQUE 25
Ainsi le théorème de Gauss s'écrit :

QD
φΣ (E) = (2.32)
0
où φΣ (E) est le ux de E à travers la surface Σ et QD la charge totale contenue dans le domaine
D. Dit autrement, le ux du champ électrique à travers une surface (dit surface de Gauss) est
égal à la charge contenue dans le volume délimité par la surface divisé par la permittivité du
vide.
Ce théorème est très utile pour calculer des champs électriques mais il nécessite que le problème possède
au moins deux symétries dans le problème pour pouvoir calculer le ux de E sans ambiguïté. Il nous faut
donc voir maintenant en quoi les symétries du problème permettent de tels calculs.
2.6 Structure du champ électrostatique

2.6.1 Invariance et principe de Curie
Les deux équations de Maxwell de l'électrostatique
ρ(r)
∇.E(r) =
0
∇ ∧ E(r) = 0
ne jouent pas le même rôle. En eet la première relie le champ électrique à sa source, à savoir une
distribution de charges dans l'espace. À contrario, la seconde équation ne possède que le champ électrique
et ainsi nous renseigne uniquement sur la structure de celui-ci.
Principe de Curie :
Les invariants d'une distribution de sources laissent invariant la distribution des eets. Les mêmes
causes produisent les mêmes eets et toute transformation laissant invariant les causes laisse
invariant les eets.
En pratique pour l'électrostatique les causes sont les distributions de charges et toute modication
continue qui laisse invariant ρ laisse invariant le champ électrique E. Ceci permet ainsi de réduire le
nombre de variables dont dépend le champ, uniquement grâce à l'équation de Maxwell-Gauss.
Exemples :
• La translation le long d'un l inni laisse invariant le problème donc si on prend l'axe x le long du
l, il n'y aura pas de dépendance en x.
• La rotation autour d'un l autour de son axe permet de dire qu'il y a une symétrie axiale et donc
une indépendance avec l'angle de rotation.
Remarque : Dans le cas général le principe de Curie n'est vrai que en moyenne. En eet on peut prendre
deux exemples où il faut moyenner pour ne pas mettre en défaut le principe de Curie :
• Lorsqu'une bille se situe en haut d'une barrière de potentiel à une dimension, elle se trouve dans un
état métastable. Elle peut tomber soit à gauche ou soit à droite quand on la lâche en haut.
• Lorsqu'une tige rigide d'une certaine longueur est xée dans un mur, si on appuie sur la tige dans
son axe alors la tige va amber. Elle peut amber dans n'importe qu'elle direction de l'espace, c'est
à dire sous n'importe quel angle entre 0 et 2π .
Dans ce deux cas on a une symétrie du système et si on fait une expérience le système va "choisir" une
direction et donc briser cette symétrie. Cependant si on répète un grand nombre de fois l'expérience et
qu'on trace la distribution de probabilité alors celle-ci va posséder les mêmes symétries que le système. Le
principe de Curie est alors vérié.
Une autre manière de le dire est qu'il faut retrouver le principe de Curie sur l'ensemble des solutions
et non sur une solution particulière.
2.6.2 Exemples d'invariants

En électromagnétisme il est possible de rencontrer plusieurs types de symétries. Voici des exemples :
• Invariant par translation : si le système est invariant par translation selon l'axe x alors le champ
électrostatique ne dépendra pas de la coordonnée x.
• Invariant par rotation : si le système est invariant par rotation autour d'un certain axe dont la rotation
est représenté par la coordonnée θ alors le champ électrostatique ne dépend pas de la coordonnée θ.
• Invariant cylindrique : si le système est invariant par rotation autour de l'axe z ainsi que par transla-
tion selon l'axe z alors le champ électrostatique ne dépendra pas de ces deux coordonnées c'est-à-dire
qu'il ne dépend que de la distance à l'axe.
• Invariant sphérique : si le système est invariant selon deux rotations de coordonnées θ et φ c'est-à-dire
invariant par rotation autour d'un certain point O alors le champ électrique ne dépend pas de ces
deux coordonnées c'est-à-dire qu'il ne dépend que de la distance r au point O.
2.6.3 Propriétés de symétrie des champs électriques

Soit D une distribution de charges. Soit Π+ un plan de symétrie de la distribution de charges. Ainsi
pour tout M ∈ D alors si M 0 est le symétrique du point M par rapport au plan Π+ qu'on notera
M 0 = symΠ+ (M ) alors on a nécessairement ρ(M ) = ρ(M 0 ). De même si Π− est un plan d'antisymétrie de
la distribution de charge alors si M 0 = symΠ− (M ) alors on a ρ(M ) = −ρ(M 0 ).
Si on se place dans le cas où la distribution possède une symétrie Π+ . Soit M et M 0 deux points
symétriques l'un par rapport à l'autre par Π+ et soit P un point de D. Calculons alors le champ créé au
point M par la distribution D.
ˆ
1 3 ρ(P )
E(M ) = dEP (M ) avec dEP (M ) = d r PM (2.33)
D 4π0 |PM|3
On peut alors calculer le champ créé au point M par les points P et P 0 symétriques par rapport au plan
Π+ . Par superposition on a :
ρ(P 0 ) 0

1 3 ρ(P )
dE P,P 0 (M ) = d r PM + 0 3 P M (2.34)
4π0 |PM|3 |P M|
On peut ensuite calculer le champ créé au point M 0 par les mêmes points P et P 0 . On exprime alors PM0
et P0 M0 comme les symétriques de P0 M et PM. On remarque également par symétrie que les distances
|PM0 | et |P0 M| sont égales ainsi que les distances |P0 M0 | et |PM|. Enn par hypothèse on a ρ(P 0 ) = ρ(P ).
ρ(P 0 ) 0 0

0 1 3 ρ(P ) 0
dEP,P 0 (M ) = d r PM + 0 0 3 P M (2.35)
4π0 |PM0 |3 |P M |
ρ(P 0 )

1 3 ρ(P ) 0
= d r symΠ+ (P M) + symΠ+ (PM) (2.36)
4π0 |P0 M|3 |PM|3
ρ(P 0 )

1 3 0 ρ(P )
= d r symΠ+ (P M) + symΠ+ (PM) (2.37)
4π0 |P0 M|3 |PM|3
= symΠ+ (dEP,P 0 (M )) (2.38)
Ainsi chaque champ élémentaire créé au point M 0 est symétrique par rapport au plan Π+ au champ
élémentaire créé au point M . Et par linéarité de l'intégrale on a :
E(M ) = symΠ+ (E(M 0 )) (2.39)
La conséquence de cette égalité est que si M ∈ Π+ alors on applique le résultat précédent avec M = M 0
et donc on a E(M ) = symΠ+ (E(M )). Ceci a pour conséquence que dans ce cas on a E⊥ (M ) = 0.
Pour une distribution de charges donnée, le champ électrique est contenu dans les plans de
symétrie. Pour xer la direction, il est alors nécessaire que le système possède deux plans de
symétrie qui s'intersectent.
De la même manière on peut montrer que si on a un plan d'antisymétrie Π− de la distribution de charge

et M ∈ Π− alors on trouve que E(M ) ⊥ Π− .
Pour une distribution de charges donnée, le champ électrique est perpendiculaire aux plans
d'antisymétrie. Pour xer la direction, il ne faut qu'un plan d'antisymétrie.
Ainsi les symétries du problème peut permettre de xer la direction du champ E (ou d'éliminer une
composante) en certains points de l'espace.
2.6.4 Calculs de champ électrostatique

• Champ de Coulomb :
Soit une charge ponctuelle q placé en O. On utilise les coordonnées sphériques pour se repérer dans
l'espace (O, r, θ, φ). Considérons un point M de l'espace.
La première chose à faire lorsqu'on aborde un problème d'électrostatique est d'étudier séparément
les symétries et les invariances du problème.
Concernant les invariances : le système est un invariant sphérique et ainsi le champ électrique re-
cherché ne dépend que de la distance à la charge : E(r, θ, φ) = E(r).
Concernant les symétries : Tout plan passant par O et contenant M est un plan de symétrie, ainsi
le champ électrique recherché doit être contenu dans ces deux plans, ainsi le champ électrique est
radial. Donc au nal :
E = E(r)er (2.40)
Appliquons à présent le théorème de Gauss. On choisit une surface, dit de Gauss, Σ à travers
laquelle le ux du champ électrique est facilement calculable, c'est-à-dire une surface parallèle ou
perpendiculaire au champ électrique en tout point. Ici la surface de Gauss naturelle est une sphère
de rayon r. Dans ce cas, en tout point le champ est perpendiculaire à la surface de Gauss et le ux
élémentaire vaut simplement E(r)dS . On pourra alors sortir E(r) de l'intégrale qui sera constant
sur la sphère. La charge contenue dans cette sphère est q quelque soit le rayon de celle-ci. Donc
l'application du théorème donne :
QD
φΣ (E) = (2.41)
0
ˆ
q
E.dS = (2.42)
Σˆ 0
q
E(r) dS = (2.43)
Σ 0
q
E(r)4πr2 = (2.44)
0
1 q
E(r) = (2.45)
4π0 r2
Et donc :
1 q
E(r) = er (2.46)
4π0 r2
On retrouve bien le champ de Coulomb sans trop de surprise.
• Fil inni :
On considère un l inni, c'est-à-dire un cylindre inniment long selon la direction z et donc on néglige
le rayon. Ce solide est uniformément chargé et on note λ la densité linéique. Dans ce problème, la
distribution de charges est invariante par translation suivant l'axe z et invariant par rotation autour
de l'axe z . On a une symétrie cylindrique donc E(ρ, θ, z) = E(ρ). Soit un point M de l'espace. Le
plan contenant l'axe de symétrie du l (c'est-à-dire le l) ainsi que le point M est un axe de symétrie
des charges, de même que le plan perpendiculaire au l et contenant le point M . Ainsi le champ
appartient à l'intersection de ces deux plans et est selon eρ au point M . Donc on a :
E = E(ρ)eρ (2.47)
La surface de Gauss Σ naturelle pour ce problème est un cylindre de hauteur h et de rayon r. Le

champ électrique est alors parallèle sur les bords supérieur et inférieur ce qui signie que le ux du
champ électrique à travers ces surfaces est nul. Le champ est perpendiculaire sur la surface latérale
et donc le calcul est le même que même que précédemment avec la surface qui vaut 2πρh cette fois :
QD
φΣ (E) = (2.48)
0
ˆ
λh
E(ρ) dS = (2.49)
Σ 0
λh
E(ρ)2πρh = (2.50)
0
λ
E(ρ) = (2.51)
2πρ0
Et nalement on trouve :
λ
E(r) = eρ (2.52)
2πρ0
• Plan uniformément chargé :
On considère un plan inni selon l'axe x et y et chargé uniformément avec une densité surfacique σ . La
distribution de charges est invariante par translation selon les axes x et y . Ainsi on a E(x, y, z) = E(z).
Soit un point M de l'espace. Tout plan orthogonale au plan chargé et contenant le point M est plan
de symétrie des charges ainsi on a un champ selon ez . Donc :
E = E(z)ez (2.53)
De plus le plan est lui même plan de symétrie des charges et donc on a que E(z) = −E(−z).
A présent il nous est possible d'appliquer le théorème de Gauss sur des boites bien choisies. Tout
d'abord si on prend n'importe qu'elle boite qui se trouve complètement dans la partie z > 0 ou z < 0.
Par exemple pour les z > 0, on prend une boite de surface S et d'ordonnées z1 et z2 > z1 . Le champ
est parallèle aux cotés de la boite donc le ux est nulle. Sur les faces supérieure et inférieure le ux
total vaut E(z2 )S − E(z1 )S . De plus il n'y a pas de charge à l'intérieur de la boite donc le théorème
de Gauss permet d'écrire E(z2 )S − E(z1 )S = 0 donc :
∀z1 , z2 > 0 E(z1 ) = E(z2 ) (2.54)

Il est possible de montrer la même chose pour les z < 0.

Nous avons montré que le champ est constant dans les deux parties de l'espace délité par le plan
chargé, on peut à présent calculer le champ. On prend alors une boite, toujours de surface S dont la
surface supérieure est en z et la surface inférieure est en −z . Le ux à travers cette boite est alors
E(z)S − E(−z)S = 2SE(z). La charge totale dans la boite vaut alors σS et donc le théorème de
Gauss donne :
σS
2SE(z) = (2.55)
0
Et donc :
σ
E(z) = ez , z>0 (2.56)
20
σ
E(z) = − ez , z<0 (2.57)
20
On remarque alors qu'on retrouve la relation de passage d'un plan chargé uniformément en surface :
σ
E(0+ ) − E(0− ) = ez (2.58)
0
• Disque uniformément chargé
On considère un disque uniformément chargé en surface de densité surfacique σ mais de taille ni, de
rayon R comme sur la Figure 2.2. Dans ce cas là, l'application du théorème de Gauss n'est
pas possible. En eet le système ne possède pas assez de symétries et d'invariances, il en faudrait
au moins 2 symétries pour xer la direction du champ partout et deux invariances pour que le champ
ne dépende que d'une coordonnée. Ici le système n'a qu'une seule invariance et symétrie de charge.
Dans ce cas là, le seul moyen de pouvoir calculer le champ est de repasser par le calcul intégral.
dS
P u
R r
α M
O θ
dE
z
Figure 2.2 Disque uniformément chargé.

Calculer le champ dans tout l'espace n'est faisable (facilement) que dans des cas très particuliers.
Ici, le champ sur l'axe est facilement calculable. La calcul consiste à calculer le champ innitésimal
dE produit par une surface innitésimale dS situé au point P sur le disque. On a :
1 σdS
dE = u = dE u (2.59)
4π0 r2
Cependant si on se place en un point M sur l'axe de symétrie z du disque alors tout plan contenant
l'axe est plan de symétrie des charges et donc le champ sur l'axe z est selon l'axe z : E = Eez . Ainsi
on va utiliser le principe de superposition en sommant les champs innitésimaux mais comme on sait
que le champ ne sera que selon l'axe z (car en réalité quand on va sommer les autres composantes
vont s'annuler) on va projeter tout d'abord sur l'axe avant le sommer.
1 σdS cos α
dEz = (2.60)
4π0 r2
Il reste maintenant à intégrer et les astuces de calculs sont alors les bienvenues dans ces cas là.
Tout d'abord, plaçons nous dans les coordonnées cylindriques, ainsi l'élément de surface s'écrit (à
justier) : dS = ρdρdθ = rdrdθ (car on a r2 = ρ2 + z 2 et on se place à z xé pour l'intégration, d'où
l'égalité des diérentielles) et donc l'intégration devient :
ˆ 2π ˆ R ˆ R
σ cos α σ cos α
Ez = dθ dr = dr (2.61)
4π0 0 0 r 20 0 r
Il s'agit à présent que d'un problème d'intégration. Eectuons un changement de variable adéquate
2
en sachant qu'on a cos α = zr donc −dα sin α = − zdr
r2
. Donc on a dr = dα sin α rz et zr = cos1 α . De plus
on appelle α0 l'angle sous lequel le point M voit le disque et donc il s'agit de la borne d'intégration
du changement de variable :
ˆ α0
σ r2 cos α
Ez = dα sin α (2.62)
20 0 z r
ˆ α0
σ r
= dα sin α cos α (2.63)
20 0 z
ˆ α0
σ 1
= dα sin α cos α (2.64)
20 0 cos α
ˆ α0
σ
= dα sin α (2.65)
20 0
σ
= [− cos α]α0 0 (2.66)
20
σ
= (1 − cos α0 ) (2.67)
20
Et donc on trouve le champ sur l'axe du disque :

σ σ z
E= (1 − cos α0 )ez = 1− √ ez (2.68)
20 20 z + R2
2
On peut vérier que notre résultat est cohérent avec ce qu'on connait déjà 3 . Si on fait tendre la taille
du disque à l'inni ce qui revient à faire tendre R → ∞ ou l'angle α0 → π2 alors on trouve :
σ
ER→∞ = ez (2.69)
20
Ce qui est en accord avec le résultat trouvé plus haut. Cependant il n'est pas simple de connaitre le
champ partout dans l'espace pour un disque de taille nie.
2.7 Potentiel électrostatique

Nous avons vu que les calculs d'intégrales étaient parfois compliqués et même dans les cas plus faciles
comme le calcul sur l'axe de symétrie du disque chargé. Il est possible d'utiliser un outil qui servira
d'intermédiaire de calcul mais à qui on va donner également une signication physique : le potentiel
électrostatique.
2.7.1 Dénition
L'équation de structure de Maxwell-Faraday de l'électrostatique s'écrit
∇∧E=0 (2.70)
3. C'est un réexe à toujours avoir en physique !
2.7. POTENTIEL ÉLECTROSTATIQUE 31
Dans ce cas là a qu'il existe un champ scalaire V , appelé potentiel électrostatique, tel que :
E(r) = −∇V (r) (2.71)
a. Nous n'avons pas eectué la démonstration mais il est nécessaire que le champ électrique s'annule à l'inni
mais dans toutes les situations physiques rencontrées c'est le cas.
Le signe "-" est pur convention et permet de faire l'analogie avec les autres cas de diusion. Ce potentiel
peut être vu comme un intermédiaire de calcul, même si nous verrons sa signication physique plus tard.
En eet, un champ scalaire est toujours plus facile à calculer qu'un champ vectoriel (3 fois moins de
composantes) et le passage en prenant le gradient n'est rarement un problème. Cette facilité de calcul est
d'autant plus vrai lorsque les symétries du système sont complexes.
2.7.2 Exemples de calculs

Reprenons l'exemple de la charge ponctuelle qui est le système le plus simple en électrostatique. Le
champ électrique a été vu plus tôt :
1 q r
E(r) = (2.72)
4π0 ||r||2 ||r||
et le potentiel électrostatique correspondant est simplement :
q
V (r) = + cte (2.73)
4π0 r
Il reste une constante à déterminer dû à l'intégration du champ électrique. Cette constante peut prendre
n'importe quelle valeur de manière arbitraire sans modier le champ électrique. Cependant un choix
naturel, dit du potentiel absolu, consiste de prendre cette constante nulle c'est à dire de prendre le potentiel
nulle à l'inni là où il n'y a pas de charge électrique.
Si on généralise, on peut écrire le potentiel pour une distribution quelconque de N charges ponctuelles
qi :
N
X qi
VN (r) = + cte (2.74)
4π0 ||r − ri ||
i=1
Pour la distribution continue de charges D, le travail a déjà été fait précédemment dans l'équation (2.25)
en écrivant le champ électrique sous forme d'un gradient d'une expression et donc dans ce cas :
ˆ
1 ρ(r0 )
VD (r) = dr0 (2.75)
4π0 D ||r − r0 ||
2.7.3 Signication physique

Le potentiel électrostatique a une signication physique assez directe. En eet, on prend deux points
de l'espace M1 et M2 ainsi qu'un chemin quelconque Γ entre les deux. Calculons le travail de la force
électrostatique qE quand on déplace la charge q dans un champ électrique E de M1 à M2 :
ˆ M2
W12 = qE.dl
M1
ˆ M2
= −q (∇V ).dl
M1
Et donc :
W12 = q[V (M1 ) − V (M2 )] (2.76)
Le potentiel est directement lié à l'énergie potentielle électrostatique. On peut également remarquer que le
résultat nal ne dépend que du point de départ et d'arrivée du trajet de la charge q et donc que le chemin
suivi Γ n'a pas d'importance pour le calcul. On retrouve la propriété que le champ électrostatique est un
champ à circulation conservative.
2.7.4 Développement multipolaire

Un des intérêts du potentiel est que le théorème de superposition fonctionne encore sans qu'on ait à se
préoccuper de la direction ou du sens du vecteur qu'on somme puisque le potentiel est un champ scalaire.
Si on prend par exemple un dipôle comme sur la Figure 2.3, c'est-à-dire deux charges q et −q séparées
M1 O M2
q -q
r1 r r2
Figure 2.3 Dipôle électrostatique.

par une distance a et qu'on regarde le champ loin de ce dipôle |OM| = r a alors on peut écrire :

q 1 1 1 p.r
Vdipole (r) = V (q, r1 ) + V (−q, r2 ) = − ' p = 2qaex (2.77)
4π0 |r1 − r| |r2 − r| 4π0 r3
où l'expression nale est obtenue de manière classique par un développement limité. Il est également
possible de généraliser le dipôle électrostatique. En eet, si on prend une distribution D de charges dans
une région de l'espace et qu'on regarde cette distribution "de loin". Pour être plus précis, si O est le
barycentre de la distribution de charges et a sa taille caractéristique alors on se place pour des distances
|OM| = r a. On peut écrire :
ˆ
1 ρ(r0 )
VD (r) = dr0
4π0 D ||r − r0 ||
ˆ
1 ρ(r0 )
= dr0
4π0 D [(r − r0 )(r − r0 )]1/2
ˆ
1
= dr0 ρ(r0 )(r2 + r02 − 2r.r0 )−1/2
4π0 D
ˆ 0
0 2 !−1/2
1 0 ρ(r ) 2r.r0 r
= dr 1− 2 +O
4π0 D r r r
ˆ 0
0 2 !
1 0 ρ(r ) r.r0 r
' dr 1+ 2 +O
4π0 D r r r
ˆ 0 2
QD 1 0 0 0 r
' + 3
dr ρ(r )r.r + O
4π0 r 4π0 r D r
ˆ 0 2
QD r. r
' + 3
dr0 ρ(r0 )r0 + O
4π0 r 4π0 r D r
On trouve donc : ˆ
QD P.r
VD (r) = + + reste P= dr0 ρ(r0 )r0 (2.78)
4π0 r 4π0 r3 D
2.8. ÉQUATION DE POISSON 33
On obtient alors les deux premiers termes du développement multipolaire. Le premier terme représente
le champ créé par un monopole, c'est à dire qu'à l'ordre 0 le champ vu au loin est juste le champ d'une
particule ponctuelle de charge QD qui est la charge totale de la distribution D. Pour le dipôle précédent
(q ,−q ) la charge totale étant nulle, ce terme n'apparaît pas. Le second terme représente le champ d'un
dipôle électrostatique caractérisé par la grandeur P . Le développement pourrait se poursuivre pour avoir
les autres termes et notamment le terme suivant "quadripolaire". Tout ceci n'est possible que si on connait
la distribution D.
2.8 Équation de Poisson

2.8.1 Mise en équation
Si à présent on prend l'équation de Maxwell-Gauss et qu'on remplace E = −∇V on trouve :
ρ
∇.E = ∇.(−∇V ) = −∇2 V = (2.79)
0
L'équation de Poisson s'écrit alors :

ρ
∆V + =0 (2.80)
0
Cette équation permet de calculer V si on connait la densité de charge, donc cela n'a pas l'air de nous
apporter plus d'information car on peut le faire avec le champ de Coulomb. En réalité cette équation a
son utilité, tout d'abord dans certains cas particuliers. Par exemple si on a une situation physique où on
impose un potentiel ou une diérence de potentiel (comme un condensateur). Dans ce cas là on ne connait
pas la densité de charge mais on a des conditions aux limites pour le potentiel et une équation diérentielle
linéaire.
2.8.2 Application à l'écrantage de Debye [Exercice]

Les plasmas sont des milieux macroscopiquement neutre, partiellement ou totalement ionisés. Naturels
ou articiels, on les rencontre sous de nombreuses formes : gaz lors du passage d'un arc électrique, gaz
dans les étoiles, ionosphère...
On considère ici, un plasma d'argon contenant en moyenne et par unité de volume, ne électrons libres
de masse me et de charge −e, ni = ne ions Ar+ de masse mi , et no atomes Ar de masse m0 . On dénit le
degré d'ionisation de ce plasma par
ne
α= (2.81)
ne + n0
Le plasma est supposé en équilibre thermodynamique local et on note T sa température.
Soit un ion argon Ar+ en particulier, placé en O, et pris comme origine. Du fait de l'attraction cou-
lombienne, on observe au voisinage de cet ion un surplus de charge négative, responsable d'un écart local
d'électroneutralité globale du plasma. Soit V (r) le potentiel qui règne en un point M situé à la distance r
de l'ion Ar+ (l'origine des potentiels est prise à l'inni). Les densités volumiques d'ions et d'électrons en
M sont données par la loi de Boltzmann :
eV (r) eV (r)
−
n+ = ne e kB T
et n− = ne e kB T (2.82)
1. Donner l'expression de la densité volumique totale de charges ρ(r) pour r 6= 0.

2. Déterminer l'équation diérentielle satisfaite par V (r).
3. On se place dorénavant dans l'hypothèse eV (r) kB T .
(a) Simplier l'équation obtenue à la question précédente et trouver une équation diérentielle en
u(r) = rV (r).
(b) Résoudre l'équation diérentielle en introduisant deux constantes d'intégration A et B .
4. On admet que V (∞) = 0 et qu'au voisinage de l'ion Ar+ , l'inuence de sa charge, supposée ponc-
tuelle, l'emporte sur celle des charges électroniques distribuées en volume
(a) Déterminer les deux constantes d'intégration A et B .
(b) Donne l'expression du potentiel V (r) en fonction de e, 0 , r et d'une distance caractéristique
λD appelé longueur de Debye, à expliciter en fonction de 0 , kB , T , ne et e.
(c) Commenter le résultat obtenu et interpréter λD .
5. Pour ce plasma d'argon, ne = 3.0 1021 m−3 . Calculer la valeur numérique de λD pour des tempéra-
tures de 1000 et 10000 K.
Correction :
1. La densité volumique de charges se calcule comme
X
ρ(r) = nk (r)qk (2.83)
k
où nk et qk désignent respectivement la densité volumique de porteurs de charge de type k t la charge

individuelle de ce porteur. Ici :
ρ(r) = (ni (r) − ne (r))e (2.84)
Á la température du plasma, les densités de porteurs suivent la loi donnée dans l'énoncé et donc :

−
eV (r) eV (r) eV (r)
ρ(r) = ne e e kB T
−e kB T
⇒ ρ(r) = −2ne esh (2.85)
kB T
2. Le potentiel électrostatique V (r) est relié à ρ(r) par l'équation de Poisson
ρ
∆V + =0 (2.86)
0
qui compte tenu de l'expression du laplacien :
∂ 2 V (r) 2 ∂V (r) 1 ∂ 2 rV (r)

2ne e eV (r)
2
+ = 2
= sh (2.87)
∂r r ∂r r ∂r 0 kB T
3. Cas eV (r) kB T
(a) Dans le cas où l'agitation thermique l'emporte à haute température devant l'énergie d'interaction
électrostatique (tel que eV (r) kB T ), on peut linéariser la fonction sh(x) = x+o(x), et donc :
∂ 2 rV (r) 2ne re eV (r)
2
= (2.88)
∂r 0 kB T
En posant u(r) = rV (r), on trouve alors :
∂ 2 u(r) ne e 2
− 2 u(r) = 0 (2.89)
∂r2 0 kB T
En raisonnant sur l'équation, il est clair que 2n
0 kB T
e 2e
2 a la dimension d'une longueur au carré, on
dénit donc : r
0 kB T
λD = (2.90)
2ne e2
Et donc on obtient :
∂ 2 u(r) u(r)
− 2 =0 (2.91)
∂r2 λD
2.8. ÉQUATION DE POISSON 35
(b) La solution générale de cette équation est :

− λr r
u(r) = Ae D + Be λD (2.92)
où A et B sont des constante d'intégration. On aurait aussi pu l'écrire sous la forme d'une
somme de cosinus et sinus hyperbolique. Dans le cas d'un potentiel qui s'annule à l'inni, la
première forme est plus pratique.
4. Solution :
(a) Comme le potentiel doit rester nul à l'inni il faut que u(r) soit nul à l'inni et donc B = 0.
Au voisinage immédiat de l'ion Ar+ situé à l'origine r = 0, le potentiel électrostatique tend vers
le potentiel de Coulomb créé par une charge ponctuelle de charge e soit :
e
r∼ (2.93)
4π0 r
Donc nécessairement : e
A= (2.94)
4π0
(b) Compte tenu des résultats précédents :
r
e − r 0 kB T
V (r) = e λD avec λD = (2.95)
4π0 r 2ne e2
Le potentiel électrostatique créé par l'atome d'argon dans le plasma décroît plus rapidement
(c)
que le potentiel Colombien,V (r) = 4πe 0 créé par un ion Ar+ seul en r = 0. En notant VC le
r
potentiel colombien et VD le potentiel de Debye, le rapport est VD /VC = e− λD . Il s'atténue sur
la distance caractéristique λD qui est d'autant plus grande que la température est élevée. Pour
r λD , on a VD /VC → 0, ce qui signie que les charges distantes de r λD ne ressentent
plus l'eet électrostatique de la charge placé en r = 0. On dit que la charge centrale est écrantée.
Ce phénomène d'écrantage a été expliqué par le physicien américain Peter Debye.
5. Avec les valeurs données on trouve λD = 2.8 10−8 m pour T = 1000 Ket λD = 8.9 10−8 m pour T
= 10000 K. Dans le second cas,λD représente seulement une dizaine de rayons atomiques.
2.8.3 Retour à la forme intégrale : Fonction de Green

L'intérêt de la fonction de Green ici est de retrouver les anciens résultats. Pour cela nous devons utiliser
les résultats de la fonction de Green rappelé dans le chapitre préliminaire. Dans le cas de l'équation de
Poisson on a :
−∇2 G(r, r0 ) = δ(r − r0 ) (2.96)
Si on dénit la transformée de Fourier tel que :
ˆ
f˜(k) = f (r)e−ik.r dr (2.97)
On prend alors la transformée de Fourier de l'équation précédente on a alors :

0
e−ik.r
G̃(k, r0 ) = (2.98)
k2
Il faut à présent inverser la TF :
ˆ 0
0 1 eik.(r−r )
G(r, r ) = dk (2.99)
(2π)3 k2
Pour eectuer l'intégration, passons en coordonnées sphériques, où nous prenons l'axe kz parallèle à (r−r0 ).
Dans ce cas là on a k.(r − r0 ) = k||r − r0 || cos θ et dk = k 2 sin θdkdθdφ. L'intégrale s'écrit alors :
ˆ 2π ˆ π ˆ ∞
0 1 0
G(r, r ) = dφ dθ dk sin θeik||r−r || cos θ (2.100)
(2π)3 0 0 0
L'intégration selon φ nous donne directement un facteur 2π et par la suite on pose u = cos θ :
ˆ ∞ ˆ 1
0 1 0
G(r, r ) = 2
dk eik||r−r ||u du (2.101)
(2π) 0 −1
On intègre alors sur u : ˆ ∞

0 1 sin k||r − r0 ||
G(r, r ) = 2 dk (2.102)
2π 0 k||r − r0 ||
Un changement de variable évident nous donne :
ˆ ∞
0 1 sin k
G(r, r ) = 2 dk (2.103)
2π ||r − r0 || 0 k
Nous avons donc bien mis en évidence la dépendance en 1/|r − r0 |. L'intégrale n'est qu'une constante qu'il
est possible de calculer avec la théorie des fonctions analytiques et vaut π/2. Ceci nous donne alors :
1
G(r, r0 ) = (2.104)
4π||r − r0 ||
ρ(r)
Par la suite il reste plus qu'à ajouter le terme de source qui vaut S(r) = 0 et donc :
ˆ
1 ρ(r0 )
V (r) = dr0 (2.105)
4π0 ||r − r0 ||
On retrouve alors l'expression du potentiel d'une densité de charge que nous avons exprimé auparavant.
2.9 Énergie des distributions de charges discrètes

Le potentiel est relié à l'énergie électrostatique du système, essayons d'approfondir un peu cette notion
d'énergie électrostatique en comparant ce qu'il se passe pour des distributions discrètes et continues.
2.9.1 Théorème de l'extrémum du potentiel
Théorème :
Il n'existe pas d'extrémum du potentiel électrostatique V (r) dans une région de l'espace vide de
charge.
Démonstration :
Supposons qu'au voisinage d'un point M de l'espace il n'y a pas de charge. Soit Σ une surface fermée
entourant M . Raisonnons alors par l'absurde. Supposons alors que en M le potentiel V soit extremum.
Ainsi je peux trouver une surface Σ tel que le champ électrique E soit entrant (ou sortant) de cette surface
car on a E = −∇V . Par conséquence le ux de φ de E à travers Σ est non nul. Si c'est le cas il faudrait
que le volume délimité par Σ contienne une charge, or ceci est exclu par hypothèse. Donc V ne peut être
extrémum en M .
2.9. ÉNERGIE DES DISTRIBUTIONS DE CHARGES DISCRÈTES 37
Le corollaire direct est qu'on ne peut pas piéger une particule soumise uniquement à un champ électrosta-
tique puisqu'il ne peut pas avoir d'équilibre stable.
Théorème de Erauskaw :
Il n'existe aucune distribution de charges ponctuelles stable par pur eet électrostatique.
Démonstration :
Soit D une distribution discrète de charges. Pour avoir stabilité de la distribution il faut que pour toute
charge qi de la distribution, tout mouvement de cette charge entraine une force qui la ramène vers sa
position d'équilibre. Soit qi xé en ri , on doit avoir ED autour de ri qui est purement entrant ou sortant
suivant le signe de qi . Donc VD (ri ) est extrémum. Ceci n'est pas en contradiction car il y a une charge
qi . Cependant comme il n'y a pas d'interaction de i sur i on a VD (ri ) = VD−qi (ri ) or dans ce cas là il n'y
a pas de charge dans le volume autour de qi donc il n'y a pas d'extrémum possible donc pas de stabilité
possible.
2.9.2 Énergie d'une distribution discrète

Considérons une distribution de charges ponctuelles D composée de charges {qj }. Nous savons que si
la distribution est bornée le champ électrique créé par cette distribution est le champ électrostatique de
Coulomb ED . Prenons maintenant une autre charge qi à l'extérieur de D au point r. Nous connaissons la
force de la distribution sur la charge qi :
FD/qi = qi ED (r)
= −qi ∇(VD (r))
 
X q j
= −qi ∇  
4π0 ||r − rj ||
j
 
X
= −qi ∇  Vj (r)
j
Supposons maintenant que la charge qi soit amenée en r depuis l'inni en partant et en arrivant avec une
vitesse nulle. Soit WD le travail exercée sur qi par la distribution D. Soit Wop le travail que doit fournir
un opérateur extérieur pour eectuer l'opération. D'après le théorème de l'énergie cinétique on a :
∆EC (qi ) = 0 = Wop + WD (2.106)
On peut donc écrire :
Wop = −WD
ˆ r
= − Felecqi .dl
∞
ˆ r
 
X
= qi ∇ Vj (r) dr
∞ j
X
= qi (Vj (r) − Vj (∞))
j
Comme la distribution est bornée, le champ électrique décroit en 1/r2 et donc on peut choisir le potentiel
nul à l'inni Vj (∞) et donc : X
Wop = qi Vj (r) = qi VDj6=i (r) (2.107)
j
représente le travail par unité de charge à fournir par un opérateur extérieur pour
VD (r)
amener une charge depuis l'inni à une position donnée.
On peut à présent calculer l'énergie nécessaire pour créer D. Pour faire cela, on amène les charges les
unes après les autres. On note Wop
i le travail fourni pour amener la charge i depuis l'inni à sa position.
(1)
Wop =0 pas de charge présente (2.108)
(1/2) q2 q1
(2)
Wop = −Welec = (2.109)
4π0 ||r2 − r1 ||

q1 q2
(3)
Wop = q3 + (2.110)
4π0 ||r3 − r1 || 4π0 ||r3 − r2 ||
Si on considère N charges alors pour i ≤ N on a :
X qj
(i)
Wop = qi (2.111)
4π0 ||ri − rj ||
j<i
On en déduit alors l'énergie totale à apporter :
N
X
tot
Wop = i
Wop (2.112)
i=1
N
X X qj
= qi (2.113)
4π0 ||ri − rj ||
i=1 j<i
N
1 X qi qj
= (2.114)
2 4π0 ||ri − rj ||
i6=j
On vient de calculer Wop tot l'énergie d'interaction de la distribution D . Il faut noter que le signe de cette
quantité n'est pas xé. En eet si on empile une charge positive et une charge négative alternativement
on peut avoir Woptot < 0. On rappelle qu'on ne compte pas ici l'énergie d'interaction de i avec i.
2.10 Énergie des distributions continues

Cette fois-ci la distribution D est une distribution de charges continue. On transforme l'expression
précédente dans la limite continue et on introduit :
ˆ
1
ED = d3 r0 ρ(r0 )V (r0 ) (2.115)
2 D
où V (r0 ) est le potentiel électrostatique créé par toute la distribution. On peut se demander si ce passage
à la limite continue est légitime notamment en regardant si cette intégrale converge. Le potentiel étant en
1/r0 et l'élément de volume en r03 donc le produit sous l'intégrale est d'ordre de grandeur r02 et comme la
distribution est de taille nie l'intégral converge.
Localisation de l'énergie :
Si on regarde l'expression de ED on pourrait penser que l'énergie est seulement localisée dans le volume
occupé par la distribution. Qu'en est-il vraiment ? Pour répondre à cela, prenons un volume V délimité
2.11. STOCKAGE DE L'ÉNERGIE ÉLECTROSTATIQUE : LES CONDENSATEURS 39
par la surface Σ contenant la distribution de charges continue D. On a :

ˆ
1
ED = d3 r0 ρ(r0 )V (r0 )
2 D
ˆ
1
= d3 r0 ρ(r0 )V (r0 )
2 V
ˆ
1
= d3 r0 0 (∇r0 E(r0 ))V (r0 )
2 V
ˆ
0
= d3 r0 (∇r0 (E(r0 )V (r0 )) − E(r0 )∇r0 V
2 V
ˆ ˆ
0 3 0 0 2 0
= d r ||E(r )|| + d3 r0 (∇r0 (E(r0 )V (r0 ))
2 V 2 V
Prenons la seconde intégrale :

ˆ ‹
0 3 0 0 0 0
I2 = d r (∇r0 (E(r )V (r )) = d2 S.E(r0 )V (r0 ) (2.116)
2 V 2 Σ
Si on regarde la dépendance en r0 de l'argument de l'intégrale on a l'élément de surface en r02 , le champ

électrique en 1/r02 et le potentiel en 1/r0 . Ainsi l'argument de l'intégrale est en 1/r0 . Si on fait tendre la
vaille du volume V vers l'inni cette intégrale tend vers 0. Ainsi on a :
ˆ
0
ED = d3 r0 ||E(r0 )||2 (2.117)
2 R3
Á partir de cette expression, il est possible de faire plusieurs remarques :

• L'énergie de la distribution est l'intégrale sur tout l'espace car on a fait tendre le volume d'intégration
vers l'inni pour se débarrasser de la seconde intégrale. Ainsi l'énergie n'est pas seulement localisée
dans la distribution mais dans tout l'espace.
• L'énergie de la distribution est positive quelque soit la distribution que l'on considère. En eet, on
pourra essayer comme dans le cas ponctuel d'empiler des régions positives et négatives, on n'obtiendra
jamais une énergie négative. La diérence fondamentale entre les deux approches (ponctuelle et
continue) est que dans le cas de la densité continue, on considère l'interaction de la charge en un
point M avec elle même.
• On peut alors dénir une densité d'énergie électrostatique :
0
wel = ||E(r0 )||2 (2.118)
2
Ceci est un choix naturel au vu de l'expression mais il faut être conscient qu'il s'agit d'une possibilité
parmi d'autres de dénition de la densité d'énergie.
2.11 Stockage de l'énergie électrostatique : les condensateurs

Si on prend un condensateur plan composé de deux conducteurs plans de surface S distants de e (tel que
e S pour négliger les eets de bords). Les deux plans sont chargés avec des charges opposées ±Q. Si on
néglige les eets de bords, on décrit les plaques comme étant chargées uniformément en surface σ = ±Q S .
Le champ entre les armatures est donnée par :
σ Q
E= n= n (2.119)
0 0 S
où n est le vecteur normal aux plans du condensateur de l'armature chargé + vers l'armature chargée -.
On calcule alors la circulation de E pour avoir la diérence de potentiel et la capacité qu'on dénit comme
le rapport entre la charge d'une armature sur la diérence de potentiel :
Q
∆V = e (2.120)
0 S
La capacité d'un condensateur plan s'écrit alors :

Q 0 S
C= = (2.121)
∆V e
avec S la surface des armateurs et e l'espacement entre ceux-ci.
Sur des conducteurs chargés à l'équilibre, l'apport de nouvelles charges dQi de l'inni est associé à un
travail des forces électrostatiques :
X
dWel = − Vi dQi = −dWop (2.122)
i
Prenons un condensateur dans un circuit avec une armature au potentiel V1 avec une charge Q2 et l'autre
au potentiel V2 avec une charge Q2 = −Q1 .
dWel = V1 dQ1 + V2 dQ2 = (V1 − V2 )dQ1 (2.123)
Or on a 4 Q1 = C(V1 − V2 ) donc :
1
dWel = Q1 dQ1 (2.124)
C
Donc :
1 Q2 1
U= = CV 2 (2.125)
2 C 2
Pour un condensateur plan :
1 0 S 2
U= V (2.126)
2 e
Or E = V /e donc :
1
U = 0 (Se)E 2 (2.127)
2
Donc par unité de volume :
1
u = 0 E 2 (2.128)
2
Tout se passe comme si l'espace entre les plaques stockait une quantité d'énergie u = 12 0 E 2 par unité de
volume (vrai quelque soit la géométrie).
4. sous certaines conditions : voir BFR Tome 1 Chapitre 7

Chapitre 3
Magnétostatique
3.1 Introduction et rappels

Dans ce chapitre nous voulons traiter l'équivalent de l'électrostatique mais pour le magnétisme. On a
vu que dans le cas des régimes stationnaires (et non statique car sans courant il ne peut avoir de champ
magnétique) les équations pour le champ électrique et magnétiques étaient découplées. On a alors comme
équations locales pour le champ magnétique :
∇.B(r) = 0 (3.1)
∇ ∧ B(r) = µ0 j(r) (3.2)
Et également la force de Lorentz :
F(r) = qv ∧ B(r) = j ∧ B(r) (3.3)
Tout d'abord si on prend la divergence de l'équation de Maxwell-Ampère, le terme de gauche est nul
avec les relations vectorielles et on trouve que :
∇.j = 0 (3.4)
Ainsi la divergence du courant est nul ce qui signie que le courant est à ux conservatif en magnétostatique.
On aurait pu retrouver cette relation avec l'équation de conservation de la charge tout en sachant qu'en
régime stationnaire on a dρ
dt = 0.
Un peu d'histoire du magnétisme ainsi que quelques ordres de grandeur :

• Les roches possédant des propriétés magnétiques comme la magnétite (F e3 O4 ) sont connues depuis
l'âge de fer (il y a 2550-3000 ans environ).
• Jusqu'en 1800 il y a eu beaucoup d'expériences et de résultats qualitatifs (répultion entre pôles
identiques...) avant une étude plus systématique.
• En 1819 l'expérience d'Oersted met en évidence qu'un courant électrique peut inuencer un ensemble
alimentée comme une aiguille.
• Il s'en suit en 1822 les expériences de Faraday et Ampère qui donneront la loi de Faraday et le
théorème d'Ampère.
• Dans le même temps Biot et Savart propose la formule du champ magnétique créé par un courant.
La loi de Biot et Savart s'écrit

ˆ
µ0 j ∧ (r − r0 )
B(r) = d3 r0 (3.5)
4π D ||r − r0 ||3
avec j le courant volumique et µ0 la perméabilité du vide
41
42 CHAPITRE 3. MAGNÉTOSTATIQUE
Quelques ordres de grandeur en magnétisme :
? 1 Weber correspond au ux permettant par induction de créer une force électromotrice de 1 Volt si
on annule de façon entrante 1 champ magnétique en 1 seconde dans une spire.
? 1 Tesla correspond à un champ uniforme qui permet d'avoir un ux de 1 Weber par mètre carré.
? 1 Gauss = 10−4 T
? Champ magnétique terrestre : la composante horizontale vaut environ 0,2 G et la verticale 0,4 G
? Dans l'entrefer d'un électroaimant entre : 0,1-2 T
? Aimant permanent au Néodyme (Nd) : 1,25 T
? IRM : 6 T
? Champ pulsé : 100 T maximum sans destruction et jusqu'à 2800 T avec explosion
? Pulsar : 4.108
3.2 Lois locales et lois intégrales

3.2.1 Lois locales
Premièrement grâce aux équations de Maxwell :
• Comme la divergence de B est nulle on peut dire qu'il n'existe pas de monopôle magnétique.
D'après les relations vectorielles du première chapitre, il existe un certain champ vectoriel A tel
que
B=∇∧A (3.6)
Ce champ A est appelé potentiel vecteur et est déni à un gradient près car A et A + ∇ξ donnent
le même champ magnétique.
• Avec les équations de Maxwell on a également le lien entre le champ magnétique et la densité de
courant (en stationnaire) :
∇ ∧ B = µ0 j (3.7)
De plus on peut relier le potentiel vecteur à la densité de courant en remplaçant dans l'équation précédente
B par ∇ ∧ A :
∇ ∧ (∇ ∧ A) = ∇(∇.A) − ∆A = µ0 j (3.8)
Or à ce stade là nous n'avons pas imposé de condition sur le potentiel vecteur.
Le choix qui est fait est de retrouver une équation de Poisson pour le potentiel vecteur (comme
pour la densité de charge en électrostatique) et pour cela on impose la Jauge de Coulomb :
∇.A = 0 (3.9)
Ainsi le lien entre le potentiel vecteur et la densité de courant est l'équation de Poisson :
∆A + µ0 j = 0 (3.10)
Il faut noter que cette condition impose une condition sur le champ ξ qu'on peut ajouter au potentiel
vecteur :
∆ξ = 0 (3.11)
3.2. LOIS LOCALES ET LOIS INTÉGRALES 43
3.2.2 Conséquences : lois globales

En partant de l'équation de Poisson pour le potentiel vecteur, on peut eectuer le même calcul que pour
le potentiel électrostatique et on trouve en remplaçant ρ par j et 1/0 par µ0 (il sut de faire le calcul
analogue pour les trois composantes) :
ˆ
µ0 j(r0 )
A(r) = d3 r0 (3.12)
4π D ||r − r0 ||
Maintenant on peut calculer le champ magnétique à partir du potentiel vecteur :

ˆ
j(r0 )

µ0 3 0
B(r) = ∇r ∧ A(r) = d r ∇r ∧ (3.13)
4π D ||r − r0 ||
D'après le chapitre de rappels mathématiques on a ∇ ∧ (V.f ) = (∇f ) ∧ V + f (∇ ∧ V). Cependant le

second terme disparait dans le développement car le courant ne dépend que de r0 et que le ∇ agit sur r.
Donc :
j(r0 ) r − r0

1 0
∇r ∧ = ∇ r ∧ j(r ) = − ∧ j(r0 ) (3.14)
||r − r0 || ||r − r0 || ||r − r0 ||3
Ceci permet de retrouver la loi de Biot et Savart pour le champ magnétique B :
ˆ
µ0 j(r0 ) ∧ (r − r0 )
B(r) = d3 r0 (3.15)
4π D ||r − r0 ||3
La question qu'on doit se poser est de savoir si l'expression de A vérie la Jauge de Coulomb ? Calculons :
ˆ
j(r0 )

µ0 0
∇r .A(r) = dr ∇r . (3.16)
4π D ||r − r0 ||
On peut alors développer la divergence : ∇.(f V) = f (∇.V) + (V.∇)f ce qui donne :

ˆ
µ0 0 0 1
∇r .A(r) = dr j(r ).∇r (3.17)
4π D ||r − r0 ||
A présent nous allons eectuer des manipulations mathématiques pour arriver à nos ns. Tout d'abord
nous allons nous restreindre à une domaine de l'espace D délimité par une surface Σ. D est choisi tel que
de ce domaine il ne sort pas de courant, ainsi pour tout vecteur normal on a j.n = 0. On parle alors de
guidage de courant. De plus nous allons changer le gradient de r vers r0 en eectuant le changement de
variable r0 → r et r → r0 ce qui donne un signe "-" dans l'expression globale. On a donc :
ˆ
µ0 0 0 1
∇r .A(r) = − dr j(r ).∇r 0
4π D ||r − r0 ||
ˆ
j(r0 )

µ0 0
= − dr ∇r0 . car ∇.j = 0
4π D ||r − r0 ||
˛
µ0 j(r0 )
= − dΣ.
4π Σ ||r − r0 ||
˛
µ0 j(r0 ).n
= − dΣ
4π Σ ||r − r0 ||
= 0 par guidage
Donc dans le cas de la magnétostatique et du guidage du courant, le potentiel vecteur vecteur obéit bien
à la Jauge de Coulomb. Il faut noter cependant que contrairement au potentiel électrique V , le potentiel
vecteur A n'ore pas vraiment d'avantage technique de calcul.
On peut alors vérier la valeur de ∇ ∧ B :
∇∧B = ∇∧∇∧A
= ∇(∇.A) − ∇2 A
= −∇2 A
ˆ
µ0 0 0 2 1
= − dr j(r )∇r
4π D ||r − r0 ||
ˆ
µ0
= − dr0 j(r0 )(−4πδ(r − r0 ))
4π D
= µ0 j(r)
On retrouve bien l'équation de Maxwell Ampère en stationnaire.
3.2.3 Théorème d'Ampère

On part de l'équation de Maxwell-Ampère ∇ ∧ B = µ0 j. Prenons une courbe fermée orientée C et une
surface Σ sur laquelle s'appuie cette courbe. Calculons la circulation du champ magnétique le long de la
courbe C :
ˆ ˆ
B.dl = (∇ ∧ B).dΣ par Stokes
C Σˆ
= µ0 j.dΣ par Maxwell-Ampère

Σ
= µ0 I
Le théorème d'Ampère s'écrit alors :

ˆ
B.dl = µ0 I (3.18)
C
avec C une courbe fermée et I la somme des courants algébriques à travers n'importe quelle surface
s'appuyant sur C . Il faut noter qu'il est nécessaire d'orienter la courbe pour bien appliquer le
théorème : la somme des courants est algébrique.
Autrement dit : la circulation du champ magnétique sur une courbe fermée C est égale au produit
de la perméabilité magnétique du vide µ0 par l'intensité totale qui traverse toute surface Σ
s'appuyant sur la courbe.
Comme pour le théorème de Gauss, l'utilisation du théorème d'Ampère suppose un choix judicieux du
contour le long duquel on calcule la circulation de B et susamment de symétries/invariances.
3.3 Symétries et calculs de champs

3.3.1 Symétries du champ magnétique
Dans le chapitre sur l'électrostatique on a déjà vu les syémtries du champ électrique qui est un "vrai"
vecteur. Cependant pour le champ magnétique ce n'est pas pareil. En eet, si on regarde la loi de Biot et
Savart pour le potentiel vecteur A, on voit que A est un vrai vecteur car il est proportionnel au courant
(et donc à une vitesse) qui est un vrai vecteur.
3.3. SYMÉTRIES ET CALCULS DE CHAMPS 45
Le lien entre le champ magnétique et le potentiel vecteur est B = ∇ ∧ A. Or l'opérateur nabla est un
vrai vecteur également car la transformation (x, y, z) → (−x, −y, −z) transforme ∇ en −∇. Ainsi comme
le champ magnétique est le produit vectoriel entre deux vrais vecteurs, il est un pseudo vecteur 1 .
En conséquence, les propriétés de symétrie du champ magnétique sont "opposés" à celles du

champ électrique qui est un vrai vecteur.
• En tout point M d'un plan de symétrie Π+ de la distribution de courant on a B(M ) ⊥ Π+
• Si Π− est un plan d'antisymétrie de la distribution des courants alors B(M ) ∈ Π− pour
tout point M ∈ Π− .
3.3.2 Exemples de calculs du champ B

Fil rectiligne inni par théorème d'Ampère
On se place dans les coordonnées cylindriques. Tout plan contenant le l et le point M d'observation
est plan de symétrie donc le champ magnétique est orthoradial (selon eθ ). Par la suite on se place à une
distance r du l et on applique le théorème d'Ampère sur un cercle d'axe le l et de rayon r. On a alors :
ˆ
B.dl = B2πr = µ0 I (3.19)
C
Et donc :
µ0 I
B= eθ (3.20)
2πr
Fil rectiligne inni par Biot et Savart

Essayons maintenant de faire le calcul directement avec la loi de Biot et Savart. On doit donc calculer
le champ à partir de la formule pour un circuit liforme :
ˆ
µ0 r − r0
B= Idl ∧ (3.21)
4π C |r − r0 |3
On ne refait pas les symétries qui sont les mêmes que précédemment. Considérons la Figure 3.1 avec
les notations. Calculons le champ magnétique au point M en calculant tout d'abord le champ innitésimal
créé par le point P . Considérons que le point O est la référence pour notre altitude (z = 0) et que le point
P est à l'altitude z . Le l est parcouru pour un courant I . Exprimons à présent quelques éléments de
longueur :
R
|r − r0 | = (3.22)
cos θ
OP = R tan θ = z (3.23)
Rdθ
dz = (3.24)
cos2 θ
On peut alors exprimer :
dl ∧ (r − r0 ) = Rdz (3.25)
Le calcul donne alors : ˆ π/2
µ0 I µ0 I π/2
B= cos θdθ = [sin θ]−π/2 (3.26)
4π R −π/2 4π R
Et donc on trouve le même résultat :
µ0 I
B= eθ (3.27)
2πR
1. Un pseudo-vecteur est un vecteur dont le sens dépend du trièdre de référence choisi.Il se comporte comme un vecteur
lors d'une rotation (il conserve les angles orientés) mais change lors d'une isométrie indirecte comme une inversion de tous
les axes par rapport à un plan.
dl = dz ez
P
θ
O M
R
Figure 3.1 Champ magnétique créé par un l inni.
Spire circulaire - Champ sur l'axe

Considérons une spire C de rayon R et d'axe Oz parcouru par un courant I comme sur la Figure 3.2
Calculons le champ magnétique crée par cette spire sur son axe de révolution.
z
B
M
α
I
r
O
R
P
dl
Figure 3.2 Champ magnétique créé par une spire sur son axe de révolution.
Débutons par les symétries. Tout plan contenant l'axe de révolution (et donc le point M ) est un plan
d'antisymétrie du courant et donc le champ magnétique appartient à ce plan et donc B = Bez . Ici on est
obligé d'eectuer le calcul avec la loi de Biot et Savart avec r = PM :
˛ ˛
µ0 dl ∧ r Iµ0
B(r) = I 3
= dl ∧ r (3.28)
4π C r 4πr3
Et donc :
µ0 I R 2 µ0 I
B(r) = ez = sin3 αez (3.29)
2 r3 2R
Grâce à ce résultat on peut calculer le champ magnétique au point M créé par un solénoïde d'une
longueur L comprenant N spires. Pour cela découpons le solénoïde en spires innitésimales de longueur
dz parcourue par un courant dI = nI dz où n = N/L la densité linéïque de spires. Le champ au point M
sur la spire créée par cette spire est, d'après le résultat précédent est égale à
µ0 dI µ0 nIdz
dB = sin3 α = sin3 α (3.30)
2R 2R
avec α l'angle sous lequel on voit la spire. On a la relation tan α = R
z qu'on peut diérentier :
dα dz −Rdα −Rdα
= −R 2 ⇒ dz = = (3.31)
cos2 α z cos2 α tan2 α sin2 α
On peut alors intégrer le champ sur l'ensemble du solénoïde avec α1 et α2 les deux angles sous lesquelles
on voit le solénoïde : ˆ α2
Iµ0 1
B= n(− sin α)dα = µ0 nI (cos α2 − cos α1 ) (3.32)
α1 2 2
Et donc :
1
B = µ0 nI (cos α2 − cos α1 ) (3.33)
2
1
B = µ0 nI (cos α2 − cos α1 )ez (3.34)
2
Conducteur cylindrique
Soit un conducteur cylindrique de rayon R parcouru par un courant uniforme j. Calculons le champ
magnétique en tout point de l'espace. Comme pour le l inni, tous les plan contenant le point M et l'axe
de révolution du cylindre est plan de symétrie donc le champ magnétique est orthoradial. L'invariance par
rotation autour de l'axe de révolution et selon le l (coordonnée z ) nous informe que le champ ne dépend
que de la distance à l'axe noté r.
Appliquons alors le théorème d'Ampère sur un cercle de rayon r et de centre passant par l'axe de
révolution du l : ˆ
B.dl = 2πrB = µ0 Ie (3.35)
C
où Ie est le courant enlacé qui va dépendre de où on se situe dans l'espace. Si r > R alors le courant
2
enlacé est le courant total qui vaut I = jπR2 et donc le champ vaudra µ0 R2r j . Si r < R alors le courant
enlacé est une fraction du courant total et vaut I = jπr2 et donc le champ vaut µ0 2r j . Au nal, le champ
magnétique créé vaut :
r
B = µ0 jeθ r<R (3.36)
2
R2
= µ0 jeθ r>R (3.37)
2r
On peut également s'entrainer à partir de l'expression précédente à calculer le potentiel vecteur. Pour cela
il faut reprendre l'expression du rotationnel en cylindrique, remarquer que le champ magnétique n'est que
selon eθ et que A ne peut dépendre de z , alors on a :
∂Az
B=∇∧A=− eθ (3.38)
∂r
Il sut donc d'intégrer les deux relations précédentes pour r < R et r > R pour trouver les relations
suivantes, en imposant un potentiel vecteur nul en r = R :
R2 r
A = µ0 jln ez r>R (3.39)
2 R
j
= µ0 (R2 − r2 )ez r<R (3.40)
4
On remarquera dans cette expression que le potentiel vecteur ne tend pas vers 0 à l'inni. Cela devrait
poser problème mais en réalité comme la distribution de courant n'est pas borné le potentiel vecteur n'a
aucune raison de décroître à l'inni.
Champ magnétique dans un câble coaxial [Exercice]

On considère un câble coaxial (câble de sortie d'un générateur basse fréquence) constitué d'un cylindre
métallique central plein de rayon R1 , et d'une couche cylindrique périphérique de rayon interne R2 et de
rayon externe R3 .
Entre R1 et R2 se trouve une matière isolante assimilable à du vide du point de vue électromagnétique. Ce
câble est rectiligne d'axe Oz et considéré comme inniment long. Sa partie conductrice centrale transporte
une intensité i constante, dirigé selon +ez et sa partie périphérique transporte une intensité i dans le sens
−ez . Dans chacune des deux parties conductrices, la densité de courant est supposée uniforme. Soit M un
point où on calcul le champ magnétique créé : on note r la distance entre M et l'axe de l'ensemble.
1. Montrer que le champ magnétique créé par le câble en un point M quelconque se met sous la forme
B(M ) = B(r)u (3.41)
où u est un vecteur à préciser.

2. Déterminer l'expression de B(r) pour r ∈ [0; +∞[. Tracer B(r).
Dans la suite on considère le cas où les deux intensités sont surfacique, cantonnées dans de très
faibles épaisseurs au voisinage de R1 et R2 .
3. Que devient le graphe de B(r) ? Commenter les discontinuités qui apparaissent.
4. Après avoir rappelé l'expression de la densité volumique d'énergie magné"tique du champ, exprimer
l'énergie magnétique Um contenue dans une portion de longueur ` du câble.
5. En déduire le coecient d'auto-inductance (inductance propre) L de cette portion, puis l'inductance
linéïque propre du câble :
L
Λ= (3.42)
`
6. Calculer numérique L pour R1 = 1.0 cm, R2 = 2.0 cm et h = 1.0 m. Commenter.
Correction :
1. Soit un point M où on souhaite exprimer le champ. Le plan qui contient M et l'axe Oz est plan
de symétrie pour la distribution de courant, donc le champ est perpendiculaire à ce champ. Il est
donc dirigé selon eθ . La distribution est invariante par translation selon z et rotation selon θ donc
le champ ne dépend que de r. Ainsi :
B(M ) = B(r)eθ (3.43)
2. Les lignes de champ sont des cercles d'axe Oz sur lesquels B(r) est constante. On peut calculer B(r)
en utilisant le théorème d'Ampère. Le contour fermé utilié est une ligne de champ magnétique de
rayon r, orienté arbitrairement dans le sens trigonométrique :
˛
B(r).dl = B(r)2πr = µ0 Ienlace (3.44)
C
Il reste à déterminer l'intensité enlacée par le contour en fonction du rayon r. Par dénition cette
intensité traverse la surface Σ s'appuyant sur le contour C . Donc :
¨ ˆ r
Ienlace = j(r).dS = j(r)2πrdr (3.45)
Σ 0
j(r) est selon z et dépend de l'endroit considéré.

• Cas 1 : si r ∈ [0; R1 ], j est tel que i = jπR12 , donc j(r) = et Ienlace = i Rr 2 donc :
i 2
πR12 1
µ0 ir
B(r) = ∀r ∈ [0; R1 ] (3.46)
2πR12
• Cas 2 : r ∈ [R1 : R2 ], j(r) = 0 et donc l'intensité enlacée est simplement l'intensité circulant
dans le cylindre central i. Donc :
µ0 i
B(r) = ∀r ∈ [R1 ; R2 ] (3.47)
2πr
• Cas 3 : r ∈ [R2 : R3 ], j est tel que −i = jπ(R32 − R22 ) donc j = −i
π(R32 −R22 )
et donc :
ˆ r
−i
Ienlace = i + 2 2πrdr (3.48)
R2 π(R3 − R22 )
On trouve donc :
r2 − R22

µ0 i
B(r) = 1− 2 ∀r ∈ [R2 ; R3 ] (3.49)
2πr R3 − R22
• Cas 4 : si r ∈ [R3 ; +∞], j = 0. L'intensité enlacé est égale à i − i = 0 donc le champ magnétique
est nul à l'extérieur du câble coaxial.
Le graphique donne d'abord une croissance linéaire jusqu'au point R1 , 2πR . Puis la courbe décroit

µ0 i
1
en 1/r jusqu'au point R2 , 2πR et enn décroit selon la fonction calculée jusqu'à 0 en R3 où il le

µ0 i
2
reste nul jusqu'à l'inni. Une représentation est faite ci-après :
3. Dans le cas où l'intensité du câble est cantonnée au voisinage de R1 et de R2 alors on a plus la

croissance linéaire et la seconde décroissance mais à chaque fois des sauts dans la valeur du champ
magnétique. La décroissance entre R1 et R2 reste inchangé.
4. La densité volumique magnétique d'écrit um = B2

2µ0 et s'exprime en J.m−3 . Elle est non nulle que
dans l'intervalle [R1 ; R2 ] où elle vaut :
2
1 µ0 i
um = (3.50)
2µ0 2πr
Pour avoir l'énergie magnétique dans la zone considérée, on intègre sur le volume entre z = 0 et
z=`: ˆ ` ˆ R2
1 R2
Um = dz dr2πrum ⇒ Um = µ0 `i2 ln (3.51)
0 R1 4π R1
En identiant les relation on obtient L = µ0 R2
2π `ln R1 . En divisant par la longueur ` :
µ 0 R2
Λ= ln (3.52)
2π R1
5. Avec les valeurs proposées on trouve L = 1.4 10−7 H.
3.4 Dipôle magnétique

Pour le dipôle magnétique, procédons par analogie avec le dipôle électrostatique et résumons les résultats
dans un tableau :
Dipôle en O Électrostatique Magéntique

Moment dipolaire p = pez M = M ez
Constante 1
0 µ0
Champ E B
2p cos θ 2µ0 M cos θ
Er = Br =
4π0 r3 4πr3
p sin θ µ0 M sin θ
Champ créé Eθ = Eθ =
4π0 r3 4πr3
µ0 3(M.r)r

1 3(p.r)r
E= − p E= −M
4π0 r3 r2 4πr3 r2
Énergie potentielle (champ extérieur) Ep = −p.Eext Ep = −M.Bext
∂Eext ∂Bext
Fx = p. Fx = M.
∂x ∂x
Force résultante dans un champ F = ∇(p.Eext )p=cte F = ∇(M.Bext )M =cte
F = (p.∇)Eext F = (M.∇)Bext
Moment résultant Γ(O) = p ∧ Eext Γ(O) = M ∧ Bext
Chapitre 4
Relations de passage
Une propriété qu'on peut déduire du régime stationnaire est l'expression des relations de passage des
champs lors du passage d'une interface. En eet, dans une description volumique, les champs et potentiels
sont continus. Les problèmes arrivent quand on cherche à passer à une description surfacique dans laquelle
l'épaisseur de la nappe de courant ou de charge tend vers 0.
Dans les deux cas on considérera une nappe d'épaisseur a tendant vers 0 et où la charge volumique ρ et
le courant volumique j est non nul. Le plan xy contient la nappe alors que l'axe z est perpendiculaire à
celle-ci.
4.1 Relation de passage pour le champ électrique

La relation de Maxwell-Gauss s'écrit en coordonnées cartésiennes :
ρ
∂x Ex + ∂y Ey + ∂z Ez = (4.1)
0
On considère une nappe chargée, innie selon les directions x et y et qu'on fait tendre vers 0 selon z .
Intégrons l'équation selon l'épaisseur z . On considère que la nappe à une épaisseur a et on fera tendre
cette quantité vers 0 :
ˆ a/2 ˆ
1 a/2
dz (∂x Ex + ∂y Ey + ∂z Ez ) = dzρ (4.2)
−a/2 0 −a/2
Les variations du champ électrique selon x et y ne seront pas aecté par le fait qu'on fasse tendre a vers
0. Ainsi leur intégrale tend vers 0 lorsque a tend vers 0. Par contre les variations selon z peuvent être très
intense dans une nappe inniment ne et ainsi son intégrale peut prendre une valeur nie lorsque a tend
vers 0. Ainsi on a :
σ
Ez (a/2) − Ez (−a/2) = a→0 (4.3)
0
qu'on écrit parfois en nommant les espaces 1 et 2 de par et d'autre de la surface :
σ
Ez2 − Ez1 = (4.4)
0
On dénit alors proprement la densité surfacique σ comme :

ˆ a/2
σ= dzρ a→0 (4.5)
−a/2
51
52 CHAPITRE 4. RELATIONS DE PASSAGE
Concernant les composantes tangentielles, il faut utiliser la relation de Maxwell-Faraday qui s'écrit en
coordonnées cartésiennes :
∂y Ez − ∂z Ey = −∂t Bx
∂z Ex − ∂x Ez = −∂t By
∂x Ey − ∂y Ex = −∂t Bz
Lors de l'intégration, les dérivées en dehors de z tendent vers 0 pour les mêmes raisons que précédemment
(même pour le temps) et on trouve avec les deux premières équations que :
Ey2 = Ey1 (4.6)

Ex2 = Ex1 (4.7)
La relation de passage pour le champ électrique s'écrit :
σ
E2 − E1 = n12 (4.8)
0
où n12 est le vecteur unitaire perpendiculaire à la nappe orienté de 1 vers 2
4.2 Relation de passage pour le champ magnétique

On procède de la même façon pour le champ magnétique. Comme la divergence du champ magnétique
est nul, on a en coordonnées cartésiennes :
∂x Bx + ∂y By + ∂z Bz = 0 (4.9)
Et par intégration sur la coordonnées z , avec l'épaisseur du plan qui tend vers 0 on obtient que :
Bz2 = Bz1 (4.10)
et donc la composante normale du champ magnétique est continue à la traversée d'une nappe de courant.
Avec la l'équation de Maxwell-Ampère on a :
1
∇ ∧ B = µ0 j + ∂t E (4.11)
c2
1
∂y Bz − ∂z By = µ0 jx + ∂t Ex
c2
1
∂z Bx − ∂x Bz = µ0 jy + 2 ∂t Ey
c
1
∂x By − ∂y Bx = µ0 jz + 2 ∂t Ez
c
Lorsqu'on intègre, toutes les dérivées en dehors de z tendent vers 0. Ainsi la dernière équation ne nous
apporte pas d'information par contre la première nous donne :
ˆ a/2
By2 − By1 = µ0 jx dz (4.12)
−a/2
4.2. RELATION DE PASSAGE POUR LE CHAMP MAGNÉTIQUE 53
En dénissant la composant selon x du courant surfacique comme :

ˆ a/2
jSx = jx dz a→0 (4.13)
−a/2
On trouve que :
By2 − By1 = µ0 jSx (4.14)
De même avec la deuxième relation sur le rotationnel on trouve que :
Bx2 − Bx1 = µ0 jSy (4.15)
La relation de passage pour le champ magnétique s'écrit alors :
B2 − B1 = µ0 jS ∧ n12 (4.16)
avec la même dénition de n12 .
En résumé on a :
À la traversée d'une surface chargée, la composante tangentielle de E est continue et sa compo-

sante normale subit une discontinuité égale à σ0 n où σ est la densité surfacique de charges.
À la traversée d'une nappe de courant, la composante normale de B est continue et sa composante

tangentielle subit une discontinuité égale à µ0 jS ∧ n où jS est la densité surfacique de courant.
54 CHAPITRE 4. RELATIONS DE PASSAGE
Chapitre 5
Approximation des régimes
quasi-stationnaires
Lorsqu'on ne se place plus dans le régime stationnaire et qu'il n'est plus possible de découpler le champ
magnétique et électrique alors la détermination des champs à partir des sources n'est plus aussi facile.
En eet une chose à prendre en compte lorsque nous ne sommes plus en stationnaire est la propagation
de l'information (et des champs) provenant de la source. Le champ champ varie avec un certain retard
si la source varie à cause de la distance entre les deux. L'Approximation des Régimes Quasi-Stationnaire
(ARQS) consiste à négliger la propagation et le retard de l'établissement. De plus, nous avons recours
parfois à négliger une partie des champs créé en fonction du régime d'ARQS dans lequel on se place :
magnétique ou électrique. Ce chapitre détaille les hypothèses et les conséquences de l'ARQS et met au
clair des choses qui ont été omises en première cycle universitaire.
5.1 Équations de Maxwell et potentiels

5.1.1 Pourquoi le courant de déplacement ?
La grande contribution de Maxwell aux équations de l'électromagnétisme est l'ajout d'un terme créant du
champ magnétique par variation temporelle du champ électrique : le courant de déplacement. La question
est de savoir ici pourquoi ce courant de déplacement est nécessaire dans ces équations.
Considérons un condensateur dans un circuit électrique qui est traversé par un courant I(t) comme sur
la Figure 5.1. Alors ce courant crée un champ magnétique B. La circulation de ce champ magnétique le
long d'une courbe C qui entoure le l alimentant le condensateur est donnée par le théorème d'Ampère et
vaut µ0 I . Mais ceci est le cas uniquement si on ne tient pas compte du courant de déplacement 0 ∂E
∂t .
Σ1
I(t) Σ2
Figure 5.1 Pourquoi le courant de déplacement ? Une explication avec un condensateur.

55
56 CHAPITRE 5. APPROXIMATION DES RÉGIMES QUASI-STATIONNAIRES
Si à présent nous prenons deux surfaces s'appuyant sur la courbe C . La première Σ1 coupe également le
l transportant le courant alors que le la seconde Σ2 passe entre les deux armatures du condensateur. Si
on calcule le courant traversant les deux surfaces, on trouve alors I(t) pour Σ1 et 0 pour Σ2 , alors qu'on
devrait trouver le même résultat !
Il manque donc quelque chose dans notre raisonnement ou dans les équations et on résout le problème
en considérant le courant de déplacement. En eet le champ électrique varie entre les armatures du
condensateur car à cause du courant I(t) le condensateur est entrain de se charger. Si on recalcule la
circulation pour les deux surfaces on trouve :

1 ∂E ∂E
µ0 I = S 2
= Sµ0 0
(5.1)
c ∂t ∂t
Et donc :
∂E
I(t) = S0
(5.2)
∂t
Le courant de déplacement permet de fermer le circuit. Si on se place en régime stationnaire le courant
de déplacement est nul car la variation de E est nul mais comme le condensateur est chargé on a aussi
I(t) = 0.
5.1.2 D'où proviennent 0 et µ0 ?

À première vue, on a besoin de 4 constantes pour décrire à la fois la relation entre E et B et la relation
entre le champ électromagnétique et les sources ρ et j :
∂B
∇ ∧ E = −α1
∂t
∇.E = α2 ρ
∇.B = 0
∂E
∇ ∧ B = α3 j + α4
∂t
Commençons par prendre le divergence de l'équation de Maxwell Ampère :

∂∇.E ∂ρ
∇.(∇ ∧ B) = 0 = α3 ∇.j + α4 = α3 ∇.j + α2 α4 (5.3)
∂t ∂t
Si on veut une conservation de la charge, on doit imposer alors α2 α4 = 1.
Prenons à présent une région de l'espace tel que j = 0 et prenons le rotationnel de l'équation de Maxwell-
Ampère :
∂E ∂2B
∇ ∧ (∇ ∧ B) = −∆B = α3 α4 = −α1 α3 α4 2 (5.4)
∂t ∂t
Et donc en remplaçant α4 par 1/α2 on obtient l'équation de propagation du champ magnétique :
α1 α3 ∂ 2 B
∆B − =0 (5.5)
α2 ∂t2
On obtient alors l'équation de d'Alembert de propagation d'une onde (on trouve la même équation pour
le champ électrique). Notons c la vitesse de la propagation de cette onde et on a alors la relation :
α1 α3 1
= 2 (5.6)
α2 c
5.1. ÉQUATIONS DE MAXWELL ET POTENTIELS 57
On a alors 4 constantes mais 2 relations entre elles il manque alors plus que deux constantes à xer et on
prend alors :
1
α3 = µ 0 = 0 (5.7)
α2
Par choix d'un système d'unités nous pouvons xer la valeur d'un des coecients. Choisissons alors α1 = 1
et on obtient alors la relation c2 0 µ0 = 1.
Il faut savoir que la forme des équations de Maxwell qu'on connait dépend uniquement du système
d'unité qu'on a choisi. En eet si on prend par exemple le système d'unités anglais, appelé système d'unité
gaussien (basé sur le cgs pour centimètre-gramme-seconde) les équations de Maxwell prenne une forme
diérente :
1 ∂B
∇∧E = −
c ∂t
∇.E = 4πρ
∇.B = 0
1 ∂E
∇∧B = 4πj +
c ∂t
où on a au nal :
1 4π
0 = µ0 = (5.8)
4πc c
5.1.3 Équations de propagation des potentiels et des champs

On part des équations de Maxwell dans le vide sans approximation. On a alors :
ρ(r, t)
∇.E(r, t) = (5.9)
0
∇.B(r, t) = 0 (5.10)
∂B(r, t)
∇ ∧ E(r, t) = − (5.11)
∂t
1 ∂E(r, t)
∇ ∧ B(r, t) = µ0 j(r, t) + (5.12)
c2 ∂t
Il faut à présent redénir les potentiels électromagnétiques car à présent les champs électriques et
magnétiques sont couplés. On ne change pas la dénition du potentiel vecteur A car l'équation de Maxwell-
Thomson reste inchangé. Par contre l'équation de Maxwell-Faraday change et on a :

∂A
∇∧ E+ =0 (5.13)
∂t
Et donc on déni le potentiel électrique par :
∂A
E = −∇V − (5.14)
∂t
De plus on choisit toujours la même convention pour qu'en stationnaire on retrouve le potentiel électro-
statique.
Maintenant, nous pouvons établir les équations pour les potentiels. Pour le potentiel vecteur il faut dé-
marrer avec l'équation de Maxwell-Ampère :
1 ∂E
∇ ∧ B = µ0 j +
c2 ∂t
1 ∂ ∂A
∇ ∧ (∇ ∧ A) = µ0 j + 2 −∇V −
c ∂t ∂t
2
2 1 ∂ A ∂V
∇(∇.A) − ∇ A = µ0 j − 2 +∇
c ∂t2 ∂t
en réorganisant les termes on trouve l'équation de propagation du potentiel vecteur A :
1 ∂2A

1 ∂V
∇2 A − = −µ 0 j + ∇ ∇.A + (5.15)
c2 ∂t2 c2 ∂t
Pour le potentiel scalaire V on utilise l'équation de Maxwell-Gauss :

ρ
∇.E =
0

∂A ρ
∇. −∇V − =
∂t 0
∂ ρ
− ∇.A − ∇2 V =
∂t 0
On pourrait s'arrêter ici mais pour avoir la même forme d'équation pour les deux potentiels on ajoute des
2
deux côtés le terme c12 ∂∂tV2 et en réorganisant les termes on trouve :
1 ∂2V

ρ ∂ 1 ∂V
2
∇ V − 2 2 =− − ∇.A + 2 (5.16)
c ∂t 0 ∂t c ∂t
Nous avons deux équations de propagation pour les potentiels électromagnétiques qui sont cou-
plés. Cependant nous n'avons pas encore choisi de conditions de jauge pour ceux-ci. Au lieu de
garder la jauge de Coulomb du cas stationnaire nous allons choisir judicieusement une jauge pour
découpler les deux équations. Nous avons déjà mis les équations dans la bonne forme pour qu'on
impose comme jauge :
1 ∂V
∇.A + 2 =0 (5.17)
c ∂t
Cette jauge est appelé jauge de Lorenz a .
a. Il s'agit du physicien danois Ludvig Lorenz et non le physicien néerlandais Hendrik Lorentz auteur des
transformations de Lorentz.
On trouve donc pour équations de propagation :
1 ∂2A
∇2 A − = −µ0 j (5.18)
c2 ∂t2
1 ∂2V ρ
∇2 V − =− (5.19)
c2 ∂t2 0
5.1. ÉQUATIONS DE MAXWELL ET POTENTIELS 59
Qu'en est-il des équations de propagation pour les champs ? Prenons le rotationnel de l'équation de
Maxwell Faraday :

∂B
∇ ∧ ∇ ∧ E = −∇ ∧
∂t
∂
∇(∇.E) − ∇2 E = − (∇ ∧ B)
∂t
∇ρ 2 ∂ 1 ∂E
−∇ E = − µ0 j + 2
0 ∂t c ∂t
Soit :
1 ∂2E ∂j ∇ρ
∇2 E − 2 2
= µ0 + (5.20)
c ∂t ∂t 0
On peut également montrer que le champ magnétique obéit à l'équation de propagation
1 ∂2B
∇2 B − = −µ0 ∇ ∧ j (5.21)
c2 ∂t2
Ainsi que ce soit pour A, B, E ou V (en jauge de Lorenz) on trouve toujours des équations de la
forme :
1 ∂2

2
∇ − 2 2 ψ=S (5.22)
c ∂t
où S sont les sources connues.
Les remarques qu'il est possible de faire à ce niveau là sont :

• Si on se place en stationnaire, les dérivées temporelles disparaissent et on retrouve les équations de
Poisson pour les potentiels
• La jauge de Lorenz en stationnaire sans la dérivée temporelle permet de retrouver la jauge de Coulomb
• Une image plus complète de l'électromagnétisme commence à apparaitre où les champs électrique et
magnétique sont couplés et donc indissociables.
5.1.4 Solutions
Pour trouver les solutions une méthode est d'utiliser les fonctions de Green. Si on a un domaine des
charges et des courants qui est localisé alors les champs décroissent à l'inni et ainsi admettent une
transformée de Fourier. La solution par la fonction de Green est alors unique.
On doit alors résoudre :
1 ∂2G
∇2 G − 2 2 = δ(r − r0 )δ(t − t0 ) (5.23)
c ∂t
Pour cela on passe dans l'espace de Fourier avec pour dénition de la transformée de Fourier :
ˆ
0 0
G(r, r0 , t, t0 ) = dqdνG(q, ν)e2iπ[q.(r−r )−ν(t−t )] (5.24)
on obtient :
ν2

2
−q + 2 G(q, ν) = 1 (5.25)
c
Ceci donne par transformée de Fourier inverse 1 :

||r − r0 ||

0 0 −1 0
G(r, r , t, t ) = δ t − t− (5.26)
4π||r − r0 || c
Formulé ainsi, on a les solutions avancées et retardées. Seule la solution retardée respecte le principe de
causalité et on l'écrit :
−cΘ(t − t0 )
G(r, r0 , t, t0 ) = δ ||r − r0 || − c|t − t0 | (5.27)

0
4π||r − r ||
L'expression de la solution pour le potentiel retardé est donné par la convolution de la fonction de Green
avec la source. On obtient alors :
ˆ 0 ||

ρ r0 , t − ||r−r
c
V (r, t) = dr0 (5.28)
D 4π0 ||r − r0 ||
ˆ ||r−r0 ||

µ0 j r 0 , t − c
A(r, t) = dr0 0
(5.29)
D 4π||r − r ||
Par la suite pour remonter au champ électromagnétique il sut d'utiliser les expressions
∂A(r, t)
B(r, t) = ∇ ∧ A(r, t) puis E(r, t) = −∇V (r, t) − (5.30)
∂t
Á ce stade, on peut faire plusieurs remarques :
• Même si les équations de Maxwell couplent les champs E et B, on voit que le choix judicieux de
la jauge permet à V d'être toujours une expression de la densité de charge ρ et A toujours une
expression de la densité de courant j. Le couplage entre les deux potentiels est ici est caché à la fois
dans la jauge de Lorenz qu'on a choisi mais également dans l'expression de E
• Les expressions des potentiels retardés satisfont la jauge de Lorenz. On pourra vérier en exercice
que c'est bien le cas.
• En principe si on connait parfaitement les distributions des sources ρ et j on connait tout grâce aux
expressions des potentiels retardés qui permettent de remonter aux champs. Mais le problème qui va
se poser est un problème mathématique avec des calculs qui deviennent rapidement infaisables...
Essayons de nous placer dans le cas où les toutes les grandeurs varient lentement. Prenons une certaine
distribution de charges et de courants D et plaçons nous en un point M situé à une distance caractéristique
L de la source. On note τ le temps caractéristique de variation des sources ρ et j. Sous quelles conditions
peut-on négliger le retard dans le potentiel ?
Le retard est majoré par L/c et si on a un temps caractéristique τ de variation on a que si τ Lc alors
on peut écrire
||r − r0 ||

ρ r, t − ' ρ(r, t) (5.31)
c
L'Approximation des Régimes Quasi-Stationnaire (ARQS) consiste à se placer dans le cas tel que :
L
=1 (5.32)
τc
Il est important de s'interroger sur deux points essentiels dans l'ARQS :

1. le calcul de la transformée de Fourier inverse est assez complexe nécessitant des outils d'analyse mathématique sur les
fonctions complexes. Le calcul en lui même ne présente pas d'intérêt pour la préparation à l'agrégation c'est pourquoi il n'est
pas présenté ici.
5.2. DÉMARCHE NAÏVE DE L'ARQS 61
• Comment faire cette approximation de manière cohérente et contrôlée ?

• Quelles sont les conséquences pour les champs électromagnétiques ?
5.2 Démarche naïve de l'ARQS

Essayons tout d'abord de montrer de manière naïve que dans l'ARQS on a toujours ∇ ∧ B = µ0 j. Pour
cela essayons deux méthodes 2 .
5.2.1 Méthode 1 : en utilisant le calcul de magnétostatique

On a montré que dans l'ARQS que la propagation des champs étaient négligé et donc que le potentiel
vecteur pouvait s'écrire :
˚
µ0 j(r0 , t)
A(r, t) = d3 r0 = ABS (r, t) (5.33)
4π R3 ||r − r0 ||
Par la suite il sut de prendre le rotationnel du potentiel vecteur :
˚
µ0 j(r0 , t) ∧ (r − r0 )
∇∧A= d3 r0 (5.34)
4π R3 ||r − r0 ||3
Il sut alors d'eectuer les mêmes calculs que ceux menés dans le chapitre de magnétostatique pour
obtenir :
∇.B = 0 et ∇ ∧ B = µ0 j (5.35)
Donc nous venons de montrer qu'en ARQS on a toujours
∇ ∧ B = µ0 j (5.36)
5.2.2 Méthode 2 : avec des ordres de grandeur

Soit E l'ordre de grandeur de la norme de E et B celui de la norme de B. Les champs varient sur une
échelle d avec d < L (on rappelle que L est la distance d'observation caractéristique). Si on écrit en ordre
de grandeur l'équation de Maxwell-Faraday ∇ ∧ E = − ∂B ∂t on obtient :
E B E d L
' ⇒ ' < (5.37)
d τ B τ τ
An de pouvoir négliger le courant de déplacement dans l'équation de Maxwell-Ampère, il est nécessaire
de comparer le terme de dérivée temporel de E avec le rotationnel de B (comme on ne suppose pas avoir
un courant pour le moment), on a alors :

1 ∂E 1 E
c2 ∂t c2 τ E d E L E
' B = 2τ
< 2τ
= 1 (5.38)
∇ ∧ B B c B c cB

d
En première approximation (c'est vrai pour une onde plane progressive monochromatique) E et cB sont du
même ordre de grandeur et donc on peut négliger le terme électrique dans l'équation de Maxwell Ampère
et on obtient toujours en ARQS :
∇ ∧ B = µ0 j (5.39)
2. On pourra consulter le BUP par A.Domps de Février 2003, Volume 97(2) Remarques sur l'approximation des régimes
ainsi que le BUP de G. Roussaux et A.Domps de Novembre 2004, Volume 98(2)
quasi-stationnaires en électromagnétisme
Remarques supplémentaires sur l'approximation des régimes quasi-stationnaires en électromagnétisme
5.2.3 Contre exemple et erreurs commises

Cette approche naïve à l'air de fonctionner au premier abord et surtout dans tous les cas de gure.
Cependant si on prend un condensateur plan circulaire 3 en régime forcé de pulsation ω on montre qu'à
une distance r du centre on a : rω
E = E0 J0 ez (5.40)
c
iE0 rω
B= J1 eθ (5.41)
c c
où J0 et J1 sont des fonctions de Bessel. En régime de basses fréquences, avec les équivalents pour les
fonctions de Bessel, on trouve :
E ' E0 ez (5.42)
E0 irω
B' eθ (5.43)
c 2c
On remarque que E est similaire au champ statique mais que B est non nul alors que dans ce cas on a
j = 0. Donc on n'a pas toujours ∇ ∧ B = µ0 j comme on a pu le montrer précédemment. Cela est aussi
cohérent avec l'introduction du chapitre montrant que le courant de déplacement permettait de fermer un
circuit avec un condensateur à l'intérieur.
Où sont les erreurs qu'on a commises ?

• Pour la première méthode, il faut refaire le calcul entièrement pour montrer que l'expression du
potentiel vecteur ne redonne pas l'équation de Maxwell Ampère sans le terme électrique. En eet
dans le cas de l'ARQS, on est avec la jauge de Lorenz et non de Coulomb or on a supprimé un terme
dans le développement du double produit vectoriel de A qui ne disparait pas dans le cas de la jauge
de Lorenz.
• Lorsqu'on écrit les termes en ordre de grandeur il faut faire attention. Quand on dit que ||∇ ∧ E|| '
E/d il faudrait plutôt dire que ||∇ ∧ E|| . E/d car on majore la norme du produit vectoriel. Dans ce
cas là on obtient :
E d
& (5.44)
B τ
Et donc on ne peux plus dire que :
E L
. (5.45)
B τ
Ainsi les démarches qu'on a faite ici sont des démarches qu'on pourrait être tenté de faire mais qui sont
naïves et non rigoureuses et amène ainsi à un résultat faux. Nous allons maintenant dénir proprement
l'ARQS et les conditions d'application.
5.3 Démarche rigoureuse de l'ARQS

5.3.1 Introduction
Les équations de Maxwell contiennent deux dérivées temporelles :
∂B
∇∧E = − (5.46)
∂t
1 ∂E
∇ ∧ B = µ0 j + (5.47)
c2 ∂t
et elles sont responsables des phénomènes de propagation qu'on néglige dans l'ARQS.
3. Cela fonctionne avec un condensateur qui n'est pas forcément circulaire mais l'expression avec les fonctions de Bessel
sont plus facile à manipuler pour l'exemple ici.
5.3. DÉMARCHE RIGOUREUSE DE L'ARQS 63
L'idée directrice de l'ARQS est de retrouver une relation entre les champs et leurs sources identiques
à celle qu'on trouve en statique. Le champ suit les variations lentes des sources : le temps typique de
variation τ est long devant le temps de propagation L/c
On souhaite comparer l'importance des diérentes quantités du problèmes pour faire correctement les
approximations voulues dans les équations de Maxwell. On a les quantités sans dimension suivantes :
ρc cτ cB cA
, , ,
j L E V
On introduit ρ̃ et j̃ les ordres de grandeur relatif à ρ et à j dans le système considéré. On suppose que ρ,
j, E et B varient spatialement sur une échelle caractéristique 4 d.
On introduit alors la quantité :

j̃
ξ= (5.48)
ρ̃c
qui va rendre compte si le système est plutôt un type conducteur pour ξ 1 ou type condensateur
ξ1
Pour les potentiels si on appelle Ω la volume caractéristique de la distribution de charges et de courants

alors on a pour les ordres de grandeur des potentiels calculés en un point M qui est dans l'ARQS tel que
cτ 1 :
L
µ0 j̃ 1 Ωρ̃
A= Ω , V= (5.49)
4π L 4π0 L
On peut maintenant comparer ces deux grandeurs et plus particulièrement cA avec V qui ont les mêmes
dimensions :
cA µ0 c Ωj̃ 4π0 L j̃
' ' =ξ (5.50)
V 4π L Ωρ̃ ρ̃c
On voit alors que le paramètre ξ permet à la fois de comparer les sources et les potentiels des champs
associés et ainsi pouvoir préciser dans quel régime on se place en fonction de la valeur de ξ .
5.3.2 Régime quasi-stationnaire électrique ( 1, ξ 1)

Plaçons nous dans le cas où :
L cA j̃
=1 et = =ξ1 (5.51)
τc V ρ̃c
Ceci signie que les charges dominent sur les courants. Prenons le champ électrique qui s'écrit :
∂A
E=− − ∇V (5.52)
∂t
et comparons les deux termes qui apparaissent :
||∇V || V 1 V τc V τc 11
∂A ' = > = 1 (5.53)

∂t
d A/τ cA d cA L ξ
Donc la variation temporelle du potentiel vecteur est négligeable par rapport au gradient du potentiel V
et on peut écrire :
E ' −∇V ⇒ ∇∧E'0 (5.54)
4. Si les échelles de variation sont diérentes pour chaque grandeur on majorera ou on minorera ces échelles pour eectuer
les ordres de grandeur. Donc supposer une seule grandeur caractéristique n'enlève rien à la généralité du problème.
Ainsi cette approximation revient à négliger le terme ∂B

∂t dans l'équation de Maxwell-Faraday.
Qu'en est-il des champs électrique et magnétique ? Calculons le rapport entre les champs électrique et
magnétique :
cB cA/d cA
≤ = =ξ1 (5.55)
E V/d V
Ainsi dans la force de Lorentz on ne garde que le terme électrique. :
F = qE (5.56)
On en déduit alors en terme d'énergie :
B2
cB E ⇒ c2 B 2 E 2 ⇒ E2 (5.57)
µ0 0
Et donc :
1 1 2
0 E 2 B (5.58)
2 2µ0
Dans l'ARQS électrique, l'énergie est principalement stocké sous forme électrique.
Il nous reste une équation que nous n'avons pas essayer de simplier, l'équation de Maxwell-Ampère. Si
on compare les deux courants qui apparaissent :
||µ j|| µ0 j̃ µ0 j̃ 0 τ c2 j̃ τ c j̃ τ c ξ
1 0∂E ' dρ̃ = = > = (5.59)
2
c ∂t 2
ρ̃ d ρ̃c d ρ̃c L
0 τ c
On sait que ξ et sont très petits devant 1 mais on ne peut pas les comparer donc on ne peut pas conclure
si un des deux termes est prépondérant. Nous sommes obligés dans ce cas là de garder les deux termes et
de ne pas simplier l'équation de Maxwell-Ampère.
La jauge de Lorenz reste inchangée.
Dans une situation quasi-stationnaire où les charges sont prépondérentes sur les courants, le
champ éclectique se calcule de la même manière qu'en électrostatique. Dans les équations de
Maxwell on néglige le terme de la dérivée temporel du champ magnétique dans Maxwell-Faraday.
Dans la force de Lorentz, le terme électrique domine. Par contre dans l'équation de Maxwell
Ampère pour le calcul du champ magnétique le courant de déplacement n'est pas négligeable.
L'énergie est stocké principalement sous forme électrique.
5.3.3 Régime quasi-stationnaire magnétique ( 1, ξ 1)

Plaçons nous dans le cas où :
L cA j̃
=1 et = =ξ1 (5.60)
τc V ρ̃c
Ceci signie que les courants dominent sur les charges comme dans le cas d'un conducteur avec électro-
neutralité.
Prenons l'équation de Maxwell-Ampère en remplaçant le champ électrique :
∂2A

1 ∂ ∂A 1 ∂V
∇ ∧ B = µ0 j + 2 −∇V − = µ0 j + 2 −∇ − (5.61)
c ∂t ∂t c ∂t ∂t2
5.3. DÉMARCHE RIGOUREUSE DE L'ARQS 65
Comparons alors le terme de courant avec celui du potentiel vecteur :
||µ j|| µ j̃ µ0 j̃ 4πLc2 τ 2 1 4πL3 1

0 ' A0 = = = 2 ∼ 2 1 (5.62)
1 ∂ 2 A
c2 τ 2
µ0 Ωj̃ Ω Ω
c2 ∂t2 4πLc2 τ 2
Ainsi le terme de variation temporelle du potentiel vecteur est négligeable. Comparons maintenant le
second terme :
||µ0 j|| dj̃τ j̃ dτ L2 c ξ dL2
' 4π0 L = = 1 (5.63)
1 ∂(∇V )
c2 ∂t 0 ρ̃Ω ρ̃c LΩ Ω
En eet le premier terme est grand devant 1 et le second terme est un nombre ni comme rapport de 2
volumes nis. On se retrouve donc avec un courant de déplacement négligeable et donc :
∇ ∧ B ' µ0 j ⇒ ∇.j = 0 (5.64)
La conservation de la charge se traduit alors par un courant à ux conservatif dans ce cas là. On aurait
pu également le montrer en considérant que :
||∇.j|| j̃ cτ j̃ cτ ξ
' > = 1 (5.65)
∂t ρ ρ̃c d ρ̃c L
On voit alors que on peut négliger la variation de ρ.
D'un autre côté les trois autres équations inchangés. Le champ magnétique est alors déterminé comme en
magnétostatique avec les mêmes équations.
Considérons maintenant la Jauge de Lorenz et regardons l'ordre de grandeur :
||∇.A|| µ0 j̃Ω 4πdτ µ0 j̃Ω 4πdτ ξ

1 ' 2
> 2
= 1 (5.66)
∂V
c2 t
4π d µ0 ρ̃Ω 4π L µ0 ρ̃Ω
Ainsi la Jauge de Lorenz se transforme en jauge de Coulomb pour l'ARQS magnétique. De plus on voit que
si nous sommes en ARQS électrique, on ne peut pas conclure sur quel terme peut être grand devant l'autre
et on garde alors la Jauge de Lorenz. Cette méthode n'est pas très rigoureuse car nous concluons qu'un
terme est négligeable devant un terme qu'on prend à 0 par la suite. Pour montrer plus rigoureusement
que la Jauge de Lorentz se transforme en Jauge de Coulomb, il faut reprendre les équations de Maxwell
en ARQS magnétique et refaire les calculs de magnétostatique.
On peut également montrer que :
||∇V ||
∂A ' ξ (5.67)

∂t
et donc qu'il n'y a pas de partie dominante entre les deux.
Si on compare maintenant les deux champs électrique et magnétique :
cB cµ0 j̃d j̃
' 0 = =ξ1 (5.68)
E ρ̃d ρ̃c
Donc le champ magnétique est prépondérant sur le champ électrique dans la force de Lorentz qui s'écrit
alors :
F = qv ∧ B (5.69)
Et en terme d'énergie on a directement :
||B||2 1
0 ||E||2 (5.70)
2µ0 2
Presque toute l'énergie est stocké sous forme magnétique.
Lorsque les courants sont prépondérants sur les charges alors le champ magnétique est prépon-
dérant et est de type statique. On peut négliger le courant de déplacement dans l'équation de
Maxwell-Ampère. Le terme de couplage dans Maxwell Faraday est toujours là et sera responsable
de l'induction. L'énergie est stockée sous forme magnétique.
5.3.4 Remarques supplémentaires

Quand on parle d'ARQS traditionnellement on parle d'ARQS magnétique où on néglige le courant de
déplacement. Mais cela ne rend pas compte de toutes les situations (seulement de 90%). Il est nécessaire
de distinguer les deux types d'ARQS et on peut dire que :
• ARQS magnétique : c'est l'ARQS des bobines. C'est celle qui est faite pour l'induction. Elle s'applique
à la physique des bobines en électrocinétique à basse fréquence.
• ARQS électrique : c'est l'ARQS des condensateurs. En électrocinétique à basse fréquence. Ce régime
décrit la physique des condensateurs qui stockent d'importantes charges tout en étant isolant, ce qui
assure de faibles courants.
Il est également possible de faire une analogie de l'ARQS magnétique avec la mécanique des uides
où le courant j a pour homologue la vitesse v en mécanique des uides et ainsi la nullité de leur
divergence respective indique une non accumulation de charge ou l'incompressibilité du uide.
La Table 5.1 récapitule tous les résultats et les équations simpliées dans les deux types d'ARQS.
Type d'ARQS Magnétique Électrique

Courants/Charges Les courants dominent Les charges dominent
L L
= 1 = 1
Hypothèses τc τc
cA j̃ cA j̃
= =ξ1 = =ξ1
V ρ̃c V ρ̃c
ρ ρ
∇.E = ∇.E =
0 0
Équations de Maxwell ∇.B = 0 ∇.B = 0
∂B ∇∧E = 0
∇∧E = −
∂t 1 ∂E
∇ ∧ B = µ0 j ∇ ∧ B = µ0 j +
c2 ∂t
Jauge Jauge de Coulomb : ∇.A = 0 Jauge de Lorenz : ∇.A + 1 ∂V

c2 ∂t
=0
Conservation de la charge Courants à ux conservatifs : ∇.j = 0 Pas de changement : ∇.j + ∂ρ

∂t =0
B2 0 E2
Énergie Forme magnétique : 2µ0 Forme électrique : 2
Électrocinétique Bobines à basses fréquences Condensateurs à basses fréquences
Table 5.1 Récapitulatif des deux types d'ARQS.

Chapitre 6
Conducteurs
Ce chapitre a pour but de montrer comment les champs électromagnétiques se comportent dans les
conducteurs électriques. Il sera nécessaire de modéliser la conduction électrique dans ceux-ci puis de
distinguer diérentes gammes de fréquences où les champs ne vont pas avoir le même comportement.
Certaines applications des conducteurs seront abordé, notamment la notion de magnétorésistance et d'eet
Hall.
6.1 Régime statique - Conducteurs à l'équilibre
On parle de conducteur à l'équilibre si les charges en son sein n'ont pas de mouvement d'en-
semble. Dans ce cas là, seul reste leur mouvement désordonné à l'échelle microscopique. En par-
ticulier, on suppose que le conducteur est homogène et possède une température uniforme.
Les conséquences de cette dénition sont :

• La densité locale de charges totale (électrons + ions du réseau) à l'intérieur du conducteur est nulle :
ρint = 0
• Par conséquent, dans le cas où on ne impose pas de champ électrique permanent à l'intérieur (ou
une diérence de potentiel), le champ électrique est nul à l'intérieur du conducteur : Eint = 0
• Le fait que le champ soit nul à l'intérieur du conducteur impose également une condition sur le
potentiel qui est constant partout à l'intérieur du conducteur : Vint = cte
• On en déduit également qu'au voisinage de la surface, le champ électrique est normal au conducteur.
Tout ceci n'est valable que pour l'intérieur du conducteur c'est-à-dire loin de la surface en principe. En
eet, dans ces conditions le conducteur porte une charge surfacique σ liée au champ perpendiculaire à la
surface. Proche de la surface le champ est donnée par le théorème de Coulomb :
σ
E= n (6.1)
0
On retrouve alors la discontinuité du champ électrique au passage d'une surface chargée.
6.2 Régime dynamique - Conduction

La dénition d'un conducteur au sens le plus large englobe bien entendu les métaux et les semi-
conducteurs mais également les solutions ioniques ainsi que les plasmas. Dans les deux derniers cas, les
porteurs (mobiles) de charges sont les ions en solution et les ions et les électrons dans le plasma.
67
68 CHAPITRE 6. CONDUCTEURS
Le temps de retour typique à l'équilibre lorsqu'il y a une inhomogénité de charge est de l'ordre de
10−14 s pour un bon conducteur (au mieux le temps de relaxation devient beaucoup plus long lorsque la
conductivité diminue). En conséquence la présence d'un courant dont le temps caractéristique d'évolution
est long devant ce temps de relaxation τrelax nous permet de dire que la densité de charges volumique reste
nul en très bonne approximation.
Modélisons le conduction électrique dans les conducteurs. La densité de courant s'écrit :

X
j= ρi vi (6.2)
i
avec ρi la densité des porteurs de charges de type i et vi leur vitesse.

On peut dire que dans un milieu isotrope et en l'absence de champ magnétique la vitesse des porteurs de
charges est, en régime linéaire, proportionnel au champ électrique E :
vi = µi (||E||)E (6.3)
où µi est la mobilité des porteurs du type i qui dépend à priori de la norme du champ électrique et donc
on a : !
X
j= ρi µi (||E||) E (6.4)
i
En très bonne approximation dans la limite des champs peu intenses, la mobilité est indépendante
de l'intensité du champ électrique et on a donc :
!
X
j= ρi µi E = σE (6.5)
i
On a donc, sous ces hypothèses, une proportionnalité entre le courant j et le champ électrique E.
Le coecient de proportionnalité σ s'appelle la conductivité du matériau et ne dépend que du
matériau (type de porteurs de charges et leur mobilité). Cette relation est la loi d'Ohm locale
car son intégration sur un conducteur cylindrique redonne la loi d'Ohm U = RI .
Dans le cas où on se place dans un milieu anisotrope mais toujours dans le régime linéaire, on trouve la
même expression mais σ devient une matrice qu'on peut diagonaliser avec des coecients diérents sur la
diagonale caractérisant l'anisotropie.
Cette relation nous informe juste sur la proportionnalité entre le courant et le champ électrique mais
non pas sur le comportement en fréquence de la conductivité et donc de la réponse du conducteur. On se
propose ici de développer un modèle microscopique pour en sortir une dépendance en la fréquence.
Modèle microscopique
Le modèle de Drude 1 est le modèle de conduction pour les métaux le plus simple. Elle a été proposé par P.
Drude en 1900 et suppose que les électrons, responsable de la conduction électrique, se comporte comme
un gaz. La théorie cinétique des gaz suppose que les molécules de gaz sont des sphères dures se déplaçant
en ligne droite jusqu'à rentrer en collision les unes avec les autres (ou les parois mais ce n'est pas le cas
des électrons). Le temps d'une seule collision est négligeable et les forces qui agissent lors des collisions
sont les seules forces prises en compte. Dans le cas des électrons dans un gaz on a comme hypothèses :
1. Voir le Cap prépa ou le Ashcroft de physique des solides.

6.2. RÉGIME DYNAMIQUE - CONDUCTION 69
1. Entre deux collisions toute interaction avec un autre électron ou un ion est négligée. Donc en absence
de champ électromagnétique externe, les électrons se déplacent en mouvement rectiligne uniforme.
Si on applique un champ extérieur, alors la dynamique est donné par le principe fondamentale
de la dynamique avec ces champs. On appelle cela approximation des électrons indépendants pour
négliger l'interaction électron-électron et approximation des électrons libres pour négliger l'interaction
électron-ion.
2. Les collisions sont des événements instantanés qui change de manière abrupte la vitesse des électrons.
Drude supposait que les électrons rebondissait sur le c÷ur des atomes, donc sur les ions. Ceci se
rapproche de la réalité car la diusion électron-électron est loin d'être prépondérante sauf dans
certains cas.
3. On suppose qu'un électron prend part à une collision avec une probabilité de 1/τ . τ est appelé
temps de relaxation, temps de collision ou temps de libre parcours moyen. C'est lui qui va piloter le
mécanisme de conduction électrique. Dans le modèle le plus simple, ce temps est indépendant de la
position et de la vitesse de l'électron car on prend un temps moyenné.
4. On suppose que les électrons établissent un équilibre thermique avec leur entourage uniquement avec
les collisions. On a un équilibre thermodynamique local : immédiatement après chaque collision,
l'électron émerge avec une vitesse indépendante de sa vitesse immédiatement avant la collision, mais
ayant une direction aléatoire et une vitesse correspondant à la température dominante au lieu où la
collision s'est produite.
On peut trouver une conductivité dans le cas stationnaire, mais le calcul en dynamique permet de directe-
ment le retrouver. Supposons alors qu'on applique un champ électrique sinusoïdale sur notre conducteur
pour voir la réponse en fréquence. Pour trouver une dépendance de la conductivité en fonction de la
pulsation du champ électrique σ(ω), on applique la modélisation microscopique que l'on vient d'exposer.
Pour cela il faut appliquer le PFD aux électrons mais comment modéliser la force de collision ? Soit p la
quantité de mouvement par électron. Soit p(t) la quantité de mouvement à l'instant. Calculons la quantité
de mouvement à un instant t + dt soit p(t + dt). La probabilité de collision est dt/τ et la probabilité que
rien ne se passe est (1 − dt/τ ). S'il n'y a pas de collision, l'électron est soumis à la force électromagnétique
qu'on appellera f (t) et acquiere donc une quantité de mouvement supplémenatire f (t)dt + O(dt)2 . La
contribution des électrons qui ont subi une collision dans l'intervalle dt est d'ordre (dt)2 et ne contribue
pas au premier ordre. Donc :

dt
p(t + dt) = 1− [p(t) + f (t)dt + O(dt)2 ] (6.6)
τ
dt
= p(t) − p(t) + f (t)dt + O(dt)2 (6.7)
τ
Donc la variation de quantité de mouvement s'écrit :
dt
p(t + dt) − p(t) = p(t) + f (t)dt + O(dt)2 (6.8)
τ
et donc :
dp(t) p(t)
=− + f (t) (6.9)
dt τ
Cela montre simplement que l'eet des collisions individuelles des électrons est d'introduire un terme
d'amortissement dans l'équation du mouvement de la quantité de mouvement par électron.
Prenons alors un champ électrique sinusoïdal en notation complexe. Le PFD s'écrit avec cette nouvelle
force :
me
me v̇ = −eE − v (6.10)
τ
avec E = E0 e−iωt et v = v0 e−iωt on trouve :

me
−me iωv0 = −eE0 − v0 (6.11)
τ
On trouve alors :
−eE0 /me
v0 = (6.12)
1/τ − iω
et donc :
ne2 τ /me
j0 = −nev0 = E0 (6.13)
1 − iωτ
On trouve alors la réponse du courant dans un conduction sous l'action d'un champ sinusoïdal
à la pulsation ω :
σ0 ne2 τ
j0 = E0 avec σ0 = (6.14)
1 − iωτ me
où σ0 est la valeur de la conductivité en statique.
On voit alors que tant que ωτ 1, la conductivité est proche de la conductivité statique (ce qui est le
cas jusqu'aux fréquences optiques pour un bon conducteur). Dans le régime opposé ωτ 1 on a une
conductivité imaginaire pur et le résultat se rapproche au résultat obtenu dans un plasma :
ne2 ωp2
σ(ω) = i = i0 (6.15)
me ω ω
On dénit alors une pulsation caractéristique appelé pulsation plasma qui vaut :
ne2
ωp2 = (6.16)
0 me
On peut eectuer une mesure de conductivité statique sur un métal (59.6 106 S/m pour le cuivre) et nous
permet de trouver un temps τ ' 2.10−14 s et donc un libre parcours moyen l ' 380.10−10 m.
6.3 Électroneutralité d'un conducteur

Considérons un conducteur initialement neutre partout dans l'espace. Cela se traduit que pour tout
point M du conducteur on a ρ(M ) = 0. À t = 0 on suppose que l'on perturbe l'état d'équilibre initial du
système en imposant une charge ρ(M0 ) = ρ0 6= 0. Cet excès local de charge est, bien entendu, compensé
par un décit de charges dans le reste du conducteur pour garder la neutralité. La question est de savoir
au bout de combien de temps voit-on restauré l'état d'électroneutralité local ? La réponse est loin d'être
triviale et trois contraintes sont à prendre en compte :
1. La relation local de la charge :
∂ρ
∇.j + =0 (6.17)
∂t
qui montre que toute variation temporelle de ρ est la source d'une densité de ux de courant j.
2. La loi d'Ohm locale :
j = σ(ω)E (6.18)
6.3. ÉLECTRONEUTRALITÉ D'UN CONDUCTEUR 71
3. Les équations de Maxwell et notamment celle de Maxwell-Gauss :

ρ
∇.E = (6.19)
0
Pour connaitre le temps de relaxation du système vers un état électriquement neutre il faut établir
l'équation diérentielle satisfaite par la densité de charge avec la condition initiale ρ(M0 ) = ρ0 . On exploite
ici l'équivalence entre les représentations temporelle et fréquentielle. Pour cela, on cherche le comportement
forcé d'un conducteur soumis à un champ électrique sinusoïdal de pulsation ω :
E(M, t) = E(M, ω)e−iωt (6.20)
La loi d'Ohm locale injectée dans la conservation de la charge donne :
∇.σE = iωρ ⇒ σ∇.E = iωρ (6.21)
On déduit en utilisant l'équation de Maxwell-Gauss que :
ρ
σ = iωρ (6.22)
0
Or on a vu que σ = σ0
1−iωτ donc :
σ0
−τ ω 2 ρ − iωρ + ρ=0 (6.23)
0
En repassant dans le domaine temporel on aboutit à :
∂ 2 ρ 1 ∂ρ ne2
+ + ρ=0 (6.24)
∂t2 τ ∂t me 0
On trouve alors une équation diérentielle sur la densité de charge du métal :

s
2
∂ ρ ωp ∂ρ ne2
+ + ω p ρ = 0 ω p = Q = ωp τ (6.25)
∂t2 Q ∂t me 0
On reconnaît alors l'équation caractéristique d'un oscillateur libre amorti de pulsation propre ωp
(pulsation plasma du métal) et de facteur de qualité Q. Avec des grandeurs typiques de métaux
on trouve que τ ∼ 10−14 s et un facteur de qualité Q ∼ 100. On en déduit donc que :
La relaxation vers l'état d'électroneutralité est réalisée au terme d'un régime pseudo-
périodique à la pulsation ωp ∼ 1016 rad.s−1 sur un temps caractéristique τ ∼ 10−14 s.
Ainsi tant que ωτ 1 on peut considérer que le conducteur est localement neutre
en tout temps.
Remarque : Cette démonstration suppose qu'on a déjà fait le modèle microscopique pour trouver la dé-
pendance en fréquence. On peut faire la démonstration en supposant qu'on a uniquement une conductivité
statique qui ne dépende pas de la fréquence et on trouve ainsi un temps de relaxation de τrelax = 10−18 s.
Cependant comme la loi d'Ohm locale n'est valable que si la durée caractéristique de variation du champ
électrique est bien supérieure à 10−14 s, on retrouve au nal le temps caractéristique τ .
6.4 Retour sur l'ARQS

On a vu que l'ARQS consistait à négliger le temps de propagation du champ électromagnétique à l'échelle
du système. On a vu que l'ARQS n'était valable que à basses fréquences mais ici on peut déterminer une
fréquence caractéristique qui délimite le régime de l'ARQS.
Dans un conducteur, ce sont les courants qui vont dominer par rapport aux charges et on se retrouve
dans l'ARQS magnétique. Dans ce cas là, le courant de déplacement est négligeable dans l'équation de
Maxwell-Ampère. Il faut donc trouver une condition sur la pulsation telle que :

∂E σ0 σ0
σ(ω)E 0 ⇒ |iω0 | ⇔ ω0 (6.26)
∂t |1 + iωτ | |1 + iωτ |
On peut passer l'équation au carré et donc la pulsation limite obéit à l'équation :
Ω σ02
Ω2 + − =0 Ω = ω2 (6.27)
τ 2 20 τ 2
On pourrait résoudre cette équation complètement an de trouver la pulsation caractéristique mais une
méthode plus simple est possible. Elle consiste à simplier l'étude et à rechercher directement la pulsation
critique ωc dans l'intervalle [1/τ ; +∞[. Si on trouve une solution ωc dans cet intervalle (et susamment
loin de 1/τ ), cette solution a de bonnes chances d'être proche de l'équation bicarrée précédente. Ainsi dans
l'approximation haute fréquence ω 1/τ on a :
σ0 σ0
' (6.28)
|1 + iωτ | ωτ
Et donc la condition précédente pour négliger le courant de déplacement devient :

r
σ0 σ0
0 ω ⇒ω = ωp (6.29)
ωτ 0 τ
On montre ainsi que la pulsation caractéristique recherchée s'indentie à la pulsation plasma ωp du

conducteur. On vérie a posteriori que cette pulsation est bien supérieure à 1/τ et se trouve bien dans
l'intervalle de départ. On a ωp ∼ 1016 rad.s−1 et 1/τ ∼ 1014 rad.s−1 .
Un conducteur se trouve en approximation des régimes quasi-stationnaires (ARQS) si
ω ωp ∼ 1016 rad.s−1 (6.30)
Ainsi, en résumé, on a trois gammes de fréquences où le conducteur va se comporter diéremment :
• Cas où ω 1014 rad.s−1 : l'électroneutralité et l'ARQS est vérié donc on a ρ = 0 ainsi

que le courant de déplacement négligeable.
• Cas où ω ∈ [1014 rad.s−1 ; 1016 rad.s−1 ] : L'électroneutralité n'est pas vériée de sorte que
l'équation de Maxwell-Gauss ne peut être simplié. L'hypothèse de l'ARQS peut théori-
quement être proposé. On notera cependant que la bande de fréquence de ce domaine est
très réduit.
• Cas où ω 1016 rad.s−1 : Dans ce régime, ni l'ARQS ni l'électroneutralité ne sont permises.
Aucune simplication dans équations de Maxwell sont possible.
6.5. ONDES ET MÉTAUX CONDUCTEURS : EFFET DE PEAU 73
6.5 Ondes et métaux conducteurs : eet de peau

A priori, le problème de la propagation d'une onde électromagnétique dans un métal est relativement
compliqué au vu des eets du champ sur les charges et du champ qu'elles créent en retour. On se limitera
alors à une description macroscopique du problème dans les cas assez simples : une Onde Plane Progressive
Harmonique en incidence normale sur un conducteur imparfait (conductivité nie) avant de discuter la
réexion d'une onde sur un métal parfait.
6.5.1 Métal imparfait

Prenons le cas d'un métal de conductivité nie σ qui occupe le demi-espace inni x > 0. L'onde incidente
normale est polarisée selon z , ne dépend que de la coordonnée x et est en notation complexe de la forme :
E = E0 f (x)eiωt ez (6.31)
L'ensemble des autres grandeurs seront donc de la même forme. Avec l'équation de Maxwell-Faraday on
trouve que :
−E0 f 0 (x)ey = −iωB (6.32)
1
B = E0 f 0 (x)eiωt ey (6.33)
iω
où la divergence de E et de B sont bien nulles.
Pour un conducteur, le courant de déplacement est négligeable devant le courant de conduction (ARQS).
Typiquement le rapport vaut environ 109 dans le domaine micro-ondes, et donc :
∇ ∧ B = µ0 j = µ0 σE (6.34)
où on va supposer que la conductivité σ est statique prend sa valeur statique σ0 pour éviter une dépendance
en ω supplémentaire. On a alors l'équation sur f :
1 00
f (x) = µ0 σ0 f (x) (6.35)
iω
f 00 (x) − iωµ0 σ0 f (x) = 0 (6.36)
Cette équation a pour solution une somme de deux exponentielles, une positive et une négative. Nous
éliminons l'exponentielle positive pour ne pas avoir un champ qui diverge car cela n'est pas physique. On
a également : r
µ0 σ 0 ω
(6.37)
p
iωµ0 σ0 = ±(1 + i)
2
On pose alors la quantité : r

2
δ= (6.38)
µ 0 σ0 ω
homogène à une longueur appelé épaisseur de peau
On a alors :
x x
E = E0 e− δ ei(ωt− δ ) ez (6.39)
On a alors un champ qui se propage dans le métal mais en étant atténué exponentiellement sur une
distance caractéristique δ . On obtient une onde amortie 2 . Cette épaisseur est d'autant plus faible que la
conductivité du matériau est grand et que la fréquence de l'onde est élevé.
2. Il ne s'agit pas d'une onde évanescente car il faut que la direction de l'atténuation et celle de propagation ne soit pas
la même.
L'eet de peau est d'autant plus marqué que la profondeur de peau est faible, c'est-à-dire que le conduc-
teur s'oppose ecacement à la pénétration du champ. L'expression de δ montre qu'il s'agit d'une fonction
décroissante de ω . Plus les champs varient rapidement dans le temps plus l'eet de peau est marqué. Il
s'agit là d'un eet inductif. Plus ω est grand, plus le champ électromagnétique varie rapidement dans le
temps. D'après la loi de Lenz, le milieu conducteur devient alors le siège d'une loi de modération visant à
s'opposer à cette variation de ux. Il y a création dans le conducteur d'un courant électrique et d'un champ
électrique induit se superposant au champ électrique incident et de sens opposé. Ce phénomène d'induction
justie l'atténuation du champ électrique dans le conducteur. La nature inductive de l'eet de peau ne doit
pas surprendre. Le phénomène d'induction est directement inscrit dans l'équation de Maxwell-Faraday qui
dans l'ARQS est la seule équation rendant compte du couplage du champ électromagnétique.
L'eet de peau saura d'autant plus marqué que le mécanisme d'induction sera ecace. La réponse
au conducteur aux variations du ux de B sera d'autant plus ecace que celui-ci présentera une bonne
aptitude à développer des courants induits à partir du champ électrique induit.
La puissance dissipée par eet Joule par unité de volume dans le conducteur est égale à :
dP
= j.E = σE2 (6.40)
dτ
Cela permet de calculer la puissance dissipée sur chaque tranche de l'onde.
6.5.2 Limite du métal parfait

Dans la limite d'un métal parfait, c'est-à-dire la conductivité qui tend vers l'inni, la longueur de peau
tend vers 0 et l'onde électromagnétique ne pénètre pas dans le conducteur. On peut également voir que
pour la puissance dissipée, si σ → ∞ alors il faut nécessairement que E = 0 pour éviter la divergence de
la dissipation.
De plus on peut montrer que au voisinage du conducteur on a :
σS
E= n B = µ0 jS ∧ n (6.41)
0
Le fait que le champ électromagnétique à l'intérieur du conducteur parfait soit nul lui confère des
propriétés intéressantes en réexion et donc des propriétés intéressantes pour du guidage des ondes élec-
tromagnétiques.
6.5.3 Force exercée sur le conducteur par le champ

On considère le cas d'une onde plane qui arrive sur une plaque métallique qu'on considère comme un
conducteur parfait. La question est de savoir quelle est la force exercée par l'onde incidente sur la plaque
métallique.
Comme précédemment, la plaque occupe le plan yOz . On a E selon ez et B selon ey Le champ va

engendrer des courants surfaciques qui sont dirigés selon ez . Ces courants subissent à leur tour l'action du
champ. On revient alors à une description volumique de la plaque à laquelle on donne une épaisseur a.
Dans cette couche on a :
B = By (x)ey (6.42)
En négligeant le courant de déplacement on a :
1 ∂By
j= ez (6.43)
µ0 ∂x
6.6. MAGNÉTORÉSISTANCE ET EFFET HALL 75
Un volume innitésimale dans la plaque est donc traversé par un courant :

1 ∂By
di = dxdy (6.44)
µ0 ∂x
qui se traduit par l'élément de courant :
1 ∂By
didl = dxdydzez (6.45)
µ0 ∂x
Cet élément subit alors une force :
dxdydz ∂By2
d2 F = didl ∧ B = − ex (6.46)
2µ0 ∂x
En intégrant sur l'épaisseur a de la plaque :
dydz 2
dF = − (By (x = a) − By2 (x = 0))ex (6.47)
2µ0
On a un champ mégnétique nul en x = a car pas de pénétration du champ et en x = 0 on a By (0) =
− 2Ec 0 cos(ωt) pour une onde plane. Et donc :
2E02
dF = cos2 (ωt)dydzex (6.48)
µ0 c2
La force est dirigée vers l'intérieur du métal, normale à sa surface et proportionnel à l'élément de surface :
on a une force de pression. La moyenne temporelle de cette pression s'écrit :
hdFi = P dS (6.49)
E02
P = = 0 E02 (6.50)
µ0 c2
qui est la pression de radiation.
Remarque : Cette expression permet de "sentir" intuitivement l'interprétation du vecteur de Poynting
comme quantité de mouvement véhiculée par le champ.
6.6 Magnétorésistance et eet Hall

6.6.1 Magnétorésistance
La question qu'on se pose dans ce paragraphe est de savoir ce qu'il se passe dans un conducteur si on
applique un champ magnétique dans celui-ci. On se passe ici en régime statique ou basse fréquence de tel
sorte que la conductivité ne dépend pas de la fréquence.
Dans le cas de l'application d'un champ magnétique supplémentaire dans le conducteur, la force de
Lorentz est modiée de cette façon :

j
−eE → −e(E + hvi ∧ B) = −e E − ∧B (6.51)
ne
Et donc dans la loi d'Ohm on peut également faire la substitution pour obtenir :

j
j = σ0 E − ∧B (6.52)
ne
Il faut remarquer que la champ magnétique est la somme du champ extérieur qu'on applique ainsi que le
champ propre induit qu'on négligera par la suite. Sous cette forme, la loi d'Ohm met en évidence qu'en
présence d'un champ magnétique dans un conducteur les lignes de champ électrique et de courant ne sont
pas confondues.
An de comprendre un peu cet eet, nous allons simplier un peu le problème en considérant un repère
cartésien et en plaçant le champ magnétique constant selon l'axe Oz : B = Bez . On peut alors écrire la
loi d'Ohm sous forme matricielle :
       
jx Ex jx 0
σ
jy  = σ0 Ey  − 0 jy  ∧  0  (6.53)
ne
jz Ez jz B
On trouve alors les trois composantes de j :

B
jx = σ0 Ex − jy
ne

B
jy = σ0 Ey + jx
ne
jz = σ0 Ez
Si on pose la pulsation cyclotron caractéristique ωcy = eB m et qu'on remplace la conductivité statique

ne2 τ
par son expression σ0 = m alors on trouve :
−ωcy τ
 
1
0
  
jx 1+(ωcy τ )2 1+(ωcy τ )2 Ex
ω τ
 jy  = σ0  1+(ωcycy τ )2 1+(ω1cy τ )2 0   Ey  (6.54)

jz 0 0 1 Ez
On obtient alors la loi d'Ohm matricielle. Les composantes de la densité de courant sont des fonctions
linéaires de celles du champ électrique mais il n'y a plus proportionnalité.
Si à présent on simplie encore le problème en imposant uniquement un champ électrique (une diérence
de potentiel) selon une seule direction, disons selon l'axe Ox. Alors à ce moment là, la loi d'Ohm s'écrira :
σ0 σ0 eB
jx = E < σ0 avec ωcy = (6.55)
1 + (ωcy τ )2 1 + (ωcy τ )2 m
La conductivité diminue lors de la présence d'un champ magnétique non aligné avec le champ
électrique. On parle alors de magnétorésistance.
On remarquera que plus le champ B est intense plus la conductivité diminue. Il en est de même pour
le rapport e/m. Ainsi, plus il est aisé de mettre en mouvement des charges (soit une charge de particule
plus grande ou une masse de particule plus faible) plus l'eet du champ magnétique sera important. Les
particules au lieu de se propager en ligne droite voient leur trajectoire courbée par le champ magnétique.
6.6.2 Eet Hall

L'eet Hall n'est qu'un cas particulier de la magnétorésistance mais dans un cas conné. On considère
un barreau métallique de taille a selon la direction y . Le barreau a une taille b selon l'axe z et un champ
magnétique uniforme est toujours placé dans cette direction. On fait alors parcourir un courant I selon l'axe
x et ainsi à cause du champ magnétique il apparait une force sur un porteur de charges fm = −e hvi ∧ B.
6.7. PLASMA 77
On choisit le barreau susamment étroit selon l'axe y (donc a faible) pour que les charges qui voit leur
trajectoire courbée rencontrent directement le bord du conducteur et il apparait une accumulation de
charges négative sur une face (et donc un défaut de charges dans l'autre). Il se forme donc un champ
électrique dit de Hall noté EH . En régime stationnaire, le nombre de charges sur les bords ne varie pas et
donc le force due au champ de Hall compense exactement la force magnétique et on a alors :
1
EH = − hvi ∧ B = j∧B (6.56)
ne
où on note RH = 1
ne la résistance de Hall.
La loi d'Ohm matricielle permet de retrouver directement ce résultat en prenant jy = jz = 0 soit :

Ex − ωcy τ Ey
jx = σ0 2 τ2
1 + ωcy
0 = ωcy τ Ex + Ey
0 = Ez
Et donc on retrouve la loi d'Ohm :
jx = σ0 Ex (6.57)
ainsi que le champ de Hall :
ωcy τ B
Ey = − jx = − jx (6.58)
σ0 ne
Le champ de Hall est responsable de la tension de Hall qui apparaît entre les faces inférieure et
supérieure du barreau :
IB
∆V = (6.59)
neb
Toute le raisonnement est transposable pour des semi-conducteurs où les porteurs de charges peuvent
être positifs mais sont surtout en nombre plus important car les semi-conducteur sont dans ce cas là dopé.
On s'en sert pour mesurer des champs magnétiques une fois que la sonde a été calibrée.
6.7 Plasma
Un autre type de conducteur qu'on peut étudier est le plasma. Il s'agit d'un ensemble macroscopiquement
neutre de charges positives (ions positifs de masse M et de charge +e) et d'électrons (de masse m et de
charges −e) de densité commune n. On peut alors étudier la propagation d'une onde plane progressive
monochromatique (OPPH) dans un tel milieu. On supposera l'onde polarisée rectilignement.
Pour établir la loi de conduction on applique PFD aux électrons et aux ions. La seule force non négligeable
appliquée sur un électron ou un ion est la force de Lorentz. Si v+ et v− sont les vitesses d'un ion et d'un
électron alors les équations sont :
dv− dv+
m = −e(E + v− ∧ B) M = e(E + v+ ∧ B) (6.60)
dt dt
Dans le plasma les normes des champs sont liés par la relation E = vϕ B avec vϕ la vitesse de phase. Or
comme les vitesses des ions et des électrons sont très inférieurs à la vitesse de phase (hypothèse à vérier
a postériori) alors la composante magnétique de la force de Lorentz est négligeable et donc en se plaçant
en notation complexe on trouve les équations :
ie ie
v− = − E v+ = E (6.61)
mω Mω
La densité de courant s'écrit :

ne2

1 1
j = ρ− v− + ρ+ v+ = ne(v+ − v− ) = i + E (6.62)
ω m M
2
On obtient alors une loi d'Ohm complexe avec la conductivité σ = i neω m 1 1
. La conductivité est

+M
imaginaire pur dans ce cas là et le courant est en retard de π/2 par rapport au champ électrique. Comme
la masse d'un proton est 1836 fois celle d'un électron on peut négliger 1/M devant 1/m et donc en posant
la pulsation plasma s
ne2
ωp = (6.63)
0 m
on obtient comme conductivité :
ωp2
σ(ω) = i0 (6.64)
ω
D'après la conservation de la charge, l'équation de Maxwell-Gauss ainsi que la loi d'Ohm on obtient :
ρ
iσk.E − iωρ = 0 ik.E = (6.65)
0
soit :
(ωp2 − ω 2 )ρ = 0 (6.66)
Si la pulsation est diérente de la pulsation plasma alors le plasma est neutre, sinon il n'y a pas de
raison d'avoir électroneutralité mais on verra que dans ce cas le vecteur d'onde est nul et qu'il n'y a pas
de propagation possible.
Établissons à présent la relation de dispersion à partir des équations de Maxwell en se plaçant en notation
complexe E = E0 ei(k.r−ωt) :
ik ∧ E = iωB
ik.B = 0
ik.E = 0
iω
ik ∧ B = µ0 σ − 2 E
c
Ainsi on a : !
iω ωp2
2
ik E = ik ∧ (k ∧ E) = 2 −ω E (6.67)
c ω
Donc :
ω 2 − ωp2
k2 = (6.68)
c2
Ainsi en fonction de la valeur de ω on peut avoir un k 2 positif ou négatif et donc un nombre d'onde
réel ou imaginaire. Un k réel permettra une propagation d'une onde plane normalement, on parle alors de
bande passante (si ω > ωp ). On obtient une onde qui se propage et tel que :
1q 2 c
k=± ω − ωp2 vϕ = ± q (6.69)
c ωp2
1− ω2
Dans le cas où k est immaginaire pur alors le champ électrique décroit exponentiellement et on a plus de
prorogation d'onde. Dans cette gamme de pulsation appelée bande interdite (si ω < ωp ) on obtient une
6.7. PLASMA 79
onde évanescente et en réalité comme il ne peut avoir propagation dans le plasma, l'onde incidente est
rééchie.
On remarquera que dans la bande passante on a vϕ > c donc notre hypothèse au départ sur la vitesse
des électrons faible devant vϕ est justié si les électrons sont non-relativistes. Ceci n'est pas anormal car
même si c constitue la vitesse maximale de propagation des interactions et de l'information. Pour une onde
plane qui occupe tout l'espace un observateur la voit passer indéniment devant lui mais rien ne change au
cours du temps à part l'oscillation temporelle. Une telle onde est une idéalisation et ne transporte aucune
information et ne peut servir qu'à indiquer le début ou la n d'un message. Il faudra utiliser la vitesse de
groupe caractériser la vitesse de propagation de l'énergie et donc de l'information.
Une application possible est la propagation dans la ionosphère qui est un plasma de pulsation plasma
de l'ordre de 10 MHz. Elle constitue une couche supérieure de l'atmosphère. Ainsi en se plaçant à des
pulsations inférieures à la pulsation plasma on aura réexion sur celle ci et elle peut se comporter comme
un guide d'onde pour des ondes type radio qui vont se rééchir sur la ionosphère puis sur la Terre. Il
faut en revanche faire attention lorsqu'on veut envoyer des signaux dans l'espace à utiliser des ondes de
pulsations plus grandes que la pulsation plasma.
Chapitre 7
Guidage d'ondes
Le guidage d'ondes électromagnétiques est omniprésente dans le transport de l'information. Nous verrons
dans ce chapitre les bases du guidage d'ondes avec l'utilisation de conducteurs "parfaits" pour leur réexion
totale.
7.1 Réexion sur un métal parfait : nécessité d'une onde rééchie

Reprenons le cas d'un onde plane arrivant sur un métal parfait occupant le demi-espace x > 0. On a :
ex ex
Ei = E0 cos(ωt − kx) Bi = ∧ Ei = ∧ E0 cos(ωt − kx) (7.1)
c c
L'onde rééchie est de la forme :
Er = E00 cos(ω 0 t + k 0 x + ϕ) (7.2)
avec une pulsation et un nombre d'onde a priori diérents, une phase supplémentaire et une propagation
selon les x décroissants.
En x = 0, la condition aux limites s'écrit :
Ei (0, t) + Er (0, t) = 0 ⇒ E0 cos(ωt) + E00 cos(ω 0 t + ϕ) = 0 ∀t (7.3)
Comme la condition aux limites est vrai pour tout temps alors on a nécessairement on a ω = ω 0 et donc
par la relation de dispersion k = k 0 . De plus, pour la phase :
ϕ=π E0 = E00 ou ϕ = 0 E0 = −E00 (7.4)
Donc :
ex ex
Er = −E0 cos(ωt + kx) Br = −
∧ Er = ∧ E0 cos(ωt + kx) (7.5)
c c
Au voisinage du métal on a B = 2Bi qui vérie bien les conditions aux limites.
Ainsi la réexion sur un métal parfait a pour eet un déphasage de π pour le champ électrique
et une simple réexion sans déphasage pour le champ magnétique.
L'onde résultante vaut donc :

E = E0 [cos(ωt − kx) − cos(ωt + kx)] (7.6)
Le passage en complexe permet d'écrire :

−iωt 2π
E = −2iE0 sin(kx)e E = 2E0 sin x sin(ωt) (7.7)
λ
81
82 CHAPITRE 7. GUIDAGE D'ONDES
De même on trouve :
2ex 2π
B= ∧ E0 cos x cos(ωt) (7.8)
c λ
Les champs sont de la forme f (x)g(t), ce sont donc des ondes stationnaires. Les n÷uds vérient :

2π λ
sin x =0 ⇒ x = −q q∈N (7.9)
λ 2
Les ventres vérient :

2π 1 λ
sin x = ±1 ⇒ x=− q+ q∈N (7.10)
λ 2 2
Le vecteur de Poynting s'écrit alors :

4π
Π= 0 cE02 sin(2ωt) sin x ex (7.11)
λ
La puissance transportée à travers S est alors :

4π
pt = 0 ScE02 sin(2ωt) sin x ex (7.12)
λ
Et la puissance moyenne vaut :

P = h||pt ||i = 0 (7.13)
On a également que la charge surfacique est nulle car le champ transverse est nul. Le courant surfacique
vaut :
jS = 20 cE0 cos(ωt) (7.14)
7.2 Réexion sous incidence oblique

On considère une OPPM polarisée rectilignement au plan d'incidence. Elle arrive sur une métal parfait
avec un vecteur k faisant un angle i avec la normale au métal.
Ei = E0 cos(ωt − k.r)ey Er = E00 cos(ω 0 t − k0 .r + ϕ)ey (7.15)
Au voisinage du métal, le champ est nul donc comme précédemment, on a :
ω = ω0 k.r = k0 .r ⇒ (k − k0 ).r = 0 (7.16)
Donc k − k0 est orthogonal au plan d'incidence et au plan métallique. Comme les deux vecteurs d'onde ont
la même norme d'après la relation de dispersion alors k0 est le symétrique de k par rapport à la normale
au plan métallique et au plan d'incidence. On retrouve la loi de la réexion. La conclusion sur la phase
est la même que précédemment. On trouve alors :
E = −2E0 sin(kx cos i) sin(ωt − kz sin i)ey (7.17)
By = 0 (7.18)
2E0
Bx = sin i sin(kx cos i) sin(ωt − kz sin i) (7.19)
c
2E0
Bz = − cos i cos(kx cos i) cos(ωt − kz sin i) (7.20)
c
7.3. GUIDE D'ONDES 83
En espace limité, l'onde résultante est non plane, progressive et se propage selon z . On note :
ω ω c
K = k sin iez vϕ = = = >0 (7.21)
K k sin i sin i
On a une onde stationnaire selon x avec des plan n÷udaux :
λ
x=q q∈N (7.22)
2 cos i
Et des ventres tel que :

1 λ
x= q+ q∈N (7.23)
2 2 cos i
On montre toujours que la charge surfacique est nulle et le courant surfacique vaut :
jS = 20 cEi cos iy (7.24)
7.3 Guide d'ondes

7.3.1 Plan-plan
Deux plans métalliques, parfaits, illimités, parallèles, distants de a. Les plans étant innis selon les
directions y et z on va se placer dans un plan perpendiculaire aux deux plans du guide d'onde et propager
une onde selon z . On place un plan en x = 0 et un plan en x = a.
Modélisation et solutions
On envoie une OPPM dans ce guide d'ondes de pulsation ω et de vecteur d'onde K faisant un angle i avec
les plans. L'onde résultante est propagative selon z avec un vecteur d'onde :
K = k sin i ez (7.25)
et on prend une onde est de la forme :
E = E(x) cos(ωt − Kz)ey (7.26)
qui est un mode transverse électrique (TE). Les modes transverses magnétiques fonctionnent de la même
façon mais avec le champ magnétique.
Les conditions aux limites vont donner la fonction E(x) : les plans métalliques doivent être des plans
nodaux pour E. La distance d'un plan nodal au plan x = 0 est :
λ λ
x=q q∈N ⇒ a=n n ∈ N∗ (7.27)
2 cos i 2 cos i
Pour chaque entier n on a un mode de propagation avec :
2π π
k= =n (7.28)
λ a cos i
avec le champ :
nπx
E = −2E0 sin cos(ωt − Kz)ey (7.29)
a
Relations de dispersion
Pour trouver les relations de dispersion pour chaque mode, injectons la solution dans l'équation de
d'Alembert. On trouve :
ω2 nπ 2
= K 2
+ (7.30)
c2 a
On obtient donc un ltre passe haut de pulsation de coupure ωcn = nπc
a .
• Si ω < ωcn alors K est immaginaire et on a une onde évanescente.
• Si ω > ωcn alors on a propagation à la vitesse de phase :
c
vϕ = q >c vg vϕ = c2 (7.31)
2
ωcn
1− ω2
Les modes transverse électrique (TE) et transverse magnétique (TM) de rang n, noté TEn et TMn sont
entièrement caractérisés par un unique entier n non nul.
On voit donc que l'onde est loin d'être plane au sens strict du terme. En eet son amplitude varie
sinusoïdalement dans les plans x = cte qui seraient les plans d'ondes naturels. Il existe des n÷uds et des
ventres dans les plans équiphases. Pour une pulsation ω donnée, il existe plusieurs modes de propagations,
caractérisés par n, qui vérient la relation de dispersion.
On remarque que pour une valeur de ω donné tout les n ne sont pas accessibles. Une onde de fréquence
donnée ne peut se propager que dans un nombre ni de modes.
Calculons maintenant le champ magnétique associé. Pour cela il faut revenir aux équations de Maxwell.
En eet les relations de structures ne sont pas valable dans le cas présent, l'onde n'étant pas une onde
plane ! En utilisant l'équation de Maxwell-Faraday, il vient que :
∂B
= ∇∧E
∂t
∂Ey ∂Ey
= ex − ez
∂x ∂x
nπ nπx nπx
= 2E0 cos cos(ωt − Kz)ez + K2E0 sin sin(ωt − Kz)ex
a a a
En intégrant par rapport au temps :
nπ nπx K nπx
B = E0,n cos sin(ωt − Kz)ez − E0,n sin cos(ωt − Kz)ex (7.32)
aω a ω a
L'onde se propage dans le vide mais pour autant elle n'est pas transverse ! Ce sont les conditions aux
limites qui modient considérablement la structure de l'onde. Celle-ci ne dépend donc pas uniquement du
milieu de propagation.
Décomposition en ondes planes

Remarquons que l'onde TEn précédente peut se mettre sous la forme d'une superposition de deux ondes
planes :
nπx
En = −2E0 sin cos(ωt − Kz)ey
h a nπx nπx i
= −E0 sin ωt − Kz + + sin ωt − Kz − ey
a a
= −E0 [sin(ωt − k1,n .r) + sin(ωt − k2,n .r)] ey
= E1,n + E2,n
On y reconnaît la somme de deux vraies OOPH de vecteurs d'onde k1,n et k2,n de même composante
selon z mais de composante opposée selon z . Il n'est pas étonnant que la solution se décompose sous forme
d'OOPH (puisque ce sont une base des solutions de l'équation de d'Alembert). Nous verrons par la suite
pourquoi cette décomposition est intéressante. Ces ondes par contre vérient la relation de structure et la
relation des ondes électromagnétiques dans le vide.
Mode TM et mode TEM

Les calculs sont analogues à ceux menés dans la partie précédente dans le cas du mode TM et donne la
même relation de dispersion. Un point important est l'existence d'une onde supplémentaire qui n'existe pas
dans le cas des ondes TE. En eet dans le cas des ondes TM il existe une onde transverse électromagnétique
(TEM). En eet, le groupe TE ne contient aucun mode TEM, puisque n ∈ N∗ et Bz est non nul dès que
n est non nul. C'est diérent pour le groupe TM, puisque :
B = B0 ei(ωt−Kz) ey (7.33)
E = cB0 ei(ωt−Kz) ex (7.34)

dénit bien une onde TEM qui est bien du groupe TM et qui de plus compatible avec les conditions aux
limites. En eet ce conditions impliquent la nullité de Bx , Ey et Ez mais rien sur les autres composantes en
raison de l'existence de charges et de courants surfaciques. En un sens, les ondes TM dans cette géométrie
sont moins contraintes que les ondes TE. Cette onde est plane et vérie la relation de dispersion ω = Kc.
Il existe un certain nombre de résultats sur les ondes TEM :

• Un guide d'onde électromagnétique ne peut propager d'onde TEM que si il est constitué d'au moins
deux conducteurs diérents (nous aurons l'occasion de revenir sur ce résultat).
• Un mode TEM vérie toujours la relation de dispersion des ondes planes dans le milieu illimité,
même si cette onde n'est pas une onde plane au sens strict.
• Par opposition aux modes TE et TM, lorsque le mode TEM existe il est unique.
Propagation des modes et dispersion

Pour le mode TEn comme pour le mode TMn nous avons obtenu la relation de dispersion :
ω 2 nπ 2
K2 = 2 − (7.35)
c a
C'est une relation de dispersion de Klein-Gordon qui dépend de l'indice n du mode. On est donc amené à
identier :
pπc
ωc,n = (7.36)
a
qui est la fréquence de coupure du mode n. En réécrivant cette relation de dispersion, on obtient :
cK 2 = ω 2 − ωc,n
2
(7.37)
s
cK ω2
= 2 −n
2 (7.38)
ωc,1 ωc,1
comme tracé sur la Figure 7.1.
Si on envoie une onde de pulsation ω celle-ci va eectivement se déformer. En eet, l'onde va être trans-
portée par diérents modes qui ont des célérités diérentes. On peut faire la construction comme sur la
Figure 7.2 en exprimant les vitesses de phase et de groupe de chaque mode :
r ω 2
c c,n
vϕ,n = q et vg,n = c 1 − (7.39)
ωc,n 2 ω
1− ω
14
12
10
0
0 2 4 6 8 10 12 14
Figure 7.1 Relation de dispersion du guide d'ondes électromagnétiques plan-plan. Pour une pulsation
donnée, on va avoir plusieurs mode de propagation possible.
Pour un paquet d'onde à spectre fréquentiel étroit, les vitesses de propagation accessibles sont fortement
élargies par la dispersion de mode, bien plus que la dispersion intramodale. Cela conduit à des fortes
déformations voir même des dédoublements du signal. Les dimensions des guides d'ondes sont adaptés à
la propagation de ces ondes pour s'assurer que ceux-ci restent monomodes.
On appelle dispersion intermodale ou dispersion de mode la diérence de vitesse eective de

propagation de deux ondes harmoniques de même pulsation mais portées par des modes diérents.
Elle est due aux conditions aux limites imposées à l'onde.
La dispersion intramodale caractérise la diérence de vitesse eective de propagation entre
deux harmoniques portées par le même mode (n xé) mais de pulsation diérentes. Elle est due
aux propriétés intrinsèques du milieu (milieu dispersif) et aussi aux conditions limites.
La dispersion de mode est d'importance technologique car elle limite grandement l'information que peut
transporter un guide d'onde, notamment les bres optiques. Dans le cas illustré sur la Figure 7.3, on
préféra utiliser un guide monomode.
En conclusion, si on considère l'ensemble des modes TE, la dispersion dans le guide d'onde est imposé
par l'équation de Klein-Gordon, dupliquée autant de fois qu'il y a de modes. Si une OPPH du groupe TE
de pulsation ω est forcée à l'entrée du guide, sa propagation par le guide dépend de sa fréquence :
• Si ω < ωc,1 , tout les modes sont amortis et l'onde ne se propage pas dans le guide.
• Si ωc,1 < ω < ωc,2 l'onde ne se propage que dans son mode fondamental, avec dispersion intramodale.
La propagation se fait sans atténuation et le guide est monomode.
• Si ωc,n < ω < ωc,n+1 , l'onde est propagée par n modes, sans atténuation mais avec une double
dispersion, à la fois intermodale et intramodale. Le guide d'onde est multimode.
Pour qu'une OPPH puisse se propager dans le guide, il faut que :
πc
ω > ωc,1 =
a
2πc πc
>
λ a
λ < 2a
Figure 7.2 Pour une valeur de ω donnée, il existe plusieurs modes de propagation possibles (ici trois)
correspondants à des valeurs de k diérentes et donc à des valeurs de vitesses de phase et groupe diérentes.
ceci induit de la dispersion intermodale.
Figure 7.3 Illustration de la dispersion intermodale : dédoublement d'un paquet d'ondes dans un guide
pour les modes TE1 et TE2 avec ω = 23 ωc,1 .
On retrouve le critère géométrique sur la longueur d'onde de l'OPPH qui est très général. Le préfac-
teur dépend du guide d'onde et de sa géométrie. Attention cependant à ne pas interpréter ce critère
comme une longueur qui rentre ou pas dans le guide d'onde puisque λ correspond à la longueur
d'onde dans la direction de propagation.
7.3.2 Extension au guide rectangulaire

On ajoute deux plans de métal parfait parallèles distants de b. Une onde TE se propage selon z , soit E x ,
E y , B x , B y , B z les composantes complexes cherchées. On cherche Bz sous la forme :
B z = f (x)g(y)ei(Kz−ωt) (7.40)
On a les équations de Maxwell :
iω
∇.E = 0 ∇.B = 0 ∇ ∧ E = iωB ∇∧B=− E (7.41)
c
On trouve alors les relations :
K K
Bx = − E By = E (7.42)
ω y ω x
iω ∂Bz iω ∂Bz
− 2
Ex = = −iKyy − 2
Ey = − = iKBx (7.43)
c ∂y c ∂x
On pose :
ω2
kc2 = − K2 (7.44)
c2
d'où :
iω ∂Bz iω ∂Bz
Ex = Ey = − (7.45)
kc2 ∂y kc2 ∂x
Pour avoir une onde TEM il faut que Bz soit égale à 0 mais dans ce cas on a toutes les autres
composantes nulles et donc il n'y a pas de telle onde qui se propage.
Ce résultat est indépendant de la forme de la section droite du guide à condition qu'elle soit
délimitée par une courbe fermée ce qui n'est pas le cas du guide plan-plan où un mode TEM est
possible.
On a par ailleurs :
∂ 2 Bz ∂ 2 Bz f 00 g 00
+ + k 2
c zB = 0 ⇒ + + kc2 = 0 (7.46)
∂x2 ∂y 2 f g
Or, la somme d'une fonction de x et d'une fonction de y ne peut être constante quelque soit x et y sauf si
p, q ∈ R :
f 00 g 00
= −p2 = −q 2 p2 + q 2 = kc2 (7.47)
f g
Alors :
f (x) = A cos(px) + B sin(px) g(x) = C cos(qy) + D sin(qy) (7.48)
Donc :
Bz = (A cos(px) + B sin(px))(C cos(qy) + D sin(qy))ei(Kz−ωt) (7.49)
Conditions aux limites pour le champ électrique :
π
Ey (0 ou a, y, z, t) = 0 ∂x Bz |x=0,a = ∂x f |x=0,a = 0 ⇒ B=0 p=m m ∈ Z (7.50)
a
π
Ey (x, 0 ou b, z, t) = 0 ∂y Bz |y=0,b = ∂y g|y=0,b = 0 ⇒ D=0 q=n m ∈ Z (7.51)
b
Donc p et q sont quantiés et ne peuvent être nuls simultanément sinon les composantes transverses des
champs seraient nulles.
Chaque couple (n, m) détermine un mode de propagation noté T Emn de champs :
mπx nπy
Bz = B0 cos sin cos(ωt − kz) (7.52)
a b
ω nπ mπx nπy K
Ex = − 2 B0 cos sin sin(ωt − kz) By = Ex (7.53)
kc b a b ω
ω mπ mπx nπy K
Ey = 2 B0 cos sin sin(ωt − kz) Bx = − Ey (7.54)
kc a a b ω
Les conditions sur B sont alors vériées. La relation de dispersion s'écrit :
ω2
2
n2

2 m
−K =π2
+ 2 (7.55)
c2 a2 b
Chaque mode (n, m) est un ltre passe haut de pulsation de coupure :
r
m2 n2
ωcn,m = πc + (7.56)
a2 b2
Dans la bande passante on a :
c
vϕ = q 2 >c (7.57)
1 − ωcnm
ω
Condition monomode : pour a > 2b, le mode T E10 est le seul mode à se propager si ωc10 < ω < ωc20 , le
guide est monomode si a < λ < 2a. En ordre de grandeur, pour une bre standard on a a = 2.25b = 2.28
cm soit 2.28 cm< λ < 4.56 cm. Pour être sur d'être monomode on prend 2.42 cm< λ < 4.35 cm
Chapitre 8
Induction électromagnétique
L'induction est un exemple très important d'un phénomène électromagnétique. En eet, il s'agit d'un
exemple typique des eets du couplage entre E et B. De plus, c'est un peu "l'antichambre" de l'électroma-
gnétisme relativiste. Enn, ces phénomènes physiques ont une importance conceptuelle et pratique avec
une foule d'applications.
8.1 Énoncé historique de la loi de Faraday

On doit les lois de l'induction à Michel Faraday. La loi qui porte son nom date de 1831 :
Toute variation du ux du champ magnétique à travers un circuit mobile se traduit

par l'apparition d'un champ électrique à rotationnel non nul vériant :
˛
dφ
e= Ei .dl = −k
C dt
où :
• e est la force électromotrice
• Ei est le champ électrique vu dans le référentiel du circuit (pas forcément Galiléen)
• k est une constante qui dépend du système d'unités choisi.
• φ est le ux du champ magnétique vu dans le référentiel du laboratoire à travers le circuit.
Le signe "-" de la formule est la signature de la loi de Lenz qui stipule que le champ électromoteur créé par
la variation de ux crée un champ magnétique qui tend à s'opposer au champ qui lui a donné naissance.
On parle donc d'induction lorsqu'il y a variation du champ magnétique. Pour cela deux façons de faire :
• On peut faire bouger le circuit avec le champ magnétique constant : Induction de Lorentz
• Le circuit ne bouge pas et le magnétique varie dans le temps : induction de Neumann
Les deux façons de faire sont complètement équivalentes et peuvent être uniées grâce au principe de
relativité, qui souligne le fait que l'entité physique interessante n'est ni E ni B mais le couple (E,B) qui
constitue le champ électromagnétique.
Dans le cas de l'induction on se place toujours dans le cadre de l'ARQS magnétique et donc on néglige
le courant de déplacement dans l'équation de Maxwell Ampère, cependant l'équation de Maxwell-Faraday
reste inchangé.
89
90 CHAPITRE 8. INDUCTION ÉLECTROMAGNÉTIQUE
8.2 Cas de l'induction de Neumann

L'équation de Maxwell-Faraday nous permet d'écrire le champ électrique comme un potentiel :
∂A
E = −∇V − (8.1)
∂t
En calculant la circulation du champ électrique le long d'un circuit xe et fermé C :
˛ ˛ ˛
∂A
E.dl = − ∇V.dl − .dl (8.2)
C C C ∂t
Le premier terme est nul car la courbe C est fermé et comme C est xe on peut sortir la dérivée temporelle
de l'intégrale :
˛ ˛ ¨ ¨
d d d dφ(B)Σ
E.dl = − A.dl = − (rotA).dS = − B.dS = − (8.3)
C dt C dt Σ dt Σ dt
avec Σ une surface s'appuyant sur le contour C . La seconde égalité est obtenue grâce au théorème de
Stokes (cf chapitre 1). On désigne par φ(B)Σ le ux de B à travers la surface Σ. De plus, la circulation
de E s'identie à celle de la force électromotrice qui met les charges en mouvement (à la charge près).
On retrouve bien la loi de Faraday. Ici le champ E n'a pas de circulation nulle le long d'un courbe fermée
puisqu'il n'est plus le champ statique.
La séparation entre les 2 termes de E avec l'un à circulation nulle et l'autre pas est a priori arbitraire.
Cependant en se plaçant dans l'ARQS magnétique on la Jauge de Coulomb ∇.A = 0 et les équations
donnant V et A sont séparées et on peut identier le champ induit (celui à circulation non nul) avec ∂A∂t :
∂A
= Einduit : champ de Neumann (8.4)
∂t
8.3 Induction de Lorentz

On lit souvent que dans le cas où le circuit se déplace on a :
v = ve + vr (8.5)
avec ve la vitesse d'entraînement du circuit et vr la vitesse relative des porteurs par rapport au circuit.
Dans le circuit règne un champ magnétique qu'on suppose permanent donc on a E = −∇V . Le seul champ
de force appliqué est la force de Lorentz :
F = q(E + v ∧ B) (8.6)
Et donc la fem du circuit s'écrit comme le travail par unité de charge de la force de Lorentz :
˛ ˛
1
e= F.dl = (E + v ∧ B).dl (8.7)
q C C
En remplaçant : ˛ ˛ ˛
e=− ∇V.dl + (ve ∧ B).dl + (vr ∧ B).dl (8.8)
C C C
La première intégrale est nulle car la courbe C est fermé et la dernière également car vr est parallèle à dl.
Donc : ˛
e = (ve ∧ B).dl (8.9)
C
8.4. TRANSFORMATIONS DES CHAMPS E ET B SUIVANT LES RÉFÉRENTIELS 91
ainsi ve ∧ B se substitue à ∂A
∂t . On peut réécrire cette expression à l'aide d'un ux en identiant que ce
produit vectoriel correspond à un ux coupé par unité de temps et donc d'une dérivée temporelle par
rapport au temps du ux de B. On parle alors de la règle du ux coupé et on a :
dφ
e=− (8.10)
dt
avec φ le ux de B coupé par le circuit. Cette approche est celle faite en classe prépa et en licence qui
omet beaucoup de choses en adoptant ce point de vue. Essayons de clarier les choses.
8.4 Transformations des champs E et B suivant les référentiels

An de pouvoir préciser l'induction et plus particulièrement le cas de Lorentz, il faut s'intéresse à la
transformation du champ électromagnétique lors d'un changement de référentiel car dans ce cas le circuit
est en mouvement par rapport au référentiel du laboratoire où on eectue la mesure.
8.4.1 Cadre de l'étude

Nous nous placerons dans le cas non relativiste où la vitesse de déplacement du circuit v a une norme
bien inférieur à la vitesse de a lumière c. Nous supposons que les axiomes de la théorie de Newton sont
vraies et entre autre que la force est invariante par changement de référentiel. On va supposer que la
physique est Galiléenne.
Soit une charge q qui se déplace à une vitesse v par rapport à R0 et soit R1 en translation rectiligne
uniforme par rapport à R0 . Dans R0 on a des champs E0 et B0 qui agissent sur q :
F/q,R0 = q(E0 + v ∧ B0 ) (8.11)
Dans R1 on a un champ E1 et B1 et la force qui agit sur q est :
F/q,R1 = q(E1 + v1 ∧ B1 ) = q(E0 + v ∧ B0 ) (8.12)
où la deuxième égalité est donné par l'axiome de Newton et la vitesse v1 est simplement une composition
de vitesse :
d d
v1 = (O1 M) = (O1 O0 + O0 M) = −vR1 /R0 + vq/R0 (8.13)
dt dt
Donc :
q(−vR1 /R0 ∧ B1 + v ∧ B1 + E1 ) = q(v ∧ B0 + E0 ) (8.14)
Ce qui mène à (après avoir remarqué que vR1 /R0 = −vR0 /R1 ) :
E0 = E1 + vR0 /R1 ∧ B1
B0 = B1
Transformations T1
Le problème c'est que ces transformations sont FAUSSES ! En eet ce n'est pas compatible avec :
• Une particule avec une vitesse rectiligne uniforme. En eet, dans son référentiel propre, elle crée un
champ électrostatique mais pas de champ magnétique. Dans le référentiel du laboratoire, elle crée
pourtant un champ magnétique non nul ce qui est en contradiction avec la transformation T1 .
• Si on considère maintenant un conducteur inni chargé uniformément et qu'on le déplace à une

vitesse v1 . On a montré que le champ magnétique s'écrit en dehors du conducteur :
R2
B = µ0
jeθ (8.15)
2r
avec R le rayon du conducteur, r la distance à laquelle on se place. On a de plus j = vR1 /R0 ρ =
ρR 2
vR1 /R0 ρez avec ρ la densité de charges uniforme. De plus, on montre que E = 2r 0
er . Ainsi le champ
magnétique créé par le conducteur en mouvement s'écrit B = c2 vR1 /R0 ∧ E. Alors on en déduit que
1
la transformation s'écrit dans ce cas là :

E1 = E 0
B0 = B1 + 1
v ∧ E0
Transformations T2
c2 R1 /R0
8.4.2 Transformation des champs

Pour eectuer proprement la transformation des champs il faut se placer en relativité restreinte en posant
le quadrivecteur vC comme vecteur vitesse qui s'identiera à la vitesse de translation d'un référentiel par
rapport à l'autre vR1/R0 . Soit V = (a, b1 , b2 , b3 ) un quadrivecteur. On peut alors écrire :
V 0 = ΛV (8.16)
On a alors la transformation :
vC .b
a0 = γa − β (8.17)
c
vC (γ − 1)vC .(vC − b)
b0 = b−γ a+ (8.18)
c c2
avec :
1 vC
γ=p β= (8.19)
1 − β2 c
On peut prendre comme quadrivecteur : (cρ, j), (ct, r), (E, pc), (V, Ac) ou encore (E, cB) qu'on va expli-
citer.

0 vC ∧ cB vC (vC .E)
E = γ E+ + (γ − 1) 2
(8.20)
c {z c
| {z } | }
1 2
h vC i vC (vC .cB)
cB0 = γ cB − ∧ E + (γ − 1) 2
(8.21)
| {zc } | {z c }
3 4
Chirons les diérents termes tel que β = vC

c 1 (Galiliéen).
1 ' βcB
!
2 1
2 ' (γ − 1)β E = p − 1 β2E ' β4E
1−β 2
3 ' βE
4 ' β 4 cB
Si on prend brutalement c → ∞ alors :
β2

E0 = E + v C ∧ B + E + O(β 4 E, β 3 cB) (8.22)
2
2
E β
B0 = B − vC ∧ 2 + 4
Bc + O(β Bc, β E) 3
(8.23)
c 2
Si on ne prend pas le terme d'ordre 2 en β on a pas tout à fait du type T1 ou T2 . Cela pose également
deux problèmes :
• Pas stable car la transformation R0 → R2 est diérent de R0 → R1 → R2
• N'assure pas la covariance
En résumé, si on garde les termes en β dans le développement on suppose qu'on a pas de comportement
du type cB E ou cB E . Si on a aucun comportement pathologique :
E0 = E + vC ∧ B (8.24)
vC ∧ E
B0 = B − (8.25)
c2
Mais ce n'est pas stable et n'assure pas la covariance des équations de Maxwell. Peut-on tout de même
trouver une transformation de champ acceptable ?
Si on se place dans le cas E cB alors on a :
E0 = E (8.26)
vC ∧ E
B0 = B− (8.27)
c2
Et on retrouve la transformation de type T2 . Dans le cadre de l'ARQS électrique on retrouve une trans-
formation précédente.
Si à présent on se place dans l'ARQS magnétique et qu'on prend E cB alors les transformation
s'écrivent :
E0 = E + v C ∧ B (8.28)
0
B = B (8.29)
Et on retrouve la transformation T1 .
À présent plaçons dans l'ARQS et regardons comme se transforme le quadrivecteur (ρc, j) :

γvC .j
ρ0 c = γρc − (8.30)
c
(γ − 1)vC (vC .j) vC
j0 = γj − 2
− γρc (8.31)
c c
Dans le cas où ||vC || c, donc γ → 1 alors :
vC .j v 2
C
ρ0 c = ρc − +O (8.32)
c c
vC v 2
C
j0 = j − ρc +O (8.33)
c c
On a alors en ARQS magnétique que :
vC .j
ρ0 c = ρc − (8.34)
c
j0 = j (8.35)
Et en ARQS électrique que :
ρ0 = ρ (8.36)
0
j = j − ρvC (8.37)
On peut montrer maintenant que dans l'ARQS les équations de Maxwell sont invariantes par transforma-
tion de Galilée de la forme :
r 0 = r − vC t (8.38)
0
t =t (8.39)
Tout d'abord on a que :

∂ ∂x0 ∂ ∂t0 ∂ ∂
= 0
+ 0
= (8.40)
∂x ∂x ∂x ∂x ∂t ∂x0
Et donc ainsi on a ∇r = ∇r0 . Et pour le temps :
∂ ∂x0 ∂ ∂y 0 ∂ ∂z 0 ∂ ∂t0 ∂ ∂ ∂ ∂ ∂
= 0
+ 0
+ 0
+ 0
= −vCx 0 − vCy 0 − vCz 0 + 0 (8.41)
∂t ∂t ∂x ∂t ∂y ∂t ∂z ∂t ∂t ∂x ∂y ∂z ∂t
Soit,
∂ ∂
= 0 − (vC .∇r0 ) (8.42)
∂t ∂t
Cas ARQS électrique :

Supposons que dans le référentiel R, en se plaçant en ARQS électrique, on a les équations de Maxwell
suivantes :
ρ
∇.E =
0
∇.B = 0
∇∧E = 0
1 ∂E
∇ ∧ B = µ0 j +
c2 ∂t
Et les transformations des champs s'écrivent :
E0 = E
vC ∧ E
B0 = B −
c2
Si on transforme l'équation de Maxwell-Faraday, comme ∇r = ∇r0 et E0 = E, on a :
∇r ∧ E = ∇r0 ∧ E0 = 0 (8.43)
De même pour l'équation de Maxwell-Gauss comme on a ρ = ρ0 alors :
ρ ρ0
∇r .E = ⇒ ∇r0 .E0 = (8.44)
0 0
Regardons maintenant l'équation de Maxwell-Ampère :
v C ∧ E0

0 vC ∧ E 0 0 0
∇r ∧ B = ∇r0 ∧ B + = ∇r ∧ B + ∇r ∧
c2 c2
1 ∂E 1 ∂E 0 (vC .∇0r )E0
= µ0 j + 2 = µ0 (j0 + ρvC ) + 2 0 −
c ∂t c ∂t c2
On utilise la relation vectorielle :
∇ ∧ (a ∧ b) = a(∇.b) − ∇(a.b) (8.45)

vC 0 0 v
C

∇r ∧ B = ∇0r ∧ B0 + (∇ r .E ) − .∇ 0
r E0
c2 c2
vC ρ0 vC 0 0
= ∇0r ∧ B0 + 2 − .∇ E
c 0 c2 r
v
C
= ∇0r ∧ B0 + µ0 vC ρ0 − .∇ 0
E0
c2 r
Donc :
v
C
1 ∂E0 (vC .∇0r )E0
∇0r ∧ B0 + µ0 vC ρ0 − .∇0 E0 = µ0 (j0 + ρvC ) + − (8.46)
c2 r c2 ∂t0 c2
Donc :
1 ∂E0
∇0r ∧ B0 = µ0 j + (8.47)
c2 ∂t0
De même on trouve que :
∇0r .B0 = 0 (8.48)
Cas ARQS magnétique :

On peut eectuer la même chose pour l'ARQS magnétique. Supposons que dans le référentiel R, en se
plaçant en ARQS magnétique, on a les équations de Maxwell suivantes :
ρ
∇.E =
0
∇.B = 0
∂B
∇∧E = −
∂t
∇ ∧ B = µ0 j
Et les transformations des champs s'écrivent :
E0 = E + vC ∧ cB
B0 = B
De plus :
vC .j
ρ0 c = ρc − (8.49)
c
j0 = γj (8.50)
Prenons les équations dans l'ordre :
∇r0 .E0 = ∇r .E0

= ∇r .E + ∇r .(vC ∧ B)
ρ
= + (∇r ∧ vC ).B − vC .(∇r ∧ B)
0
ρ
= + µ0 vC .j
0
ρ vC .j
= +
0 0 c2
ρ0
=
0
L'équation de Maxwell-Gauss est invariante par transformation de Galilée.

∇r0 ∧ E0 = ∇r ∧ E0
= ∇r ∧ E + ∇r ∧ (vC ∧ B)
∂B
= − + (∇r .B)vC − (∇r .vC )B + (B.∇r )vC − (vC .∇r )B
∂t
∂B
= − − (vC .∇r )B
∂t
∂
= − + vC .∇r B
∂t
∂B
= − 0
∂t
∂B0
= − 0
∂t
L'équation de Maxwell-Faraday est invariante par transformation de Galilée.
On a immédiatement que ∇r0 .B0 = ∇r .B et de même ∇r0 ∧ B0 = µ0 j0 . Maxwell-Thomson et Maxwell-

Ampère sont invariantes par transformation.
Ainsi l'ARQS (dans le cas électrique ou magnétique) est la limite dans laquelle les transformations
galiléennes des champs font sens.
8.5 Cadre général de l'induction de Lorentz

Si on se place dans l'ARQS magnétique alors on a les transformations :
E1 = E0 + v1/0 ∧ B0 (8.51)
B1 = B0 (8.52)
Avec R0 et R1 deux référentiels galiléen et on prend un contour C(t) dans le référentiel R1 .
Expérimentalement on a : ˛
dφ0B
E0 .dl = −k (8.53)
C(t) dt
On sait que si v1/0 = 0 alors on a :
˛
dφB
E1 .dl = −k (8.54)
C(t) dt
On fait toujours l'hypothèse de la physique galiléenne où¸les forces sont identiques dans tous les référentiels.
Ainsi le travail de la force sur le contour C(t) qui vaut C(t) F.dl est indépendant du référentiel choisi.
Force électromotrice : ˛ ˛
1
e= E.dl = qE.dl (8.55)
C(t) C(t) q
est indépendant du référentiel choisi. Dans le référentiel R1 , si la vitesse de translation est constante :
˛ ˛ ˛ ˛
e= E1 .dl = [E0 + v1/0 ∧ B0 ].dl = E0 .dl + (v1/0 ∧ B0 ).dl (8.56)
C(t) C(t) C(t) C(t)
8.5. CADRE GÉNÉRAL DE L'INDUCTION DE LORENTZ 97
On a également : ˛ ¨
∂B0
E0 .dl = − .dΣ (8.57)
C(t) ∂t
Soit : "¨ ˛ #
∂B0
e(t) = − .dΣ − (v1/0 ∧ B0 ).dl (8.58)
∂t C(t)
On va relier ceci à la variation de ux de champ magnétique. Considérons la courbe C(t) à un instant t et
à un instant t+dt ainsi que deux points A et B sur cette courbe. On a :
¨ ¨
dφB
φB,t+dt − φB,t = dt φB,t+dt = B(r, t + dt).dΣ φB,t = B(r, t).dΣ (8.59)
dt S(t+dt) S(t)
Comme la divergence de B est nulle alors on a formellement que :

¨ ¨
− B(r, t).dΣ + B(r, t).dΣ + φlatéral = 0 (8.60)
S(t) S(t+dt)
Soit :
¨ ¨ ¨
∂B ∂B
B(r, t) + (r, t)dt .dΣ = B(r, t).dΣ + (r, t)dt.dΣ + O(dt2 ) (8.61)
S(t+dt) ∂t S(t+dt) S(t) ∂t
Donc :
¨ ¨ ¨
∂B
φ(t + dt) − φ(t) = B(r, t).dΣ + (r, t)dt.dΣ − B(r, t).dΣ + O(dt2 ) (8.62)
S(t+dt) S(t) ∂t S(t)
Soit : ¨
dφ ∂B φlatéral
= (r, t).dΣ − (8.63)
dt S(t) ∂t dt
Le problème est maintenant de savoir comment calculer le ux sur la surface latérale. Si le circuit est
indéformable il est facile d'exprimer ce terme. Prenons un élément de la courbe C(t) qui est en translation
à la vitesse v1/0 alors l'élément de surface latéral est perpendiculaire à dl et v1/0 et vaut dlv1/0 dt. Ainsi
la variation de ux latéral vaut :
ˆ
dφlatéral (t)
= B.(dl ∧ v1/0 ) (8.64)
dt C(t)
Ceci est vrai que si le circuit est indéformable et en translation. Soit :

¨
dφ ∂B φlatéral
− = − (r, t).dΣ +
dt S(t) ∂t dt
¨ ˆ
∂B
= − (r, t).dΣ + B.(dl ∧ vC )
S(t) ∂t C(t)
¨ ˆ
∂B
= − (r, t).dΣ + dl.(vC ∧ B)
S(t) ∂t C(t)
Et donc si vC est uniforme :

¨ ˆ
∂B dφ
e=− (r, t).dΣ + dl.(vC ∧ B) = − (8.65)
S(t) ∂t C(t) dt
Donc : ˛ ¨ ˆ
∂B dφ
e= E.dl = − (r, t).dΣ + dl.(v(t) ∧ B) = − (8.66)
S(t) ∂t C(t) dt
Et cette quantité est égale à − dφ

dt si le circuit est indéformable et en translation.
Attention : Pour appliquer e = − dφdt il faut que le circuit reste bien déni, c'est-à-dire que
la constitution du circuit ne doit pas changer. Par exemple si on prend un rhéostat manuel,
constitué d'un bobine avec un contact qui peut bouger le long des spires, si on place un champ
magnétique permanent dans la bobine et qu'on fait bouger le contacteur alors on ne voit pas de
courant apparaitre alors qu'il y a variation du ux du champ magnétique. En eet ici comme
B ne varie pas dans le temps on a pas d'induction de Neumann et les spires ne bougent pas
donc pas d'induction de Lorentz. Dans ce cas-là il vaut mieux revenir aux expressions locales a .
L'exemple classique de l'induction est l'expérience des rails de Laplace pour illustrer la formule
e = − dφ
dt qui n'a aucune raison de marcher car le circuit est déformable et pourtant c'est le cas.
Il s'agit simplement d'une coïncidence et surtout d'un système avec une géométrie très simple
qui se déforme tel que la loi de Faraday est quand même vériée.
a . Un autre exemple qu'on peut trouver est le disque de Faraday qu'on peut retrouver dans le BFR3.
Enn on peut réécrire le loi d'Ohm dans le cas où on a un circuit rigide qui a un mouvement à la vitesse
vC dans le référentiel du laboratoire. Dans le référentiel du circuit R0 les porteurs de charges ont une
vitesse d'ensemble hvR i et sont soumis à la force de Lorentz, la force électrique et éventuellement d'autre
actions et donc on peut écrire la loi d'Ohm dans le référentiel du circuit :

∂A
j = σ −∇V − + vC ∧ B + hvR i ∧ B (8.67)
∂t
Une démonstration plus complète de cette loi demande de faire un bilan des actions s'exerçant sur chaque
porteur de charge avant de faire une moyenne d'ensemble.
8.6 Courants de Foucault

On considère un conducteur de volume V dans un champ magnétique uniforme Ba , créé par des sources
extérieures. Ce champ magnétique induit un champ électrique E1 tel que :
∂Ba
∇ ∧ E1 + =0 (8.68)
∂t
Si le conducteur a une conductivité nie σ , alors des courants, appelés courants de Foucault,
vont se former tel que :
∂Ba
jf = σE1 ∇ ∧ E1 = − (8.69)
∂t
On peut analyser cela dans un géométrie simple : un conducteur cylindrique de rayon R et de hauteur
h. On applique un champ magnétique Ba (t) parallèle à l'axe de révolution. En supposant que le champ
magnétique s'identie au champ appliqué alors on a en un point M distant de r du centre :
1 dBa σr dBa
E1 = − r eθ jf = − eθ (8.70)
2 dt 2 dt
8.6. COURANTS DE FOUCAULT 99
Ainsi les courants de Foucault sont plus intenses en périférie du conducteur et leur direction obéit à la loi
de Lenz : le champ magnétique qu'il engendrent tend à s'opposer à Ba .
Pour le cas d'un champ sinusoïdal pour voir la dépendance en fréquence :
Ba = B0 cos(ωt) ez (8.71)
on trouve :
σω
jf = rB0 sin(ωt) eθ (8.72)
2
La puissance dissipée par eet Joule s'écrit :
ˆ ˆ
D E
j2f V 2 2 2
2
(8.73)

hP i = dτ jf .E1 = dτ = σω B0 R sin (ωt)
σ 8
Donc :
hP i σω 2 B02 R2
= (8.74)
V 16
La puissance dissipée par les courants de Foucault est proportionnelle à ω 2 et donc possède une forte
dépendance en la fréquence du champ magnétique appliqué. Ceci permet soit d'en tirer prot soit de
tenter de l'éviter en fonction des applications.
Il faut faire attention car ici nous nous sommes contenter que de l'ordre 1 des équations de Maxwell. En
eet nous nous sommes limité au cas où le champ magnétique se confondait avec le champ magnétique
appliqué, ce qui n'est pas strictement le cas puisque les courant induits vont créer un champ magnétique qui
va lui aussi induire un courant, qui vont eux même créer un nouveau champ magnétique, etc... Dans le cas
où B = Ba n'est plus vérié, typiquement pour les hautes fréquences il vaut mieux résoudre directement
les équations de Maxwell en prenant en compte la loi d'Ohm.
Exemple de la chute d'un aimant dans un tube conducteur 1 [Exercice]

Soit un tuyau cylindrique conducteur de rayon moyen a, d'épaisseur e et de conductivité σ . On a e a.
L'extrémité supérieure est en z = 0 et possède une longueur L susamment grande devant le régime
transitoire. L'axe z est dirigé vers le bas. On fait tomber un aimant A, supposé ponctuel, sans vitesse
initiale en z = 0. Son moment magnétique µ = µez est au cours du mouvement toujours vertical et dirigé
vers le bas. On note sa position za (t).
Le but de l'exercice est de trouver quelle est la force d'exercice le tuyau sur l'aimant par l'intermédiaire
des courants induits par la chute de l'aimant. On raisonnera sur un circuit C de cote z constitué d'un
tronçon de tuyau de hauteur dz . Un point M de cette boucle est repéré par ses coordonnées sphérique
d'origine l'aimant et l'orientation positive de la courbe est choisi selon +eϕ .
1. À partir du potentiel vecteur A créé par l'aimant, déterminer le champ électromoteur Em lié à sa
chute au niveau de C , puis le courant di induit dans C en notant v = dzdt la vitesse de chute de
a
l'aimant.
2. En déduire à l'aide d'une intégrale portant sur z , la projection Fz de la force exercée par le tuyau
sur l'aimant en évaluant son opposé Fz0 = Fz , c'est-à-dire la force exercée par l'aimant sur le tuyau.
3. Sachant qu'on donne la valeur de l'intégrale
ˆ +∞
x2 dx 5π
= (8.75)
−∞ (1 + x2 )5 128
calculer la force de freinage sur l'aimant et la mettre sous la forme F = −αv où α est indépendant
de za .
1. Se référer au Garing de magnétisme.
4. Intégrer l'équation diérentielle du mouvement et exprimer la vitesse v(t) et la cote za (t) de l'aimant
à l'aide d'un temps caractéristique τ du régime transitoire et de v` la vitesse limite atteinte.
5. Application numérique avec un tuyau de cuivre. a = 12 mm, e = 2 mm, L = 1 m, σ = 5.8 107 S/m,
m = 6 g et µ = 0.6 A.m2 .
1. Appelons θ l'angle sous lequel l'aimant voit le point M (angle des coordonnées sphériques). Le po-
tentiel vecteur A créé en un point M du circuit C par le dipôle magnétique est, avec sin θ = a/r avec
r = AM :
µ0 µ ∧ er µ0 µ sin θ µ0 µa
A(M ) =. = . eϕ = eϕ (8.76)
4π r2 4π r2 4π (a2 + (z − za )2 )3/2
De par la chute du petit aimant, za est fonction du temps et donc le potentiel vecteur aussi. Le champ
électromoteur du type Neummann au niveau de la courbe C s'écrit alors :
∂A(M ) ∂A(M )
Em (M ) = − =− veϕ (8.77)
∂t ∂za
avec v = dza
dt qui est toujours positif. Et donc :
3µ0 µa(z − za )v
Em (M ) = − . 2 eϕ (8.78)
4π (a ) + (z − za )2 )5/2
Le courant di induit par ce champ électromoteur dans le circuit C est :
di = j.dS = (σEm )(e dz eϕ ) = σEm e dz (8.79)
Et donc :
3µ0 σeµa(z − za )v
di = − dz (8.80)
4π (a2 ) + (z − za )2 )5/2
Tous ces calculs suppose que les champs sont uniformes sur l'épaisseur du tuyau et donc qu'on a bien
e a.
2. La force dF qu'exerce le circuit C sur l'aimant est l'opposé de la force dF0 qu'exerce l'aimant su le
circuit C . Ainsi, sur une longueur élémentaire dl du circuit au niveau duquel le champ B de l'aimant
est donné en coordonnées sphériques :
d2 F = didl ∧ B = di(dleϕ ) ∧ (Br er + Bθ eθ ) = didl(Br eθ − Bθ er ) (8.81)
Pour le circuit entier, la résultante des forces donne une composante suivant z par symétrie :
dFz0 = di2πa(−Br sin θ − Bθ cos θ) (8.82)
Or on rappelle les expressions du champ magnétique d'un dipôle :

µ0 2µ cos θ µ0 µ sin θ
Br = . , Bθ = . (8.83)
4π r3 4π r3
et donc :
3µ0 µdia cos θ sin θ
dFz0 = − (8.84)
2 r3
Or avec sin θ = a
r et cos θ = z−za
r et l'expression du courant di de la question précédente on trouve :
9µ20 σeµ2 a3 v(z − za )2 dz
dFz0 = (8.85)
8π(a2 + (z − za )2 )5
Et donc l'action du tuyau sur l'aimant est donné par l'intégration :

ˆ L
Fz = − dFz0 (8.86)
0
8.6. COURANTS DE FOUCAULT 101
3. An de calculer l'intégrale précédente il faut eectuer le changement de variable Z = z−z a

a
et on
trouve :
ˆ ˆ (L−za )/a
9µ20 σeµ2 a3 v L (z − za )2 9µ20 σeµ2 a3 v 1 Z2
Fz = − 2 2 5
dz = − dZ (8.87)
8π 0 (a + (z − za ) ) 8π a7 −za /a (1 + Z 2 )5
Or l'intégrale qui est donné dans l'énoncé possède des bornes innies. On ne peut remplacer −za /a
par −∞ dès que za > 0 (c'est-à-dire juste après le début de la chute) et (L − za )/a par +∞ tant que
za < L c'est-à-dire tant que l'aimant est dans le tube. Il n'est pas nécessaire que a za et a (L−za )
pour faire tendre les bornes vers l'inni mais il faut remarquer que la fonction est équivalente à 1/Z 8
pour des Z grands et tend donc très rapidement vers 0. Ainsi, étendre les bornes de l'intégrale ne
change pas la valeur de celle-ci car l'aire sous la courbe pour cette partie est négligeable. On trouve
alors une force qui vaut :
45µ20 σeµ2
Fz = − v = −αv (8.88)
1024a4
On trouve bien une force de type frottement uide avec un coecient α indépendant de za .
4. On applique le principe fondamentale de la dynamique à l'aimant :
dv
m = mg − αv (8.89)
dt
Ainsi, la vitesse limite vaut vl = mg
α et est atteinte au bout d'un temps de l'ordre de quelques τ = mα .
On a :
1024a4 mg
vl = (8.90)
45µ20 σeµ2
L'équation diérentielle projeté sur l'axe Oz donne :
dv v m
+ =g , τ= (8.91)
dt τ α
La solution, avec un vitesse initiale nulle est :
t

v(t) = vl 1 − e− τ (8.92)
et donc la position au cours du temps vaut (avec za (0) = 0) :

t
za (t) = vl t + vl τ e− τ − 1 (8.93)
5. On trouve une vitesse limite de vl = 0.42 m/s et un τ = 42 ms soit une position za (τ ) = 1.8 cm. Le
temps de chute tc (en considérant que l'exponentielle vaut 1) est L ' gτ tc et donc tc = 2.37 s. On
voit alors que la vitesse limite est très rapidement atteinte car τ tc et que le mouvement rectiligne
uniforme concerne la majeure partie de la chute car za (τ ) L.
Il est possible de réaliser l'expérience assez facilement mais il faut être conscient de plusieurs choses pour
faire du quantitatif :
• Connaître le moment dipolaire d'un aimant est compliqué.
• L'aimant ne tombe pas bien droit et son moment dipolaire fait un angle avec la verticale.
• La dépendance de σ avec la température n'est peut-être pas linéaire à basse température (azote
liquide). Il est par exemple possible de faire l'expérience avec du cuivre et de l'aluminium et d'avoir
un rapport des conductivité facilement. Pour mesurer le temps entre deux positions, il sut de placer
deux uxmètres à distance xe reliés à un oscilloscope. Le temps entre les deux pics et la distance
entre les uxmètres donne directement la vitesse limite si on ne se place pas en haut du tuyau.
Chapitre 9
Rayonnement dipolaire
Lorsqu'une charge est accélérée elle rayonne de l'énergie sous la forme d'un champ électromagnétique.
Ce rayonnement a mis en défaut le modèle planétaire de Bohr mais permet de créer des champs électro-
magnétique à partir de courant sinusoïdaux : ce sont des antennes. Le calcul du champ électromagnétique
est assez fastidieux et nécessite certaines hypothèses qu'il est important de préciser et développer.
9.1 Cadre de l'étude et hypothèses

Considérons le système suivant : un domaine D limité dans l'espace et en son sein des charges et des
courants. Quels sont les propriétés d'un tel système ? On notera les charges qi qui auront une vitesse vi .
Trois hypothèses sont nécessaires dans le cadre du rayonnement dipolaire :
• Hypothèse 1 : Les charges sont non relativistes et on peut alors écrire le potentiel vecteur pour des
charges continues ou en discrètes sous la forme :
ˆ

j r 0 , t − ||PM|||
µ0 c
A(M, t) = dτ (9.1)
4π P ∈D ||PM||

q v r 0 , t − ||Pi M|||
µ0 X i i c
A(M, t) = (9.2)
4π ||Pi M||
i∈D
De plus nous avons la Jauge de Lorenz qui est valable ∇.A + c12 ∂V
∂t = 0 ainsi si A est connu alors V
l'est également avec : ˆ t
V (M, t) = V (M, 0) − c2
∇.A(M, t0 )dt0 (9.3)
0
Le problème est que dans le cas général, même dans une géométrie simple on ne peut pas calculer
A. Une simplication est alors nécessaire.
• Hypothèse 2 : Si a est la taille typique de D et O le barycentre des charges de la distribution alors
on se place en un point M de l'espace susamment loin, c'est-à-dire tel que :
||OM|| a (9.4)
l'hypothèse dipolaire.
c'est
• Hypothèse 3 : On suppose que les retards entre les diérents points des sources sont identiques.
Prenons deux points P1 et P2 de la distribution distant de a. La diérence de retard entre les deux
points est a/c. Si les charges et les courants varient sur un temps typique τ alors on doit avoir :
a
τ ⇒ a τc (9.5)
c
c'est l'ARQS des sources.
103
104 CHAPITRE 9. RAYONNEMENT DIPOLAIRE
Si on a les trois hypothèses vériées alors il est possible de faire quelques calculs. Dans l'expression de A
l'hypothèse 2 nous permet de remplacer le ||PM|| 1
par ||OM||
1
et l'hypothèse 3 nous permet de remplacer
dans la dépendance temporelle de jvi de ||PM|| par ||OM||. Et donc :
ˆ

j r 0 , t − ||OM|||
µ0 c
A(M, t) = dτ (9.6)
4π P ∈D ||OM||

q v r 0 , t − ||OM|||
µ0 X i i c µ0 1 X
0 ||OM|||
A(M, t) = = qi v i r , t − (9.7)
4π ||OM|| 4π ||OM|| c
i∈D i∈D
||OM ||
Dénissons t0 =t− c alors on a :

||OM||| d ||OM|||
vi t − = 0 ri t − (9.8)
c dt c
Soit : " #
µ0 1 d X ||OM|||
A(M, t) = qi ri r0 , t − (9.9)
4π ||OM|| dt0 c
i∈D
Or on peut dénir le moment dipolaire de la distribution par :
X
P(t0 ) = qi ri (t0 ) (9.10)
i
Et la limite continue donne : ˆ

0
P(t ) = dτ ρ(P, t0 )OP (9.11)
Et donc on a pour la description discrète ou continue :
µ0 Ṗ(t0 ) dP(t0 )
A(M, t) = avec Ṗ(t0 ) = (9.12)
4π ||OM|| dt0
9.2 Expression du champ électromagnétique

Si le point d'observation M est situé à une grande distance de la distribution alors la connaissance du
potentiel vecteur est susante pour déterminer le champ électromagnétique. En eet comme B = ∇ ∧ A
on peut calculer B et l'intégration de l'équation de Maxwell-Ampère permet de remonter à E. Pour le
calcul on xe l'axe z selon le moment dipolaire de la distribution :
P(t0 ) = P (t0 )ez (9.13)
Ceci n'est pas une perte de généralité car si on a un mouvement quelconque alors on pourra toujours par
superposition linéaire exprimer le moment dipolaire de la distribution comme la somme de trois moments
dipolaires de directions xées. On a :
µ0 Ṗ (t0 )
A(M, t) = ez (9.14)
4π ||OM||
On se place en coordonnées cylindrique (ρ, ϕ,z ) 1 . On note r = OM et θ l'angle que fait le vecteur r avec
l'axe z . Alors en utilisant l'expression du rotationnel en cylindrique :
∂Az
B(M, t) = ∇ ∧ A(M, t) = − eϕ (9.15)
∂ρ
1. On ne pourra pas confondre avec la densité de charge car elle n'apparaitra pas.
9.2. EXPRESSION DU CHAMP ÉLECTROMAGNÉTIQUE 105
et comme on a r = z 2 + ρ2 , t0 = t − r/c et sin θ = ρ/r :

p
∂r ρ ∂ ṗ ∂ ṗ ∂r ∂(t − r/c) ∂r p̈ sin θ

=p = sin θ ⇒ = = p̈(t0 ) =− (9.16)
∂ρ z 2 + ρ2 ∂ρ ∂r ∂ρ ∂ρ ∂ρ c
Soit :
µ0 ∂ ṗ µ0 ṗ ∂(1/r) µ0 p̈ ṗ
B(M, t) = − eϕ − eϕ = + 2 sin θeϕ (9.17)
4πr ∂ρ 4π ∂ρ 4π cr r
On peut le réécrire comme :

µ0 p̈ ṗ
B(M, t) = + ∧ er (9.18)
4π cr r2
On voit que la contribution majoritaire du champ magnétique provient majoritairement, à grande distance,
aux modications temporelles du dipôle électrique de la distribution.
Dans la zone où est calculé le champ on a par hypothèse que le courant est nul. On a donc par Maxwell-
Ampère :
ˆ t
E(M, t) = E(M, 0) + c2
dt00 ∇ ∧ B(M, t00 ) (9.19)
0
On peut eectuer le calcul et on trouve 2 :

1 p ṗ Qtot sin θ p ṗ p̈
E(M, t) = 2 cos θ + + 2 er + + + eθ (9.20)
4π0 r3 r2 c r 4π0 r3 r2 c rc2
où Qtot est la charge totale de la distribution qui est nécessairement invariante au cours du temps pour
assurer la conservation globale de la charge électrique.
Si à présent on revient sur les hypothèses qu'on a faites, à savoir l'observation à grande distance et négliger
le retard entre deux points de la sources, on peut regarder en ordre de grandeur :
ṗ(t0 ) p p̈(t0 ) p
2
' 2 2
' 2 2 (9.21)
r c r cτ rc rc τ
On peut alors envisager deux limites :
• Zone de champ proche (a r cτ ) : même si la distance est grande devant la taille de la
distribution elle reste petite devant la distance parcourues par l'information sur une durée typique
de variation. Dans ce cas on a :
µ0 ṗ(t) ∧ r 3(p(t).er )er − p(t) Qtot

B(M, t) ' , E(M, t) ' + er (9.22)
4π r3 4π0 r3 4π0 r2
et les champs perçus sont alors conformes à ceux calculés dans une limite quasi-stationnaire. Ils
décroissent au moins proportionnellement à l'inverse du carré de la distance aux sources. On note de
plus que l'état électrique du système de sources est perçu sans retard (p(t) et non plus p(t0 )).
• Zone de rayonnement (a cτ r) : la distance d'observation est cette fois-ci très grande et le
retard induit par la vitesse nie de propagation de l'information est primordial. On a alors :
µ0 p̈(t0 ) sin θ p̈(t0 ) sin θ

B(M, t) ' eϕ , E(M, t) ' eθ (9.23)
4πcr 4π0 c2 r
2. On pourra trouver certains détails de calculs dans le Pérez d'électromagnétisme.
et les champs perçus décroissent proportionnellement à l'inverse de la distance aux sources. Dans la zone
de rayonnement on a donc des champs qui décroissent lentement.
Pour la majorité des applications, les situations physiques rencontrées sont celles de la zone de rayonnement
et il convient de détailler ce cas. Les champs perçus à grande distance possèdent plusieurs caractéristiques :
1. Ils décroissent lentement (en 1/r)
2. Ils sont anisotropes et sont nuls en tout point situé suivant l'axe du moment dipolaire
3. La dépendance temporelle en t0 = t − r/c indique que le champ créé à grande distance se propage
suivant er à une vitesse c
4. Localement, les deux champs sont perpendiculaires entre eux et le champ électromagnétique a une
structure d'onde plane.
5. Le champ électromagnétique à grande distance est généré par les dérivées temporelles secondes du
moment dipolaire, c'est-à-dire qu'il n'existe que si les charges en mouvement sont accélérées.
9.3 Comportement énergétique du rayonnement

9.3.1 Le vecteur de Poynting
2
En statique on veut en général interpréter 20 E 2 et 2µ
B
0
comme des termes d'énergie volumique. Soit V un
volume de l'espace et S la surface fermée qui le délimite. Notons U (t) l'énergie électromagnétique stockée
au temps t dans le volume. On a alors :
U (t + dt) = U (t) + δU (t) (9.24)
où on a :
δU (t) = δUsurface (t) + δUvolume (t) (9.25)
Soit Π(r, t) le vecteur dont le ux à travers dS indique la puissance transmise à V à travers cette surface :
δe = −Π(r, t).dS (9.26)
Soit σ(r, t) la source volumique d'énergie dans le volume dτ et pendant le temps dt :
δe = σ(r, t)dτ dt (9.27)
Le bilan s'écrit alors :

˚ ‹
U (t + dt) = U (t) + d3 r0 σ(r0 , t)dt − Π(r0 , t)dSdt (9.28)
V S
Donc : ˚ ˚
dU 3 0 0
= d r σ(r , t) − ∇r0 .Π(r0 , t)d3 r0 (9.29)
dt V V
Ceci est vrai pour tout volume V et donc en notant u la densité d'énergie volumique alors :
∂u(r, t)
= σ(r, t) − ∇r .Π(r, t) (9.30)
∂t
Il s'agit maintenant de trouver les valeurs des grandeurs que nous avons introduites. Pour cela il faut
utiliser les équations de Maxwell :
1 1 ∂E 2

∂ 1
(∇ ∧ B).E = µ0 j.E + 2 = µ0 j.E + µ0 0 E 2
(9.31)
c 2 ∂t ∂t 2
9.3. COMPORTEMENT ÉNERGÉTIQUE DU RAYONNEMENT 107
Or on a ∇.(A ∧ B) = A.(∇ ∧ B) + B.(∇ ∧ A) donc :

∂B
E.(∇ ∧ B) = B.(∇ ∧ E) − ∇.(E ∧ B) = −∇.(E ∧ B) + B. − (9.32)
∂t
Au nal on trouve :
0 E 2

E∧B ∂ 1 2
−∇. = j.E + + B (9.33)
µ0 ∂t 2 2µ0
Il faut comparer ça à :
∂u(r, t)
−∇r Π(r, t) = −σ(r, t) + (9.34)
∂t
Le choix naturel pour est le choix de Poynting :
E∧B 0 E 2 1 2
Π= σ = −j.E u= + B (9.35)
µ0 2 2µ0
Mais on aurait pu prendre :
E∧B ∂Y 0 E 2 1 2
Π= +∇∧X+ σ = −j.E u= + B − ∇.Y (9.36)
µ0 ∂t 2 2µ0
Qu'on appelle jeu de Mac Arthur.
9.3.2 Cas du rayonnement dipolaire

Dans le cas du rayonnement dipolaire à grande distance (c'est-à-dire dans la zone de rayonnement) on
trouve directement que :
p̈2 (t0 ) sin2 θ
Π(r, t) = er (9.37)
16π 2 0 c3 r2
On voit que dans cette expression l'anisotropie car il n'y a pas d'énergie rayonné le long de l'axe du moment
dipolaire. La dépendance en 1/r2 montre la conservation de l'énergie.
Le diagramme de rayonnement consiste à tracer la fonction :
||Π(r, θ, ϕ, t0 )||
F (θ, ϕ) = (9.38)
||Πmax (r, θmax , ϕmax , t0 )||
où Πmax est la valeur maximale à r et t0 xé. Dans le cas présent on a :
F (θ, ϕ) = sin2 θ (9.39)
Ce qui correspond à deux lobes symétriques.
On peut regarder la puissance rayonnée par une sphère :

‹ ˆ π ˆ 2π ˆ π
0 0 p̈2 (t) sin2 θ p̈2 (t0 )
P(t) = Π(r , t ).dS = dθ 2
dϕr sin θ = dθ sin3 θ (9.40)
sphère θ=0 ϕ=0 16π 2 0 c3 r2 8π0 c3 θ=0
or l'intégrale vaut 4/3 donc :

p̈2 (t0 )
P(t0 ) = (9.41)
6π0 c3
On voit alors que la puissance varie comme le carré de p̈ et également qu'elle ne dépend pas de r.
On peut alors traiter deux cas :

1. Cas d'un charge isolée : Dans ce cas là on a P(t) = OP(t)q et donc en prenant l'accélération selon
une direction on a : P̈(t) = qz̈(t)ez = qa(r, t0 )ez . On a donc :
a2 (r, t0 )q 2
P(t) = (9.42)
6π0 c3
Ainsi une charge accélérée rayonne de la puissance et c'est pourquoi le modèle planétaire de Bohr
ne peut pas marcher dans ce cas là.
2. Dipole de Hertz qui vérie :
P(t0 ) = qa cos ωt0 (9.43)
On se place alors dans le cas non relativiste, dipolaire ainsi que dans l'ARQS des sources c'est-à-dire
que aω c. On peut alors calculer la puissance rayonné sur une période :
ˆ T
1 a2 ω 4 q 2
< P >T = P (t0 )dt0 = (9.44)
T 0 12π0 c3
9.4 Rayonnement d'antenne

9.4.1 Antenne électrique
Une antenne est un morceau de conducteur avec des charges en mouvement. On notera l la taille du
conducteur et S la section du conducteur. On a des charges qi possédant chacune une vitesse vi alors on
a:
X X qi vi
P= qi r i j= (9.45)
lS
i i
Et donc : P
i qi v i
I = jS = (9.46)
l
Soit : !
dP(t0 ) d X X
0
= 0 qi ri (t0 )) = qi vi (t0 )) = I(t0 )l (9.47)
dt dt
i i
Soit :
P̈ (t0 ) = I(t
˙ 0 )l (9.48)
La puissance rayonnée est donc donnée par :
I˙2 (t0 )l2

P= (9.49)
6π0 c3
Ainsi pour gagner de la puissance à intensité xé il faut augmenter la longueur du conducteur l mais si on
augmente l on va perdre la condition d'ARQS l cτ . De plus si l'ARQS des sources est vériée on a :

dI 1
E(r, t) = l sin θ eθ (9.50)
dt 4π0 c2
9.5. RAYONNEMENT DE LA MATIÈRE 109
9.4.2 Antenne réelle

Si on prend une antenne rectiligne réelle, le problème est que celle-ci n'est pas dans l'ARQS. En revanche,
on peut la décomposer en morceaux élémentaires qui le sont. On est toujours dans l'approximation dipolaire
et donc pour une portion d'antenne dz située au point P :
dI ˆ
dz sin θ dt 0 (z) 0
dEP (M ) = eθ E(M, t ) = dEP (M ) (9.51)
4π0 c2 r P
On fait alors une TF temporelle du signal I(z, t) :

ˆ
1
I(z, t) = dωejωt I(z, ω) (9.52)
(2π)1/2
Pour une pulsation ω xé ou si on considère que I est monochromatique alors on aurait I(z, ω) = f (z).
On a alors :
I(t, z) = f (z)ejωt I(t0 , z) = f (z)ejω(t−r/c) (9.53)
On a alors : ˆ
dz sin θf (z)jωejω(t−P M/c)
E(M, t) = eθ (9.54)
antenne 4π0 c2 r
On est dans la zone de rayonnement et on a :
OM.OP 1/2

||PM|| = ||PO + OM|| = ||OM|| 1 − 2 ' r − z cos θ (9.55)
||OM||2
ˆ
jω r z
E(M, t) = dz sin θf (z)ejω(t− c + c cos θ) eθ (9.56)
4π0 c2 r antenne
ˆ ˆ
jω jωt0 jω zc cos θ) jω jωt0 2πz
E(M, t) = 2
e dz sin θf (z)e eθ = 2
e dz sin θf (z)ej λ cos θ) eθ
4π0 c r antenne 4π0 c r
| antenne {z }
F (θ)
Ainsi le champ électrique est proportionnel à la transformée de Fourier de la répartition de courant f (z).
Par exemple pour une antenne rectiligne on a :

2π l
f (z) = I0 sin |z| − (9.57)
λ 2
où λ est la longueur d'onde associé à la pulsation ω dans le vide.
cos πl πl

λ − cos λ cos θ
F(θ) = (9.58)
sin2 θ
Ainsi pour avoir la résonance d'antenne il faut avoir une longueur d'antenne l = 2 .
nλ
Le lecteur pourra se demander comme faire pour avec 4 antennes on peut émettre dans une direction
donnée.
9.5 Rayonnement de la matière

9.5.1 Expérience introductive
On place du lait (entier de préférence) dans une cuve à eau. On fait traverser un rayon, on constate :
• La lumière au bout est plus rouge
• On constate un eet de polarisation
• Une diusion bleuté
9.5.2 Un modèle simple

On considère un atome soumis à un champ E. On suppose que :
• Le noyau est très lourd par rapport aux électrons : pas de mouvement du noyau
• On suppose que les électrons répartis de façon sphérique et on note Z le numéro atomique de l'atome
donc le nombre d'électron ou de protons.
Au repos les barycentres des charges positives et négatives sont confondus. Avec un champ électrique
extérieur alors les barycentres ne sont plus confondus et on note r = PN avec P le barycentre des charges
positives et N le barycentre des charges négatives. Appliquons le PFD pour la sphère électronique :
Zme r̈ = FP/N + (−Ze)Eext (r, t) + Fmagnétique + fdissip (9.59)
Regardons les termes un à un :

• Le second terme correspond à la force électrique du champ sur les électrons. Si la longueur d'onde
λ du rayonnement est grande devant la taille du nuage électronique a alors on peut supposer que le
champ est uniforme sur la sphère électronique et donc Eext (r, t) = Eext (t) dans l'équation
• Le troisième terme est la partie magnétique de la force de Lorentz qui est négligeable si ||qv ∧ B||
||qE||. Or on peut écrire que ||E|| ' vϕ ||B|| avec vϕ la vitesse de phase de l'onde. On suppose alors
que la vitesse des particules est très petite devant vϕ ce qu'il faudra vérier par la suite.
• Le quatrième terme est un terme de dissipation dont on ne connait pas a priori la forme. On prend
alors une forme linéaire et simple fdissip = − mτe v ce qui donne une puissance dissipée en ω 2
• Pour exprimer le premier terme qui est la force électrostatique du noyau positif sur la sphère négative
on va commencer par utiliser le principe de réciprocité et dire que FP/N = −FN/P . Or comme on
considère le noyau comme ponctuel on a :
FN/P = Esphère (P ).Ze (9.60)
Or on montre facilement que :
ρr −Zer −ZeNP
Esphère = er = 3
er = (9.61)
30 4πa 0 4πa3 0
Et donc :
(Ze)2 NP −(Ze)2 PN
FP/N = −FN/P = = (9.62)
4πa3 0 4πa3 0
Donc le PFD s'écrit alors :
Z 2 e2 Zme
Zme r̈ = − 3
r − ZeEext − ṙ (9.63)
4πa 0 τ
Si on introduit alors le moment dipolaire P = −Zer = −ZeNP alors :
Ze2 Ze2 Ṗ
P̈ = − 3
P − Eext (t) + (9.64)
4π0 a me me τ
On peut alors distinguer deux cas :

• Si on prend un champ électrique constant E = E0 et qu'on prend P = 0 α0 E0 , où α0 ' 4πa3 qui est
à peut près le volume de la sphère. Si à t = 0 on coupe le champ alors l'équation de la forme :
d2 P 1 dP Ze2
+ + ω02 P = 0 ω02 = (9.65)
dt2 τ dt 4π0 a3 me
Et on a l'équation diérentielle de relaxation de la polarisation et on voit apparaitre τ comme le

temps de relaxation.
• Si maintenant on se place en notation complexe pour voir la réponse en fréquence du système :

iω Ze2
−Pω ω 2 + Pω + ω02 Pω = Eω (9.66)
τ me
Et donc :
Ze2
me
Pω = Eω (9.67)
ω0 − ω 2 + iω
2
τ
Puis pour remonter à P(t) on peut eectuer la TF inverse. On peut alors poser :
α0 Ze2
α(ω) = 2 α0 = (9.68)
me 0 ω02

ω
1− ω0 + i τ ωω2
0
Et on obtient la polarisabilité du matériau en fonction de la fréquence et on peut tracer son module

sur la Figure 9.1 qui a un comportement compliqué avec une résonance :
1 ω
α0 Q2 ω0
|α(ω)| = " #1/2 tan Φ = 2 (9.69)
2 2
2 ω
− 1
1 − ωω0 + 1
Q2
ω
ω0
ω0
0
102
-1
α/α0
100
Φ
-2
-3
-2
10
-4
10-2 10-1 100 101 102 10-1 100 101
ω/ω0 ω/ω0
Figure 9.1 Variation du module et de la phase de la polarisabilité normalisée αα0 en fonction de la
pulsation normalisée ωω0 . La courbe rouge correspond à un Q = 10 et la courbe noire à Q = 106 .
On peut alors distinguer 3 régimes particuliers :

• Régime de Rayleigh (ω ω0 ) : Dans ce cas là : on trouve :
2 !
ω 1
|α(ω)| ' α0 1 + 1− Φ=0 (9.70)
ω0 2Q2
Lorsque les fréquences d'excitation sont faibles devant la fréquence propre de vibration de l'atome, le
dipôle généré voit son amplitude contrôlée majoritairement par la susceptibilité statique et celui-ci
oscille en phase avec le champ inducteur.
• Régime de résonance (ω ' ω0 ) : La résonance a lieu lorsque |α(ω)| est maximum c'est à dire
pour :
1
ω = ωr = ω0 1 + (9.71)
2Q2
Or en pratique le facteur de qualité de l'oscillateur est très grand (108 pour le césium par exemple) et
la fréquence de résonance s'identie à la fréquence propre de vibration du nuage électronique. Dans
le régime de résonance on a alors :
α0 Q π
|α(ω)| ' 2 1/2 Φ=− (9.72)
2
ω−ω0
1 + 4Q2 ω0
On remarque alors que pour |ω − ω0 | ω0 , la polarisabilité de l'atome évolue de façon lorentzienne.

• Régime de Thomson (ω ω0 )
ω 2 ω 2 1

0 0
|α(ω)| ' α0 1+ 1− Φ = −π (9.73)
ω ω 2Q2
Lorsque les fréquences d'excitation sont beaucoup plus grandes que la fréquence propre de vibration
de l'atome alors l'amplitude du dipôle induit décroît proportionnellement au carré de la pulsation.
Le dipôle induit est alors en opposition de phase avec le champ inducteur.
On peut maintenant calculer la puissance totale rayonnée par un dipôle :

ω 4 |Eω |2 |α(ω)|2
Pray = (9.74)
12π0 c3
Dans l'expérience abordée en début de chapitre, le champ électrique employé appartient au domaine visible
ω ' 1014 rad.s−1 et les particules qui vont se polariser sous l'action du champ électrique sont des molécules
lipidiques du lait. Si la quantité de lait introduite dans l'eau est faible alors ces molécules ne se regroupent
pas sous forme de micelles géantes et leurs tailles restent petites devant les longueurs du spectre visible.
On a alors ω0 ' 1016 rad.s−1 et nous sommes dans le cas ω ω0 . La puissance rayonnée devient alors :
ω 4 |Eω |2 α02
PrayRayleigh = (9.75)
12πc3 0
La puissance rayonnée par les dipôles, proportionnelle à ω 4 , est donc d'autant plus importante que la
fréquence est élevée. Le rayonnement bleuté observé vient précisément de cela : les dipôles absorbent
ne partie du rayonnement incident et le réémettent sous forme de rayonnement dipolaire et dans tout
l'espace. C'est pour les pulsations les plus importantes (correspondant au bleu) que ce phénomène est
le plus important. Le faisceau de lumière blanche initial va donc au cours de son parcours dans la cuve,
s'appauvrir petit à petit en rayonnement haute fréquence. Tout ceci est bien entendu à relier avec la couleur
bleue du ciel ainsi que la couleur rouge du coucher de soleil. On peut également remarquer que la lumière
diusée est polarisée et ici encore les caractéristiques du rayonnement dipolaire permettent d'expliquer
cette propriété : la lumière émise initialement est non polarisée et le champ électrique oscillateur est
transverse. Si la direction de propagation est selon l'axe z alors des moments dipolaires sont excités de
façon équiprobable selon la direction x ou y . Si la direction d'observation se fait selon l'axe x alors aucune
puissance des dipôles ne nous parvient des dipôles excités suivant cette direction et donc le champ électrique
reçu par l'observateur ne provient que des dipôles polarisés selon l'axe y qui émettent un champ électrique
polarisé suivant cette direction : on a donc bien polarisation de la lumière.
9.5.3 Auto-consistance du modèle

On a pris comme force de dissipation une force proportionnelle à la vitesse des électrons. Est-ce que ceci
est bien consistant avec la puissance rayonnée ?
A priori non car pour une excitation sinusoïdale on a la puissance mécanique dissipée en ω 2 alors que la
puissance rayonnée par un dipôle est en ω 4 . On peut alors considérer deux situations physiques simples :
1. Le cas d'une charge uniformément accéléré. On peut alors comparer l'énergie mécanique de la par-
ticule à l'énergie rayonnée pendant un temps T :
1 1
Eciné = mv 2 = m(aT )2 (9.76)
2 2
q 2 a2 T τc q2
Eray = 3
= Eciné τc = ' 10−14 s (9.77)
6π0 c T 6π0 c3 m
Si l'accélération est constante sur un temps T τc alors l'énergie rayonnée est une faible quantité
d'énergie cinétique.
2. Sur un mouvement périodique de la forme :
d = d0 cos ω0 t (9.78)
On a alors comme T = 2π/ω0 :
q 2 d20 ω04 2π τ
Eméca = mω02 d2 Eray = 3
' Eméca (9.79)
12π0 c ω0 T
Ainsi pour des systèmes usuels, on a que le rayonnement est une faible perturbation.
On peut prouver que à faible perturbation la forme qu'on a prise est la bonne au premier ordre.

Electromagnetisme Dans Le Vide

Transféré par

Droits d'auteur :

Formats disponibles

Electromagnetisme Dans Le Vide

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Electromagnetisme Dans Le Vide

Transféré par

Droits d'auteur :

Formats disponibles

Électromagnétisme dans le vide

Préparation à l'agrégation de l'ENS de Lyon 2020-2021

Bibliographie pour aller plus loin

Bibliographie pour méthodes mathématiques

Épreuves d'agrégation lié à l'électromagnétisme

2.7.4 Développement multipolaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

5 Approximation des régimes quasi-stationnaires 55

6.6 Magnétorésistance et eet Hall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75

9 Rayonnement dipolaire 103

1.1 Analyses vectorielles

On parle de champ uniforme lorsqu'il ne dépend pas de l'espace et permanent ou stationnaire

On appelle ligne de champ, à un instant donné t, la courbe tangente en chaque point M de

Notons deux choses :

Le sens du vecteur de surface élémentaire dS se déduit de l'orientation de la courbe par la règle de la

Figure 1.1  Calcul du ux sortant d'un cube innitésimal.

De même pour la face 2 :  

où Vy = V.n est la composante du champ selon ey .

où ∇.V est appelé divergence du champ vectoriel V.

Figure 1.2  Calcul de la circulation en découpant la courbe en deux.

Figure 1.3  Calcul de la circulation sur un carré.

où n est le vecteur unitaire normal en tout point à la surface S .

1.1.5 Laplacien scalaire

∂2U ∂2U ∂2U

1.1.6 Laplacien vectoriel

1.2 Relations vectorielles

1.2.1 Théorème de Stokes

où dS est le vecteur normal à la surface élémentaire et proportionnel à la dite surface.

1.2.2 Théorème de Green-Ostrogradsky

1.2.4 Égalités vectorielles

1.3 Diérents systèmes de coordonnées

Figure 1.4  Système de coordonnées cartésien.

Figure 1.5  Système de coordonnées cylindrique.

Figure 1.6  Système de coordonnées sphérique.

1.4 Fonction de Green

2.1 Introduction et historique

La contribution de Maxwell consiste essentiellement à l'ajout du terme c12 ∂E(r,t)

2.2 Cadre de l'électrostatique

2.3 Approche intégrale et loi de Coulomb

Figure 2.1  Deux charges ponctuelles en interaction.

La force électrostatique (loi de Coulomb) s'écrit alors :

2.3.2 Champ de Coulomb

2.3.3 Répartition continue de charges

2.4 Liens avec les équations de Maxwell

Or comme ∇r ∧ (∇r f ) = 0 pour tout f cela signie que :

Et donc ED (r) vérie bien la seconde équation de l'électrostatique.

2.4.2 Équivalence des deux approches

∇.V = a(r) (2.28)

ik.V = A(k) (2.30)

2.5 Théorème de Gauss

Ainsi le théorème de Gauss s'écrit :

2.6 Structure du champ électrostatique

2.6.2 Exemples d'invariants

2.6.3 Propriétés de symétrie des champs électriques

De la même manière on peut montrer que si on a un plan d'antisymétrie Π− de la distribution de charge

2.6.4 Calculs de champ électrostatique

La surface de Gauss Σ naturelle pour ce problème est un cylindre de hauteur h et de rayon r. Le

∀z1 , z2 > 0 E(z1 ) = E(z2 ) (2.54)

Il est possible de montrer la même chose pour les z < 0.

Figure 2.2  Disque uniformément chargé.

2.7 Potentiel électrostatique

6.6 Magnétorésistance et eet Hall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75

Figure 1.1 Calcul du ux sortant d'un cube innitésimal.

De même pour la face 2 :

Figure 1.2 Calcul de la circulation en découpant la courbe en deux.

Figure 1.3 Calcul de la circulation sur un carré.

1.3 Diérents systèmes de coordonnées

Figure 1.4 Système de coordonnées cartésien.

Figure 1.5 Système de coordonnées cylindrique.

Figure 1.6 Système de coordonnées sphérique.

Figure 2.1 Deux charges ponctuelles en interaction.

Or comme ∇r ∧ (∇r f ) = 0 pour tout f cela signie que :

Et donc ED (r) vérie bien la seconde équation de l'électrostatique.

Figure 2.2 Disque uniformément chargé.

2.7.3 Signication physique

Figure 2.3 Dipôle électrostatique.

On peut alors vérier la valeur de ∇ ∧ B :

Fil rectiligne inni par Biot et Savart

Figure 3.1 Champ magnétique créé par un l inni.

On dénit alors proprement la densité surfacique σ comme :

En dénissant la composant selon x du courant surfacique comme :

avec la même dénition de n12 .