第10讲 时域有限差分法 I.ppt

电磁波时域有限差分方法 Finite-difference time-domain method for electromagnetic wave 目录 引言 Maxwell方程及其FDTD形式 数值稳定性 吸收边界 激励源 近－远场外推 应用算例 其它参考书 第一章参考文献 [1] Yee(1966) 第一篇FDTD论文 [36] Kunz(1993) [44] Taflove(1995, 2000 Second Ed.) [46] Sullivan(2000) [42] 王长清（1994） [43] 高本庆（1995） 电磁场计算方法：解析方法 典型目标 散射问题：球，圆柱，劈 辐射问题：平面波，线电流，电偶极子 用途 实际问题的近似 导出新的物理概念和方法（如GTD） 其它计算方法的验证算例 电磁场计算方法：解析方法 用途 广泛应用于电磁散射和雷达截面的计算、波导与谐振腔系统、辐射天线分析、周期结构分析、电子封装和电磁兼容分析、核电磁脉冲的传播和散射、微光学元器件中光的传播和衍射特性、电磁波生物效应、微波及毫米波集成电路分析、超高速集成电路互连封装电磁特性的分析、双负介质以及各向异性介质中的电磁波传播、逆散射与遥感、地下电磁探测和电磁成像等方面。 电磁场计算方法：高频技术 目标电磁散射(如RCS)特性的理论建模与分析是雷达共性基础研究中不可或缺的重要课题。在雷达常用的波段上，大多数的目标既表现出宏观的电大尺寸，同时又不可避免地具有细节上的复杂结构。 长期以来，人们普遍应用解析解方法和各种高频近似方法(如几何光学(GO)法，物理光学(PO)法，几何绕射理论(GTD)，物理绕射理论(PTD)，复射线(CR)法或弹射射线(SBR)法)来分析各类电大尺寸目标的电磁散射特性。 这些方法的主要优点是简单明晰，容易掌握，计算方便，甚至可“实时”显示近似计算结果。但其普遍的缺点则是理论模型粗糙，计算精度太低。 电磁场计算方法：高频技术 几何光学（GO）方法，射线方法，弹跳射线方法——频域 几何绕射理论（GTD），一致性几何绕射方法(UTD)，物理绕射理论(PTD) ——频域 物理光学（PO）方法——频域，时域 等效边缘电磁流（EEC）方法——频域，时域 电磁场计算方法：数值方法 矩量法（MoM），快速多极子方法FMM） ——频域，时域 有限元方法（FEM） ——频域，时域 边界元方法（BEM） ——频域 有限差分（FD）方法 时域有限差分（FDTD）方法——时域 第一章 引言 FDTD的发展 FDTD的应用 FDTD的基本点 1.1 FDTD的发展：提出 Yee（1966年）[1]首先提出Maxwell方程的差分离散方式，并用来处理电磁脉冲的传播和反射问题。 FDTD的发展(续)：推向应用 Taflove 等（1975年）[5]用FDTD计算非均匀介质在正弦波入射时的时谐场（稳态）电磁散射，讨论了时谐场情况的近－远场外推，以及数值稳定性条件。 Holland（1977年）[6]和Kunz（1978年）[7]用FDTD计算F117飞机这种复杂目标的电磁脉冲散射。 FDTD的发展(续)：时域外推 Britt （1989年）[21]首次给出时域远场结果，但论文未给出外推具体方法。 Yee 等（1991年）[22]和Luebbers等（1991年）[23]提出了三维FDTD时域近－远场外推方法，随后Luebbers等（1992年）[24]提出二维FDTD时域近－远场外推方法。 FDTD的发展(续)：吸收边界1 Mur（1981年）[8]提出在计算区域截断边界处的一阶和二阶吸收边界条件及其在FDTD的离散形式。这是FDTD的一种十分有效的吸收边界条件，获得广泛应用。 FDTD的发展(续) ：吸收边界2 Berenger （1994，1996年）[30-32]提出将麦克斯韦方程扩展为场分量分裂形式，并构成完全匹配层（PML)。 Sacks 等（1995年）[33]和Gedney （1996年）[34]提出各向异性介质的PML，其支配方程是各向异性介质麦克斯韦方程。在FDTD计算中这两种PML作为吸收边界已得到广泛应用。 1.2 FDTD的应用 天线辐射的分析 散射和雷达截面计算 微波器件和导行波结构的研究 周期结构分析 电子封装，电磁兼容分析 光学元器件中光的传播和衍射特性 同轴线馈电时无限大理想导体地面上轴对称天线 天线辐射场 金属方柱散射场 介质圆柱散射场 用于FDTD的“战斧”导弹模型 1.3 FDTD的基本点(1)：Yee元胞 最初由Yee提出（1966） E、H场分量节点在空间和时