添加各向异性参数修正ADIFDTD算法中的数值色散研究.pdfVIP

添加各向异性参数修正ADI—FDTD 算法中的数值色散 郑宏兴郁道银 天津大学精密仪器与光电子工程学院，天津300072 Email：hxzheng@tju．edu．cn 摘要：对变换方向隐式(ADD时域有限差分(FDTD)算法提出了一种改进方法，即通过添加各向异性 参数来修正相速度误差，从而减少数值色散，提高计算精度。给出了改进后的ADI--FDTD计算 公式，并对算法的稳定性、色散关系以及各向异性参数的选取进行了研究，数值结果证明了所给 方法的有效性。 关键词：时域有限差分各向异性相速度数值色散 过程～： 一、fIu吾 为了减小时域有限差分法(FDTD)的数值 FP72 。x 色散，Juntunen等采用在计算空间添加各向异 li+1／2，J，I 色‰_。 性参数的手段【l】。实现这种方法较容易，也不 E =C：E、，f” bn+l／：2亿。 ’li,j+l／2．女 增加计算机资源——cPU时间和内存的占用。 F[n+1／2 “z lf，j，k+112 跳m州： T．Namiki提出的变换方向隐式(ADD的FDTD 方法拉1能够克服Coumt Friedrich。Levy(CVL)条 件【3】对时间步长的限制。但是随着时间步长的 增大，ADI—FDTD方法的数值色散效应也会增 捆：『 a日，r2 8，ay 大。因此数值色散成为影响时间步长选取的主 lⅢ，2．m exOz[f+1，2．m 要约束。本文将Juntunen的方法运用到 +砬塑I ～亟I ADI．FDTD方法中，给出了一种改进方法。即 sy吼州，2tg，孤|f√+1／2．女 通过添加各向异性参数来修正相速度误差，从 而减少数值色散，提高ADI．FDTD方法的计算 汨，I“ 弛r心 精度和效率。 ￡，叭_州，。吼_㈧，： 二、改进的ADI．FDTD算法 日P72 ⅣJ” 一。li,j+l／2．k+1／2 一xli,j+l／2．k+1／2 在三维各向同性介质中，介电常数岛和胁 Hn+l／2 H～r Y =GY h+l／2。』，k+l／2 1i+1／2，，，k+l／2 分型乘以各向异性参数8=diag[Ex，￡y，g：】阿一72 Ⅳl“ 一z 一z[i+1／2．j+l／2，k I／+1／2。j+112．k 和p=diag[儿，p，，p：】，这些添加的参数 与电磁场各分量对应。不失一般性，我们假设