一种简化的时域非连续伽略金阻抗边界算法

杨谦; 魏兵; 李林茜; 邓浩川

doi:10.7498/aps.72.20222104

摘要

针对时域非连续伽略金(discontinuous Galerkin time-domain, DGTD)算法中的阻抗边界条件问题开展研究. 阻抗边界条件中的频域算符${\rm{j}}\omega $一般在根号内部, 其在时域数值算法中的实现有一定难度. 另一方面, DGTD算法中数值通量表达式也含有阻抗边界条件, 这也进一步增加了频时转换难度. 为了能给出简化的DGTD阻抗边界算法, 本文首先针对数值通量表达式进行推导, 得到一个特定函数$ {\tilde Z_R} $, 该函数包含频域算符${\rm{j}}\omega $, 函数以外表达式不含频域算符${\rm j}\omega$, 这样就可以仅处理$ {\tilde Z_{\rm{R}}} $的频时转换问题. 由于$ {\tilde Z_{\rm{R}}} $形式复杂, 对$ {\tilde Z_{\rm{R}}} $进行矢量匹配处理, 得到关于${\rm{j}}\omega$的一阶有理分式, 进而得到其时域迭代式. 这一过程简明、易于实施, 还可避开矩阵计算. 本文方案经一维及三维算例验证, 精度很好, 针对特定电磁问题如涂覆层问题可大幅降低计算时间.

关键词:

Abstract

Large-size conductive targets or coated targets are difficult problems in computational electromagnetics. In general, these problems can be classified as multi-scale problems. Multi-scale problems usually consume a large quantity of computational resources. A lot of efforts have been devoted to seeking for fast methods for these problems. When the skin depth is less than the size of a conductive target, the tangential component of the electric field and magnetic field over the surface of the target can be correlated by the surface impedance $ \tilde Z $. The $ \tilde Z $ is usually a complex function of the frequency, and it can be used to formulate an impedance boundary condition (IBC) to describe iterative equations in time domain methods, avoiding the volumetric discretization of the target and improving computational efficiency. This condition is commonly known as the surface impedance boundary condition (SIBC). Similarly, for a conductor whose thickness is in the order of skin depth or less, it also has high resource requirements, if the target is of direct volume discretization. The transmission impedance boundary condition (TIBC) can be utilized instead of a coated object to reduce resource requirements. Therefore, there is no need to discretize volume.There are few studies on the IBC scheme by using the discontinuous Galerkin time-domain (DGTD) method. Li et al. (Li P, Shi Y, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 5686; Li P, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 3065 ; Li P, Jiang L J, Bağcι H 2018 IEEE Trans. Antennas Propag. 66 3590 ) discussed the IBC scheme by using the DGTD, which involves complex matrix operations in the processing of IBC. In the DGTD method, numerical flux is used to transmit data between neighboring elements, and the key to the IBC scheme in DGTD is how to handle numerical flux. We propose a DGTD method with a simple form and matrix-free IBC scheme. The key to dealing with IBC in DGTD is numerical flux. Unlike the way in the literature, the impedance $ \tilde Z $ is not approximated by rational functions in our study. A specfic function $ {\tilde Z_R} $ obtained after the derivation in this work is approximated by rational functions in the Laplace domain through using the vector-fitting (VF) method, and its time-domain iteration scheme is given. This approach avoids matrix operations. The TIBC and SIBC processing schemes are also given. The advantage of the proposed method are that the upwind flux’s standard coefficients are retained and the complex frequency-time conversion problem is implemented by the vector-fitting method. The one-dimensional and three-dimensional examples also show the accuracy and effectiveness of our proposed method in this work.

Keywords:

impedance boundary condition /
vector-fitting method /
discontinuous Galerkin time domain method

作者及机构信息

1.
西安电子科技大学物理学院, 西安　710071

2.
北京环境特性研究所, 电磁散射重点实验室, 北京　100854

通信作者: 杨谦, qyang@xidian.edu.cn

基金项目: 国家自然科学基金(批准号: 61901324, 62001345)、中国博士后科学基金(批准号: 2019M653548, 2019M663928XB)、电波环境特性及模化技术国防科技重点实验室基金(批准号: 201903002)、电磁散射重点实验室基金(批准号: 61424090111)和中央高校基本科研业务费(批准号: XJS200501, XJS200507, JB200501)资助的课题.

Authors and contacts

1.
School of Physics, Xidian University, Xi’an 710071, China

2.
Science and Technology on Electromagnetic Scattering Laboratory, Beijing Institute of Environmental Features, Beijing 100854, China

