时域有限差分法(FDTD算法)的基本原理及仿真.pdf

时域有限差分法〔FDTD算法〕时域有限差分法是1966年发表在AP上的一篇论文建立起来的，后被称为Yee网格空间离散方式。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解,通过建立时间离散的递进序列,在相互交织的网格空间中交替计算电场和磁场。FDTD算法的基本思想是把带时间变量的Maxwell旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应。需要考虑的三点是差分格式、解的稳定性、吸收边界条件。有限差分通常采用的步骤是：采用一定的网格划分方式离散化场域；对场内的偏微分方程及各种边界条件进行差别离散化处理，建立差分格式，得到差分方程组；结合选定的代数方程组的解法，编制程序，求边值问题的数值解。1.FDTD的基本原理FDTD方法由Maxwell旋度方程的微分形式出发，利用二阶精度的中心差分近似，直接将微分运算转换为差分运算，这样到达了在一定体积内和一段时间上对连续电磁场数据的抽样压缩。Maxwell方程的旋度方程组为：EHHEEH〔1〕ttm在直角坐标系中，〔1〕式可化为如下六个标量方程：HHEEzyxEEHzyxyztxHyztmxHHEEEHxzyE，xzyH〔2〕zxtymyzxtHHEEEHyxzEyxzHxytzmzxyt上面的六个偏微分方程是FDTD算法的基础。Yee首先在空间上建立矩形差分网格，在时刻nt时刻，F(x,y,z)可以写成F(x,y,z,t)F(ix,jy,kz,nt)Fn(i,j,k)〔3〕用中心差分取二阶精度：对空间离散：F(x,y,z,t)Fn(i12,j,k)Fn(i12,j,k)Ox2xxixxF(x,y,z,t)Fn(i,j12,k)Fn(i,j12,k)Oy2yyjyyF(x,y,z,t)Fn(i,j,k12)Fn(i,j,k12)Oz2zzkzz对时间离散：F(x,y,z,t)Fn12(i,j,k)Fn12(i,j,k)Ot2〔4〕ttnttYee把空间任一网格上的E和H的六个分量，如下列图放置:zEyExHzExEyEzEzEzHyHxEyExoyx图1Yee氏网格及其电磁场分量分布在FDTD中，空间上连续分布的电磁场物理量离散的空间排布如下图。由图可见，电场和磁场分量在空间交叉放置，各分量的空间相对位置也适合于Maxwell方程的差分计算，能够恰当地描述电磁场的传播特性。同时，电场和磁场在时间上交替抽样，抽样时间间隔相差半个时间步，使Maxwell旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求逆运算。因此，由给定相应电磁问题的初始条件，FDTD就可以逐步推进地求得以后各个时刻空间电磁场的分布。根据这一原则可以写出六个差分方程:(i1/2,j,k)t12(i1/2,j,k)En1(i1,j,k).En(i1/2,j,k)x(i1/2,j,k)tx12(i1/2,j,k)t1〔5〕..(i1/2,j,k)t(i1/2,j,k)12(i1/2,j,k)Hn1/2(i1/2,j,k)Hn1/2(i1/2,j1/2,k)Hn1/2(i1/2,j,k1/2)Hn1/2(i1/2,j,k1/2)[zzyy]yz其余的也如法可以写出，每个网格点上的个场分两的新值依赖于该点在前一时间步长时刻的值机该点周围的临近点上另一场量在早半个时间步长时的值。因此任一时刻可一次算出一个点，并行算法可计算出多个点。通过这些运算可以交替算出电场磁场在各个时间步的值。根据上述FDTD差分方程组可得出计算电磁场的时域推进计算方法，如图2ttnt已知100时刻空间各处的电磁场初始所示。计算ttt/2时刻空间各处的磁场值21循环n次ttt/2计算12时刻空间各处的电场值图2FDTD在时域的交叉半步逐步推进计算时间步长t，空间步长x，y，z必须满足一定的关系，否则就使得数值表现不稳定，表现为：随着计算步数的增加，计算场量的数值会无限的增大，这种增大不是由于误差积累造成的，而是由于电磁波的传播关系被破坏造成的。所以t，x，y，z必须满足一定的关系以保证稳定性。Taflove等在1975年对Yee氏差分格式的稳定性进行了讨论，并导出了对时间步长的限制条件。数值解是否稳定主要取决于时间步长t与空间步长x、y