Chapitre II plasma en équilibre thermodynamique

Chapitre II

Plasma en équilibre thermodynamique

II.1. Introduction :

L’équation de Boltzmann (1872) est une équation intégro-différentielle de la théorie

cinétique qui décrit l'évolution d'un gaz peu dense hors d'équilibre. Elle permet notamment
d'étudier la relaxation du gaz d'un état d'équilibre local vers l'équilibre global caractérisé par
la distribution de Maxwell des vitesses.

II.2. Equation de Boltzmann non-relativiste :

L’équation de Boltzmann est une équation de continuité de la fonction de distribution

 
dans l’espace des phases d r d v .

On peut dire que la variation de la fonction de distribution dans l’élément de volume

dans l’espace des phases est due au flux de particules qui entrent et qui sortent et qui
provoquent des collisions dans cet élément de volume. Ce que l’on exprime en disant que la
variation de la fonction de distribution est due seulement aux collisions.

       
df  r , v , t   f  r , v , t  
     (II-1)
dt  t 
 
  collisions

  
df  r , v , t 
   f  f r  f v
dt t r t v t (II-2)

f   F   f 
  v . r f  . v f   
t m  dt col

26
Chapitre II plasma en équilibre thermodynamique

 Signification des différents termes

 
  
 v . r f  r , v , t  : exprime l’influence des phénomènes de diffusion ; par exemple :
 
considérons un nuage de particules électrons et aucune particule à l’extérieur e nuage. Du
faite de la vitesse des électrons cet état va évaluée les molécules vaut diffusées vers les
régions vides et au bout de certain temps l’ensemble tendras vers une situation homogène.

F    
 .  v f  r , v , t  : exprime l’action des forces appliquées.
m  

F
 : est l’action imposée aux particules par les forces d’origine
m
extérieure (champ électrique ou magnétique par exemple).
 f 
   Terme de collisions.
  t  col
II.3. Etablissement des équations de transport :
Le modèle hydrodynamique que nous utilisons pour analyser la décharge dans un gaz
impose au plasma d’agir de façon globale ou encore de se comporter comme un fluide ±
conducteur et non pas par l’intermédiaire des particules individuelles. La condition de validité
de cette approximation est que la distance entre les particules soit faible par rapport aux
distances où agissent les forces inter particules qui se traduit par l’équation :

n3D  1 (II-3)

 D : La longueur de Debye.

n : la densité des particules.

et que le libre parcourt moyen ( l..P.m.) Soit faible par rapport à l’échelle des variations de
quantités macroscopiques.

27
Chapitre II plasma en équilibre thermodynamique

  
II.3.1. équation de transport d’une quantité : A  r , v , t 
 

Ces conditions en générale sont vérifier dans les plus parts des décharges électriques
dans les gaz et pratiquement la majorité des processus de base produisent dans les décharges

peuvent être traités sous la base des modèles hydrodynamique.

  
Les équations macroscopiques se déduisent de l’équation de Boltzmann. Soit f  r , v , t  la
 
fonction de distribution qui obéit à l’équation de Boltzmann :


f   F   f 
 v .r f  .v f    (II-4)
t m  dt  col

Qui s’écrit aussi

f f F f  f 
  .vi  i .   (II-5)
t i xi i m  vi  dt  col

i =1, 2, 3.

  
Soit la quantité A  r , v , t  une fonction de la vectrice vitesse, de la vectrice position et
 

  
du temps. Multiplions l’équation (II-5) par A  r , v , t  , et effectuant une intégration sur tout
 
l’espace des vitesses. Ce calcul est développé dans [39, 40] ; il conduit au résultat suivant :


n  A   n  A  . n  vA  n  v.A  n.  F A  C  A (II-6)
t t m v

28
Chapitre II plasma en équilibre thermodynamique

  
Dans lequel on a introduit la valeur moyenne A  r , v , t  de la grandeur A définit par
 
l’équation :

  
 
