一种基于波动方程的新高阶FDTD方法.pdfVIP

第 19卷 第l2期 强 激 光 与 粒 子 束 Vo1．19，NO．12 2007年 l2月 HIGH POWER LASER AND PARTICLE BEAMS Dec．。2007 文章编号： i001-4322(2007)12—2069—05 一 种基于波动方程的新高阶FDTD方法 王 烁 ， 邵振海 ， 文光俊 ， 李世贵。， 谢 康 (1．电子科技大学通信与信息工程学院射频集成电路实验室，成都610054； 2．Micro Design Technology Pte Ltd，Singapore 117674 3．重庆科技学院数理学院，重庆400042； 4．电子科技大学 宽带光纤传输与通信网技术教育部重点实验室，成都 610054) 摘 要： 提出了一种基于二阶波动方程的(2M，4)高阶时域有限差分(FDTD)方法，通过使用辛积分传播 子(SIP)在时域上获得 4阶精度，使用离散奇异卷积(DSC)方法在空域上达到 2M阶精度。与已有的(2M，4) 阶FDTD方法相比，虽然两者都采用SIP和DSC方法，但是此二者的不同点在于：第一，新方法基于二阶波动 方程；第二，在离散计算空间时使用单一网格而不是传统的Yee网格；第三，单独计算某一场分量从而节约内存 并减少计算量。数值计算结果表明，与传统高阶算法相比，基于波动方程的高阶FDTD方法耗费的机时只有 它的5O ，内存消耗下降 1O ，而两者的计算结果之间相对误差小于5‰。 关键词： 时域有限差分； 波动方程； 离散奇异卷积； 拉格朗日多项式； 辛积分传播子 中图分类号： TN0l1； TP3l1．11 文献标识码： A 传统时域有限差分(FDTD)法[1 在解决实际电磁场问题时，需要将6个场分量同时进行交叉迭代运算。 即使在2维TM／TE问题中，基于Yee网格的传统FDTD法至少需要3个分量。而从二阶波动方程出发的 2D-FDTDc。 只需要计算3个场分量中的一个，在导波结构场分析等应用中具有较好的灵活性。另一方面，尽 管FDTD方法有着广泛的应用，但其精度却无法与以格林函数为基础的积分方法相比。近年来，学者们提出 了许多方案来提高FDTD的精度，例如具有(4，4)阶精度的基于龙格一库塔方法的FDTD[3]、基于辛积分传播子 (SIP)方法的FDTD[4 ]、基于拉格朗日多项式插值与SIP方法[E 以及基于离散奇异卷积(DSC)与龙格一库塔方 法的高阶FDTD【7 等。然而，从已发表的文献来看，以上这些方案都是Yee网格的FDTD法。 本文提出一种基于波动方程的具有(2M，4)阶精度的FDTD算法。其特点如下：首先，波动方程是一个场 分量的方程，故该方案在计算时只需要考虑单个场分量，减少了计算场量，从而节约了计算量和存储量。又因 为在离散计算区域时只使用单一网格而不需要引入形如Yee空问的交错网格，因此空间结构得以简化。其 次，该(2M，4)方案在空域和时域上都采用高阶算法，其精度远远高于标准的FDTD。在空域上，采用基于拉格 朗日 核的DSC方法，时域上采用4阶精度的SIP方法。选择适当的计算带宽M值能够大大提高计算精度。 本文介绍基于波动方程的高阶FDTD方法的基本理论，以导波系统为例对所提出的方法进行了计算精度 验证并讨论了该方法的优势。 1 基本原理 1．1 空域2M阶DSC方法 文献[7]提出了一种DSC时域方法，其在计算电磁场问题时，引用DSC方法离散计算空问，构造高阶插值 方程从而提高计算精度。 假设 T表示系统的特征函数，厂(z)为一个测试函数，则．厂(z)的奇异卷积表示为 广+∞ F(￡)=T(￡)*厂(￡)一I 了、(￡一z)厂(z)dx (1) J 一 ∞ 其离散的奇异卷积(DSC)形式为 F。(￡)一 T。(￡一zt)厂(zt)