Décharge Couronne-Application Et Modélisation

DAILL Guillaume
FOURNAUD Benot
Projet Utilisation Rationnelle de lnergie :

DECHARGE COURONNE,
APPLICATIONS ET MODELISATION

Encadrants: Hamou Sadat
Nicolas Dubus
nergtique Industrielle Anne 2004/2005
Rsum :
Ces dernires annes, la communaut scientifique sintresse de plus en plus aux
plasmas froids. En effet, les plasmas froids sont une technologie prometteuse qui
tend pntrer un nombre croissant de domaines dapplications. Ils peuvent tre
gnr par plusieurs procds. Dans ce projet, seul les plasmas froids gnrs par
dcharges couronnes ont t tudi. Dans un premier temps, les configurations et
les mcanismes de formation et de propagation du steamer sont exposs. Ensuite,
un inventaire non exhaustif des champs dapplication des plasmas froid gnrs par
dcharge couronne est effectu, parmi lesquels les traitements des effluents gazeux
et les traitements de surfaces. Dans une dernire partie, un tat de lart est
rapidement prsent sur les travaux visant modliser les phnomnes se droulant
lors de la cration dun plasma froid.

Abstract:

These last years, the scientific community is interested more and more in cold
plasmas. Indeed, cold plasmas are a promising technology which tends to penetrate
a growing number of applicability. They can be generated by several processes. In
this project, only the cold plasmas generated by corona discharges were studied.
Initially, the configurations and the mechanisms of formation and propagation of the
steamer are exposed. Then, a non exhaustive inventory of the applications fields of
the cold plasmas generated by corona discharge is carried out, among which off-gas
treatments and surface treatments. In a last part, a state of the art is quickly
presented on work aiming at modelling the phenomena proceeding during the
creation of cold plasma.

http://kronos.fbhost.be
Dcharge couronne, applications et modlisation 1
SOMMAIRE
Introduction................................................................................................................. 4
I. Les dcharges couronnes : ................................................................................. 5
A. Les plasmas..................................................................................................... 5
B. Configuration des dcharges couronnes : ....................................................... 6
C. Formation des dcharges dans les gaz........................................................ 6
1. Mcanismes dactivation du gaz et paramtres fondamentaux.................... 6
2. Dcharge de Townsend et critre dauto-entretien de la dcharge.............. 8
3. Loi de Paschen............................................................................................. 9
D. Proprits spcifiques aux dcharges couronnes........................................ 9
1. Caractristique courant-tension et rgimes de dcharge............................. 9
2. Champ lectrique et topographie de la dcharge....................................... 11
E. Les dcharges filamentaires .......................................................................... 12
1. Dcharges pointe positive-plan ou streamers .......................................... 13
2. Dcharges pointe ngative-plan................................................................. 15
3. Dcharges alternatives............................................................................... 15
4. Vent lectrique et consquences hydrodynamiques .................................. 15
II. Application des dcharges couronnes............................................................... 17
A. Traitement des effluents gazeux.................................................................... 17
1. Les prcipitateurs lectrostatiques............................................................. 17
2. Destruction des polluants par dcharges couronne ................................... 18
B. Gnrateurs d'ozone ..................................................................................... 19
1. Gnrateur avec Dcharge Couronne Oxygne pure (CDOX) ............ 19
2. Gnrateur avec dcharge Couronne coefficient rduit (RCCD) ......... 20
C. Traitement des surfaces............................................................................. 20
1. Nettoyage et rosion de surface................................................................. 20
2. Dpt de films ............................................................................................ 21
3. Modification de la chimie de surface .......................................................... 21
4. Conclusion ................................................................................................. 22
D. Autres applications..................................................................................... 22
1. limination de llectricit statique.............................................................. 22
2. Contrle actif dun coulement ................................................................... 23
E. Conclusion..................................................................................................... 24
III. Modlisation................................................................................................... 25
A. Mise en quation............................................................................................ 25
B. Les diffrentes mthodes de rsolutions ....................................................... 26
1. Mthode des lments frontires ............................................................... 26
2. Mthode des lments finis........................................................................ 27
3. Mthode des caractristiques..................................................................... 28
C. Exemples de modlisation ......................................................................... 29
1. Exemple faisant intervenir la BEM et la FEM
[3]
........................................... 29
2. Exemple faisant intervenir la mthode des diffrences finies
[2]
.................. 30
3. Conclusion ................................................................................................. 35
Conclusion................................................................................................................ 36
Bibliographie............................................................................................................. 37
INDEX DES FIGURES
Figure I-1 Systmes dlectrodes couramment utiliss pour la production de
dcharges............................................................................................................ 6
Figure I-2 Tension de claquage pour diffrents gaz en fonction du produit pression,
distance interlectrode. ....................................................................................... 9
Figure I-3 Caractristique courant-tension dune dcharge couronne dans lair
pression atmosphrique (configuration pointe-plan, distance interlectrode d=13
mm). .................................................................................................................. 10
Figure I-4 Lignes isopotentiel avant (a) et 10 Os aprs (b) passage du streamer. .... 11
Figure I-5 Topographie de la dcharge couronne en configuration pointe plan........ 12
Figure I-6 Mcanisme de dveloppement du streamer propos par Townsend. ...... 13
Figure I-7 Mcanisme de propagation du streamer par photoionisation................... 14
Figure I-8 Evolution du champ rduit E/N au passage du streamer. ........................ 14
Figure I-9 Visualisation de lhydrodynamique par introduction de fume (a) et par
simulation (b). .................................................................................................... 16
Figure II-1 Principe de fonctionnement d'un ioniseur................................................ 17
Figure II-2 Les trois fonctions dun catalyseur dNOx standard. .............................. 19
Figure II-3 Eliminateur d'lectricit statique. ............................................................. 22
Figure II-4 Principe de contrle d'un coulement par dcharge couronne................ 24
Figure III-1 Dcharge Couronne grce un dispositif double; (1) masse, (2) lectrode
ionisante, (3) support mtallique. ...................................................................... 31
Figure III-2 Modlisation relle (a) et simplifie (b) de la gomtrie de l'lectrode. .. 31
Figure III-3 Schma suivi lors de la rsolution numrique. ....................................... 33
Glossaire
Pre disruptif : Etat prcdant la disruption.

Disruption : Claquage lectrique.

Steamer : Zone de fort champ et de faible conductivit o se forment les espces
excites et les radicaux.

Enduction : Action denduire un support textile dun produit destin lui confrer des
qualits particulires, en modifier laspect, etc.

C.O.V. : Composs Organiques Volatils.

C.E.A. : Commissariat lnergie Atomique.

B.E.M. : Mthode des lments frontires (Boundary Element Method).

F.E.M. : Mthode des lments finis (Finite Element Method).

F.D.M. : Mthode des diffrences finies (Finite Difference Method).

M.O.C. : Mthode des caractristiques (Method of Caracteristic).

I.A.A. : Industrie Agro Alimentaire.

Gap : Distance inter lectrode.

Introduction

De nos jours, les plasmas sont de plus en plus utiliss dans lindustrie. Il existe deux
catgories de plasma, les plasmas dits thermiques et les plasmas dits froids. Les
dcharges couronnes sont un des procds pouvant conduire la cration de ces
derniers.
Bien que relativement jeune, lutilisation des techniques faisant intervenir les
dcharges couronnes tend prendre de limportance. En effet, elles sont dhors et
dj largement utilise dans les domaines de traitement des effluents gazeux et des
surfaces. Aussi, afin de minimiser les cots de dveloppement, des rcentes
recherches essayent de modliser les phnomnes mis en jeu.

I. Les dcharges couronnes :

A. Les plasmas

Dans la nature, tous les corps sont constitus partir d'atomes (ou d'assemblages
d'atomes, appels molcules). A son tour, un atome est compos :

d'un noyau central, charg positivement, qui est un assemblage de protons et
de neutrons.
d'un nuage priphrique compos d'un cortge d'lectrons, chargs
ngativement, qui tournent des vitesses trs leves autour du noyau.

A l'tat solide, les atomes sont fermement emprisonns dans un rseau rigide
(comme dans la glace par exemple). Lorsque la temprature monte, le solide se
liqufie et passe ltat liquide, o les atomes peuvent glisser les uns par rapport
aux autres, ce qui permet au liquide d'pouser la forme d'un rcipient. En continuant
de chauffer, l'tat gazeux est atteint : les atomes se dplacent alors librement,
indpendamment les uns des autres (l'eau s'est transforme en vapeur). Enfin, de
trs hautes tempratures (plusieurs millions de degrs), les constituants de l'atome
se sparent, noyaux et lectrons se dplacent indpendamment et forment un
mlange globalement neutre : c'est un plasma.

Un plasma est un gaz ionis obtenu en fournissant de lnergie la matire par
adjonction dun champ lectrique ou magntique, ou sous forme de chaleur.