1
n
  
 A  r , v , t    A r , v, t f r , v, t d v (II-7)

et le terme :

 f  
C  A    A  d v (II-8)
 t  col

L’équation (II-8) est l’équation de transport de la grandeur A. En remplacent A par

diverses fonctions de plus en plus compliquées de la vitesse, nous obtiendrons une suite
d’équations hydrodynamiques.

 
Lorsque A est fonction seulement de v (mais non de r ni de t) l’équation de transport
(II-8) s’écrit plus simplement :


n  A   . X  A  n  F A  C  A (II-9)
t m v

Avec :


X  A  n v A (II-10)

Sous cette forme on connaît une équation de conservation de grandeur A : n A et X (A)

sont la densité et le flux de cette grandeur en un point du fluide. Les deux termes sources qui
figurent au second membre représentent la création éventuelle d’une certaine quantité de cette
grandeur par unité de temps et de volume respectivement par la force extérieure et par les
interactions entre particules.

29
Chapitre II plasma en équilibre thermodynamique

II.3.2 Les équations fondamentales de conservation :

II.3.2.1.L’équation de continuité :

On pose dans l’équation (II-9) A≡1 et C (1)=0. On obtient l’équation de conservation

des particules :

n 
 .n v  0 (II-11)
t

Dans un gaz réactif il y aurait des collisions inélastiques pouvant créer ou détruire des
particules et on aurait alors au deuxième membre un terme de source positif ou négatif.

II.3.2.2.équation de conservation du moment :


  
Posons maintenant dans (II-9) A  m v et C  m v   p  . . On obtient :
 


    F     
 nm   .nm  v v   n   m v   C m v  (II-12)
t   m 
v    

 
En tenant compte de C  m v  et de l’identité :
 

   

  m v   m  (II-13)
v 

On obtient l’équation de conservation de la quantité de mouvement :

 
   
 nm    .nm  v v  n  F    (II-14)
t  

30
Chapitre II plasma en équilibre thermodynamique

II.3.2.3.équation de conservation de l’énergie cinétique :

1
Posons enfin A  mv 2 dans l’équation (II-9) ; on obtient l’équation de conservation
2
de l’énergie cinétique :


U K   . S K  n  F    C  1 mv 2 
  
(II-15)
t 2 


S K est le vecteur flux d’énergie cinétique total. Au deuxième membre apparaissent
deux termes de source : le premier est le travail moyen effectué par la force extérieure (la

valeur moyenne marquée sur  F  est en fait inutile ; pour les forces indépendantes de la
 
vitesse on a  F   F ; les seules forces dépendantes de la vitesse sont les forces de Lorentz

qui, étant perpendiculaire à v , apportent une contribution nulle.), le deuxième terme
représente l’effet des interactions.

II.4.Formalisme hydrodynamique :
II.4.1 Formalisme hydrodynamique du gaz ionisé :
L’intégration de l’équation générique

avec donne l’équation de transport de quantité de mouvement :

(II-16)

est le tenseur de pression cinétique décrivant la densité d’énergie d’agitation thermique.

Si on néglige le terme de dérivée temporelle de l’énergie de dérive devant l’énergie

thermique, on obtient alors l’équation de conservation de la quantité de mouvement simplifiée
définissant l’espace vitesse des espèces chargées

31
Chapitre II plasma en équilibre thermodynamique

(II-17)

Où est le flux de particules ; le terme de conduction électrique ( étant la

vitesse de dérive de l’espèce) ; est le terme de diffusion ( étant le tenseur de diffusion).

L’hypothèse de champ local permet, en plus de préciser la dépendance vis-à-vis du champ.

Finalement, l’équation (II-19) couplée à l’équation de conservation de la densité des

électrons et des ions, et à l’équation du champ électrique, définissent le modèle fluide d’ordre
un utilisé dans ce travail, résumé par :

(II-18)

Avec

