FDTD原理及例子.ppt

单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,FDTD,数值分析法,目录,1.,麦克斯韦方程的基础知识,2.,一维和三维,Maxwell,方程的,Yee,算法,3.,数值稳定性分析,4.,吸收边界条件,5.,波源的设置,6.,编程思路,基础知识,麦克斯韦方程微分形式：,BPM,方式,基础知识,麦克斯韦方程微分形式：,FDTD,方式将时间进行差分，并且磁场与电场交替迭代更新,时谐场形式：,对于有耗媒质：,一维,Maxwell,方程的,Yee,算法,一维,Maxwell,方程，介质参数和场量均与,x,，,y,无关,(,无损，电导率和磁导率为,0),利用一阶导数的二阶中心差分近似，上面的方程变为,一维,Maxwell,方程的,Yee,算法,采用归一化磁场,使得电场与归一化磁场有相同的数量级，于是可以得到,FDTD,迭代公式为,式中，为自由空间中的光速。,一维,Maxwell,方程的,Yee,算法,用计算机语言表示的,FDTD,公式,式中，时间变量已隐含在迭代公式中，以及,只要给定了所有空间点上电,/,磁场的初值，就可以一步一步地求出任意时刻所有空间点上的电,/,磁场值。,一维,Maxwell,方程的,Yee,算法,0,1,2,3,三维,Maxwell,方程的,Yee,算法,E,y,(i,j,k),x,y,z,E,z,E,y,E,y,E,z,E,z,H,x,E,x,E,x,E,x,H,z,H,y,FDTD,离散中的,Yee,元胞,1.,每一个磁场分量由四个电场分量环绕；每一个电场分量由四个磁场分量环绕,2.,这种空间取样方式符合法拉第感应定律和安培环路定律,3.,电场和磁场在时间上交替抽样，抽样时间彼此相差半个时间步,三维,Maxwell,方程的,Yee,算法,Maxwell,旋度方程为,三维,Maxwell,方程的,Yee,算法,采用时间平均近似,最后得：,可以看到有个,2,阶小量，这个要忽略掉，因此是近似的,三维,Maxwell,方程的,Yee,算法,同理，可以得到其他,2,个磁场分量的,FDTD,方程,三维,Maxwell,方程的,Yee,算法,利用对偶原理：,，并注意到,E,与,H,在时间上差半个步长，可以直接从磁场,FDTD,公式得到电场的,FDTD,公式。如：,三维,Maxwell,方程的,Yee,算法,媒质参数赋值,在所有空间点给电磁场分量赋初值,求所有空间离散点上,n+1,时间步的磁场,求所有空间离散点上,n+1,时间步的电场,n=n+1,nn,max,结 束,No,Yes,三维,Maxwell,方程的,Yee,算法,介绍了求解矢量,Maxwell,方程的,FDTD Yee,算法，归纳起来，,Yee,算法的主要特点有：,1,）,Yee,算法采用耦合的,Maxwell,旋度方程，同时在时间和空间求解电场和磁场，而不是采用波动方程只求解电场或磁场。,2,）,Yee,网格在三维空间这样安排,E,和,H,分量，使得每一个,E,或,H,分量由四个,H,或,E,循环的分量所环绕。,3,）,Yee,算法以蛙跳算法在时间上安排,E,和,H,分量。在某一时刻，使用前一时刻的,E,数据计算所有,H,分量。然后，再使用刚计算的,H,数据计算所有的,E,分量。如此循环，直至完成时间步进过程。,数值稳定性问题,（1,）,FDTD,计算中每一步都是有误差的，随着时间步进，误差会不断积累。如果误差的积累不会造成总误差的增加，就成,FDTD,法是稳定的，否则成为不稳定的。数值不稳定性会造成计算结果随时间步进无限增加。,（2,）,FDTD,法是有条件稳定的，即：时间步必须必须小于一定值以避免数值不稳定性。,（,3,）数值稳定性分析方法是建立在,Courant,等人几十年前提出的经典方法基础上。这种方法首先把有限差分算法分解为相互分离的时间和空间本征值问题。,数值稳定性问题,二维的,Yee,算法数值,Courant,稳定性条件,三维的,Yee,算法数值,Courant,稳定性条件,便于理解，当,Yee,元胞为立方体时，此两式表明时间间隔必须等于或小于波以光速通过,yee,元胞对角线长度,1/3,或,1/2,所需时间,差分近似所带来的数值色散,一维标量波动方程为例,设在离散空间点,离散行波解为,将上式代入差分方程,得：,最后得色散关系,我们想要,k,与没加元胞之前的,k,一样，但是差分近似后,k,变的不一样,差分近似所带来的数值色散,将平面波带入差分方程所出现的稳定性、色散、以及各向异性，这些特性并非由介质的物理特性所引起，而是数值计算中的差分近似所致，在,FDTD,数值计算中，稳定性、色散、各向异性将影响计算精度。,吸收边界条件,理论上说，求解空间是无限大的，但是由于计算的数据容量问题，需要在有限空间的周围做特殊处理，使得向边界面行进的波在边界处保持“外向行进”，无明显的反射现象，并且不会使内部空间的场产生畸变。,1,）,Mur,吸收边界条件,2,）,Berenger,完美匹配层（,PML,）,3,）各向异性匹配层吸收边界条件,Berenger,完美匹配层（,PML,）,只考虑二维,TE,情况，对于二维,TM,和三维情况可采用类似方法进行分析。分析范围在,PML,层内,分别表示自由空间中的电导率和磁阻率,分别表示自由空间中的电导率和磁阻率,Berenger,为了引入规定损耗和阻抗匹配的新自由度，将,Hz,分裂为两个分量 和 ，即 同时引入了新的电导率 和磁损耗 ，并规定,TE,情形的四个场分量（而不是通常的,3,个）由下列方程耦合在一起：,Hz,Hz,Berenger,完美匹配层（,PML,）,换句话说，,Berenger,构造了一般新的非物理媒质（称为,PML,媒质），在该媒质中场满足左端的方程（并不一定是,Maxwell,方程）,经过一系列公式推导：最后得出无反射匹配条件,表明,PML,媒质中的波阻抗与入射角无关，与真空中波阻抗相同，这意味着从真空中的任意角度入射到,PML,媒质交界面时将会无反射地进入,PML,媒质中，并在,PML,媒质中衰减地传播,Berenger,完美匹配层（,PML,）,1.,若分界面垂直于,x,轴（,y,轴），要求二者具有相同的横向电导率和磁导率，且横向和纵向电导率、磁导率均满足阻抗匹配条件。,Berenger,完美匹配层（,PML,）,如何计算：任何普通介质均可以视为特殊的,PML,介质，但麦克斯韦方程只需计算三个分量，而不是,PML,的四个分量。,另外由于在介质中电磁波衰减的很快，常规,FDTD,中,Yee,的差分格式已不再适用，需要指数差分,波源的设置,为了用,FDTD,法模拟电磁场工程问题，必须在,FDTD,网格中引入电磁波激励源。,常用的波源种类：,平面波源,用于电磁散射问题；,导波源,用于微波网络参数计算；,电流源或电压源,用于微波电路或天线的激励。,常用的设置方式：,初始条件；,硬源；,电流源；,总场,/,散射场公式,硬波源,硬波源可以通过规定,FDTD,空间网格中电场或磁场分量满足所希望的时间函数进行简单地设置。如：在网格点（,i,s,j,s,k,s,),建立,E,z,的硬源为,:,硬波源相当于点源，模拟从源点向外辐射的具有与源函数性质相同的数值波。如果材料结构放在离波源一定距离的地方，辐射的数值波最终传播到这个结构，一部分传输过去，一部分反射。在原理上，时间步进可以连续进行，直到所有的瞬态波消失。,硬波源,常用源波形有以下三种,正弦波源为时谐源，提供了频率为,f,0,的正弦波。,FDTD,模拟时，一般采用,Gaussian,脉冲信号激励，因为它可以提供宽频带特性。,Gaussian,脉冲信号的频谱也是,Gaussian,函数。,硬波源,但是，在连续时间步进时，当反射波到达波源位置，由于波源处总场已被规定，并没有考虑网格中可能的反射波（所以称为硬源），硬源对这些反射波就会造成寄生的非物理的再次反射。当规定了一个表面上的总场而没有考虑反射场值时，寄生反射就一定会发生。,解决这一问题的一个简单方法是当反射波到达波源位置时去掉波源。这就要求源离散射体足够远，使散射波到达源位置时源入射波已基本为零，所以这种方法不适合正弦入射波源。,编程思路（二维,TE+PML,边界条件）,第一步：基本参数设置，包扩介电系数，磁导系数，频率，波长等,编程思路（二维,TE+PML,边界条件）,第二步：网格参数设置,编程思路（二维,TE+PML,边界条件）,第三步：,材料属性以及激励源设置,编程思路（二维,TE+PML,边界条件）,第四步：,初始场设置（无激励），矩阵形式,编程思路（二维,TE+PML,边界条件）,第五步：中心场区的系数设置,编程思路（二维,TE+PML,边界条件）,第六步：,PML,区的设置，指数差分系数设置,以,Front,层为例，按照相似的办法定义,back,，,left,，,right,，,PML,边界层,编程思路（二维,TE+PML,边界条件）,第七步：,movie,初始化,编程思路（二维,TE+PML,边界条件）,第八步：电磁场随时间迭代运算，之前都是初始化,从电场迭代到磁场迭代过程，中间还需要将,ex,在,PML,层更新,编程思路（二维,TE+PML,边界条件）,第九步：将电场在,PML,层迭代，以,Ex,在,PML,边界迭代为例,编程思路（二维,TE+PML,边界条件）,第九步：将电场在,PML,层迭代，以,Ex,在,PML,边界迭代为例,编程思路（二维,TE+PML,边界条件）,第十步：将磁场在,PML,层迭代，分,Hzx,和,Hzy,第十一步：将,n/4,为整数的时间帧保存到,M(:,nn)=getframe(gcf,rect);,最后,for,循环结束，用函数,movie(gcf,M,0,10,rect);,得到电磁场在以,PML,边界的区域内的随时间变化的仿真分析,预祝大家儿童节快乐,报告人：邵海峰,