Approximate Crank–Nicolson Algorithm with Higher-Order PML Implementation for Plasma Simulation in Open Region Problems

Abstract

By incorporating the higher-order concept with the perfectly matched later (PML) scheme, unconditionally stable approximate Crank–Nicolson algorithm is proposed for plasma simulation in open region problems. More precisely, the proposed implementation is based on the CN Direct-Splitting (CNDS) procedure for the finite-difference time-domain (FDTD) unmagnetized plasma simulation. The unmagnetized plasma can be regarded as frequency-dependent media which can be calculated by the piecewise linear recursive convolution (PLRC) method. The proposed implementation shows the advantages of higher-order concept, CNDS procedure, and PLRC method in terms of improved absorbing performance, enhanced computational efficiency, and outstanding calculation accuracy. Numerical examples are introduced to indicate the effectiveness and efficiency. It can be concluded from results that the proposed scheme shows considerable efficiency, accuracy, absorption, and unconditional stability.