基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构

计算物理 ›› 2016, Vol. 33 ›› Issue (3): 333-340.

基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构

叶珍宝, 周海京

北京应用物理与计算数学研究所, 北京 100094

收稿日期:2015-03-18 修回日期:2015-06-24 出版日期:2016-05-25 发布日期:2016-05-25
通讯作者: 周海京,E-mail:zhou_haijing@iapcm.ac.cn
作者简介:叶珍宝(1982-),女,博士,助理研究员,主要研究方向为电磁场时域数值计算和区域分解方法及其应用研究,E-mail:ye_zhenbao@iapcm.ac.cn
基金资助:
国家重点基础研究发展计划(2013CB328904),国家自然科学基金重点项目(61431014)资助

E-H Time-Domain Finite-Element Method with Hierarchical Vector Basis Functions for Analysis of Cavity and Waveguide Structures

YE Zhenbao, ZHOU Haijing

Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2015-03-18 Revised:2015-06-24 Online:2016-05-25 Published:2016-05-25

摘要/Abstract

摘要： 将高阶叠层矢量基函数用于E-H时域有限元方法,电场和磁场用相同的基函数展开并同时求解,时间离散采用Crank-Nicolson差分格式使得时间步长的选取摆脱稳定性条件的限制,同时采用完美匹配层来截断计算区域.对三维谐振腔及波导结构进行数值模拟与分析,结果表明,相较于低阶基函数,高阶叠层矢量基函数可以有效提高E-H时域有限元方法的计算精度.

关键词: 叠层矢量基函数, 高阶矢量基函数, E-H时域有限元方法, Crank-Nicolson差分格式

Abstract: We present E-H time-domain finite-element method using hierarchical high-order vector basis functions. Electric fields and magnetic fields are computed together and Crank-Nicolson difference scheme is implemented to obtain an unconditionally stable algorithm. Meanwhile, perfectly matched layers are used for truncation of unbounded regions. Three-dimensional cavity and waveguide structures are simulated. It shows that accuracy is improved by E-H TDFEM method with higher-order basis functions.

Key words: hierarchical basis functions, higher-order basis functions, time-domain finite-element method, Crank-Nicolson scheme

中图分类号:

叶珍宝, 周海京. 基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构[J]. 计算物理, 2016, 33(3): 333-340.

YE Zhenbao, ZHOU Haijing. E-H Time-Domain Finite-Element Method with Hierarchical Vector Basis Functions for Analysis of Cavity and Waveguide Structures[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(3): 333-340.

参考文献

[1] YEE K S. Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media[J]. IEEE Trans Antennas Propagat, 1966, 14:302-307.
[2] ZHOU H, LIU Y, LI H, et al. Computational electromagnetic and applications in numerical simulation of electromagnetic environment effects and development tendency[J].Chinese J Comput Phys, 2014, 31(4):379-389.
[3] KUANG X, WANG D, ZHANG L, et al. Finite difference time-domain method based on high order compact scheme[J]. Chinese J Comput Phys, 2014, 31(1):91-95.
[4] LEE J F. WETD-a finite element time-domain approach for solving Maxwell's equations[J]. IEEE Microwave and Guided Wave Letters, 1994, 4:11-13.
[5] GEDNEY S D, NAVSARIWALA U. An unconditionally stable finite element time-domain solution of the vector wave equation[J]. IEEE Microwave and Guided Wave Letters, 1995, 5:332-334.
[6] DONDERICI B, TEIXEIRA F L. Mixed finite-element time-domain method for transient Maxwell equations in doubly dispersive media[J]. IEEE Trans Microwave Theory Tech, 2008, 56(1):113-120.
[7] JIN J M. The finite element method in electromagnetics[M]. 3rd ed. New York:Wiley, 2014.
[8] WEBB J P, FORGHANI B. Hierarchal scalar and vector tetrahedral[J]. IEEE Trans Magn, 1993, 29:1495-1498.
[9] ANDERSEN L S, VOLAKIS J L. Hierarchical tangential vector finite elements for tetrahedral[J]. IEEE Microwave Guided Wave Lett, 1998, 8:127-129.
[10] SAVAGE J S, PETERSON A F. Higher-order vector finite elements for tetrahedral cells[J]. IEEE Trans Microwave Theory Tech, 1996, 44:874-879.
[11] WEBB J P. Hierarchal vector basis functions of arbitrary order for triangular and tetrahedral finite elements[J]. IEEE Transactions on Antennas and Propagation, 1999, 47(8):1244-1253.
[12] ZHU Y, CANGELLARIS A. Multigrid finite element methods for electromagnetic field modeling[M]. Hoboken, NJ:Wiley-IEEE, 2006.
[13] YANG Y, CHEN R S,WANG D X, YUNG E K N. Unconditionally stable Crank-Nicolson finite-difference time-domain method for simulation of 3-D microwave circuits[J]. IET Microw Antennas Propag, 2007, 1(4):937-942.
[14] GEDNEY S D. An anisotropic perfectly matched layer absorbing media for the truncation of FDTD lattices[J]. IEEE Transactions on Antennas and Propagation, 1996, 44(12):1630-1639.
[15] JIAO D, JIN J M, MICHIELSSEN E, et al. Time-domain finite-element simulation of three-dimensional scattering and radiation problems using perfectly matched layers[J]. IEEE Transactions on Antennas and Propagatiion, 2003,51(2):296-305.
[16] JN W P Carpes, PICHON L, RAZEK A. Efficient analysis of resonant cavities by finite element method in the time domain[J]. IEE Proc-Microw Antennas Propagation, 2000, 147(1):53-56.
[17] MOVAHHEDI M, ABDIPOUR A. Alternation-direction implicit formulation of the finite-element time-domain method[J]. IEEE Trans Microwave Theory Tech, 2007, 55:1322-1331.
[18] CHEN R S, YUNG E K N, CHAN C H, WANG D X, FANG D G. Application of the SSOR preconditioned CG algorithm to the vector FEM for 3-D full-wave analysis of electromagnetic-field boundary-value problems[J]. IEEE Trans Microwave Theory Tech, 2002, 50(4):1165-1172.