computational electrodynamics The Finite Difference Time Domain Method 3rd edition.docVIP

完全匹配层吸收边界条件 Stephen Gedney 7.1 介绍 无界空间中电磁波相互作用问题有效且准确的解决方案一直是FDTD（时域有限差分）方法最大的挑战。对于这些问题，必须引入一种位于外层网格边界的吸收边界条件（ABC），以仿真至无限远的网格的存在。正如我们在第六章回顾过的，已经引入大量的分析技术手段来实现这一目标。 一种实现ABC的代替方案是在吸收材料媒质中的空间网格的外层边界处终结。这与吸波暗室墙壁的物理方法类似。理想情况下，吸收媒质只有几个网格胞元的厚度，对传播的波在其全频谱无反射，高吸收，且在源或散射体近场处效率高。早期的一个实现这样的吸收材料边界条件的尝试是由一荷兰人提出的，他利用了一种传统的有损耗无色散的吸收媒质。此方式的困难在于这样一种吸收层只能与垂直入射的平面波相匹配。因此，有损耗ABC的整个范畴在当时的应用仅仅被限定在电磁领域内。 自1994年起，Berenger提出的一种其命名为完全匹配层（PML）的高效吸收材料的ABC在这个领域中掀起一股新的热潮。Berenger的PML的创新在于任意入射角度、极化和频率的的平面波都可以在边界处匹配。或许具有同等重要的意义的是PML也可以被用作为一种吸收边界来终结由不均匀、色散、各向异性以及甚至非线性媒质组成的作用域，这在之前仅由解析导出的ABC是不可能实现的。 在他开创性的工作中，Berenger得到了一个新的分裂场形式的Maxwell方程组，在这个方程组中每个矢量场分量被分为两个正交分量。Maxwell的旋度方程也是通过适当地分离，得到一组十二元耦合的一阶偏微分方程组。之后，通过选择与非色散媒质相一致的损耗参数，便得到了一个完全匹配的平面。 在连续空间中，PML吸收体和基本媒质是完全匹配的。然而，在离散FDTD网格中，电与磁的材料参数是由一种分段常数的方式来表示并且在空间上交错的。这导致了能够降低PML的理想表现的离散化误差。Berenger提出了一种PML参数的空间尺度来减少与材料交界面有关的离散化误差。他也提出了减少对于二阶FDTD算法的离散化误差以及减小PML需要厚度的方法。 随着Berenger性的，许多论文证实了他的研究，并应用了使用PML媒质的FDTD方法。为了强化PML的表现，人们也对其进行了一些修改。原始的分离场PML概念也以一种延伸坐标的形式进行了重申。这种形式拓展了PML对于通过度量的坐标映射的其他正交坐标系统的应用。分离场PML也能够直接用于FDTD离散化，这需要使用一种叫做卷积完全匹配层（CPML）的极为有效的算法，这种算法是由Roden和Gedney提出的。CPML适应了更普遍的能够使对缓慢变化的衰逝波的吸收提高的度量张量系数。 也有观点将分离场PML作为一种同时具有电、磁介电常数张量的单轴各向异性媒质。将损耗引入张量，得到了一种完全匹配的吸收媒质。单轴完全匹配层（UPML）最早是由Gedney在FDTD方法中应用的。UPML吸引人的地方在于它是以Maxwell方程而非数学模型为基础的。事实上，实现物理上的UPML媒质的尝试已经有人提出过了。 这一章将会回顾Berenger的分离场PML、它的延伸坐标形式、UPML以及CPML。本书将着重讲解这些方法间的关系。针对对笛卡尔坐标中的FDTD空间网格的UPML和CPML，本书会进行有效的补充。附加的内容包括了多种PML的多离散化误差和本构参数与空间分段函数的最优选择。我们发现PML在含有通常情况会从空间网格终止的一些均匀、导电以及色散媒质区域的FDTD模型中是可以应用的。 7.2 平面波垂直射入部分有损耗空间 为了建立一个我们讨论PML吸收体的基础，这部分回顾了一个正弦平面波以相对轴任意角度射入一个常规的有损耗材料的物理过程。特别地，我们考虑一波由无损耗材料区域1射入有损耗材料区域2，区域1为，区域2（）材料的电导率与磁损耗参量分别为、。故，入射波磁场为 (7.1) 这里的上圆弧表示矢量。则，区域1内的总场可以表示为， (7.2a) (7.2b) 入区域2内的场分量可以表示为 (7.3a) (7.3b) 、分别为磁场的反射、系数， (7.4b) 其中中，i=1,2。为了加强外围场区在区域1与区域2于交界x=0处的连续性，使，以及 (7.5a,b) 总之，对应任意入射角度。然而，对于垂直入射=0的特殊情形有， (7.6) 其中区域1与区域2的波阻抗可表示为 (7.7a,b) 然后我们令，更进一步加强条件 (7.8) 故有，。由此可得，即一种垂直入射时区域1与区域2的无反射交界的情形。由(7.4b)，我们也可以找到对应这种情形， (7.9) 得到下面区域2中的透射场量： (7.10a,b) 我们发现，对于垂直