fdtd基本原理解析.ppt

( 44a) ( 44b) ( 44c) ( 44d) ( 44e) ( 44f) 在上组公式中，系数 分别表示各场量位置处的媒质参数， 为各场量所对应的数组下标。 ( 45) 上述时域有限差分方程表明，任何时刻的电磁场取决于上一时间步的电磁场、与此电磁场正交的面上前半个时间步相邻的磁电场以及媒质参数。由于采用了中心差分近似，时域有限差分方程在空间和时间上具有二阶精度。 作为时域方法，时域有限差分法把所有研究的电磁问题作为初值问题，初始时刻模拟区域内的电磁场为零，在源激励下，以蛙跳的方式迭代时域有限差分方程，在时间上逐步向前推进电场和磁场。随着时间的发展，在有限计算区域内，时间和空间上离散取样电磁场量，数值模拟电磁波传播以及与媒质间的相互作用，近似实际连续的电磁波，获得整个计算区域内时域电磁信息。 为了方便编程计算，用计算机语言表示上述六个差分方程，得到： ( 46a) ( 46b) ( 46c) ( 46d) ( 46e) ( 46f) ?图9形象的给出了“计算空间的Yee网格离散化——电磁场分量空间分布——编程实现时各分量表示”这一实现过程 有助于我们清晰的理解程序上FDTD的实现方法。 图9三维空间内Yee网格的划分 依据这些三维空间内差分方程组可得出计算电磁场的时域推进计算方法，其流程图可以参考图6给出的流程图。 3. 媒质不连续性处理 介质交界面的边界条件要求电场E和磁场H在穿过边界时连续。假设计算目标非磁性材料，即 ，因此交界面总位于电场分量交界面。相应的电场迭代方程也应该做些修正，以考虑相邻网格不同的介电常数和电导率。 接下来我们推导有四种不同介质网格包围的电场Ex分量的迭代公式。假设环线C围绕 分量。环线C包含四个路径， ,每一个路径落入一个网格内，令S1处于网格Cube(i,j,k)内，S2处于网格Cube(i,j-1,k)内，S3处于网格Cube(i,j-1,k-1)内，S4处于网格Cube(i,j,k-1)内，如10图所示： 图 10 媒质不连续性处理 沿着路径C应用积分形式安培－麦克斯韦电路定律 ，可得： ( 47a) 利用时间导数的中心差分近似，这个方程可以近似为： ( 47b) 重新组织上式各项，可以得到由四种不同媒质包围的电场Ex电场分量的一般迭代公式： ( 47c) 其中： 和。 类似可得Ey与Ez分量。概括而言，我们可以通过介电常数和电导率的平均值来修改电磁场量的迭代公式。(47c)可以自动保证切向电场连续，因为穿过交界面时磁导率没有变化，这也保证了磁场的连续性。 不连续的磁导率可以按照类似于介电常数和电导率按照平均的方法处理，但必须指出：目前的空间网格结构不能容许介电常数与磁阻率同时不连续（1994 ----Hockanson）。 FDTD基本原理 直角坐标系中的时域有限差分方法 简单地说，将麦克斯韦方程组中的旋度方程分解为6个（电磁场各三个）标量方程以后，运用差分近似来代替各微分运算，即可得到麦克斯韦方程的时域有限差分计算格式。 首先回顾麦克斯韦方程组： 1. 一维情形 (1) 其旋度方程可分别写为： ( 2) 相应的本构关系为： ( 3) 则可以得到: ( 4) 特殊地，真空中的本构关系为： 则可知： ( 6) ( 5) 这种电磁场之间的耦合关系可以用图1表示为： 图1 电磁场耦合变化关系示意图 为简单起见，我们首先考虑1维情况下麦克斯韦旋度方程以及其FDTD差分计算。由真空中的本构关系可得到： ( 7) 可得到1维麦克斯韦方程组中的旋度方程为： 设： 对上式中的空间变量以及时间变量分别利用具有二阶精度的中心差分近似： 则由(8a)式，将磁场分量与电场分量在不同位置处离散，如图2所示。 (a) 电场与其周围磁场 (b) 磁场与其周围电场 图 2 电磁场空间离散位置示意图 (c)电磁场分量在空间上的分布 =1,