Cours d’électromagnétisme classique

Valentin Clarisse

Dans ce cours, on adopte une approche historique de la théorie de l’électromagnétisme.

Dans un premier temps, on considère la loi de Coulomb (1785), donnant le champ

électrique créé par une charge immobile, comme un postulat, et on en déduit toute la
théorie de l’électrostatique. On terminera cette section par le dipôle électrostatique.

Dans un deuxième temps, on partira de la loi de Biot et Savart (1820), pour

construire la théorie de la magnétostatique.

Dans un troisième temps, on verra comment concilier les deux théories précédentes
en une théorie de l’électromagnétisme en régime variable, où l’équation des ondes
joue un rôle central.

Nous terminerons en étudiant quelques milieux : le vide, les conducteurs, et les

plasmas, avec une ouverture sur l’optique ondulatoire.

Note temporaire : ce cours n’est qu’une ébauche en construction. Certaines parties

et points de vue développés ici ne sont pas au programme de CPGE.

1
Table des matières
1 Électrostatique 5
1.1 Notion de charge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 De la loi de Coulomb au théorème de Gauss . . . . . . . . . . . . . . 7
1.3 Invariances et symétries . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3.1 Invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3.2 Symétries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4 Exemples de calculs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.4.1 Plaque infinie avec charge volumique . . . . . . . . . . . . . . 11
1.4.2 Boule avec charge volumique . . . . . . . . . . . . . . . . . . . 13
1.4.3 Cylindre avec charge volumique . . . . . . . . . . . . . . . . . 14
1.4.4 Plan infini avec charge surfacique . . . . . . . . . . . . . . . . 15
1.4.5 Sphère avec charge surfacique . . . . . . . . . . . . . . . . . . 16
1.4.6 Cylindre avec charge surfacique . . . . . . . . . . . . . . . . . 17
1.4.7 Fils fini et infini . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.5 Potentiel scalaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.6 Principe de superposition et applications . . . . . . . . . . . . . . . . 20
1.6.1 Enoncé du principe . . . . . . . . . . . . . . . . . . . . . . . . 20
1.6.2 Condensateur plan . . . . . . . . . . . . . . . . . . . . . . . . 21
1.6.3 Condensateur sphérique . . . . . . . . . . . . . . . . . . . . . 22
1.6.4 Condensateur cylindrique . . . . . . . . . . . . . . . . . . . . 24
1.7 Relation de passage pour le champ électrostatique . . . . . . . . . . . 26
1.7.1 Discontinuité de la composante normale . . . . . . . . . . . . 26
1.7.2 Continuité de la composante tangentielle . . . . . . . . . . . . 27
1.7.3 Formulations faible et forte de l’électrostatique . . . . . . . . . 28
1.8 Développement multipolaire du potentiel scalaire . . . . . . . . . . . 29
1.9 Dipôle électrostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.9.1 Calcul du champ électrique . . . . . . . . . . . . . . . . . . . 31
1.9.2 Interaction d’un dipôle avec un champ électrique . . . . . . . 31
1.10 Équations de l’électrostatique, existence et unicité du champ électro-
statique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33

2 Magnétostatique 34
2.1 Notion de courant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.2 De la loi de Biot et Savart au théorème d’Ampère . . . . . . . . . . . 36
2.3 Invariances et symétries . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.3.1 Invariances . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.3.2 Symétries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.4 Exemples de calculs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.4.1 Plaque infinie avec courant volumique . . . . . . . . . . . . . . 39
2.4.2 Cylindre avec courant volumique vertical . . . . . . . . . . . . 41
2.4.3 Cylindre avec courant volumique ortho-radial . . . . . . . . . 42
2.4.4 Plan avec courant surfacique . . . . . . . . . . . . . . . . . . . 43
2.4.5 Cylindre avec courant surfacique vertical . . . . . . . . . . . . 44
2.4.6 Cylindre avec courant surfacique ortho-radial . . . . . . . . . 45

2
 2.4.7 Fils fini et infini . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.4.8 Spire de courant et solénoïdes fini et infini . . . . . . . . . . . 47
2.5 Potentiel vecteur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.6 Relation de passage pour le champ magnétostatique . . . . . . . . . . 50
2.6.1 Continuité de la composante normale . . . . . . . . . . . . . . 50
2.6.2 Discontinuité de la composante tangentielle . . . . . . . . . . 51
2.6.3 Formulations faible et forte de la magnétostatique . . . . . . . 52
2.7 Développement multipolaire du potentiel vecteur . . . . . . . . . . . 53
2.8 Dipôle magnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.8.1 Calcul du champ magnétique . . . . . . . . . . . . . . . . . . 55
2.9 Équations de la magnétostatique, existence et unicité du champ ma-
gnétostatique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56

3 Induction 57
3.1 Champ électromoteur . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.2 Loi de Lenz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.3 Induction de Lorenz . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4 Induction de Neumann . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.5 Inductance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

4 Équations de Maxwell 59
4.1 Conservation de la charge . . . . . . . . . . . . . . . . . . . . . . . . 59
4.2 Formulation locale des équations de Maxwell . . . . . . . . . . . . . . 60
4.3 Relations de passage en régime variable . . . . . . . . . . . . . . . . . 62
4.3.1 Relation de passage pour le champ électrique . . . . . . . . . . 62
4.3.2 Relation de passage pour le champ magnétique . . . . . . . . 63
4.4 Des équations de Maxwell à l’équation des ondes . . . . . . . . . . . . 64
4.5 Potentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.6 Modes de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.7 ARQS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
4.8 Aspects énergétiques du champ électromagnétique . . . . . . . . . . . 69
4.9 Formulation des équations de Maxwell comme un problème d’évolution 70

5 Conducteurs 72
5.1 Modèle de Drude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
5.2 Loi d’Ohm locale et effet de peau . . . . . . . . . . . . . . . . . . . . 73
5.3 Conducteur parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
5.4 Plasmas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
5.5 Champs moyennes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
5.6 Milieux diélectriques . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
5.7 Milieux diamagnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . 77

6 Reflexion/refraction/transmission 77
6.1 Vide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
6.2 Plasmas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
6.3 Conducteurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
6.4 Diélectriques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77

3
7 Vers l’optique ondulatoire 77
7.1 Equation eikonale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77

8 Appendice 78
8.1 Analyse vectorielle . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
8.1.1 Opérateurs différentiels . . . . . . . . . . . . . . . . . . . . . . 78
8.1.2 Lemme de Poincaré . . . . . . . . . . . . . . . . . . . . . . . . 80
8.1.3 Théorème de Stokes . . . . . . . . . . . . . . . . . . . . . . . . 81
8.1.4 Théorèmes de transport . . . . . . . . . . . . . . . . . . . . . 81
8.2 Théorie du signal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
8.2.1 Produit de convolution . . . . . . . . . . . . . . . . . . . . . . 82
8.2.2 Transformation de Fourier . . . . . . . . . . . . . . . . . . . . 82
8.2.3 Equation des ondes . . . . . . . . . . . . . . . . . . . . . . . . 83
8.3 Equation de Poisson . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
8.3.1 Existence et unicite de la solution avec condition de Neumann 84
8.3.2 Principe du maximum . . . . . . . . . . . . . . . . . . . . . . 85
8.3.3 Fonction de Green . . . . . . . . . . . . . . . . . . . . . . . . 85
8.4 Equation des ondes . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
8.4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
8.4.2 Fonctions de Green du d’alembertien . . . . . . . . . . . . . . 87

4
1 Électrostatique
1.1 Notion de charge
L’ensemble de la physique moderne repose sur quatre interactions dites fondamen-
tales :
— L’interaction gravitationnelle, qui est responsable de l’attraction à grande
échelle des corps massifs, chaque corps massif possédant une grandeur scalaire
intrinsèque strictement positive, appelée masse, dont l’unité est le kilogramme
kg.
— L’interaction nucléaire faible, qui a lieu entre leptons et quarks à l’échelle du
noyau atomique. Elle explique une partie de la cohésion de la matière.
— L’interaction nucléaire forte, qui a lieu entre quarks et gluons, à très courte
portée, qui assure la cohésion des nucléons.
— L’interaction électromagnétique, qui explique l’interaction entre certaines par-
ticules et un champ dit électromagnétique. Ces dernières particules possèdent
une grandeur scalaire intrinsèque, appelée charge, dont l’unité est le Coulomb
C.
Dans ce paragraphe, nous étudierons comment donner un sens mathématique exploi-
table à cette notion de charge.

On se place dans le cadre classique (c’est-à-dire non quantique et non relativiste), à

un instant donné. Soit M ∈ R3 . Un volume δV de R3 contenant M contient une
certaine charge δQ, obtenue en sommant les charges de tous les porteurs de charges
présents dans δV . On définit alors la densité volumique de charge en M par :
δQ
ρ(M ) = lim .
δV →0 δV

L’unité de ρ est le C · m−3 . Il s’agit d’une grandeur intensive en tant que quotient
de deux grandeurs extensives. Ainsi, la charge Q contenue dans un volume V peut
s’écrire : Z
Q= ρ(M )dV.
M ∈V

Pour des raisons pratiques, on sera parfois amené à manipuler des densités de charge
surfacique et linéique (autrement dit des charges par unité de surface ou de longueur).

5
On considère une surface S de R3 et un point M ∈ S. On note δS une surface de
S contenant le point M , et δQ la charge contenue dans δS. On définit la densité
surfacique de charge en M par :
δQ
σ(M ) = lim .
δS→0 δS

L’unité de σ est le C · m−2 . Il s’agit d’une grandeur intensive en tant que quotient
de deux grandeurs extensives. Ainsi, la charge Q contenue dans une surface S de S
peut s’écrire : Z
Q= σ(M )dS.
M ∈S

On considère une courbe C de R et un point M ∈ C. On note δℓ un élément de

3

courbe de C contenant le point M , et δQ la charge contenue dans δC. On définit la

densité linéique de charge en M par :
δQ
λ(M ) = lim .
δS→0 δS

L’unité de λ est le C · m−1 . Il s’agit d’une grandeur intensive en tant que quotient
de deux grandeurs extensives. Ainsi, la charge Q contenue dans une courbe C de C
peut s’écrire : Z
Q= λ(M )dℓ.
M ∈C

On retiendra les égalités symboliques suivantes :

δQ(M ) = ρ(M )dV = σ(M )dS = λ(M )dℓ.

On conclut ce paragraphe en donnant un sens mathématique précis à la notion de

densité de charge sur un volume V de R3 .

Une densité de charge sur V est une distribution ρ ∈ D′ (V, R). Elle est dite :
— volumique si ρ est une fonction de L1loc (V, R) ;
— surfacique s’il existe une surface S de V et une fonction σ ∈ L1loc (S, R) telles
que :
ρ = σdS
avec dS la mesure de surface de S ;
— linéique s’il existe une courbe C de V et une fonction λ ∈ L1loc (C, R) telles
que :
ρ = λdℓ
avec dℓ la mesure de longueur de C.

6
1.2 De la loi de Coulomb au théorème de Gauss
On rappelle que, d’après la loi de Coulomb, le champ électrique créé par une charge
(P, q) en un point M ∈ R3 est :

1 PM
E(M ) = q .
4πε0 |PM|3

On cherche maintenant à déterminer le champ électrique créé par une densité

volumique de charges ρ. On postule que le champ électrique vérifie le principe de
superposition. Étant donné un point P de l’espace, il porte la charge ρ(P )dV où
dV est un volume élémentaire entourant P . D’après la loi de Coulomb, le champ
électrique créé par le point (P, ρ(P )dV ) évalué au point M est :

ρ(P )dV PM
dEP (M ) = .
4πε0 |PM|3

Ainsi, d’après le principe de superposition (dont on reparlera plus tard), le champ

électrique total créé par ρ au point M est :
Z
1 Z PM
E(M ) = dEP (M ) = ρ(P ) dV
P ∈R 3 4πε0 P ∈R 3
|PM|3

Cette forme explicite du champ électrique est cependant peu propice à un calcul
explicite de E. Elle servira plus tard pour établir certaines preuves générales, mais
ne sera pas utile pour réaliser des calculs effectifs.

Une première étape dans la réalisation de calculs effectifs consiste à trouver une
équation plus simple (mais plus implicite) liant E et ρ. On commence par observer
que le champ électrique créé par une charge ponctuelle décroît en inverse du carré de
la distance, ce qui traduit une conservation du flux du champ électrique à travers
les sphères centrées en la charge dans le cas d’une charge ponctuelle. Ainsi, le flux
du champ électrique apparaît comme une quantité intéressante à calculer. On note
toujours E le champ électrique créé par ρ, et on se donne V un volume dont on note
∂V la surface de bord, le vecteur normal n étant orienté vers l’extérieur de V :
Z
1 Z Z
PM
E · ndS = ρ(P ) · n(M )dV dS
∂V 4πε0 M ∈∂V P ∈R 3
|PM|3
1 Z Z
PM
= ρ(P ) · n(M )dSdV
4πε0 P ∈R M ∈∂V
3
|PM|3
1 Z
!
Z
PM
= ρ(P ) 3 · n(M )dS dV.
4πε0 P ∈R3 M ∈∂V |PM|

7
On admet la formule de la moyenne suivante :

Z
PM 4π si P ∈ V
· n(M )dS = 0 sinon
,
M ∈∂V |PM|3

ce qui donne : Z
Q
E · ndS =
∂V ε0
Z
où Q = ρ(P )dV est la charge contenue dans V . Il s’agit du théorème de Gauss.
P ∈V

Le prochain paragraphe est destiné à la notion d’invariance et de symétrie : le membre

de gauche du théorème de Gauss est considérablement simplifié dans des situations
à haute symétrie.

8
1.3 Invariances et symétries
1.3.1 Invariances
On se donne R une transformation orthogonale affine, c’est-à-dire une application
affine R : R3 → R3 dont la partie linéaire R est dans O(3).

La distribution de charges est dite invariante sous la transformation R si :

∀M ∈ R3 , ρ(M ) = ρ(R(M )).

Effectuons le changement de variable R dans l’expression explicite du champ élec-

trique E créé par ρ :

1 Z R(PM)
∀M ∈ R , E(R(M )) =
3
ρ(R(P )) dV
4πε0 P ∈R 3 |R(P )R(M )|3

Or ∀P ∈ R3 , ρ(R(P )) = ρ(P ) et, puisque la transformation R conserve les distances,

∀P, M ∈ R3 , |R(P )R(M )| = |PM|, donc :

∀M ∈ R3 , E(R(M )) = R(E(M )),

donc le champ électrique est invariant sous la transformation R.

Voyons deux exemples majeurs :

— On suppose disposer d’un repère cartésien (O, ux , uy , uz ), et que ρ est invariante
par translation d’axe z (autrement dit, ρ ne dépend que des coordonnées x, y).
Alors E ne dépend que de x, y.
— On suppose disposer de coordonnées polaires (r, θ, φ) de centre O, et que ρ
est invariant par rotation autour de O (autrement dit, ρ ne dépend que de la
coordonnée r). Alors E ne dépend que de r.

9
1.3.2 Symétries
Un plan affine Π de R3 est un plan de symétrie de ρ si :

∀M ∈ R3 , ρ(SΠ (M )) = ρ(M )

où SΠ est la réflexion orthogonale par rapport à Π. Le calcul précédent donne que :

∀M ∈ R3 , E(SΠ (M )) = SΠ (E(M )).

Ainsi, Π est également un plan de symétrie pour le champ E.

Remarquons que, dans le cas où M ∈ Π :

E(M ) = SΠ (E(M )),

donc E(M ) appartient à la direction de Π.

Un plan affine Π⋆ de R3 est un plan d’antisymétrie de ρ si :

∀M ∈ R3 , ρ(SΠ⋆ (M )) = −ρ(M )

où SΠ⋆ est la réflexion orthogonale par rapport à Π⋆ . Le calcul précédent donne que :

∀M ∈ R3 , E(SΠ⋆ (M )) = −SΠ⋆ (E(M )).

Ainsi, Π⋆ est également un plan d’antisymétrie pour le champ E.

Remarquons que, dans le cas où M ∈ Π⋆ :

E(M ) = −SΠ⋆ (E(M )),

donc E(M ) est orthogonal à la direction de Π⋆ .

10
1.4 Exemples de calculs
1.4.1 Plaque infinie avec charge volumique
On munit R3 d’un repère cartésien (O, ux , uy , uz ).

On considère la densité volumique de charge ρ définie par :

ρ = ρ0 1R2 ×[−ℓ,ℓ]

avec ρ0 ∈ R et ℓ > 0. Autrement dit, on considère une plaque infinie de demi-épaisseur

ℓ uniformément chargée.

Invariances : ρ est invariante par translations de vecteurs ux , uy , donc E est invariant

par translations de vecteurs ux , uy .

Symétries : soit M ∈ R3 . Les plans (M, ux , uz ), (M, uy , uz ) sont des plans de

symétrie de ρ, donc :

E(M ) ∈ Vect(ux , uz ) ∩ Vect(uy , uz ) = Ruz .

