第2讲-Maxwell方程Yee算法

第2讲-Maxwell方程Yee算法第一页，共29页。2.1一维Maxwell方程的Yee算法（1）一维Maxwell方程利用一阶导数的二阶中心差分近似，上面的方程变为第一页第二页，共29页。2.1一维Maxwell方程的Yee算法（2）采用归一化磁场使得电场与归一化磁场有相同的数量级，于是可以得到FDTD迭代公式为式中，为自由空间中的光速。第二页第三页，共29页。2.1一维Maxwell方程的Yee算法（3）用计算机语言表示的FDTD公式式中，时间变量已隐含在迭代公式中，以及

只要给定了所有空间点上电/磁场的初值，就可以一步一步地求出任意时刻所有空间点上的电/磁场值。第三页第四页，共29页。2.1一维Maxwell方程的Yee算法（3）电场与磁场分量的空间-时间分布图0123第四页第五页，共29页。2.1一维Maxwell方程的Yee算法（4）Mainloopin1DFDTDC-program:

for(k=0;k<=kmax;k++){Hy[k]=0;Ex[k]=0;}for(n=1;n<=nmax;n++){Ex(0)=Source(n);for(k=0;k<kmax;k++){Hy[k]=Hy[k]-ca[k]*(Ex[k+1]-Ex[k]);}for(k=1;k<kmax;k++){Ex[k]=Ex[k]-cb[k]*(Hy[k]-Hy[k-1]);}Ex(kmax)=Boundary;}第五页第六页，共29页。2.2三维Maxwell方程的Yee算法（1）

考虑非时变、线性、各向同性媒质填充的无源区域，Maxwell旋度方程为第六页第七页，共29页。2.2三维Maxwell方程的Yee算法（2）

以及第七页第八页，共29页。2.2三维Maxwell方程的Yee算法（3）Yee首先将空间按立方体分割，电磁场的六个分量在空间的取样点分别放在立方体的边沿和表面中心点上，电场与磁场分量在任何方向始终相差半个网格步长。Ey(i,j,k)xyzEzEyEyEzEzHxExExExHzHy第八页第九页，共29页。2.2三维Maxwell方程的Yee算法（4）在时间上，Yee把电场分量与磁场分量也差半个步长取样。xEzt=2tEzEzEzEzEzEzEzEzEzEzEzHyHyHyHyHyHyt=1.5tt=tt=0.5tt=03x2xx0第九页第十页，共29页。2.2三维Maxwell方程的Yee算法（5）于是，利用一阶导数的二阶中心差分近似就可以导出旋度方程的FDTD公式。如：第十页第十一页，共29页。2.2三维Maxwell方程的Yee算法（6）

将上述公式代入相应的方程，得第十一页第十二页，共29页。2.2三维Maxwell方程的Yee算法（7）

采用时间平均近似最后，忽略高次项，得第十二页第十三页，共29页。2.2三维Maxwell方程的Yee算法（8）同理，可以得到其他2个磁场分量的FDTD方程第十三页第十四页，共29页。2.2三维Maxwell方程的Yee算法（9）

