伊犁师范学院硕士研究生————期末考核科目：电磁波有限时域差分方法姓名：***学号：*************学院：电子与信息工程学院专业：无线电物理时域有限差分法1 选题背景在多种可用的数值方法中，时域有限差分法(FDTD)是一种新近发展起来的可选方法。

1966年，K.S.Yee 首次提出电磁场数值计算的新方法—时域有限差分法(Finite Difference- Time Domain ，简称FDTD)。

经历了二十年的发展FDTD 法才逐渐走向成熟。

上世纪80年代后期以来FDTD 法进入了一个新的发展阶段，即由成熟转为被广泛接受和应用的阶段。

FDTD 法是解决复杂问题的有效方法之一，是一种直接基于时域电磁场微分方程的数值算法，它直接在时域将Maxwell 旋度方程用二阶精度的中心差分近似，从而将时域微分方程的求解转换为差分方程的迭代求解。

是电磁场和电磁波运动规律和运动过程的计算机模拟。

原则上可以求解任意形式的电磁场和电磁波的技术和工程问题，并且对计算机内存容量要求较低、计算速度较快、尤其适用于并行算法。

现在FDTD 法己被广泛应用于天线的分析与设计、目标电磁散射、电磁兼容、微波电路和光路时域分析、生物电磁剂量学、瞬态电磁场研究等多个领域[1]。

2 原理分析2.1 FDTD 的Yee 元胞E,H 场分量取样节点在空间和时间上采取交替排布，利用电生磁，磁生电的原理t t ∂∂=∂∂=⨯∇E D H ε t t ∂∂-=∂∂-=⨯∇HB E μ图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 的差分格式麦克斯韦第一、二方程 ⎪⎪⎩⎪⎪⎨⎧⎪⎭⎫⎝⎛+∂∂-=⨯∇+∂∂=⨯∇m t t J B E J D H （1）式中，J 时电流密度，反映电损耗，mJ 是磁流密度，单位2m V /，反映磁损耗。

主要与上式对应。

各向同性介质中的本构关系：H JE J H B E D m mγγμε==== （2）其中m γ是磁阻率，计算磁损耗的。

以H E ,为变量，在直角坐标中，展开麦克斯韦第一、二方程，分别为⎪⎪⎪⎩⎪⎪⎪⎨⎧+∂∂=∂∂-∂∂+∂∂=∂∂-∂∂+∂∂=∂∂-∂∂zz x y y y zx x xy z Et E y H x H E t E x H z H E t E z H yH γεγεγε （3） ⎪⎪⎪⎩⎪⎪⎪⎨⎧-∂∂-=∂∂-∂∂-∂∂-=∂∂-∂∂-∂∂-=∂∂-∂∂zm z x y y m y zx x m xy z Ht H y E x E H t H x E z E H t H z E yE γμγμγμ （4） 令()t ,z ,y ,x f 代表H E,在直角坐标中的任何一个分量，离散符号取为()()()k ,j ,i ft n ,z k ,y j ,x i f t z y x f n=∆∆∆∆=,,, （5）()t ,z ,y ,x f 关于时间和空间的一阶偏导数取中心差分近似为()()[]()()[]()()[]()()⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡-∆≈∂∂--+∆≈∂∂--+∆≈∂∂--+∆≈∂∂+∆=∆=∆=∆=k ,j ,i f k ,j ,i f t 1t f k ,j ,i f k ,j ,i f z 1z f k ,j ,i f k ,j ,i f y 1y f k ,j ,i f k ,j ,i f x 1x f 21-n 21n t n t 21n 21n z k z 21n21n yj y 21n21n xi x （6） 可以看出，每一节点上沿某一方向场分量的一阶偏微分可以用在该方向上相邻两点的一阶中心差商来描述，将式(1)用一阶中心差商方程取代，整理后便得到一阶差分方程，它具有二阶精度[3]。

Yee 元胞如图1所示，规定为1）剖分节点与场分量所在棱边中点不同，场分量的位置，即H E,节点是Yee 元胞节点的相对位置，不需要单独编码；2）当空间存在媒质分界面时，场量自动满足场的连续性条件，2t1t 2t 1t H H ,E E ==电磁分量的取样方式不仅符合法拉第电磁感应定律和安培环路定律的自然结构，也符合麦克斯韦方程的差分计算。

其次，时间步长可以取为电磁波传播一个空间步长所需时间的一半，因此E 与H 在时间顺序上交替抽样，时间间隔相差半个时间步长。

