电子科技大学：《计算电磁学 Computational Electronmagentics》课程教学资源（课件讲稿，有限差分法）第4章 时域有限差分法 I 4.1 FDTD基本原理 4.2 解的稳定性条件 4.3 非均匀网格 4.4 共形网格 4.5 半解析数值模型 4.6 良导体中的差分格式

4.1FDTD基本原理 ◆ Yee的差分格式 >每一个网格点上，各场分量的新值依赖于该点在前一时间步长时刻的值及该点周围邻近点上另一场量的场 分量早半个时间步长时刻的值 >在任一给定时刻，场分量的计算可一次算出一个点，或者采用P个并行处理器一次算P个点（并行算法） >蛙跳格式：通过FDTD差分格式，逐个时间步长对模拟区域各网格点的电、磁场交替进行计算，在执行到适 当的时间步数后，即可获得需要的时域数值结果。 HO. H12 30 >显式差分格式：由{信息直接逐点计算，无需矩阵求逆 7

7 4.1 FDTD基本原理  Yee的差分格式  每一个网格点上，各场分量的新值依赖于该点在前一时间步长时刻的值及该点周围邻近点上另一场量的场 分量早半个时间步长时刻的值  在任一给定时刻，场分量的计算可一次算出一个点，或者采用P个并行处理器一次算P个点(并行算法)  蛙跳格式：通过FDTD差分格式，逐个时间步长对模拟区域各网格点的电、磁场交替进行计算，在执行到适 当的时间步数后，即可获得需要的时域数值结果。  显式差分格式：由 信息直接逐点计算 ，无需矩阵求逆

4.1FDTD基本原理 956 环路积分解释 >积分形式的时域麦克斯韦方程：安培定律和法拉第定律 H 6+k+ 安培定律 E 环路C (,j,k+ ■场量在环路每一边中点的值等于场量在该边的平均值 H 右边=（一+A+H(+分k+4 -++》-,-+}y气e i-,k+ ■E,k+ 等于Ez在小面元S1的平均值：时间中点取在（n+1/2)△t 东边=+引别 fa-生- 2k+ 2k+2 8

8 4.1 FDTD基本原理  环路积分解释  积分形式的时域麦克斯韦方程：安培定律和法拉第定律  安培定律：  场量在环路每一边中点的值等于场量在该边的平均值  等于Ez在小面元S1的平均值; 时间中点取在（n+1/2)t 1 1 1 1 S C d d t        D S H l 环路 E i j k z , ,        1 2 0 1 1 1 , , , , 2 2 n n z z x y E i j k E i j k t                            左边 1 1 1 1 2 2 2 2 1 0 1 1 1 1 1 1 1 1 , , , , , , , , 1 1 2 2 2 2 2 2 2 2 , , , , 2 2 n n n n yyxx n n z z H i j k H i j k H i j k H i j k t E i j k E i j k  x y                                                                    1 1 1 1 , , , , 2 2 2 2 1 1 1 1 , , , , 2 2 2 2 x y x y H i j k x H i j k y H i j k x H i j k y                                         右边 H x 1 1 ( , , ) 2 2 i j k   H x 1 1 ( , , ) 2 2 i j k   H y H y 1 1 ( , , ) 2 2 i j k   1 1 ( , , ) 2 2 i j k   x y z E y E y 1 ( , , ) 2 i j k  1 ( , , 1) 2 i j k   E z 1 ( , , ) 2 i j k  E z 1 ( , 1, ) 2 i j k   C1 1 S

4.1FDTD基本原理 966 环路积分解释 +1+ >积分形式的时域麦克斯韦方程：安培定律和法拉第定律 E 法拉第定律： (i+ 环路C E ■场量在环路每一边中点的值等于场量在该边的平均值 H 右边=-+2-E(+l1+2Ay 5,j+1k +++E(*小 H++人等于Hz在小面元S2的平均值；时间中点取在nAt 另一条获得差分格式的途径 可将环路路径选取来与弯曲表 左边：(++小(++ 面、槽缝等结构共形相配 tic -E;

9 4.1 FDTD基本原理  环路积分解释  积分形式的时域麦克斯韦方程：安培定律和法拉第定律  法拉第定律：  场量在环路每一边中点的值等于场量在该边的平均值  等于Hz在小面元S2的平均值; 时间中点取在 nt 2 2 2 2 S C d d t         B S E l 环路 1 1 , , 1, , 2 2 1 1 , 1, , , 2 2 x y x y E i j k x E i j k y E i j k x E i j k y                                        右边 1 1 , , 2 2 H i j k z         1 1 0 2 2 1 1 1 1 , , , , 2 2 2 2 n n z z x y H i j k H i j k t                               左边 1 1 2 2 0 1 1 1 1 , 1, , , , , 1, , 1 1 1 1 2 2 2 2 , , , , 2 2 2 2 n n n n x x y y n n z z E i j k E i j k E i j k E i j k t H i j k H i j k  y x                                                                     E x 1 ( , , ) 2 i j k  E x 1 1 ( , 1, ) 2 2 i j k    H y H y 1 ( , , ) 2 i j k  1 ( 1, , ) 2 i j k   x y E y z E y 1 1 ( , 1, ) 2 2 i j k    H z 1 1 ( , , ) 2 2 i j k   H z C2 1 ( , 1, ) 2 i j k   1 3 ( , , ) 2 2 i j k   2 S  另一条获得差分格式的途径  可将环路路径选取来与弯曲表 面、槽缝等结构共形相配

4.2解的稳定性条件 956 ◆不稳定性表现：随着时间步数的增加，计算结果将无限制地增加 ◆稳定性条件：△t和X、△y和△z必须满足一定的关系 ◆稳定性条件推导思路 >任何波都能展开为平面波谱的迭加 >考察算法对一平面波是否稳定 >把有限差分算式分解为时间本征值问题和空间本征值问题 >将平面波本征模带入空间本征值方程→本征值入 >入应在由时间本征值方程确定的稳定区间 10

10 4.2 解的稳定性条件  不稳定性表现：随着时间步数的增加，计算结果将无限制地增加  稳定性条件：t和x、y和z必须满足一定的关系  稳定性条件推导思路  任何波都能展开为平面波谱的迭加  考察算法对一平面波是否稳定  把有限差分算式分解为时间本征值问题和空间本征值问题  将平面波本征模带入空间本征值方程本征值  应在由时间本征值方程确定的稳定区间