Corresponding author: Yang Qian, qyang@xidian.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 61901324, 62001345), the China Postdoctoral Science Foundation (Grant Nos. 2019M653548, 2019M663928 XB), the Foundation of National Key Laboratory of Electromagnetic Environment, China (Grant No. 201903002), the Foundation of the Science and Technology on Electromagnetic Scattering Laboratory, China (Grant No. 61424090111), and the Fundamental Research Funds for the Central Universities of China (Grant Nos. XJS200501, XJS200507, JB200501).

文章全文 : translate this paragraph

参考文献

[1]	Senior T 1960 Appl. Sci. Res. 8 418 Google Scholar
[2]	Feliziani M 2011 IEEE Trans. Electromagn. Compat. 54 299 Google Scholar
[3]	Nayyeri V, Soleimani M, Ramahi O M 2013 IEEE Trans. Electromagn. Compat. 56 385 Google Scholar
[4]	Hanson G W 2008 J. Appl. Phys. 103 064302 Google Scholar
[5]	Nayyeri V, Soleimani M, Ramahi O M 2013 IEEE Trans. Antennas Propag. 61 6107 Google Scholar
[6]	Shapoval O V, Gomez-Diaz J S, Perruisseau-Carrier J, Mosig J R, Nosich A I 2013 IEEE Trans. Terahertz Sci. Technol. 3 666 Google Scholar
[7]	da Costa K Q, Dmitriev V, Nascimento C M, Silvano G L 2014 Microwave Opt. Technol. Lett. 56 1019 Google Scholar
[8]	Beggs J H, Luebbers R J, Yee K S, Kunz K S 1992 IEEE Trans. Antennas Propag. 40 49 Google Scholar
[9]	oh K S, Schutt-Aine J E 1995 IEEE Trans. Antennas Propag. 43 660 Google Scholar
[10]	Kobidze G 2010 IEEE Trans. Antennas Propag. 58 2394 Google Scholar
[11]	Gustavsen B, Semlyen A 1999 IEEE Trans. Power Delivery 14 1052 Google Scholar
[12]	Yi M, Ha M, Qian Z, Aydiner A, Swaminathan M 2013 IEEE Trans. Microwave Theory Tech. 61 4029 Google Scholar
[13]	Glisson A W 1992 Radio Sci. 27 935 Google Scholar
[14]	Yan S, Jin J M 2013 IEEE Trans. Antennas Propag. 61 5533 Google Scholar
[15]	Ylä-Oijala P, Kiminki S P, Järvenpää S 2010 IEEE Trans. Antennas Propag. 58 3997 Google Scholar
[16]	Gyselinck J, Dular P, Geuzaine C, Sabariego R 2009 IEEE Trans. Magn. 45 1280 Google Scholar
[17]	Li P, Shi Y, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 5686 Google Scholar
[18]	Li P, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 3065 Google Scholar
[19]	Li P, Jiang L J, Bağcı H 2018 IEEE Trans. Antennas Propag. 66 3590 Google Scholar
[20]	Wang P, Shi Y, Tian C Y, Li L 2018 IEEE Antennas Wirel. Propag. Lett. 17 2169 Google Scholar

施引文献

图 1 半空间示意图

Fig. 1. Plane wave incident on a lossy dielectric half-space.

下载: 全尺寸图片幻灯片

图 2 虚拟单元

Fig. 2. Virtual element.

下载: 全尺寸图片幻灯片

图 3 SIBC($ {\tilde Z_R} $)

Fig. 3. SIBC ($ {\tilde Z_R} $).

下载: 全尺寸图片幻灯片

图 4 薄涂层示意图

Fig. 4. A conductive object with thin sheet.

下载: 全尺寸图片幻灯片

图 5 TIBC ($ {\tilde Z_R} $)

Fig. 5. TIBC ($ {\tilde Z_R} $).

下载: 全尺寸图片幻灯片

图 6 对于一维SIBC算例, 反射波及入射波时域数据

Fig. 6. Incident wave and reflected wave for one-dimensional SIBC example.

下载: 全尺寸图片幻灯片

图 7 对于一维SIBC算例, 反射系数幅值

Fig. 7. Magnitude of the reflection coefficient for one-dimensional SIBC example.

下载: 全尺寸图片幻灯片

图 8 对于一维SIBC算例, 反射系数相位

Fig. 8. Phase of the reflection coefficient for one-dimensional SIBC example.

下载: 全尺寸图片幻灯片

