时域有限差分法(姚伟)介绍.doc

伊犁师范学院硕士研究生————期末考核科目：电磁波有限时域差分方法姓名：姚伟学号：学院：电子与信息工程学院专业：无线电物理时域有限差分法1 选题背景在多种可用的数值方法中，时域有限差分法 (FDTD)是一种新近发展起来的可选方法 1966 年，首次提 出电 磁场 数值计算 的新 方法 — 时 域有 限差分法 (Finite Difference- Time Domain ，简称 FDTD)经历了二十年的发展 FDTD法才逐渐走向成熟上世纪 80 年代后期以来 FDTD法进入了一个新的发展阶段，即由成熟转为被广泛接受和应用的阶段 FDTD法是解决复杂问题的有效方法之一， 是一种直接基于时域电磁场微分方程的数值算法，它直接在时域将 Maxwell 旋度方程用二阶精度的中心差分近似，从而将时域微分方程的求解转换为差分方程的迭代求解是电磁场和电磁波运动规律和运动过程的计算机模拟原则上可以求解任意形式的电磁场和电磁波的技术和工程问题，并且对计算机内存容量要求较低、计算速度较快、尤其适用于并行算法现在 FDTD法己被广泛应用于天线的分析与设计、目标电磁散射、电磁兼容、微波电路和光路时域分析、生物电磁剂量学、瞬态电磁场研究等多个领域 [1] 。

2 原理分析2.1 FDTD 的 Yee 元胞E,H 场分量取样节点在空间和时间上采取交替排布，利用电生磁，磁生电的原理DEBHHtEttt图1 Yee模型如图 1 所示， Yee单元有以下特点 [2] ：1）E 与 H 分量在空间交叉放置，相互垂直；每一坐标平面上的E 分量四周由H 分量环绕， H分量的四周由 E 分量环绕；场分量均与坐标轴方向一致2）每一个 Yee 元胞有 8 个节点， 12 条棱边， 6 个面棱边上电场分量近似相等，用棱边的中心节点表示，平面上的磁场分量近似相等，用面的中心节点表示3）每一场分量自身相距一个空间步长， E 和 H相距半个空间步长4）每一场分量自身相距一个时间步长， E 和 H相距半个时间步长， 电场取 n 时刻的值，磁场取 n+0.5 时刻的值；即：电场 n 时刻的值由 n-1 时刻的值得到，磁场 n+0.5 时刻的值由 n-0.5 时刻的值得到；电场 n 时刻的旋度对应 n+0.5 时刻的磁场值，磁场 n+0.5 时刻的旋度对应 (n+0.5)+0.5 时刻的电场值，逐步外推5） 3 个空间方向上的时间步长相等，以保证均匀介质中场量的空间变量与时间变量完全对称。

应用这种离散方式，将含时间变量的 Maxwell 方程转化为一组差分方程，并在时间轴上逐步推进地求解空间电磁场由电磁问题的初值和边界条件，就可以逐步推进地求解以后各时刻空间电磁场分布2.2 Maxwell 方程 FDTD的差分格式HDJ麦克斯韦第一、二方程t（1）BEJ mt式中， J 时电流密度，反映电损耗，J m 是磁流密度，单位 V/ m2，反映磁损耗主要与上式对应各向同性介质中的本构关系：D E B H J EJmm H（ 2）其中 m 是磁阻率，计算磁损耗的以 E , H 为变量，在直角坐标中，展开麦克斯韦第一、二方程，分别为H zH yExExyztH xH zEyE y（ 3）zxtH yH xEzEzxytEzEyH xm H xyztE xEzH ym H y（4）zxtE yExH zm H zxyt令 f x, y, z,t 代表 E, H 在直角坐标中的任何一个分量，离散符号取为f x , y , z,t f ix, j y,k z,n t f n i , j ,k（ 5）f x, y, z,t 关于时间和空间的一阶偏导数取中心差分近似为fxfyx i xy j y1 fx1fynni21 , j ,kfi , j21 ,kfnni21 , j ,ki , j21 ,k（ 6）fzft1fni , j ,k1fni , j ,kzkzz2111fni , j ,kfn-i , j ,k22tntt12可以看出，每一节点上沿某一方向场分量的一阶偏微分可以用在该方向上相邻两点的一阶中心差商来描述，将式 (1) 用一阶中心差商方程取代，整理后便得到一阶差分方程，它具有二阶精度 [3] 。

Yee元胞如图 1 所示，规定为1）剖分节点与场分量所在棱边中点不同，场分量的位置，即E, H 节点是 Yee 元胞节点的相对位置，不需要单独编码；2）当空间存在媒质分界面时，场量自动满足场的连续性条件，E1tE2t , H 1t H 2t电磁分量的取样方式不仅符合法拉第电磁感应定律和安培环路定律的自然结构，也符合麦克斯韦方程的差分计算其次，时间步长可以取为电磁波传播一个空间步长所需时间的一半，因此 E 与 H 在时间顺序上交替抽样，时间间隔相差半个时间步长2.3 一维问题均匀平面波（ TEM波）是一维问题，设电磁波沿z 轴方向传播，则 Ez0, H z 0 ，0,0场量和介质参数均与 x，y 无关，即xy，麦克斯韦方程为H yεE xγEztx（7）E xμH yγm H yzt和H xεEyγEzty（8）E yH xμγmH xzt旋转坐标轴后可以只保留一组公式[4]，设保留（ 7）Yee元胞如图 2 所示E xH yz图 2 一维 Yee 元胞差分格式为n1n1n11n11k CB mH y2 kH y2 kE xk CA m Exz22（9）n1n112121nnH yCP m H yCQ mkk 2k 2Exk 1 E xz（ 10）如果介质无损耗，则0,m02.4 二维问题三维通常是散射问题，二维是TE、 TM波问题，一维是 TEM波问题。

在二维场中，所有物理量与Z 坐标无关，既/ z0 于是在 TE和 TM波的表达式分别为H zExExytTE 波（ E z0）H zEyEy（ 11）xtEyExH zm H zxytEzH xm H xytTM波（ H z0）EzH ym H y（ 12）xtH yH xE zEz。