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.

Corrections on radial edge errors of charge and current density deposition for two-dimension radial-axial particle-in-cell simulations

Correction of Particle-Moving Algorithm for Particle-in-cell Simulation of Plasma

A Volume-Weighting Cloud-In-Cell Model for Particle Simulation of Axially Symmetric Plasmas

A VOLUME-WEIGHTING CIC METHOD FOR NUMERICALSIMULATION OF AXIAL SYMMETRICAL PLASMAS

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

Geometric Electrostatic Particle-In-Cell Algorithm on Unstructured Meshes

Particle-Moving Algorithm for Particle-in-Cell Simulation of Plasmas in Polar Coordinates System

Exploring the Microscopic Physical Processes of Z-pinch by a Modified Electrostatic Direct Implicit Particle-in-Cell Algorithm

On a 1D-electrostatic test problem for the PIC method

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

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

Note on particle balance in Particle-In-Cell/Monte Carlo model and its implications on the steady-state simulation

A generalized external circuit model for electrostatic particle-in-cell simulations

EZ: An Efficient, Charge Conserving Current Deposition Algorithm for Electromagnetic Particle-In-Cell Simulations

3D PIC-MCC Simulation of Corona Discharge in Needle-Plate Electrode with External Circuit

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

Monte Carlo Approach to Calculate Ionization Dynamics of Hot Solid-Density Plasmas Within Particle-In-Cell Simulations

An Improved Particle Distribution and Correction Procedure for Particle Level Set Method

Monte-Carlo approach to calculate the ionization of warm dense matter within particle-in-cell simulations

An Exactly Energy-conserving Electromagnetic Particle-in-cell Method in Curvilinear Coordinates