A finite element procedure for time-dependent radio-frequency sheaths based on a two-dimensional microscale fluid model

• Published in: Computer physics communications Vol. 291; no. C; p. 108841
• Main Authors: Kohno, H., Myra, J.R.
• Format: Journal Article
• Language: English
• Published: Netherlands Elsevier B.V 01.10.2023; Elsevier
• Subjects: Admittance; Finite element method; Magnetic confinement fusion; Plasma heating; RF sheath; Finite element method; RF sheath; Admittance; Plasma heating; Magnetic confinement fusion
• DOI: 10.1016/j.cpc.2023.108841
• ISSN: 0010-4655

In this paper, we present a finite element scheme for the analysis of time-dependent radio-frequency (RF) sheath behavior in a plasma-filled domain bounded by periodically curved plates. This numerical scheme is based on a two-dimensional (2D) microscale model in which the time-dependent cold-ion fluid equations, the Maxwell–Boltzmann relation for electrons, and Poisson's equation are solved subject to periodic boundary conditions (BCs) and conducting-wall BCs. The continuity of the total current is employed in localized regions to provide a constraint on the reference sheath potential. The primary purpose of this work is to understand 2D dynamic sheath behavior in order to improve predictive capabilities for RF wave interactions in magnetic fusion experiments. In particular, this work treats cases where the local radius of curvature of the wall surface is comparable to the non-neutral sheath width. Using the developed numerical code, the dependences of the ion, electron, and displacement admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle are investigated. It is shown that the ion and electron admittances are nearly unchanged over wide ranges of the bump height and ion magnetization. In addition, the 2D sheath effects are assessed through the spatial distributions of various quantities such as the electron density, electrostatic potential, ion velocity, and surface wall current. •New finite element scheme based on a 2D microscale radio-frequency sheath model.•Analysis of dynamic behaviors of non-neutral sheaths and magnetic presheaths.•Use of a field aligned mesh and flux tube segments for the continuity of current.•Dependence of the surface-integrated admittances on various parameters.