Modeling low pressure collisional plasma sheath with space-charge effect

The present work develops a computationally efficient one-dimensional subgrid embedded finite element formulation for plasma-sheath dynamics. The model incorporates space-charge effect throughout the whole plasma and the sheath region using multifluid equations. Secondary electron emission is not considered. A third-order temperature dependent polynomial is used to self-consistently calculate the rate of ionization in the plasma dynamic equations. The applications include dc and rf sheath inside a glow discharge tube where the noble gas is immobile, and a partially ionized plasma sheath inside an electric propulsion thruster channel in which the gas flows. The electron and ion number densities of the numerical solution decrease in the sheath region as expected. The ion velocity and electron temperature profiles also exhibit the expected behavior. The computed sheath potential compares well with the available experimental data.

[1]  J. Boeuf,et al.  Two-dimensional model of a capacitively coupled rf discharge and comparisons with experiments in the Gaseous Electronics Conference reference reactor. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[2]  A. J. Baker,et al.  NONLINEAR SUBGRID EMBEDDED FINITE-ELEMENT BASIS FOR STEADY MONOTONE CFD SOLUTIONS, PART II: BENCHMARK NAVIER-STOKES SOLUTIONS , 1998 .

[3]  W. Lotz Electron Impact Ionization Cross Sections and Ionization Rate Coefficients for Atoms and Ions , 1967 .

[4]  Sternberg,et al.  Dynamic model of the electrode sheaths in symmetrically driven rf discharges. , 1990, Physical review. A, Atomic, molecular, and optical physics.

[5]  I. Langmuir,et al.  CONTACT POTENTIAL MEASUREMENTS WITH ADSORBED FILMS , 1929 .

[6]  B. P. Pandey,et al.  Numerical investigation of a Hall thruster plasma , 2002 .

[7]  Joseph Shang,et al.  RECENT RESEARCH IN MAGNETO-AERODYNAMICS , 2001 .

[8]  D. Graves,et al.  Matching an RF sheath model to a bulk plasma model , 1995 .

[9]  Valery Godyak,et al.  Approximation of the bounded plasma problem by the plasma and the sheath models , 1996 .

[10]  V. I. Lagutin,et al.  Aerodynamic drag reduction by plasma and hot-gas injection , 2000 .

[11]  B. P. Pandey,et al.  Plasma–wall interaction inside a Hall thruster , 2002, Journal of Plasma Physics.

[12]  D. Wolff,et al.  Modelling of sheaths, presheaths, current-voltage characteristics and double probes in low-pressure plasmas in coaxial and concentric geometries , 1995 .

[13]  Lawrence J. Overzet,et al.  Comparison of electron‐density measurements made using a Langmuir probe and microwave interferometer in the Gaseous Electronics Conference reference reactor , 1993 .

[14]  H. Valentini Bohm criterion for the collisional sheath , 1996 .

[15]  W. Lotz,et al.  Electron-impact ionization cross-sections and ionization rate coefficients for atoms and ions from hydrogen to calcium , 1968 .

[16]  Self-consistent modelling of helium discharges: investigation of the role of He2+ ions , 2001 .

[17]  A. W. Trivelpiece,et al.  Introduction to Plasma Physics , 1976 .

[18]  A. M. Bishaev,et al.  Local plasma properties in a Hall-current accelerator with an extended acceleration zone , 1978 .

[19]  Iu. P. Raizer Gas Discharge Physics , 1991 .

[20]  M. Mitchner,et al.  Partially ionized gases , 1973 .

[21]  A. Lichtenberg,et al.  Principles of Plasma Discharges and Materials Processing , 1994 .

[22]  E. Ahedo,et al.  Model of the plasma discharge in a Hall thruster with heat conduction , 2002 .

[23]  N. Hershkowitz,et al.  First experimental measurements of the plasma potential throughout the presheath and sheath at a boundary in a weakly collisional plasma. , 2002, Physical review letters.