Où Se est le taux de création ou de disparition de particules (suite à l’ionisation,

l’attachement, etc.).
Remarque : pour des raisons de simplicité, nous avons choisi un modèle du premier ordre
pour les particules chargées, qu’un modèle du second ordre incluant en plus de l’équation de
l’énergie des électrons. Le modèle d’ordre deux est plus précis [41 ,42] mais en raison surtout
du manque de données de base et de l’allongement des temps de calcul, nous avons préféré
l’utilisation d’un modèle d’ordre beaucoup plus souple lors d’un couplage de modèle gaz
ionisé-gaz neutre.

II.4.2 Formalisme hydrodynamique du gaz neutre :

II.4.2.1.Equation de continuité :
Le gaz considéré est faiblement ionisé impliquant donc que l’ionisation, la
recombinaison, etc. ne modifient pas de manière perceptible la population des neutres (avec N
la densité des particules neutres), dans ce cas, l’équation de conservation de densité s’écrit :

32
Chapitre II plasma en équilibre thermodynamique

(II-19)

II.4.2.2.Equation de conservation de la quantité de mouvement :

L’équation de la conservation de la quantité de mouvement, du fait que les forces
appliquées sur les neutres soient nulles, superposée aux conditions de fermeture peut se
simplifier à :

) (II-20)

Le terme de collision représente les transferts de quantité de mouvement lors du choc

dus aux composantes du gaz ionisé (S : électrons et ions) avec les particules neutres, ie,
. Les théories classiques des collisions permettent d’écrire :

Où est la fréquence de collision élastique entre les neutres et l’espèce (électron ou ion).
Cette fréquence représente le nombre de fois qu’une particule neutre donnée, entre en
collision avec une particule de l’espèce chargée s par unité de temps. Si N T représente le
nombre d’interactions par unité de temps entre le gaz de particules neutres et le gaz de
particules chargées, alors

Où est la fréquence de collisions élastiques entre l’espèce s et les particules

neutres, impliquant que le terme de collision devient :

Les fréquences de collisions peuvent s’écrire différemment en faisant appel à la

définition de la mobilité et la diffusion

33
Chapitre II plasma en équilibre thermodynamique

avec

Impliquant donc

L’équation de conservation des neutres comprend :

(i) un terme de quantité de mouvement dû au choc direct avec les particules chargées.
(ii) un terme des effets visqueux.
La notion de viscosité, apparentée au frottement interne du fluide, n’apparaît que
lorsque les différentes portions de celui-ci s’écoulent avec des vitesses différentes et se
déplacent donc les unes par rapport aux autres.
La viscosité se manifeste par le transport supplémentaire de la quantité de mouvement,
à partir des endroits de grandes magnitudes vers ceux des plus petites. Cela implique donc que
le tenseur des contraintes visqueuses doit dépendre des dérivées de la vitesse par rapport aux
coordonnées spatiales. Si les gradients ne sont pas de trop grandes amplitudes, la dépendance
de ce terme de transport de quantité de mouvement se fait par rapport aux dérivées premières
de la vitesse d’après L.Landau [43]. est le tenseur des tensions visqueux dont la

composante (par exemple en cordonnées cartésiennes) pour un fluide compressible s’écrit

de la forme générique

Le coefficient d’échange moléculaire

II.4.2.3.Equation de conservation de l’énergie :

Le terme collisionnel de l’équation de l’énergie sert à accroître directement l’énergie
de translation d’agitation thermique dans le cas d’un gaz neutre[44,45].
L’énergie ε des neutres est composée d’une partie d’agitation (due au mouvement
d’ensemble des particules) et d’une fraction thermique (due à la vitesse d’agitation
thermique)
34
Chapitre II plasma en équilibre thermodynamique

et c-à-d

(II-21)

