FDTD原理及例子解析.ppt

编程思路（二维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边界的区域内的随时间变化的仿真分析 预祝大家儿童节快乐 报告人：邵海峰 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算法 Ey ( i, j, k ) x y z Ez Ey Ey Ez Ez Hx Ex Ex Ex Hz Hy 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 nnmax 结 束 No Yes 三维Maxwell方程的Yee算法 介绍了求解矢量Maxwell方程的FDTD Yee算法，归纳起来，Yee算法的主要特点有： 1）Yee算法采用耦合的Maxwell旋度方程，同时在时间和空间求解电场和磁场，而不是采用波动方程只求解电场或磁场。 2） Yee网格在三维空间这样安排E和H分量，使得每一个E或H分量由四个H或E循