一种粗网格快速时域有限差分方法.pdf

《一种粗网格快速时域有限差分方法.pdf》由会员分享，可在线阅读，更多相关《一种粗网格快速时域有限差分方法.pdf（16页完整版）》请在专利查询网上搜索。

1、(10)申请公布号 CN 103514143 A (43)申请公布日 2014.01.15 CN 103514143 A (21)申请号 201310433245.9 (22)申请日 2013.09.22 G06F 17/11(2006.01) (71)申请人 西安交通大学 地址 710049 陕西省西安市咸宁西路 28 号 (72)发明人 陈娟 王建国 张安学 田春明 (74)专利代理机构 西安通大专利代理有限责任 公司 61200 代理人 汪人和 (54) 发明名称 一种粗网格快速时域有限差分方法 (57) 摘要 本发明公开了一种粗网格快速时域有限差分 方法， 属于电磁场数值计算领域， 本。

2、发明方法的时 间步长 t 与空间网格长度 z 无关， 空间网格 长度x只需小于等于模拟频段最小波长的1/2。 本发明能够同时解决传统时域有限差分方法的两 大限制条件， 即 Courant-Friedrich-Levy 时间稳 定性条件和空间离散间隔限制条件， 能够在降低 时间稳定性条件的同时， 改善波长对空间网格长 度的限制， 本发明方法适用于模拟同时具有精细 结构和电大尺寸结构的复杂目标， 相比于传统时 域有限差分方法， 该方法具有计算效率高、 计算所 需内存少两大优点。 (51)Int.Cl. 权利要求书 3 页 说明书 9 页 附图 3 页 (19)中华人民共和国国家知识产权局 (12)。

3、发明专利申请 权利要求书3页 说明书9页 附图3页 (10)申请公布号 CN 103514143 A CN 103514143 A 1/3 页 2 1. 一种粗网格快速时域有限差分方法， 其特征在于， 包括以下步骤 ： 1） 对待求电磁目标模型进行空间离散 ： 磁场节点和电场节点的空间排布采用 Yee 元 胞， 电场节点 Ex、 Ey和 Ez位于元胞的各个棱上， 磁场节点 Hx垂直于元胞的 yz 平面， 磁场节 点 Hy与电场节点 Ez的空间位置重合， 磁场节点 Hz与电场节点 Ey的空间位置重合 ； 2） 对待求电磁目标模型进行时间取样 ： 电场分量的时间步取值为 n 和 n+1 时刻， 磁。

4、场 分量的时间步取值也为 n 和 n+1 时刻 ； 3） 对 Maxwell 方程中的空间求导项采用混合时 间步法进行时间离散 ； 4） 对得到的求导项采用傅立叶变换求解， 其 余空间求导项采用二阶中心差分近似 ； 5） 利用公式 （1） 和 （2） 求解 n+1 时刻的电场分量和磁场分量 6） 利用公式 （3） 求解 n+1 时刻的电场分量 权 利 要 求 书 CN 103514143 A 2 2/3 页 3 7） 利用公式 （4） 求解 n+1 时刻的电场分量 8） 利用公式 （5） 和 （6） 求解 n+1 时刻的磁场分量和 权 利 要 求 书 CN 103514143 A 3 3/3 。

5、页 4 以上各式中， 表示傅立叶变换，表示逆傅立叶变换 ； 9） 令 n=n+1， 重复执行步骤 5） 8） 直至迭代完成。 2. 根据权利要求 1 所述的一种粗网格快速时域有限差分方法， 其特征在于， 步骤 3） 所 述的对 Maxwell 方程中的空间求导项采用混合时间 步法进行时间离散 ； 所述的混合时间步法是指在对场量 Ex、 Ey、 Hx和 Hy， 在时刻的选取上， 包 含了未知的 n+1 时刻。 3. 根据权利要求 1 所述的一种粗网格快速时域有限差分方法， 其特征在于， 时间步长 t 与空间网格长度 z 无关， 时间稳定性条件满足 ：空间离 散间隔条件满足 ： 沿 x 方向的空间。

