Bandpass approximate Crank-Nicolson implementation for anisotropic gyrotropic plasma open region simulation

For overcoming the Courant-Friedrich-Levy (CFL) stable condition and maintaining the stability of the algorithm, unconditionally stable implementation is proposed for the simulation of anisotropic gyrotropic plasma in bandpass open region problems. To be more accurately, the proposed implementation is based on the complex envelope method, higher order concept, perfectly matched layer formulation and Crank-Nicolson Direct-Splitting procedure which shows advantages over them including accuracy, efficiency and absorption. The gyrotropic plasma with unique anisotropic characteristic can be simulated by the modified approach. For the further illustration of effectiveness, several numerical examples are introduced. Through the results and comparison with previous work, it can be concluded that the entire performance which includes accuracy, efficiency and absorption can be significantly improved. In addition, the mesh size can be no longer limited by the CFL condition rather than the computational accuracy which shows considerable potential in numerical simulation of fine structures. Meanwhile, the time step can be obtained according to the maximum frequency rather than bandwidth resulting in the exponential decrement of running time.