Kinetic simulations of low-pressure inductively coupled plasma: an implicit electromagnetic PIC/MCC model with the ADI-FDTD method

Abstract

Low-pressure inductively coupled plasma (ICP) is promising for space electric propulsion. For the first time, an implicit electromagnetic particle-in-cell Monte Carlo collision model based on the alternating-direction-implicit finite-difference time-domain (ADI-FDTD) method is developed to investigate low-pressure xenon plasma characteristics of a miniature ICP source. The induced simulated electric field is well consistent with that calculated by the finite element method, indicating that this method can provide an accurate estimation of the electromagnetic field. The simulation time step used in the ADI-FDTD method is no longer restricted by the Courant–Friedrichs–Lewy constraints. Compared with the FDTD method, the ADI-FDTD method increases the size of the time step and significantly improves computational efficiency. The method is validated by comparing the simulated and measured electron density and plasma potential profile and reasonable agreement is reached. Therefore, the model is used to investigate the temporal and spatial distribution of plasma properties and the influence of the current amplitude of radio frequency (RF) coil, applied frequency of RF coil and neutral gas pressure on the plasma dynamics in the ionization chamber of a miniature gridded RF ion thruster. To explain the influence of the operating parameters, a concept called ‘the energy relaxation characteristics of electrons in response to the change of electric field’ is proposed and verified. The simulations also find that the oscillation frequency of plasma properties is twice the applied frequency of RF coil. The oscillation characteristics reveal the dynamic energy balance in the ICP. The experiment on the gridded RF ion thruster BHRIT-4 confirms the oscillation by measuring the plasma sheath potential.