时域有限差分法fdtd算法的基本原理及仿真.docVIP

时域有限差分法（FDTD算法） 时域有限差分法是1966年K.S.Yee发表在AP上的一篇论文建立起来的，后被称为Yee网格空间离散方式。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解, 通过建立时间离散的递进序列, 在相互交织的网格空间中交替计算电场和磁场。 FDTD算法的基本思想是把带时间变量的Maxwell旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应。需要考虑的三点是差分格式、解的稳定性、吸收边界条件。有限差分通常采用的步骤是：采用一定的网格划分方式离散化场域；对场内的偏微分方程及各种边界条件进行差分离散化处理，建立差分格式，得到差分方程组；结合选定的代数方程组的解法，编制程序，求边值问题的数值解。 FDTD的基本原理 FDTD方法由Maxwell旋度方程的微分形式出发，利用二阶精度的中心差分近似，直接将微分运算转换为差分运算，这样达到了在一定体积内和一段时间上对连续电磁场数据的抽样压缩。 Maxwell方程的旋度方程组为： （1） 在直角坐标系中，（1）式可化为如下六个标量方程： ， （2） 上面的六个偏微分方程是FDTD算法的基础。 Yee首先在空间上建立矩形差分网格，在时刻时刻，F(x,y,z)可以写成 （3） 用中心差分取二阶精度： 对空间离散： 对时间离散： （4） Yee把空间任一网格上的E和H的六个分量，如下图放置: 图1 Yee氏网格及其电磁场分量分布 在FDTD中，空间上连续分布的电磁场物理量离散的空间排布如图所示。由图可见，电场和磁场分量在空间交叉放置，各分量的空间相对位置也适合于Maxwell方程的差分计算，能够恰当地描述电磁场的传播特性。同时，电场和磁场在时间上交替抽样，抽样时间间隔相差半个时间步，使Maxwell旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求逆运算。因此，由给定相应电磁问题的初始条件，FDTD就可以逐步推进地求得以后各个时刻空间电磁场的分布。根据这一原则可以写出六个差分方程: （5 （5） 其余的也如法可以写出，每个网格点上的个场分两的新值依赖于该点在前一时间步长时刻的值机该点周围的临近点上另一场量在早半个时间步长时的值。因此任一时刻可一次算出一个点，并行算法可计算出多个点。通过这些运算可以交替算出电场磁场在各个时间步的值。 已知 0时刻空间各处的电磁场初始值根据上述FDTD 已知 0时刻空间各处的电磁场初始值 计算 时刻空间各处的磁场值 计算 时刻空间各处的磁场值 计算 时刻空间各处的电场值 循环n 计算 时刻空间各处的电场值 图2 FDTD在时域的交叉半步逐步推进计算 2.数值稳定性条件 时间步长，空间步长，，必须满足一定的关系，否则就使得数值表现不稳定，表现为：随着计算步数的增加，计算场量的数值会无限的增大，这种增大不是由于误差积累造成的，而是由于电磁波的传播关系被破坏造成的。所以，，，必须满足一定的关系以保证稳定性。Taflove等在1975年对Yee氏差分格式的稳定性进行了讨论，并导出了对时间步长的限制条件。数值解是否稳定主要取决于时间步长与空间步长、、的关系。对于非均匀媒质构成的计算空间选用如下的稳定性条件： （6） 若采用均匀立方体网格：， （7） 而一般取：，c为光速。 当，，不相等时， （8） 3.数值色散 FDTD网格中，会导致数字波模在网格中发生改变，这种改变是由于计算网格本身引起的，而非物理因素，所以必须考虑。即在FDTD网格中，电磁波的相速与频率有关，电磁波的相速度随波长、传播方向及变量离散化的情况不同而改变。色散将导致非物理因素引起的脉冲波形畸变、人为的各向异性和虚假折射等现象。显然，色散与空间、时间的离散间隔有关，如下式所示： （9） 与数值色散关系相对应，在无耗介质中的单色平面波，色散解析关系是： （10） 由式（9）可知，当式（9）中的、、、均趋于零时，它就趋于式（10）。也就是说数值色散是由于用