Chapitre 3 : Théorie fluide et cinétique du plasma

Au chapitre 2, nous avons analysé le mouvement d’une particule chargée isolée, c’est-à-dire
sans interaction aucune avec d’autres particules, et soumise à des champs électrique E et
magnétique B appliqués de l’extérieur. Dans le présent chapitre, nous abordons un modèle qui
considère un ensemble de particules.
Cependant, ce que nous observons ce sont paramètres macroscopiques moyens tels que la
densité de courant j, la densité d’électrons ne, la pression P, la température des ions Ti, etc…
Ces paramètres sont des moyennes macroscopiques de la fonction de distribution des vitesses
et/ou de positions. Cette fonction de distribution s’obtient dans le cadre de la théorie cinétique
des gaz.
Dans ce chapitre, nous allons:

 Introduire les concepts de fonction de distribution fα(r, v, t) pour une espèce

de plasma.
 Déterminer les équations d’évolutions temporelles de la fonction de distribution
(équation de Vlasov/Boltzmann).
 Déterminer les différents moments de l’équation de Vlasov.

I. Eléments de la théorie fluide

La théorie fluide considère le plasma comme un fluide se mouvant avec une vitesse u.
Prenons un petit volume V, dans un plasma de masse volumique La masse totale et la
quantité de mouvement relié à ce volume sont:
(1)

Les équations fluides apparaissent en explicitant les lois de conservation de la masse et de

l'impulsion.

(2)
F représente les forces externes (pesanteur, pression, forces électriques) exercées sur le
volume.

Figure 1: Evolution d’un élément de surface

L’évaluation de la variation temporelle de  f (t )dV

V
fait apparaitre deux termes:


 L’un du à la variation temporelle de f, soit  t f (t )dV .
V
 
 L’autre du au déplacement de la surface qui entoure le volume V, soit
V
 f u dS .
Nous obtenons alors, à titre d'exemple, pour la dérivée totale sur la masse l'expression
suivante:

(3)
En utilisant le théorème de la divergence nous obtenons ensuite:
(4)

Le volume V étant quelconque, la dernière équation est satisfaite si:

(5)

C'est l'équation de continuité pour un fluide, qui représente la conservation de la masse.

La même démarche sur l'équation de Newton, suivie par l'application du théorème du
gradient, nous donne, en l'absence des forces électromagnétiques:

(6)
L'opérateur s'appelle dérivée convective et on le note souvent, par
commodité, d/dt.
Si nous considérons aussi les forces électromagnétiques, chose tout à fait naturelle et
indispensable dans un plasma, il faut rajouter un terme à l'équation de Newton.
Avec les forces de Lorentz l'équation devient, pour l'espèce :
(7)

Dans notre système d'équation, Newton + continuité, nous avons cinq inconnues (u, p, ,E et
B), il faut donc le compléter avec trois équations supplémentaires pour que le problème soit
défini. On utilise alors une équation d'état pour relier pression et densité, et les équations de
Maxwell pour B et E.
Nous obtenons finalement, pour la théorie fluide, le système d'équations suivant, où =i,e
dénote l'espèce de particule:
CONTINUITE:

NEWTON:

ETAT:

MAXWELL:
Les deux autres équations de Maxwell sont:

(où est maintenant la densité de charge) n'apportent pas une information supplémentaire; on
peut en effet les considérer comme conditions initiales des deux premières.
Les densités de charge et de courant sont facilement exprimables en fonction des densités des
particules et des vitesses:
II. Théorie cinétique
II.1 Espace de phase
Considérons une particule d’espèce α (α=e ou α=i). Elle peut être repérée par sa position
r  xi  y j  z k et le vecteur vitesse v  vx i  v y j  vz k . Les coordonnées ( r, v ) définissent la
position de la particule dans l’espace de phase.
Pour un système de plusieurs particules, on introduit la fonction de distribution fα( r, v ,t) tel
que: fα( r, v ,t) d r d v = dN( r, v ,t) est le nombre de particules dans l’élément de volume
dV= d r d v de l’espace de phase, avec d r ≡ d 3r ≡ dxdydz et d v ≡ d 3v ≡ dvx dv y dvz .
fα( r, v ,t) est une fonction finie, positive qui décroit vers zéro lorsque v   .

Figure 2: Un élément de l’espace de phase de largeur dx et de hauteur dv.