2.3 一维问题均匀平面波（TEM 波）是一维问题，设电磁波沿z 轴方向传播，则00==z z , H E ，场量和介质参数均与x ，y 无关，即0y ,0x =∂∂=∂∂，麦克斯韦方程为ym y xxxyH γt H μz E γE t E εz H +∂∂=∂∂-+∂∂=∂∂-（7）和xm x y yy x H γtH μz E γE tE εz H +∂∂=∂∂-+∂∂=∂∂- （8）旋转坐标轴后可以只保留一组公式[4]，设保留（7） Yee 元胞如图2所示E xH yz图2 一维Yee 元胞差分格式为()()()()()() k H k H z 1m CB k E m CA k E2121n y2121n y nx1n x⎪⎪⎭⎫ ⎝⎛--+∆-=+++ （9）()()()()()()[]k E 1k E z1m CQ k H m CP k H n x n x2121n y2121n y-+∆-+=+-+（10）如果介质无损耗，则0 ,0m ==γγ2.4 二维问题三维通常是散射问题，二维是TE 、TM 波问题，一维是TEM 波问题。

在二维场中，所有物理量与Z 坐标无关，既0z /=∂∂。

于是在TE 和TM 波的表达式分别为TE 波（0E z =） ⎪⎪⎪⎩⎪⎪⎪⎨⎧-∂∂-=∂∂-∂∂+∂∂=∂∂-+∂∂=∂∂H t H y E x E E t E x H E tE y H zm z x y y y zx x z γμγεγε （11）TM 波（0H z =） ⎪⎪⎪⎩⎪⎪⎪⎨⎧+∂∂=∂∂-∂∂-∂∂-=∂∂--∂∂-=∂∂zz x y y m y zxm x zEt E y H x H H t H x E H t H y E γεγμγμ （12）图3分别给出了TM 波和TE 波的Yee 元胞图图3 TM 波的Yee 元胞 图4 TE 波的Yee 元胞对于TE 波，只要令0=z E ，在z ∆上，yx H H , 不随z 变化，m 中去掉k 即可得到：()()()()()()j ,i H j ,i H y 1m CB j ,i E m CA j ,i E212121n z212121n z 21n x211n x⎪⎪⎭⎫ ⎝⎛-+-++∆++=++++ 式中：j ,i m 21+= （13）()()()()()()⎪⎪⎭⎫ ⎝⎛+--++∆-+=++++212121n z212121n z 21n y211n yj ,i H j ,i H x 1m CB j ,i E m CA j ,i E式中：21j ,i m += （14）()()()()()()()()⎥⎥⎦⎤⎢⎢⎣⎡∆+-++-∆+-++⋅-++=++-+yj ,i E 1j ,i E x j ,i E j ,1i E m CQj ,i H m CP j ,i H 21n x21n x21n y21n y212121n z212121n z式中，2121j ,i m ++= （15） 对TM 波，只要令0=z H ，在z ∆上，yx E E , 不随z 变化，m 中去掉k ，即可得到：()()()()()()[] j ,i E 1j ,i E y1m CQ j ,i H m CP j ,i H n z n z2121n x2121n x-+∆-+=+-+式中，21j ,i m += （16） ()()()()()()[]j ,i E j ,1i E x1m CQ j ,i H m CP j ,i H n z nz 2121n y2121n y-+∆++=+-+式中，j ,i m 21+= （17）()()()()()()()()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆--+-∆--++=+++++yj ,i H j ,i H x j ,i H j ,i H m CBj ,i E m CA j ,i E 2121n x 2121n x 2121n y 2121n y n z 1n z式中：j ,i m = （18）为了编写统一的TE 和TM 波二维FDTD 程序，可将描述TE 波差分公式(13)～(15)中相应的标号整体移动1/2，即坐标（x,y ）分别沿x 和y 轴方向移动半个网格，并将离散时间也移动半个时间步长，式(13)～(15)可以重新写为()()()()()()[] j ,i H 1j ,i Hy 1m CB j ,i E m CA j ,i E n z n z2121n x2121n x-+∆++=+-+式中：21j ,i m += （19）()()()()()()[]j ,i H j ,1i H x1m CB j ,i E m CA j ,i E n z n z 2121n y2121n y-+∆-+=+-+式中：j ,i m 21+= （20）()()()()()()()()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∆--+-∆--+-=+++++yj ,i E j ,i E x j ,i E j ,i E m CQj ,i H m CP j ,i H 2121n x 2121n x 2121n y 2121n y n z 1n z式中，j ,i m += （21） 可以看出，TE 波的FDTD 公式(19)～(21)与TM 波的FDTD 公式(16)～(18)形式相同，给编程带来极大方便。