Modélisation et simulation numérique de la génération de plasma ...

6. Description de la méthode numérique et validationun nuage d’ions positifs. Entre ce nuage et l’anode, le champ faiblit. Cette diminution entraînel’extinction de l’avalanche électronique. Quelques électrons s’attachent aux neutres pour formerdes ions négatifs, mais en quantité très faible, Figure 6.11.Les ions positifs (O + 2 et N+ 2 ) sont répartis selon des créneaux successifs s’éloignant de l’anode.A chaque pulse de courant correspond un créneau d’ions positifs, créneaux dont les densitéss’homogénéisent par la suite pour former un nuage uniforme dont la densité est de quelque 10 10particules par cm 3 . C’est cette charge d’espace continue et importante qui provoque la diminutiondu champ maximal de 210 kV/cm initialement à 95 kV/cm lors de la phase établie de pulses.Lorsque le créneau issu de la dernière avalanche électronique est suffisamment éloigné de l’anode,le champ électrique retrouve un niveau suffisant permettant la création d’un nouveau pulse. Cenouveau pulse est également rendu possible par le détachement électronique qui, même si lesions négatifs sont en faible nombre, fournit un apport électronique suffisant, comme le rappellela Figure 6.7. On observe donc un régime établi constitué de pulses de courant successifs.Du point de vue numérique, on vérifie a posteriori que le champ électrique ne varie pas defaçon trop abrupte entre deux mailles de calcul voisines. Cela assure ainsi que les coefficients detransport sont bien estimés aux interfaces par la méthode décrite au paragraphe 6.2.2, évitantainsi de soumettre le problème à une source éventuelle d’instabilité numérique.6.4 ConclusionLa méthode numérique a été testée et validée sur une géométrie d’électrodes de type Sphère -Sphère. Elle est suffisamment précise pour simuler le développement des pulses de courant d’unecouronne positive. Cette étude sert donc de base aux travaux présentés dans les chapitres suivantset qui ont pour objet de simuler les décharges couronnes se produisant autour d’électrodes filaires.Dans un premier temps, on isole les couronnes positives des couronnes négatives afin d’étudierle comportement du modèle physique (et numérique) sur des cas plus simples que le dispositifexpérimental constitué de deux fils. Le prochain chapitre a donc pour objectif de cerner lesprincipales caractéristiques des décharges couronnes positives puis négatives afin de préparer lasimulation du cas expérimental.80