II.2 Fonctions de distribution et valeurs moyennes

En un temps donné, chaque particule a une position et une vitesse spécifique. On peut alors
caractériser la configuration temporelle d’un grand nombre de particules en spécifiant la
densité des particules à chacun des points r et v de l’espace de phase.
La fonction permettant de décrire l’évolution instantanée de la densité des particules dans
l’espace de phase est appelée fonction de distribution et noté f( r, v ,t).
Ainsi, f( r, v ,t) d r d v est le nombre de particules à l’instant t, ayant une position comprise
 
entre r et r  dr , et une vitesse appartenant à l’intervalle v et v  dv .
La valeur moyenne < X > d’une quantité X est définie à partir de la fonction de distribution :
  

X (r , t )   X f  ( r , v , t ) d 3v
  .
3
 f ( r , v , t ) d v
(8)
Cette quantité permet de calculer les grandeurs physiques du plasma, telles que la densité de
particules, la pression, la vitesse moyenne, etc…
II-3 Moments de la fonction de distribution
Les variables macroscopiques (quantités mesurables) sont obtenues en considérant les
moments des fonctions de distribution de vitesses appropriés.
D’une manière générale, M (n ) le moment d’ordre n de la fonction de distribution fα( r, v ,t) est
défini par :
    
M ( n ) (r )   v n f (r , v , t ).d 3v , (9)
où n est un entier.
On définit les grandeurs fluides suivantes:
-Densité des particules (nombre de particules par unité de volume) :
La première grandeur physique accessible est la densité de particules n( r, t ) liée au moment
d’ordre 0 :
   

n(r , t )  f (r , v , t ).d 3v. (10)


-La vitesse moyenne, ou vitesse fluide noté u est lié au moment d’ordre 1 de f:
En appliquant à v la définition de la valeur moyenne donnée par l’équation (8), on obtient :
       
 
u (r , t )   v f (r , v , t ).d 3v  v f (r , v , t ).d 3v
    .
3
 f (r , v , t ).d v n( r , t )
(11)
– Pression cinétique
Une autre quantité importante pouvant être extraite de la fonction de distribution des
particules est la pression cinétique. Cette pression cinétique a pour origine, pour un élément

de volume se déplaçant à la vitesse moyenne u , les mouvements aléatoires des particules qui

entrent et sortent de cet élément de volume. La vitesse fluctuante w est définie comme la
vitesse d’une particule dans le référentiel de la vitesse moyenne : w  v  u.
Dans le cas d’un milieu isotrope, la pression cinétique est alors définie par :
      
P(r , t )  m w2 f (r , v , t ) d 3v  m v  u  f (r , v , t ) d 3v ,
2

     
 P(r , t )  m v 2 f (r , v , t ).d 3v  m n(r , t ) u 2 (r , t ).
(12)

II-4 Equation de Vlasov (Boltzmann)

Il a été vu dans la section précédente que la fonction de distribution permet de calculer les
grandeurs caractéristiques du plasma. Il reste maintenant à trouver l’équation permettant de
déterminer l’évolution de f( r, v ,t) en fonction des conditions physiques auxquelles sont
soumis les particules.
Cette évolution temporelle de f donne une description du système beaucoup plus détaillée que
la description fluide, mais moins détaillée que si l’on suivait la trajectoire individuelle de
chaque particule. Donc l’évolution de f pour caractériser le système ne permet pas de suivre la
trajectoire individuelle des particules, mais plutôt de caractériser une classe de particules
 
ayant la même position r et la même vitesse v .
A l’instant t, les particules se trouvent dans un volume D à
 
la position ( r , v ) de l’espace de phase.
En l’absence de collisions, l’écoulement dans l’espace des
phases est celui d’un fluide incompressible de sorte que, à
un temps t+dt ultérieur, ces mêmes particules de l’élément
de volume d r d v vont se retrouver au point r  v dt ,
 
 
v  d v / dt dt de l’espace des phases où d v / dt  F / m
représente l’accélération produite par les forces imposées de l’extérieur. En développant en
série de TAYLOR limitée au premier ordre, la fonction de distribution en ce nouveau point
par rapport au point initial, nous obtenons:
 
  f f Fi f
3 3
F
f  r  v dt , v  dt , t  dt   f r , v, t   vi dt  dt  dt , (13)
 m  i 1 xi i 1 vi m t