De plus, le plan (O, ux , uy ) est un plan de symétrie de ρ, donc, pour tout (x, y) ∈ R2 ,
z 7→ E(x, y, z) est impaire.

Ainsi, il existe une fonction E : R → R impaire telle que :

∀x, y, z ∈ R, E(x, y, z) = E(z)uz .

On se donne z ⩾ ℓ.

On considère le volume V = [0, L] × [0, L] × [z, z + L] avec L > 0, de sorte que :

Z
E · ndS = (E(z + L) − E(z))L2 .
∂V

De plus, ρ|V = 0, donc, d’après le théorème de Gauss :

E(z + L) − E(z) = 0.

En faisant tendre L vers 0, on trouve E ′ (z) = 0. On en déduit que E est constante

sur le demi-espace z ⩾ ℓ.

11
On se donne z ∈]0, ℓ].

On considère le volume V = [0, L] × [0, L] × [0, z] avec L > 0, de sorte que :

Z
E · ndS = (E(z) − E(0))L2 .
∂V
Z
De plus, ρdV = ρ0 L2 z, donc, d’après le théorème de Gauss :
V

ρ0
E(z) − E(0) = z.
ε
ρ0
Par imparité de E, E(0) = 0, donc : E(z) = z.
ε0

En appliquant à nouveau l’imparité de E, on trouve :

 ρ0

 z si |z| ⩽ ℓ



 ε0





ρ0
ℓ si z ⩾ ℓ

∀z ∈ R, E(z) = .



ε0



ρ0


ℓ si z ⩽ −ℓ

−



ε0

ρ0 E(z)
ℓ
ε0

z
−ℓ ℓ

ρ0
− ℓ
ε0

12
1.4.2 Boule avec charge volumique
On munit R3 de coordonnées sphériques (r, θ, φ). On considère la densité volumique
de charge ρ définie par : ρ = ρ0 1B(O,R) , avec ρ0 ∈ R et R > 0. Autrement dit, on
considère une boule de rayon R uniformément chargée.

Invariances : ρ est invariante par rotations de centre O, donc E est invariant par
rotations de centre O.

Symétries : soit M ∈ R3 . Les plans (M, ur (M ), uθ (M )), (M, ur (M ), uφ (M )) sont

des plans de symétrie de ρ, donc : E(M ) ∈ Rur (M ).

Ainsi, il existe une fonction E : R∗+ → R telle que :

