Author:
Xie Guo-Da,Hou Gui-Lin,Niu Kai-Kun,Feng Nai-Xing,Fang Ming,Li Ying-Song,Huang Zhi-Xiang, ,
Abstract
Dispersive media refer to a class of natural substances, including living organisms, composite materials, plasma and water, with diverse applications in areas such as biomedicine, microwave sensing, electromagnetic protection, and stealth technology. In the pursuit of investigating the electromagnetic properties of these media, time-domain numerical methods, including finite difference in time domain (FDTD), finite element method, and time domain boundary integral equation method, have been widely utilized. Time-domain numerical methods are preferred to their frequency-domain counterparts owing to their ability to handle nonlinear and wideband problems, as well as various material properties. The FDTD method, in particular, is a highly adaptable, robust, and easy-to-use numerical method that directly solves the Maxwell equations while also simulating the reflection, transmission, and scattering of electromagnetic waves in complex dispersion media. Nonetheless, the traditional FDTD method suffers low computational efficiency arising from the Courant-Friedrichs-Lewy (CFL) stability condition. To solve the problem of low computational efficiency, a new method, the complying divergence implicit finite-difference time-domain (CDI-FDTD) method with a one-step leapfrog scheme, is introduced for lossy Debye dispersive media. The Maxwell equations in the frequency domain form a starting point, and the Fourier transform is utilized to transform the electromagnetic field components from the frequency domain to the time domain. To approximate the integral terms arising from the frequency-to-time domain transformation, a recursive integration (RI) method is employed. Subsequently, the time-domain Maxwell equations and auxiliary variables are discretized with a one-step leapfrog implicit scheme. The iterative formula of the RI-CDI-FDTD algorithm for lossy Debye dispersive media is then derived. The RI-CDI-FDTD method does not change the formulas of the traditional CDI-FDTD method while only requiring to add auxiliary variables for updating field components to the dispersive medium region. The numerical implementation is straightforward, and the electromagnetic modeling is flexible. Moreover, the unconditional stability of the RI-CDI-FDTD algorithm is proven by using the von Neumann method. Finally, some numerical examples are presented to demonstrate the effectiveness and efficiency of the proposed method. In conclusion, our work contributes a crucial numerical simulation tool to accurately modeling complex dispersive media while providing a systemic stability analysis method for time-domain numerical methods.
Publisher
Acta Physica Sinica, Chinese Physical Society and Institute of Physics, Chinese Academy of Sciences
Subject
General Physics and Astronomy
Reference30 articles.
1. Cheng X, Shao W, Wang K, Wang B Z 2019 IEEE Antennas Wirel. Propag. Lett. 18 1931
2. Chakarothai J 2018 IEEE Trans. Antennas Propag. 67 6076
3. Fang Y, Liu J F, Jiao Z H, Bai G H, Xi X L 2018 IEEE Trans. Plasma Sci. 47 173
4. Tian H M, He S M, Wang M Y, Li G P, Yu M X, Zhang S T, Xu J 2021 IEEE Antennas Wirel. Propag. Lett. 20 2392
5. Chen B Y, Zhang Y R, Wang L, Wang F F 2016 Acta Phys. Sin. 65 144101
陈碧云, 张业荣, 王磊, 王芳芳 2016 物理学报 65 144101
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献