FDTD的PML边界条件.docVIP

FDTD的PML边界条件 　　周期结构 FDTD 算法的各向异性 PML 边界条件 　　康利鸿 1 　　, 刘克成 　　2 　　(1. 北京遥感信息研究所 北京 , 100011; 2. 长沙国防科大 403教研室 , 长沙 410073) 　　摘 　 要 :　 本文建立了一种新的适用于周期 FDTD 计算 的边界条件 — 周 期 /各 向异性 P M L 混合边 界条件 . 计 算验 　　证表明 , 本文的方法改善了原来单向波边 界的缺点 , 具有较高的精度 , 使周期 FDTD 算法进一步走向成熟 . 　　关键词 :　 周期结构 ; FDTD ; 边界条件 　　中图分类号 :　 TN015　　　 文献标识码 :　 A 　　　 文章编号 :0372-2112(2000) 12-0111-02 　　Anisotropic PML B . C . in FDTD Calculation of Periodic Structures 　　KANG Li -hong 1, LIU Ke -cheng 2 　　(1. Beij in Remote Sens ing Institute , Beij ing 100011, China ; 　　2. Dept . of El ectr onic Te chnical , National Univer sity De fenc e of Tec hnic al , C hangsha 410073, C hina ) 　　Abstract :　 In this paper , a new algorithm , periodic /Anisotropic P M L B . C . for periodic FDTD calculation , is introduced . The re -sult show that the new method is better in precision than the old one and eliminate the shortcoming in the old . This work makes period -ic FDTD method more convenient and accurate in use . 　　Key words :　 periodic structures ; FDTD ; boundary condition 　　1　 引言 　　在处理电磁问题时 , 经常遇到周期结构问题 , 即同样的结 构向两侧无限延伸 . 周期 结构的 电磁计 算一直 得到广 泛的关 注和研究 [1~4] 　　, 研究周期结构算法有助于简化目标的计 算 , 得 到更精确的结果 . 　　最早对周期结构的计算用频域 方法 :利用 Floquet 模的原 理 , 建立起积分方程进 行求解 . 为了 得到时 域的响 应 , 需 要求 出大量频率 点上的值进行 傅里叶变 换 , 这是非 常低效 率的方 法 . 1991年 , C . H . Chen 等把 FDTD 方法 引入周 期结 构电 磁问 题 [5], 求解时域 响 应 , 得 到 了 令 人满 意 的 结 果 . 图 1是周 期 FDTD 算法的计算周期结构电磁问题的示意图 , 只取了无限周 期中的一个进行计算 , 边界情况如图 1所示 . 　　同期结构 FDTD 算法 和非周 期结 构的 FDTD 算法有 明显 的差异 . 在用周期 FDTD 算法解周期问题 时 , 为 了从已 知场推 出未知场的迭代式能够顺 利进行 , 需 要用 E 、 H 的替换量来代 替 E 、 H 本 身参 加迭 代 , 并把 原来 的方 程改 变为 替换 量 的方 程 . 原来对应 E 、 H 建立的 吸收边界条件要改为 替换量的吸收 边界条件 . 而能否找到新的适于周期算法的边界条件 , 是周期 FDTD 算法所面临的问题 , 新边 界条件引起误差的大小将决定 算法的精度 . 文献 [6]中对周期边界的处理是采用 变形的单向 波边界条件 , 它具有一定的精 度 , 但是 如同经 典 FDTD 方法中 的单向波边 界一样 , 精 度有限 , 尤其 在三维 问题中 , 有时 会引 起不小的边 界反 射 , 给计 算带 来很 大误 差 , 甚至 于使 计 算发 散 . 这在很大程度上限制了周期算法在实际中 的应用 . 　　本文所 做 的 工作 就 是根 据 FDTD 算 法 中各 向 异性 PML 　　(完全匹配层 ) 边界及 周期算 法的原理 , 对 原来的 周期迭 代式 进行 分裂 、 变形 后 , 得 到了 周期 /各 向异 性 PML 混 合 边界 条 件 . 并在文章的最后进行了数值验证 . 　　2