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.

One of the main factors that limit the lifetime of the Hall effect Thrusters (HETs) is the erosion of the acceleration channel caused by the flux of energetic ions. The magnetic field that is curved and convex towards the anode has been widely used in HETs because of its role in reducing the divergence angle of the ion beam and the channel wall erosion. However, the mechanism of the influence of the magnetic field line curvature on the wall erosion is not clear. Therefore, in this paper, a 2D3V numerical model based on the immersed-finite-element and particle-in-cell (IFE-PIC) method is established to simulate the radial-azimuthal plane near the exit of the Hall thruster. The effect of the tilt angle of the magnetic field line on the wall sputtering erosion rate is analyzed. The results show that compared to the case with the electric field E perpendicular to the magnetic field B, the energy of the ions hitting the channel wall is smaller and the wall erosion is weaker when the magnetic field lines are convex to the anode. As the tilt angle of the magnetic field lines increases from 0° to 60°, the erosion rate is reduced by 90%. Conversely, when the magnetic field lines are convex to the exit plane of the channel, the wall erosion is much more serious compared to the case with the orthogonal electric field E and the magnetic field B. As the tilt angle of the magnetic field line changes from 0° to 60°, the erosion rate is enhanced by 171%. The results in this paper are instructive for the design and optimization of the magnetic field of the HETs.