WE—FDTD中近似吸收边界条件的角点计算.pdfVIP

2014年第4期 文章编号：1009—2552(2014)04—0042—02 中图分类号：TN911 文献标识码：A WE—FDTD中近似吸收边界条件的角点计算 赵长青 (河海大学数理部，江苏 常州213022) 摘 要：从线性插值思想出发研究波动方程时域有限差分法 (WE—FDTD)中二维二阶近似吸 收边界条件差分格式下矩形网格空间角点的场量计算方法。提出确定线性插值点位置的原则和 线性插值点场量的迭代方法，考虑衰减因素进而计算角点的场量。 关键词：波动方程时域有限差分法；线性插值；近似吸收边界条件；角点；差分格式 Thecornercalculation0fW E．FDTD inabsorbingboundary condition ZHA0 Chang—qing (DepartmentofMathematicalPhysics，HohmUninersity，Changzhou213022，JiangsuProvince，China) Abstract：Thispaperstartsbasedonthewaveequationfinite—differencetime—domainme~odfrom linear interpolationme~odintwodimensionalorderabsorbingboundaryconditiondifferencefieldrectangular gridspatialformatanglecalculationme~od．Itpresentstodeterminetheprincipleoflinearinteprolation pointandlinearinteprolationmethod，consideringtheattenuationfactor． Keywords：waveequationthefinite—differencetime—domainmethod；linearinteprolation；approximation oftheabsorbingboundary conditions；comer；differencescheme 0 引言 有限的网格空间来模拟电磁波在无界空间中的传 由波动方程时域有限差分法(WE—FDTD)的基 播。加于边界场的适合上述要求的算法称为吸收边 本原理可知，这种算法的一个重要特点是，在需要计 界条件(AbsorbingBoundaryCondition，ABC) 。 算电磁场的全部区域建立Yee氏网格计算空间。于 1 二维精确吸收边界条件 是对于像辐射、散射等这类开放问题，所需要的网格 如图1所示，在平面二维问题中，左右边界的精 空间成为无限大。任何计算机的存储空间都是有限 确吸收边界条件为 j： 的，因此在无限大网格空间中计算电磁场是根本不 可能的。在实际计算 中总是在某处把网格空间截 断，使之成为有限的。这样一来，在网格空间的截断 处就会出现非物理的电磁波的反射，这将严重地影 响计算的精度，必须设法消除这种反射现象 j。另 一 方面，中心差分格式的时域有限差分方程 由于需 要截断处边界外场的信息用于边界网点上场的计 图1 二维矩形计算网格空间 算，故也需要适合于截断处边界网格点计算的算 法 J。总体来讲，需要一种截断边界网格点处场的 收稿 日期：2013—05—20 基金项 目：国家 自然科学基金项 目；中央高校基本科研 特殊计算方法，它不仅能保证边界场计算的必要精 业务费专项资金项 目(2010B23414) 度，而且还能大大消除非物理因素引起的入射到截