图 9 对于一维TIBC算例, 反射波及入射波时域数据

Fig. 9. Incident wave and reflected wave for one-dimensional TIBC example.

下载: 全尺寸图片幻灯片

图 10 对于一维TIBC算例, 反射系数幅值

Fig. 10. Magnitude of the reflection coefficient for one-dimensional TIBC example.

下载: 全尺寸图片幻灯片

图 11 对于一维TIBC算例, 反射系数相位

Fig. 11. Phase of the reflection coefficient for one-dimensional TIBC example.

下载: 全尺寸图片幻灯片

图 12 单站RCS结果

Fig. 12. Monostatic RCS.

下载: 全尺寸图片幻灯片

图 13 涂覆球

Fig. 13. Coated sphere.

下载: 全尺寸图片幻灯片

图 14 涂覆球单站RCS结果

Fig. 14. Monostatic RCS.

下载: 全尺寸图片幻灯片

表 1 数值通量系数表达式

Table 1. Coefficients of numerical flux.

	通量系数
	$ k_E^e $	$ k_H^e $	$ v_H^e $	$ v_E^e $
迎风通量表达式(upwind flux)	$\dfrac{{{{\tilde Y}^{e + }}}}{{{Y^e} + {{\tilde Y}^{e + }}}}$	$\dfrac{{{{\tilde Z}^{e + }}}}{{{Z^e} + {{\tilde Z}^{e + }}}}$	$\dfrac{1}{{{Y^e} + {{\tilde Y}^{e + }}}}$	$\dfrac{1}{{{Z^e} + {{\tilde Z}^{e + }}}}$

下载: 导出CSV

表 2 计算参数

Table 2. Calculation parameters.

	离散尺度/m	四面体数量	时间步/(10^–11 s)	计算时间/min	均方根误差/dB
DGTD (SIBC)	0.05	402094	0.328972	54	0.6237
DGTD (standard)	0.05	487677	0.328447	61	0.5836

下载: 导出CSV

表 3 计算参数

Table 3. Calculation parameters.

离散尺度/m	四面体数量	时间步/(10^–11 s)	计算时间/min	均方根误差/dB
0.05	402094	0.328972	55	0.8782

下载: 导出CSV

[1]	Senior T 1960 Appl. Sci. Res. 8 418 Google Scholar
[2]	Feliziani M 2011 IEEE Trans. Electromagn. Compat. 54 299 Google Scholar
[3]	Nayyeri V, Soleimani M, Ramahi O M 2013 IEEE Trans. Electromagn. Compat. 56 385 Google Scholar
[4]	Hanson G W 2008 J. Appl. Phys. 103 064302 Google Scholar
[5]	Nayyeri V, Soleimani M, Ramahi O M 2013 IEEE Trans. Antennas Propag. 61 6107 Google Scholar
[6]	Shapoval O V, Gomez-Diaz J S, Perruisseau-Carrier J, Mosig J R, Nosich A I 2013 IEEE Trans. Terahertz Sci. Technol. 3 666 Google Scholar
[7]	da Costa K Q, Dmitriev V, Nascimento C M, Silvano G L 2014 Microwave Opt. Technol. Lett. 56 1019 Google Scholar
[8]	Beggs J H, Luebbers R J, Yee K S, Kunz K S 1992 IEEE Trans. Antennas Propag. 40 49 Google Scholar
[9]	oh K S, Schutt-Aine J E 1995 IEEE Trans. Antennas Propag. 43 660 Google Scholar
[10]	Kobidze G 2010 IEEE Trans. Antennas Propag. 58 2394 Google Scholar
[11]	Gustavsen B, Semlyen A 1999 IEEE Trans. Power Delivery 14 1052 Google Scholar
[12]	Yi M, Ha M, Qian Z, Aydiner A, Swaminathan M 2013 IEEE Trans. Microwave Theory Tech. 61 4029 Google Scholar
[13]	Glisson A W 1992 Radio Sci. 27 935 Google Scholar
[14]	Yan S, Jin J M 2013 IEEE Trans. Antennas Propag. 61 5533 Google Scholar
[15]	Ylä-Oijala P, Kiminki S P, Järvenpää S 2010 IEEE Trans. Antennas Propag. 58 3997 Google Scholar
[16]	Gyselinck J, Dular P, Geuzaine C, Sabariego R 2009 IEEE Trans. Magn. 45 1280 Google Scholar
[17]	Li P, Shi Y, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 5686 Google Scholar
[18]	Li P, Jiang L J, Bağcι H 2015 IEEE Trans. Antennas Propag. 63 3065 Google Scholar
[19]	Li P, Jiang L J, Bağcı H 2018 IEEE Trans. Antennas Propag. 66 3590 Google Scholar
[20]	Wang P, Shi Y, Tian C Y, Li L 2018 IEEE Antennas Wirel. Propag. Lett. 17 2169 Google Scholar

