PML吸收边界条件下TM波时域有限差分法分析.pdfVIP

PML 吸收边界条件下TM 波的时域有限差分法分析 1 2 2 马晓晖 ，刘丽卓 ，赵维 1.北京邮电大学光通信中心，北京（100876） 2.陕西渭南师范学院计算机科学系，陕西渭南（714000） E-mail ：mxhui123@ 摘 要：介绍了时域有限差分法的基本原理，通过场区的FDTD 公式分析了PML 边界及其 相关参数设置的要求。最后以自由空间中TM 波传播为实例，利用C 语言编写了FDTD 程 序。并对结果进行了详细分析。 关键词：完全匹配层，TM波，时域有限差分法 中图分类号：TN252 1. 引 言 目前，电磁场的时域计算方法越来越引人注目。时域有限差分（Finite Difference Time Domain，FDTD ）法作为一种主要的电磁场时域计算方法，最早是在 1966 年由 K. S. Yee 提出的。这种方法通过将 Maxwell 旋度方程转化为有限差分式而直接在时域求解，通过建立 时间离散的递进序列，在相互交织的网格空间中交替计算电场和磁场。经过三十多年的发展， [1] 这种方法已经广泛应用到各种电磁问题的分析之中 。 时域有限差分法的主要思想是把Maxwell 方程在空间、时间上离散化，用差分方程代替 一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。FDTD 空间网格单元上电 场和磁场各分量的分布如图 1 所示。 图 1 FDTD 网格中电磁场分量的分布 　 电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元 每个面的中心，每个磁场（电场）分量都有 4 个电场（磁场）分量环绕。这样不仅保证了介 质分界面上切向场分量的连续性条件得到自然满足，而且还允许旋度方程在空间上进行中心 差分运算，同时也满足了法拉第电磁感应定律和安培环路积分定律，也可以很恰当地模拟电 [2] 磁波的实际传播过程 。 2. 场区 FDTD 公式的讨论 ∂E z ∂H x ⎫ − − m x ∂y µ ∂ t σ H ⎪⎪ ∂E z ∂ H y + m y ⎪ TM 波 (1) µ σ H ⎬ ∂x ∂ t ⎪ ∂H y ∂ H x ∂ E z ⎪ ∂x − ∂y ε ∂ t + σ E z ⎭⎪ - 1 - 对于 TM 波，Ex = Ey = Hz =0 ，FDTD 公式为： n ( , 1) n ( , ) n +1/ 2 n −1/ 2 E z i j + −E z i j