Discrete unified gas kinetic scheme for electrostatic plasma and its comparison with the particle-in-cell method.

Abstract

In this paper, we present a finite-volume direct kinetic method, the so-called discrete unified gas kinetic scheme (DUGKS), for electrostatic plasma. One key feature of this method is the semi-implicit unsplitting treatment of particle transport and collision, and thus the time step in current DUGKS is not limited by the particle collision time. In addition, a fourth-order compact MUSCL scheme with a positivity preserving limiter is implemented in the interface reconstruction, which enables present DUGKS to preserve the favorable conservative property and positivity of distribution function. Combined with this MUSCL method, the semi-Lagrangian scheme is used for the particle transport in velocity space to remove Courant-Friedricks-Lewy restriction induced by the large electric force. As a result, the proposed DUGKS becomes an efficient and stable multiscale scheme. Several numerical experiments, including plasma sheath, linear Landau damping, collisional nonlinear Landau damping, and plasma ion acceleration, are performed to validate current DUGKS. A comparative study of current DUGKS with a general particle in cell (PIC) method which could handle particle collision in a conservative way is also presented. Theory and numerical experiments demonstrate that DUGKS is preferred for plasma flows involving small electrostatic perturbation and high collision regimes, while the PIC method is desired for the field- dominated plasma flows involving a wide range of velocities.