上面公式之间有明显的规律，便于记忆，如：系数Da和Db在空间的位置就是方程左边项的场分量的空间位置；右边第一项的场分量与左边的相同，但为n时间步，而左边场分量的时间步为n+1；右边第二、第三项的场分量与左边的相反（电场与磁场）三者的坐标分量满足循环关系：x-y-z-x;右边第二、第三项为空间差分形式。右边第二项分子上场量的坐标分量与分母上空间步长的坐标分量也满足x-y-z-x的循环关系。而右边第三项分子上场量的坐标分量和空间步长坐标分量与第二项恰好对调。第三项符号为负。第十四页第十五页，共29页。2.2三维Maxwell方程的Yee算法（10）利用对偶原理：，并注意到E与H在时间上差半个步长，可以直接从磁场FDTD公式得到电场的FDTD公式。如：第十五页第十六页，共29页。2.2三维Maxwell方程的Yee算法（11）给出n=0时刻电磁场的初值和媒质参数；由磁场FDTD公式，根据n时间步的电场值和磁场值求得n+1时间步空间所有点的磁场分量；由电场FDTD公式，根据n时间步的电场值和n+1时间步的磁场值求得n+1时间步空间所有点的电场分量；如此迭代，可获得任何时刻空间所有点的电磁场值。每一过程常称为蛙跳法(leapfrog)。对于连续变化的媒质，FDTD法需储存的量有n时间步和n+1时间步的六个场分量，六个D参数和六个C参数，所以总储存量近似为24N,N为空间网格数。但对于均匀媒质，D参数和C参数为常数，故总储存量减少为12N。第十六页第十七页，共29页。2.2三维Maxwell方程的Yee算法（12）媒质参数赋值在所有空间点给电磁场分量赋初值求所有空间离散点上n+1时间步的磁场求所有空间离散点上n+1时间步的电场n=n+1n>nmax结束NoYes第十七页第十八页，共29页。2.3以积分形式的Faraday和Ampere定理解释Yee算法（1）上面介绍的FDTD算法是从点的观点对Maxwell方程微分形式中的两个旋度方程直接进行导数二阶中心差分近似得到的。这种观点对理解FDTD如何模拟波在媒质中的传播是有用的。但是，当模拟细几何结构如导线、槽和曲面时，点的观点对于指导为了获得适当解需要作怎样的算法修正却帮助甚少。为了解决这一问题，我们从积分形式的Ampere和Faraday定理出发来解释Yee算法。第十八页第十九页，共29页。2.3以积分形式的Faraday和Ampere定理解释Yee算法（2）仅讨论自由空间的情况。考虑右图中实线网格的y-z平面上包围面积S1的矩形围线C1。沿C1应用Faraday定理:(i+1/2,j+1/2,k-1/2)C2S2Ey(i,j,k)C1S1EzEyEzHxHxHxHzHz第十九页第二十页，共29页。2.3以积分形式的Faraday和Ampere定理解释Yee算法（3）假设在围线任一边上和围线所围的面积上相应的场值不变，时间导数采用中心差分近似，则Faraday定理近似为整理后便可得到自由空间中FDTD公式。第二十页第二十一页，共29页。2.3以积分形式的Faraday和Ampere定理解释Yee算法（4）以相同的方式，把Ampere定理应用于图中虚线网格的x-z平面上包围面积S2的矩形围线C2，并作类似的假设，也可以得到相应的FDTD公式。所以，FDTD公式既是微分形式的Maxwell旋度方程的中心差分近似，也自然满足积分形式的Ampere和Faraday定律。第二十一页第二十二页，共29页。

2.4Yee算法的无散性（1）对于无源区域，满足Maxwell两个旋度方程的场也一定满足Maxwell的两个散度方程或它们的积分形式

下面证明对于从旋度方程近似而来的FDTD公式也满足两个散度方程。第二十二页第二十三页，共29页。

2.4Yee算法的无散性（2）在自由空间的一个Yee网格上考虑，有利用磁场分量的FDTD公式，把与磁场分量时间导数相关的电场空间有限差分代入上式中各项，可得第二十三页第二十四页，共29页。

2.4Yee算法的无散性（3）于是设初始时磁场为零，则

所以，对于无源区域，FDTD公式满足磁场Gauss定理，即对于磁场是无散的。同理可以证明，对于电场，FDTD公式也满足Gauss定理，即电场也是无散的。第二十四页第二十五页，共29页。结论2（1）本讲介绍了求解矢量Maxwell方程的FDTDYee算法，归纳起来，Yee算法的主要特点有：Yee算法采用耦合的Maxwell旋度方程，同时在时间和空间求解电场和磁场，而不是采用波动方程只求解电场或磁场。同时使用E和H信息比只使用其中一个的优点是获得的解更稳固（robust），即算法可以适用非常广泛的电磁波物理结构，并且电场和磁场的特性可以用更直接的方式模拟。如果同时使用电场和磁场，每一种场的独立特性，如边沿和角处切向磁场的奇异性、细线附近磁场的奇异性以及靠近点、边沿和细导线处径向电场的奇异性就能够独立地模拟。第二十五页第二十六页，共29页。结论2（2）Yee网格在三维空间这样安排E和H分量，使得每一个E或H分量由四个H或E循环的分量所环绕。这提供了一幅三维空间中由相互交链的Faraday定理和Ampere定理围线阵列构成的优美而简单的图画。保证了Yee算法同时模拟了Maxwell方程点意义上的微分形式和宏观的积分形式。后者对于处理边界条件和奇异性是极其有用的。旋度算子中空间导数的差分公式是二阶精度的中心差分。如果不同材料的交界面平行于Yee网格的一个坐标轴，在交界面上切向E和H的连续性自然保持。在Yee算法隐含地执行了两个高斯定律。所以，同时保证了无源区域中电磁场的无散性。第二十六页第二十七页，共29页。结论2（3