第七讲吸收边界条件——局域abc.docVIP

第5讲 吸收边界条件 当利用FDTD模拟开区域的电磁场问题时，由于计算机内存有限，必须对计算区域截断，使计算区域有限。 在截断边界上必须满足吸收边界条件(ABC)，以模拟电磁波无反射地通过截断边界，向无限远处传播。理想的ABC是难以实现的，通常只能采用近似的ABC。 对于近似ABC，要求 能够模拟向外传播的波； 引入的反射应足够小，对计算结果的影响可忽略； 保证算法稳定。 FDTD差分方程是不能用于截断边界的，因为它需要边界外面的点。根据实现ABC的途径，ABC分为如下三大类： 全域ABC 根据Green定理可知，时刻边界上的场可以用期间边界内部的所有点的场表示。这种ABC是严格的，但它涉及空间所有点在整个时间进程中的场值，即是“全域”的，所以，计算复杂，占内存大。实际模拟中很少采用。 局域ABC 局域ABC是利用边界附近点相邻几个时刻内的场来构造ABC，这种方法计算方便，占内存少。局域ABC的最大难点是如何建立截断边界上所满足的波以任何角度无反射地向外传播的方程。由于方程不同，而出现了各种类型的局部ABC，局部 ABC是当前FDTD模拟的主流。 损耗媒质ABC 在截断边界外部加上一层虚拟的损耗媒质，调整媒质参数和层厚度，就可以使波无反射地通过ABC，理想匹配层（PML）是这种ABC的典型代表，也是目前最好的ABC。 5.1 局域ABC 5.1.1 单向波动方程 单向波动方程: 波只沿一个方向传播的偏微分方程。如： (5-1) 表示波沿与轴夹角为的方向传播的波动方程。 上式并不能直接用于截断边界，因为在截断边界处通常并不知道波的传播方向，而且在边界上不同处入射角也不相同。 Engquist和Majda给出了适合于直角FD-TD网格吸收边界条件的单向波动方程理论。他们的理论可以用偏导数算子的因式分解来解释。 以直角坐标系中的二维波动方程为例 (5-2) 其中，偏微分算子L为 (5-3) (5-2)的解为 式中、、、表示沿四个不同方向传输的波，以及 利用算子的因式分解，可以得到 （5-4） 式中, (5-5a) 考虑图5-1所示的以任意角度从的区域入射到的边界上的波，则 因此，方程可以完全地吸收以任意角度从的区域入射到的边界上的波。同理可以证明，可以完全地吸收以任意角度从区域入射到的边界上的波。所以，和就是单向波动方程。 遗憾的是中方根的存在使其成为一类伪微分算子，妨碍了作为吸收边界条件的直接数值实现。不过，我们可以利用近似方根函数，产生一系列正常的偏微分方程，并进行数值实现，从而构成ABC。 这些近似使方程变得不精确，因而，当数值波通过截断边界时将产生一定的反射，但是这种反射可以减少到足够小。 近似方根函数的方法很多，常用是Taylor级数近似。设外向波沿方向变换缓慢，将（5-5）改写为 (5-5b) 式中，。 一阶Taylor近似 (5-6) 适用于非常小的情况。非常小意味非常小，即外向波沿方向变化很小，等效地，波以基本上垂直于外边界的角度入射到外边界上。 (5-6)代入(5-5b)，得 (5-7) 即一阶精度的单向波动方程 (5-8) 这正是沿方向传播的波方程。 为了检验（5-8）的吸收特性，设 代入（5-8），得 于是，反射系数为 式中表示入射角度。可见时反射为零，越大反射越大，时，全反射。 二阶Taylor 近似 (5-9) (5-10) 将两边同乘以，则有 即 （5-11） 将(5-8)和(5-11)进行数值离散，就构成了Mur一阶和二阶ABC。但是，并非展开项越多越精确，当三阶展开后，离散化会出现数值不稳定。 采用与前面类似的方法，可以得到反射系数为 同样，当时反射为零，越大反射越大，时，全反射。 但二阶近似比一阶近似反射特性有所改善。 同理，我们可以得到三维情况的二阶单向波动方程 (5-12) 5.1.2 Mur有限差分方案 考虑(5-11)，以网格边界为例。波从区域入射到的边界，设表示位于处并与Yee网格边