第10讲时域有限差分法I

,Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,*,电磁波时域有限差分方法,Finite-difference time-domain method for electromagnetic wave,1,目录,引言,Maxwell,方程及其FDTD形式,数值稳定性,吸收边界,激励源,近远场外推,应用算例,2,其它参考书,第一章参考文献,1 Yee(1966),第一篇FDTD论文,36 Kunz(1993),44 Taflove(1995, & 2000 Second Ed.),46 Sullivan(2000),42 王长清（1994）,43 高本庆（1995）,3,电磁场计算方法：解析方法,典型目标,散射问题：球，圆柱，劈,辐射问题：平面波，线电流，电偶极子,用途,实际问题的近似,导出新的物理概念和方法（如GTD）,其它计算方法的验证算例,4,电磁场计算方法：解析方法,用途,广泛应用于电磁散射和雷达截面的计算、波导与谐振腔系统、辐射天线分析、周期结构分析、电子封装和电磁兼容分析、核电磁脉冲的传播和散射、微光学元器件中光的传播和衍射特性、电磁波生物效应、微波及毫米波集成电路分析、超高速集成电路互连封装电磁特性的分析、双负介质以及各向异性介质中的电磁波传播、逆散射与遥感、地下电磁探测和电磁成像等方面。,5,电磁场计算方法：高频技术,目标电磁散射(如RCS)特性的理论建模与分析是雷达共性基础研究中不可或缺的重要课题。在雷达常用的波段上，大多数的目标既表现出宏观的电大尺寸，同时又不可避免地具有细节上的复杂结构。,长期以来，人们普遍应用解析解方法和各种高频近似方法(如几何光学(GO)法，物理光学(PO)法，几何绕射理论(GTD)，物理绕射理论(PTD)，复射线(CR)法或弹射射线(SBR)法)来分析各类电大尺寸目标的电磁散射特性。,这些方法的主要优点是简单明晰，容易掌握，计算方便，甚至可“实时”显示近似计算结果。但其普遍的缺点则是理论模型粗糙，计算精度太低。,6,电磁场计算方法：高频技术,几何光学（GO）方法，射线方法，弹跳射线方法频域,几何绕射理论（GTD），一致性几何绕射方法(UTD)，物理绕射理论(PTD) 频域,物理光学（PO）方法频域，时域,等效边缘电磁流（EEC）方法频域，时域,7,电磁场计算方法：数值方法,矩量法（MoM），快速多极子方法FMM）,频域，时域,有限元方法（FEM） 频域，时域,边界元方法（BEM） 频域,有限差分（FD）方法,时域有限差分（FDTD）方法时域,8,第一章 引言,FDTD的发展,FDTD的应用,FDTD的基本点,9,1.1 FDTD,的发展：提出,Yee（1966年）1首先提出Maxwell方程的差分离散方式，并用来处理电磁脉冲的传播和反射问题。,10,FDTD,的发展(续)：推向应用,Taflove 等（1975年）5用FDTD计算非均匀介质在正弦波入射时的时谐场（稳态）电磁散射，讨论了时谐场情况的近远场外推，以及数值稳定性条件。,Holland（1977年）6和Kunz（1978年）7用FDTD计算F117飞机这种复杂目标的电磁脉冲散射。,11,FDTD,的发展(续)：时域外推,Britt （1989年）21首次给出时域远场结果，但论文未给出外推具体方法。,Yee 等（1991年）22和Luebbers等（1991年）23提出了三维FDTD时域近远场外推方法，随后Luebbers等（1992年）24提出二维FDTD时域近远场外推方法。,12,FDTD,的发展(续)：吸收边界1,Mur（1981年）8提出在计算区域截断边界处的一阶和二阶吸收边界条件及其在FDTD的离散形式。这是FDTD的一种十分有效的吸收边界条件，获得广泛应用。,13,FDTD,的发展(续) ：吸收边界2,Berenger （1994，1996年）30-32提出将麦克斯韦方程扩展为场分量分裂形式，并构成完全匹配层（PML)。,Sacks 等（1995年）33和Gedney （1996年）34提出各向异性介质的PML，其支配方程是各向异性介质麦克斯韦方程。在FDTD计算中这两种PML作为吸收边界已得到广泛应用。,14,1.2 FDTD,的应用,天线辐射的分析,散射和雷达截面计算,微波器件和导行波结构的研究,周期结构分析,电子封装，电磁兼容分析,光学元器件中光的传播和衍射特性,15,同轴线馈电时无限大理想导体地面上轴对称天线,16,天线辐射场,脉冲越过天线锥体,顶部棱角后的波面,沿顶面传播到另一棱,角后辐射的新波面,17,金属方柱散射场,幅值 相位,18,介质圆柱散射场,幅值 相位,19,用于FDTD的,“战斧”,导弹模型,复杂目标建模步骤,20,1.3 FDTD,的基本点(1)：,Yee元胞,最初由Yee提出（1966）,E、H场分量节点在空间和时间上采取交替排布，每一个E（或H）场分量周围有四个H（或E）场分量环绕,将Maxwell旋度方程转化为一组差分方程，并在时间轴上逐步推进地求解；由电磁问题的初始值及边界条件逐步推进地求得以后各时刻空间电磁场分布,21,Yee元胞,分量节点位置见p.10,表2-1,E取n时刻，H取n1/2时刻,22,电场z分量被H分量环绕,23,磁场z分量被E分量环绕,24,FDTD,的基本点(2)：,FDTD区的划分,对于散射问题，通常在FDTD计算区域中引入总场边界（即连接边界），如图1-2-1所示。FDTD计算区域划分为总场区和散射场区。这样做的好处是：一、应用惠更斯（Huygens）原理, 可以在连接边界处设置入射波，使入射波的加入变得简单易行；二、可以在截断边界（即吸收边界）处设置吸收边界条件，利用有限计算区域就能够模拟开域的电磁散射过程；三、根据等效原理，应用数据存储边界（即输出边界）处的近区场便可以实现远场的外推计算。对于辐射问题，激励源直接加到辐射天线上，整个FDTD计算区域为辐射场，如图1-2-2所示，不再区分总场区和散射场区。,25,FDTD,的基本点(2)：,FDTD区的划分,对于散射问题，划分为总场区和散射场区。,26,FDTD,的基本点(2)：,FDTD区的划分,对于辐射问题，激励源直接加到辐射天线上，整个FDTD计算区域为辐射场区,27,FDTD,的基本点(3)：,吸收边界条件,为了在有限计算区域模拟无界空间中的电磁问题，必须在计算区域的截断边界上设置吸收边界条件。,吸收边界从开始简单的插值边界，已经发展了多种吸收边界条件。目前比较广泛采用的有Mur吸收边界；以及近几年发展的完全匹配层(PML)吸收边界。,28,FDTD,的基本点(4)：,近远场变换,FDTD的模拟只能限于有限空间，为了获得计算域以外的散射或辐射场，必须借助等效原理应用计算区域内的近场数据实现计算区域以外远场的外推。,对于时谐场和瞬态场分别采用不同的外推方法。,29,近远场外推,30,第二章 Maxwell方程及其FDTD形式,麦克斯韦方程和Yee元胞,直角坐标中的FDTD：三维情形,直角坐标中的FDTD：二维情形,直角坐标中的FDTD：一维情形,柱坐标中的FDTD（7.6节）,球坐标中的FDTD（7.7节）,31,2.1 Maxwell方程和Yee元胞,麦克斯韦旋度方程,32,Maxwell,旋度方程直角分量式,电场,x,分量,磁场,x,分量,33,中心差分近似,在时间和空间域中的离散取以下符号：,中心差分近似,34,FDTD离散中的Yee元胞,分量节点位置见p.10,表2-1,E取n时刻，H取n1/2时刻,35,2.2 直角坐标中的FDTD：三维情形,电场时间推进,计算式,36,三维FDTD系数,37,三维FDTD公式（续）,磁场时间推进,计算式,38,三维FDTD系数,39,电场,z,分量被H分量环绕,40,磁场z分量被E分量环绕,41,FDTD在时域的交叉半步逐步推进计算流程,42,2.3 直角坐标中的FDTD：二维情形（TE波）,43,二维TE和TM波Yee元胞,(a)TM波 (b) TE波,44,二维FDTD时域推进公式,45,二维FDTD时域推进公式（续）,46,TE和TM波之间的对偶关系,可以编写统一适用于TE和TM波,情况的二维FDTD计算程序,47,2.4 直角坐标中的FDTD：一维情形（TEM波）,48,一维电磁场分量节点取样,49,一维FDTD离散公式,50,FDTD,方法特点,随时间推进计算，无需矩阵求逆,电磁场节点当前值只与前一时刻节点值相关，只需存储当前时间步电磁场值,电磁场节点只与相邻节点相关，便于并行计算,通过空间节点给介质参数赋值，便于处理复杂目标；推进计算程序可以用典型目标检验,可以分析开域（引入吸收边界）和闭域多种电磁问题,51,