摘要
为了克服在传统时域有限差分算法中时间步长受网格尺寸大小的限制,提出一种用于精确模拟波在色散金属结构传播的3维交替方向隐式时域有限差分算法.首先,利用广义关键点色散模型描述金属材料在光学频率下的电磁特性;然后,采用2个1阶微分方程代替2阶广义辅助微分方程的方法,避免时间步长引起的不稳定,降低算法的复杂度和应用难度;最后,用交替方向隐式时域有限差分算法研究1维等离子体模型和3维金属谐振腔结构,将交替方向隐式时域有限差分算法的结果与传统时域有限差分算法进行对比.结果表明二者吻合较好,但交替方向隐式时域有限差分算法的效率更高.
In order to overcome the time step size that was limited by the size of the grid in traditional finite difference time domain algorithm,a unified three-dimensional alternating direction implicit finite difference time domain algorithm was presented for accurately simulating the propagation of waves in a dispersion metal structure.First,the generalized critical points dispersion model was used to describe the electromagnetic properties of metallic materials at optical frequencies.Then,two first-order differential equations were used instead of the second-order generalized auxiliary differential equations,which avoided the instability caused by the time step component,and reduced the complexity and application difficulty of the algorithm.Finally,the alternating direction implicit finite difference time domain algorithm was used to study the one-dimensional plasma model and the three-dimensional metal resonator structure.The results of the alternate direction implicit finite difference time domain algorithm were compared with the traditional finite difference time domain algorithm.The results showed that the two agreed well,and the implicit finite difference time domain algorithm with alternating directions had higher efficiency.
作者
高林
谢国大
黄志祥
许杰
吴先良
GAO Lin;XIE Guoda;HUANG Zhixiang;XU Jie;WU Xianliang(Key Laboratory of Computational Intelligence and Signal Processing, Ministry of Education, Anhui University, Hefei 230039, China)
出处
《安徽大学学报(自然科学版)》
CAS
北大核心
2020年第4期52-59,共8页
Journal of Anhui University(Natural Science Edition)
基金
国家优秀青年科学基金资助项目(61722101)
国家自然科学基金资助项目(61871001,61471001)。
关键词
时域有限差分法
交替方向隐式时域有限差分法
金属结构
关键点模型
辅助微分方程
finite difference time domain
alternating direction implicit finite difference time domain
metal structure
critical points model
auxiliary differential equation