FDTD原理及例子课件.ppt

《FDTD原理及例子课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、FDTD数值分析法目录1.麦克斯韦方程的基础知识2.一维和三维Maxwell方程的Yee算法3.数值稳定性分析4.吸收边界条件5.波源的设置6.编程思路基础知识麦克斯韦方程微分形式：BPM方式基础知识麦克斯韦方程微分形式：FDTD方式将时间进行差分，并且磁场与电场交替迭代更新时谐场形式：对于有耗媒质：一维Maxwell方程的Yee算法一维Maxwell方程，介质参数和场量均与x，y无关(无损，电导率和磁导率为0)利用一阶导数的二阶中心差分近似，上面的方程变为一维Maxwell方程的Yee算法采用归一化磁场使得电场与归一化磁场有相同的数量级，于是可以

2、得到FDTD迭代公式为式中，为自由空间中的光速。一维Maxwell方程的Yee算法用计算机语言表示的FDTD公式式中，时间变量已隐含在迭代公式中，以及只要给定了所有空间点上电/磁场的初值，就可以一步一步地求出任意时刻所有空间点上的电/磁场值。一维Maxwell方程的Yee算法0123三维Maxwell方程的Yee算法Ey(i,j,k)xyzEzEyEyEzEzHxExExExHzHyFDTD离散中的Yee元胞1.每一个磁场分量由四个电场分量环绕；每一个电场分量由四个磁场分量环绕2.这种空间取样方式符合法拉第感应定律和安培环路定律3.电场和磁场在时间

3、上交替抽样，抽样时间彼此相差半个时间步三维Maxwell方程的Yee算法Maxwell旋度方程为三维Maxwell方程的Yee算法采用时间平均近似最后得：可以看到有个2阶小量，这个要忽略掉，因此是近似的三维Maxwell方程的Yee算法同理，可以得到其他2个磁场分量的FDTD方程三维Maxwell方程的Yee算法利用对偶原理：，并注意到E与H在时间上差半个步长，可以直接从磁场FDTD公式得到电场的FDTD公式。如：三维Maxwell方程的Yee算法媒质参数赋值在所有空间点给电磁场分量赋初值求所有空间离散点上n+1时间步的磁场求所有空间离散点上n+1

4、时间步的电场n=n+1n>nmax结束NoYes三维Maxwell方程的Yee算法介绍了求解矢量Maxwell方程的FDTDYee算法，归纳起来，Yee算法的主要特点有：1）Yee算法采用耦合的Maxwell旋度方程，同时在时间和空间求解电场和磁场，而不是采用波动方程只求解电场或磁场。2）Yee网格在三维空间这样安排E和H分量，使得每一个E或H分量由四个H或E循环的分量所环绕。3）Yee算法以蛙跳算法在时间上安排E和H分量。在某一时刻，使用前一时刻的E数据计算所有H分量。然后，再使用刚计算的H数据计算所有的E分量。如此循环，直至完成时间步进过程。数

5、值稳定性问题（1）FDTD计算中每一步都是有误差的，随着时间步进，误差会不断积累。如果误差的积累不会造成总误差的增加，就成FDTD法是稳定的，否则成为不稳定的。数值不稳定性会造成计算结果随时间步进无限增加。（2）FDTD法是有条件稳定的，即：时间步必须必须小于一定值以避免数值不稳定性。（3）数值稳定性分析方法是建立在Courant等人几十年前提出的经典方法基础上。这种方法首先把有限差分算法分解为相互分离的时间和空间本征值问题。数值稳定性问题二维的Yee算法数值Courant稳定性条件三维的Yee算法数值Courant稳定性条件便于理解，当Yee元胞

6、为立方体时，此两式表明时间间隔必须等于或小于波以光速通过yee元胞对角线长度1/3或1/2所需时间差分近似所带来的数值色散一维标量波动方程为例设在离散空间点,离散行波解为将上式代入差分方程得：最后得色散关系我们想要k与没加元胞之前的k一样，但是差分近似后k变的不一样差分近似所带来的数值色散将平面波带入差分方程所出现的稳定性、色散、以及各向异性，这些特性并非由介质的物理特性所引起，而是数值计算中的差分近似所致，在FDTD数值计算中，稳定性、色散、各向异性将影响计算精度。吸收边界条件理论上说，求解空间是无限大的，但是由于计算的数据容量问题，需要在有限空

7、间的周围做特殊处理，使得向边界面行进的波在边界处保持“外向行进”，无明显的反射现象，并且不会使内部空间的场产生畸变。1）Mur吸收边界条件2）Berenger完美匹配层（PML）3）各向异性匹配层吸收边界条件Berenger完美匹配层（PML）只考虑二维TE情况，对于二维TM和三维情况可采用类似方法进行分析。分析范围在PML层内分别表示自由空间中的电导率和磁阻率分别表示自由空间中的电导率和磁阻率Berenger为了引入规定损耗和阻抗匹配的新自由度，将Hz分裂为两个分量和，即同时引入了新的电导率和磁损耗，并规定TE情形的四个场分量（而不是通常的3个）

8、由下列方程耦合在一起：HzHzBerenger完美匹配层（PML）换句话说，Berenger构造了一般新的非物理媒质（称为