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.

Analysis of Dielectric Relaxation Effect in Fluid Modeling of AC PDP

A Fast Time Domain PMCHW Equations for Predicting Transient Responses of Dielectric Objects

A Stabilized and Well-Conditioned Marching-on-in Degree Method for Transient Analysis of Dielectric Objects

One-step leapfrog ADI-FDTD method for simulating electromagnetic wave propagation in general dispersive media.

A Gradually Implicit Method with Deferred Correction and Adaptive Relaxation Factors for Simulation to Nonlinear Heat Transfer.

A Fast Time Domain Electric Field Integral Equations for Predicting Transient Responses of Dielectric Objects

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

Discontinuous Galerkin Time Domain Electric Field Integral Equation Method For Solving Non-Conformal Structures

Studies on Fluid Model for Numerical Simulation of Gas Discharges in Color Plasma Displays

Influence of Parameters on Discharge Characteristics of Alternating Current Plasma Display Panel Cell

Improved HIE-FDTD Method for Simulating Thin Layer Structure in Linear Dispersive Media

Numerical Simulation of the Discharge Efficiency in Five-Electrode Ac Pdp

Simulation of Partial Discharge Sequences Using Fluid Equations

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

Numerical Analysis of the Discharge Characteristics of the Surface Charge AC-PDP

Direct Implicit and Explicit Energy-Conserving Particle-in-Cell Methods for Modeling of Capacitively-Coupled Plasma Devices

Energy stable finite element method for an electrohydrodynamic model with variable density

An Efficient Numerical Technique for the Simulation of Charge Transport in Polymeric Dielectrics

Novel Scheme for GPU‐Accelerated Finite‐Difference Time‐Domain Simulation of Electromagnetic Wave Interaction With Magnetic Plasma

Optimized Ade-Fdtd Method In Unmagnetized Plasma

Reduction of the grid effects in simulation plasmas