Ce quatrime tat de la matire, qui se retrouve dans les toiles et le milieu
interstellaire, constitue la majorit de notre univers (autour de 99 %). Il nexiste pas
sur Terre l'tat naturel, si ce n'est dans les clairs ou les aurores borales, mais il
est produit artificiellement en appliquant des champs lectriques ou magntiques
suffisamment puissants pour sparer le noyau de ses lectrons dans les gaz.

On distingue 2 types de plasmas :

Les plasmas thermiques pour lesquels les ions, les lectrons et les molcules
ont des tempratures sensiblement identiques, typiquement de lordre de
5 000 K 50 000 K (utilis notamment pour la fusion thermonuclaire).
Les plasmas hors quilibre (ou froid ) pour lesquels la temprature
lectronique (1000 K 10000K) est suprieure la temprature du gaz
ambiant (gnralement infrieure 500 K).

Les plasmas froids peuvent tre forms par des dcharges couronnes. Ces dernires
sont particulirement bien adaptes la formation de plasmas froid car aucun autre
type de dcharges ne permet en effet une localisation et une distribution en nergie
des espces (lectrons, molcules neutres ou ions) aussi bien dfinie.

B. Configuration des dcharges couronnes :

Les dcharges couronnes sont des dcharges lectriques pr-disruptives qui se
dveloppent dans lintervalle gazeux entre deux lectrodes fortement asymtriques.
Parmi les configurations trs varies de systmes dlectrodes utiliss, les plus
courantes sont les systmes pointe/plan, multipointe/plan, couteau/cylindre,
fil/cylindre ou cylindre/cylindre. La distance inter-lectrodes est habituellement
infrieure au centimtre.

Figure I-1 Systmes dlectrodes couramment utiliss pour la production de dcharges.

Du fait de cette dissymtrie, le champ lectrique entre les deux lectrodes est
htrogne. Les dcharges couronne peuvent sappliquer des milieux gazeux
varis des pressions allant de quelques centaines de pascals quelques
atmosphres. En revanche, ces mmes pressions, les plasmas de gaz rares
exempts dimpurets volueront pratiquement instantanment en arc, entranant le
passage au plasma thermique. En outre, les dcharges couronnes ont pour
avantage dtre faciles produire dans de bonnes conditions de stabilit et de
reproductibilit.
Suivant les applications envisages, la haute tension, qui est gnralement
applique llectrode active, peut tre pulse, alternative ou continue de polarit
positive ou ngative. Lutilisation dune alimentation pulse permet de gnrer de
forts courants pendant des temps trs et donc dobtenir, sans passage larc, des
lectrons de haute nergie, adapts au traitement de certains polluants. Les tensions
de travail habituellement utilises varient de quelques kilovolts une trentaine de
kilovolts environ tandis que les puissances dpassent rarement quelques kilowatts.

C. Formation des dcharges dans les gaz
1. Mcanismes dactivation du gaz et paramtres
fondamentaux

La formation dun plasma est due un transfert dnergie cintique par collision entre
un lectron acclr par un champ lectrique et les molcules neutres du gaz. Les
collisions sont de deux types :
les collisions lastiques
les collisions inlastiques.
Les premires entranent le transfert dune certaine quantit de mouvement des
lectrons vers les molcules neutres qui se trouvent alors acclres. Ces collisions
ne modifient donc pas lnergie interne des molcules en question, mais sont
responsables de phnomnes thermiques (conduction et diffusion de chaleur au sein
du gaz) et hydrodynamiques.
Les collisions inlastiques sont quant elles responsables de la fragmentation des
molcules impactes ou de la modification de la distribution de lnergie (modification
de ltat vibrationnel, translationnel ou lectronique). Si son nergie cintique est
suffisante, un lectron peut, lors dune collision inlastique, ioniser ou exciter la
molcule neutre rencontre.

Pour obtenir une description plus formelle de ces processus, il est intressant de
dterminer le nombre de chocs que subira llectron par unit de temps. Nous
supposerons que la distribution de vitesse des molcules du gaz obit la loi de
Maxwell-Boltzmann :

En admettant que llectron se dplace une vitesse moyenne v pendant un temps
Yt, il parcourt alors la distance v.Yt et balaie un domaine despace de volume Z.v.Yt,
Z reprsentant la section efficace de collision entre llectron et une molcule. Le
nombre de collision est donc gal :
Z.v.Yt.N (1)
avec N = densit du gaz en nombre de molcules par unit de volume.
Le nombre de collisions par unit de temps se dduit alors de lexpression
prcdente :
n = Z.v.N (2)

La notion de libre parcours moyen, qui reprsente la distance moyenne que peut
parcourir un lectron entre deux chocs successifs, dcoule naturellement de la
relation prcdente :

(3)

Cette notion est intressante car si llectron acclr par le champ lectrique
acquiert une nergie cintique suffisante, il pourra alors ioniser ou exciter la molcule
neutre rencontre. La ractivit du milieu plasmagne dpendra donc des
paramtres suivants :

Composition du gaz introduit : la section efficace Z varie en fonction de la
nature du gaz et de son niveau dexcitation lectronique (pour des molcules
dans leur tat fondamental, nous avons par exemple Z(N
2
)=0,43; Z(CO
2
)=
0,52 et Z(C
6
H
6
)= 0,88)
Pression et temprature du gaz introduit
Tension applique (responsable de lintensit du champ lectrique)
Gomtrie des lectrodes (responsable de la forme des lignes de champ).

2. Dcharge de Townsend et critre dauto-entretien de la
dcharge

Comme nous venons de le voir, lnergie cintique dun lectron dtermine sa
capacit exciter ou dissocier les molcules du gaz plasmagne. Lnergie
acquise par un lectron plac entre deux lectrodes dpend de son libre parcours
moyen et du champ lectrique auquel il est soumis. Lobtention dun nombre
dlectrons suffisant pour exciter le gaz de faon significative exige un mcanisme de
multiplication. Si lnergie dun lectron est suprieure lnergie dionisation dune
molcule M, alors un choc inlastique est susceptible de librer un lectron, ce que
traduit lexpression suivante : M + e
-
]M+ + 2e
-
Lefficacit dionisation est estime 10 % pour un lectron possdant une nergie
trois fois suprieure au potentiel dionisation de lespce considre (par exemple,
les nergies dionisations de certains gaz couramment utiliss sont : 15,6 eV pour
N2, 12,2 eV pour O2, 12,6 eV pour H2O ou 13,7 eV pour CO2).
Les deux lectrons obtenus par ionisation sont alors acclrs par le champ
lectrique et peuvent leur tour ioniser de nouvelles molcules.

