时域有限差分法中几种吸收边界条件的比较和数值验证.pdfVIP

时域有限差分法中几种吸收边界条件的比较与数值验证’ 马双武高攸纲 (北京邮电大学171信箱100876) 擅叠本文对时域有限差分法(FDTD)中的Mur近似吸收边异条件、超吸收技术、PML技术等有影响 的吸收边界备件进行了比较，进行了数值试验，分析和总结了几点结论． 关销≈时域有限差分法吸收边界条错数值试验 1，引言 们的关注并得到了快速的发展．FDTD的优越性能，使它获得了广泛的应用，几乎被用到 了电磁场工程中的各个方面。利用FDTD计算电磁场，需要在全部场域内划分网格并实旖 计算。有些问题．如辐射、散射等，属于开放问题，其场域一般是无限的。由于计算机资 源的有限性，在一个无限大的场域中划分网格并实施计算是根本不可能的。因此，在实际 计算中．必须将网格在某处截断，使之成为有限空间。但这样的截断边界一般会带来对电 磁波的非物理性的反射，从而严重影响计算精度。所以，为了能够利用有限的网格空间模 拟电磁波在无界空间中的传播，必须对截断边界实施某种特殊的算法，这种算法一方面要 使得FDTD算法在有限的网格空间得以实施，另一方面要能最大限度地减少边界对波的反 射，从而提高计算精度。这种作用于截断边界上的特殊算法一般称之为吸收边界条件 (ABC)。 在FDTD的发展中，人们提出了多种吸收边界条件【2】-[10l，其中比较有影响的有【5j 【6】【7】等。Mur(1981)【5】通过对波动方程的因式分解获得了单行波方程，并将之作为 吸收边界条件。由于精确的单行波方程无法实施数值计算，Mur讨论了它的各阶近似算法 及相应的差分格式，其二阶近似ABC是目前应用非常广泛的吸收边界条件。Mei等(1992) 【6]对人们提出的吸收边界条件进行了分析，提出了一种修正的算法(超吸收技术)．提高 PML技术。这一技术很快就得到了重视和发展【8】[9】【lO】。 本文对上述几种有影响的吸收边界条件的吸收性能进行了比较，并进行了数值试验， 在此基础上给出了一些建议和结论。 2．理论简述 1) MIll二阶近似吸收边界条件 以直角坐标系为例。设场量々满足波动方程：Lqb=O 其中 如 + + 一 生妒 生劳 ●一矿 竺酽 将算子三分解为 悦生舻口 ·车工作获国家臼然科学基金资助(项目编号 ．390． ”{ ∥毋翥 其中 p=；±三昙麝 …恸舻(妄／甜 垂直x轴的两边界面x=xmh和#x。：：上的精确吸收边界条件为： L-4,=0 (x嚷．) L++--0(H) 对于垂直Y，z轴的其它边界面，容易得到相应的结果。 由于上述公式中的根式不适合直接数值计算。Mur对之进行了近似处理。其二阶近似边界 条件是在以下近似的基础上得到的： 2 √l—sz1一S212 其二维空间中垂直x轴的二阶近似边界条件为： 苏西v甜2 却2J 一 ft盖丢等十三2甜=。c踽。， ‘～ 2 V L(旦OxOt七罢Ot上2豺=。c—劫2J’ 其它边界面上的二阶近似边界条件可以同样求出。 Mur二阶近似边界条件的反射系数R与波的入射角e有关，为： 2)超吸收技术