时域有限差分法

时域有限差分法第1页，共23页。引言引言11966年,K.S.Yee美籍香港人首先提出了Finite-Difference Time-Domain Method,并用于柱形金属柱电磁散射分析。由于当时计算机技术还比较落后，这一方法并未引起重视。1972年,A.Taflovey应用FDTD研究了UHF和微波对人类眼睛的穿透，以了解“微波白内障的成因。Taflove成功地应用和开展了Yee的FDTD算法。80年代后期，随着高速大容量计算机的普及，FDTD法得到了迅速开展。如今已应用于涉及波动现象的任何领域。至今，FDTD法的研究与应用仍方兴未艾。第2页，共23页。引言引言2 本本课课程程采采用用研研讨讨班班形形式式。教教师师讲讲授授FDTDFDTD的的根根本本知知识识，学生针对某一方向进展较深入的研究。学生针对某一方向进展较深入的研究。本本讲讲我我们们考考虑虑描描述述波波动动现现象象的的最最根根本本偏偏微微分分方方程程：一一维维标标量量波波动动方方程程的的数数值值FDTDFDTD解解，为为以以后后二二维维、三三维维MaxwellMaxwell方程的方程的FDTDFDTD分析奠定根底分析奠定根底 课程内容取自以下的参考书和近年来相关的一些文献课程内容取自以下的参考书和近年来相关的一些文献 1A.Taflove,Computational 1A.Taflove,Computational Electrodynamics Electrodynamics The The Finite-Difference Finite-Difference Time-Domain Time-Domain Method,Method,Artech Artech HourseHourse，1995.1995.22高高本本庆庆，时时域域有有限限差差分分法法，国国防防工工业业出出版版社社，1995.1995.33葛葛德德彪彪，闫闫玉玉波波，电电磁磁场场时时域域有有限限差差分分法法，西西电电出版社，出版社，20022002第3页，共23页。1.1 1.1 差分近似差分近似1 1一维标量波动方程 1-1上式的解为 1-2 采用Taylor 展开1-3第4页，共23页。1.1 1.1 差分近似差分近似2 2于是，有 (1-4同理，有 (1-5上式称为二阶偏导数的二阶中心差分格式。将它们代入(1-1)，得1-6忽略高次项，便可得到求解的差分迭代公式。第5页，共23页。1.1 1.1 差分近似差分近似3 3NoYesn=0在所有空间点给uin,uin-1(i=1:imax)赋初值 n=n+1由(1-6)在所有空间点求uin+1(i=1:imax)结 束nnmax?图1.1 一维波动方程FDTD流程图第6页，共23页。1.1 差分近似4 应当注意，在一般情况下(1-6)对时间或空间具有二阶精度。但对 于 的特殊情况，根据解(1-2)，可以证明 于是 所以，1-6中的两个剩余项抵消，得到了准确的数值差分公式 1-7正因为有这样的奇妙特性，为“魔时间步(Magic time step).第7页，共23页。1.2 1.2 数值色散关系数值色散关系(1)(1)色散关系定义为行波的波长随频率的变化关系。为方便起见，色散关系也常表示为行波的波数关于角频率的变化关系。考虑(1.1)的正弦行波解 代入(1-1)得 即 (1-8)上式便是一维标量波动方程的色散关系。由上式得相速度 1-9可见，相速与频率无关，称为非色散。非色散意味着对于具有任意调制的包络或脉冲形状的波传播任意距离后波形保持不变。进一步由(1-8)可以得到群速关系 (1-10)这种情况下，群速也是与频率无关。第8页，共23页。1.2 1.2 数值色散关系数值色散关系(2)(2)上述过程也可用于一维标量波动方程差分近似的数值色散分析。设在离散空间点 ,离散行波解为 ，式中，为存在于有限差分网格中的数值正弦波的波数。一般情况下，不同于连续物理波的波数。正是这种不同导致了数值相速和群速偏离了准确解。进而导致了数值色散误差。将上式代入差分方程(1-6),得 (1-11)重新组合并应用 Euler恒等式，最后得到数值色散关系为 (1-12)第9页，共23页。1.3 1.3 数值相速数值相速1 1 类似于(1-9)，定义数值相速为 由1-12可得1-13可见数值相速与频率有关。因此，由FDTD得到的数值波是色散的。取 那么数值相速为 。相对误差为-1.27%。如果物理波传播了 距离(100空间格)时，数值模拟波只传播了0。取 那么 。这时数值相速的相对误差为0.31，减少了4倍。同样，当物理波传播了同样的 时(200空间格)，数值模拟传播了0，也减少了4倍。误差减少了4倍反映了差分算法是二阶精度的。第10页，共23页。1.3 1.3 数值相速数值相速2 2情况情况1 1：非常细网格：非常细网格 根据根据 ，数值色散关系，数值色散关系(1-12)(1-12)变为变为即，即，最后得，最后得 ，于是有，于是有 。所以，在非常细的网格条件下，差分解逼近准确解。所以，在非常细的网格条件下，差分解逼近准确解。情况情况2 2：魔时间步魔时间步 (1-12)(1-12)变为变为 ,即即 。所以，所以，。可见，魔时间步下差分解与准确解一样。可见，魔时间步下差分解与准确解一样。第11页，共23页。1.4 1.4 数值群速数值群速 定义数值相速为 1-14 情况1 非常细网格 利用正弦函数的一阶Taylor展开，可得 1-15 所以，群速与相速一样，在细网格条件下趋近准确解。这证明了 当空间步长和时间步长趋于零时，数值解变得准确。情况2 魔时间步 将魔时间步条件和波数代入(1-14)，得 1-16 再次验证了魔时间步下数值解等于准确解。第12页，共23页。1.5 1.5 数值稳定性数值稳定性1 1FDTD计算中每一步都是有误差的，随着时间步进，误差会不断积累。如果误差的积累不会造成总误差的增加，就成FDTD法是稳定的，否那么成为不稳定的。数值不稳定性会造成计算结果随时间步进无限增加。FDTD法是有条件稳定的，即：时间步必须必须小于一定值以防止数值不稳定性。本节的数值稳定性分析方法是建立在Courant等人几十年前提出的经典方法根底上。这种方法首先把有限差分算法分解为相互别离的时间和空间本征值问题。第13页，共23页。1.5 1.5 数值稳定性数值稳定性1 1时间本征值问题时间本征值问题 (1-17)差分近似，得 (1-18)定义不变增长因子(1-19)第14页，共23页。1.5 1.5 数值稳定性数值稳定性2 2将(1-19)代入(1-18)，有 ，于是 算法稳定性要求 。如果 ，那么总有 ，于是 ，满足稳定性要求。这样可得1-20这就是稳定的数值差分解所要求的时间本征值谱。第15页，共23页。1.5 1.5 数值稳定性数值稳定性3 3空间本征值问题空间本征值问题 (1-21)代入中心差分公式，得 (1-22)令 ，Eular公式可得因为 ，所以 (1-23)上式给出了差分网格中任意空间Fourier模的本征值谱。第16页，共23页。1.5 1.5 数值稳定性数值稳定性4 4稳定性稳定性 为了保证任何空间模式的数值稳定性，(1-23)给出的空间模式的本征值范围必须完全落在(1-20)所给出的时间本征值的稳定范围内，于是即 (1-24)可见，时间步长 必须是有界的。上式称为Courant稳定性条件。有趣的是其上界恰好是魔时间步。第17页，共23页。1.6 1.6 鼓励源的设置鼓励源的设置 在FDTD模拟电磁波传播时需要设置初始条件和鼓励源。最简单的源设置方法是“硬源,即在鼓励源的位置令 u满足ui=f(n),常用的有 正弦函数 ui=sin(nt+)高斯函数 ui=exp-(n-n0)2/T2 阶跃函数 ui=0 nn1 =(n-n1)/(n2-n1)n1nn2 “硬源设置简单，但当反射波回到“硬源位置时，会引起寄生反射，所以，要在这之前“关掉源。以后会有有关源设置的更详细讨论。第18页，共23页。1.7 1.7 吸收边界条件吸收边界条件 由于计算机容量所限，计算域必须是有限的。对于理想电壁或磁壁的边界条件的设置是直接的。但如果模拟的是“开问题，就要设置截断边界。在截断边界上要设置吸收边界条件，使得电磁波可以被完全吸收，模拟波无反射的通过吸收边界。对于一维问题，采用单向波方程 于是利用单向差分近似得到吸收边界条件，详细讨论见后面章节。第19页，共23页。结结 论论 1 1 本讲介绍了一维标量波动方程的FDTD求解过程：利用Taylor级数展开方法获取空间/时间导数的二阶中心差分近似，从而得到具有二阶精度的方程数值解的时间步进迭代公式。一般情况下，数值解引入了寄生的数值色散。当空间步长和时间步长非常小时，数值解逼近准确解。当时间步长满足魔时间步条件时，数值解等于准确解。空间步长和时间步长必须满足Courant稳定性条件才能保证数值解的稳定性。第20页，共23页。习习 题题 1 11.1 利用Taylor级数展开方法分别推导一阶导数 的二阶和四阶精度中心差分近似。1.2 利用数值相速和群速公式分别画出数值相速和群速在 ，和 条件下关于网格空间分辨率 的曲线，并进展相应的讨论。1.3 编写本讲介绍的一维标量波动方程FDTD求解程序。设在网格左边界波源为以下脉冲函数：a.高斯函数 b.矩形函数 在网格右边界 。时间总步数为5000,空间总步为200。在 ，和 条件下每隔1000个时间步画出每个传播脉冲关于位置的分布曲线。第21页，共23页。谢谢谢谢欣欣欣欣赏赏第22页，共23页。人有了知识，就会具备各种分析能力，明辨是非的能力。所以我们要勤恳读书，广泛阅读，古人说“书中自有黄金屋。”通过阅读科技书籍，我们能丰富知识，培养逻辑思维能力；通过阅读文学作品，我们能提高文学鉴赏水平，培养文学情趣；通过阅读报刊，我们能增长见识，扩大自己的知识面。有许多书籍还能培养我们的道德情操，给我们巨大的精神力量，鼓舞我们前进。第23页，共23页。