时域有限差分中各向异性完全匹配层的实现方法.docVIP

时域有限差分中各向异性完全匹配层的实现方法 　　摘要：本文讨论了如何在给定边界层数的条件下，通过调整各介质层吸收层的介质参数来实现入射波的最佳吸收，并利用自适应遗传算法对算法中的边界参数进行优化。 　　关键词：时域有限差分算法；各向异性；完全匹配层 　　中图分类号：O441文献标识码：A文章编号：1009-3044(2008)20-30234-03 　　 　　The Implementation Methods of PML in the Finite Difference Time Domain 　　SUN Yun-he, LI Jian-xiong 　　(Tianjin Industry University, Information and Communication Engineering Institute, Tianjin 300160, China) 　　Abstract: This article discusses how a given number of boundary conditions, by adjusting the dielectric layer absorption of the medium parameters to achieve the best incident wave absorption and use of adaptive algorithms in the genetic algorithm to optimize the parameters of the border. 　　Key words: FDTD method; anisotropy; PML 　　 　　1 引言 　　 　　时域有限差分算法(FDTD)是一种近年来应用最广泛、发展最迅速的电磁场仿真算法，用该算法计算电磁场散射问题时需要设置吸收边界条件来模拟无限空间中的电波传播。时域有限差分(FDTD)方法已逐渐成为解决电磁散射和电磁波传输等问题的有力工具。研究FDTD方法的核心问题是寻求一种理想的吸收边界，使截断面反射最小。 　　各向异性介质在微波元件、微带电路中得到广泛的应用，特别是随着微波与毫米波集成技术的发展,集成电路由混合集成发展为单片集成，各向异性的介质基片得到越来越广泛的应用。若在电路特性参数的分析中，忽略介质基片的各向异性而作为各向同性介质处理，随着频率的升高，将产生较大的误差。因此，分析各向异性介质问题的数值方法得到人们的重视。完全匹配层PML吸收边界是目前最常用、边界处理比较方便且吸收效果比较理想的一种吸收边界条件。与传统的Mur二阶近似吸收边界条件相比，PML吸收边界条件可提高精度40dB左右。完全匹配层参数的设置合适与否直接影响电磁场模拟仿真的精度，为了提高介质层吸收效果，一般采用增加吸收层的物理厚度。由于吸收层位于网格最外层，因此每增加一层都会增加计算机的内存需求和计算时间。 　　吸收边界的处理是FDTD算法的一个关键问题，吸收边界条件的效果直接关系到FDTD计算的正确性和精确度，是影响FDTD计算品质的决定因素。J.P.Berenger于1994年首先提出了理想匹配层(PML)吸收边界条件的概念。1996年，S.D.Gedney提出了各向异性介质PML理论并用于FDTD区域的吸收边界。与传统的Mur二阶近似吸收边界条件相比[7]，PML吸收边界条件可提高精度40dB左右，是目前最好的吸收边界条件。PML吸收边界条件是针对矢量FDTD算法提出的，并不能直接用于标量的WE-FDTD算法。最近，M.A.Alsunaidi导出了PML媒质中波动方程的表达式，但它忽略了电导率是空间坐标的函数。D.Zhou基于伸张坐标系下的标量波动方程导出了一维PML的差分格式，并用于光器件的仿真中，获得了非常好的结果。本文实现了波动方程FDTD算法的PML吸收边界条件(WE-PML)，并利用二维空间中的TM波进行了数值验证，同时对PML媒质参数进行了优化和分析。 　　 　　2 WE-PML吸收边界条件 　　 　　PML边界是在FDTD区域边界处设置一种特殊介质层，该介质层的波阻抗与相邻真空的波阻抗完全匹配，即在连续和离散情况下都满足阻抗匹配条件。理论上希望垂直入射时入射波能无反射地进入PML区域，即反射系数为零。但PML区域不可能是无限大，因此在区域外围要用截断边界（反射系数为――1）或其他吸收边界。Berenger通过数值实验发现，真空一PML交界面电导率的突变是导致较大数值反射的主要原因。因此，良好PML参数设置的目的是：①尽量用少量的吸收网格使计算区域中入射波能无反射地进入PML区域；②尽量减少吸收网格之间的电导率变化速率；③尽量减少到达截断边界的入射波强度。