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

时域有限差分法 (FDTD 算法 ) 的基本原理及仿真 时域有限差分法（ FDTD 算法） 时域有限差分法是 1966 年 K.S.Yee 发表在 AP 上的一篇论文建立起来的，后被称为 Yee网格空间离散方式。这种方法通过将 Maxwell 旋度方程转化为有限差分式而直接在时域求解 , 通过建立时间离散的递进序列 , 在相互交织的网格空间中交替计算电场和磁场。 FDTD 算法的基本思想是把带时间变量的 Maxwell 旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应。需要考虑的三点是差分格式、解的稳定性、吸收边界条件。有限差分通常采用的步骤是：采用一定的网格划 分方式离散化场域；对场内的偏微分方程及各种边界条件进行差分离散化处理，建立差分格式，得到差分方程组；结合选定的代数方程组的解法，编制程序，求边值问题的数值解。 FDTD 的基本原理 FDTD 方法由 Maxwell 旋度方程的微分形式出发，利用二阶精度的中心差分近似，直接将微分运算转换为差分运算，这样达到了在一定体积内和一段时间上对连续电磁场数据的抽样压缩。 Maxwell 方程的旋度方程组为： H E E E H mH t t 在直角坐标系中，（1）式可化为如下六个标量方程： H z H y E x E x Ez E y H x y z t y z t H x H z E y E y Ex Ez H y z x t ， x t z H y H x E z E z E y E x H z x y t x y t  （ 1） m H x m H y （2） m H z 上面的六个偏微分方程是 FDTD 算法的基础。 Yee 首先在空间上建立矩形差分网格，在时刻 n t 时刻， F(x,y,z) 可以写成 F (x, y, z,t ) F (i x, j y, k z, n t ) F n (i , j ,k) （3） 用中心差分取二阶精度： 对空间离散： F ( x, y, z,t) F n (i 1 2, j ,k) F n (i 1 2, j, k) 2 x i x x O x x F ( x, y, z,t ) F n (i, j 1 2,k) F n (i, j 1 2, k) 2 y j y y O y y F ( x, y, z, t ) F n (i, j , k 1 2) F n (i, j ,k 1 2) 2 z k z z O z z 对时间离散： F ( x, y, z, t) F n 1 2 (i , j , k) F n 1 2 (i, j ,k ) t 2 t n t t O （4） t Yee 把空间任一网格上的 E 和 H 的六个分量，如下图放置 : z Ey Ex Hz Ex Ey Ez Ez Ez Hy Hx Ey Ex o y x 图 1 Yee氏网格及其电磁场分量分布 在 FDTD 中，空间上连续分布的电磁场物理量离散的空间排布如图所示。由图可见，电场和磁场分量在空间交叉放置，各分量的空间相对位置也适合于 Maxwell 方程的差分计算， 能够恰当地描述电磁场的传播特性。 同时，电场和磁 场在时间上交替抽样，抽样时间间隔相差半个时间步，使 Maxwell 旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求 逆运算。因此，由给定相应电磁问题的初始条件， FDTD 就可以逐步推进地求得 以后各个时刻空间电磁场的分布。根据这一原则可以写出六个差分方程 : Exn 1(i (i  1, j,k) t 1/ 2, j ,k)  1 (i 1/ 2, j,k) t 2 (i 1/ 2, j,k) n 1 (i 1/ 2, j, k) t .Ex (i 1 / 2, j ,k) 2 (i 1/ 2, j,k) . 1 t . （5） (i 1/ 2, j, k) 1 2 (i 1/ 2, j,k) n 1/2 n 1/ 2 n 1/2 n 1/ 2 (i 1/ 2, j,k 1/ 2) [ Hz (i 1/ 2, j,k) Hz (i 1/ 2, j 1/ 2, k) H y (i 1/ 2, j, k 1/ 2) Hy z ] y 其余的也如法可以写出，每个网格点上的个场分两的新值依赖于该点在前一时间步长时刻的值机该点周围的临近点上另一场量在早半个时间步长时的值。因此任一时刻可一次算出一个点，并行算法可计算出多个点。通过这些运算可以交替算出电场磁场在各个时间步的值。 根据上述 FDTD 差分方程组可得出计算电磁场的时域推进计算方法， 如图 2 所示。 已知 t1 t0 n t 0 时刻空间各处的电磁场初始 计算 t 2 t1 t / 2