TDEFIE的奇异性积分的精确快速计算.pdfVIP

第19卷增刊 电波科学学报 V01．19，Sup． RADIO8CIflENCE October，2004 2004年10月aⅡ憾JOURN．Mt．,OF TDEFIE的奇异性积分的精确快速计算+ 赵延文武胜波聂在平徐建华 (电子科技大学电子工程学院，ywzhao@u鹤tc．edu．cn，I匹llJl成都610054) 摘要利用时域积分方程(TDIE)法模拟时域散射现象越来越受到学者们的关注。 用于求解TDIE的时间步进(MoT)算法的后时不稳定性严重阻碍了TDIE法的广 泛应用和发展。本文首先利用参数坐标和Duffy坐标变换将TDEFIE的奇异性积 分转换成非奇异性积分，然后根据时间基函数的特点将该积分转换成可以快速精确 计算的分区域积分。数值结果表明，该方法大幅提高了利用MOT算法求解TDE— FIE的MOT算法的后时稳定性和计算精度。 关键词 时域积分方程，时间步进算法，奇异性积分 生感应电流J(，，￡)。利用理想导体表面上切向电 1 引 言 场为零的边界条件，可得感应电流J(r，￡)所满足的 近年来，利用时域方法模拟电磁现象越来越受 时域电场积分方程为 到人们的重视。在分析均匀无限大介质中的理想导 体或均匀介质体的表面时域散射问题时，采用时域 积分方程(TDIE)方法比采用时域微分方程 (TDDE)方法更具优势：TDIE法的未知量仅仅位于 n …～I。、-… 散射目标的表面；通过格林函数自动满足辐射条件， 其中，下标“tan”表示切向分量，推迟时间r=t—R／ 无需强加吸收边界条件；网格色散误差小。然而，基 r—r，l，￡和卢 C，为场点r和源点r7间的距离R—l 于时间步进(MOT)算法求解时域电场积分方程 分别为散射体周围介质的介电常数和磁导率，c为 (TDEFIE)所出现的时间后期不稳定性严重地阻碍 电磁波在其中的传播速度。 了TDIE法的广泛应用和发展。有许多学者对 利用矩量法求解方程(1)，可以求得电流J(r， TDIE的时间后期稳定性问题作了广泛深入的研 ￡)。为此，将感应电流J(r，￡)用空间基函数^(，)(咒 究：例如，利用混合场积分方程(CFIE)法来克服内 谐振所导致的不稳定性[13；采用时间平均等滤波方 法来克服高频成份的影响[2]以及通过避免直流成 展开为 MNr 份在时间迭代过程的积累来提高时间后期稳定 J(r，￡)一∑∑P．Tf(t)，。(r)(2) 性[33等等。实际上，MOT算法出现时间后期不稳 定性主要是由于在离散时域积分方程时采用了不精 其中E为与时空基函数f。(r)正(t)相关的待求加 确的数值计算方法和不恰当的解析近似所导致的， 权系数。在本文中，选取时间基函数Tl(t)=T(￡一 精确计算TDEFIE的时域阻抗矩阵元素(特别是自 lAt)(At为时间步长)为三角形函数、空间基函数^ 阻抗元素)可以大幅度提高利用MOT算法求解 TDEFIE的计算精度和时间后期稳