第一部分 有限差分时域法课件.pptVIP

* * **满足CFL条件只能保证数值计算是稳定的, 并不能保证计算结果的精度. * **满足CFL条件只能保证数值计算是稳定的, 并不能保证计算结果的精度. * * * **满足CFL条件只能保证数值计算是稳定的, 并不能保证计算结果的精度. * * * * 由于采取了近似，采用Mur吸收边界条件进行FDTD运算时，只有当入射角θ较小时才能取得较为精确的结果。 为了克服这一缺点，Berenger(1994年)提出在截断边界处设置一特殊介质层（即完全匹配层，Perfectly Matched Layer），其特点为： PML为虚拟的有一定厚度（一般设为3-5个网格尺寸）的有耗介质，透射波会迅速衰减为零。 PML的波阻抗与相邻介质（自由空间）波阻抗完全匹配，这样入射波在任意角度下可无反射地进入PML。 3.2 完全匹配层PML(Perfectly Matched Layer) Berenger证明当介质的电导率σ和导磁率σ满足以下关系时 该介质为PML。 (3-20) 从源所产生的电磁场结构来分类，主要分为近场源（如电压源、电流源）和远场源两大类。我们这里只讨论远场源。 对激励源的模拟是FDTD方法分析电磁问题的一个重要任务。 第四章 激励源 从源随时间变化看，主要有两类激励源：一类是随时间周期变化的时谐场源，一类是对时间呈脉冲函数形式的波源。 从空间分布来看，有面源、线源和点源等。 为了用FDTD方法来计算时谐场情况下的电磁问题，激励源定义为： 即为从t=0开始的半无限正弦波。 4.1 随时间变化的源 4.1.1 时谐场源 (4-1) 高斯脉冲函数的时域表达式为 1.高斯脉冲 (4-2) 4.1.2 脉冲源 通过对(5-2)式进行傅立叶变换，可得高斯脉冲的频谱特性。 (4-3) 双指数脉冲的时域表达式为 2.双指数脉冲 (4-4) 当A0=5.25×10^4V/m, α=4×10^6/s, β=4.76×10^8/s时被称为Bell波形，它被经常用在核电磁脉冲和雷电脉冲的研究中。 (4-5) (4-4)的傅立叶变换为 为FDTD仿真构造的激励波形是时间的函数，而FDTD仿真主要输出的也是时域数据。使用傅立叶变换，输入和输出的时间函数就可以变换到频域，从而得到系统对各个频率的响应。 (4-6) 上式的逆傅立叶变换为 4.2 时域与频域的相互转换 (4-7) 1、考虑函数f(x)=cos(x)e^(-0.1x), 其一阶导数为f’(x)=-sin(x)e^(-0.1x)-0.1cos(x)e^(-0.1x)。试用matlab编程，做出(1)当取样周期Δx=π/5,在一张图中绘出一阶导数曲线，以及由前向差分、后项差分和中心差分来计算导数的近似曲线。(2)当取样周期Δx=π/10,在一张图中绘出一阶导数曲线，以及由前向差分、后项差分和中心差分来计算导数的近似曲线。(3)在同一张图中画出Δx=π/5时,前向差分、后项差分和中心差分来计算导数的误差曲线。(4)在同一张图中画出Δx=π/10时,前向差分、后项差分和中心差分来计算导数的误差曲线。 作业： 2: 两个相距1米的平行无限大金属平面位于x=0和x=1处. 在x=0.5处有一高斯激励源,,产生电场分量为Ez,磁场分量为Hy的TEM波.金属板中间为空气,试用matlab编写FDTD程序,计算并显示电磁场在金属板内的变化情况.(**此例不能直接用(2-31)(2-32)式,因为这两个式子适用电场分量为Ex,磁场分量为Hy的TEM波. 编写此例运算程序时,应按照(2-2)麦克斯韦标量方程推导出适用电场分量为Ez,磁场分量为Hy的一维TEM波FDTD方程) 电磁场分析方法一般分为解析法、近似解析法、数值方法 *解析解即为精确解。解析法包括分离变量法和变换数学法。前者针对微分方程，后者针对积分方程。解析法能得到带求函数的闭式解，但仅能解决几种简单经典的微波结构，如矩形、六面体、圆柱等。 **近似解析法包括变分法、微扰法、高频和低平近似法、直线法等。高频近似法如PO、GO、GTD等用于电大结构的电磁散射辐射。低频近似在电小结构时，用静电场近似简化。 * * * * * * * **满足CFL条件只能保证数值计算是稳定的, 并不能保证计算结果的精度. * **满足CFL条件只能保证数值计算是稳定的, 并不能保证计算结果的精度. * * * * * * * * * * * * * * * (2-15) (2-14) 磁场的相应FDTD公式为 (2-17) (2-16) 其中 (2-18) 通常采用迭代的方式来利用上述FDTD差分方程计算电磁场 若已知t1=t0=nΔt时刻空间各处E 的值 计算t2=t1+Δt/2时刻空间各处H 的值