一种新的基于递归卷积方法的ADIFDTD+PML边界条件.pdfVIP

2005。全国微波毫米波会议论文集 一种新的基于递归卷积方法的 ADI．FDTD PML边界条件 王林年，梁昌洪，李龙 (西安电子科技大学电子工程学院，陕西西安710071) 摘要：本文提出了一种新的基于递归卷积方法的理想匹配层(PML)边界条件，适用于使 用交替方向隐式(ADI)时域有限差分法(FDTD)分析开区域电磁问题。这是一种基于复频 移(CFS)坐标伸缩变换的PML方法，其优点是可以吸收低频电磁波。数值结果表明， 和已有的类似方法相比，新方法可以将截断边界的反射误差减小30分贝左右，且明显 减少了计算时间。 关键词：交替方向隐式(ADI)；时域有限差分法(FDTD)；理想匹配层(PML)；复频移 (CFS) 1． 引言 交替方向隐式(ADI)时域有限差分法(FDTD)是一种无条件稳定的FDTD分析方法 【1]，其时间步的大小不再受至IJCFL稳定条件的限制，因而可以根据电磁信号的周期大小 或者根据对计算精度的要求来确定。 但是ADI．FDTD方法对计算机资源的要求增加了。当被分析对象具有过于细腻的几何 特征，因而空间网格不得不划分的非常细密，此时使用ADI．FDTD方法可以节约计算时 间。这就要求有一种适用ADI方法的、可以吸收低频电磁成分的边界条件。一种基于复 频移(CFS)坐标伸缩变换的PML边界条件被证明可以有效地吸收低频电磁成分[2，3]， 在文献[4]中，这一方案已被嫁接至UADI．FDTD方法中。 本文改进了文献[4]的方法，有关时间离散的过程作了适当的调整。数值结果表明， 改进后的方法可以大幅度地减小反射误差，提高计算速度。 2．ADI差分公式 为简单起见，本文以二维陋“隋况为例，且使用了归一化的电磁场量和Maxwell 方程。在PML区域以下方程代替了Maxwell方程组[3】 jcoc～E，=Sy-10。H： (1a) ycoc～E，=一Sx-lO。H： (1b) H jcoc一1 x—s Oc) z=sv-18vE x-18xEv 其中s。为坐标伸缩因子 (叩=x，Y) (2a) S口=％+盯口／(aq+joeo)， 因而 s_～=1／k”+％I(jco+dq) (2b) 其中b口=一仃口k々～／e。， [5】，其中增加了口。选项。这一方法被称为复频移．(CFS)坐标伸缩变换PML边界条件。 1257 当取％0时，PML边界条件可以吸收低频电磁波。 。 个莶盏簇皋1雾-荨、器I[骺1]o撒贼替方向隐舫法要批标准FDTD方法一 第一分步： 一 ，= 吖‘lb №协卜 日 №‰ m恸 (3a) f| ¨‘h¨ 光巴 √ E。洲 “ (3b) E髟蛾 ：Il 引销皿 + 圯r¨m忪 篇麓：W 。洲，川k 肌2 川忱刎暇70 墩。毗也 ∽吼儿