李林茜 魏兵 杨谦 葛德彪 王飞

(1.西安电子科技大学物理与光电工程学院,西安 710071;2.西安电子科技大学信息感知技术协同创新中心,西安 710071)



二维TM波时域非连续伽略金算法理论数值通量研究

李林茜1,2魏兵1,2杨谦1,2葛德彪1,2王飞1,2

(1.西安电子科技大学物理与光电工程学院,西安 710071;2.西安电子科技大学信息感知技术协同创新中心,西安 710071)

采用数值通量的方式进行场量交互是时域非连续伽略金(Discontinuous Galerkin Time Domain,DGTD)算法区别于时域有限元(Finite Element Time Domain, FETD)方法的主要方面.从二维TM情形弱解方程出发,讨论了当前三角形单元和相邻单元进行场量交互时数值通量物理意义和不同形式.结合数值通量和弱解方程得到了DGTD算法的迭代计算式.给出了线元辐射和双线元干涉的数值算例,算例结果表明了文中方法的正确性.

时域非连续伽略金算法;算法值通量;结点基函数

DOI 10.13443/j.cjors.2015111702

引 言

时域电磁仿真方法因其能够通过一次时域计算结合傅里叶变换即可得到宽频带信息的特点,近年来受到人们广泛关注.时域有限差分(Finite Difference Time Domain, FDTD)方法[1-2]是目前电磁学领域内被人们广泛、深入地研究,并取得巨大成功应用的方法.该方法原理直观、编程简便、实用性强,一直在时域方法中占主导地位.但由于FDTD方法采用规则网格剖分,其建模的能力弱,对弯曲表面的阶梯近似,严重限制了将其应用于复杂几何结构时的计算精度.另外一种时域算法为时域有限元(Finite Element Time Domain, FETD)方法,该算法采用非结构单元拟合(如四面体),与真实目标的差异……