6、网格长度 x 只需小于等于模拟频段最小波长的 1/2， x /2。 权 利 要 求 书 CN 103514143 A 4 1/9 页 5 一种粗网格快速时域有限差分方法 技术领域 0001 本发明属于电磁场数值计算领域， 具体涉及一种粗网格快速时域有限差分方法。 背景技术 0002 目前， 随着科学技术的发展， 电磁波频谱的应用范围已从微波波段扩展至光波波 段， 器件尺寸也缩小至纳米级范围。越来越多的电磁目标呈现出同时具有精细结构和电大 尺寸结构的复杂特征， 如太赫兹波段下的大规模集成电路器件、 大型复杂载体上的发射天 线、 核电磁脉冲作用下的电子信息系统以及薄层涂覆电大目标等等。这些同时具有。

7、精细结 构和电大尺寸结构的复杂目标， 对电磁场数值计算方法的计算速度以及计算机的内存提出 了新的挑战。 0003 众所周知， 时域有限差分方法是当今计算电磁学中应用最为广泛的数值模拟 方法之一。然而， 采用该方法来模拟同时具有精细结构和电大尺寸结构的复杂目标， 却 面临着很大的困难。首先， 该方法需满足 Courant-Friedrich-Levy 时间稳定性条件 ： 即时间步长 t 受空间最小网格尺寸的限制。该限制条 件使得时域有限差分方法对具有精细结构的时域电磁问题的模拟存在计算效率较低的问 题 ； 其次， 为了减小差分近似所带来的数值色散误差， 该方法还需满足空间离散间隔限制条 件， 即。

8、空间网格长度不得大于模拟频段最小波长的 1/10， 该限制条件使得时域有限差分方 法对电大尺寸结构的模拟存在内存需求过大的问题。针对精细结构问题， 可采用弱条件稳 定时域有限差分方法和交变方向隐式时域有限差分方法进行解决。 而对于电大尺寸目标的 电磁模拟， 目前， 也有一些方法可以解决。如 M.Krumpholz 结合小波技术提出的时域多分辨 小波方法和柳清伙教授提出的伪谱时域差分方法。虽然经过研究者们的不断努力与完善， 具有精细结构的时域电磁模拟和电大尺寸目标的模拟问题都分别得到了解决， 但是， 如果 模拟目标同时具有精细结构和电大尺寸， 则无论是传统时域有限差分方法， 还是弱条件稳 定时域。

9、有限差分方法， 以及时域多分辨小波方法或伪谱时域差分方法均没有足够的分析能 力。 发明内容 0004 为了解决上述现有技术存在的缺陷， 本发明的目的在于提供了一种粗网格快速时 域有限差分方法， 该方法能够在降低时间稳定性条件的同时， 改善波长对空间网格长度的 限制， 计算效率高、 内存需求低。 0005 本发明的目的是通过以下技术方案来实现的 ： 0006 一种粗网格快速时域有限差分方法， 包括以下步骤 ： 0007 1） 对待求电磁目标模型进行空间离散 ： 磁场节点和电场节点的空间排布采用 Yee 元胞， 电场节点 Ex、 Ey和 Ez位于元胞的各个棱上， 磁场节点 Hx垂直于元胞的 yz 。

10、平面， 磁场 节点 Hy与电场节点 Ez的空间位置重合， 磁场节点 Hz与电场节点 Ey的空间位置重合 ； 说 明 书 CN 103514143 A 5 2/9 页 6 0008 2） 对待求电磁目标模型进行时间取样 ： 电场分量的时间步取值为 n 和 n+1 时刻， 磁场分量的时间步取值也为 n 和 n+1 时刻 ； 0009 3） 对 Maxwell 方程中的空间求导项采用混 合时间步法进行时间离散 ； 0010 4） 对得到的求导项采用傅立叶变换求 解， 其余空间求导项采用二阶中心差分近似 ； 0011 5） 利用公式 （1） 和 （2） 求解 n+1 时刻的电场分量和磁场分量 0012。

