电磁场与电磁波 PPT-32.ppt

电磁场与电磁波 时域有限差分方法（FD?TD） 时域有限差分（FD?TD）方法，不仅在三维空间对场量离散，而且在时间轴上也对场量离散，并用差分代替微分，得到麦克斯韦方程在一定边值与初值条件下的近似解。 差分法的优点是，解题有很大的灵活性，缺点是计算工作量大。 如果在x、y、z方向各取Nx、Ny、Nz个离散点，在（t1–t2）时间再取Nt个时间间隔，对于一个具体的场问题就要计算（Nx?Ny?Nz?Nt）个离散未知量，其计算工作量之大可想而知。 近代高速大容量数字计算机的发展，使时域有限差分法的应用成为可能，因而自20世纪80年代以来，时域有限差分法得到了极大的发展。在计算电磁学中得到广泛应用。 基本方程式（适合处理电各向异性介质） 对于单轴晶体 网格划分 E，D和H在空间位置上相差 半个差分步长，以后对时间 差分上E和H也相差半个步长。 这样保证了差分方程的形式 为中心差分，有二阶差分精 度，同时可以获得蛙跳式差 分格式。 定义： ?x，?y，?z 分别表示在空间x, y, z 三个方向上的差分步长，?t 表示时间上的差分间隔。Fn(i, j, k)表示在（i?x, j?y, k?z）位置上n?t时刻的场分量。 时间和空间上用差分代替微分 用差分代替微分 时间和空间上用差分代替微分 用差分代替微分 时域有限差分法基本算法 求得D分量后，E分量可以利用 得出 上面诸式给出了从n?t 时刻的场计算（n+1）?t 时刻的场，它们是FD-TD法的基本算式。 稳定性条件 对于任何一个差分方程组，为了使每一步差分求解所引起的误差不会积累，也就是使差分方程有稳定的解，差分步长?x, ?y, ?z和?t必须满足一定的关系，这种关系称为差分方程的稳定性条件。 从偏微分方程差分解的稳定性条件可以导出： 式中vp为所分析波导的最大相速。 边界条件—等效电壁、磁壁和介质交界面 电壁：对于切向电场或法向磁场，其值为零； 对于法向电场与切向磁场，法线方向导数为零 磁壁：对于切向磁场或法向电场，其值为零； 法向磁场与切向电场，法线方向导数为零 介质分界面： 相加取平均可得： 同样可得Dy，Dz，经组合便得到 由此，在求边界上的场量时，可引入一等效介电常数 在介质交界面上从D求E时，用 ‘ 代替 即可。 边界条件—开放边界 对于开放边界，严格来说只有向外传播的波，因而这类边界条件也可称之为辐射条件，或者称之为匹配边界条件。 原则上只要使开放的边界在无穷远处，就不会引起反射，但这在计算机模拟时是不可能实现的。 处理方法之一： 开放边界问题转化成为电壁、磁壁类边界问题。 在其他条件完全相同的情况下，一个边界等效开路（?0=1）与另一个边界等效短路（?s=–1）的场分布的叠加可以得到一个匹配边界（?=0）条件下的场分布，即 这样就把原来难于处理的开放边界问题转化为电壁、磁壁类边界问题。 这一方法只能用于线性系统。而且当结构具有n个开放边界时，必须计算2n个等效结构的场分布。 边界条件—开放边界 处理方法之二： （n+1）?t 时刻，开放边界（i, j, k）上的值可用n?t 时刻开放边界上的值以及边界内（i, j, k–1）网格上（n+1）?t 时刻的场来近似。 假设开放边界在z=k?z处，且在这 一平面内只有向正z方向传播的波 vz表示成： 取一阶近似 将上式代入式（1）可得： 这就是开放边界的一阶近似表达式。由此把F 看成切向电场并对式（2）进行差分近似，可得到： 上式即为开放边界条件的一阶近似差分表达式。 边界条件—开放边界 处理方法之三：是对一阶近似差分表达式的改进 c0用vp来代替， 假设 可以通过选择?t与?eff使 M为整数。对式