∀(r, θ, φ) ∈ R∗+ × [0, π] × [0, 2π[, E(r, θ, φ) = E(r)ur (θ, φ).

Soit (r, θ, φ) ∈ R∗+ × [0, π] × [0, 2π[.

Z
On considère le volume V = B(O, r), de sorte que : E · ndS = E(r)4πr2 .
∂V

Z
4
De plus : ρdV = ρ0 π min(R, r)3 , donc, d’après le théorème de Gauss :
V 3
ρ0
E(r)r2 = min(R, r)3 .
3ε0
Ainsi :  ρ0
 r si r ⩽ R
3ε0





∀r ∈ R+ , E(r) = .
3
ρ0 R


si r ⩾ R



3ε0 r2


E(r)
ρ0
R
3ε0

r
R

13
1.4.3 Cylindre avec charge volumique
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité volumique
de charge ρ définie par : ∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, ρ(r, θ, z) = ρ0 1[0,R] (r), avec
ρ0 ∈ R et R > 0. Autrement dit, on considère un cylindre de rayon R uniformément
chargé.

Invariances : ρ est invariante par rotation d’axe (Oz) et par translation de vecteur
uz , donc E est invariant par rotation d’axe (Oz) et par translation de vecteur uz .

Symétries : soit M ∈ R3 . Les plans (M, ur (M ), uθ (M )), (M, ur (M ), uz ) sont des

plans de symétrie de ρ, donc : E(M ) ∈ Rur (M ). Ainsi, il existe une fonction
E : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, E(r, θ, z) = E(r)ur (θ).

On considère le volume V = D(O, r) × [0, L] avec r, L > 0, de sorte que :

Z
E · ndS = E(r)2πrL.
∂V
Z
De plus : ρdV = ρ0 π min(R, r)2 L, donc, d’après le théorème de Gauss :
V

ρ0
E(r)r = min(R, r)2 .
2ε0
Ainsi :  ρ0
 r si r ⩽ R
2ε0





∀r ∈ R∗+ , E(r) = .
2
ρ0 R


si r ⩾ R



2ε0 r


E(r)
ρ0
R
2ε0

r
R

14
1.4.4 Plan infini avec charge surfacique
On munit R3 d’un repère cartésien (O, ux , uy , uz ).

On considère la densité surfacique de charge σ uniforme de valeur σ0 ∈ R portée par

la surface R2 × {0} (en coordonnées cartésiennes). Autrement dit, on considère un
plan infini uniformément chargé.

La même étude des invariances et des symétries que dans le cas volumique donne
l’existence d’une fonction E : R → R impaire telle que :

∀(x, y, z) ∈ R3 , E(x, y, z) = E(z)uz .

En procédant comme dans le cas volumique, on montre que E est constante sur
R2 × R∗+ et R2 × R∗− .

Soit L > 0. On applique le théorème de Gauss au volume [−L, L]3 , ce qui donne :
σ0 2
(E(L) − E(−L))L3 = L,
ε0
σ0
donc, par imparité de E : E(L) = . Ainsi :
2ε0
σ0
∀z ∈ R, E(z) = sgn(z).
2ε0

σ0 E(z)
2ε0

σ0
−
2ε0

Le champ électrique n’est pas défini sur la plaque : c’est un défaut de modélisation
du fait du caractère surfacique de la densité de charge. Graphiquement, on retrouve
σ0
le même résultat en prenant ℓ → 0 et ρ0 ℓ → .
2

15
1.4.5 Sphère avec charge surfacique
On munit R3 de coordonnées sphériques (r, θ, φ).

On considère la densité surfacique de charge σ uniforme de valeur σ0 ∈ R portée par

S(O, R) avec R > 0. Autrement dit, on considère une sphère de rayon R uniformément
chargée.

La même étude des invariances et des symétries que dans le cas volumique donne
l’existence d’une fonction E : R∗+ → R telle que :

∀(r, θ, φ) ∈ R∗+ × [0, π] × [0, 2π[, E(r, θ, φ) = E(r)ur (θ, φ).

On applique le théorème de Gauss au volume B(O, r) avec r > 0 :

σ0
E(r)4πr2 = 4πR2 1[R,+∞[ (r),
ε0
donc :
0 si r < R






∀r ∈ R∗+ , E(r) = .
 σ0 R2
si r ⩾ R



ε0 r 2


E(r)
σ0
ε0

r
R

Le champ électrique n’est pas défini sur la sphère : c’est un défaut de modélisation
du fait du caractère surfacique de la densité de charge.

16
1.4.6 Cylindre avec charge surfacique
On munit R3 de coordonnées cylindriques (r, θ, z).

On considère la densité surfacique de charge σ uniforme de valeur σ0 ∈ R portée par

D(0, R) × R avec R > 0 (en coordonnées cartésiennes). Autrement dit, on considère
un cylindre de rayon R dont la surface est uniformément chargée.

La même étude des invariances et des symétries que dans le cas volumique donne
l’existence d’une fonction E : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, E(r, θ, z) = E(r)ur (θ).

On applique le théorème de Gauss au volume D(O, r) × [0, L] avec r, L > 0 :

σ0
E(r)2πrL = 2πRL1[R,+∞[ (r),
ε0
donc : 

 0 si r < R
∀r ∈ R∗+ , E(r) =  σ0 R .
 si r ⩾ R
ε0 r

E(r)
σ0
ε0

r
R

Le champ électrique n’est pas vraiment défini sur le cylindre : c’est un défaut de
modélisation du fait du caractère surfacique de la densité de charge.

17
1.4.7 Fils fini et infini
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité de charge
ρ suivante :
ρ = λ0 dC
avec λ0 ∈ R et C = {0}2 × [−ℓ, ℓ] (où ℓ > 0). Autrement dit, on considère un fil
rectiligne de longueur 2ℓ uniformément chargé.

Soit M ∈ R3 , de coordonnées cylindriques (r, θ, z). Le champ électrique créé au point

M par un point P du fil est :
1 PM
dEP (M ) = λ0 dzP .
4πε0 |PM|3

On note zP l’altitude de P , de sorte que PM = OM − OP = rur (θ) + (z − zP )uz .

D’après le principe de superposition, le champ E total créé au point M est :

λ0 Z ℓ rur (θ) + (z − zP )uz

E(M ) = dzP .
4πε0 −ℓ (r2 + (z − zP )2 )3/2

Or :
Z ℓ
r 1 Z (z+ℓ)/r 1
dzP = 2 du
−ℓ (r + (z − zP ) )
2 2 3/2 r (z−ℓ)/r (1 + u2 )3/2
#(z+ℓ)/r
1
"
u
= 2 √
r 1 + u2 (z−ℓ)/r
 
 z +ℓ z−ℓ 
1
 
= 2
 r r 
!2 − v
v 
r u
z+ℓ
!2 
t1 + z − ℓ
u u 
1+
t u 
r r

et :
Z ℓ
z − zP 1 Z (z+ℓ)/r u
dzP = du
−ℓ (r + (z − zP ) )
2 2 3/2 r (z−ℓ)/r (1 + u2 )3/2
#(z+ℓ)/r
1 1
"
= −√
r 1 + u2 (z−ℓ)/r
 
 
1 1 1
 
= 

v − v .
r u 2 2
t1 + z + ℓ
! !
z−ℓ
u 
1+
u
t u 
r r

18
Ainsi :
 
 z+ℓ z−ℓ 
λ0 
 
E(M ) = v r − r  ur (θ)

4πε0 r  u
v
2 2
z+ℓ
u ! !
z−ℓ 
u
1+ 1+
t u 
t
r r
 
 
λ0  1 1
 
+

v −v !2  uz

4πε0 r 
u z+ℓ
!2 u
z −ℓ 
ut1 + t1 +
u 
r r

Dans le cas où ℓ → +∞, on trouve le champ électrique créé par un fil infini
uniformément chargé :
λ0
E(M ) = ur (θ).
2πε0 r
On retrouve le même résultat pour le fil infini en utilisant le théorème de Gauss
(alors qu’on manque de symétries pour l’appliquer dans le cas du fil fini).

En reprenant le cylindre infini de rayon R portant une charge volumique uniforme

λ0
ρ0 , on retrouve, lorsque R → 0 et ρ0 R2 → , le fil infini de charge linéique λ0 .
π

19
1.5 Potentiel scalaire
Le but de ce paragraphe est de donner du sens à la notion de potentiel, d’où découle
la notion de tension. On commence par reprendre l’expression explicite de E pour
trouver :
∇ ∧ E = 0.
Ainsi, d’après le lemme de Poincaré, il existe une fonction numérique lisse sur R3
telle que :
E = −∇Φ.
Une telle fonction Φ est appelée potentiel scalaire. Une telle fonction n’est pas unique.
En effet, si Φ0 est une constante et si Φ est un potentiel scalaire de E, alors Φ + Φ0
est encore un potentiel scalaire de E.

Enfin, on remarque que l’expression explicite de E permet de trouver un potentiel

scalaire explicite Φ défini par :
1 Z 1
∀M ∈ R3 , Φ(M ) = ρ(P ) dV.
4πε0 P ∈R3 |PM|

1.6 Principe de superposition et applications

1.6.1 Enoncé du principe
On rappelle l’expression du champ électrique E créé par une densité de charges ρ :
1 Z PM
∀M ∈ R3 , E(M ) = ρ(P ) dV.
4πε0 P ∈R3 |PM|3

On remarque que cette expression est linéaire en ρ. Ainsi, si ρ1 , ρ2 sont deux distribu-
tions de charge donnant des champs électriques E1 , E2 , alors le champ électrique créé
par ρ1 + ρ2 est E1 , E2 . Cette linéarité est aussi valide pour les potentiels scalaires : si
Φ1 , Φ2 sont deux potentiels scalaires créés par ρ1 , ρ2 , alors Φ1 + Φ2 est un potentiel
scalaire créé par ρ1 + ρ2 . De façon symbolique, cette linéarité est aussi valide en
superposant une infinité (potentiellement « continue » dans le cas d’une intégrale)
de densités de charges.

Cette remarque anodine permet de calculer séparément les champs électriques créés
par certaines parties de l’espace (en appliquant le théorème de Gauss par exemple),
le champ électrique total étant obtenu en sommant les champs électriques calculés.

20
1.6.2 Condensateur plan
On se donne un repère cartésien (O, ux , uy , uz ).

Un condensateur plan est modélisé par la distribution de charge ρ :

ρ = σ0 dS+ − σ0 dS−

avec S± = R2 × {±ℓ} où ℓ ∈ R∗+ est la demi-épaisseur du condensateur et σ0 ∈ R

est la densité surfacique de charge de la plaque supérieure.

Autrement dit, on superpose une plaque en z = ℓ de charge surfacique σ0 , et une

plaque en z = −ℓ de charge surfacique −σ0 .

Le principe de superposition donne le champ électrique créé par la superposition des

deux plaques :

0 si |z| ⩾ ℓ





∀(x, y, z) ∈ R3 , E(x, y, z) =  .
σ0
− uz si |z| < ℓ



ε0
Notons que E est dirigé des charges positives vers les charges négatives.

On peut aussi calculer un potentiel Φ :

 σ0

 ℓ si z ⩾ ℓ



 ε0





σ0
z si |z| < ℓ

∀(x, y, z) ∈ R3 , Φ(x, y, z) =  .


ε0



σ0


ℓ si z ⩽ −ℓ

−



ε0

σ0 Φ(0, 0, z)
ℓ
ε0

z
−ℓ ℓ

σ0
− ℓ
ε0

21
On peut retrouver σ0 à partir d’une différence de potentiel entre les deux plaques :
ε0
σ0 = U
2ℓ
où U = Φ(0, 0, ℓ) − Φ(0, 0, −ℓ). On suppose que les deux plaques ont une surface S
finie suffisamment importante pour que les expressions précédentes soient approxi-
mativement vraies.

On note Q = σ0 L2 la charge portée par la plaque supérieure, de sorte que :

Q ≃ CU
ε0 S
avec C = la capacité du condensateur.
2ℓ

1.6.3 Condensateur sphérique

On munit R3 de coordonnées sphériques (r, θ, φ). On considère la densité surfacique
de charges définie par :
R− 2
ρ = σ0 dS− − σ0 dS+
R+ 2
avec S± = S(O, R± ) σ0 ∈ R et 0 < R− < R+ . Autrement dit, il s’agit de deux
sphères concentriques de rayons R− , R+ portant des charges surfaciques, de charges
totales opposées. Une telle distribution de charges est qualifiée de condensateur
sphérique.

Le principe de superposition donne :

si r ⩽ R− ou r ⩾ R+



0


∀(r, θ, φ) ∈ R∗+ × [0, π] × [0, 2π[, E(r, θ, φ) =  .
σ0 R− 2
ur (θ, φ) si R− < r < R+




ε0 r 2
Notons que E est dirigé des charges positives vers les charges négatives.

On en déduit également un potentiel :

0 si r ⩾




R+





2
1 1
 !
 σ0 R−


si r ∈]R− , R+ [

−
∀(r, θ, φ) ∈ R∗+ × [0, π] × [0, 2π[, Φ(r, θ, φ) = ε0 r R+ .






 σ0 R− 2 1 1

 !


 − si r ⩽ R−
ε0 R− R+


22
 Φ(r, 0, 0)

σ0 R− 2 1 1
!
−
ε0 R− R+

r
R− R+

On peut retrouver σ0 à partir de la différence de potentiel entre les deux sphères :

ε0
σ0 = !U
1 1
R− 2 −
R− R+

où U = Φ(R− , 0, 0) − Φ(R+ , 0, 0). On note Q = 4πR− 2 σ0 la charge de la sphère

intérieure, de sorte que :
Q = CU
4πε0
avec C = 1 1 la capacité du condensateur.
−
R− R+

23
1.6.4 Condensateur cylindrique
On munit R3 de coordonnées cylindriques (r, θ, z).

On considère la densité surfacique de charges ρ définie par :

R−
ρ = σ0 dS− − σ0 dS+
R+

avec S± = D(O, R± ) × R, σ0 ∈ R et 0 < R− < R+ . Autrement dit, il s’agit de deux

cylindres co-axiaux portant des charges surfaciques uniformes opposées. Une telle
distribution de charges est qualifiée de condensateur cylindrique.

Le principe de superposition donne :




 0 si r ⩽ R− ou r ⩾ R+


∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, E(r, θ, z) =  .
σ0 R−
ur (θ) si R− < r < R+



ε0 r


Notons que E est dirigé des charges positives vers les charges négatives.

On en déduit également un potentiel :




0 si r ⩾ R+






 σ0 R− R+

  
ln si R− < r < R+

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, Φ(r, θ, z) = ε0 r .





 !
σ0 R− R+


ln si r ⩽ R−





ε0 R−

24
 Φ(r, 0, 0)
!
σ0 R− R+
ln
ε0 R−

r
R− R+

On peut retrouver σ0 à partir de la différence de potentiel entre les deux cylindres :

ε0
σ0 = !U
R+
R− ln
R−

où U = Φ(R− , 0, 0) − Φ(R+ , 0, 0). On suppose que les deux cylindres ont une
hauteur H finie suffisamment importante pour que les expressions précédentes soient
approximativement vraies.

On note Q = σ0 2πR− H la charge portée par la plaque intérieure, de sorte que :

Q ≃ CU
2πHε0
avec C = ! la capacité du condensateur.
R+
ln
R−

25
1.7 Relation de passage pour le champ électrostatique
Dans le cas de la plaque infiniment fine uniformément chargée, on a vu que la
σ
composante normale à la plaque du champ électrique subissait un saut de , tandis
ε0
que sa composante tangentielle était continue à travers la plaque.

Il est possible de généraliser cette identité à une distribution de charge surfacique

(S, σ) quelconque.

Étant donné M ∈ S, on note E(M ± ) = E(M + 0± n(M )).

1.7.1 Discontinuité de la composante normale

On considère un point M ∈ S, et un cylindre Σ de hauteur h et de rayon r, d’axe
dirigé par la normale n(M ) à S en M .

Écrivons le flux de E à travers Σ. On note Σ± les deux faces supérieure et inférieure

de Σ et Σlat sa surface latérale. On a, en notant ν la normale à Σ :
Z
E · νdS = E(M ± ) (±n(M )) πr2 + O(r3 ) + O(h)
Σ± (r,h)→(0,0)

et : Z
E · νdS = O(h),
Σlat (r,h)→(0,0)

donc :
Z  
E · νdS = E(M + ) − E(M − ) · n(M )πr2 + O(r3 ) + O(h).
Σ (r,h)→(0,0)

Or la charge enfermée dans le volume de bord Σ est :

σ(M )πr2 + O(r3 ),

donc, d’après le théorème de Gauss :

  σ(M ) 2
E(M + ) − E(M − ) · n(M )r2 = r + O(r3 ) + O(h).
(r,h)→(0,0) ε0

On fait d’abord tendre h vers 0, puis on divise par r2 et on fait tendre r vers 0 dans
la relation précédente, ce qui donne :
  σ(M )
E(M + ) − E(M − ) · n(M ) = .
ε0
Il s’agit de la relation de passage pour la composante normale du champ électrique.

26
1.7.2 Continuité de la composante tangentielle
On remarque que l’expression explicite du champ électrique donne :

∇ ∧ E = 0,

ce qui donne, d’après le théorème de Stokes, pour toute surface Σ de R3 :

Z
E · τ dℓ = 0
∂Σ

où τ est le vecteur tangent à ∂Σ.

On considère un point M ∈ S, t un vecteur tangent à S en M , et un cadre

rectangulaire C dont deux cotés sont portés par t. Les deux cotés portés par t sont
de longueur r, et les deux autres cotés ont pour longueur h. On note C ± les deux
cotés portés par t, et Clat l’union des deux autres cotés. On a, en notant τ le vecteur
tangent à C :
Z
E · τ dℓ = ±E(M ± )r + O(r2 ) + O(h)
C± (r,h)→(0,0)

et : Z
E · τ dℓ = O(h),
Clat (r,h)→(0,0)

donc :
Z  
E · τ dℓ = E(M + ) − E(M − ) · tr + O(r2 ) + O(h).
C (r,h)→(0,0)

On fait tendre h vers 0, puis on divise par r et on fait tendre r vers 0 :

 
E(M + ) − E(M − ) · t = 0.

La relation précédente étant vraie pour tout vecteur t tangent à S en M :

 
E(M + ) − E(M − ) ∧ n(M ) = 0.

On a donc continuité de la composante tangentielle de E à travers S.

27
1.7.3 Formulations faible et forte de l’électrostatique
Ce paragraphe nécessite des rudiments de théorie des distributions. Il peut être omis
en première lecture.

Soit V un volume de R3 et ρ ∈ D′ (V, R). On appelle équations de l’électrostatique le

système : 
ρ
∇ · E =

ε0
∇ ∧ E = 0


où l’inconnue est E ∈ D′ (V, R3 ).

Soit S une surface de V telle que V \ S est l’union disjointe de deux ouverts de R3
V± tels que le vecteur n normal a S soit dirige vers V+ .

On suppose ρ de la forme ρ = ρ0 + σdS avec ρ0 localement intégrable sur V , dS la

mesure de surface de S, et σ : S → R localement intégrable.
 
Soit E : V → R3 tel que E|V± admet un prolongement de classe E± ∈ C 1 V± , R3 .
On a équivalence entre :

le champ E vérifie, au sens faible sur V :


ρ
∇ · E
 =
ε0
∇ ∧ E = 0


et :

les champs E± vérifient, sur V± :


ρ0
∇ · E±
 =
ε0
∇ ∧ E± =0


et :
  σ
 E+ |S
 − E− |S · n =
  ε0 .
 E+

|S − E− |S ∧ n =0

Ainsi, la relation de passage pour le champ électrostatique est encodée dans les
équations de l’électrostatique sous forme faible.

28
1.8 Développement multipolaire du potentiel scalaire
Le but de cette section est d’étudier le potentiel scalaire à longue distance créé par
une densité volumique de charge à support compact ρ. On admet l’égalité suivante,
qui peut se prouver à l’aide du théorème des résidus :
1 +∞
= Pℓ (x)z ℓ
X
√
1 − 2xz + z 2 ℓ=0
avec (Pℓ )ℓ∈N la suite des polynômes de Legendre. On se donne des coordonnées
sphériques (r, θ, φ) d’origine O. On a :
1 1 +∞
X |OM|ℓ OP · OM
!
=q = ℓ+1
Pℓ .
|PM| |OM|2 + OP2 − 2OP · OM ℓ=0 |OP| |OP||OM|

Ainsi, en remplaçant dans le potentiel scalaire :

1 +∞ 1
!
OM
∀M ∈ R \ {O}, Φ(M ) =
3
X
Mℓ
4πε0 ℓ=0 |OM|ℓ+1 |OM|
avec : !
Z
OP
∀ℓ ∈ N, ∀r ∈ S , Mℓ (r) =
2
|OP| Pℓ · r ρ(P )dV
ℓ
P ∈R 3 |OP|
le ℓ-ième moment multipolaire. On peut calculer les trois premiers moments multipo-
laires. On rappelle que :
3 1
P0 = 1, P1 = X et P2 = X 2 − ,
2 2
donc :
∀r ∈ S2 , M0 (r) = Q, M1 (r) = r · p et M2 (r) =
X
Qij ri rj
i,j∈{x,y,z}
Z Z
où Q = ρ(P )dV est la charge totale, p = OPρ(P )dV est le moment
P ∈R3 P ∈R3
dipolaire par rapport à O (indépendant de O si Q = 0), et :
3 1
Z  
Qij = OPi OPj − δij |OP|2 ρ(P )dV
P ∈R3 2 2
est le moment quadripolaire. On a la description suivante :
— Si Q ̸= 0, alors :
Q
Φ(M ) ∼ .
M →∞ 4πε0 |OM|

On retrouve à grande distance le potentiel créé par une charge ponctuelle Q en

O.
— Si Q = 0 et p ̸= 0, alors
p · OM
Φ(M ) ∼ .
M →∞ 4πε0 |OM|3
Il s’agit par exemple du potentiel scalaire créé par une molécule d’eau à grande
distance.

29
 — Si Q = 0, p = 0 et Q =
̸ 0, alors :

1
! !
OM OM
Φ(M ) ∼
X
Qij .
M →∞ 4πε0 |OM|3 i,j∈{x,y,z}
|OM| i
|OM| j

Il s’agit par exemple du champ créé par un carré de quatre charges deux à
deux opposées (d’où le nom de quadripôle), ou encore par une molécule de
dioxyde de carbone.
On peut développer les moments multipolaires en utilisant les harmoniques sphériques.
On rappelle l’identité suivante :

4π X ℓ
∀ℓ ∈ N, ∀x, y ∈ S2 , Pℓ (x · y) = Yℓ,m (x)Yℓ,m (y)⋆ .
2ℓ + 1 m=−ℓ

On peut ainsi écrire sous une forme plus explicite la dépendance angulaire :

1 +∞ ℓ
1
!
OM
∀M ∈ R \ {O}, Φ(M ) =
3
X X
Mℓ,m Yℓ,m
4πε0 ℓ=0 m=−ℓ |OM|ℓ+1 |OM|
avec :
!⋆
4π Z OP
2
∀ℓ ∈ N, ∀m ∈ J−ℓ, ℓK, ∀r ∈ S , Mℓ,m = |OP|ℓ Yℓ,m ρ(P )dV
2ℓ + 1 P ∈R3 |OP|

les moments multipolaires sphériques.

30
1.9 Dipôle électrostatique
On considère une distribution de charge ρ à support compact. L’objectif de ce
paragraphe est d’étudier Zle champ électrique à grande distance créé par ρ dans le cas
où la charge totale Q = ρ(P )dV est nulle. On choisit une origine O.
P ∈V

1.9.1 Calcul du champ électrique

D’après le paragraphe précédent, le potentiel scalaire s’écrit :
p · OM
Φ(M ) ∼ ,
M →∞ 4πε0 |OM|3
donc le champ électrique vérifie :
3(p · OM)OM − |OM|2 p
E(M ) ∼ .
M →∞ 4πε0 |OM|5

1.9.2 Interaction d’un dipôle avec un champ électrique

On place la distribution de charges ρ (de charge totale toujours supposée nulle) en
présence d’un champ électrique extérieur E.

On note V + = {P ∈ V, ρ(P ) > 0} et V − = {P ∈ V, ρ(P ) < 0}, de sorte à définir les

barycentres des charges positives et négatives P + , P − tels que :
Z Z
ρ(P )P PdV =
+
ρ(P )P− PdV = 0
P ∈V + P ∈V −

Calculons l’énergie potentielle du dipôle. On a :

Z
Ep = ρ(P )Φ(P )dV
P ∈V

Φ(P ) = Φ(P + ) + ∇Φ(P + ) · P+ P + o(|P+ P|),

donc :
Z Z
ρ(P )Φ(P )dV = Φ(P + ) ρ(P )dV
P ∈V + P ∈V +
Z
+ ∇Φ(P + ) · ρ(P )P+ PdV + o(diam(V ))
P ∈V +
Z
= Φ(P ) +
ρ(P )dV + o(diam(V )).
P ∈V +

De même :
Z Z
ρ(P )Φ(P )dV = Φ(P ) −
ρ(P )dV + o(diam(V ),
P ∈V − P ∈V −

donc :

Ep = q(Φ(P + ) − Φ(P − )) + o(diam(V )) = −qP− P+ · E0 + o(diam(V ))

31
 Z
où q = ρ(P )dV est la charge positive et E0 est le champ électrique évalué en
P ∈V +
un point de V . De plus, par définition de p :
Z Z Z
p= ρ(P )P PdV =
−
ρ(P )P PdV =
−
ρ(P )P− P+ dV = qP− P+ ,
P ∈V P ∈V + P ∈V +

donc :
Ep = −p · E0 + o(diam(V )).
On en déduit que le dipôle subit une force F :

F = ∇(p · E0 ) + o(diam(V )).

Ainsi, lorsque le dipôle est placé dans un champ électrique uniforme, la résultante
des forces qui s’exercent sur lui est nulle. Sinon, le dipôle est attiré vers les régions
de champ E intense, et a tendance à s’aligner dans la même direction que E0 .

Calculons le moment résultant MO du champ électrique sur le dipôle où O est une

origine arbitraire. On a :
Z
MO = OP ∧ ρ(P )E(P )dV
ZP ∈V
= OPρ(P )dV ∧ E0 + o(diam(V ))
P ∈V
= p ∧ E0 + o(diam(V )).

On remarque que le moment ne dépend pas du point O.

32
1.10 Équations de l’électrostatique, existence et unicité du
champ électrostatique
Ce paragraphe nécessite des connaissances mathématiques en théorie des distributions.
Il peut être omis en première lecture.

On se donne un volume V borne de R3 , ρ ∈ L2 (V, R) et En ∈ L2 (∂V, R) vérifiant

les relations de compatibilité :
1 Z Z
ρ(M )dV − En (M )dS = 0.
ε0 M ∈V M ∈∂V

On cherche à résoudre le problème :

Trouver E ∈ H 1 (V, R3 ) tel que :

 ρ
 ∇·E =
ε0



 ∇ ∧ E = 0 .

E|∂V · n = En



Le lemme de Poincaré donne :

 Z 
{E ∈ H (V, R ), ∇ ∧ E = 0} = −∇Φ, Φ ∈ H (V, R) et
1 3 2
Φ(M )dV = 0 .
M ∈V

Le théorème d’existence et d’unicité de la solution à l’équation de Poisson avec

condition de Neumann (annexe) donne que le problème précédent admet une unique
solution.

33
2 Magnétostatique
2.1 Notion de courant
Il existe plusieurs façons de modéliser le courant. On présente ici une modélisation
statistique. A l’échelle microscopique, les porteurs de charges (par exemple les
électrons, les ions) ont un mouvement désordonné à cause de l’agitation thermique.
Ce fait nécessite d’introduire une densité de particules et de moyenner les différentes
grandeurs physiques pertinentes selon cette densité. Ainsi, pour toute position
M ∈ R3 , pour toute vitesse v ∈ R3 , pour toute charge q ∈ R, on note :

f (M, t, v, q)

le nombre de particules présentes dans un volume dM dt dv dq de R3 × R × R3 × R

centre en (M, t, v, q). L’unité de f (M, v, q) est m−3 · s−1 · (m · s−1 )−3 · C−1 .

On retrouve la densité volumique de charge avec l’expression :

Z
∀(M, t) ∈ R3 × R, ρ(M, t) = qf (M, t, v, q)dvdq.
R3 ×R

La vitesse de groupe est alors définie par :

Z
∀(M, t) ∈ R × R, v(M, t) =
3
vf (M, t, v, q)dvdq,
R3 ×R

tandis que la densité volumique de courant est définie par :

Z
∀(M, t) ∈ R × R, j(M, t) =
3
qvf (M, t, v, q)dvdq.
R3 ×R

Par exemple, si une seule particule de charge q0 et de position t 7→ M (t) est présente,
alors :
dOM
!!
∀(M, t, v, q) ∈ R ×R×R ×R, f (M, t, v, q) = δ (M, v, q), M (t),
3 3
(t), q0 ,
dt

ce qui donne :
dOM
∀(M, t) ∈ R3 × R, j(M, t) = q0 δ(M, M (t)) (t).
dt
Etant donnée une surface S de R3 de vecteur normal n, l’intensité traversant S est
définie comme étant le flux de j à travers S :
Z
∀t ∈ R, I(t) = j(M, t) · n(M )dS.
S

34
Il est plus délicat de définir une bonne notion de « charge ponctuelle » pour les courants
électriques, dans la mesure où il n’existe pas de courant ponctuel (nous en reparlerons
ultérieurement). Il est cependant possible de définir un objet mathématique analogue,
l’élément de courant, qui permettra de définir les notions de courants surfacique et
linéique. Pour tout (M, t) ∈ R3 × R, on définit l’élément de courant par :
δC(M, t) = j(M, t)dV.
Supposons que l’élément de volume dV est de la forme dS dz avec z ∈] − h, h[ où S
est une surface de R3 et h > 0. On note, si elle existe js (M, t) = lim jh (M, t)dz la
h→0
densité surfacique de courant, de sorte que : δC(M, t) = js (M, t)dS.

Le cas où l’écriture précédente est licite correspond à la présence de courant sur la

surface S, sachant qu’un voisinage de S privé de S est dépourvu de courant.

Supposons que l’élément de surface dS est de la forme dℓ dy avec y ∈] − h, h[, où dℓ

est un élément de longueur élémentaire d’une courbe C de vecteur tangent t et h > 0.
On note, si elle existe, I(M, t) = lim js h (M, t)dy la densité linéique de courant, de
h→0
sorte que : δC(M, t) = I(M, t)dℓ.

Le cas où l’écriture précédente est licite correspond à la présence de courant traversant

la courbe C, sachant qu’un voisinage de C privé de C est dépourvu de courant.

Justifions la notation I. Au voisinage de M , on épaissit la courbe C en un volume

Ch , et on note Sh (M ) la section de Ch passant par M , dont le vecteur normal n
vérifie n(M ) = t(M ).
Z
On a : I(M, t) · t(M ) = lim ISh (M ) (t) où ISh (M ) (t) = j(P, t) · n(P )dP est
h→0 Sh (M )
l’intensité passant par Sh (M ).

On retiendra les égalités symboliques suivantes :

δC(M, t) = j(M, t)dV = js (M, t)dS = I(M, t)dℓ.
On conclut ce paragraphe en donnant un sens mathématique précis à la notion de
densité de courant sur un volume V de R3 .

Une densité de courant sur V est une distribution j ∈ D′ (V, R3 ). Elle est dite :
— volumique si j est une fonction de L1loc (V, R3 ) ;
— surfacique s’il existe une surface S de V et une fonction js ∈ L1loc (S, R3 ) telles
que :
j = js dS
avec dS la mesure de surface de S ;
— linéique s’il existe une courbe C de V et une fonction I ∈ L1loc (C, R3 ) telles
que :
j = Idℓ
avec dℓ la mesure de longueur de C.
35
2.2 De la loi de Biot et Savart au théorème d’Ampère
On rappelle que, d’après la loi de Biot et Savart, le champ magnétique créé en un
point M ∈ R3 par un élément de courant (P, C) est :
µ0 PM
B(M ) = C∧ .
4π |PM|3
On cherche maintenant à déterminer le champ magnétique créé par une densité
volumique de courant j. On postule que le champ magnétique vérifie le principe de
superposition. Étant donne un point P ∈ R3 , il porte un élément de courant j(P )dV ,
où dV est un volume élémentaire entourant P .

D’après la loi de Biot et Savart, le champ magnétique créé par le point (P, j(P )dV )
évalué au point M est :
µ0 PM
dBP (M ) = j(P )dV ∧ .
4π |PM|3
Ainsi, d’après le principe de superposition, le champ magnétique total créé par j au
point M est :
Z
µ0 Z PM
B(M ) = dBP (M ) = j(P ) ∧ dV.
P ∈R3 4π P ∈R3 |PM|3
Comme dans le cas du champ électrique, cette expression est trop complexe pour
calculer B explicitement à partir de j. Les expériences d’Ampère montrent qu’une
quantité intéressante à considérer est la circulation de B le long d’un contour fermé.
On se donne ainsi S une surface orientée. On a :
!
Z
µ0 Z Z
PM
B · tdℓ = j(P ) ∧ dV · t(M )dℓ
∂S 4π M ∈∂S P ∈R3 |PM|3
!
µ0 Z Z
PM
= ∧ t(M )dℓ · j(P )dV
4π P ∈R3 M ∈∂S |PM|3

Or une formule de la moyenne donne :


Z
PM 4πn(P ) si P ∈ S
∧ t(M )dℓ = 0 sinon
,
M ∈∂S |PM|3

donc : Z
B · tdℓ = µ0 I
∂S
Z
où I = j · ndS est l’intensité du courant j à travers S. Il s’agit du théorème
S
d’Ampère.

Le prochain paragraphe est destiné à la notion d’invariance et de symétrie : le membre

de gauche du théorème d’Ampère est considérablement simplifié dans des situations
à haute symétrie.

36
2.3 Invariances et symétries
2.3.1 Invariances
On se donne R une transformation orthogonale affine, c’est-à-dire une application
affine R : R3 → R3 dont la partie linéaire R est dans O(3).

La distribution de courant j est dite invariante sous la transformation R si :

∀M ∈ R3 , j(R(M )) = R(j(M ))

où R est la partie linéaire de R.

Effectuons le changement de variable R dans l’expression explicite du champ magné-

tique B créé par ρ :

µ0 Z R(PM)
∀M ∈ R3 , B(R(M )) = j(R(P )) ∧ dV
4π P ∈R3 R(P )R(M )3
µ0 Z R(PM)
= R(j(P )) ∧ dV.
4π P ∈R3 |PM|3

Or ∀x, y ∈ R3 , R(x) ∧ R(y) = det(R)R(x ∧ y), donc :

B(R(M )) = det(R)R(B(M )),

donc le champ magnétique est invariant sous R si R est directe, sinon le champ
magnétique est changé en son opposé. Cette propriété donne que B est un pseudo-
vecteur.

Voyons deux exemples majeurs :

— On suppose disposer d’un repère cartésien (O, ux , uy , uz ), et que j est invariante
par translation d’axe z (autrement dit, j ne dépend que des coordonnées x, y).
Alors B ne dépend que de x, y.
— On suppose disposer de coordonnées polaires (r, θ, φ) de centre O, et que j
est invariant par rotation autour de O (autrement dit, j ne dépend que de la
coordonnée r). Alors B ne dépend que de r.

37
2.3.2 Symétries
Un plan affine Π de R3 est un plan de symétrie de j si :

∀M ∈ R3 , j(SΠ (M )) = SΠ (j(M ))

où SΠ est la réflexion orthogonale par rapport à Π. Le calcul précédent donne que,

puisque les réflexions sont des isométries indirectes :

∀M ∈ R3 , B(SΠ (M )) = −SΠ (B(M )).

Ainsi, Π est un plan d’antisymétrie pour le champ B.

Remarquons que, dans le cas où M ∈ Π :

B(M ) = −SΠ (B(M )),

donc B(M ) est orthogonal à la direction de Π.

Un plan affine Π⋆ de R3 est un plan d’antisymétrie de j si :

∀M ∈ R3 , j(SΠ⋆ (M )) = −SΠ⋆ (j(M )).

Le calcul précédent donne :

∀M ∈ R3 , B(SΠ⋆ (M )) = SΠ⋆ (B(M )).

Ainsi, Π⋆ est plan de symétrie pour le champ B.

Remarquons que, dans le cas où M ∈ Π⋆ :

B(M ) = SΠ⋆ (B(M )),

donc B(M ) appartient à la direction de Π⋆ .

38
2.4 Exemples de calculs
2.4.1 Plaque infinie avec courant volumique
On munit R3 d’un repère cartésien (O, ux , uy , uz ).

On considère la densité volumique de courants j définie par :

j = j0 1R2 ×[−ℓ,ℓ] ux

avec j0 ∈ R et ℓ > 0. Autrement dit, on considère une plaque infinie de demi-épaisseur

ℓ parcourue par un courant volumique uniforme.

Invariances : j est invariante par translations de vecteurs ux , uy , donc B est invariant

par translations de vecteurs ux , uy .

Symétries : soit M ∈ R3 . Le plan (M, ux , uz ) est un plan de symétrie de j, donc :

B(M ) ∈ Vect(ux , uz )⊥ = Ruy .

De plus, le plan (O, ux , uy ) est un plan de symétrie de j, donc, pour tous (x, y) ∈ R2 ,
z 7→ B(x, y, z) est impaire.

Ainsi, il existe une fonction B : R → R impaire telle que :

∀x, y, z ∈ R, B(x, y, z) = B(z)uy .

On se donne z ⩾ ℓ.

On considère la surface S = {0} × [0, L] × [z, z + L] avec L > 0, de sorte que :

Z
B · tdℓ = (−B(z + L) + B(z))L.
∂S

De plus, comme j|S = 0, on a, d’après le théorème d’Ampère :

B(z + L) − B(z) = 0.

En faisant tendre L vers 0, on trouve B ′ (z) = 0, donc B est constante sur le

demi-espace z ⩾ ℓ.

On se donne z ∈]0, ℓ], et on considère la surface S = {0} × [0, L] × [0, z] avec L > 0,
de sorte que : Z
B · tdℓ = (−B(z) + B(0))L.
∂S

39
 Z
De plus, j · ndS = j0 zL, donc, d’après le théorème d’Ampère :
S

−B(z) + B(0) = µ0 j0 z,

donc, puisque B(0) = 0 par imparité de B : B(z) = −µ0 j0 z. On applique à nouveau

l’imparité de B :
−µ0 j0 z si |z| ⩽ ℓ









∀z ∈ R, B(z) = −µ0 j0 ℓ si z ⩾ ℓ .







µ0 j0 ℓ si z ⩽ −ℓ



B(z)
µ0 j0 ℓ

z
−ℓ ℓ

−µ0 j0 ℓ

40
2.4.2 Cylindre avec courant volumique vertical
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité volumique
de courant j définie par : ∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, j(r, θ, z) = j0 1[0,R[ (r)uz avec
j0 ∈ R et R > 0. Autrement dit, on considère un cylindre de rayon R parcouru par
des courants verticaux.

Invariances : j est invariante par rotation d’axe (Oz) et par translation de vecteur
uz , donc B l’est aussi.

Symétries : soit M ∈ R3 . Le plan (M, ur (M ), uz ) est un plan de symétrie de j,

donc (M, ur (M ), uz ) est un plan d’antisymétrie de B, donc :

B(M ) ∈ Vect(ur (M ), uz )⊥ = Ruθ (M ).

Ainsi, il existe une fonction B : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, B(r, θ, z) = B(r)uθ (θ).

On considère la surface S = D(O, r) × {0} avec r > 0, de sorte que :

Z
B · tdℓ = B(r)2πr.
∂S
Z
De plus : jdS = j0 π min(r, R)2 , donc, d’après le théorème d’Ampère :
S

2B(r)r = µ0 j0 min(r, R)2 .

Ainsi :
µ0 j0

r si r ⩽ R


 2



∀r ∈ R∗+ , B(r) = .
µ j R2


 0 0 si r ⩾ R



2 r

B(r)
µ0 j0
R
2

r
R

41
2.4.3 Cylindre avec courant volumique ortho-radial
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité volumique
de courant j définie par : ∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, j(r, θ, z) = j0 1[0,R[ (r)uθ (θ) avec
j0 ∈ R et R > 0. Autrement dit, on considère un cylindre de rayon R parcouru par
un courant ortho-radial uniforme.

Invariances : j est invariante par rotation d’axe (Oz) et par translation de vecteur
uz , donc B l’est aussi.

Symétries : soit M ∈ R3 . Le plan (M, ur (M ), uθ (M )) est un plan de symétrie de

j, donc (M, ur (M ), uθ (M )) est un plan d’antisymétrie de B, donc B(M ) ∈ Ruz ).
Ainsi, il existe une fonction B : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, B(r, θ, z) = B(r)uz .

On considère la surface (cartésienne) S = {0} × [0, r] × [0, L] avec r, L > 0, de sorte

que : Z
B · tdℓ = (B(r) − B(0))L.
∂S
Z
De plus : j · ndS = j0 L min(r, R), donc, d’après le théorème d’Ampère :
S

B(r) − B(0) = µ0 j0 min(r, R).

Des considérations de symétrie donnent B(0) = 0, ce qui fournit :

si r ⩽ R

µ0 j0 r


∀r ∈ R∗+ , B(r) = .

µ0 j0 R si r ⩾ R



B(r)

µ0 j0 R

r
R

42
2.4.4 Plan avec courant surfacique
On munit R3 d’un repère cartésien (O, ux , uy , uz ).

On considère la densité surfacique de courant js uniforme de valeur js 0 ux portée par

la surface R2 × {0}.

La même étude des symétries et invariances que dans le cas volumique montre qu’il
existe B : R → R impaire telle que :

∀x, y, z ∈ R, B(x, y, z) = B(z)uy .

On montre de même que B est constante sur R2 × R∗+ et sur R2 × R∗− .

Soit L > 0. On considère la surface S = {0} × [0, L] × [−L, L], ce qui donne :
Z
B · tdℓ = (−B(L) + B(−L))2L.
∂S
Z
De plus, j · ndS = js 0 L, donc, d’après le théorème d’Ampère :
S

2(−B(L) + B(−L)) = µ0 js 0 .

L’imparité de B donne :
µ0 js 0
∀z ∈ R, B(z) = − sgn(z).
2

µ0 js 0
2 B(z)

µ0 js 0
−
2

Le champ magnétique n’est pas défini sur la plaque : c’est un défaut de modélisation
du fait du caractère surfacique de la densité de courant. Graphiquement, on retrouve
js
le même résultat en prenant ℓ → 0 et j0 ℓ → 0 .
2

43
2.4.5 Cylindre avec courant surfacique vertical
On munit R3 de coordonnées cylindriques (r, θ, z).

On considère la densité surfacique de courant js uniforme de valeur j0 s uz portée par

la surface D(0, R) × R.

La même étude des symétries et invariances que dans le cas volumique montre qu’il
existe B : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, B(r, θ, z) = B(r)uθ (θ).

On considère la surface S = D(O, r) × {0} avec r > 0, de sorte que :

Z
B · tdℓ = B(r)2πr.
∂S

Le théorème d’Ampère donne :

µ0 2πRj0 s si r ⩾ R



B(r)2πr = .

0 si r < R



Ainsi : 
R
µ0 j0 s r si r ⩾ R




∀r ∈ R∗+ , B(r) = .


0 si r < R



B(r)

µ0 j0 s

r
R

44
2.4.6 Cylindre avec courant surfacique ortho-radial
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité surfacique
de courant js = j0 s uθ (θ) portée par D(0, R) × R.

La même étude des symétries et invariances que dans le cas volumique montre qu’il
existe B : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, B(r, θ, z) = B(r)uz .

On considère S = {0} × [0, r] × [0, L] avec r, L > 0, de sorte que :

Z
B · tdℓ = (B(r) − B(0))L.
∂S

Des considérations de symétrie donnent B(0) = 0, ce qui fournit, avec le théorème

d’Ampère :
0 si r < R



∀r ∈ R∗ , B(r) =
+
.

µ0 j0 s si r ⩾ R



B(r)

µ0 j0 s

r
R

45
2.4.7 Fils fini et infini
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité de courant
j suivante :
j = I0 uz dC
avec I0 ∈ R et C = {0}2 × [−ℓ, ℓ] (où ℓ > 0). Autrement dit, on considère un fil
rectiligne de longueur 2ℓ parcouru par un courant uniforme.

Soit M ∈ R3 , de coordonnées cylindriques (r, θ, z). Le champ magnétique créé au

point M par un point P du fil est :
µ0 PM
dBP (M ) = I0 uz dℓ ∧ .
4π |PM|3

On note zP l’altitude de P , de sorte que PM = rur (θ) + (z − zp )uz .

Le champ B total créé au point M est :

µ0 I0 Z ℓ uθ (θ)
B(M ) = r 2 dzP .
4π −ℓ (r + (z − zP )2 )3/2

Or :
 
 z +ℓ z−ℓ 
Z ℓ
1 1
 
dzP = 2
 r r 
!2 − v !2  ,
v 
−ℓ (r + (z − zP ) )
2 2 3/2 r u
z+ℓ t1 + z − ℓ
u u 
1+
t u 
r r

donc :  
 z+ℓ z−ℓ 
µ0 I0 1 
 
B(M ) =  r −v r  uθ (θ).

4π r  u
v
z+ℓ
!2 !2
t1 + z − ℓ
u u 
1+
t u 
r r
Dans le cas où ℓ → +∞, on retrouve le champ magnétique créé par un fil infini
parcouru par une intensité uniforme I0 :
µ0 I0 1
B(M ) = uθ (θ).
2π r
On obtient le même résultat si on effectue le calcul avec le théorème d’Ampère.

En reprenant le cylindre de rayon R portant une densité volumique de courant

uniforme j0 , on retrouve, lorsque R → 0 et j0 πR2 → I0 , le fil infini parcouru par une
intensité uniforme I0 .

46
2.4.8 Spire de courant et solénoïdes fini et infini
On munit R3 de coordonnées cylindriques (r, θ, z). On considère la densité linéique
de courant I = I0 uθ (θ) portée par C(0, R) × {0}. Autrement dit, on considère une
spire de courant.

Calculons le champ magnétique créé par la distribution précédente sur l’axe (Oz).
On se donne un point M ∈ (Oz) d’ordonnée z. Le champ magnétique créé au point
M par un point P de la spire de coordonnée angulaire θ est :
µ0 PM
dBP (M ) = I0 uθ (θ)dℓ ∧ .
4π |PM|3

Or PM = −Rur (θ) + zuz , donc le champ B total créé au point M est :

µ0 I0 Z 2π Ruz + zur (θ)

B(M ) = Rdθ,
4π 0 (R2 + z 2 )3/2

donc :
µ0 I0 R2
B(M ) = uz .
2 (R2 + z 2 )3/2
On considère maintenant un empilement de spires de −ℓ à ℓ comprenant n spires par
unité de longueur, chacune parcourue par un courant d’intensité I0 uniforme. Un tel
dispositif est appelé solénoïde fini. Le principe de superposition fournit :

µ0 nI0 Z ℓ R2
B(M ) = dh.
2 −ℓ (R2 + (z − h)2 )3/2

Or :
#(ℓ−z)/R
du
"
Z ℓ
R2 Z (ℓ−z)/R
u
dh = = √
−ℓ (R + (z − h) )
2 2 3/2 (−ℓ−z)/R (1 + u )
2 3/2
1 + u2 (−ℓ−z)/R
 
 z+ℓ z−ℓ 
µ0 nI0
 
B(M ) =
 R −v R 
!2  uz .
v 
2 u
z+ℓ
!2
t1 + z − ℓ
u u 
1+
t u 
R R
Lorsque ℓ tend vers +∞, on trouve :

B(M ) = µ0 nI0 uz .

47
On considère toujours un solénoïde infini d’axe (Oz). La distribution de courant est
invariante par rotation d’axe (Oz) et par translation de vecteur uz , donc B l’est
aussi.

Soit M ∈ R3 . Le plan (M, ur (M ), uθ (M )) est un plan de symétrie de la distribution

de courant, donc (M, ur (M ), uθ (M )) est un plan d’antisymétrie de B, donc :

B(M ) ∈ Ruz .

Ainsi, il existe une fonction B : R∗+ → R telle que :

∀(r, θ, z) ∈ R∗+ × [0, 2π[×R, B(r, θ, z) = B(r)uz .

Le théorème d’Ampère applique à un cadre de côté (Oz) à l’intérieur du solénoïde

infini donne que B est constant à l’intérieur du solénoïde. De même, on montre que
B est constant à l’extérieur du solénoïde. En prenant un cadre traversant le solénoïde,
et en appliquant le théorème d’Ampère, on montre que :

si r ⩽ R

µ0 nI0


∀r ∈ R∗+ , B(r) = .

0 si r < R



B(r)

µ0 nI0

r
R

48
2.5 Potentiel vecteur
Au début de ce chapitre, nous avons calculé la circulation de B le long d’un contour
fermé, au lieu du flux de B à travers une surface fermée comme dans le cas du champ
électrique. Le but de ce paragraphe est de faire ce calcul pour le champ magnétique,
et d’en déduire un analogue du potentiel scalaire pour le champ magnétique.

Soit V un volume. On a :
!
Z
µ0 Z Z
PM
B · ndS = j(P ) ∧ dV · n(M )dS
∂V 4π M ∈∂V P ∈R 3
|PM|3
!
µ0 Z Z
PM
= 3 ∧ n(M )dS · j(P )dV.
4π P ∈R3 M ∈∂V |PM|

Une propriété de la moyenne donne :

Z
PM
∀P ∈ R3 , ∧ n(M )dS = 0,
M ∈∂V |PM|3

donc : Z
B · ndS = 0
∂V
Il s’agit de l’équation de Maxwell-Thomson. Par analogie avec le champ électrique,
on en déduit qu’il n’existe pas de monopole magnétique, autrement dit pas de charge
magnétique. On peut réécrire cette équation sous forme locale. Pour tout volume V :
Z
∇ · BdV = 0,
V

donc ∇ · B = 0. D’après le lemme de Poincaré, il existe A tel que :

B = ∇ ∧ A.

Un tel champ de vecteurs A est appelé potentiel vecteur. Il n’est pas unique.

En effet, si on prend χ une fonction numérique lisse sur R3 , on trouve que, si A est
un potentiel vecteur, alors A + ∇χ est encore un potentiel vecteur.

La forme explicite du premier paragraphe donne un potentiel vecteur A de B de la

forme :
µ0 Z 1
∀M ∈ R , A(M ) =
3
j(P ) dV.
4π P ∈R 3 PM

49
2.6 Relation de passage pour le champ magnétostatique
Dans le cas de la plaque infiniment fine parcourue par un courant uniforme, on a vu
que la composante tangentielle à la plaque du champ magnétique subissait un saut
de µ0 js , tandis que sa composante normale était continue à travers la plaque.

Il est possible de généraliser cette remarque à une distribution de courant surfacique

(S, js ).

Étant donné M ∈ S, on note B(M ± ) = B(M + 0± n(M )).

2.6.1 Continuité de la composante normale

On considère un point M ∈ S et un cylindre Σ de hauteur h et de rayon r, d’axe
dirigé par la normale n(M ) à S en M .

Écrivons le flux de B à travers Σ. On note Σ± les deux faces supérieure et inférieure

de Σ, et Σlat sa surface latérale. On a, en notant ν la normale à Σ :
Z
B · νdS = B(M ± ) (±n(M )) πr2 + O(r3 ) + O(h)
Σ± (r,h)→0

et : Z
B · νdS = O(h),
Σlat (r,h)→(0,0)

donc :
Z  
B · νdS = B(M + ) − B(M − ) · n(M )πr2 + O(r3 ) + O(h).
Σ (r,h)→(0,0)

On applique l’équation de Maxwell-Thomson, on fait tendre h vers 0, puis on

divise par r et on fait tendre r vers 0 :
 
B(M + ) − B(M − ) · n(M ) = 0.

On a donc continuité de la composante normale du champ magnétostatique à travers

S.

50
2.6.2 Discontinuité de la composante tangentielle
On considère un point M ∈ S, t un vecteur tangent à S en M , et un cadre
rectangulaire C dont deux cotés sont portés par t. Les deux cotés portés par t sont
de longueur r, et les deux autres cotés ont pour longueur h. On note C ± les deux
cotés portés par t, et Clat l’union des deux autres cotés. On a, en notant τ le vecteur
tangent à C :
Z
B · τ dℓ = ±B(M ± )r + O(r2 ) + O(h)
C± (r,h)→(0,0)

et : Z
B · τ dℓ = O(h),
Clat (r,h)→(0,0)

donc :
Z  
B · τ dℓ = B(M + ) − B(M − ) · tr + O(r2 ) + O(h).
C (r,h)→(0,0)

De plus, l’intensité passant par le cadre est :

js (M ) · tr + O(r2 ),

donc, d’après le théorème d’Ampère :

 
B(M + ) − B(M − ) · tr = µ0 js · tr + O(r2 ) + O(h).
(r,h)→(0,0)

On fait tendre h vers 0, puis on divise par r et on fait tendre r vers 0, ce qui donne :
 
B(M + ) − B(M − ) · t = µ0 js (M ) · t.

La relation précédente étant vraie pour tout vecteur tangent t à S en M :

 
B(M + ) − B(M − ) ∧ n(M ) = µ0 js (M ) ∧ n(M ).

Il s’agit de la relation de passage pour la composante tangentielle du champ magné-

tique.

51
2.6.3 Formulations faible et forte de la magnétostatique
Ce paragraphe nécessite des rudiments de théorie des distributions. Il peut être omis
en première lecture.

Soit V un volume de R3 et j ∈ D′ (V, R3 ). On appelle équations de la magnétostatique

le système : 
∇ · B = 0
∇ ∧ B = µ0 j

où l’inconnue est B ∈ D′ (V, R3 ).

Soit S une surface de V telle que V \ S est l’union disjointe de deux ouverts de R3
V± tels que le vecteur n normal à S soit dirigé vers V+ .

On suppose j de la forme j = j0 + js dS avec j0 localement intégrable sur V , dS la

mesure de surface de S, et js : S → R3 localement intégrable.
 
Soit B : V → R3 tel que B|V± admet un prolongement de classe B± ∈ C 1 V± , R3 .
On a équivalence entre :

le champ B vérifie, au sens faible sur V :


∇ · B =0
∇ ∧ B = µ0 j

et :

les champs B± vérifient, sur V± :


∇ · B =0
±
∇ ∧ B ± = µ0 j

et :
 
 B+
|S − B− |S · n = 0
  .
 B+
|S − B− |S ∧ n = µ0 js ∧ n

Ainsi, la relation de passage pour le champ magnétostatique est encodée dans les
équations de la magnétostatique sous forme faible.

52
2.7 Développement multipolaire du potentiel vecteur
Le but de ce paragraphe est d’étudier le potentiel vecteur créé à longue distance par
une densité volumique de courant j. Puisque j est conservé : ∇ · j = 0. On a :

µ0 +∞ 1
!
OM
∀M ∈ R \ {O}, A(M ) =
3
X
Mℓ
4π ℓ=0 |OM|ℓ+1 |OM|
avec : !
Z
OP
∀ℓ ∈ N, ∀r ∈ S , Mℓ (r) =2
|OP| Pℓ ℓ
· r j(P )dV
P ∈R3 |OP|
le ℓ-ième moment multipolaire. On peut calculer les deux premiers moments multi-
polaires. On rappelle que :
P0 = 1 et P1 = X,
donc :
Z
M0 (r) = j(P )dV
3
ZP ∈R
M1 (r) = (OP · r) j(P )dV
P ∈R3

On rappelle la formule du double produit vectoriel :

∀x, y, z ∈ R3 , (x ∧ y) ∧ z = (x · z)y − (y · z)x.

Ainsi :
(OP · r)j(P ) = (OP ∧ j(P )) ∧ r + (j(P ) · r)OP.
Or :
Z Z Z
j(P )j OPi = ∇ · (OPj j)OPi dV − OPi OPj ∇ · j(P )dV
P ∈R3 P ∈R3 P ∈R3
Z Z
=− OPj (j · ∇)(OPi )dV − OPi OPj ∇ · j(P )dV
P ∈R3 P ∈R3
Z Z
=− j(P )i OPj − OPi OPj ∇ · j(P )dV,
P ∈R3 P ∈R3

donc :
1Z 1Z
M1 (r) = (OP ∧ j(P )) ∧ rdV − (OP · r) (∇ · j) (P )OPdV.
2 P ∈R3 2 P ∈R3
1Z
On pose : µ = OP ∧ j(P )dV le moment dipolaire magnétique, de sorte que,
2 P ∈R3
puisque ∇ · j = 0 :
M1 (r) = µ ∧ r.

53
On a la description suivante :
Z
— Si j(P )dV ̸= 0, alors :
P ∈R3

µ0 Z
A(M ) ∼ j(P )dV.
M →∞ 4π P ∈R3
Z
— Si j(P )dV = 0 et µ ̸= 0, alors :
P ∈R3

µ0 µ ∧ OM
A(M ) ∼ .
M →∞ 4π |OM|3

Il s’agit par exemple du potentiel vecteur créé par une boucle de courant
parcourue par une intensité constante. Notons que, dans ce cas, µ = ISn où
I est l’intensité, S la surface de la boucle, et n le vecteur normal à la surface
d’orientation compatible avec I.
On peut développer les moments multipolaires en utilisant les harmoniques sphériques.
On rappelle l’identité suivante :

4π X ℓ
∀ℓ ∈ N, ∀x, y ∈ S , Pℓ (x · y) =
2
Yℓ,m (x)Yℓ,m (y)⋆ .
2ℓ + 1 m=−ℓ

On peut ainsi écrire sous une forme plus explicite la dépendance angulaire :

µ0 +∞ ℓ
1
!
OM
∀M ∈ R \ {O}, A(M ) =
3
X X
Mℓ,m Yℓ,m
4π ℓ=0 m=−ℓ |OM|ℓ+1 |OM|
avec :
!⋆
4π Z OP
∀ℓ ∈ N, ∀m ∈ J−ℓ, ℓK, ∀r ∈ S , Mℓ,m 2
= |OP|ℓ Yℓ,m j(P )dV
2ℓ + 1 P ∈R3 |OP|

les moments multipolaires sphériques.

54
2.8 Dipôle magnétostatique
On considère une distribution de courant j à support compact telle que ∇ · j = 0.
L’objectif de ce paragraphe est d’étudier le champ magnétique à grande Z distance
créé par j dans le cas où il n’y a pas de courant de fuite sortant de V : jdV = 0.
V
On choisit de plus une origine O.

2.8.1 Calcul du champ magnétique

D’après le paragraphe précédent, le potentiel vecteur s’écrit :

µ0 µ ∧ OM
A(M ) ∼ ,
M →∞ 4π |OM|3
donc le champ magnétique vérifie :

µ0 3(µ · OM)OM − |OM|2 µ

B(M ) ∼
M →∞ 4π |OM|5

TODO : a finir

55
2.9 Équations de la magnétostatique, existence et unicité
du champ magnétostatique
Ce paragraphe nécessite des connaissances mathématiques en théorie des distributions.
Il peut être omis en première lecture.

On se donne un volume V borné de R3 , j ∈ L2 (V, R3 ) et Bt ∈ L2 (∂V, R3 ) vérifiant

les relations de compatibilité :
Z Z
µ0 jdV + Bt dS = 0.
V ∂V

On cherche à résoudre le problème :

Trouver B ∈ H 1 (V, R3 ) tel que :

∇ · B =0





∇∧B = µ0 j .
B|∂V ∧ n = Bt



Le lemme de Poincaré donne :

{B ∈ H 1 (V, R3 ), ∇ · B = 0} = {∇ ∧ A, A ∈ H 2 (V, R3 ) et ∇ · A = 0}.

Une version modifiée du théorème d’existence et d’unicité de la solution à l’équation

de Poisson avec condition de Neumann donne l’existence et l’unicité de la solution
au problème suivant :

Trouver A ∈ H 2 (V, R3 ) tel que :

= µ0 j

−∆A




∇·A =0 ,
(∇ ∧ A)|∂V ∧ n = Bt



ce qui établit l’existence et l’unicité de la solution au problème initial sur B.

56
3 Induction
Dans cette section, on s’intéresse au mouvement des charges à l’intérieur d’un
circuit électrique filiforme, potentiellement déformable, en présence d’un champ
électromagnétique.

On se donne un circuit filiforme ∂S en présence d’un champ électromagnétique (E, B)

écrit dans le référentiel du laboratoire, supposé inertiel.

3.1 Champ électromoteur

On considère une charge ponctuelle (M, q) se déplaçant dans le circuit. La force de
Lorentz qui s’exerce sur elle s’écrit, dans le référentiel du laboratoire :

∀t ∈ R, F(t) = qE(M (t), t) + q (v(M (t)) + V(M (t))) ∧ B(M (t), t),

où V est la vitesse du circuit dans le référentiel du laboratoire, et v est la vitesse des

charges dans le référentiel du circuit.

On en déduit le travail élémentaire de la force de Lorentz, en remarquant que v et

t sont colinéaires :

∀t ∈ R, δW (t) = q (E(M (t), t) + V(M (t)) ∧ B(M (t), t)) · t(M (t))dℓ.

Autrement dit, formellement, on peut réécrire :

∀t ∈ R, δW (t) = qEm (M (t), t) · t(M (t))dℓ

avec ∀t ∈ R, ∀M ∈ ∂S(t), Em (M, t) = E(M, t) + V(M ) ∧ B(M, t) le champ électro-

moteur.

3.2 Loi de Lenz

On se donne un circuit ∂S (mobile ou immobile), placé dans un champ électroma-
gnétique (E, B). Une série d’expériences due à Lenz montre que la variation du flux
de B à travers S modifie l’intensité parcourant le circuit suivant la loi suivante, dite
de Lenz-Faraday :
dΦ
e=−
dt
Z Z
avec e = Em · dℓ (où Em est le champ électromoteur) et Φ = B · ndS est le
∂S S
flux du champ magnétique à travers S. Il s’agit d’une loi dite de modération.

La présence du signe − indique que le comportement du circuit en réaction à la

variation du flux du champ B va avoir tendance à minimiser le flux du champ
magnétique. Dans le pratique, si le circuit est plongé dans un champ magnétique
Bext extérieur de flux variable, un courant (dit induit) apparaît dans le circuit de
sorte à créer un champ magnétique induit Bind qui va s’opposer à la variation du
flux de Bext par rétroaction, de sorte que B = Bext + Bind .

57
3.3 Induction de Lorenz
Dans ce paragraphe, on suppose que le circuit est mobile ou déformable, mais que le
champ magnétique B est constant.

3.4 Induction de Neumann

3.5 Inductance

58
4 Équations de Maxwell
L’objectif de ce chapitre est de coupler les effets électriques et magnétiques étudiés
dans les chapitres précédents. On y établit et étudie les équations fondamentales de
l’électromagnétisme : les équations de Maxwell.

4.1 Conservation de la charge

On commence par étudier le lien entre la densité de charge ρ et la densité de courant
j. Le point clé est la conservation de la charge. Plus précisément, on se donne V un
volume quelconque, jouant le rôle d’un volume de particules chargées à un instant
initial t0 . On note V (t) le volume obtenu à l’instant t, et v le champ de vitesse des
particules chargées. La conservation de la charge s’écrit :

d Z
ρ(M, t)dV = 0.
dt V (t) t=t0

On note Φ le flot de v, de sorte que :

Z Z
ρ(M, t)dV = ρ(Φt (M ), t)|det(dΦt (M ))|dV,
M ∈V (t) M ∈V

donc :
d Z
ρ(M, t)dV
dt V (t) t=t0
Z
= (∂t ρ(M, t0 ) + v(M, t0 ) · ∇ρ(M, t0 ) + ρ(M, t0 )∇ · v(M, t)) dV
ZM ∈V
= (∂t ρ(M, t0 ) + ∇ · j(M, t0 )) dV
M ∈V

car j = ρv. La relation précédente étant valable pour tout volume V de R3 et pour
tout temps t0 , on en déduit :
∂t ρ + ∇ · j = 0.
Il s’agit de l’équation de conservation de la charge.

59
4.2 Formulation locale des équations de Maxwell
On postule que le théorème de Gauss reste valide en régime variable. Ainsi, pour
tout volume V de R3 : Z
Q
E · ndS =
∂V ε0
Z
où Q = ρdV est la charge contenue dans V . Un théorème de Stokes donne :
V
Z Z
E · ndS = ∇ · EdV,
∂V V

donc : Z 
ρ

∇·E − dV = 0.
V ε0
La relation précédente étant vraie pour tout volume V , on a :
ρ
∇·E = .
ε0
Il s’agit de l’équation de Maxwell-Gauss.

Passons à l’étude du théorème d’Ampère, qu’on suppose temporairement vrai. Pour

toute surface S de R3 : Z
B · tdℓ = µ0 I
∂S
Z
où I = j · ndS. Un théorème de Stokes donne :
S
Z Z
B · tdℓ = ∇ ∧ B · ndS,
∂S S

donc : Z
(∇ ∧ B − µ0 j) · ndS = 0,
S
donc, puisque S est quelconque :

∇ ∧ B − µ0 j = 0,

donc ∇ · j = 0, ce qui contredit en général l’équation de conservation de la charge.

Ainsi, le théorème d’Ampère est faux en régime variable. Afin de satisfaire la
conservation de la charge, on postule que l’équation suivante doit être vérifiée :

∇ ∧ B = µ0 (j + ε0 ∂t E) .

Il s’agit de l’équation de Maxwell-Ampère. Vérifions qu’elle satisfait la conservation

de la charge. En prenant la divergence, on trouve :

∇ · j + ε0 ∂t ∇ · E = 0,

donc, d’après l’équation de Maxwell-Gauss :

∂t ρ + ∇ · j = 0.

60
Ainsi, la charge est bien conservée.

On peut également réécrire l’équation de Maxwell-Ampère sous forme globale en

l’intégrant sur une surface S :

d Z
Z !
B · tdℓ = µ0 I + ε0 E · ndS
∂S dt S
Z
d Z
où I = ρdV . Le terme E · ndS est appelé courant de dérive. Notons que,
V dt S
dans le cas statique, on retrouve bien le théorème d’Ampère.

Quant à l’équation de Maxwell-Thomson, elle est postulée vraie en régime

variable : à ce jour, on n’a jamais observé de monopole magnétique.

Passons à l’équation de Lenz-Faraday : étant donnée une surface S (fixe), on a

dΦ
e=−
dt
Z Z
avec e = E · tdℓ et Φ = B · ndS. On applique un théorème de Stokes :
∂S S
Z
(∇ ∧ E + ∂t B) · ndS = 0,
S

donc, puisque S est quelconque :

∇ ∧ E = −∂t B.

Il s’agit de l’équation de Maxwell-Faraday.

61
4.3 Relations de passage en régime variable
On considère une surface S (fixe ou variable) portant des densités surfaciques de
charge et de courant σ, js . Le but de ce paragraphe est de montrer que les relations de
passage pour (E, B) sont encore valides en régime variable. Soit t ∈ R et M ∈ S(t).

4.3.1 Relation de passage pour le champ électrique

Le théorème de Gauss reste encore vrai en régime variable, donc :
  σ(M, t)
E(M + , t) − E(M − , t) · n(M, t) = .
ε0
Soit t un vecteur tangent à S(t) en M , et un cadre rectangulaire C dont deux cotés
sont portés par t. Les deux cotés portés par t sont de longueur r, et les deux autres
cotés ont pour longueur h. On a, en notant τ le vecteur tangent à C :
Z  
E · τ dℓ = E(M + , t) − E(M − , t) · tr + O(r2 ) + O(h),
C (r,h)→(0,0)

donc, d’après l’équation de Maxwell-Faraday, en notant Σ le rectangle dont le

bord est C :
  Z
E(M + , t) − E(M − , t) · tr + O(r2 ) + O(h) = − ∂t B( · , t) · νdS
(r,h)→(0,0) Σ

où ν est la normale à Σ. Or :
Z
∂t B( · , t) · νdS = O(h),
Σ (r,t)→(0,0)

donc, en faisant tendre h vers 0, puis en divisant par r et en faisant tendre r vers 0 :
 
E(M + , t) − E(M − , t) · t = 0.

Or t est un vecteur tangent quelconque à S(t) en M , ce qui donne :

 
E(M + , t) − E(M − , t) ∧ n(M, t) = 0.

62
4.3.2 Relation de passage pour le champ magnétique
La loi de Maxwell-Thomson reste encore vraie en régime variable, donc :
 
B(M + , t) − B(M − , t) · n(M, t) = 0.

On note t un vecteur tangent à S(t) en M , et un cadre rectangulaire C dont deux

cotés sont portés par t. Les deux cotés portés par t sont de longueur r, et les deux
autres cotés ont pour longueur h. On a :
Z  
B · τ dℓ = B(M + ) − B(M − ) · tr + O(r2 ) + O(h),
C (r,h)→(0,0)

donc, d’après l’équation de Maxwell-Ampère, en notant Σ le rectangle dont le

bord est C :
 
B(M + , t) − B(M − , t) · tr
Z
= µ0 js (M, t) · tr + µ0 ε0 ∂t E( · , t) · ν + O(r2 ) + O(h)
(r,h)→(0,0) Σ

où ν est la normale à Σ. Or :
Z
∂t E( · , t) · ν = O(h),
Σ (r,h)→(0,0)

donc, en faisant tendre h vers 0, puis en divisant par r et en faisant tendre r vers 0 :
 
B(M + , t) − B(M − , t) · t = µ0 js (M, t) · t.

La relation précédente étant vraie pour tout vecteur tangent t à S(t) en M :

 
B(M + , t) − B(M − , t) ∧ n(M, t) = µ0 js (M, t) ∧ n(M, t).

Il s’agit de la relation de passage pour la composante tangentielle du champ magné-

tique.

63
4.4 Des équations de Maxwell à l’équation des ondes
Le but de ce paragraphe est de mettre en évidence la nature ondulatoire du champ
électromagnétique.

On commence par calculer le laplacien vectoriel du champ électrique :

∆E = ∇(∇ · E) − ∇ ∧ (∇ ∧ E).

On applique ensuite les équations de Maxwell-Faraday et de Maxwell-Gauss :

1
∆E = ∇ρ + ∂t ∇ ∧ B.
ε0
L’équation de Maxwell-Ampère donne :
1
∆E = ∇ρ + ∂t (µ0 j + µ0 ε0 ∂t E),
ε0
1 2 1
ce qui se réécrit, en posant □ = 2
∂t − ∆ avec c = √ :
c ε0 µ 0

1
□E = −µ0 ∂t j − ∇ρ.
ε0
Il s’agit clairement d’une équation des ondes de célérité c. Montrons que le champ
magnétique vérifie aussi une équation des ondes de célérité c. On a :

∆B = ∇(∇ · B) − ∇ ∧ (∇ ∧ B).

Les équations de Maxwell-Ampère et de Maxwell-Thomson donnent :

∆B = −∇ ∧ (µ0 j + µ0 ε0 ∂t E) = −µ0 ∇ ∧ j − µ0 ε0 ∂t ∇ ∧ E

On applique ensuite l’équation de Faraday :

∆B = −µ0 ∇ ∧ j + µ0 ε0 ∂t2 B.

Ainsi :
□B = µ0 ∇ ∧ j.
Ainsi, le champ (E, B) est une onde de célérité c.

64
4.5 Potentiels
Le but de ce paragraphe est de réécrire les équations de Maxwell en termes de
potentiels, et de voir leur application à la résolution explicite des équations de
Maxwell.

On part de l’équation de Maxwell-Thomson : ∇ · B = 0, donc il existe un champ

de vecteurs A tel que : B = ∇ ∧ A. On réinjecte dans l’équation de Maxwell-
Faraday :
∇ ∧ (E + ∂t A) = 0,
donc il existe une fonction scalaire Φ telle que :

E + ∂t A = −∇Φ.

Un tel champ A est appelé potentiel vecteur, et un tel champ Φ est appelé potentiel
scalaire. Ils vérifient par définition les relations :

B = ∇ ∧ A et E = −∇Φ − ∂t A.

Une telle paire (Φ, A) n’est pas unique. En effet, en partant d’un potentiel vecteur,
on remarque que, pour toute fonction scalaire χ, si (Φ, A) est un potentiel pour
(E, B), alors (Φ − ∂t χ, A + ∇χ) est encore un potentiel pour (E, B).

La transformation (Φ, A) → (Φ − ∂t χ, A + ∇χ) est appelée transformation de jauge.

Toutes les quantités physiques calculées à partir des potentiels doivent être laissées
invariantes lorsqu’on applique une transformation de jauge.

On peut réécrire les équations de Maxwell avec sources en termes de potentiels :

ρ
— équation de Maxwell-Gauss : −∆Φ − ∂t ∇ · A = .
ε0
— équation de Maxwell-Ampère : ∇ ∧ (∇ ∧ A) = µ0 j − µ0 ε0 ∂t (∇Φ + ∂t A).
Voyons comment écrire l’équation des ondes à l’aide des potentiels. La deuxième
équation peut se réécrire :

∇ ∧ (∇ ∧ A) = ∇(∇ · A) − ∆A

□A = µ0 j − ∇ (µ0 ε0 ∂t Φ + ∇ · A) ,
ce qui invite à choisir (Φ, A) de sorte que :

µ0 ε0 ∂t Φ + ∇ · A = 0.

Une telle contrainte sur le choix de (Φ, A) est appelée jauge d’onde. Dans cette jauge,
les équations de Maxwell avec sources se réécrivent :
ρ
□Φ = et □A = µ0 j.
ε0
Il s’agit de deux équations des ondes à la célérité c. Montrons qu’on peut effectivement
faire un tel choix de jauge. On se donne (Φ0 , A0 ) un potentiel quelconque, et χ un

65
champ scalaire. Sous la transformation de jauge (Φ0 , A0 ) → (Φ0 − ∂t χ, A0 + ∇χ),
on a :
µ0 ε0 ∂t Φ0 + ∇ · A0 → µ0 ε0 ∂t Φ0 + ∇ · A0 − □χ.
On choisir donc χ tel que :

□χ = µ0 ε0 ∂t Φ0 + ∇ · A0 ,

de sorte que le potentiel (Φ, A) = (Φ0 − ∂t χ, A0 + ∇χ) vérifie la jauge d’onde. Notons
que la jauge d’onde est invariante sous les transformations de jauge associées aux
champs scalaires de d’alembertien nul.

On peut trouver explicitement des solutions aux équations de Maxwell. Puisque

(E, B) est une onde de célérité c, une approche possible consiste à prendre les
potentiels de l’électrostatique et de la magnétostatique, puis de les écrire sous forme
retardée :
!
|PM|
ρ P, t −
1 Z c
∀M ∈ R , ∀t ∈ R, Φ(M, t) =
3
dV
4πε0 P ∈R3 |PM|
et : !
|PM|
j P, t −
µ0 Z c
∀M ∈ R , ∀t ∈ R, A(M, t) =
3
dV.
4π P ∈R3 |PM|
On vérifie facilement que (Φ, A) satisfont les équations de Maxwell avec sources en
termes de potentiels, ce qui donne une expression de (E, B) satisfaisant les équations
de Maxwell :
! !
|PM| |PM|
 
 ρ P, t − ∂t ρ P, t −
1 c c
Z  
E(M, t) = +  PM dV
  

4πε0 P ∈R 
3

|PM| 3 c|PM| 2 




!
|PM|
∂t j P, t −
1 Z c
− dV
4πε0 P ∈R3 c |PM|
2

et :
! !
|PM| |PM|

 j P, t − ∂t j P, t −
µ0 c c
Z 
B(M, t) = +
 
  ∧ PMdV
4π P ∈R 
3

|PM| 3 c|PM| 2 


Il s’agit des équations de Jefimenko.

On remarque qu’on retrouve les expressions déduites de lois de Coulomb, Biot et

Savart dans les cas statiques.

66
De plus, les invariances et symétries de B en fonction de celles de j sont les mêmes que
dans le cas statique. Cependant, le terme de courant dans E fait que les invariances
et symétries de E sont modifiées :
— si ρ, j sont invariantes sous une transformation R de déterminant positif, alors
E est encore invariant sous R. Cette conclusion est fausse en général dans le
cas où R est négative.
— si Π (respectivement Π⋆ ) est un plan de symétrie (respectivement d’antisymétrie)
de ρ et j, alors Π (respectivement Π⋆ ) est un plan de symétrie (respectivement
d’antisymétrie) de E.

67
4.6 Modes de Fourier
Les équations de Maxwell peuvent se réécrire à l’aide de la transformée de Fourier,
en choisissant une origine spatiale O ∈ R3 . On note k, ω les variables de Fourier
duales associées a OM et t En prenant les transformées de Fourier en temps et en
espace, on a :

ρ(k, ω)
∀R3 , ∀ω ∈ R, −ik · E(k, ω) =
ε0
−ik ∧ B(k, ω) = µ0 j(k, ω) − iµ0 ε0 ωE(k, ω)
k · B(k, ω) = 0
k ∧ E(k, ω) = ωB(k, ω)

Commentons les expressions précédentes.

— La troisième équation (celle de Maxwell-Thomson) donne que B est trans-
verse.
— La première équation (celle de Maxwell-Gauss) donne que E est transverse
si et seulement si ρ = 0.
— La quatrième équation (celle de Maxwell-Faraday) donne la relation de
structure des ondes électromagnétiques :
k
∀k ∈ R3 , ∀ω ∈ R, B(k, ω) = ∧ E(k, ω).
ω

Nous étudierons plus en détails l’équation de Maxwell-Ampère dans une future

section, consacrée a la loi d’Ohm locale.

On peut aussi réécrire l’équation des ondes pour le champ électromagnétique dans le
vide : !
|k|2  
ω − 2 E(k, ω) = 0 et ω 2 − c2 |k|2 B(k, ω) = 0.
2
c
Ainsi, le support du champ électromagnétique en variables (k, ω) est inclus dans le
cône d’équation ω = c|k| : cette dernière relation est la relation de dispersion du
vide.

4.7 ARQS

68
4.8 Aspects énergétiques du champ électromagnétique
On cherche à donner du sens à la notion d’énergie du champ électromagnétique,
et d’établir une équation de conservation de cette dernière. On peut chercher à la
décomposer en deux parties : une électrique, et une magnétique, chacune proportion-
nelle à |E|2 et |B|2 . Par analyse dimensionnelle, on définit ainsi la densité d’énergie
du champ électromagnétique (E, B) par :
1 1
u = ε0 |E|2 + |B|2 .
2 2µ0
On peut calculer la dérivée temporelle de u :
1
∂t u = ε0 E · ∂t E + B · ∂t B.
µ0

Les équations de Maxwell-Ampère et de Maxwell-Faraday donnent :

1 1 1
∂t u = E · (∇ ∧ B − µ0 j) − B · ∇ ∧ E = (E · ∇ ∧ B − B · ∇ ∧ E) − j · E,
µ0 µ0 µ0
1
donc, en posant Π = E ∧ B le vecteur de Poynting, on a :
µ0
∂t u + ∇ · Π = −j · E.

Il s’agit de l’équation de Poynting, écrite sous forme locale. On peut l’écrire sous
forme globale en l’intégrant sur un volume V de R3 :
dU Z Z
+ Π · ndS = − j · EdV.
dt ∂V V
Z
avec U = udV . Ainsi, Π peut s’interpréter comme le vecteur courant d’énergie
V
électromagnétique. Une partie de cette dernière est dissipée via le terme de droite,
qui n’est autre qu’une manifestation de l’effet Joule.

69
4.9 Formulation des équations de Maxwell comme un pro-
blème d’évolution
Cette section nécessite des rudiments de théorie des distributions, et peut être
omis en première lecture. On cherche à donner un sens mathématique plus précis
aux équations de Maxwell. Une bonne formulation consiste à les voir comme un
problème d’évolution, et à remplacer les dérivées temporelles en utilisant la formule
des sauts. Soit V un volume de R3 . On se donne ρ0 ∈ S ′ (R3 , R), E0 ∈ S ′ (R3 , R3 )
vérifiant les équations de l’électrostatique :

ρ0
∇ · E0
 =
ε0
∇ ∧ E0 =0


et j0 ∈ S ′ (R3 , R3 ), B0 ∈ S ′ (R3 , R3 ) vérifiant les équations de la magnétostatique :


∇ · B =0
0
.
∇ ∧ B 0 = µ0 j 0

On dit que des distributions (ρ, j) ∈ S ′ (R3 × R, R) × S ′ (R3 × R, R3 ) et (E, B) ∈

S ′ (R3 × R, R3 ) × S ′ (R3 × R, R3 ) à supports dans R3 × R+ (elles sont dite à support
dans le futur) sont solutions aux équations de Maxwell si :
 ρ
∇ · E
 =


 ε0
∇·B = 0


.
∇ ∧ E = −∂ B + B ⊗ δ

t 0 0




∇ ∧ B = µ0 (j + ε0 (∂t E − E0 ⊗ δ0 ))



En prenant la divergence de l’équation de Maxwell-Faraday et en utilisant

l’équation de Maxwell-Gauss, on trouve l’équation de conservation de la charge :

∂t ρ + ∇ · j = ρ0 ⊗ δ0 .

Le même calcul que précédemment donne les équations des ondes suivantes :
∇ρ

□E =− − µ0 (∂t j − j0 ⊗ δ0 ) + ε0 µ0 E0 ⊗ δ0′

ε0
□B = µ ∇ ∧ j + ε µ B ⊗ δ ′

0 0 0 0 0

L’équation de Maxwell-Thomson donne l’existence de A ∈ S ′ (R3 × R, R3 ) telle

que B = ∇ ∧ A. De plus, il existe A0 ∈ S ′ (R3 , R3 ) tel que B0 = ∇ ∧ A0 . Ainsi :

∇ ∧ (E + ∂t A − A0 ⊗ δ0 ) = 0,

donc il existe Φ ∈ S ′ (R3 × R, R) tel que :

E = −∇Φ − ∂t A + A0 ⊗ δ0 .

70
En transformée de Fourier, les équations de Maxwell se réécrivent comme suit :


 ik · B(k, ω) =0
1



ik ∧ E(k, ω) = −iωB(k, ω) + √ B0 (k)


2π .
1

 !!
= µ0 j(k, ω) + ε0 iωE(k, ω) − √ E0 (k)

ik ∧ B(k, ω)


2π


TODO

71
5 Conducteurs
5.1 Modèle de Drude
Le but de ce paragraphe est de présenter un modèle simplifié donnant une relation
entre j et E. On suppose que B = 0 dans le référentiel du laboratoire, et que l’espace
ne contient qu’une seule charge ponctuelle en mouvement (M, q, m). Le principe
fondamental de la dynamique donne, les champs étant écrits au point M (t) :

dv
m (t) = qE(M (t), t) − αv(t)
dt
avec α > 0 modélisant l’interaction entre M et le milieu dans lequel elle évolue. On
peut réécrire l’équation précédente en termes de courant :
m dj α
(t) = qE(t) − j(t),
q dt q
donc :
j(ω) = γ(ω)E(ω)
γ0 q2 m
avec γ(ω) = où γ0 = et τ0 = . La relation précédente peut se réécrire
1 + iτ0 ω α α
en notation réelle :
γ(t)
j(t) = ⋆ E(t).
(2π)1/2
Notons que, si le spectre de E est concentré en basse fréquence, on a :

j(t) ≃ γ0 E(t).

72
5.2 Loi d’Ohm locale et effet de peau
Plus généralement, la loi d’Ohm locale postule que, pour une large classe de conduc-
teurs, pour les champs électriques faibles, il existe une fonction γ : R3 × R → C
telle que :
∀M ∈ R3 , ∀ω ∈ R, j(M, ω) = γ(M, ω)E(M, ω).
Les conducteurs vérifiant la loi précédente sont dits ohmiques. Dans la suite du
paragraphe, on se place dans un conducteur ohmique.

On suppose que γ ne dépend pas de la position. On rappelle l’équation de conservation

de la charge :
∂t ρ + ∇ · j = ρ0 ⊗ δ0 ,
donc, en prenant la transformée de Fourier en temps et appliquant la loi d’Ohm :
1
iωρ(M, ω) + γ(ω)∇ · E(M, ω) = ρ0 (M ) √
2π
Le théorème de Gauss fournit :
γ(ω) 1
iωρ(M, ω) + ρ(M, ω) = ρ0 (M ) √ ,
ε0 2π
donc :
1 ε0 ρ0 (M )
ρ(M, ω) = √ .
2π γ(ω) + iε0 ω
Regardons ce qui se passe dans le cas particulier d’une conductivité de Drude : il
γ0
existe γ0 , τ0 ∈ R∗+ tels que ∀ω ∈ R, γ(ω) = .
1 + iτ0 ω
s
γ0 γ0
Ainsi, en posant ω0 = et Q = :
ε0 τ 0 ε0 ω0

i ω
1+
ε0 Q ω0
∀M ∈ R3 , ∀ω ∈ R, ρ(M, ω) = √ 2 ρ0 (M ).
2πγ0 i ω ω

1+ −
Q ω0 ω0
Notons que le milieu agit comme un filtre passe-bande. On peut calculer la densité
en temps :
i ω
1+
ε0 ρ0 (M ) Z
Q ω0
∀M ∈ R3 , ∀t ∈ R, ρ(M, t) =  2 e dω
iωt
2πγ0 R i ω ω
1+ −
Q ω0 ω0
i
1+ x
ε0 ω0 ρ0 (M ) Z
Q
= eixω0 t dx.
2πγ0 R i
1+ x−x 2
Q

73
Or :
i
1+ x
Q
Z
eixω0 t dx
R i
1 + x − x2
Q
1
 !
2π −1


2Q2 1
 ! !!

t t
exp − − exp − 1R+ (t) si Q <



1 2
 s
τ+ τ−


4

−





 Q 2
1


=  2πe−ω0 t (2 − ω0 t)1R+ (t) si Q =
2
1

  

1−


1 1

t  2Q2

   
2π exp  cos (ω1 t) − s sin (ω  1 + (t) si Q >
 
− t)


1 R
1 2


 τ Q 
1−

 
4Q2



!−1
1 2Q 1
s s
ω0
où τ± = ∓ ω0 −1 ,τ= et ω1 = 1− ω0 . Ainsi :
2Q 4Q2 ω0 4Q2

∀M ∈ R3 , ∀t ∈ R, ρ(M, t)
1
 !
−1


ε0 ω0 ρ0 (M ) 2Q2 1
 ! !!
t t


exp − − exp − 1R+ (t) si Q <



1 2
 s

γ0 τ+ τ−
4

−



Q2



ε0 ω0 ρ0 (M ) −ω0 t 1



= e (2 − ω0 t)1R+ (t) si Q =

 γ0 2
1

  

1−


ε0 ω0 ρ0 (M ) 1 2Q2 1

t 

   
exp −  cos (ω1 t) − sin (ω1 t)
 1R+ (t) si Q >

 

2

1
s
γ0 τ Q


1−

 

4Q

 2

On retrouve les trois régimes de l’oscillateur harmonique amorti de pulsation propre

ω0 et de facteur de qualité Q.

Dans le cas où τ0 = 0 (ie la conductivité est réelle, ie le courant j est en phase avec
le champ E), en reprenant le calcul précédent on a :
γ0
 
∀M ∈ R , ∀t ∈ R , ρ(M, t) = ρ0 (M ) exp − t 1R+ (t).
3 +
ε0
ε0
Le temps est appelé temps de relaxation diélectrique. Il donne l’ordre de grandeur
γ0
que met la densité volumique de charge à s’annuler dans le conducteur.

74
Dans la suite de ce paragraphe, on reprend une conductivité γ quelconque et on
poursuit notre étude autour du champ électrique. D’après les équations de Maxwell-
Ampère et de Maxwell-Faraday :
iωµ0 ω2
−|k|2 E(k, ω) = γ(k, ω) ⋆k E(k, ω) − E(k, ω),
(2π)3/2 c2
Il est difficile de tirer de cette expression une relation de dispersion. Si γ est indépen-
dante de k, alors :
ω2
|k|2 = µ0 ℑ(γ(ω)) + et ℜ(γ(ω)) = 0.
c2
Cette relation contraint beaucoup γ, qui doit être imaginaire pur.

Étudions ce qui se passe pour une autre géométrie où on ne peut pas utiliser la
transformée de Fourier en espace : le cas du demi-espace R2 × R+ . On suppose
que le champ électrique est de la forme E(M, ω) = f (z)E0 (ω) et que γ ne dépend
que de ω. On a :
ω2
f ′′ (z)E0 (ω) = iωγ(ω)E0 (ω)f (z) − E0 (ω)f (z).
c2
Soit ω ∈ R tel que E0 (ω) ̸= 0. On a :
!
ω2
∀z ∈ R , f (z) = iωγ(ω) − 2 f (z).
+ ′′
c

ω2 ω2
Après calcul, les racines carrées de iωγ − = −ωℑ(γ(ω)) − + iωℜ(γ(ω)) sont :
c2 c2
i
k+
δ
avec :
vv
uu !2
ω2 ω2
uu
ωℑ(γ(ω)) + 2 + (ωℜ(γ(ω)))2 − ωℑ(γ(ω)) −
ut
u
u c c2
k(ω) = ±
t
2
et :
vv
uu !2
ω2 ω2
uu
ωℑ(γ(ω)) + 2 + (ωℜ(γ(ω)))2 + ωℑ(γ(ω)) +
ut
u
1 c c2
u
= ± sgn (ωℜ(γ(ω)))
t
δ(ω) 2
(avec, par convention, δ = +∞ si ωℜ(γ(ω)) = 0).

Ainsi, il existe f± (ω) ∈ C tels que :

∀z ∈ R+ , f (z) = f+ (ω) exp(i(k(ω) + i/δ(ω))z) + f− (ω) exp(−i(k(ω) + i/δ(ω))z).

75
On choisit le signe tel que δ > 0, et on interdit les solutions divergentes dans le
conducteur), ce qui donne f− (ω) = 0. Ainsi, quitte à absorber le facteur f+ (ω) dans
E0 (ω), on a : !
z
E(M, ω) = exp − exp (ik(ω)z) E0 (ω).
δ(ω)
Le champ E pénètre dans le conducteur sur une épaisseur caractéristique δ(ω),
appelée épaisseur de peau.

TODO : tracer les graphes de δ et k pour un profil de Drude.

5.3 Conducteur parfait

Dans ce paragraphe, on se place dans le cadre où γ = ∞. Dans ce cas, on a E = 0
et, pour le cas du demi-plan : δ = 0.

5.4 Plasmas
Parler de la conductivité, parler du filtrage

5.5 Champs moyennes

A REPRENDRE On étudie le champ electromagnetique a une echelle de volume δV .
Plus precisement, on convole toutes les quantites physiques par une certaine fonction
densité d’ecart-type δV . On peut prendre par exemple la gaussienne GδV . Le choix
de la densité n’a pas d’importance physique dans la limite où δV tend vers 0. Plus
precisement, on remplace toute quantite physique g par sa moyenne suivant GδV :
Z
⟨g⟩(M, t) = g(P, t)GδV (MP)dV,
P ∈R3

de sorte que ∂i G = ∂i ⟨G⟩ et ∂t G = ∂t ⟨G⟩. On peut reecrire les equations de Maxwell

pour les grandeurs moyennees :

⟨ρ⟩
∇ · ⟨E⟩ =
ε0
∇ ∧ ⟨B⟩ = µ0 (⟨j⟩ + ε0 ∂t ⟨E⟩)
∇ · ⟨B⟩ = 0
∇ ∧ ⟨E⟩ = −∂t ⟨B⟩.

On peut chercher a comparer la difference ⟨ρ⟩ − ρ :

dGδV
⟨ρ⟩ − ρ = ρ ⋆ GδV − ρ = ρ ⋆ δV + o(δV )
dδV δV =0

La distribution ⟨ρ⟩ − ρ s’apparente a un moment dipolaire, ce qui invite a poser la

dGδV
densité volumique dipolaire P = −ρ ⋆ .
dδV δV =0

76
On peut faire de même pour le courant. On a :

dGδV
⟨j⟩ − j = j ⋆ δV + o(δV ).
dδV δV =0

Or :
dGδV Z
dGδV
j⋆ = j(P ) (PM)dV =
dδV δV =0 P ∈R3 dδV δV =0

1
!
x2
GδV (x) = exp −
(2π)3/2 δV 1/3 2δV 2/3
!
−1 x2
G′δV (x) = δV −4/3
exp −
3(2π)3/2 2δV 2/3

5.6 Milieux diélectriques

On se place dans l’approximation des regimes quasi-stationnaire électrique. On se
donne un milieu qui contient une densité volumique de dipoles p. Autrement dit, il
existe une echelle de volume On suppose que le milieu possede une densité

5.7 Milieux diamagnétiques

6 Reflexion/refraction/transmission
On considere la situation ou une onde Ei arrive a l’interface entre deux conducteurs

6.1 Vide
6.2 Plasmas
6.3 Conducteurs
6.4 Diélectriques

7 Vers l’optique ondulatoire

7.1 Equation eikonale
Faire Snell-Descartes, angle de Brewster et polarisation

77
8 Appendice
8.1 Analyse vectorielle
8.1.1 Opérateurs différentiels
On note Ω un ouvert de R3 .
Définition 1 (Gradient). Soit f ∈ C 1 (Ω, R). Le gradient de f est le champ vectoriel
∇f ∈ C 0 (Ω, R3 ) tel que :
∀M ∈ Ω, ∀h ∈ R3 , df (M ) · h = ∇f (M ) · h.
Proposition 2. L’application C 1 (Ω, R) → C 0 (Ω, R3 ), f 7→ ∇f est linéaire, et vérifie
la règle de Leibniz :
∀f, g ∈ C 1 (Ω, R), ∇(f g) = f ∇g + (∇f )g.
Proposition 3 (Gradient en coordonnées cartésiennes). Soit (ux , uy , uz ) une base
orthonormée de R3 et f ∈ C 1 (Ω, R). On a :
∇f = (∂x f )ux + (∂y f )uy + (∂z f )uz .
Proposition 4 (Gradient en coordonnées cylindriques). Soit (r, θ, z) des coordonnées
cylindriques et f ∈ C 1 (Ω, R). On a, partout où les coordonnes cylindriques sont bien
définies :
1
∇f = (∂r f )ur + (∂θ f )uθ + (∂z f )uz .
r
Proposition 5 (Gradient en coordonnées sphériques). Soit (r, θ, φ) des coordonnées
sphériques et f ∈ C 1 (Ω, R). On a, partout où les coordonnées sphériques sont bien
définies :
1 1
∇f = (∂r f )ur + (∂θ f )uθ + (∂φ f )uφ .
r r sin(θ)
Définition 6 (Divergence). Soit V ∈ C 1 (Ω, R3 ). On définit la divergence de V
∇ · V ∈ C 0 (Ω, R) par :
∀M ∈ Ω, ∇ · V(M ) = Tr(dV(M )).
Proposition 7. L’application C 1 (Ω, R3 ) → C 0 (Ω, R), V 7→ ∇ · V est linéaire, et
vérifie la règle de Leibniz :
∀V ∈ C 1 (Ω, R3 ), ∀f ∈ C 1 (Ω, R), ∇ · (f V) = (∇f ) · V + f ∇ · V.
Proposition 8 (Divergence en coordonnées cartésiennes). Soit (ux , uy , uz ) une base
orthonormée de R3 et V ∈ C 1 (Ω, R3 ). On a :
∇ · V = ∂x Vx + ∂y Vy + ∂z Vz .
Proposition 9 (Divergence en coordonnées cylindriques). Soit (r, θ, z) des coordon-
nées cylindriques et V ∈ C 1 (Ω, R3 ). On a, partout où les coordonnées cylindriques
sont bien définies :
1 1
∇ · V = ∂r (rVr ) + ∂θ Vθ + ∂z Vz .
r r

78
Proposition 10 (Divergence en coordonnées sphériques). Soit (r, θ, φ) des coordon-
nées sphériques et V ∈ C 1 (Ω, R3 ). On a, partout où les coordonnées sphériques sont
bien définies :
1  2  1 1
∇·V = ∂ r r Vr + ∂θ (sin(θ)Vθ ) + ∂φ Vφ
r 2 r sin(θ) r sin(θ)

Définition 11 (Rotationnel). Soit V ∈ C 1 (Ω, R3 ). Le rotationnel de V est le champ

vectoriel ∇ ∧ V ∈ C 0 (Ω, R3 ) tel que :

∀X ∈ R3 , ∀M ∈ Ω, ∇ · (X ∧ V) = −X · (∇ ∧ V(M )) .

Proposition 12. L’application C 1 (Ω, R3 ) → C 0 (Ω, R3 ), V 7→ ∇ ∧ V est linéaire et

vérifie la règle de Leibniz :

∀V ∈ C 1 (Ω, R3 ), ∀f ∈ C 1 (Ω, R), ∇ ∧ (f V) = f ∇ ∧ V + (∇f ) ∧ V.

Proposition 13 (Rotationnel en coordonnées cartésiennes). Soit (ux , uy , uz ) une

base orthonormée directe de R3 et V ∈ C 1 (Ω, R3 ). On a :

∇ ∧ V = (∂y Vz − ∂z Vy ) ux + (∂z Vx − ∂x Vz ) uy + (∂x Vy − ∂y Vx ) uz .

Proposition 14 (Rotationnel en coordonnées cylindriques). Soit (r, θ, z) des coor-

données cylindriques et V ∈ C 1 (Ω, R3 ). On a, partout où les coordonnées cylindriques
sont bien définies :
1 1
 
∇∧V= ∂θ Vz − ∂z Vθ ur + (∂z Vr − ∂r Vz ) uθ + (∂r (rVθ ) − ∂θ Vr ) uz .
r r
Proposition 15 (Rotationnel en coordonnées sphériques). Soit (r, θ, φ) des coor-
données sphériques et V ∈ C 1 (Ω, R3 ). On a, partout où les coordonnées sphériques
sont bien définies :
1
∇∧V= (∂θ (sin(θ)Vφ ) − ∂φ Vθ ) ur
r sin(θ)
1 1
!
+ ∂φ Vr − ∂r (rVφ ) uθ
r sin(θ) r
1
+ (∂r (rVθ ) − ∂θ Vr ) uφ
r
Définition 16 (Laplacien scalaire). Soit f ∈ C 2 (Ω, R). Le laplacien de f est la
fonction scalaire ∆f ∈ C 0 (Ω, R) définie par :

∀M ∈ Ω, ∆f (M ) = ∇ · (∇f ).

Proposition 17. L’application C 2 (Ω, R) → C 0 (Ω, R), f 7→ ∆f est linéaire.

Proposition 18 (Laplacien en coordonnées cartésiennes). Soit (ux , uy , uz ) une base
orthonormée de R3 . On a :

∆f = ∂x2 f + ∂y2 f + ∂z2 f.

79
Proposition 19 (Laplacien scalaire en coordonnées cylindriques). Soit (r, θ, z)
des coordonnées cylindriques et f ∈ C 2 (Ω, R). On a, partout où les coordonnées
cylindriques sont bien définies :
1 1
∆f = ∂r (r∂r f ) + 2 ∂θ2 f + ∂z2 f.
r r
Proposition 20 (Laplacien scalaire en coordonnées sphériques). Soit (r, θ, φ) des
coordonnées sphériques et f ∈ C 2 (Ω, R). On a, partout où les coordonnées sphériques
sont bien définies :
1 1 1
∆f = ∂r (r2 ∂r f ) + 2 ∂θ (sin(θ)∂θ f ) + 2 2 ∂φ2 f.
r 2 r sin(θ) r sin (θ)

Définition 21 (Laplacien vectoriel). Soit V ∈ C 2 (Ω, R3 ). On définit le laplacien

vectoriel de V ∆ ∈ C 0 (Ω, R3 ) par :

∆V = ∇(∇ · V) − ∇ ∧ (∇ ∧ V).

Proposition 22 (Laplacien vectoriel en coordonnées cartésiennes). Soit (ux , uy , uz )

une base orthonormée directe de R3 et V ∈ C 2 (Ω, R3 ). On a :

∆V = (∆Vx )ux + (∆Vy )uy + (∆Vz )uz .

Proposition 23. On a :

∀f ∈ C 2 (Ω, R), ∇ ∧ (∇f ) = 0

et :
∀V ∈ C 2 (Ω, R3 ), ∇ · (∇ ∧ V) = 0.

8.1.2 Lemme de Poincaré

On remarque les deux faits suivants, vrais pour tout ouvert Ω de R3 :

∀f ∈ C 2 (Ω, R), ∇ ∧ (∇f ) = 0

et :
∀V ∈ C 2 (Ω, R3 ), ∇ · (∇ ∧ V) = 0.
Les réciproques sont vraies dans le cas d’un ouvert simplement connexe.
Lemme 24 (Lemme de Poincaré). Soit Ω un ouvert simplement connexe de R3
et V ∈ C 1 (Ω, R3 ).
— Si ∇ ∧ V = 0, alors il existe f ∈ C 2 (Ω, R) tel que V = ∇f .
— Si ∇ · V = 0, alors il existe A ∈ C 2 (Ω, R3 ) tel que V = ∇ ∧ A.

80
8.1.3 Théorème de Stokes
Le but de cette section est de décliner la célèbre formule de géométrie différentielle :
Z Z
dα = α
X ∂X

à l’aide des opérateurs vectoriels définis précédemment.

Proposition 25 (Formule de la divergence). Soit V un volume de R3 et un champ
vectoriel V ∈ Cc1 (R3 , R3 ). On a :
Z Z
V · ndS = ∇ · VdV.
∂V V

Proposition 26 (Formule du rotationnel). Soit S une surface de R3 et un champ

vectoriel V ∈ Cc1 (R3 , R3 ). On a :
Z Z
V · tdℓ = (∇ ∧ V) · ndS.
∂S S

Proposition 27 (Intégration par parties). Soit V un volume de R3 , V ∈ Cc1 (R3 , R3 )

et f ∈ Cc1 (R3 , R). On a :
Z Z Z
V · ∇f dV = f V · ndS − f ∇ · V.
V ∂V V

Proposition 28 (Theoreme du rotationnel). Soit V un volume de R3 et V ∈

Cc1 (R3 , R3 ). On a : Z Z
∇ ∧ VdV = − V ∧ ndS.
V ∂V

8.1.4 Théorèmes de transport

Soit Ω0 un ouvert borné de R3 , V ∈ C 1 (R3 × R, R3 ) un champ de vitesses dont le
flot φ est défini globalement. On note ∀t ∈ R, Ωt = φt (Ω0 ).
Proposition 29 (Transport d’une quantité scalaire). Soit ρ ∈ C 1 (R3 × R, R).

On a :
d Z Z
ρ(M, t)dV = (∂t ρ(M, t) + ∇ · (ρV)(M, t)) dV.
dt M ∈Ωt M ∈Ωt

Proposition 30 (Transport d’une quantité vectorielle). Soit A ∈ C 1 (R3 × R, R3 ).

On a :
d Z Z
A(M, t)dV = (∂t A(M, t) + ∇ · (V ⊗ A)(M, t)) dV.
dt M ∈Ωt M ∈Ωt

Proposition 31 (Transport d’un flux). Soit A ∈ C 1 (R3 × R, R3 ) et Σ0 une surface

bornée de R3 . Pour tout t ∈ R, on note Σt = φ(Σ0 ). On a :
d Z
A(M, t) · n(M, t)dS
dt ZM ∈Σt
= (∂t A(M, t) − ∇ ∧ (V ∧ A)(M, t)) · n(M, t)dS.
M ∈Σt

81
8.2 Théorie du signal
8.2.1 Produit de convolution
8.2.2 Transformation de Fourier
Dans cette section, pour plus de simplicité, on énonce les résultats usuels de la théorie
de Fourier pour les fonctions de la classe de Schwartz :
( )
S(R ) = f ∈ C (R , C), ∀α, β ∈ N , sup x ∂ f (x) < +∞ .
d ∞ d d α β
x∈Rd

On peut montrer que S(Rd ) est une sous-algèbre de CR stable par produit de
d

convolution.
Définition 32 (Transformée de Fourier). Soit f ∈ S(Rd ). La transformée de
Fourier de f , notée fˆ, est définie par :
1
∀k ∈ Rn , fˆ(k) = f (x)e−ik · x dx.
Z

(2π)d/2 Rd

Proposition 33. L’application S(Rd ) → S(Rd ), f 7→ fˆ est bien définie et C-linéaire.

On l’appelle transformation de Fourier.
Proposition 34 (Formule d’inversion de Fourier). Soit f ∈ S(Rd ). On a :
1
fˆ(k)eik · x dk.
Z
∀x ∈ Rd , f (x) =
(2π)d/2 Rn

Proposition 35 (Formule de Plancherel). Soit f, g ∈ S(Rd ). On a :

fˆ(k)ĝ(k)dk.
Z Z
f (x)g(x)dx =
Rd Rd

Proposition 36 (Transformée de Fourier et dérivation). Soit f ∈ S(Rd ) et α ∈ Nd .

On a :
α f (k) = (−ik)α fˆ(k).
∀k ∈ Rd , ∂d
Remarque 37. Ainsi, en un sens, la transformée de Fourier « diagonalise la
dérivation ». C’est pour cette raison qu’elle est particulièrement utilisée dans la
résolution d’équations différentielles.
Proposition 38 (Transformée de Fourier et convolution). Soit f, g ∈ S(Rd ).

On a :
⋆ g = (2π)d/2 fˆĝ.
f[

82
8.2.3 Equation des ondes
Dans ce paragraphe, on se donne c ∈ R∗+ .
Proposition 39 (Résolution explicite de l’équation des ondes en dimension 1). Soit
f ∈ C 0 (R2 , C) et u0 ∈ C 0 (R, C).

On se donne u ∈ C 2 (R)

83
8.3 Equation de Poisson
8.3.1 Existence et unicite de la solution avec condition de Neumann
Soit Ω un ouvert borne de Rd , f ∈ L2 (Ω) et g ∈ L2 (∂Ω) verifiant la relation de
compatibilite : Z Z
f+ g = 0.
Ω ∂Ω
On considere le probleme suivant :

Trouver u ∈ H 2 (Ω) telle que :


−∆u =f
∂n u =g

Il est clair que, si u est solution, alors, pour tout c ∈ R, u + c est encore solution.
On cherche donc les solutions dans un supplementaire dans H 1 (Ω) des fonctions
constantes :  Z 
V = u ∈ H (Ω), u = 0 .
1
Ω
Ainsi, le probleme a resoudre est :

Trouver u ∈ H 2 (Ω) telle que :

=f



 −∆u
=g

Z
∂n u

=0


 u
Ω

Unicite : soit u1 , u2 ∈ H 2 (Ω) deux solutions. On a, en notant u = u2 − u1 :

Z
|∇u|2 = 0,
Ω

donc u est constante, puis u = 0.

Existence : on considere le probleme mis sous forme variationnelle

Trouver u ∈ V tel que :

Z Z Z
∀v ∈ V, ∇u · ∇v = fv + vg
Ω Ω ∂Ω

L’inegalite de Poincare-Wirtinger et le theoreme de Lax-Milgram donnnent

l’existence et l’uniciteZde la solution au probleme precedent, notee u. Soit φ ∈ D(Ω).
1
On note v = φ − φ, de sorte que, d’apres la relation de compatibilite :
|Ω| Ω
1 Z
Z Z Z Z Z  Z Z Z
∇u · ∇φ = f φ + gφ + φ f+ g = fφ + gφ = f φ,
Ω Ω ∂Ω |Ω| Ω Ω ∂Ω Ω ∂Ω Ω

84
donc −∆u = f , puis, d’apres un argument de regularite elliptique non detaille ici,
u ∈ H 2 (Ω). En refaisant le calcul precedent avec φ ∈ D(Ω), on a :
Z Z
φ∂n u = φg,
∂Ω ∂Ω

donc ∂n u = g, ce qui conclut.

8.3.2 Principe du maximum

TODO

8.3.3 Fonction de Green

On appelle fonction de Green de −∆ une distribution G ∈ D′ (Rd ) telle que :
−∆u = δ0 .
Le but de ce paragraphe est de mettre en évidence une fonction de Green de −∆,
et de voir son application a la résolution de l’équation de Poisson.

On cherche les solutions fortes radiales a l’équation de Laplace :

∆u = 0.
On rappelle que, si u ∈ C 2 (Rd , R) est radiale, alors :
d−1 ′
∀r ∈ R∗+ , ∆u(r) = u′′ (r) + u (r).
r
Ainsi, si u ∈ C 2 (Rd , R) vérifie ∆u = 0, alors :





∃A, B ∈ R, ∀r ∈ R∗+ , u(r) = Ar + B si d = 1
∃A, B ∈ R, ∀r ∈ R∗+ , u(r) = A ln(r) + B si d = 2 .


∃A, B ∈ R, ∀r ∈ R+ , u(r) = A + B si d ⩾ 3



 ∗
rd−2
Un calcul non détaillé ici donne des fonctions de Green de −∆ :
1

G = − ln(|x|) si d = 2

2π

1
G = si d ⩾ 3

(d − 2)ωd |x|d−2


2π d/2
ou ωd = est l’aire de Sd−1 .
Γ(d/2)

Il est clair que G ∈ L1loc (Rd ) et que G|Rd \{0} ∈ C ∞ (Rd \ {0}, R).
Théorème 40. On suppose que d ⩾ 3. Soit f ∈ E ′ (Rd ). La distribution G ⋆ f est
une fonction lisse sur Rd \ {0} tend vers 0 à l’infini

A FINIR

Commencer par les solutions a symetrie spherique, puis montrer une solution de
Green, et finir par l’unicite en cas de distribution a support compact

85
8.4 Equation des ondes
8.4.1 Introduction
Le but de ce paragraphe est d’etablir le theoreme suivant :
Théorème 41. Soit S ∈ E ′ (Rd × R∗+ ) et u0 , u1 ∈ E ′ (Rd ). Il existe une unique
distribution u ∈ S ′ (Rd × R) telle que :
1 2

∂ u − ∆u = S + u0 ⊗ δ0′ + u1 ⊗ δ0


c 2 t .
supp(u) ⊂ Rd × R+



Démonstration. Unicite : soit u, v ∈ S ′ (Rd ×R) deux solutions. On note w = u−v,

de sorte que :
 1 ∂ 2 w − ∆w = 0


c2 t .
supp(w) ⊂ Rd × R+


On peut prendre la transformée de Fourier en espace :

1 2
∂ Fx w − |k|2 Fx w = 0,
c2 t
1 1
  
donc ∂t − |k| ∂t + |k| Fx w = 0, donc e2it|k| ∂t (e−it|k| Fx w) est constante par
c c
rapport a t. Or supp(w) ⊂ Rd × R+ , donc e−it|k| Fx w est constante par rapport a t,
donc, puisque supp(w) ⊂ Rd × R+ , Fx w = 0, donc w = 0.

Existence : il suffit de montrer l’existence d’une fonction de Green G de □ a support

dans le futur. Si on dispose d’une telle distribution G, alors G ⋆ (S + u0 ⊗ δ0′ + u1 ⊗ δ0 )
est solution.

86
8.4.2 Fonctions de Green du d’alembertien
1
On cherche une fonction de Green G de 2 ∂t2 − ∆. Par définition, en prenant la
c
transformée de Fourier en temps et en espace :
1
!
ω2
− |k|2 Ft Fx G = −
c 2 (2π) (d+1)/2

Notons que plusieurs choix sont possibles pour (2π)(d+1)/2 Ft Fx G :

1 1 1 1
2 |k| 2 |k| 2 |k| 2 |k|
ω + ω , ω + ω
|k| − + i0+ |k| + + i0+ |k| − + i0+ |k| + + i0−
c c c c
et leurs conjugués. On a donc, pour (2π) d/2+1
Fx G, après un calcul de résidus :
ei|k|ct e−ic|k|t ei|k|ct e−i|k|ct
πic 1R− (t) − πic 1R− (t) et πic 1R− (t) + πic 1R+ (t)
|k| |k| |k| |k|
les conjugués donnant :
ei|k|ct e−ic|k|t ei|k|ct e−i|k|ct
−πic 1R+ (t) + πic 1R+ (t) et − πic 1R+ (t) − πic 1R− (t)
|k| |k| |k| |k|
Dans la suite, on notera G+ , G− , GF , G⋆F définies par :
ei|k|ct e−i|k|ct sin(|k| ct)
!
(2π) d/2+1
Fx G+ = πic − + 1R+ (t) = 2πc 1R+ (t)
|k| |k| |k|
ei|k|ct e−i|k|ct sin(|k| ct)
!
(2π) d/2+1
Fx G− = πic − 1R− (t) = −2πc 1R− (t)
|k| |k| |k|
ei|k|ct e−ic|k|t
!
(2π) d/2+1
Fx GF = πic 1R (t) +
− 1R+ (t)
|k| |k|
ei|k|ct e−ic|k|t
!
(2π) d/2+1
Fx G⋆F = −πic 1R+ (t) + 1R− (t)
|k| |k|
Les deux dernières distributions ne sont pas temperees (nous y reviendrons plus
tard), ce qui nous laisse comme fonctions de Green possibles G± . On remarque
qu’elles sont a support dans le passe et le futur respectivement : supp(G± ) ⊂ R × R.

Dans la suite, on retiendra seulement la fonction de Green a support dans le futur

G+ . On vient d’etablir l’existence de la solution presentee dans l’introduction. On se
donne S, u0 , u1 comme precedemment, et on note u = G+ ⋆ (S + u0 ⊗ δ0′ + u1 ⊗ δ0 ).

On peut expliciter G+ et u lorsque d = 1, 2, 3, S = 0 et u0 , u1 ∈ C 0 (R).

sin(|k|ct)
Dimension 1 : pour tout t > 0, la transformée de Fourier inverse de k 7→
|k|
√ 1
est 2π 1[−ct,ct] . Ainsi :
2
c
∀(x, t) ∈ R2 , G+ (x, t) = 1[−ct,ct] (x)1R+ (t).
2

87
Par consequent :

c Z Z ct′
∀(x, t) ∈ R × R∗+ , G+ ⋆ (u0 ⊗ δ0′ )(x, t) = u0 (x − x′ )δ0′ (t − t′ )dx′ dt′
2 R+ −ct′
c
= (u0 (x − ct) + u0 (x + ct))
2
et
cZ
∀(x, t) ∈ R × R∗+ , G+ ⋆ (u1 ⊗ δ0 )(x, t) = 1[−ct,ct] (x′ )u1 (x − x′ )dx′
2 R
c Z ct
= u1 (x − x′ )dx′
2 −ct
c Z x+ct
= u1 (x′ )dx′ ,
2 x−ct
donc :
c c Z x+ct
∀(x, t) ∈ R × R∗+ , u(x, t) = (u0 (x − ct) + u0 (x + ct)) + u1 (x′ )dx′ .
2 2 x−ct
sin(|k|ct)
Dimension 2 : pour tout t > 0, la transformée de Fourier inverse de k 7→
|k|
1
est x 7→ q 1[0,ct] (|x|), donc :
(ct)2 − |x|2

c 1
∀(x, t) ∈ R2 × R, G+ (x, t) = 1[0,ct] (|x|)1R+ (t).
2π (ct)2 − |x|2
q

Par conséquent :
c Z Z 1
∀(x, t) ∈ R2 × R∗+ , G+ ⋆ (u0 ⊗ δ0′ ) = u0 (x − x′ )δ0′ (t − t′ )dx′ dt′
2π R+ |x′ |⩽ct′ (ct′ )2 − |x′ |2
q

Parler de la formulation au sens des distibutions temperee, moyennes spheriques,

Huygens 2d, resolution en 3d, en Fourier, etc

88