时域有限差分法(姚伟)介绍 下载本文

伊犁师范学院硕士研究生

科目:电磁波有限时域差分方法姓名:姚伟学号:学院:电子与信息工程学院专业:无线电物理

————期末考核

1076411203009

时域有限差分法

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场分量取样节点在空间和时间上采取交替排布,利用电生磁,磁生电的原理

?D?E?B?H??H?????E??????t?t ?t?t

图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的差分格式

?D???H??J ??t?麦克斯韦第一、二方程 ? (1)

?B????E????Jm?????t??式中,J时电流密度,反映电损耗,Jm是磁流密度,单位V/m,反映磁损耗。主要

2与上式对应。各向同性介质中的本构关系:

D??E B??H J??E Jm??mH (2) 其中?m是磁阻率,计算磁损耗的。

以E ,H为变量,在直角坐标中,展开麦克斯韦第一、二方程,分别为

??Hz?Hy?Ex?????Ex??y?z?t??Ey??Hx?Hz?????Ey (3) ??z?x?t???Hy??Hx???Ez??Ez??x?y?t???Ez?Ey?Hx??????mHx??z?t??y?Hy??Ex?Ez ???????mHy (4)

?z?x?t???Ey??Ex????Hz??Hmz??x?y?t?令f?x,y,z,t?代表E,H在直角坐标中的任何一个分量,离散符号取为

f?x,y,z,t??f?i?x,j?y,k?z,n?t??fn?i,j,k? (5)

f?x,y,z,t?关于时间和空间的一阶偏导数取中心差分近似为

??f1n11n?fi?,j,k?fi?,j,k??x22?x?x?i?x??f1n11n?fi,j?,k?fi,j?,k??y22?y?y?j?y (6) ?1n11??f?fi,j,k?2?fni,j,k?2??zz?k?z?z?11n-?n?2??f12??????fi,j,k?fi,j,k?? ????tt?n?t?t???????????????????可以看出,每一节点上沿某一方向场分量的一阶偏微分可以用在该方向上相邻两点的

一阶中心差商来描述,将式(1)用一阶中心差商方程取代,整理后便得到一阶差分方程,

[3]

它具有二阶精度。

Yee元胞如图1所示,规定为

1)剖分节点与场分量所在棱边中点不同,场分量的位置,即E,H节点是Yee元胞节点的相对位置,不需要单独编码;

2)当空间存在媒质分界面时,场量自动满足场的连续性条件,E1t?E2t, H1t?H2t 电磁分量的取样方式不仅符合法拉第电磁感应定律和安培环路定律的自然结构,也符合麦克斯韦方程的差分计算。其次,时间步长可以取为电磁波传播一个空间步长所需时间的一半,因此E与H在时间顺序上交替抽样,时间间隔相差半个时间步长。 2.3 一维问题

均匀平面波(TEM波)是一维问题,设电磁波沿z轴方向传播,则Ez?0, Hz?0,

???0, ?0?y场量和介质参数均与x,y无关,即?x,麦克斯韦方程为

?Ex?γEx?z?t (7)

?Hy?Ex??μ?γmHy?z?t??ε?Hy和

?Ey?Hx??ε?γEy?z?t (8)

?Ey?Hx??μ?γmHx?z?t旋转坐标轴后可以只保留一组公式[4],设保留(7) Yee元胞如图2所示

Ex

Hy 图2 一维Yee元胞

z

差分格式为

En?1x?n?2?k??CA?m?E?k??CB?m?1?Hyk???z?nx1?12? ??H?k????n?y1212? (9)

Hyn?121?E?k???CP?m?H?k?? ?CQ?m??z12n?y1212nx?k?1??Exn?k?? (10)

如果介质无损耗,则??0, ?m?0 2.4 二维问题

三维通常是散射问题,二维是TE、TM波问题,一维是TEM波问题。

在二维场中,所有物理量与Z坐标无关,既?/?z?0。于是在TE和TM波的表达式分别为

?Ex??Hz????Ex ??y?t??Ey??Hz????Ey (11) TE波(Ez?0) ? ??x?t?? ?Ey??Ex????Hz??H mz??x?y?t??Hx??Ez?????mHx??y?t??Hy??Ez?????mHy (12) TM波(Hz?0) ???x?t?? ?Hy??Hx???Ez??Ez??x?y?t?图3分别给出了TM波和TE波的Yee元胞图

(i, j+1/2)HxHyEyHz(i, j+1/2)(i+1/2, j+1/2)Ex(i, j)(i+1/2, j)Ez

图3 TM波的Yee元胞 图4 TE波的Yee元胞