Abstract:A wide variety of plasma phenomena have been investigated during the past decades using the particle-in-cell/Monte Carlo collisions (PIC/MCC) method. As an important component of the PIC/MCC method, solving Poisson's equation is crucial for the accuracy and efficiency of calculations. Different acceleration techniques for solving finite difference discretization Poisson's equation are investigated and compared, including direct method, iterative method, multigrid (MG) method, parallel computing and inherited initial value. The charge density distribution with a known analytical solution is used to validate the algorithm and code. The optimal relaxation factor for the successive over-relaxation (SOR) method in 2D Poisson's equation with unequal grid node numbers in different dimensions is derived, which is only related to the dimension with the largest grid number. Although there will be a "more optimal" relaxation factor deviated from in some simulation cases, selecting the optimal relaxation factor derived always leads to a not slow solving speed. However, when SOR is used in MG for smoothing, the optimal relaxation factor will shift to 0.5-1.2 from the theoretical optimal value derived with the increase of MG levels. By comparing the convergence order under different relaxation factors and MG levels, the suitable MG level is proposed as log2[min(Nx, Ny)]-2. Combining the optimal SOR relaxation factor, MG, parallel computing and inherited initial values, the computational cost may decrease by 5 orders of magnitude than that by the simple Gaussian elimination (GE). Based on the optimal acceleration techniques mentioned above, a benchmark simulation case electron cyclotron drift instability (ECDI) in magnetized plasmas was run to further validate the developed PIC/MCC code. The distributions of electric field in the x-direction, electron density and electron temperature are all consistent with the literatures. This paper provides a reference for the acceleration strategy selection for solving Poisson's equation quickly in plasma simulations.

The PIC-DSMC numerical simulation for vacuum plume of MPT

Study on hybrid DSMC/EPSM method for chemical reacting gas flow

The study of hybrid DSMC/EPSM method of supersonic flat plate

Numerical investigation of plasma properties for the microwave discharge ion thruster μ10 using PIC-MCC simulation

Numerical Simulation of Plasma Plume with an Insulated Wall in Sheath Region Using a DSMC-PIC Parallel Algorithm

The Overall Numerical Simulation of Both High-Altitude Plume Flowfield and Inflow Region with DSMC Method

The Hypersonic Chemical Reaction Flow of Flat-Nosed Cylinder with Hybrid DSMC/EPSM Method

Particle-in-cell plus direct simulation Monte Carlo (PIC-DSMC) approach for self-consistent plasma-gas simulations

Numerical study on wave attenuation via 1D fully kinetic electromagnetic particle-in-cell simulations

Self consistent kinetic simulations of SPT and HEMP thrusters including the near-field plume region

Kinetic, 3-D, PIC-DSMC Simulations of Ion Thruster Plumes and the Backflow Region

Numerical Analysis of a Miniature Microwave-discharge Ion Thruster Using Water as the Propellant

Coupling of Fluid and Particle-in-Cell Simulations of Ambipolar Plasma Thrusters

Study of scaling law for particle-in-cell/Monte Carlo simulation of low-temperature magnetized plasma for electric propulsion

Numerical Research on Plasma Generation and Expansion into Vacuum

Three-Dimensional Pic-Mc Modeling for Relativistic Electron Beam Transport Through Dense Plasma

A Comparative Study on the Acceleration Techniques for Solving Finite Difference Discretization Poisson's Equation in the PIC/MCC Method

Three-Dimensional Direct-Implicit Particle-in-Cell Model Using Trilinear Anisotropic Immersed-Finite-Element for Plasma Propulsion

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

Monte Carlo Simulation for Aerodynamic Coefficients of Satellites in Low-Earth Orbit

PIFE-PIC: Parallel Immersed-Finite-Element Particle-In-Cell For 3-D Kinetic Simulations of Plasma-Material Interactions