Computationally Efficient CN-PML for EM Simulations

• Published in: IEEE transactions on microwave theory and techniques 2019-12, Vol.67 (12), p.4646-4655
• Main Authors: Jiang, Hao Lin, Wu, Li Ting, Zhang, Xin Ge, Wang, Qiang, Wu, Pei Yu, Liu, Che, Cui, Tie Jun
• Format: Article
• Language: English
• Subjects: Bilinear transform (BT); Computational efficiency; Computational modeling; Computer simulation; Crank–Nicolson (CN); Finite difference methods; Finite difference time domain method; finite-difference time domain (FDTD); Iterative methods; Lattices (mathematics); Numerical stability; perfectly matched layer (PML); Perfectly matched layers; Run time (computers); Stability; Stability criteria; Time-domain analysis
• ISSN: 0018-9480; 1557-9670
• DOI: 10.1109/TMTT.2019.2946160

Since nearly all studies concerning the approximate Crank-Nicolson perfectly matched layer (CN-PML) are limited to 2-D cases, a computationally efficient implementation that can be used to truncate 3-D finite-difference time-domain (FDTD) lattices is presented in this article. More precisely, it is based on the CN direct-splitting (DS) scheme and the bilinear transform (BT) method. This article can fully exploit the unconditional stability of the standard CN-FDTD method and can be free from the Courant-Friedrich-Lewy (CFL) limit; hence, it is especially suitable for situations where space discretization step is much smaller than 1/10th or 1/20th of the smallest wavelength of interest. Aiming at further reducing the requirement of the computer resources, this new implementation can be reformulated in more simple forms if proper auxiliary variables are introduced. It therefore shows a higher iteration speed than other published unconditionally stable PMLs as fewer numbers of arithmetic operations are involved. Finally, three numerical examples, including scatting, transmission, and radiation, are also provided to validate its running time, unconditional stability, and absorption characteristic.