11、 0013 0014 6） 利用公式 （3） 求解 n+1 时刻的电场分量 0015 说 明 书 CN 103514143 A 6 3/9 页 7 0016 7） 利用公式 （4） 求解 n+1 时刻的电场分量 0017 0018 8） 利用公式 （5） 和 （6） 求解 n+1 时刻的磁场分量和 说 明 书 CN 103514143 A 7 4/9 页 8 0019 0020 0021 以上各式中， 表示傅立叶变换，表示逆傅立叶变换 ； 0022 9） 令 n=n+1， 重复执行步骤 5） 8） 直至迭代完成。 0023 步 骤 3）所 述 的 对 Maxwell 方 程 中 的 空 间 求。

12、 导 项 采用混合时间步法进行时间离散 ； 所述的混合时间步法是指在对场量 Ex、 Ey、 Hx和 Hy， 在时刻的选取上， 包含了未知的 n+1 时刻。 0024 时 间 步 长 t 与 空 间 网 格 长 度 z 无 关，时 间 稳 定 性 条 件 满 足 ： 空间离散间隔条件满足 ： 沿 x 方向的空间网格长度 x 只需小 于等于模拟频段最小波长的 1/2， x /2。 0025 与现有技术相比， 本发明具有以下有益效果 ： 0026 本发明在对同时具备精细结构和电大尺寸的复杂模型进行模拟分析时， 能够同 时解决传统时域有限差分方法的两大限制条件， 即 Courant-Friedrich。

13、-Levy 时间稳定性 条件和空间离散间隔限制条件， 能够在降低时间稳定性条件的同时， 改善波长对空间网格 长度的限制 ； 本发明方法的时间步长 t 与空间网格长度 z 无关， 时间稳定性条件为 ： 同时， 本发明方法沿 x 方向的空间网格长度只需小于等于模拟 频段最小波长的 1/2 ： x /2。本发明具有以下优点 ： 0027 1、 本发明提供的粗网格快速时域有限差分方法， 其时间步长 t 不受空间网格长 度 z 的限制， 因此， 在模拟沿 z 方向具有精细结构的时域电磁问题时， 相比于传统时域有 限差分方法， 计算效率大大提高。 说 明 书 CN 103514143 A 8 5/9 页 。

14、9 0028 2、 本发明提供的粗网格快速时域有限差分方法， 空间网格长度 x 只需小于等于 模拟频段最小波长的 1/2。因此， 在模拟沿 x 方向具有电大尺寸结构的时域电磁问题时， 相 比于传统时域有限差分方法， 计算所需内存大大减少。 0029 3、 本发明提供的粗网格快速时域有限差分方法， 在模拟同时具有精细结构 （沿 z 方向） 和电大尺寸结构 （沿 x 方向） 的复杂目标时， 相比于传统时域有限差分方法， 具有计算 效率高、 计算所需内存少两大优势。 附图说明 0030 图 1 为本发明所述的电场节点和磁场节点的空间排布示意图 ； 0031 图 2 为本发明的流程图 ； 0032 图。

15、 3 为本发明具体实施例提供的开缝金属板结构图 ； 0033 图 4 为本发明具体实施例提供的开缝金属板的透射场。 具体实施方式 0034 下面结合具体的附图及实施例对本发明做进一步的详细说明， 所述是对本发明的 解释而不是限定。 0035 本发明提供的粗网格快速时域有限差分方法， 对待求电磁目标模型进行空间离散 时， 磁场节点和电场节点的空间排布采用 Yee 元胞， 各节点的空间排布示意图如图 1 所示， 电场节点 Ex、 Ey和 Ez的位于元胞的各个棱上， 磁场节点 Hx垂直于元胞的 yz 平面， 磁场节点 Hy与电场节点 Ez的空间位置重合， 磁场节点 Hz与电场节点 Ey的空间位置重合。

16、。 0036 本发明提供的粗网格快速时域有限差分方法， 对待求电磁目标模型进行时间取样 时， 电场分量的时间步取值为 n 和 n+1 时刻， 磁场分量的时间步取值也为 n 和 n+1 时刻。 0037 本发明提供的粗网格快速时域有限差分方法对 Maxwell 方程中的空间求导 项采用混合时间步法进行时间离散 ； 对得到的 求导项采用傅立叶变换求解， 其余空间求导项采 用二阶中心差分近似 ； 0038 本发明提供的粗网格快速时域有限差分方法采用以下步骤进行求解 ： 0039 1） 利用公式 （1） 和 （2） 直接求解电场分量和磁场分量 0040 0041 说 明 书 CN 103514143 。

