-
混合显-隐式时域有限差分方法的时间步长只由空间两个方向上的网格大小决定, 因而该方法被广泛应用于沿一个方向具有精细结构的电磁目标的模拟. 本文通过对传统时域有限差分方法基本公式进行近似, 给出适用于线性、非色散空间, 一维精细结构电磁目标模拟的通用混合显-隐式时域有限差分方法. 在通用混合显-隐式时域有限差分方法中, 通过媒质编号索引和迭代系数空间定位实现三对角矩阵的自动构建, 引入等效参数方法解决麦克斯韦方程微分形式在不同介质分界面失效的问题, 引入卷积完全匹配层边界条件对计算区域进行截断. 此外, 本文提出了一种缩减三对角矩阵数量的方法, 该方法可有效降低程序运行所需内存, 提高程序运行效率. 应用通用混合显-隐式时域有限差分方法仿真了平面波照射介质板和双频微带倒F天线两个模型, 数值计算结果与传统时域有限差分方法和CST软件的计算结果一致, 但与传统时域有限差分方法相比, 通用混合显-隐式时域有限差分方法的计算效率大大提高.
-
关键词:
- 电磁目标 /
- 混合显-隐式时域有限差分 /
- 等效参数 /
- 卷积完全匹配层
The hybrid implicit-explicit finite-difference time-domain (HIE-FDTD) method is a weakly conditionally stable finite-difference time-domain (FDTD) method. The time step size of the HIE-FDTD method is only confined by the two coarse spatial cell sizes, so it is widely used in the simulation of electromagnetic targets with fine structures along one direction. In this work, the basic iterative formulations of the HIE-FDTD method are proposed by approximating the formulations of the FDTD method. In these formulations, the iterative coefficients are marked with spatial grid numbers. Therefore, the coefficients can be calculated automatically according to the medium parameters of the spatial cells by indexing the grid numbers. Since the triangular matrix which is used to calculate the electric field is based on the iterative coefficients, the triangular matrix can also be updated automatically. In addition, a method to reduce the number of tridiagonal matrices is proposed in this work, which can effectively reduce the calculation memory and improve the calculation efficiency. In the proposed HIE-FDTD method, equivalent parameters are employed at the interfaces of different media and the convolution perfectly matched layer boundary condition is used to truncate the computational region. Based on the proposed HIE-FDTD method, a series of programs are implemented, which can simulate arbitrary electromagnetic targets with fine structure in one dimension in linear and non-dispersive space. A dielectric plate irradiated by planar wave and a dual-frequency microstrip inverted F antenna are simulated by using these programs. The numerical results are in good agreement with those from the traditional FDTD method and CST software, and the computational efficiency of the proposed HIE-FDTD method is greatly improved in comparison with that of the traditional FDTD method. This study provides a reliable simulation tool for the wide application of the HIE-FDTD method.-
Keywords:
- electromagnetic target /
- hybrid implicit-explicit finite-difference time-domain /
- equivalent parameter /
- convolution perfectly matched layer
[1] Yee K S 1966 IEEE Trans. Antennas Propag. 14 302
Google Scholar
[2] Taflove A, Hagness S C 2005 Computational Electrodynamics: The Finite-Difference Time-Domain Method (3rd Ed.) (Norwood: Artech House) pp107–167
[3] Kim Y H, Choi H, Cho J, Jung K Y 2020 J. Electr. Eng. Technol. 15 1281
Google Scholar
[4] Zhang Y X, Feng N X, Wang L X, Guan Z, Liu Q H 2020 IEEE Trans. Antennas Propag. 68 366
Google Scholar
[5] Li B H, Du Y P, Chen M L, Li Z 2021 IEEE Trans. Electromagn. Compat. 63 762
Google Scholar
[6] Kim Y J, Jung K Y 2021 IEEE Trans. Antennas Propag. 69 6600
Google Scholar
[7] 杨利霞, 刘超, 李清亮, 闫玉波 2022 物理学报 71 064101
Google Scholar
Yang L X, Liu C, Li Q L, Yan Y B 2022 Acta Phys. Sin. 71 064101
Google Scholar
[8] 叶志红, 张杰, 周健健, 苟丹 2020 物理学报 69 060701
Google Scholar
Ye Z H, Zhang J, Zhou J J, Gou D 2020 Acta Phys. Sin. 69 060701
Google Scholar
[9] 方云 2020 博士学位论文 (西安: 西安理工大学)
Fang Y 2020 Ph. D. Dissertation (Xi’an: Xi’an University of Technology) (in Chinese)
[10] Huang B K, Wang G, Jiang Y S, Wang W B 2003 Microwave Opt. Technol. Lett. 39 97
Google Scholar
[11] Chen J, Wang J G 2006 Microwave Opt. Technol. Lett. 48 2291
Google Scholar
[12] Chen J, Wang J G 2007 Microwave Opt. Technol. Lett. 49 1001
Google Scholar
[13] Chen J, Xu N, Zhang A X, Guo J Y 2016 IEEE Trans. Antennas Propag. 64 3011
Google Scholar
[14] Zhai M L, Peng H L, Li D M, Wang D W, Xie H, Yin W Y 2017 Int. J. Numer. Modell. 30 e2250
Google Scholar
[15] Xu N, Chen J, Wang J G, Qin X J, Shi J P 2017 IET Microwaves Antennas Propag. 11 92
Google Scholar
[16] Chen J, Guo J Y, Tian C M 2018 IEEE Trans. Electromagn. Compat. 60 362
Google Scholar
[17] Zhang K L, Wang L, Wang M J, Zheng H X, Li E 2020 Int. J. RF Microwave Comput. Aided Eng. 30 e22151
Google Scholar
[18] Chen J, Wang J G 2007 IEEE Trans. Antennas Propag. 55 3603
Google Scholar
[19] Cui Y L, Chen B, Xiong R, Mao Y F 2013 J. Electromagn. Waves Appl. 27 529
Google Scholar
[20] Unno M, Asai H 2011 IEEE Microwave Wireless Compon. Lett. 21 453
Google Scholar
[21] Ge D B, Yan Y B 2005 Finite-Difference Time-Domain Method for Electromagnetic Waves (3rd Ed.) (Xian: Xidian University Press) pp9–27 (in Chinese) [葛德彪, 闫玉波 2005 电磁波时域有限差分方法 (第三版) (西安: 西安电子科技大学出版社) 第9–27页
[22] Atef Z, Elsherbeni, Demir V 2015 The Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations (2nd Ed.) (Scotland: SciTech Publishing) pp43–48
[23] Chen J, Wang J G, Xu N 2016 Weakly Conditionally Stable Finite-Difference Time-Domain Method (Beijing: Science Press) pp63–69 (in Chinese) [陈娟, 王建国, 许宁 2016 弱条件稳定时域有限差分方法 (北京: 科学出版社) 第63—69页]
[24] Chen J, Wang J G 2013 Appl. Comput. Electromagnet. 28 680
[25] Roden J A, Gedney S D 2000 Microwave Opt. Technol. Lett. 27 334
Google Scholar
[26] 王秀芝, 高劲松, 徐念喜 2013 物理学报 62 237302
Google Scholar
Wang X Z, Gao J S, Xu N X 2013 Acta Phys. Sin. 62 237302
Google Scholar
-
图 5 观察点处Ex的时域波形 (a) 观察点1; (b) 观察点2
Fig. 5. Time domain waveforms of Ex at the observation point: (a) Observation point one; (b) observation point two.
图 6 双频微带倒F天线结构图
Fig. 6. Structure of the dual-frequency microstrip inverted-F antenna.
图 8 天线在2.43 GHz的辐射方向图 (a) XOY面; (b) XOZ面
Fig. 8. Radiation patterns of the inverted-F antenna at 2.43 GHz: (a) XOY plane; (b) XOZ plane.
图 9 天线在4.96 GHz的辐射方向图 (a) XOY面; (b) XOZ面
Fig. 9. Radiation patterns of the inverted-F antenna at 4.96 GHz: (a) XOY plane; (b) XOZ plane.
表 1 直接法和缩减法对比结果
Table 1. Comparison results of the direct method and the reduction method.
计算方法 直接法 缩减法 1×1 5×5 10×10 1×1 5×5 10×10 三对角矩阵数量 684 11100 41760 12 12 12 程序运行时间/s 0.35 84.78 1299.4 0.02 0.20 0.51 程序所需内存/
MB1.40 24.14 90.80 0.98 15.77 59.27 
下载: 导出CSV
-
[1] Yee K S 1966 IEEE Trans. Antennas Propag. 14 302
Google Scholar
[2] Taflove A, Hagness S C 2005 Computational Electrodynamics: The Finite-Difference Time-Domain Method (3rd Ed.) (Norwood: Artech House) pp107–167
[3] Kim Y H, Choi H, Cho J, Jung K Y 2020 J. Electr. Eng. Technol. 15 1281
Google Scholar
[4] Zhang Y X, Feng N X, Wang L X, Guan Z, Liu Q H 2020 IEEE Trans. Antennas Propag. 68 366
Google Scholar
[5] Li B H, Du Y P, Chen M L, Li Z 2021 IEEE Trans. Electromagn. Compat. 63 762
Google Scholar
[6] Kim Y J, Jung K Y 2021 IEEE Trans. Antennas Propag. 69 6600
Google Scholar
[7] 杨利霞, 刘超, 李清亮, 闫玉波 2022 物理学报 71 064101
Google Scholar
Yang L X, Liu C, Li Q L, Yan Y B 2022 Acta Phys. Sin. 71 064101
Google Scholar
[8] 叶志红, 张杰, 周健健, 苟丹 2020 物理学报 69 060701
Google Scholar
Ye Z H, Zhang J, Zhou J J, Gou D 2020 Acta Phys. Sin. 69 060701
Google Scholar
[9] 方云 2020 博士学位论文 (西安: 西安理工大学)
Fang Y 2020 Ph. D. Dissertation (Xi’an: Xi’an University of Technology) (in Chinese)
[10] Huang B K, Wang G, Jiang Y S, Wang W B 2003 Microwave Opt. Technol. Lett. 39 97
Google Scholar
[11] Chen J, Wang J G 2006 Microwave Opt. Technol. Lett. 48 2291
Google Scholar
[12] Chen J, Wang J G 2007 Microwave Opt. Technol. Lett. 49 1001
Google Scholar
[13] Chen J, Xu N, Zhang A X, Guo J Y 2016 IEEE Trans. Antennas Propag. 64 3011
Google Scholar
[14] Zhai M L, Peng H L, Li D M, Wang D W, Xie H, Yin W Y 2017 Int. J. Numer. Modell. 30 e2250
Google Scholar
[15] Xu N, Chen J, Wang J G, Qin X J, Shi J P 2017 IET Microwaves Antennas Propag. 11 92
Google Scholar
[16] Chen J, Guo J Y, Tian C M 2018 IEEE Trans. Electromagn. Compat. 60 362
Google Scholar
[17] Zhang K L, Wang L, Wang M J, Zheng H X, Li E 2020 Int. J. RF Microwave Comput. Aided Eng. 30 e22151
Google Scholar
[18] Chen J, Wang J G 2007 IEEE Trans. Antennas Propag. 55 3603
Google Scholar
[19] Cui Y L, Chen B, Xiong R, Mao Y F 2013 J. Electromagn. Waves Appl. 27 529
Google Scholar
[20] Unno M, Asai H 2011 IEEE Microwave Wireless Compon. Lett. 21 453
Google Scholar
[21] Ge D B, Yan Y B 2005 Finite-Difference Time-Domain Method for Electromagnetic Waves (3rd Ed.) (Xian: Xidian University Press) pp9–27 (in Chinese) [葛德彪, 闫玉波 2005 电磁波时域有限差分方法 (第三版) (西安: 西安电子科技大学出版社) 第9–27页
[22] Atef Z, Elsherbeni, Demir V 2015 The Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations (2nd Ed.) (Scotland: SciTech Publishing) pp43–48
[23] Chen J, Wang J G, Xu N 2016 Weakly Conditionally Stable Finite-Difference Time-Domain Method (Beijing: Science Press) pp63–69 (in Chinese) [陈娟, 王建国, 许宁 2016 弱条件稳定时域有限差分方法 (北京: 科学出版社) 第63—69页]
[24] Chen J, Wang J G 2013 Appl. Comput. Electromagnet. 28 680
[25] Roden J A, Gedney S D 2000 Microwave Opt. Technol. Lett. 27 334
Google Scholar
[26] 王秀芝, 高劲松, 徐念喜 2013 物理学报 62 237302
Google Scholar
Wang X Z, Gao J S, Xu N X 2013 Acta Phys. Sin. 62 237302
Google Scholar
下载:
计量
- 文章访问数: 3825
- PDF下载量: 99
- 被引次数: 0
