FDTD方法中的吸收边界条件

1、如果你遇到不懂的专业词汇，请您首先用鼠标选中该词汇，然后按住鼠标左键将该词汇拖入查询框中，再点击“查询”，就会弹出窗口显示该词汇的专业注释！关闭FDTDFDTD 方法中的吸收边界条件方法中的吸收边界条件胡来平，刘占军胡来平，刘占军（重庆邮电学院光电工程学院（重庆邮电学院光电工程学院重庆重庆400065400065）摘摘要要：介绍并分析了时域有限差分中的吸收边界条件，对各种条件的应用进行了讨论，对时域差分技术的吸收边界条件进行了一定的总结和展望。关键词关键词：时域有限差分方法；吸收边界条件；电磁散射；完全匹配层时域有限差分法（FDTD）是一种分析各种电磁问题的全波方法。用 FDTD分析电磁辐射、

2、散射等开放或者半开放性质问题时，不可能直接对无限的结构进行计算，因此必须在截断处设置适当的吸收边界条件，以便用有限网格空间模拟开放的无限空间或无限长的传输结构。理想的吸收边界条件应在截断边界上只有向外传输的波而没有向内的反射波。自从 Yee 提出 FDTD 方法以来，对 FDTD 方法中的重要组成部分-吸收边界条件的研究就一直没有停止过。目前，构造吸收边界条件的思路主要有 2种：一种是在边界上引入吸收材料，电磁波在无反射地进入吸收材料后被衰减掉，如 PML。这种方法构造复杂，内存需求较大，但在很大的入射角度上吸收效果较好。另一种是从外行波方程出发构造的透射边界条件，如 Mur 边界条件等。这种

3、类型的透射边界条件具有构造简单，内存需求小，基本上不额外消耗内存等特点。下面介绍几种应用较为广泛的吸收边界条件。1 1MurMur 吸收边界条件吸收边界条件1 1考察一维波动方程：他可分解为 2 个单向波方程：当边界上电磁场满足式（2）时，电磁场仍是单向波形式，不产生反射，这就是 Mur 一阶吸收边界条件。同法对二维情况，有二维波动方程：把式（4）根号部分进行 Taylor 展开，然后取其前 2 项，即令：这就是 Mur 所建议的具有二阶近似的，适用于二维问题的近似吸收边界条件。他在 FDTD 中有广泛地应用。Mur 吸收边界条件具有实施方便简单、吸收边界条件效果好的特点，然而在使用中注意到，

4、一阶近 Yee 网格划分，在角区域存在较大误差，而二阶近似尽管就算精度较高，但编程复杂，且对三维情况还可能出现结果发散的现象2。2 2廖氏吸收边界条件廖氏吸收边界条件廖氏吸收边界条件比同阶的 Mur 吸收边界条件反射小约一个数量级，并且各阶吸收边界条件可用统一的公式表示。由于推导繁琐，这里直接给出其吸收边界条件公式：其中：CjN为组合数，N 表示廖氏吸收边界条件的阶数。N1 时，给出了一阶吸收边界条件：其中：Svtx，是一个常数，052。N2 时，给出了二阶吸收边界条件：3 3超吸收边界条件超吸收边界条件Fang 和 Mei 于 1988 年第一次提出了超吸收边界条件的概念。传统的吸收边界条件

5、只在边界上给电场或磁场进行特殊处理，不同时计算二者，因为只要在边界知道了电场或磁场，则内部区域的场就能惟一地确定。在超吸收边界条件中让磁场也参与计算并用他来减少计算电场时所产生的非物理因素引起的反射，从而改善原始吸收的性能。其具体的理论推导可见文献1。超吸收边界条件的具体运用格式为：（1）在自由空间的截断边界上采用一种吸收边界M05，），重新由前，（5）利用得到的面的差分方程中 H 分量在 n05 时刻的表达式计算出）4 4完全匹配层完全匹配层1994 年 JPBerenger 提出了“完全匹配层（PML）”这种新边界，他是应用某种依赖于方向的、满足匹配条件的地导电和导磁媒质来吸收反射波。通常在计算域的截断面之外布置完全匹配层，当波进入 PML 中时，因波阻抗保持不变而无反射发生。当外行波到达 PML 的最外层时，其幅度近似衰减为 0。微弱的反射波在此产生，反射波传播到计算域之前，又要经过完全匹配层的再次衰减，这样进入计算域的反射波是十分微弱的。完全匹配层是一种非物理性的电磁波吸收层，用电导率和磁导率表征，和可分别解为x，y，。可以证明，