[1]	杨谦, 魏兵, 李林茜, 邓浩川. 一种强鲁棒性的时域非连续伽略金可穿透薄层算法. 物理学报, 2023, 72(8): 080202. doi: 10.7498/aps.72.20222230
[2]	马奥杰, 陈颂佳, 李玉秀, 陈颖. 纳米颗粒布朗扩散边界条件的分子动力学模拟. 物理学报, 2021, 70(14): 148201. doi: 10.7498/aps.70.20202240
[3]	黄燕燕, 张旭琳, 杨伟, 王笑冰, 雷蕾, 彭文达, 徐平. 微光学元件宽带光源照明下无色散边界条件探讨. 物理学报, 2019, 68(22): 224203. doi: 10.7498/aps.68.20190716
[4]	林晨森, 陈硕, 肖兰兰. 适用复杂几何壁面的耗散粒子动力学边界条件. 物理学报, 2019, 68(14): 140204. doi: 10.7498/aps.68.20190533
[5]	岳松, 张兆传, 高冬平. 阻抗匹配条件下磁控管的注入锁频. 物理学报, 2013, 62(17): 178401. doi: 10.7498/aps.62.178401
[6]	杨利霞, 马辉, 施卫东, 施丽娟, 于萍萍. 基于表面阻抗边界条件的等离子体薄涂层电磁散射的时域有限差分分析. 物理学报, 2013, 62(3): 034102. doi: 10.7498/aps.62.034102
[7]	颛孙旭, 马西奎. 一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法. 物理学报, 2012, 61(11): 110206. doi: 10.7498/aps.61.110206
[8]	王龙, 钟易成, 张堃元. 金属/介质涂覆的S形扩压器电磁散射特性. 物理学报, 2012, 61(23): 234101. doi: 10.7498/aps.61.234101
[9]	张旭, 周玉泽, 闭强, 杨兴华, 俎云霄. 有边界条件的忆阻元件模型及其性质. 物理学报, 2010, 59(9): 6673-6680. doi: 10.7498/aps.59.6673
[10]	吴兆春. 导热几何形状反演的变分原理及边界条件的确立. 物理学报, 2010, 59(9): 6326-6330. doi: 10.7498/aps.59.6326
[11]	张燕萍, 赵晓鹏, 保石, 罗春荣. 基于阻抗匹配条件的树枝状超材料吸收器. 物理学报, 2010, 59(9): 6078-6083. doi: 10.7498/aps.59.6078
[12]	谭康伯, 梁昌洪. 高频条件下运动边界的电磁场边界条件. 物理学报, 2009, 58(10): 6770-6771. doi: 10.7498/aps.58.6770
[13]	李博, 王延申. 可积开边界条件下的q形变玻色子模型. 物理学报, 2007, 56(3): 1260-1265. doi: 10.7498/aps.56.1260
[14]	唐黎明, 王艳, 王丹, 王玲玲. 边界条件对介电量子波导中声子输运性质的影响. 物理学报, 2007, 56(1): 437-442. doi: 10.7498/aps.56.437
[15]	张波. 二维介质柱光子晶体波导吸收边界条件. 物理学报, 2005, 54(12): 5677-5682. doi: 10.7498/aps.54.5677
[16]	梁昌洪, 褚庆昕. 运动边界的电磁场边界条件. 物理学报, 2002, 51(10): 2202-2204. doi: 10.7498/aps.51.2202
[17]	韩一平, 吴振森. 椭球粒子电磁散射的边界条件的讨论. 物理学报, 2000, 49(1): 57-60. doi: 10.7498/aps.49.57
[18]	颜家壬. 边界条件对非传播孤子运动的影响. 物理学报, 1995, 44(10): 1571-1576. doi: 10.7498/aps.44.1571
[19]	刘登云. 含时边界条件和Berry相位. 物理学报, 1993, 42(5): 705-710. doi: 10.7498/aps.42.705
[20]	吴式玉, 周子舫. 边界条件对无序体系本征态的影响. 物理学报, 1984, 33(12): 1650-1660. doi: 10.7498/aps.33.1650

计量

文章访问数: 3348
PDF下载量: 72
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板