Le taux de gain d’énergie par unité de volume est fonction d’énergie par conservation
et conduction. Il est également dû au travail des forces de pression et de viscosité.
Dans le cas d’un gaz moléculaire, l’énergie déposée dans le milieu résistif neutre par unité de
temps est distribuée en fonction de l’état statistique du gaz et des différentes sections efficaces
de collisions dans plusieurs réservoirs d’énergie dont certains relaxent vers le réservoir
thermique en des temps assez longs (supérieurs à la durée de la décharge). Les différents
processus collisionnels élastique et inélastique ont pour conséquence de déposer l’énergie qui
est stockée sous plusieurs formes correspondantes aux processus : de translation, de rotation,
de vibration, d’excitation, de dissociation et d’ionisation. A fin de prendre en compte la
fraction d’énergie qui est effectivement relaxée sous forme thermique pendant la durée de la
décharge, le terme source d’énergie apportée aux neutres, est pondéré par un coefficient t

Ft, i.e. .
La modélisation de l’action de la décharge sur les particules neutres est composée de
trois équations de conservation (II-21 à II-23) où l’interaction entre particules neutres et
particules chargées est prise en compte par le biais des termes explicités dans le prochain
paragraphe.

II.5.Couplage de la dynamique des particules chargees a la dynamique des

particules neutres :
Le modèle auto cohérent décrivant la dynamique des particules chargées et la
dynamique des particules neutres d’une décharge comprend le système d’équations
étroitement couplées (II-20 à II-23) et l’équation du champ électrique. Le système d’équations
illustre clairement que le modèle peut être tributaire des particules chargées et des neutres à la
double hétérogénéité
(Charge d’espace et dépollution des neutres).

35
Chapitre II plasma en équilibre thermodynamique

II.6. Processus élémentaires dans un plasma :

On ne peut pas décrire l’équilibre thermodynamique sans avoir étudié tous les
processus élémentaires qui se produisent dans un plasma. En effet le plasma a le siège de
divers processus élémentaires qui décrivent les interactions entre ces particules. Ces
phénomènes sont soit collisionels soit radiatifs.

II.6.1.Processus collisionels :

En dehors des collisions élastiques qui sont les causes des transferts d’énergie entre
les particules et l’établissement d’un équilibre thermique, d’autres collisions inélastiques
peuvent se produisent dans le plasma

II.6.1.1. Excitation et désexcitation :

Sous l’effet d’un choc ou collision une particule A (neutre ou chargée) peut absorber
une énergie cédée par la particule B (qui peut être aussi un électron). Cette énergie va
modifier l’état interne de la particule A et permettre à un électron de son cortége électronique
de passer d’un niveau d’énergie En à un niveau supérieur d’énergie Em. On dit alors que
cette particule a été excitée.

B + A (n) 
 B + A*(m)

Avec

n : le niveau initial ; m : le niveau final.

B : la particule incidente.

A : la particule cible dans l’état d’énergie inférieur En qui n’est pas nécessairement l’état
fondamental

A* : la particule excité (dans l’état d’énergie Em)

36
Chapitre II plasma en équilibre thermodynamique

Inversement, un atome dans un niveau supérieur peut se désexcité spontanément

vers le niveau fondamental, en libérant un photon h = E. une telle transition est connue
sous le non de l’émission spontanée.

Par ailleurs, ces processus ne sont pas toujours réalisables car il existe des niveaux
métastables dont la durée de vie est de l’ordre de 3 10-3 s. Les atomes restent piégés jusqu'à ce
que une collision électronique arrive les désexcité vers leurs niveaux fondamental, ce
phénomène corresponds à l’émission induite. Les états métastables sont importants dans un
plasma, parce qu’ils peuvent stoker de l’énergie interne, et la livrer à une autre espèce dans le
plasma ou à une surface [46] .

II.6.1.2. Ionisation et dissociation :

Sous l’effet d’une collision un atome ou une molécule peut libérer un électron (le plus
éloigner du noyau), lorsqu’elle absorbe suffisamment d’énergie. Le phénomène observé : est
l’ionisation. Son processus inverse s’appel la recombinaison

A + e-    A+ + e- + e-
ionisation