Townsend a propos en 1914 lintroduction dun coefficient dionisation ^ (appel
premier coefficient de Townsend) pour traduire mathmatiquement ce phnomne
dionisation. Ce coefficient dpend du champ rduit E/N (avec N dsignant la densit
du gaz) et son inverse 1/ ^ reprsente le libre parcours moyen dun lectron entre
deux collisions ionisantes. Lapparition dun lectron primaire au niveau de la
cathode se traduit par larrive de exp(^d) lectrons au niveau de lanode. Ce
processus est appel avalanche lectronique. La relation de Townsend scrit
alors :
I = I
0
. exp(^d) (4)
Avec I : courant (A)
I
0
: courant initial d aux lectrons primaires (A)
d : distance interlectrode (m)
^ : premier coefficient de Townsend (m
-1
)
En dehors de ce mcanisme de multiplication des lectrons propos par Townsend,
des lectrons secondaires peuvent tre mis par impact dions positifs M
+
sur la
cathode ou par impact de photons (mis par dsexcitation radiative aprs impact
lectronique : M + e
-
] M* + e
-
] M + e
-
+ h`) sur la cathode ou les molcules
proches. La probabilit de produire ainsi des lectrons secondaires se traduit par
lintroduction dun second coefficient de Townsend : a.
Le critre dauto-entretien de la dcharge (nombre adimensionnel) scrit alors :
a.(exp(^d)-1) = 1 (5)

Comme ce mcanisme est consommateur dlectrons, il est comptitif du mcanisme
dionisation. Pour prendre en compte ce phnomne, il convient alors dintroduire
dans la relation de Townsend une probabilit dattachement lectronique b :
I = I
0
.exp( ^- b)d (6)
Le critre dauto-entretien de la dcharge scrit donc :
a.( exp( ^- b)d-1) = 1 (7)

3. Loi de Paschen

La tension dalimentation minimale pour dclencher la dcharge est appele tension
de claquage. Il existe une relation qui lie la valeur de cette tension de claquage la
pression du gaz et de la distance interlectrodes : la loi de Paschen.
Cette loi est habituellement reprsente par la courbe caractristique Uc = f(p.d), o
p est la pression et d la distance interlectrode, et dpend de la nature du gaz
introduit.
Figure I-2 Tension de claquage pour diffrents gaz en fonction du produit pression, distance
interlectrode.

Nous pouvons constater que pour un gaz pression atmosphrique, les courbes
passent par un minimum pour une distance interlectrodes denviron 10 Om (p.d =
0,5 5 environ). Cela signifie que, pour la production de dcharges pression
atmosphrique dans des conditions dutilisation usuelles (d c 1 mm), la tension de
claquage est une fonction croissante de la distance interlectrode.

D. Proprits spcifiques aux dcharges couronnes
1. Caractristique courant-tension et rgimes de dcharge

La caractrisation courant-tension de la dcharge couronne est obtenue en faisant
crotre la tension applique lintervalle gazeux. Diffrents rgimes de dcharge sont
alors parcourus (Figure 3).

Figure I-3 Caractristique courant-tension dune dcharge couronne dans lair pression
atmosphrique (configuration pointe-plan, distance interlectrode d=13 mm).

Le premier rgime (I) consiste en la collecte des espces charges prsentes
naturellement dans le milieu. Ces espces sont produites par limpact de rayons
cosmiques ou de particules issues de la dsintgration spontane de certains
atomes (radioactivit naturelle) sur le gaz prsent dans lespace interlectrodes.

Le second rgime (II), qui se produit partir dune tension seuil V
0
, traduit le
dclenchement de lavalanche lectronique de Townsend mais ne rpond pas
encore au critre dauto-entretien de la dcharge. La dcharge est qualifie de non-
autonome car elle dpend de processus dionisation externes pour donner naissance
aux lectrons germes. Pour ce rgime, le courant augmente trs fortement pour une
trs faible variation de la tension applique.

Pour le troisime rgime (III), lionisation du gaz par la dcharge est suffisante pour
ne plus ncessiter de phnomnes dionisation extrieurs : la dcharge est auto-
entretenue. Cest le rgime de la dcharge couronne, o apparaissent les dards
(aussi appels streamers) responsables de la composante impulsionnelle du
courant.

Au-del dune tension applique Vr, la temprature du canal augmente fortement, ce
qui provoque la dilatation brutale du gaz. Le champ rduit E/N est alors suffisamment
important pour que les phnomnes dionisation deviennent majoritaires par rapport
lattachement lectronique. Le rgime darc (IV) est atteint, le milieu passe de ltat
de plasma hors quilibre ltat de plasma thermique.

2. Champ lectrique et topographie de la dcharge

La dissymtrie des lectrodes entrane la non homognit du champ lectrique.
Parmi les configurations les plus tudies, nous retrouvons les configurations pointe-
plan et coaxiale (cylindre-cylindre ou fil-cylindre). La distribution spatiale du champ
lectrique est donne par les deux relations suivantes :

En gomtrie pointe hyperbolique-plan :

(8)

avec
E(x) : champ lectrique (V.m
-1
),
U : tension applique la pointe (V),
x : distance la pointe sur laxe pointe-plan (m),
r : courbure de la pointe son extrmit (m),
d : distance interlectrode (m).

En gomtrie coaxiale :
(9)
avec
E(r) : champ lectrique (V.m
-1
),
U : tension applique llectrode centrale (V),
r : distance radiale du point considr (m),
R
e
: rayon de llectrode externe (m),
R
i
: rayon de llectrode interne (m).

Des simulations ont t effectues en configuration pointe-plan pour dterminer le
champ lectrique dans lespace interlectrode avant ou aprs passage du streamer.

Figure I-4 Lignes isopotentiel avant (a) et 10 ,s aprs (b) passage du streamer.

Dans le cas dune gomtrie pointe-plan, le champ intense qui rgne proximit de
llectrode active (10
-19
< E/N < 10
-17
V.m
2
) permet lionisation du gaz. Les ions
prsents dans cette zone peuvent tre positifs ou ngatifs. Au-del de cette rgion
se trouve une zone de drive ionique o le champ lectrique est trop faible pour
permettre les ractions dionisation (10
-21
< E/N < 10
-19
V.m
2
). Des ions ngatifs sont
alors forms par attachement lectronique. Ils sont ensuite entrans le long des
lignes de champ et ragissent avec les molcules rencontres. La zone de drive
ionique, lectriquement passive, est le sige dune intense activit physico-chimique
(dispersion des espces par le vent lectrique, ractions radicalaires).
A linterface entre la zone dionisation et la zone de drive ionique, la consommation
dlectrons par attachement lectronique gale leur libration par impact
lectronique : les coefficients ^ et b sont alors gaux.

Figure I-5 Topographie de la dcharge couronne en configuration pointe plan.
E. Les dcharges filamentaires

Les phnomnes lmentaires de dveloppement et de propagation des dcharges
ont surtout t dcrits pour des configurations pointe-plan. Cette gomtrie est en
effet la plus simple car le systme est mono-impulsionnel et permet une parfaite
reproductibilit des mesures.
Deux cas principaux sont distinguer :
lapplication dune tension positive la pointe,
lapplication dune tension ngative la pointe.
Dans le cas dune alimentation lectrique alternative, les mcanismes lmentaires
sont alternativement les mmes que pour les dcharges positives ou ngatives,
condition que les charges despace aient pu se diluer dans le milieu avant
linversion de la polarit.
Nous pourrons enfin admettre que le mcanisme gnral de propagation du dard et
les ordres de grandeur des diffrents phnomnes observs sont similaires avec
dautres gomtries (fil-cylindre, fil-plan), mme si des nuances peuvent apparatre,
comme par exemple la localisation du point dancrage de la dcharge ou lorientation
du flux gazeux.

1. Dcharges pointe positive-plan ou streamers
Townsend a propos en 1915, un modle de dcharge qui sapplique aux systmes
pointe positive-plan. Dans ce modle, un lectron prsent au voisinage de la pointe
est acclr vers celle-ci par le champ lectrique. Il subit lors de son trajet des
collisions ionisantes, produisant une avalanche dlectrons attirs eux aussi vers la
pointe. Aprs leur collecte par lanode, ils laissent dans leur sillage des cations qui
forment une zone de charges despace positive (Figure 6). Celle ci modifie
localement le champ lectrique, prolongeant celui cr par lanode dans lespace
interlectrodes.
Figure I-6 Mcanisme de dveloppement du streamer propos par Townsend.

Le mcanisme de propagation du streamer a t propos par Loeb

et Meek. Cette
thorie du streamer introduit la notion de production dlectrons en amont de la
charge despace. Des photons issus de la dsexcitation des molcules vont tre
absorbs par des molcules neutres qui se sparent alors en un cation et un
lectron : cest le phnomne de photoionisation.
La photoionisation a lieu suffisamment proche de la charge despace pour bnficier
dun champ lectrique local important. Ces lectrons sont leur tour acclrs par le
champ lectrique et engendrent des avalanches secondaires provoquant la
propagation du streamer (Figure 7). Aprs passage du streamer, il se forme entre
celui-ci et lanode un plasma neutre englob par une zone de charge positive. Le
champ de ce plasma est trs faible (de lordre de 5 kV.cm
-1
), favorisant alors les
phnomnes dattachement lectronique par les espces neutres. La chute du
champ lectrique dans le plasma rend les vitesses de diffusion plus importantes que
les vitesses de convection lies directement au champ lectrique.
Figure I-7 Mcanisme de propagation du streamer par photoionisation.

Le dard, dont le diamtre mesure environ 20 Om, se propage une vitesse de 10
5
10
6
m.s
-1
. La densit lectronique en tte de streamer atteint 10
15
-10
16
lectrons par
cm
3
. Lnergie moyenne des lectrons y est denviron 15 eV contre 2 4 eV.
La temprature est de lordre de 300 350 K. Lamplitude de limpulsion de courant
due larrive des lectrons sur le plan peut atteindre quelques centaines de
milliampres.
Les mcanismes de propagation du streamer en gomtrie pointe-plan ont fait lobjet
de nombreuses tudes et modlisations permettant de dfinir linfluence de diffrents
paramtres comme la pression de travail, la composition du mlange gazeux ou la
gomtrie du racteur. La srie dimages sur la figure I-8 illustre ainsi lvolution du
champ rduit E/N au passage du streamer dans lespace interlectrode.

Figure I-8 Evolution du champ rduit E/N au passage du streamer.
2. Dcharges pointe ngative-plan

Dans le cas dune pointe alimente par une tension ngative (pointe cathodique),
lapparition dun lectron germe proximit de la pointe peut provoquer une
avalanche lectronique qui se dveloppe en champ divergent. La cathode acclre
alors les ions positifs qui la percutent et lui arrachent des lectrons secondaires. Ce
phnomne, propre lutilisation dune pointe cathodique, entrane un abaissement
de la tension de claquage de la dcharge par rapport la configuration pointe
positive-plan. La charge despace positive domine dans la zone dionisation tandis
que la zone de drive ionique possde une charge despace ngative, faible ou forte
suivant limportance des processus dattachement lectronique.

3. Dcharges alternatives
En appliquant un champ alternatif au gaz, alors le comportement du systme dpend
de la distance interlectrode et de la frquence de lalimentation lectrique. En effet,
si le temps de parcours du gap par le gaz est infrieur au temps dinversion de la
polarit, alors les charges despace cres par une dcharge auront quitt lespace
interlectrode avant lapparition de la dcharge suivante. Le comportement
lmentaire du systme consistera donc en une alternance de dcharges positives
puis ngatives indpendantes.
En augmenttant la frquence du signal lectrique, une charge despace rsiduelle
persistera dans lespace interlectrode, favorisant le dclenchement de la dcharge
suivante. Ce phnomne peut se traduire par le maintien de la dcharge avec une
tension dalimentation lgrement infrieure la tension de claquage de la dcharge.
Cette accumulation de charges entrane aussi une baisse de la tension de passage
en rgime darc.

4. Vent lectrique et consquences hydrodynamiques

Sous laction de la dcharge, de nombreuses espces charges (ions positifs,
ngatifs ou lectrons) sont cres dans lespace interlectrode. Ces espces sont
acclres par le champ lectrique et entranes vers la cathode pour les ions
positifs, et vers lanode pour les ions ngatifs et les lectrons. En chemin, elles
percutent dautres molcules auxquelles elles transmettent une partie de leur
quantit de mouvement.
Lacclration du gaz dans lespace interlectrode, sous leffet de la dcharge, cre
un vritable courant appell le vent lectrique. Il est responsable de la dispersion et
du mlange des molcules neutres, excites ou ionises.
Seule une faible partie de lnergie consomme par le systme est utilise pour la
mise en mouvement des molcules du gaz, les mcanisme de chauffage et
dexcitation reprsentant la majorit du transfert dnergie. Le caractre impulsionnel
de la dcharge agit sur lcoulement comme une srie dondes de pression dont la
frquence est identique celle des pulses. La vitesse du vent lectrique dpend de
la distance interlectrode et du courant de dcharge.
L encore, le phnomne a surtout t tudi dans le cas de configurations pointe-
plan : dans une gomtrie coaxiale, il semble difficile de raliser des mesures ou de
visualiser la dcharge car le systme est multi-impulsionnel (la dcharge ne se
produit pas toujours au mme point) et le corps du racteur (lectrode ou dilectrique
entour dune lectrode) est gnralement opaque.
De nombreuses mthodes ont t utilises pour visualiser les effets du vent
lectrique sur lhydrodynamique au sein des systmes dcharge couronne telles
que lintroduction de particules ou la simulation.

Figure I-9 Visualisation de lhydrodynamique par introduction de fume (a) et par simulation (b).

Les images de la figure I-9 illustrent parfaitement les phnomnes de recirculation,
de transfert de matire et de mlangeage qui ont lieu dans les racteurs dcharges
couronne.
II. Application des dcharges couronnes
A. Traitement des effluents gazeux

1. Les prcipitateurs lectrostatiques

Les premires applications des plasma froids ont t les prcipitateurs
lectrostatiques, qui utilisent des dcharges couronnes appliques au traitement des
effluents gazeux contenant des poussires. Nous pouvons les trouver la sortie des
gaz des centrales thermolectriques, des cimenteries, etc. Ils sont utiliss galement
pour amliorer la qualit de lair intrieur dun local (ioniseur dair : OR'IONS,
DAVIDSON).
Des tudes ont montres que la prsence dun grand nombre dions ngatifs
entrane une faible pollution. Un ioniseur dair pige en fait la micro pollution
(particule).

Figure II-1 Principe de fonctionnement d'un ioniseur.

Lioniseur met des electrons. Ces lectrons vont venir sattacher aux gros ions
positifs ( la pollution ) pour former des particules neutres. Ces particules vont
acqurir des charges ngatives en quantits suffisantes pour devenir des gros ions
ngatifs. Ces gros ions ngatifs vont ensuite se repousser mutuellement et le seul
endroit possible o ils peuvent aller (attraction lectrique) est la surface conductrice
(donc de dcharge) que reprsentent le sols et les murs : ils ne seront donc plus en
suspension dans l'air.
2. Destruction des polluants par dcharges couronne

Depuis une bonne dizaine d'annes, les milieux scientifiques d'abord, industriels
ensuite, s'intressent la possibilit d'utiliser les plasmas froids pour dtruire les
polluants et non plus les piger. L'ide de base est d'exploiter la ractivit chimique
des radicaux et espces engendrs au sein du plasma pour transformer les
molcules polluantes en molcules inoffensives, la raction se droulant dans l'air
la pression ordinaire et sans pertes thermiques. Les rsidus obtenus peuvent en
outre tre parfois valorisables. Les polluants viss dans les dveloppements actuels
sont essentiellement les oxydes de soufre et d'azote, les COV et les molcules
malodorantes.

En prenant un volume actif correspondant un volume de plasma, un facteur
volumique est dfinissable : le quotient dun volume de plasma par un volume de gaz
dans le racteur. Ce facteur est toujours plus petit que lunit, alors que des
conversions de polluant comprises entre 80 et 100 % sont aisment observes Cela
signifie que le traitement du polluant peut seffectuer hors de la zone ionise.
Sachant que la prsence du plasma est ncessaire, cela montre que la chimie est
initie par des espces produites dans le plasma, dont certaines sont actives hors de
la zone ionise. Ces constatations illustrent lintrt des dcharges lectriques qui
permettent de traiter un volume important alors que les apports dnergie restent trs
localiss. Les processus dexcitation et de dsexcitation par collisions lectroniques
trs rapides sont suivis de ractions inities par des radicaux libres. Certains
radicaux vont diffuser ou se mlanger par convection avec les autres espces du gaz
traiter avec une chelle de temps de lordre de la milliseconde.

Le traitement des effluents par plasma froid est dj entr dans son ge industriel.
C'est ainsi que la socit Paganetti, Moosch en Alsace, a dvelopp avec le
soutien du rseau Novlect son procd "Electronflux", essentiellement destin au
traitement des nuisances olfactives. Ce procd s'applique en milieu ferm, c'est
dire soit l'extrieur du local, en canalisant l'effluent traiter, soit l'intrieur, en
utilisant le local comme chambre de raction, ce qui permet d'amliorer la qualit de
l'air intrieur et du mme coup les conditions de travail. La premire ralisation a eu
lieu au Havre, la station de relevage des eaux uses et pluviales "Augustin
Normand" o l'quipement mis en place (1996) traite un dbit d'air malodorant de 7
500 m
3
/h. Dans la mme ville, Paganetti a install une autre centrale sur un bassin
orage d'une contenance de 50 000 m
3
o le dbit trait atteint 12 500 m
3
/h. Plus
rcemment, le SAN d'Evry (91) a pris la mme dcision pour viter les nuisances
olfactives manant de son unit de traitement des produits de curage qui spare les
solides selon leur nature (sable, dbris vgtaux, encombrants) de la phase liquide
charge en matires organiques. L'quipement plasma froid traite 8 500 m3/h
extraits du local de 1200 m3. Sa puissance lectrique totale est de 45 kW, dont 200
W seulement pour les gnrateurs de plasma. D'aprs la socit Paganetti, elle
consomme en moyenne 11 kW par heure et sa maintenance cote 10 k par an, y
compris les consommables.

Cependant les applications des dcharges pour lenvironnement sont encore un
domaine de recherche trs ouvert et ce, pour trois grandes raisons :
la complexit des phnomnes impliqus qui combine effets chimiques,
physiques et lectriques,
les aspects lgislatifs sur la pollution qui ne sont pas encore suffisamment
coercitifs pour inciter les industriels investir dans de nouvelles technologies,
le cot nergtique associ au procd.

De plus, des travaux rcents ont montr que les techniques de plasmas froids,
associes la catalyse, pouvaient tre performantes et conduire la rduction des
cots nergtiques. Par exemple, dans les traitements DNO
x
qui consistent
transformer les oxydes dazote en azote et oxygne, les deux premires fonctions du
catalyseur standard (figure 2) peuvent tre remplaces par le plasma. Il est ainsi
possible datteindre un taux de conversion de NO en N
2
pouvant dpasser 60 %. La
perspective daccrotre lactivit catalytique par un plasma temprature ambiante
est dont particulirement intressante, car cela permet au catalyseur de fonctionner
ds le dmarrage du moteur.

Figure II-2 Les trois fonctions dun catalyseur dNOx standard.

B. Gnrateurs d'ozone

Une autre solution pour traiter les effluents gazeux est dutiliser lozone (O
3
) : au lieu
de crer un milieu fortement ractif, une espce possdant un fort pouvoir oxydant
est introduite dans le gaz traiter. Cette espce est lozone.
Lozone est galement utilis dans le traitement de leau mais aussi dans de
nombreux domaine comme lindustrie agro-alimentaire en tant que dsinfectant ou
conservateur.
Il existe trois gnrateurs pour la production artificielle d'ozone et ils emploient des
mthodes diffrentes de production:
Gnrateur Ultraviolets (UV),
Gnrateur avec Dcharge Couronne Oxygne pure (CDOX),
Gnrateur avec Dcharge Couronne coefficient rduit (RCCD).
Le premier rplique le mcanisme de cration dozone par le soleil ; les deux autres
les phnomnes lectriques orageux.
1. Gnrateur avec Dcharge Couronne Oxygne pure
(CDOX)

Il produit lozone travers le passage d'oxygne dshumidifi, par l'intermdiaire
d'un champ lectrique trs haute intensit. Ces appareils produisent des grandes
quantits d'ozone, mais ils ncessitent toujours et seulement de l'oxygne l'tat
sec. Lavantage est quen employant seulement de loxygne, ils n'engendrent pas
de sous-produits dangereux. Ces gnrateurs sont jusqu' 10 fois plus puissants
quun gnrateur UV. Ils produisent cependant une considrable quantit de chaleur
(deux fois plus quun gnrateur UV) qui doit tre dissipe. Le cot de dissipation
doit toujours tre prit en compte dans le cot global de linstallation.
Les gnrateurs de ce type ne sont employs que dans le domaine des traitements
industriels de haut niveau, et dans les villes pour rendre potable des millions de
mtres cube d'eau l'heure.

2. Gnrateur avec dcharge Couronne coefficient rduit
(RCCD)

Il produit lozone travers le passage d'air, par l'intermdiaire un champ lectrique
trs haute intensit. Ces types d'appareils n'exigent que de l'air et ils sont tout fait
immuniss contre toute sous-production de gaz nitrognes dangereux. La capacit
de gnration est jusqu' 4 fois suprieure aux appareils UV de mme puissance
lectrique.
Il peut tre employ par exemple, dans l'aviculture o la production demande est
trs basse : Un gnrateur de type RRCD qui produit 8 milligrammes d'ozone l'heure,
suffit pour dsinfecter: couveuses, closoirs, chambres chaudes, box avec nids,
salles et chambres grandes mme jusqu' 12 mtres cube.

C. Traitement des surfaces

Le traitement des surfaces est le deuxime domaine dutilisation des plasmas froids
produits par dcharge couronne. Le premier secteur les avoir adopts couramment
est la micro-lectronique pour la fabrication des circuits intgrs et des circuits
imprims. En raison de leurs caractristiques et de leurs qualits intrinsques, ils se
diffusent dans dautres secteurs comme la mtallurgie et la mcanique (durcissement
superficiel et protection anti-corrosion des pices et outillages), la plasturgie
(prparation de la surface une enduction ultrieure, nettoyage et dgraissage,
laboration de fonctions slectives sur membranes, modification du coefficient de
frottement, cration dun effet barrire, notamment pour les emballages
alimentaires,), le textile (anti-mouillabilit, anti-glisse), ou encore loptique et lopto-
lectronique (dpts organo-mtalliques multicouches).

1. Nettoyage et rosion de surface

Cela comprend les oprations de nettoyage, dcapage, de prparation de surface
avant dpt et de gravure. En ce qui concerne lenvironnement et la mise en uvre
de ces procds, les techniques classiques sont souvent surpasses par lutilisation
de dcharges couronnes. Cela peut aussi tre plus intressant conomiquement.

Pour nettoyer une pice couverte d'huile et de particules, au lieu d'utiliser un solvant,
coteux et dangereux, il faut commencer par chasser le plus gros mcaniquement,
avec de l'air comprim ou de l'eau, puis liminer la fine couche d'huile rsiduelle
l'aide d'un plasma froid l'oxygne qui va volatiliser les polluants organiques de
surface par un effet de "combustion froide" contrle. Rsultat, de la matire est
conomise et l'environnement protg (consommation et manation de solvants)
tout en traitant la pice avec dlicatesse. Ce procd est en outre capable d'enlever
mme des substances chimiquement stables ou trs adhsives, et ce sur peu prs
n'importe quelle gomtrie de pice.

Enfin, il est possible de faire suivre le nettoyage d'une opration d'activation de la
surface sans renchrir le procd. Cette activation est souvent ncessaire en
plasturgie pour permettre l'adhrence des peintures, vernis ou colles sur les
matriaux polymres.

2. Dpt de films

Les dcharges couronnes peuvent tre utilises pour dposer des films minces ou
multicouches.

La polymrisation plasmagne est un procd de ce type. Des monomres ltat
gazeux, excits par un plasma, se dposent en couche fortement rticule sur un
support. Il faut bien sr que le gaz utilis contienne des atomes capables de former
des chanes, comme le carbone, le silicium ou le soufre, et sa composition initiale va
tre partiellement modifie au cours du processus. La structure obtenue, trs
dsordonne, et le degr de rticulation sont influencs par les paramtres de
fonctionnement (pression, quantit de gaz, nergie lectrique). Il est ainsi possible
de crer des couches en gradient , dans lesquelles le degr de rticulation varie
en fonction de lpaisseur de la couche. Les revtements obtenus prsentent une
trs forte adhsion sur presque tous les substrats et sont stables chimiquement,
thermiquement et mcaniquement. Ce procd est utilis pour dposer des
revtements anti-abrasion, anti-corrosion, anti-salissures, anti-adhrents et effet
barrire.

3. Modification de la chimie de surface

La modification de la chimie de surface peut tre la fonctionnalisation, le greffage, le
dopage par implantation, diffusion ou strilisation. Cest laide de plasmas froids
que sont ralises un grand nombre de membranes polymres

La modification des fibres cellulosiques par dcharges couronnes a connu plusieurs
succs dans le cadre dtudes sur les composites bois-polymre (proprits
mcaniques, rduction de la viscosit au point de fusion). Le traitement Corona est
une technique dactivation de surface par oxydation. Cela se traduit par une
augmentation de la concentration des groupements aldhyde lors de lactivation de
la surface du bois. La varit de modifications possibles avec les plasmas provient
du choix du gaz impliqu dans le traitement. Loxygne entranera la formation de
groupes plutt acides (accepteurs dlectrons), alors que lazote ou lammoniaque
entranera la formation de groupes plutt basiques (donneurs dlectrons). Ce type
de solution par dcharge lectrique ncessite toutefois une technologie assez
complexe.
4. Conclusion

Les plasmas au travers de ses diffrentes techniques de prparation de surface, de
modification, de polymrisation ou de greffage induit permettent denvisager
llaboration de membranes aux proprits slectives ou dlaborer un matriau
ayant en surface des proprits barrire. La simplicit de mise en uvre, le respect
de lenvironnement (technique sche, sans solvant), le cot de revient relativement
bas au regard de la chimie classique dveloppe pour la synthse de nouveaux
polymres sont aussi des facteurs favorables pour une application industrielle.

D. Autres applications

1. limination de llectricit statique

Llectricit statique est, dans lindustrie, la source de dysfonctionnements gnants.
Elle provoque des dcharges nfastes aux dispositifs lectroniques ou des collages
de poussire sur des produits souhaits purs.
La socit Valitec a dvelopp un systme utilisant l'injection supersonique d'ions
afin dliminer llectricit statique dans de gros dbits de production.
Pour viter les inconvnients dune production dions par dcharge lectrique dans
un milieu susceptible dexploser, les ions sont produits par effet couronne lintrieur
dune buse, constitue dun tube aliment en gaz pur sous pression et ouvert sur la
veine traiter travers le col dune tuyre supersonique.

Figure II-3 Eliminateur d'lectricit statique.
Lcoulement a peu dactions sur les ions, que le champ lectrique dirige vers les
parois du tube. Mais il suffit que le gaz contienne un trs faible pourcentage
dhumidit (point de rose de - 20 C, soit beaucoup plus sec que le produit traiter)
pour que chaque ion se recouvre dune pellicule de glace, sil est cr au voisinage
de la zone dexpansion supersonique du jet caractrise par un important
refroidissement. Larosol ainsi form est expuls par le jet dans la veine. En moins
dune milliseconde, il se sublime alors et lion est restitu. Cet artifice est donc
essentiellement utilis pour favoriser le transfert des charges travers lorifice de la
tuyre, ce qui permet de sparer la rgion de formation des ions de leur rgion
dutilisation.
Alors quun jet subsonique en gaz sec ne permettrait de transfrer dans la veine que
quelques pour mille du courant dions de la dcharge couronne, le jet supersonique
en gaz pur point de rose de - 21 C permet linjection de prs de 80% des ions
produits dans le tube. Il faut noter que lhumidit contenue dans le gaz est minime :
l'injection opre contribue scher et non pas humidifier le produit.
la suite dune srie de tests raliss par le LCIE, le matriel Valitec a obtenu en
2003 un agrment europen autorisant son emploi en atmosphre explosive. Il a
galement t test par le CEA pour vrification de ses missions NO
x
, qui se sont
avres ngligeables. Une troisime batterie de tests au LCIE a permis de vrifier sa
porte : il neutralise une cible une distance de trois mtres.

2. Contrle actif dun coulement

Contrler un coulement consiste modifier ses caractristiques pour lamener dans
un tat souhait. Cela peut permettre une diminution de trane, une augmentation
de la portance des profils, une augmentation du mlange ou une rduction des
bruits.
Limpact de telles recherches est immense pour lindustrie. Les enjeux conomiques
du contrle de lcoulement sont multiples dans de nombreux systmes industriels,
quils fassent intervenir des processus dchange de chaleur, despces chimiques,
de stabilisation de structure ou bien quil sagisse damliorer le mlange pour des
problmes de combustion interne ou de contrle acoustique. Par exemple, rduire la
trane de 10% dun avion civil entranerait une rduction de la consommation en
carburant de lordre de 50 millions de litres durant sa dure de vie.
De nouveaux actionneurs sont donc ltude dont lactionneur
lectroarodynamique . Il utilise le phnomne qui consiste induire un
coulement de gaz partir dune injonction de charge. La mthode choisie pour
injecter ces charges consiste crer un plasma froid laide dune dcharge de
surface de type couronne.
Lactionneur lectroarodynamique est donc constitu de deux lectrodes
places sur une surface isolante entre lesquelles est applique une diffrence de
potentiel. Des ions positifs produit lanode et attirs en direction de la cathode par
les forces de Coulomb, entrent en collision avec les particules neutres du gaz, et
crent ainsi un coulement secondaire (le vent ionique). Cet coulement augmente
ou diminue ainsi la vitesse de lcoulement le long de lobstacle et modifie ses
caractristiques.
Figure II-4 Principe de contrle d'un coulement par dcharge couronne.

Cette mthode active de contrle des coulements convertit directement lnergie
lectrique en nergie cintique, sans apport de masse. Lnorme avantage de ce
dispositif est labsence de partie mcanique mobile souvent fragile. De plus, son
contrle est totalement lectrique ce qui permet des temps de raction trs court et
une modulation en temps et en amplitude trs simple.
Cette technique et lactionneur qui en dcoule font encore lobjet de recherche mais
laissent entrevoir de nombreuses possibilits.

E. Conclusion

Les dcharges couronne sont utilises dans de nombreux domaines o elles
remplacent des techniques souvent plus polluantes (utilisation de solvants).
Cependant elles restent pour certaines applications encore trop coteuses dun point
de vue nergtique compares dautre mthode mieux matrises. Mais Les
recherches effectues sur le sujet tendent accrotre sa pntration dans le milieu
industriel et ouvrent de nouvelles possibilits. La dcharge couronne est donc une
technologie davenir quil est ncessaire de mieux comprendre et matris.

III. Modlisation
A. Mise en quation

La majeure partie des modlisations introduit les quations de Poisson, ainsi,
Armand Caron et Lucian Dascalescu
[2]
supposent :
Llectrode ionisante est un fil sur lequel la dcharge est uniformment
rpartie
Pour des tensions suprieures la tension de claquage, le champ lectrique
sur la surface de llectrode est donn par la loi de Peek
Lpaisseur de la couche dionisation est ngligeable
La zone de drive contient seulement des ions de mme signe
La diffusion thermique et la recombinaison des ions sont ngligeables
La mobilit des ions est constante

Le champ lectrique E
r
drive du potentiel lectrique scalaire V :

(10)

et satisfait lquation de Gauss :

(11)

o q est la densit de charge volumique et l la permittivit du milieu.

La combinaison des quations (10) et (11) conduisent lquation de poisson :

(12)

La densit de charge q doit satisfaire lquation de conservation :

(13)

O J
r
est le vecteur densit de courant. En utilisant lhypothse dune mobilit
ionique des ions constante (b), J
r
peut sexprimer ainsi :

(14)

La majorit des chercheurs rsolvent (12) et (13) sparment, excutant la correction
mutuelle des solutions dans un processus itratif, jusqu ce que la convergence soit
atteinte. La symtrie du systme d'lectrode les a aids trouver des solutions au
problme de champ. Malheureusement, les algorithmes dvelopps pour ces
applications ne sont pas utilisables pour des gomtries plus complexes.

V grad - E =
r
l
q
E div =
r
l
q
nV =
0 J div
Yt
Yq
= +
r
E qb J
r r
=
Certains auteurs ont essay de coupler les quations (12) (13) et (14), et de
rsoudre ainsi :

(15)

Cette approche mne parfois un comportement imprvisible du systme
numrique, cest pourquoi, lquation suivante est prfre:

(16)

B. Les diffrentes mthodes de rsolutions

Il existe plusieurs mthodes possibles pour rsoudre numriquement ces quations.
Parmi celles qui reviennent le plus souvent dans la littrature, sont distingues :
BEM (Mthode des lments frontires)
FEM (Mthode des lments finis)
MOC (Mthodes des caractristiques)
Diffrences finies
Mthodes faisant intervenir les R-fonctions

En rgle gnrale, la rsolution du problme fait intervenir plusieurs de ces
mthodes.

Le champ lectrique a deux composantes. La premire, gnre par la charge
despace, est gouverne par lquation de Poisson avec les conditions aux limites
uniformes de Dirichlet. Lautre composante est produite par la diffrence de
potentiels applique aux lectrodes, et est dcrite par les quations de Laplace avec
des conditions aux limites de Dirichlet non uniforme.
Les quations de Laplace sont rsolues par la BEM. Quand la BEM est utilise, il
suffit que la surface des lectrodes soit discrtise. Le problme est formul en
fonction de la densit de charge la surface des lectrodes, et le champ lectrique
en tout point du domaine peut tre calcul par une intgration numrique.
Pour trouver la seconde composante du champ lectrique, les quations de poissons
sont rsolues en utilisant la FEM afin de minimiser lordre du systme plutt que de
rsoudre directement les quations.
La troisime tape dans le processus de rsolution implique le calcul du flux de
charges en utilisant une intensit connue de champ lectrique. Plusieurs approches
existent, la plus usite tant la MOC.

1. Mthode des lments frontires

La BEM offre plusieurs avantages par rapport la FEM, ou la FDM (Finite Difference
Method). Un des atouts de cette mthode est un nombre plus petit dquation et une
rduction considrable du nombre de donnes ncessaires pour rsoudre un
problme. Elle sapplique trs bien aux problmes avec des domaines infinis, et est
V grad q grad
l
bq
Yt
Yq
+
|
.
|
\
|
E q grad
l
bq
Yt
Yq
r
+
|
.
|
\
|
capable doffrir des solutions continues sur le domaine. Les paramtres du problme
tel que le potentiel et lintensit du champ lectrique peuvent tre valu rapidement
en tout point.

Les solutions de Green des quations aux quations de Laplace en 2D sont :

(17)

pq 0
r
1
ln
2qq
1
Q) G(R, = Prs, Qrt (18)

o r
PQ
reprsente la distance entre les deux points P et Q, s le domaine, t la
frontire de s.
Lquation (17) est valable en tout point de lespace la condition que le paramtre ^
soit correct. En particulier, ^=1 lintrieur du domaine s. La solution en tout point
est exprime en fonction de la valeur la frontire et de sa driv. La solution peut
tre facilement dduite si ces deux conditions sont connues. Dans le cas gnral,
seulement une de ses conditions est spcifie, dans ce cas, lquation (17) devient
une quation intgrale avec une fonction inconnue.

Le long de la limite o le potentiel est connu, lquation (17) devient une quation
intgrale du premier ordre avec pour inconnue les drivs de la solution.

(18)

2. Mthode des lments finis

Bien que la BEM soit prfrable sur la FEM pour la rsolution des quations de
Laplace, les champs lectriques produits par la charge volumique sont mieux rsolus
par la FEM. En effet, dans ce cas, il ny pas de variation brutale du champ lectrique,
et donc, par voie de consquence, une discrtisation fine ne ce justifie pas. De plus,
la FEM est mieux adapte pour des problmes o la charge volumique nest pas
uniforme. Dans un tel cas, la BEM demanderait un temps de calcul plus important
pour lintgration.
Le principe de la FEM est de diviser le domaine en plusieurs mailles et ensuite,
dinterpoler par une fonction, la solution en chaque nud du maillage. Par exemple,
pour le problme considr, peuvent tre utilis :
des lments triangulaires
une interpolation linaire :

} (
c
c
c
c
=
t
dI Q) G(P,
n
V(Q) V(Q)
n
Q) G(P, ^V(P)
^(P) V(Q)dI
n
Q) G(P,
t
=
c
c
}

(19)

o y) (x, N y), (x, N y), (x, N
e
k
e
j
e
i
sont les fonctions de forme au nud i, j et k. V
i
, V
j
et V
k
sont les potentiels aux nuds. La densit de charge u(x,y) peut aussi tre interpol
de la mme manire.

Rsoudre lquation
l
q
V = V revient trouver un extremum la fonctionnelle, ce
qui peut sexprimer par :

(20)

o [S]et [T] sont obtenues par lassemblage des matrices globales. En appliquant la
condition dextremum :

0
V
F
i
=
c
c
i=1, 2, , n (21)
Nous obtenons ainsi, un jeu de n quations avec pour inconnues, la valeur du
potentiel en chaque nud.

3. Mthode des caractristiques

La MOC est base sur une technique o lquation aux drives partielles exprimant
la densit de charge, devient une quation diffrentielle classique le long dune ligne
de courant.

Les quations (10), (11), (13) et (14), deviennent en rgime stationnaire:

(22)

Le long de la ligne caractristique donne par lquation :
bE
dt
dr
= o r(t) est lquation paramtrique de la ligne caractristique.

En intgrant lquation (22):

(23)

k
e
k j
e
j i
e
i
V N V N V N V + + =
[ ][ ] [ ][ ]
(
(
(
(
(
(
(
(
=
n
2
1
n 2 1
n
2
1
n 2 1 0
u
u
u
T V ... V V
V
V
V
S V ... V V l
2
1
F
M M
0
l
q
q E = V
0
l
q
b
dt
dq
=
On peut connatre la densit de charge le long dune ligne caractristique :

(24)

o q
0
est la densit de charge au dpart de la ligne caractristique et b, la mobilit
ionique dans lair (suppose constante).
Le paramtre t peut tre interprt comme tant un temps fictif au cours duquel les
ions voyagent depuis llectrode jusqu un point donn. Cependant, cela ne veut
pas dire que la solution dpend du temps. Pour viter toute confusion, t peut tre
utilis comme un paramtre formel :
t=0 correspond la surface de llectrode. Quand t augmente, la trajectoire de
llectron est dtermine par lquation :
bE
dt
dr
= . (25)
La densit de charge diffrent point de cette trajectoire est calcule en substituant
la valeur de t en dans (24) :
La dernire quation peut tre utilise pour valuer la densit de charge dans le
domaine entier. Finalement, la densit peut tre calcule par :
qbE J = .
C. Exemples de modlisation

Au cours de nos recherches, nous avons rencontr plusieurs mthodes possibles
pour rsoudre les quations rgissant ce problme. Nous allons en prsenter deux.

1. Exemple faisant intervenir la BEM et la FEM
[3]

La technique utilise fait appel lutilisation de la BEM et la FEM pour dterminer la
distribution du champ lectrique. La MOC est quand elle utilise pour dterminer la
densit de charge.

La valeur de la densit de charge sur la surface de llectrode de travail est un autre
paramtre qui doit tre dtermin au cours du processus itratif. Ce processus
implique lutilisation dun double algorithme itratif :

1. Rsoudre les quations de Laplace par la BEM,

2. Faire une premire estimation de la densit de charge sur la surface de
llectrode de travail q
0
. Il ny a pas de mthode prcise pour cette estimation,
mais il existe des rsultats exprimentaux au sujet de la valeur du courant de
dcharge et des relations entre la densit de courant, le champ lectrique et la
densit de charge,

3. Supposer une distribution spatiale de la densit de charge q(r)

dans tout le
domaine. Une bonne estimation de cette distribution acclra la rsolution du
(
+
=
t
l
b
q
1
1
q
0 0
problme. Malheureusement, il ny a pas de mthode gnrale pour une
bonne estimation cest pourquoi, nimporte quelle technique peut tre utilise.
Par exemple, les auteurs utilisent ici une interpolation linaire entre la valeur
q
0
sur llectrode de travail et 0 sur la deuxime lectrode,

4. Calculer le champ lectrique gnr par la densit de charges en utilisant la
FEM,

5. Mettre jour la distribution de la densit de charge par la MOC,

6. Rpter les tapes 4 et 5 jusqu ce que la solution converge. Un cart de
1% sur la valeur du champ lectrique entre 2 itrations successives est utilis
comme condition darrt,

7. Comparer le champ lectrique surfacique sur llectrode de travail et le
champ lectrique dtermin par la relation de Peek,

8. Corriger lestimation de q
0
,
9. Rpter les tapes 3 8 jusqu ce que la valeur de lintensit du champ
lectrique sur llectrode de travail diffre de 3% par rapport la valeur fournie
par la relation Peek.

Ce procd itratif semble converger rapidement. Les tapes 3 8 constituent la
boucle extrieure. Quand la boucle extrieure est dans ses premires itrations, la
boucle intrieure (tape 4 6), demande gnralement 4 6 itrations pour
converger. Quand la boucle extrieure est prs de la solution, seulement 2 itrations
sont ncessaires pour que la boucle interne converge. Dans le but de diminuer le
temps CPU durant les premires boucles extrieures, la convergence de la boucle
intrieure nest pas contrle et seulement 1 itration est effectue (ltape 6 est
enleve). Quand la convergence de la boucle extrieure se fait plus prcise, ltape 6
est rintgre au processus. Bien videmment, plusieurs expriences sont
ncessaires pour estimer quand lalgorithme entier peut tre utilis. Cela a
gnralement lieu quand le nombre ditrations dans la boucle extrieure est de
lordre de 7. Gnralement, il faut entre 8 et 12 itrations lalgorithme pour
converger.

2. Exemple faisant intervenir la mthode des diffrences
finies
[2]

Le problme est formul en 2D. Le modle gomtrique est prsent en figure 2 et
correspond au modle exprimental de la figure 1. Llectrode de travail, dun
diamtre de 0,3 mm, est situe une distance de 50 mm de la surface de lautre
lectrode. Cette seconde lectrode est de forme circulaire et a un diamtre de 25
mm.
Figure III-1 Dcharge Couronne grce un dispositif double; (1) masse, (2) lectrode ionisante, (3)
support mtallique.

Figure III-2 Modlisation relle (a) et simplifie (b) de la gomtrie de l'lectrode.

t
Di
, i=1,2,3 dsigne la frontire sur laquelle est appliqus des conditions aux limites
de Dirichlet, respectivement llectrode de travail, llectrode cylindrique, et la masse.
La section circulaire de llectrode de travail est modlise par un carr de 0,3 mm
de cot. Llectrode cylindrique est modlise comme tant un octogone.
Les conditions suivantes sont imposes aux limites :

(26)

D3
D2 D1 c
t y) (x, 0 V
t t y) (x, V V
e =
e =
En accord avec lhypothse "Llectrode ionisante est un fil sur lequel la dcharge est
uniformment rpartie", la condition aux limites pour la charge q prend en compte le
fait quelle doit tre uniformment rpartie sur la surface de llectrode de travail :

(27)

o I
c
est la valeur estime du courant de dcharge par unit de longueur de
llectrode et E
c
est donn par la loi de Peek. De mme, la seconde hypothse nous
donne la valeur du champ lectrique sur llectrode de travail.

(28)

E
c
est dtermin par la relation suivante :

(29)

O k
1
et k
2
varient en fonction de la configuration tudie. Dans le cas dune
configuration fil - cylindre, k
1
=30,1x10
5
V/m et k
2
=0,0301 m
1/2
, ainsi E
c
=104kV/cm.

La mthode de rsolution utilise est de type FDM, dont le principe est de convertir
un systme non linaire (quations (1) (3) et (7)) en un systme itratif dquation :

(30)

(31)

(32)

o m est lindice ditration au temps t=nt
s
.
Lalgorithme est prsent en figure III-3.

D1
c c
c
c
t y) (x,
bE 8R
I
q q e = =
D1 c
t y) (x, E E e =
|
.
|
\
|
+ =
R
k
1 k E
2
1 c
0 E q bgrad
l
bq
Yt
Yq
n 1, m n m,
n 1, m n m,
= + +
|
|
.
|
\
|
r
l
q
nV
n m,
n m,
=
n m, n m,
gradV E =
r
Figure III-3 Schma suivi lors de la rsolution numrique.
A linstant t=0, les valeurs de q, V et E
r
en tout point de lespace D, sont dtermins
comme suit :
La densit de charge :

D1 c 0,0 D1 0,0
t y) (x, q q et t - D y) (x, 0 q e = e = (33)

Le potentiel V
0,0
est la solution de lquation de Laplace :

0 nV
0,0
= (34)

Le champ lectrique est dduit de (1) :

(35)

A linstant t=t
s
, i.e pour n=1 :

0,0 0,1 0,0 0,1
E E et q q = = et la charge q
1,1
peut tre dduite de lquation (30)
avec m=1 et n=1. Cette valeur de densit de charge est alors utilise pour
calculer V
1,1
partir de (31).
La premire itration sachve par le calcul de E
1,1
grce lquation (32).
Pour chacune des itrations suivantes, les valeurs de q
m-1,1
et E
m-1,1
sont
utilises pour calculer q
m,1
qui est alors utilis pour estimer V
m,1
et E
m,1
.
A chaque itration,les diffrences q
m,n
q
m-1,n
, V
m,n
V
m-1,n
et E
m,n
E
m-1,n
sont
compares avec les erreurs maximales admissibles (cest le critre de
convergence). Si ces valeurs sont infrieures la tolrance alors :

(36)

A linstant t = nt
s
, le mme schma est utilis.

D1 c 0,0 D1 0,0
t y) (x, E E et t - D y) (x, V V e = e =
n m, s n m, s n m, s
E ) t y, E(x, et V ) t y, V(x, , q ) t y, q(x, = = =
3. Conclusion

A la lecture des diffrents travaux que nous avons consult, il sest avr que les
rsultats de simulations obtenus sont parfois trs loin dtre valid
exprimentalement. En effet, une modlisation qui donnerait des rsultats plausibles
dans une configuration bien prcise nest jamais transportable dautres
configurations (gomtrie, pression, etc.). Il y a deux raisons cela.

Il est souvent suppos que les particules (plus particulirement les lectrons) sont en
quilibre avec le champ lectrique local hors, il existe au moins deux causes de non
quilibre:
prsence de parois
fortes variations de champ lectrique (tte de streamer, gaines)

Lhypothse dquilibre est utilise mais pas vrifie partout. Il existe cependant des
modles nutilisant pas cette hypothse. Cest le cas dans les modlisations utilisant
une quation de conservation des espces + deux quations supplmentaires (pour
la quantit de mouvement des lectrons, et pour leur nergie moyenne). Le
problme est alors de trouver des expressions des coefficients de transport et des
constantes de raction.

De plus, avant de modliser, il est ncessaire de bien comprendre les phnomnes
physiques mis en jeux. A ce jour, aucune tude na t ralise concernant le
dclenchement de deux streamer conscutifs. Et mme si la propagation du
streamer est assez bien comprise, il y a des carts entre les vitesses calcules et les
vitesses mesures.

Ainsi donc, les deux raisons principales majeures des checs dans la mise en
oeuvres de modlisation concerne donc :
Les modles trop simplistes
La mconnaissance des phnomnes mis en jeu.

Conclusion

Il existe plusieurs configurations gomtriques dlectrodes permettant la ralisation
de dcharges couronnes. Lutilisation dune configuration plutt quune autre est
conditionne par lapplication envisage. Bien souvent, les dcharges couronnes
sont employes la place de techniques plus polluantes. Leurs domaines
dapplication sont multiples et leur relative simplicit de mise en uvre leurs ouvre
de nouvelles perspectives. Cependant lutilisation des dcharges couronnes restent
conditionnes par les cot la fois financier et nergtique. En effet elles sont
encore inintressantes pour certains champs dapplication. Aussi, de nombreuses
tudes sont actuellement en cours afin de mieux matriser cette technologie. Parmi
celles ci, une large place est laisse la modlisation numrique des phnomnes
mis en jeu. Mais les rsultats obtenus ce jour ne sont concluants que dans des
configurations bien dfinies. Au-del des limitations machines, les carts entre
simulation et exprience sont bien souvent dus soit des modles trop simplistes,
soit une mconnaissance des phnomnes mis en jeu.

Bibliographie

[1] Conception et simulation dun racteur fil-cylindre dcharge couronne avec
barrire dilectrique adapt au traitement des oxydes dazote dans des effluents
marqus par un isotope, thse, Axel Vincent, Universit Paris VI, 2002
[2] Numrical modeling of combined corona-electrostatics fields, Armand Caron,
Lucian Dascalescu, IUT de Lorient, LAII-ESIP, 2003
[3] A numerical algorithm for simulation of the electric corona discharge in the triode
system, X. Deng et K. Adamiak, Department of Electrical and Computer Engineering,
University of Western Ontario, London, Ontario, Canada N6A 5B9, 1998
[4] Modlisation des dcharges filamentaires la pression atmosphrique,
prsentation, A. Bourdon (CORIA, Rouen), K. Hassouni (LIMHP, Villetaneuse), E.
Marode (LPGP-Suplec, Orsay), J. Paillol (LGE, Pau), P. Sgur (CPAT, Toulouse),
2003
[5] Hexamethyldisiloxane-Plasma coating of wood surfaces for creating water
reppelent characteristics, Denes A.R, Tshabalala M.A, Rowell R, Denes F, Young
R.A, 1999
[6] Modification of wood wettability by plasma and corona treatments, Podgorski L,
Chevet B, Onic L, Merlin A, 2000
[7] Contole actif dun coulement dair par un plasma froid surfacique, thse, Luc
Lger, LEA, universit de Poitiers, 2003
[8] Rhologie de panneaux composites bois/thermoplastiques sous chargement
thermomcanique : aptitude au postformage, thse, Michaud F, Universit
Bordeaux1, 2003
[9] Dpollution des effluents gazeux par plasma, Jean-Marie Cormier, Ahmed
Khacef, Olivier Motret, 2002

Site Web :
Date de consultation : 15 Janvier 2005
[1] http://www.paganetti.fr
[2] http://www.energie-plus.com
[3] http://www.valitec.fr
[4] http://www.lpgp.u-psud.fr
[5] http://plasmas.agmat.asso.fr
[6] http://www.maino.it
[7] http://www.ozone.ch
[8] http://www.info-systel.com
[9] http://www.plastiquesmagazine.com


ANNEXE 1
CLASSIFICATION DES PLASMAS
Il existe de nombreux types de plasmas que lon diffrencie selon diffrents critres.
Parmi ces critres nous trouvons la densit lectronique ne et la temprature
lectronique Te (ou nergie kTe) ou la longueur de Debye xD. Cette grandeur, qui
dfinit la distance sur laquelle les effets de champ lectrique dus une charge e ne
sont pas neutralis par un ensemble de charges de signe oppos, reprsente aussi
la limite de neutralit lectrique du plasma : si xD est infrieur aux dimensions du
plasma, celui-ci sera globalement neutre. xD est dfini par la relation suivante :

avec l0 (permittivit du vide) = 8,84.10-12 F.m-1
k (constante de Boltzmann) = 1,3806.10-23 J.K-1.molcule-1
La figure suivante rcapitule les diffrents types de plasmas selon les trois critres :
kTe, ne et xD :
Diffrentiation des plasmas par lnergie lectronique kTe, la densit lectronique ne
et la longueur de Debye xD.

Une autre forme de diffrentiation des plasmas est leur tat thermodynamique : on
distingue alors les plasmas lquilibre thermodynamique (ou plasmas thermiques)
des plasmas hors quilibre (plasmas froids) suivant que les tempratures qui les
caractrisent (lectronique, rotationnelle, vibrationnelle, cintique) sont gales ou
non :
Les plasmas thermiques sont gnrs par des torches plasma ou des arcs.
Les tempratures de ces plasmas sont typiquement comprises entre 5.000 K
et 50.000 K. Ces plasmas, dont la puissance varie de quelques kilowatts
quelques mgawatts, ncessitent la mise en jeu de courants importants.
Les plasmas hors quilibre sont gnralement produits par dcharge
luminescente, dcharge couronne ou courant inductif radiofrquence. Leur
temprature lectronique (1000 K 10000 K) est trs suprieure la
temprature du gaz, proche de lambiante (y 500 K). Dans ce type de plasma,
la majeure partie de lnergie nest pas convertie en chauffage du gaz mais en
ractivit chimique par production despces mtastables, dissocies et
ionises. Les puissances typiquement mises en jeu pour obtenir un plasma
hors quilibre vont de quelques watts quelques kilowatts.


ANNEXE 2
DOMAINE DUTILISATION DES
PLASMA FROIDS
POLYMERES

1) Polymrisation :
Schage de peintures, vernis, encres, adhsifs
Traitement de surfaces : papier, plastique, feuilles mtalliques,
contreplaqus, circuits imprgns...
Fabrication de matriaux composites: bois, btons imprgns, verre,
kevlar
Rsines radiodurcissables.

2) Rticulation :
Renforcement des proprits mcaniques, de la tenue en temprature,
de l'tanchit ...
Rubans, manchons..
Films thermortractables, isolants, mousses expanses ...
Vulcanisation de caoutchouc, d'lastomres...
Fabrication de disquettes et bandes magntiques.

3) Greffage:
Textiles synthtiques
Membranes permislectives
Polythylne thermoadhsif
Papier hyginique.

4) Dpolymrisation :
Tflon
Caoutchouc butyle
Polypropylne

PRODUITS MDICAUX ET PHARMACEUTIQUES

1) Strilisation de matriel mdico-chirurgical:
Articles en non-tisss (couvre-tables opratoires, masques..)
Articles en PVC
Pansements, seringues, gants en latex...

2) Strilisation de produits mdicaux:
Articles de conditionnement et d'emballages
Articles de laboratoires et de prlvement
Prothses et implants

3) Dcontamination:
Produits pharmaceutiques et cosmtologiques (poudres, talcs,
colorants, produits d'hygine...).


ENVIRONNEMENT

1) Dsinfection de boues d'gouts. Traitement de l'eau
Purification d'eau potable
Dcontamination des eaux uses
Strilisation de bouillons de culture
Dpollution chimique des eaux uses (limination des hydrocarbures
lgers, solvants pesticides).

2) Dpollution des fumes (S02, NOx, H2S).

3) Purification de l'air:
Neutralisation des odeurs
Elimination des composs organiques volants des fumes
Aseptisation d'ambiance par limination de germes.

4) Strilisation des dchets hospitaliers

5) Dpollution des sols

IAA

1) Produits et emballages :
Strilisation
Dcontamination

2) Produits :
Inhibition de la germination
Allongement de la dure de conservation
Retard de maturation
Dsinsectisation
Strilisation
Dcontamination bactrienne

Décharge Couronne-Application Et Modélisation

Transféré par

Droits d'auteur :

Formats disponibles

Décharge Couronne-Application Et Modélisation

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Décharge Couronne-Application Et Modélisation

Transféré par

Droits d'auteur :

Formats disponibles

DAILL Guillaume

Vous aimerez peut-être aussi