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.

Bandpass Approximate Crank-Nicolson Implementation for Anisotropic Gyrotropic Plasma Open Region Simulation

One-step Leapfrog ADI-FDTD Method for Anisotropic Magnetized Plasma

Finite-Element-Based Generalized Impedance Boundary Condition for Modeling Plasmonic Nanostructures

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Numerical Simulation of Heat Transfer and Fluid Flow of a Non-Equilibrium Argon Plasma Jet with Confined Wall

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

ICN‐FDTD scheme with absorption boundary condition for nonuniform rotational symmetric geometrics

A fully decoupled iterative method with three-dimensional anisotropic immersed finite elements for Kaufman-type discharge problems

An Unstructured Body-of-Revolution Electromagnetic Particle-in-Cell Algorithm with Radial Perfectly Matched Layers and Dual Polarizations

Multiscale simulation of neutral particle flows in the plasma edge

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

Analysis of electromagnetic wave propagation and scattering characteristics of plasma shealth via high order ADE-ADI FDTD

Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction

ANSYS HFSS as a new numerical tool to study wave propagation inside anisotropic magnetized plasmas in the Ion Cylotron Range of Frequencies

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

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

High-efficient Particle-in-cell/monte Carlo Model for Complex Solution Domain Andsimulation of Anode Layer Ion Source

Performance improvement of particle-in-cell method for numerical modelling of open magnetic system

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

High-order Targeted Essentially Non-Oscillatory Scheme for Two-Fluid Plasma Model

Higher-order moment convergent method in weakly anisotropic plasma and the NLVFP code for solution of the 0D-2V Vlasov-Fokker-Planck equation