La dissociation d’une molécule est un phénomène qui nécessite un apport d’énergie

qui peut être cédé à la molécule par une autre particule lors d’une collision inélastique.

AB 
 A + B
E AB

Avec E AB est l’énergie de dissociation.

II.6.2. Processus radiatifs :

Plusieurs mécanismes sont responsables de l’émission du rayonnement total, leurs

importances relatives dépendent de la température, de la fréquence et de la nature des atomes
émetteurs. Les différents processus radiatifs sont les suivants :

37
Chapitre II plasma en équilibre thermodynamique

II.6.2.1. Transition liée - liée (spectre des raies) :

Un électron dans un atome peut effectuer une transition d’un niveau d’énergie
supérieur Em vers un autre niveau d’énergie inférieur E n, en émettant un photon de fréquence
proportionnelle à la variation de l’énergie :
hνmn = Em - En

II.6.2.2. Transition libre - lié:

Un électron libre peut être capturé par un ion de charge Ze en émettant un photon hν
Il forme ainsi un atome neutre ou un ion de charge (Z-1) e se phénomène correspond à la
recombinaison radiative Son processus inverse est la photo - ionisation.

A+ + e- 
 A + hν

Si la particule formée est dans son état fondamental, l’équation de la conservation d’énergie
s’écrit sous la forme suivante :

1
me ve2  Eion  h
2

Certains atomes neutres peuvent capturer un électron et former un ion négatif. Le

champ moyen du noyau positif de l’atome neutre s’étend au delà de l’écran formé par le
cortège électronique. Les gaz rares qui ont une couche externe saturé ont un écran
particulièrement efficace et ne forme pas donc d’ion négatif (sauf l’hélium). Par contre les
halogènes (Cl, F,…) auxquels il manque un électron pour avoir une couche externe saturée
ont le plus for champ attractif et forma ainsi des ions négatifs. Le phénomène de la capture
d’électron est accompagné par l’émission d’un photon ou d’un rayonnement continue ;

A+ + e- 
 A + hν

38
Chapitre II plasma en équilibre thermodynamique

1
me ve2  E a  h
2

Ea représente l’affinité électronique de l’atome A

Le processus inverse est la photo détachement.

II-6-2-3.Transition libre – libre :

Un électron peut céder une partie de son énergie cinétique, lorsqu’il est freiné dans un
champ d’un atome ou d’un ion (en restant libre). Ce processus est accompagné d’une
émission d’un photon.

A + e-vi 
 A + e-vf + h

Avec

m
h = (Vi2 – Vf2)
2

m : la masse de l’électron.

Vi : vitesse initiale de l’électron.

Vf : vitesse finale de l’électron.

II.7. Equilibre thermodynamique complet (etc) :

Le plasma peut atteindre un état d’équilibre thermodynamique complet en absence des

gradients de température est de pression. On peut alors parler d’une température unique, le
théorème de l’équipartition est vérifier : les différentes espèces de particules du milieu ont
alors la même énergie cinétique moyenne. Il s’ajoute à cette condition un équilibre entre les
phénomènes de création et de disparition comme l’ionisation et la recombinaison, l’excitation
et la désexcitation ou l’émission et l’absorption de photon. Autrement dit chaque processus
est contrebalancé par son processus inverse. Dans ces conditions le plasma doit suivre les lois
d’équilibre [47].

39
Chapitre II plasma en équilibre thermodynamique

II.8.Les lois d’équilibre :

Nous pouvons compter cinq lois importantes dont quatre régissent l’émission du
rayonnement thermique à l’équilibre : la loi de Maxwell, la loi de Boltzmann, la loi de Saha et
la loi de Guldberg-Waage, la loi de Planck donne la distribution spectrale. Comme on peu
ajouter des lois complémentaire tel que : la loi de neutralité électrique et la loi de Dalton.

II.8.1. La loi de Maxwell :

A l’équilibre la loi de distribution de vitesse est celle de Maxwell, elle permet de

