Multidimensional Numerical Simulation of Glow Discharge by Using the N-BEE-Time Splitting Method
|
|
- Darlene Simon
- 6 years ago
- Views:
Transcription
1 Plasma Science and Technology, Vol.14, No.9, Sep Multidimensional Numerical Simulation of Glow Discharge by Using the N-BEE-Time Splitting Method Benyssaad KRALOUA, Ali HENNAD Electrical Engineering Faculty, University of Sciences and Technology (USTO-MB), Oran, Algeria Abstract In this work, a new numerical technique is proposed for the resolution of a fluid model based on three Boltzmann moments. The main purpose of this technique is to calculate electric and physical properties in the non-equilibrium electric discharge at low pressure. The transport and Poisson s equations form a self-consistent model. This equation system is written in cylindrical coordinates following the geometric shape of a plasma reactor. Our transport equation system is discretized using the finite volume approach and resolved by the N-BEE explicit scheme coupled to the time splitting method. This programming structure reduces computation time considerably. The 2D code is carried out and tested by comparing our results with those found in literature. Keywords: 2D cylindrical fluid, glow discharge, N-BEE scheme, time splitting method PACS: v, Ex, Hc DOI: / /14/9/06 1 Introduction During the past few years, glow discharge been used for a large number of important applications, such as water treatment, sterilization medical material, low energy consumption source light, screen plasma for micro and nanoelectronic industries [1 7]. In this contribution, a recent numerical scheme is proposed for the two-dimensional resolution of the transport equations of the second-order fluid model in non-equilibrium electric discharges at low pressure. Our model is based on the first three moments of Boltzmann s equations for ions and electrons, coupled to Poisson s equation. By considering the geometry of a plasma reactor, the transport equations are written in cylindrical coordinates. In this work, we consider that glow discharge arises between two metal electrodes circular, plane and parallel and radius R. These two electrodes are spaced a distance d and the wall of our reactor is a dielectric (see Fig. 1). The bidimensional transport equations are transformed into one-dimensional equations (driftdiffusion equation), while a time splitting approach [8,9] is adopted. The transport equations are then solved by using the N-BEE explicit scheme [10,11]. The numerical solution of the transport equations depends essentially on the choice of the boundary and initial conditions. On the discharge symmetrical axis, Neumann boundary conditions are mainly used for the electron and ion densities, electric potential and electron temperature. Dirichlet boundary conditions are used at the electrodes and the dielectric walls. Schematic representation of the glow discharge re- Fig.1 actor 2 Fluid model presentation In the present model, the transport equations derived from the first three moments of Boltzmann s equation are written only for electrons and positive ions. The electron and ion densities are described by continuity equations expressed in cylindrical geometry [12 14] : e t + J ez + 1 (r J er ) = S e, (1) + + J +z + 1 (r J +r ) = S +. (2) t In the model, the momentum conservation equation is replaced by the drift-diffusion approximation; hence, the transport equation is represented by two separate
2 Benyssaad KRALOUA et al.: Multidimensional Numerical Simulation of Glow Discharge terms, i.e., drift and diffusion terms. The electron and ion density fluxes are given by: J ez = µ e n e E z D e e, (3) J er = µ e n e E r D e e r, (4) J +z = +µ + n + E z D + +, (5) J +r = +µ + n + E r D + + r. (6) Here J e and J + are electron and ion density fluxes, respectively. Poisson s Eq. (7) is resolved for the determination of the electric field and the discharge potential. In this equation, V is the electric potential, e is the elementary charge, and ε 0 is the permittivity of free space. 2 V r V r + 2 V r 2 = e (n + n e ). (7) ε 0 The axial E z and radial E r components of the electric field strength are then found from the partial derivatives of the potential function, and are given by E z = V and E r = V r. (8) In our model, the ionization coefficient of neutral particles is a function of electron temperature. The third moment of Boltzmann s equation (electron energy Eq. (9)) is used to calculate the electron temperature profile). ε t + J εz + 1 (r J εr ) = S ε, (9) J εz = 5 3 µ en ε E z 5 3 D ε e, (10) J εr = 5 3 µ en ε E r 5 3 D ε e r, (11) where n ε is the electron density energy and J ε is the electronic energy density flux. The right-hand sides of Eqs. (1), (2) and (9) include all source terms, S e, S + and S ε are calculated with the corresponding coefficients rates tabulated in Table 1. Table 1. Transport parameters used in the present simulation [4,10,11] Transport parameters Values N (cm 3 ) µ e (cm 2 V 1 s 1 ) D (cm 2 s 1 ) µ (cm 2 V 1 s 1 ) E i 24 (ev) H i (ev) γ V (V) T c 0.5 (ev) K i (cm 3 s 1 ) d (cm) R 5.08 (cm) ( S e = S + = K i Nn e exp E i k b T e ( S ε = ej + E H i K i Nn e exp E i k b T e ), (12) ). (13) The pre-exponential factor, energy of process collisional activation, energy lost by ionization and density of the neutral gas are K i, E i, H i, and N, respectively. Our equation system is solved by means of two electrical and physical approximations: local electric field approximation for ions and local energy approximation for electrons [12,15]. 3 Boundary conditions The use of these boundary conditions for the above mentioned transport equations is essential for the description of our problem. As a matter of fact, we used mainly Neumann boundary conditions on the symmetrical axis: for the gradient particle densities: = 0, for the electric potential: T e r = 0. V r + r = 0 and e r = 0, for the electron temperature: The Dirichelet boundary conditions for the electric potential at electrodes are: V =0 V at the grounded electrode (anode) and V = 77.4 V at the powered electrode (cathode). At the anode; the gradient ion density + = 0, and for the electron density n 0 = 0, for the gradient of the electron temperature is T e = 2e V 5k b. At the dielectric wall; the gradient potential is supposed to be a floating potential. This potential is determined in an auto-coherent approach based on the flow of the charged particles which affects the dielectric according to the following relation: V r = e ε 0 t 0 (J + J e )dt. (14) We assume that the recombination of the charged particles on the wall is immediately completed and the conductivity of dielectric is perfect. Consequently, we suppose that the values of the electron and ion densities on the dielectric wall are nulls (n + =0 and n e =0). The gradient of the electron temperature is supposed to be floating: T e r = 2e V 5k b r. (15) At the cathode, the electron temperature is T c = 0.5 ev, the gradient ion density is + /=0 and the electron density is n e = 0. In this study, the secondary electron emission coefficient, γ 0, is equal to Fig. 2 shows the schematic representation of the different boundary conditions used in our model. 803
3 Plasma Science and Technology, Vol.14, No.9, Sep Here the flux J m i+1/2,j is: J m i+1/2,j = nm i,j (1 ν i) ( n m i+1,j n m i,j) ϕ (ri,j, ν i ). (19) The condition of Current Frederic Lewy (CFL) is: ν i = t z ν e (20) The flux limiter ϕ(r i,j, ν i ) of Eq. (17) is: [ ( ϕ(r i,j, ν i ) = Max 0, Min 1, 2r ) ( i,j, Min r i,j, ν i By using the density ratio: 2 1 ν i (21) )]. Fig.2 Schematic representation of the boundary conditions used in our model 4 N-BEE-time splitting method The numerical resolution of the transport equations of our 2D auto-coherent fluid model is laborious. In this section, we use the N-BEE algorithm coupled to the time splitting method for the numerical solution of the transport equations [8 11]. The splitting methods are well-known to solve this kind of multidimensional and multi-physical problems. In effect, the equations written initially in 2D are transformed into one-dimensional equations in each spatial direction. So by using this time-splitting method, it is possible to find n m i,j firstly for time step m with no knowledge of the numerical approximation, then for the new time step (m +1) and finally calculate Ji,j m using the solution that has already been computed for n m+1 i,j. In this section, we present in detail the N-BEE algorithm used in this work for the resolution of the transport equation cited above [16,17]. t + J + 1 (r J) = S. (16) By using the time splitting method, we have seen that the transport equations written initially in 2D can be transformed into 1D equations in each spatial direction: t + J = 0, (17.a) t + 1 (r J) = 0, (17.b) t = S. n m+1 i,j n m i,j t + J m i+1/2,j J m i 1/2, j z + r j+1/2j m i,j+1/2 r j 1/2J m i,j 1/2 r j r (17.c) = S m i,j. (18) r i,j = nm i,j nm i 1,j n m i+1,j. (22) nm i,j The solution of Eq. (17.a) is as follows: n m+1 i,j = n m i,j ν i (Ji+1/2,j m J i 1/2,j m ). (23) From Fig. 3, we have to compute the first step and then get the result of n m+1 i,j that will be considered as a further initial condition for the second step, where the n m+1 i,j solution is obtained. In the third step, we have put Si,j m as a further initial condition and get finally the result of n m+1 i,j. Fig.3 Time splitting method In this figure, the scheme represents the different steps of the transport equations resolution of the fluid model using the time splitting technique. With t, being the time step n m i,j and the discretization in time of n i,j at time m t, z and r are space steps and Ji,j m is the discretization in space following the nodes of the grid named i and j. 5 Results and interpretations A glow discharge reactor has a cylindrical symmetry and therefore the 2D geometry can be reduced to 1D for the fluid model [18]. This reduction enables a comparison with the 1D numerical simulation results obtained by LIN and ADOMAITIS [19,20] in which a good agreement can be noticed. Fig. 4 shows one-dimensional distributions of the steady-state electron and ion densities where the three different regions characterizing glow discharge (cathode 804
4 Benyssaad KRALOUA et al.: Multidimensional Numerical Simulation of Glow Discharge and anode sheaths and positive column) appear clearly. The anode region is characterized by a relatively important ion density compared to the electron density. The positive column associated with the plasma may be defined by electron and ion densities which are quasi equal characterising the plasma neutrality. Consequently, in the positive column the net space charge is negligible. The density of the species charged in this region has a maximum value. In the cathode region, the electron density is very low compared to the density of the ions. This gradient of density in this region is due to the fact that the electrons move away much more quickly from the cathode region than the ions in the presence of a potential gradient. This causes the depopulation of the electrons in this region. Fig.4 One-dimensional distributions of electron and ion densities at the At the level of the cathode sheath, the values of electron density resulting from our model and also from the model of LIN and ADOMAITIS [19 20] are different from zero, because this zone represents a source of electron production by the secondary emission process. Figs. 5, 6 and 7 represent the one-dimensional distributions of the potential, the electric field and the electron temperature in the steady state of the glow discharge. We notice in Fig. 5 an important gradient of potential in the region of the cathode sheath. This potential gradient is one of the characteristics of the glow discharge. In the positive column and the anode region, the potential is almost constant. The value of the potential in these two regions is nearly equal to the potential applied to the anode. This behaviour is normal, because of the value of the net space particle density that offers toward zero in the positive column. In Fig. 6, the variation of the field in the cathode sheath is linear because of the potential gradient. The field is quasi null in the positive column because the potential is almost unvaried. In the anode region, we can observe an inversion of the electric field because the gradient of plasma density is sufficiently important in this region. Fig. 7 shows the spatio-temporal distribution of the steady-state electron temperature over the whole distance. From this distribution, one can deduce the electron energy. We can easily see that the electron energy in the cathode region is relatively important because of presence of the relatively intense electric field. Figs. 8 and 9 represent contour lines of the electron and ion densities at a steady state where the behaviour of the luminescent discharge is correctly represented by its three regions. The drift velocity of electrons is 100 times larger than that of ions. We note that the contours level lines show the rapid movement of electrons of the cathode region. Fig.7 One-dimensional distributions of electron temperature at the Fig.5 One-dimensional distributions of electric potential at the Fig.6 One-dimensional distributions of electric field at the Fig.8 Bi-dimensional contour of electron density at the 805
5 Plasma Science and Technology, Vol.14, No.9, Sep due to the increase of electronic temperature. We also notice that the variation of the electron temperature in the positive column and the anodic region is almost independent of the value of the applied voltage. This is reasonable because of the nature of the electric field in both two regions. Fig.9 Bi-dimensional contour of ion density at the steady state (color online) The contour lines of the steady-state electric potential are illustrated in Fig. 10. One can see a potential gradient relatively important at the cathode sheath level. In the positive column region where the plasma is formed, the potential is practically constant. Fig.10 Bi-dimensional contour of electric potential at the Fig. 11 represents the contour level lines of the steady-state axial electric field, which confirms the relatively high value of the electric potential at the cathode sheath. In the plasma column and near the anode region, the axial field is almost negligible or null because of the net space charge. We can notice on contour level lines of the steadystate radial electric field in Fig. 12 that this component of the field presents an asymmetric distribution with respect to the symmetry axis. Due to this fact, the resulting radial field is null in the inter-electrodes gap. The main role of this component is to maintain the glow discharge confined around the symmetry axis. Fig. 13 represents contour lines of the steady-state electron temperature. We can see the profile of expected electron temperature because of the high electric field region near the cathode. Generally, the increase of electron density is produced by secondary emission 806 Fig.11 Bi-dimensional contour of axial electric field at the Fig.12 Bi-dimensional contour of transversal radial field at the Fig.13 Bi-dimensional contour of the steady-state electron temperature (color online)
6 Benyssaad KRALOUA et al.: Multidimensional Numerical Simulation of Glow Discharge 6 Conclusion In this paper, we have developed a bidimensional model of the glow discharge being maintained by a term source of electrons and ions uniformly produced by a secondary electron emission at the cathode. The N-BEE algorithm based on the finite difference method coupled to the time splitting technique is a good approach to resolve the bi-dimensional transport equation of the glow discharge model at low pressure. Our simulation of glow discharge with a symmetry axis are in good agreement with the 1D work of LIN and ADOMATIS [19,20]. References 1 Sugama C, Tochikubo F, Uchida S. 2006, J. Appl. Phys., 45: Arif Malik M. 2010, Plasma Chem. Plasma Process., 30: 21 3 Moisan M, Barbeau J, Crevier M C, et al. 2002, Pure Appl. Chem., 74: Lee H W, Park G Y, Seo Y S, et al. 2011, J. Phys. D: Appl. Phys., 44: Arkhipenko V I, Kirillov A A, Callegari T, et al. 2009, IEEE Transactions Plasma Science, 37: Duan Y, Huang C, Yu Q. 2005, IEEE Transactions Plasma Science, 33: Samanta M, Srivastava A K, Sharma S, et al. 2010, Electrical Discharge as an Inspection Method for Imperfect Plasma Display Cells. J. Phys. Conf. Ser., 208: Bott A. 2010, Improving the time-splitting errors of one-dimensional advection schemes in multidimensional applications. Atmospheric Research, 97: 619; Boris J P, Book D L. 1973, Flux-corrected transport. I. SHASTA, A fluid transport algorithm that works, J. Comput. Phys., 11: 38 9 Tahmouch G. 1995, A numerical multidimensionalmodelization and simulation of the interaction dynamics between charged particles and neutral particles during the transition to a point-plane transitory discharge in gas [Ph.D]. University Louis Pasteur, France 10 Bokanowski O, Zidani H. 2007, Journal of Scientific Computing, 30: 1 11 Bokanowski O, Martin S, Munos R, et al. 2006, Applied Numerical Mathematics, 56: Meyyappan M, Kreskovsky J P. 1990, J. Appl. Phys., 68: Lee C, Graves D B, Lieberman M A, et al. 1994, J. Electrochem, Soc., 141: Lee I, Graves D B, Lieberman M A. 2008, Plasma Sources Sci. Technol., 17: Kraloua B. 2010, Two dimensional modeling of the electrical discharges in low pressure: Second order fluid model [Ph.D]. University Sciences and Technology Mohamed Boudiaf Oran, Algeria 16 Kraloua B, Hennad A. 2008, Transport equations resolution by N-BEE anti-dissipative scheme in 2D model low pressure glow discharge. AIP Conf. Proc., 1047: Kraloua B, Hennad A, Yousfi M. 2009, Application anti-dissipative scheme N-BEE for 2D characterization glow discharge. IREE, 4: Kraloua B, Hennad A. 2010, Bidimensional Modelling Non-Equilibrium Fluid Model of Glow Discharge at Low Pressure. IREE, 5: Lin Y H, Adomatis R A. 1998, Physics Letters A, 243: Lin Y H. 1999, From Detailed Simulation to Model Reduction: Development of Numerical Tools for a Plasma Processing Application [Ph.D]. Faculty graduate school of the University of Maryland, USA (Manuscript received 13 January 2011) (Manuscript accepted 5 January 2012) address of corresponding anthor Ali HENNAD: ali hennad@yahoo.fr 807
PIC/MCC Simulation of Radio Frequency Hollow Cathode Discharge in Nitrogen
PIC/MCC Simulation of Radio Frequency Hollow Cathode Discharge in Nitrogen HAN Qing ( ), WANG Jing ( ), ZHANG Lianzhu ( ) College of Physics Science and Information Engineering, Hebei Normal University,
More informationThree-dimensional modelling of the DC glow discharge using the second order fluid model
Hocine TEBANI, Ali HENNAD University of Sciences and Technology of Mohamed Boudiaf USTOMB, BP 155 El M'Naouer, 31 Oran, Algeria Three-dimensional modelling of the DC glow discharge using the second order
More informationOne dimensional hybrid Maxwell-Boltzmann model of shearth evolution
Technical collection One dimensional hybrid Maxwell-Boltzmann model of shearth evolution 27 - Conferences publications P. Sarrailh L. Garrigues G. J. M. Hagelaar J. P. Boeuf G. Sandolache S. Rowe B. Jusselin
More informationSimulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 5 1 SEPTEMBER 2003 Simulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma Doosik Kim
More informationDepartment of Physics, University Abou Bekr Belkaid, Tlemcen, Algeria
Chemical Analysis of NO x Removal Under Different Reduced Electric Fields A. HADDOUCHE, M. LEMERINI Department of Physics, University Abou Bekr Belkaid, Tlemcen, Algeria Abstract This work presents a chemical
More informationCharacteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges
Plasma Science and Technology, Vol.6, No.6, Jun. 204 Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges M. M. HATAMI, A. R. NIKNAM 2 Physics Department
More informationStudy of a Micro Hollow Cathode Discharge at medium argon gas pressure
Study of a Micro Hollow Cathode Discharge at medium argon gas pressure Claudia LAZZARONI Antoine ROUSSEAU Pascal CHABERT LPP Ecole Polytechnique, Palaiseau, FRANCE Nader SADEGHI LSP Grenoble, FRANCE I-V
More informationKinetic simulation of the stationary HEMP thruster including the near field plume region
Kinetic simulation of the stationary HEMP thruster including the near field plume region IEPC-2009-110 Presented at the 31st International Electric Propulsion Conference, University of Michigan Ann Arbor,
More informationSimple examples of MHD equilibria
Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will
More informationContents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments
Contents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments IEC:! IEC at UIUC modified into a space thruster.! IEC has
More informationNonequilibrium discharges in air and nitrogen plasmas at atmospheric pressure*
Pure Appl. Chem., Vol. 74, No. 3, pp. 337 347, 2002. 2002 IUPAC Nonequilibrium discharges in air and nitrogen plasmas at atmospheric pressure* Charles H. Kruger, Christophe O. Laux, Lan Yu, Denis M. Packan,
More informationNumerical Simulation of Townsend Discharge, Paschen Breakdown and Dielectric Barrier Discharges Napoleon Leoni, Bhooshan Paradkar
Numerical Simulation of Townsend Discharge, Paschen Breakdown and Dielectric Barrier Discharges Napoleon Leoni, Bhooshan Paradkar HP Laboratories HPL-2009-234 Keyword(s): Townsend Discharge, Paschen Breakdown,
More informationSimulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 7 1 APRIL 2004 Simulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
More informationEffect of Applied Electric Field and Pressure on the Electron Avalanche Growth
Effect of Applied Electric Field and Pressure on the Electron Avalanche Growth L. ZEGHICHI (), L. MOKHNACHE (2), and M. DJEBABRA (3) () Department of Physics, Ouargla University, P.O Box.5, OUARGLA 3,
More informationElectron Transport Behavior in a Mirror Magnetic Field and a Non-uniform Electric Field
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 207 212 c International Academic Publishers Vol. 35, No. 2, February 15, 2001 Electron Transport Behavior in a Mirror Magnetic Field and a Non-uniform
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationSPECTRAL INVESTIGATION OF A COMPLEX SPACE CHARGE STRUCTURE IN PLASMA
SPECTRAL INVESTIGATION OF A COMPLEX SPACE CHARGE STRUCTURE IN PLASMA S. GURLUI 1, D. G. DIMITRIU 1, C. IONITA 2, R. W. SCHRITTWIESER 2 1 Faculty of Physics, Al. I. Cuza University, 11 Carol I Blvd., RO-700506
More informationCurrent sheath formation in the plasma focus
Plasma Science and Applications (ICPSA 2013) International Journal of Modern Physics: Conference Series Vol. 32 (2014) 1460321 (8 pages) The Author DOI: 10.1142/S2010194514603214 Current sheath formation
More informationMatti Laan Gas Discharge Laboratory University of Tartu ESTONIA
Matti Laan Gas Discharge Laboratory University of Tartu ESTONIA Outline 1. Ionisation 2. Plasma definition 3. Plasma properties 4. Plasma classification 5. Energy transfer in non-equilibrium plasma 6.
More informationCONTROL OF UNIFORMITY IN CAPACITIVELY COUPLED PLASMAS CONSIDERING EDGE EFFECTS*
CONTROL OF UNIFORMITY IN CAPACITIVELY COUPLED PLASMAS CONSIDERING EDGE EFFECTS* Junqing Lu and Mark J. Kushner Department of Electrical and Computer Engineering at Urbana-Champaign mjk@uiuc.edu, jqlu@uiuc.edu
More informationCZ České Budějovice, Czech Republic b Technical University of Liberec, Department of Materials Science, Hálkova 6, D Dresden, Germany
INVESTIGATION OF ELECTRIC CONDITIONS IN THE VICINITY OF CARBON NANOTUBES GROWN IN A DC PLASMA SHEATH J. Blažek a, P. Špatenka b, Ch. Taeschner c, A. Leonhardt c a University of South Bohemia, Department
More informationMICRODISCHARGES AS SOURCES OF PHOTONS, RADICALS AND THRUST*
MICRODISCHARGES AS SOURCES OF PHOTONS, RADICALS AND THRUST* Ramesh Arakoni a) and Mark J. Kushner b) a) Dept. Aerospace Engineering b) Dept. Electrical and Computer Engineering Urbana, IL 61801 USA mjk@uiuc.edu
More informationMWP MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES
MWP 1.9 MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES Insook Lee, D.B. Graves, and M.A. Lieberman University of California Berkeley, CA 9472 LiebermanGEC7 1 STANDING WAVES
More informationModélisation de sources plasma froid magnétisé
Modélisation de sources plasma froid magnétisé Gerjan Hagelaar Groupe de Recherche Energétique, Plasma & Hors Equilibre (GREPHE) Laboratoire Plasma et Conversion d Énergie (LAPLACE) Université Paul Sabatier,
More informationNumerical study of a Double Stage Hall Effect Thruster
Numerical study of a Double Stage Hall Effect Thruster IEPC--1 Presented at the 9 th International Electric Propulsion Conference, Princeton University, October 1 November, C. Boniface *, G.J.M Hagelaar,
More informationInfluence of vibrational kinetics in a low pressure capacitively coupled hydrogen discharge
Influence of vibrational kinetics in a low pressure capacitively coupled hydrogen discharge L. Marques 1, A. Salabas 1, G. Gousset 2, L. L. Alves 1 1 Centro de Física dos Plasmas, Instituto Superior Técnico,
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationMAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT
MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California
More informationFigure 1.1: Ionization and Recombination
Chapter 1 Introduction 1.1 What is a Plasma? 1.1.1 An ionized gas A plasma is a gas in which an important fraction of the atoms is ionized, so that the electrons and ions are separately free. When does
More informationStudies of the cathode sheath of a low pressure hollow cathode discharge
Journal of Physics: onference Series Studies of the cathode sheath of a low pressure hollow cathode discharge To cite this article: G Petraconi et al 0 J. Phys.: onf. Ser. 370 004 View the article online
More informationGeneration and diagnostics of atmospheric pressure dielectric barrier discharge in argon/air
Indian Journal of Pure & Applied Physics Vol. 55, February 017, pp. 155-16 Generation and diagnostics of atmospheric pressure dielectric barrier discharge in argon/air R Shrestha a,b *, D P Subedi a, R
More information2D Hybrid Fluid-Analytical Model of Inductive/Capacitive Plasma Discharges
63 rd GEC & 7 th ICRP, 2010 2D Hybrid Fluid-Analytical Model of Inductive/Capacitive Plasma Discharges E. Kawamura, M.A. Lieberman, and D.B. Graves University of California, Berkeley, CA 94720 This work
More informationStudy of DC Cylindrical Magnetron by Langmuir Probe
WDS'2 Proceedings of Contributed Papers, Part II, 76 8, 22. ISBN 978-737825 MATFYZPRESS Study of DC Cylindrical Magnetron by Langmuir Probe A. Kolpaková, P. Kudrna, and M. Tichý Charles University Prague,
More informationInfluence of Axial Magnetic Field on the Electrical Breakdown and Secondary Electron Emission in Plane-Parallel Plasma Discharge
Vol:5, No:8, 211 Influence of Axial Magnetic Field on the Electrical Breakdown and Secondary Electron Emission in Plane-Parallel Plasma Discharge Sabah I. Wais, Raghad Y. Mohammed, Sedki O. Yousif International
More informationThe Effect of Discharge Characteristics on Dielectric Barrier Discharges According to the Relative Permittivity
, pp.21-27 http://dx.doi.org/10.14257/astl.2017.145.05 The Effect of Discharge Characteristics on Dielectric Barrier Discharges According to the Relative Permittivity Don-Kyu Lee Electrical Engineering,
More informationDepartment of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, Texas 78712, USA
1 MAGNETIZED DIRECT CURRENT MICRODISCHARGE, I: EFFECT OF THE GAS PRESSURE Dmitry Levko and Laxminarayan L. Raja Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at
More informationKinetic Simulations of Ion Beam Neutralization
Kinetic Simulations of Ion Beam Neutralization O. Chang and J. Wang Astronautical Engineering Department University of Southern California Los Angeles, CA 90089-1192, USA Abstract. Full particle PIC simulations
More informationPlasma transport around dust agglomerates having complex shapes
JOURNAL OF APPLIED PHYSICS VOLUME 83, NUMBER 11 1 JUNE 1998 Plasma transport around dust agglomerates having complex shapes Eric R. Keiter a) and Mark J. Kushner b) Department of Electrical and Computer
More informationPlasma Optimization in a Multicusp Ion Source by Using a Monte Carlo Simulation
Journal of the Korean Physical Society, Vol. 63, No. 7, October 2013, pp. 0 0 Plasma Optimization in a Multicusp Ion Source by Using a Monte Carlo Simulation M. Hosseinzadeh and H. Afarideh Nuclear Engineering
More informationFINAL REPORT. DOE Grant DE-FG03-87ER13727
FINAL REPORT DOE Grant DE-FG03-87ER13727 Dynamics of Electronegative Plasmas for Materials Processing Allan J. Lichtenberg and Michael A. Lieberman Department of Electrical Engineering and Computer Sciences
More informationPIC-MCC simulations for complex plasmas
GRADUATE SUMMER INSTITUTE "Complex Plasmas August 4, 008 PIC-MCC simulations for complex plasmas Irina Schweigert Institute of Theoretical and Applied Mechanics, SB RAS, Novosibirsk Outline GRADUATE SUMMER
More informationNumerical modeling of a cutting torch
Journal of Physics: Conference Series OPEN ACCESS Numerical modeling of a cutting torch To cite this article: B R Mancinelli et al 2014 J. Phys.: Conf. Ser. 511 012071 View the article online for updates
More informationModeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range
Modeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range Dr.-Ing. Frank H. Scharf CST of America What is a plasma? What is a plasma? Often referred to as The fourth
More informationModeling nonthermal plasmas generated in glow discharges*
Pure Appl. Chem., Vol. 71, No. 10, pp. 1837±1844, 1999. Printed in Great Britain. q 1999 IUPAC Modeling nonthermal plasmas generated in glow discharges* I. Revel, Ph. Belenguer, J. P. Boeuf and L. C. Pitchford²
More informationMODELING OF AN ECR SOURCE FOR MATERIALS PROCESSING USING A TWO DIMENSIONAL HYBRID PLASMA EQUIPMENT MODEL. Ron L. Kinder and Mark J.
TECHCON 98 Las Vegas, Nevada September 9-11, 1998 MODELING OF AN ECR SOURCE FOR MATERIALS PROCESSING USING A TWO DIMENSIONAL HYBRID PLASMA EQUIPMENT MODEL Ron L. Kinder and Mark J. Kushner Department of
More informationEquilibrium model for two low-pressure electronegative plasmas connected by a double layer
PHYSICS OF PLASMAS 13, 093504 2006 Equilibrium model for two low-pressure electronegative plasmas connected by a double layer P. Chabert, a N. Plihon, C. S. Corr, and J.-L. Raimbault Laboratoire de Physique
More informationIntroduction. Chapter Plasma: definitions
Chapter 1 Introduction 1.1 Plasma: definitions A plasma is a quasi-neutral gas of charged and neutral particles which exhibits collective behaviour. An equivalent, alternative definition: A plasma is a
More informationKinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas
Kinetic Solvers with Adaptive Mesh in Phase Space for Low- Temperature Plasmas Vladimir Kolobov, a,b,1 Robert Arslanbekov a and Dmitry Levko a a CFD Research Corporation, Huntsville, AL 35806, USA b The
More informationTwo-dimensional Fluid Simulation of an RF Capacitively Coupled Ar/H 2 Discharge
Two-dimensional Fluid Simulation of an RF Capacitively Coupled Ar/H 2 Discharge Lizhu Tong Keisoku Engineering System Co., Ltd., Japan September 18, 2014 Keisoku Engineering System Co., Ltd., 1-9-5 Uchikanda,
More informationKinetic modelling of the jet extraction mechanism in spherical IEC devices
Kinetic modelling of the jet extraction mechanism in spherical IEC devices Type of activity: Standard study 1 Background & Study Motivation 1.1 Introduction Inertial Electrostatic Confinement (IEC) devices
More informationFLASH CHAMBER OF A QUASI-CONTINUOUS VOLUME SOURCE OF NEGATIVE IONS
FLASH CHAMBER OF A QUASI-CONTINUOUS VOLUME SOURCE OF NEGATIVE IONS P.A. Litvinov, V.A. Baturin * Institute of Applied Physics, National Academy of Science of Ukraine, 58 Petropavlovskaya St. Sumy, 40030
More informationHong Young Chang Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea
Hong Young Chang Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea Index 1. Introduction 2. Some plasma sources 3. Related issues 4. Summary -2 Why is
More informationPRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING Second Edition MICHAEL A. LIEBERMAN ALLAN J, LICHTENBERG WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION CONTENTS PREFACE xrrii PREFACE
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
More informationKINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS
KINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS Ralf Peter Brinkmann, Dennis Krüger Fakultät für Elektrotechnik und Informationstechnik Lehrstuhl für Theoretische Elektrotechnik Magnetized low
More informationEffect of an axial magnetic field on a DC argon arc
Vol 17 No 2, February 2008 c 2008 Chin. Phys. Soc. 1674-1056/2008/17(02)/0649-06 Chinese Physics B and IOP Publishing Ltd Effect of an axial magnetic field on a DC argon arc Li Lin-Cun( ) and Xia Wei-Dong
More informationSIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(0,n) MODES*
25th IEEE International Conference on Plasma Science Raleigh, North Carolina June 1-4, 1998 SIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(,n) MODES* Ron L. Kinder and Mark J.
More informationTwo-dimensional Numerical Simulation of a Planar Radio-frequency Atmospheric Pressure Plasma Source
Two-dimensional Numerical Simulation of a Planar Radio-frequency Atmospheric Pressure Plasma Source Lei Wang 1*, Gheorghe Dscu, Eusebiu-Rosini Ionita, Christophe Leys 1, Anton Yu Nikiforov 1 1 Department
More informationContents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations
PIC Method for Numerical Simulation Ninad Joshi NNP Group 1 Contents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations Motivation 3 Particle simulation Ion beams and Plasmas Accelerators
More informationNumerical Study on Influences of Barrier Arrangements on Dielectric Barrier Discharge Characteristics
504 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 31, NO. 4, AUGUST 2003 Numerical Study on Influences of Barrier Arrangements on Dielectric Barrier Discharge Characteristics Woo Seok Kang, Jin Myung Park,
More informationA New Model of Investigating the Electric Field in Dielectric Liquid for Streamer Initiation
A New Model of Investigating the Electric Field in Dielectric Liquid for Streamer Initiation E A. El-Zein and M. Talaat Electrical Power & Machines Department, Faculty of Engineering, Zagazig University,
More informationProfiling and modeling of dc nitrogen microplasmas
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 5 1 SEPTEMBER 2003 Profiling and modeling of dc nitrogen microplasmas Chester G. Wilson a) and Yogesh B. Gianchandani EECS Department, University of Michigan,
More informationEffect of small amounts of hydrogen added to argon glow discharges: Hybrid Monte Carlo fluid model
PHYSICAL REVIEW E, VOLUME 65, 056402 Effect of small amounts of hydrogen added to argon glow discharges: Hybrid Monte Carlo fluid model Annemie Bogaerts* and Renaat Gijbels Department of Chemistry, University
More informationAssessment of fluctuation-induced and wall-induced anomalous electron transport in HET
Assessment of fluctuation-induced and wall-induced anomalous electron transport in HET IEPC-2015-418 Presented at Joint Conference of 30th International Symposium on Space Technology and Science 34th International
More informationElectrical Discharges Characterization of Planar Sputtering System
International Journal of Recent Research and Review, Vol. V, March 213 ISSN 2277 8322 Electrical Discharges Characterization of Planar Sputtering System Bahaa T. Chaid 1, Nathera Abass Ali Al-Tememee 2,
More informationCandidacy Exam Department of Physics February 6, 2010 Part I
Candidacy Exam Department of Physics February 6, 2010 Part I Instructions: ˆ The following problems are intended to probe your understanding of basic physical principles. When answering each question,
More informationNumerical Study of Power Deposition, Transport and Acceleration Phenomena in Helicon Plasma Thrusters
Numerical Study of Power Deposition, Transport and Acceleration Phenomena in Helicon Plasma Thrusters M. Magarotto, M. Manente, P. de Calro, F. Trezzolani, D. Pavarin, and D. Melazzi CISAS, Padova, Italy
More informationDevelopment of a Hall Thruster Fully Kinetic Simulation Model Using Artificial Electron Mass
Development of a Hall Thruster Fully Kinetic Simulation Model Using Artificial Electron Mass IEPC-013-178 Presented at the 33rd International Electric Propulsion Conference, The George Washington University
More informationControl of gas discharge plasma properties utilizing nonlocal electron energy distribution functions DOE-OFS Plasma Science Center
1 Control of gas discharge plasma properties utilizing nonlocal electron energy distribution functions DOE-OFS Vladimir Demidov Current cooperation with: Mark Koepke (WVU), Igor Kaganovich (PPPL), Yevgeny
More informationThe Role of Secondary Electrons in Low Pressure RF Glow Discharge
WDS'05 Proceedings of Contributed Papers, Part II, 306 312, 2005. ISBN 80-86732-59-2 MATFYZPRESS The Role of Secondary Electrons in Low Pressure RF Glow Discharge O. Brzobohatý and D. Trunec Department
More informationCHARACTERIZATION OF A DC PLASMA WITH HOLLOW CATHODE EFFECT
Romanian Reports in Phisics, Vol. 56, No., P. 71-76, 004 CHARACTERIZATION OF A DC PLASMA WITH HOLLOW CATHODE EFFECT A. R. PETRE 1, M. BÃZÃVAN 1, V. COVLEA 1, V.V. COVLEA 1, ISABELLA IOANA OPREA, H. ANDREI
More informationNonlinear Diffusion in Magnetized Discharges. Francis F. Chen. Electrical Engineering Department
Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical Engineering Department PPG-1579 January, 1998 Revised April, 1998 Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical
More informationThe low-field density peak in helicon discharges
PHYSICS OF PLASMAS VOLUME 10, NUMBER 6 JUNE 2003 Francis F. Chen a) Electrical Engineering Department, University of California, Los Angeles, Los Angeles, California 90095-1597 Received 10 December 2002;
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationTheory of Gas Discharge
Boris M. Smirnov Theory of Gas Discharge Plasma l Springer Contents 1 Introduction 1 Part I Processes in Gas Discharge Plasma 2 Properties of Gas Discharge Plasma 13 2.1 Equilibria and Distributions of
More informationPlasma parameter evolution in a periodically pulsed ICP
Plasma parameter evolution in a periodically pulsed ICP V. Godyak and B. Alexandrovich OSRAM SYLVANIA, 71 Cherry Hill Drive, Beverly, MA 01915, USA The electron energy probability function (EEPF) has been
More informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
More informationarxiv:physics/ v1 [physics.plasm-ph] 5 Nov 2004
Ion Resonance Instability in the ELTRAP electron plasma G. Bettega, 1 F. Cavaliere, 2 M. Cavenago, 3 A. Illiberi, 1 R. Pozzoli, 1 and M. Romé 1 1 INFM Milano Università, INFN Sezione di Milano, Dipartimento
More informationCharacteristics and classification of plasmas
Characteristics and classification of plasmas PlasTEP trainings course and Summer school 2011 Warsaw/Szczecin Indrek Jõgi, University of Tartu Partfinanced by the European Union (European Regional Development
More informationSCALING OF HOLLOW CATHODE MAGNETRONS FOR METAL DEPOSITION a)
SCALING OF HOLLOW CATHODE MAGNETRONS FOR METAL DEPOSITION a) Gabriel Font b) Novellus Systems, Inc. San Jose, CA, 95134 USA and Mark J. Kushner Dept. of Electrical and Computer Engineering Urbana, IL,
More informationSheaths: More complicated than you think a
PHYSICS OF PLASMAS 12, 055502 2005 Sheaths: More complicated than you think a Noah Hershkowitz b University of Wisconsin-Madison, Madison, Wisconsin 53706 Received 7 December 2004; accepted 7 February
More informationFundamentals of Plasma Physics
Fundamentals of Plasma Physics Definition of Plasma: A gas with an ionized fraction (n i + + e ). Depending on density, E and B fields, there can be many regimes. Collisions and the Mean Free Path (mfp)
More informationBreakdown behavior in radio-frequency argon discharges
PHYSICS OF PLASMAS VOLUME 10, NUMBER 3 MARCH 2003 H. B. Smith, C. Charles, and R. W. Boswell Plasma Research Laboratory, Research School of Physical Sciences and Engineering, Australian National University,
More informationEffect of negative ions on the characteristics of plasma in a cylindrical discharge
Araghi and Dorranian Journal of Theoritical and Applied Physics 2013, 7:41 RESEARCH Open Access Effect of negative ions on the characteristics of plasma in a cylindrical discharge Farnaz Araghi and Davoud
More informationStatistical Instability of Barrier Micro-Discharges Operating in Townsend Regime
Statistical Instability of Barrier Micro-Discharges Operating in Townsend Regime V. P. Nagorny, V. N. Khudik Plasma Dynamics Corporation, Waterville, OH 43566 New kind of instability of a macroscopic physical
More informationBoundary Conditions for the Child Langmuir Sheath Model
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 6, DECEMBER 2000 2207 Boundary Conditions for the Child Langmuir Sheath Model Mikhail S. Benilov Abstract A collision-free space-charge sheath formed by
More informationSupporting Information for Conical Nanopores. for Efficient Ion Pumping and Desalination
Supporting Information for Conical Nanopores for Efficient Ion Pumping and Desalination Yu Zhang, and George C. Schatz,, Center for Bio-inspired Energy Science, Northwestern University, Chicago, Illinois
More informationNumerical Simulation of Trichel Pulses in a Negative Corona Discharge in Air
Proc. ESA Annual Meeting on Electrostatics, Paper K Numerical Simulation of Trichel Pulses in a Negative Corona Discharge in Air P. Sattari, G.S.P. Castle, and K. Adamiak Dept. of Electrical and Computer
More informationOpen Research Online The Open University s repository of research publications and other research outputs
Open Research Online The Open University s repository of research publications and other research outputs You cannot patch active plasma and collisionless sheath Journal Item How to cite: Franklin, R.
More informationIonization Detectors. Mostly Gaseous Detectors
Ionization Detectors Mostly Gaseous Detectors Introduction Ionization detectors were the first electrical devices developed for radiation detection During the first half of the century: 3 basic types of
More informationA Magnetohydrodynamic study af a inductive MHD generator
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan A Magnetohydrodynamic study af a inductive MHD generator Augusto Montisci, Roberto Pintus University of Cagliari, Department of Electrical
More informationarxiv: v1 [physics.plasm-ph] 10 Nov 2014
arxiv:1411.2464v1 [physics.plasm-ph] 10 Nov 2014 Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas A. Derzsi 1, I. Korolov
More informationNitrogen Glow Discharge by a DC Virtual Cathode
Nitrogen Glow Discharge by a DC Virtual Cathode Azza M. Shager a, Amany T. Sroor b, Hoda A. El Tayeb a, Hoda A. El Gamal a, and Mohamed M. Masoud a a Plasma Physics and Nuclear Fusion Department, Atomic
More informationPlasma Formation in the Near Anode Region in Hall Thrusters
41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit 10-13 July 2005, Tucson, Arizona AIAA 2005-4059 41 st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit AIAA-2005-4059 Plasma Formation
More informationIntegrated Particle Transport Simulation of NBI Plasmas in LHD )
Integrated Particle Transport Simulation of NBI Plasmas in LHD Akira SAKAI, Sadayoshi MURAKAMI, Hiroyuki YAMAGUCHI, Arimitsu WAKASA, Atsushi FUKUYAMA, Kenichi NAGAOKA 1, Hiroyuki TAKAHASHI 1, Hirohisa
More informationInfluence of surface charges on the structure of a dielectric barrier discharge in air at atmospheric pressure: experiment and modeling
Eur. Phys. J. Appl. Phys. 47, 22810 (2009) DOI: 10.1051/epjap/2009078 Regular Article THE EUROPEAN PHYSICAL JOURNAL APPLIED PHYSICS Influence of surface charges on the structure of a dielectric barrier
More informationA global (volume averaged) model of a chlorine discharge
A global (volume averaged) model of a chlorine discharge Eyþór Gísli Þorsteinsson and Jón Tómas Guðmundsson Science Institute, University of Iceland, Iceland Department of Electrical and Computer Engineering,
More informationThe Computational Simulation of the Positive Ion Propagation to Uneven Substrates
WDS' Proceedings of Contributed Papers, Part II, 5 9,. ISBN 978-8-778-85-9 MATFYZPRESS The Computational Simulation of the Positive Ion Propagation to Uneven Substrates V. Hrubý and R. Hrach Charles University,
More informationMulti-fluid Simulation Models for Inductively Coupled Plasma Sources
Multi-fluid Simulation Models for Inductively Coupled Plasma Sources Madhusudhan Kundrapu, Seth A. Veitzer, Peter H. Stoltz, Kristian R.C. Beckwith Tech-X Corporation, Boulder, CO, USA and Jonathan Smith
More information