坐标伸缩完全匹配层在二维柱坐标系下的实现与应用

系统仿真学报 ›› 2015, Vol. 27 ›› Issue (7): 1644-1651.

坐标伸缩完全匹配层在二维柱坐标系下的实现与应用

朱安石¹, 陈自力¹, 杨勇², 甄云卉³

1.军械工程学院,河北石家庄 050003;
2.轻武器研究所,北京 102202;
3.河北省军区,河北石家庄 050011

收稿日期:2014-04-29 修回日期:2014-07-26 出版日期:2015-07-08 发布日期:2020-07-31
作者简介:朱安石(1985-),男,河北承德,博士,研究方向为计算电磁学;陈自力(1964-),男,山西太原,教授,研究方向为无人机数据链抗干扰;杨勇(1983-),男,湖北宜昌,硕士,研究方向为轻武器测试技术。

Implementation and Application of Coordinate Stretched PML under Two-dimensional Cylindrical Coordinate

Zhu Anshi¹, Chen Zili¹, Yang Yong², Zhen Yunhui³

1. Mechanical Engineering College, Shijiazhuang 050003, China;
2. Small Arms Institute, Beijing 102202, China;
3. Hebei Provincial Military Command, Shijiazhuang 050011, China

Received:2014-04-29 Revised:2014-07-26 Online:2015-07-08 Published:2020-07-31

摘要/Abstract

摘要： 针对柱坐标系下电磁场的数值计算,常规的分裂场吸收边界已不再适用,将坐标伸缩完全匹配层应用到二维扩展柱坐标系中,应用辅助方程法推导了Maxwell方程在该匹配层中的时域迭代方程,仿真结果表明: 理论推导正确,且吸收边界的反射系数小于-30dB,能够满足数值计算精度需求。将该匹配层应用于单极同轴线馈电天线辐射特性的研究,得到了不同激励源下的近场和远场方向图。该研究为柱坐标系下电磁场散射与辐射的数值计算提供了一种新的思路和参考。

关键词: 时域有限差分法, 完全匹配层, 数值计算, 方向图

Abstract: According to the numerical calculation of electromagnetic field under cylindrical coordinate, the conventional split field absorbing boundary is no longer applicable. The coordinate stretched perfectly matched layers (CS-PML) were applied in the two-dimensional stretched cylindrical coordinate system. The updating equations of Maxwell curl equations in CS-PML were derived by using auxiliary equation method. The simulation results show that theoretical derivation is correct. The reflection coefficient of CS-PML is less than -30 dB, and it is able to meet the precision demand of numerical calculation. Besides, the CS-PML was adopted to investigate the radiation characteristic of the monopole antenna fed by coaxial cable. And the near-field and far-field pattern were obtained. The investigation provides a new idea and reference for the numerical calculation of electromagnetic scattering and radiation problems.

Key words: finite-difference time-domain (FDTD), perfectly matched layer (PML), numerical calculation, radiation pattern

中图分类号:

TN92

朱安石, 陈自力, 杨勇, 甄云卉. 坐标伸缩完全匹配层在二维柱坐标系下的实现与应用[J]. 系统仿真学报, 2015, 27(7): 1644-1651.

Zhu Anshi, Chen Zili, Yang Yong, Zhen Yunhui. Implementation and Application of Coordinate Stretched PML under Two-dimensional Cylindrical Coordinate[J]. Journal of System Simulation, 2015, 27(7): 1644-1651.

参考文献

[1] B Engquist, A Majda.Absorbing boundary conditions for the numerical simulation of waves[J]. Mathematics of Computation (S0025-5718), 1977, 31(139): 629-651.
[2] G Mur.Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations[J]. IEEE Trans. Electomagn. Compat.(S0018-9375), 1981, 23(4): 377-382.
[3] Liao Z P, Yang B P, Yuan Y F, et al. A transmitting boundary for transient wave analysis[J]. Scientia Sinica(series A)(S0253-5831), 1984, 27(10): 1063-1076.
[4] 王玥, 王建国, 张殿辉. 全局稳定性一致的共形卷积完全匹配层[J]. 强激光与粒子束, 2013, 25(2): 441-445.
[5] 杨利霞, 梁庆, 于萍萍. 三维新型非分裂场完全匹配层吸收边界条件[J]. 电波科学学报, 2011, 26(1): 67-71.
[6] 柴焱杰, 孙继银, 李琳琳, 等. UPML吸收边界条件统一建模与优化方法[J]. 计算物理, 2011, 28(3): 413-419.
[7] Chew W C, W H Weedon. A 3D perfectly matched medium from modified Maxwell’s equations with stretched coordinates[J]. Microwave and Optical Technology Letters (S0895-2477), 1994, 7(13): 599-604.
[8] 朱翼超, 罗根新, 方向, 等. CPML在脉冲辐射场数值模拟中的应用研究[J]. 系统仿真学报, 2012, 24(6): 1189-1191.
[9] 刘亚文, 陈亦望, 徐鑫. 基于辅助差分方程的完全匹配层在时域多分辨率分析算法中的应用与性能分析[J]. 物理学报, 2013, 62(3): 034101.
[10] Feng N X, Li J X.A Z-transform Implementation of Higher-order Perfectly Matched Layer[J]. Chinese Journal of Computational Physics (S1001246X), 2012, 29(2): 271-275.
[11] Zhang X L, Gao H T, Chen C.FDTD Analysis and Appli-cation of Sleeve Monopole Antenna[J]. Wuhan University Journal of Nature Sciences (S1007-1202), 2006, 12(2): 284-288.
[12] Elsherbeni A, Demir V.The Finite Difference Time Domain Method for Electromagnetics with MATLAB Simulations[M]. San Francisco, USA: Peter Classman, 2009.
[13] Maloney J G, Smith G S, Scott W R.Accurate computation of the radiation from simple antennas using the finite-difference time-domain method[J]. IEEE Transactions on Antennas and Propagation (S0018-926X), 1990, 38(7): 1059-1068.