déterminer le nombre de particules dN d’une espèce donnée dont la vitesse est comprise entre
v et v + dv (pour une température fixée).

 mv 2 
3/ 2
dN  m 
 f (v)dv    4v 2 exp   (II-22)
N  2KT   2 KT 

m : étant la masse de la particule et K la constante de Boltzmann.

II.8.2. La loi de Boltzmann :

Dans un système composé de plusieurs particules en équilibre thermodynamique à une

température T, la distribution canonique de l’équation de Boltzmann permet de donner les
populations des différents niveaux d’énergies.

N  E 
Nm  g m exp  m  (II-23)
U (T )  KT 

II.8.3. La loi de Saha :

Elle régie l’équilibre d’ionisation dans un plasma, elle peut donc donner le degré
d’ionisation du milieu.

40
Chapitre II plasma en équilibre thermodynamique

U  2me KT   E  E 
3/ 2
Ne Ni
2 i   exp  i  (II-24)
Na U a  h2   KT 

Ne, Ni et Na sont respectivement les densités de l’électron, de l’ion et de l’atome neutre.

Ui (T) et Ua (T) sont les fonctions de partition de l’ion et de l’atome, celle de l’électron est
égale à 2.

h est la constante de Planck.

Ei est le potentiel d’ionisation de l’espèce considérée.

∆E est l’abaissement du potentiel d’ionisation du à la présence des particules chargées, son

expression est donnée par Bruner [47] :

Ne
E  1.2110 6 3 N e  2.510 8 (II-25)
T

II.8.4. La loi de Planck :

Lorsque le plasma est en équilibre thermodynamique pour une température T,

chaque processus élémentaire est équilibré par son processus inverse. Cette condition exige
que le rayonnement émis par le plasma est totalement absorbé. Autrement dit, le rayonnement
émis subit une certaine distribution spectrale, c’est celle du corps noir donné par Planck.

2h 3 1
B0  n 2 (II-26)
c2  h 
exp  1
 KT 

B0 : puissance monochromatique émise par unité de surface, d’angle solide et de fréquence.

41
Chapitre II plasma en équilibre thermodynamique

c : vitesse de la lumière dans le vide.

h : constante de Planck

II.8.5. La loi de Guldberg – Waage :

La loi de Guldberg - Waage lie les espèces moléculaires aux espèces atomiques.
Elle exprime également l’équilibre de dissociation des espèces moléculaires. C’est une autre
forme de la loi de Saha.

3/ 2
 m A mB   2KT   E  E AB 
3/ 2
N e N i U AU B
    2  exp  AB  (II-27)
Na U AB  m AB   h   KT 

UA, UB et UAB sont les fonctions de partitions de A, B et AB, mA, mB et mAB leurs masses
respectivement. EAB est l’énergie de dissociation et ∆E AB l’abaissement de potentiel de cette
énergie de réaction.

AB 
 A  B
E AB

II.9. Fonction de partition :

La fonction de partition joue un rôle très important dans la physique des plasmas, en
dehors de son utilisation dans le calcul de la composition elle permet de déterminer plusieurs
propriétés thermodynamiques des plasmas. Cette fonction à été introduite par la statistique de
Boltzmann :

  Etot 
U tot   g tot . exp  (II-28)
 KT 

42
Chapitre II plasma en équilibre thermodynamique

La fonction de partition totale d’un système c’est le produit de plusieurs contributions :

- la fonction de partition de translation elle est donnée par la statistique classique :

 2mi KT 
3/ 2

U Tran
i (T )    (II-29)
 h 
2

-la fonction de partition interne

  En 
U iInt T    g n  exp  (II-30)
n  KT 

où U T i
Int
est la fonction de partition interne du l’espèce i. En représente l’énergie interne
(quantifiée) de la particule. Si la particule est monoatomique son énergie interne se réduit à
l’énergie d’excitation électronique, sa fonction de partition interne est égale à :

  Ek 