17、A 9 6/9 页 10 0042 2） 利用公式 （3） 和 （4） 求解电场分量和 0043 0044 说 明 书 CN 103514143 A 10 7/9 页 11 0045 3） 利用公式 （5） 和 （6） 求解磁场分量和 0046 0047 0048 以上各式中， 表示傅立叶变换，表示逆傅立叶变换 ； 0049 4） 最后， 令 n=n+1， 重复执行步骤 1） 3） 直至迭代完成。 0050 以上， 可以得出 ： 时间步长 t 与空间网格长度 z 无关， 时间稳定性条件满足 ： 空间离散间隔条件满足 ： 沿 x 方向的空间网格长度 x 只需小 说 明 书 CN 10351414。

18、3 A 11 8/9 页 12 于等于模拟频段最小波长的 1/2， x /2。 0051 上述实施步骤的流程图如图 2 所示。 0052 下面以图 3 所示的开缝金属板为例， 说明粗网格快速时域有限差分方法在计算效 率和计算所需内存方面的优势。 0053 图 3 所示的金属板尺寸为 60mm1mm， 板上开有三条细长缝隙， 尺寸均为 30mm0.2mm， 各个缝之间的距离为 0.1mm， 一沿 x 方向的电流源位于金属板正前方 1.5mm 处， 频率为 100GHz， 波长为 3mm。该模型沿 z 方向具有精细结构， 沿 x 方向具有电大尺寸结 构。 0054 采用本发明的粗网格快速时域有限差。

19、分方法计算金属板后的透射场， 所得结果如 图 4 所示。为了便于比较， 图中同时给出了传统时域有限差分方法的计算结果， 其中， 实线 表示传统时域有限差分方法的计算结果， 虚线表示粗网格快速时域有限差分方法的计算结 果。从该图可以看出， 两种方法的计算结果符合较好。 0055 完成上述模拟， 两种方法所用的空间网格尺寸 x， 时间步长 t， 计算时间和计 算所需内存如下表所示 ： 0056 0057 以上分析结果表明， 粗网格快速时域有限差分方法具有与传统时域有限差分方法 相似的计算精度， 但由于采用了较大的空间网格 x 和较大的时间步长 t， 其计算所需时 间、 所需内存均远小于传统时域有限。

20、差分方法。 0058 对 Maxwell 方程中的求导项在时刻的选取 上采用混合时间步技术， 以此消除 z 方向空间网格长度 z 对时间步长 t 的限制， 使时间 稳定性条件变为 ： 0059 0060 对 Maxwell 方程中的求导项采用傅立叶变 换求解， 而不是传统时域有限差分方法的中心差分近似， 该方法在理论上可保证沿 x 方向 一个波长仅用 2 个网格来离散。 0061 综上所述， 本发明提供的粗网格快速时域有限差分方法， 通过对 Maxwell 方程中 的求导项在时刻的选取上采用混合时间步技术， 对 说 明 书 CN 103514143 A 12 9/9 页 13 求导项采用傅立叶。

21、变换求解， 设计出了一种时间步 长t与空间网格长度z无关， 空间网格长度x只需小于等于模拟频段最小波长1/2的 新型时域有限差分方法。该方法在模拟同时具有精细结构和电大尺寸结构的复杂目标时， 相比于传统时域有限差分方法， 具有更高的计算效率和更低的内存需求。 0062 以上内容是结合具体的优选实施方式对本发明所做的进一步详细说明， 不能认定 本发明的具体实施方式仅限于此， 对于本发明所属技术领域的普通技术人员来说， 在不脱 离本发明构思的前提下， 还可以做出若干简单的推演或替换， 都应当视为属于由本发明所 提交的权利要求书确定专利保护范围。 说 明 书 CN 103514143 A 13 1/3 页 14 图 1 说 明 书 附 图 CN 103514143 A 14 2/3 页 15 图 2 说 明 书 附 图 CN 103514143 A 15 3/3 页 16 图 3 图 4 说 明 书 附 图 CN 103514143 A 16 。