Computationally Efficient Implicit ADI Theory and Its Open Region Problems for Anisotropic Gyrotropic Plasma Simulations

Abstract

With its unique anisotropic characteristic, gyrotropic plasma is extensively employed in a series of applications. As it is hard to obtain the electric behavior of anisotropic gyrotropic plasma directly, numerical simulation is regarded as one of the most important ways on the investigation. Based upon higher order formulation and an alternating‐direction‐implicit (ADI) procedure, the unconditionally stable ADI perfectly matched layer (PML) is proposed in finite‐difference time‐domain lattice to solving the open region problems for anisotropic gyrotropic plasma. To decouple the updated equations in the gyrotropic plasma regions, approximate decoupling theory is introduced during the implementations. The proposed implementation takes advantages of the higher order formulation, ADI procedure, approximate decoupling theory, and PML implementation in terms of considerable absorption, efficiency, and accuracy in open region problems. Numerical examples are carried out for the demonstration of effectiveness both in theory and experimentally. Through the results, it can be concluded that the proposed implementation is effective and accurate for terminating anisotropic gyrotropic plasma in open region problems when time step surpasses far beyond the Courant–Friedrich–Levy condition.

Efficient Enhanced Hybrid Implicit-Explicit Procedure to Gyrotropic Plasma in Open Regions With Fine Geometry Details Along Single Direction

“…Until now, several methods have been developed for the investigation of anisotropic gyrotropic plasma including the weighted laguerre polynomials, ADI, LOD, CN procedures [41][42][43][44][46][47]. All of these procedures belong to the unconditionally stable algorithms which shows significant decrement of efficiency and accuracy when simulating the fine geometry details along a single direction.…”

“…It should be noticed that they can merely solve the Maxwell's equations in two-dimensions. Recently, the ADI and CN procedures are combined with the gyrotropic plasma simulation for the three-dimensions [43][44]. Although all of them can overcome the CFL condition, the limitations with fine geometry details along single direction structures still need to be perfected.…”