U ielec T    g k  exp  (II-31)
n  KT 

S’il s’agit d’une molécule polyatomique il faut ajouter à l’énergie électronique les énergies de
vibration et de rotation. Ainsi la fonction de partition interne peut se mettre sous la forme d’un
produit de ces trois contributions [48]:

U Int  U elecU vib .U rott (II-32)

II.10. Equilibre thermodynamique local (etl) :

Les lois que nous venons de voir s’appliquent à des milieux dont la température est
uniforme. Il n’existe pas de gradient de température. Les flux sont isotropes et donc
algébriquement nuls en un point. Dans la réalité, les conditions précédentes ne sont jamais
réalisées dans le cas des plasmas d’arc. Il existe des gradients de température et de densité, il
ne peut donc y avoir micro - réversibilité des processus radiatifs. Par conséquent, la
distribution spectrale du rayonnement ne suit pas la loi de Planck.

43
Chapitre II plasma en équilibre thermodynamique

L’hypothèse de l’équilibre thermodynamique complet n’est donc plus valable pour ce

type de plasma. Mais, si les densités électroniques sont supérieures à 10 22 m-3 pour des
températures au dessus de 9000K, les travaux de Griem nous indiquent que les phénomènes
collisionnels électroniques sont prépondérants devant les phénomènes radiatifs, qu’il y a
micro - réversibilité des processus collisionnels et que certaines lois sont toujours valables (la
loi de Maxwell, la loi de Boltzmann, la loi de Saha, la loi de Guldberg - Waage)

Dans ce cas, nous pouvons supposer que l’équilibre thermodynamique est établi
localement en chaque point du milieu. Le plasma est donc considéré à l’équilibre
thermodynamique local (ou ETL).

II.11. Calcul de la composition du plasma :

Le calcul de la composition du plasma à l’équilibre du mélange O 2-N2. Ce calcule

consiste à déterminer la densité des différentes particules présentes. Le plasma est supposé en
équilibre thermodynamique locale avec une variation de température de 1000 à 20000 K, des
pressions et des pourcentages d’oxygène. Dans nos conditions, on ne tient compte que des
espèces suivantes, N, O, N2, O2, N2+, O2- O3 et les électrons. La résolution du système est
réalisée par la méthode de Newton - Raphson. Les fonctions de partitions utilisées dans le
calcul sont tirées de la littérature. Dans ces conditions, les lois d’équilibre utilisées sont la loi
de Maxwell, la loi de Saha et la loi de Guldberg - Waage. Le système doit contenir certaines
équations appelées lois annexes qui sont les suivantes :

II.11.1. La loi de Dalton :

Elle relie la densité des espèces à la pression et à la température

P   ni KT (II-33)
i

La somme porte sur toutes les espèces présentes dans le plasma. Comme exemple,
dans le cas d’un mélange de N2-O2 la loi de Dalton s’écrit sous la forme suivante :

44
Chapitre II plasma en équilibre thermodynamique

II.11.2. La loi de neutralité électrique :

Le plasma est crée à partir d’un gaz neutre, alors il est évident que cette
propriété se conserve, et on aura toujours

n e   z i ni (II-34)
i

ne : la densité électrique.

ni : la densité des ions de charge zi.

Dans notre cas cette loi a la forme suivante :

II.11.3. Loi de conservation de la quantité de la matière :

Elle exprime la conservation de nombre de moles au de la quantité de la matière dans

le mélange. Elle s’applique dans deux cas :

- si le plasma est formé de mélange de gaz, par exemple si P 1 est le pourcentage ou la fraction

molaire de l’azote dan le mélange N2-O2 on peut écrire :

1  P1 N N 2
 P1 N O2 (II-35)

II.12. Conclusion :
Nous avons présenté dans la première partie de ce chapitre les équations qui vont nous
permettre de modéliser le plasma. Dans le prochain chapitre nous exposons la méthode de
résolution de ces équations.

45