Nous pouvons donc exprimer la variation totale, df/dt, de f entre les deux points, par tranche
élémentaire de temps dt, sous la forme :
 
df 3
f 3
f Fi f    F  f
 vi    v . r f  .v f 
dt i 1 xi i 1 vi m t m t , (14)
 
où  r et  v sont respectivement les opérateurs différentiels dans l’espace des coordonnées
 

spatiales et dans l’espace des vitesses.

En l’absence de collisions, le nombre de particules contenues dans un élément de volume de
l’espace des phases est conservé et la fonction de distribution n’est donc pas modifiée, de
sorte que df/dt =0,  d’où :
f    F 
 v . r f  .v f  0.
t m (15)
C’est l’équation de VLASOV (ou équation de BOLTZMANN sans collisions).
En présence des collisions binaires, l’équation précédente devient:
f    F   f 
 v . r f  .v f    .
t m  t coll (16)
On obtient ainsi l’équation de BOLTZMANN qui décrit l’évolution de la fonction de
distribution des particules dans l’espace des phases sous l’influence, d’une part, des gradients
affectant cette distribution et, d’autre part, des forces en présence et des collisions.
En général, l’équation de Vlasov à 3D pour chaque particule est couplée aux équations de
Maxwell:

II.5 Moments de l’équation de Vlasov

A partir de l’équation de Vlasov que nous avons obtenue dans la section précédente, il est
possible de retrouver les équations fluides.
Considérons l’équation de Vlasov pour l’électron:
f
t r
 
  .v f    .
v  m

  e    

E  v  B f   0.
 (17)
 Premier moment : conservation du nombre de particules
Le premier moment de l’équation de Vlasov s’écrit:
f 3 
v t d v 
  3
v r
 .v f d v   
  e    
 v . m E  v  B f   0
(18)

    

t v
fd 3v   . v f d 3v  0.
r v
Vues les définitions ci-dessus, on a en fait montré que

(19)
c’est-à-dire l’équation de conservation du nombre d’électrons.
  Second moment : conservation de la quantité de mouvement
Pour aboutir à l’équation de conservation de la quantité de mouvement, on prend le second
moment:
 f 3 
v t
v d v 
   3
v r
v  .v f d v  v 
   e     

v v . m E  v  B f   0.
(20)
Pour le premier terme:

Pour le second terme, on cherche à faire intervenir le terme de pression:

   3  3
v v r .v f d v  r .v v f d v
2

 P 
  .  n u 2 .
r  m 
Pour le troisième terme, on utilise une intégration par partie:

On a donc montré que:

 
t
   
n u     P  n u 2   e n E  n u  B  0
r  m

 m (21)
 
 n
u 1 P

t m r
 u   n    e 
   
  n u   u    n u   n E  n u  B  0,
r

 t r  m
 n   
ce qui conduit en utilisant la relation de continuité    n u   0, à:
 t r 
 
m n
 u  u  P

 u       en E u  B .
 

 t r  r (22)
On obtient cette fois l’équation de conservation de la quantité de mouvement.
 Troisième moment : conservation de l’énergie
Nous ne le faisons pas ici, mais il est ensuite possible de prendre le troisième moment pour
aboutir à une équation de conservation de l’énergie:

etc…
III. La distribution de Maxwell-Boltzmann
Une fonction de distribution importante est la distribution de Maxwell-Boltzmann:
3/ 2 
   m    m v2     
f r , v , t     nr , t exp    , de tel sorte que :  f r , v , t  d 3v  nr , t  .
 2K BT   2 K BT  
Dans certains cas, on utilise la fonction de distribution de Maxwell-Boltzmann normalisé à
l’unité :
En utilisant , on peut vérifier que (distribution
normalisée).
Propriétés des intégrales:

Deux moyennes importantes, souvent calculées pour la distribution de Maxwell Boltzmann: la

vitesse moyenne et l'écart quadratique moyen permettant d'obtenir l'énergie:
 La vitesse moyenne v est définie par:

Comme est symétrique suivant vx, vy, vz, l’intégrale s’obtient en passant en coordonnées
sphériques

La norme de la vitesse moyenne dans une direction a une moyenne différente que la
composante dans une direction .

 L’écart quadratique moyen des vitesses à une dimension s’obtient par:

En résumé, une fonction de type Maxwell-Boltzmann (Gaussienne) possède les propriétés

suivantes:

Démontrer les relations précédentes.