FDTD原理及例子

FDTD数值分析法目录1.麦克斯韦方程的基础知识2.一维和三维Maxwell方程的Yee算法3.数值稳定性分析4.吸收边界条件5.波源的设置6.编程思路基础知识麦克斯韦方程微分形式：0DHJtBEtBDBPM方式基础知识麦克斯韦方程微分形式：0DHJtBEtBDFDTD方式将时间进行差分，并且磁场与电场交替迭代更新,,eemmmeDHJJEtBEJJsHtBD时谐场形式：00()(,,,)(,,,)(,,,)(,,,)jtjtHjjEEjHExyztExyzteHxyztHxyzte对于有耗媒质：一维Maxwell方程的Yee算法一维Maxwell方程，介质参数和场量均与x，y无关(无损，电导率和磁导率为0)利用一阶导数的二阶中心差分近似，上面的方程变为zHtEyrx01zEtHxry01zkHkHktkEkEnynyrnxnx)()()(1)()(212/1212/101zkEkEktkHkHnxnxrnyny)()1()(1)()(210212/1212/1一维Maxwell方程的Yee算法采用归一化磁场使得电场与归一化磁场有相同的数量级，于是可以得到FDTD迭代公式为式中，为自由空间中的光速。HH00~)()1()()21(~)21(~212/12/1kEkEkztckHkHnxnxrnyny)21(~)21(~)()()(2/12/11kHkHkztckEkEnynyrnxnx001c一维Maxwell方程的Yee算法用计算机语言表示的FDTD公式式中，时间变量已隐含在迭代公式中，以及只要给定了所有空间点上电/磁场的初值，就可以一步一步地求出任意时刻所有空间点上的电/磁场值。][]1[*][][][kExkExkcakHykHy]1[][*][][][kHykHykcbkExkEx)(][;)(][];[][);21(~][21kztckcbkztckcakEkExkHkHyrrxy一维Maxwell方程的Yee算法xExExEyHyH0123k0n2n2/3n1n2/1n)()1()()21(~)21(~212/12/1kEkEkztckHkHnxnxrnyny)21(~)21(~)()()(2/12/11kHkHkztckEkEnynyrnxnx三维Maxwell方程的Yee算法Ey(i,j,k)xyzEzEyEyEzEzHxExExExHzHyFDTD离散中的Yee元胞1.每一个磁场分量由四个电场分量环绕；每一个电场分量由四个磁场分量环绕2.这种空间取样方式符合法拉第感应定律和安培环路定律3.电场和磁场在时间上交替抽样，抽样时间彼此相差半个时间步三维Maxwell方程的Yee算法Maxwell旋度方程为tHHE1xzyxHyEzEtHyxzyHzExEtH11zyxzHxEyEtHxyzxEzHyHtE1yzxyExHzHtE1zxyzEyHxHtE1tEEH三维Maxwell方程的Yee算法11221111,,,,2222211,,22nnnxxijkijkxijkHHHOttt2,21,1,21,21,21,zOzEEzEnkjiynkjiynkjiy221,,21,1,21,21,yOyEEyEnkjiznkjiznkjiz采用时间平均近似22121,21,2121,21,21,21,2tOHHHnkjixnkjixnkjixynkjiEnkjiEznkjiEnkjiEkjiDnkjiHkjiDnkjiHzzyybHxxaHxx,,,,,1,,,,,1,,,,,,,,,1,,,最后得：可以看到有个2阶小量，这个要忽略掉，因此是近似的三维Maxwell方程的Yee算法同理，可以得到其他2个磁场分量的FDTD方程znkjiEnkjiExnkjiEnkjiEkjiDnkjiHkjiDnkjiHxxzzbHyyaHyy,,,,1,,,,,,,,1,,,,,,,1,,,xnkjiEnkjiEynkjiEnkjiEkjiDnkjiHkjiDnkjiHyyxxbHzzaHzz,,,,,,1,,,,,1,,,,,,,,1,,,所在空间位置wHaHwttD22所在空间位置wHbHwttD22三维Maxwell方程的Yee算法利用对偶原理：，并注意到E与H在时间上差半个步长，可以直接从磁场FDTD公式得到电场的FDTD公式。如：，，HEEH,,bbaaCDCD,ynkjiHnkjiHznkjiHnkjiHkjiCnkjiEkjiCnkjiEzzyybExxaExx1,,,1,,1,1,,,1,1,,,,,,,,,1,,,所在空间位置wEaEwttC22所在空间位置wEbEwttC22三维Maxwell方程的Yee算法媒质参数赋值在所有空间点给电磁场分量赋初值求所有空间离散点上n+1时间步的磁场求所有空间离散点上n+1时间步的电场n=n+1nnmax结束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分量。如此循环，直至完成时间步进过程。数值稳定性问题（1）FDTD计算中每一步都是有误差的，随着时间步进，误差会不断积累。如果误差的积累不会造成总误差的增加，就成FDTD法是稳定的，否则成为不稳定的。数值不稳定性会造成计算结果随时间步进无限增加。（2）FDTD法是有条件稳定的，即：时间步必须必须小于一定值以避免数值不稳定性。（3）数值稳定性分析方法是建立在Courant等人几十年前提出的经典方法基础上。这种方法首先把有限差分算法分解为相互分离的时间和空间本征值问题。数值稳定性问题22111yxct二维的Yee算法数值Courant稳定性条件三维的Yee算法数值Courant稳定性条件2221111zyxct便于理解，当Yee元胞为立方体时，此两式表明时间间隔必须等于或小于波以光速通过yee元胞对角线长度1/3或1/2所需时间差分近似所带来的数值色散一维标量波动方程为例22222xuctu设在离散空间点,离散行波解为nitx,xiktnjninietxuu~,将上式代入差分方程得：22221112122tOxOctuuuuuxtcuninininininitjxkjxkjtjeeextce22~~2最后得色散关系1cos1cos1~21ttcxxk我们想要k与没加元胞之前的k一样，但是差分近似后k变的不一样差分近似所带来的数值色散将平面波带入差分方程所出现的稳定性、色散、以及各向异性，这些特性并非由介质的物理特性所引起，而是数值计算中的差分近似所致，在FDTD数值计算中，稳定性、色散、各向异性将影响计算精度。吸收边界条件理论上说，求解空间是无限大的，但是由于计算的数据容量问题，需要在有限空间的周围做特殊处理，使得向边界面行进的波在边界处保持“外向行进”，无明显的反射现象，并且不会使内部空间的场产生畸变。1）Mur吸收边界条件2）Berenger完美匹配层（PML）3）各向异性匹配层吸收边界条件Berenger完美匹配层（PML）只考虑二维TE情况，对于二维TM和三维情况可采用类似方法进行分析。分析范围在PML层内yHEtEzxx0xHEtEzyy0xEyEHtHyxzz*0*分别表示自由空间中的电导率和磁阻率分别表示自由空间中的电导率和磁阻率Berenger为了引入规定损耗和阻抗匹配的新自由度，将Hz分裂为两个分量和，即同时引入了新的电导率和磁损耗，并规定TE情形的四个场分量（而不是通常的3个）由下列方程耦合在一起：zxHzyzxzHHHzyHHzHzzHyx,**,yxBerenger完美匹配层（PML）yHHEtEzyzxxyx0xHHEtEzyzxyxy0xEHtHyzxxzx*0yEHtHxzyyzy*0换句话说，Berenger构造了一般新的非物理媒质（称为PML媒质），在该媒质中场满足左端的方程（并不一定是Maxwell方程）经过一系列公式推导：最后得出无反射匹配条件0*00*0yyxxGHEZ000022sincosyxWWG00Z表明PML媒质中的波阻抗与入射角无关，与真空中波阻抗相同，这意味着从真空中的任意角度入射到PML媒质交界面时将会无反射地进入PML媒质中，并在PML媒质中衰减地传播Berenger完美匹配层（PML）1.若分界面垂直于x轴（y轴），要求二者具有相同的横向电导率和磁导率，且横向和纵向电导率、磁导率均满足阻抗匹配条件。Berenger完美匹配层（PML）如何计算：任何普通介质均可以视为特殊的PML介质，但麦克斯韦方程只需计算三个分量，而不是PML的四个分量。另外由于在介质中电磁波衰减的很快，常规FDTD中Yee的差分格式已不再适用，需要指数差分0011122,1,1/2,1/2,11xxnnnnttyyzzijijijijxEeEeHHx0011122,12,12,1,11yynnnnttxxzzijijijijyEeEeHHy波源的设置为了用FDTD法模拟电磁场工程问题，必须在FDTD网格中引入电磁波激励源。常用的波源种类：平面波源—用于电磁散射问题；导波源—用于微波网络参数计算；电流源或电压源—用于微波电路或天线的激励。常用的设置方式：初始条件；硬源；电流源；总场/散射场公式硬波源硬波源可以通过规定FDTD空间网格中电场或磁场分量满足所希望的时间函数进行简单地设置。如：在网格点（is,js,ks)建立Ez的硬源为:硬波源相当于点源，模拟从源点向外辐射的具