空间分割的FDTD方法研究.pdfVIP

电磁兼容论文集 空间分割的FDTD方法 薛正辉 高本庆 北京理工大学电子工程系 1000R 摘要 本文介绍了一种 “空间分割”的FDTD方法，此方法可应用于利用FDTI!方 法分析电大尺寸目标电磁兼容性的工作中、以便有效利用有限的计界机内存资源，达 到常规FDTD方法所不能企及的分析远超出个人计算机容量的目标的目的。文章以 计算实例验证了此方法的可行及计葬结果的正确性C 关健词F -I 电大尺寸目标 价宝;___ 一 、介绍 乍j“L勒 时域有限差分法(简称FDID方法)是 1966年由Yee首先提出的，该方法允许我们研究具有 任意形状和材料组成物体的电磁场分布。目前广泛应用于研究散射、吸收、生物电磁学、天线和 电磁兼容与电磁于扰等诸多领域中，取得 了令人满意的效果，其应用领域还在不断扩展，方法本 身也在不断完善和创新。 FDTD法直接用有限差分式代替麦克斯韦时域场旋度方程中的微分式，得到关于场分量的 有限差分式，用具有不同电参量的空间网格去模拟被研究体，选取合适的场初始值和计算空间 的边界条件，可以得到包括时间变量的麦克斯韦方程和四维数值解。通过傅立叶变换可求得三 维空间的频域解。这种方法的核心是将被研究媒质的只维空间沿坐标轴离散化，分成小六面体 单元网络，每个单元网络的电参量用媒质的参数(e,1,6,s?描述，用这种具有特定媒质参数的单 元网格的 堆‘砌护去模拟被研究体，并实施前述的电磁场的迭代运算 二、方法与结果 在传统的F7YID方法中，同一网点和相邻网点在上一时间步的电磁场值以及其它的边界条 件变量等是存储在计算机的内存中的，F)TD法的稳定条件又要求划分的网格尺寸要小于所研 究最小波长的十分之一，因此，被研究物体的电尺寸越大，计算所需占用的计算机内越大。举例 来说，一台具有32M内存的PC机仅能完成空间网格数大约为ROXROXRO目标的分析计算 所 以，对于分析具有电大尺寸目标的电磁问题，我们常采用专用Y十算工作站甚至巨型计算机来完 成，但是，如果我们仅有的硬件条件只是PC机和小内存工作站，或者经费不足以支持巨型计算 机的昂贵计算费用，而又必须完成电大尺寸目标的分析，找们就必须采用新的或者说改进的 FD17)方法 。 在这篇文章中，作者提出并实施 了一种可以采用现有手段解决这一问题的FDTF算法。这 421 电磁兼容论文集 种万法的核心思想是利用PC机的硬盘空间模拟内存空间，通过共用数组在应用程序级别实现 内外存交换。在这种方法中，整个计算空间沿某一坐标轴向分割成许多小段，如图I所示，每一 小段按顺序分别加载内存，实现计算机内存的复用。这一过程的细节描述如下: 图I是一个复杂形体的三维FD77〕计算空间，沿Y轴向被分割成许多小段，分别赋给标号 1,2,3n，其中侮 一小段的空间大小小于现有微机内存容量的允许范围，因此每一小段可以 装载内存完成计算。 只 ! 段 2 图1计算空间分割情况 在每一个确定的时间步中，首光，把第 一段中各网格上磁场量的初值装载于一个微机内存 中的数组，然后执行 MTDN叠代计算，完成一次扫描后，把这一段中各网格L的电磁场量的新 值存储在硬盘上;随后，相同的数组又被用来装载第三段中各网格L的电磁场量的初值，接着执 行相同的FDTD叠代ii算. 代结果又存储在硬盘 匕同样的过程一直持续到最后一段，所有的 各网格上电磁场量的叠代结果都被存储在硬盘 匕但是所需完成叠代的计算机内存却很小，完 全是P(:机所能承担的。在下一个时间步中 相邻上一个时间步存储在硬盘上的电磁场量的值 又被当作初值逐段地从硬盘上读到内存当中，与上一个时间步完全相同的计算过程又被重复执 行，一it进行到预定时fIH