FDTD办法中吸收边界条件.doc

? ? ? ? ? ? ????如果你遇到不懂的专业词汇，请您首先用鼠标选中该词汇，然后按住鼠标左键将该词汇拖入查询框中，再点击“查询”，就会弹出窗口显示该词汇的专业注释！????????????? HYPERLINK javascript:close() 关闭 ? FDTD方法中的吸收边界条件 胡来平，刘占军 （重庆邮电学院光电工程学院　重庆　400065） 　　摘　要：介绍并分析了时域有限差分中的吸收边界条件，对各种条件的应用进行了讨论，对时域差分技术的吸收边界条件进行了一定的总结和展望。??? 关键词：时域有限差分方法；吸收边界条件；电磁散射；完全匹配层 　　时域有限差分法（FDTD）是一种分析各种电磁问题的全波方法。用FDTD分析电磁辐射、散射等开放或者半开放性质问题时，不可能直接对无限的结构进行计算，因此必须在截断处设置适当的吸收边界条件，以便用有限网格空间模拟开放的无限空间或无限长的传输结构。理想的吸收边界条件应在截断边界上只有向外传输的波而没有向内的反射波。??? 自从Yee提出FDTD方法以来，对FDTD方法中的重要组成部分--吸收边界条件的研究就一直没有停止过。目前，构造吸收边界条件的思路主要有2种：一种是在边界上引入吸收材料，电磁波在无反射地进入吸收材料后被衰减掉，如PML。这种方法构造复杂，内存需求较大，但在很大的入射角度上吸收效果较好。另一种是从外行波方程出发构造的透射边界条件，如Mur边界条件等。这种类型的透射边界条件具有构造简单，内存需求小，基本上不额外消耗内存等特点。下面介绍几种应用较为广泛的吸收边界条件。1　Mur吸收边界条件［1］??? 考察一维波动方程：??? 他可分解为2个单向波方程：　　当边界上电磁场满足式（2）时，电磁场仍是单向波形式，不产生反射，这就是Mur一阶吸收边界条件。同法对二维情况，有二维波动方程：??? 把式（4）根号部分进行Taylor展开，然后取其前2项，即令：　　这就是Mur所建议的具有二阶近似的，适用于二维问题的近似吸收边界条件。他在FDTD中有广泛地应用。Mur吸收边界条件具有实施方便简单、吸收边界条件效果好的特点，然而在使用中注意到，一阶近Yee网格划分，在角区域存在较大误差，而二阶近似尽管就算精度较高，但编程复杂，且对三维情况还可能出现结果发散的现象［2］。2　廖氏吸收边界条件　　廖氏吸收边界条件比同阶的Mur吸收边界条件反射小约一个数量级，并且各阶吸收边界条件可用统一的公式表示。由于推导繁琐，这里直接给出其吸收边界条件公式：其中：CjN为组合数，N表示廖氏吸收边界条件的阶数。??? N＝1时，给出了一阶吸收边界条件：其中：S＝αvΔt／Δx，α是一个常数，0．5≤α≤2。　　N＝2时，给出了二阶吸收边界条件：3　超吸收边界条件　　Fang和Mei于1988年第一次提出了超吸收边界条件的概念。传统的吸收边界条件只在边界上给电场或磁场进行特殊处理，不同时计算二者，因为只要在边界知道了电场或磁场，则内部区域的场就能惟一地确定。在超吸收边界条件中让磁场也参与计算并用他来减少计算电场时所产生的非物理因素引起的反射，从而改善原始吸收的性能。其具体的理论推导可见文献［1］。超吸收边界条件的具体运用格式为：　　（1）在自由空间的截断边界上采用一种吸收边界　　（5）利用得到的，M－0．5，．），重新由前面的差分方程中H分量在n＋0．5时刻的表达式计算出）4　完全匹配层　　1994年J．P．Berenger提出了“完全匹配层（PML）”这种新边界，他是应用某种依赖于方向的、满足匹配条件的地导电和导磁媒质来吸收反射波。通常在计算域的截断面之外布置完全匹配层，当波进入PML中时，因波阻抗保持不变而无反射发生。当外行波到达PML的最外层时，其幅度近似衰减为0。微弱的反射波在此产生，反射波传播到计算域之前，又要经过完全匹配层的再次衰减，这样进入计算域的反射波是十分微弱的。完全匹配层是一种非物理性的电磁波吸收层，用电导率σ和磁导率σ＃表征，σ和σ＃可分别解为σx，σy，σz和。可以证明，当PML满足以下2个条件时，他对电磁波不反射［3］：　　（1）为吸收某一方向的电磁波（如z向），则σ和σ＃在其他方向上的分量均为0。??? （2）σ和σ＃与电介常数ε和磁导率μ满足下列方程：5　吸收边界条件的现状和一点展望　　时域有限差分法是最受瞩目的电磁场数值计算方法之一，现在应用于包括电磁辐射、散射等在内的几乎所有电磁问题的数值模拟中。他实质上是一种蛙跃法，具有二阶精度。　　对于开域问题，受计算机内存和计算时间限制，必须截断计算空间并设置边界条件，国内外许多人在这方面做了大量的工作，提出了各种边界条件，其目的是使有限的