交错紧凑隐式差分法的数值色散分析.pdfVIP

第 5 期 第 32卷 第 5期 佛 山科 学技术 学 院学报 （自然科 学版 ） V o l.3 2N o.5 2 0 14 年 9 月 JournalofFoshanUniversity （NaturalSciencesEdition ） S e p.2 0 14 文章编号 ：1008-0 17 1（20 14 ）05-0051-04 交错紧凑隐式差分法的数值色散分析 陈方芳 （黎 明职业大学 公共教 学部 ，福建 泉 州 362000 ） 摘要 ：首先使用空间的四阶交错 紧凑 隐式差分和时间交错后 向时间差分 ，提 出时域微分麦克斯韦方程组 的交 错紧凑隐式差分(CSBD-FDTD)方法 ，给出直角坐标系下 CSBD-FDTD 方法 的场量迭代公式 ，然后使用数学分析 方法推导 CSBD-FDTD 方法稳定性条件和数值色散关系 ，为电磁计算提供新 的途径 。 关键词 ：交错紧凑隐式差分 ；数值色散 ；稳定性条件 中图分类号 ：TN0 11 文献标志码 ：A 目前 ，关于高阶精度时域有 限差分法方面的研究得到了广泛 的关注 。 文献 ［1］提出一个新 的 SBD 理论 ，并用 SBD 理论描述时间差异 ，将空间差异与四阶中心差分近似相 比较 ，然而 由于涉及到高阶的特 征多项式 ，它很难分析理论 的稳定性与色散性。 文献 ［2-3 ］则用复杂的 Routh-Hurwitz（赫维茨 ）理论来解 [4] 决这个问题。 时域有限差分法（FDTD ）是解决电磁问题非常有效的一种数值方法 ，然而 FDTD 法属于显 性 的差分方程式 ，必须满足 Courant Friedrich Levy（CFL ）稳定性条件 ，时间步长受到限制 ，使得计算效率 不高。 1 交错紧凑隐式差分法的基本思想 为简单起见 ，二维 TE 波模式下无源、均匀 、无耗媒质 中的麦克斯韦方程可定义如下 鄣Ｖ ＝Ｒ×Ｖ ， （1） 鄣ｔ ０ 鄣 ０ ０ ０ ０ ０ 着鄣ｙ Ｅ ０ ｘ ０ 鄣 ０ Ｅ Ｒ＝ ０ ０ － 其 中， Ｖ ＝ ｙ ， ０ ０ 着鄣x ０ Ｈ ０ ｚ ０ 鄣 鄣 ０ － ０ ０滋 鄣ｙ 滋 鄣ｘ