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.

Parallelization And Optimization Of Electrostatic Particle-In-Cell/Monte-Carlo Coupled Codes As Applied To Rf Discharges

Novel Parallel Meshing Technique Implemented into Parallel FDTD Simulation Using Function Language for Solving Electrically Large and Complex EMC Problems

Massively Parallel Simulation Of Large-Scale Electromagnetic Problems Using One High-Performance Computing Scheme And Domain Decomposition Method

A Scalable HPC-Based Domain Decomposition Method for Multiphysics Modeling of RF Devices

Performance Analysis and Application of 2-D Plasma Parallel Particle-in-cell Simulations

Three-dimensional Electromagnetic Particle-in-cell Simulation by Parallel Computing

On the Energy Conservation Electrostatic Particle-In-Cell/Monte Carlo Simulation: Benchmark and Application to the Radio Frequency Discharges

Implicit and electrostatic Particle-in-cell/Monte Carlo model in two dimensional and axisymmetric geometry I: analysis of numerical techniques

Parallel 3D electromagnetic PIC simulation based on OpenMP

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

Three-dimensional Parallelized Fully-Electromagnetic Particle Simulation Code UNIPIC-3D

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

Exploiting Parallelism on Shared Memory in the QED Particle-in-Cell Code PICADOR with Greedy Load Balancing

The Plasma Simulation Code: A Modern Particle-in-cell Code with Patch-Based Load-Balancing.

Parallelization and I/O Optimization for a 3-D Plasma Simulation Code

Implicit electrostatic particle-in-cell/Monte Carlo simulation for the magnetized plasma: Algorithms and application in gas-inductive breakdown*

Implementation of parallel plasma particle-in-cell codes on PC cluster

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

VSHPIC: A Particle-In-Cell Algorithm Based On Vector Spherical Harmonics Expansion

Mitigation of Numerical Cerenkov Radiation and Instability Using a Hybrid Finite Difference-Fft Maxwell Solver and a Local Charge Conserving Current Deposit