第1课FDTD基础要点解读.ppt

Unit 1： 时域有限差分方法 FDTD简介 Mathematic Basis 差分近似 FDTD 直观 存储空间仅仅与网格总数成正比 同时获取时域信息和频域信息 可一次获得宽频带的信息 FDTD Yee cell (leapfrog) Yee网格的特点 Yee网格体系的特点是，E和H各分量在空间的取值点被交叉放置，使得在每个坐标平面上每个E分量由四个H分量环绕，同时每个H分量由四个E分量环绕。这样的电磁场空间分布符合电磁场的基本规律――法拉第电磁感应定律和安培环路定律的自然结构，即符合麦克斯韦方程的基本要求，能够恰当地描述电磁场的传播特性。 为实现空间坐标的差分计算，并考虑导电磁场在空间相互正交和交链的关系，在Yee网格中，每个坐标轴方向上场分量间相距半个网格空间步长，因而同一种场分量之间相隔一个空间步长：电场和磁场在时间顺序上相隔半个时间步长，这样，使麦克斯韦方程离散后可以构成显式差分方程，从而给出电磁问题的初始值后，可以在时间上迭代求解。 3-D Maxwell’s equations 直角坐标系中的FDTD方程 六个标量场方程 E公式推导举例 H公式推导举例 E update equations H update equations FDTD在时域的逐步推进计算 空间步长 麦克斯韦方程建立的有限差分算法方程是用对空间和时间的差分代替微分，这种替代将产生误差，能使在计算网格空间所模拟的波形产生色散。换句话说，在FDTD网格空间中存在的数值模，其相速取决于模的波长、传播方向以及网格单元的尺寸。数值相速随传播角和网格尺寸而变化的慢波效应是FDTD算法的一个固有属性。要保证算法的精度，空间步长 、 和 要取得足够小，通常取所感兴趣的最高频率所对应波长的1/10到1/20。 解的稳定性 在执行FDTD算法时，随着时间步长的增长，保证算法的稳定性是一个很重要的问题。数值解是否稳定主要取决于时间步长与空间步长间的关系。对于一般的三维FDTD网格，数值解稳定的条件（CFL条件）是： 调制高斯脉冲波形 微带线加源的方法 微带线加源的方法 时间迭代步数 散射参量的计算 散射参量的计算 入射电压 观察点1处的电压 观察点2处的电压 散射参量 观察点电压图 散射参量图 高斯脉冲源频域波形 调制高斯脉冲 调制高斯脉冲 如果要使直流分量接近为0，有效频谱中心位于 ，我们就可采用调制高斯脉冲源 调制高斯脉冲频域波形 加源的方法 有了源的形式后，源的加法要根据特定的EM结构来合适的选取。对于一般的平面微波电路和天线问题，就是要在电路或者天线的输入端口加源，使之沿馈线传播想要得到的模式。这就需要在仿真问题之前，了解其在馈线上的横向场分布。 如果被研究的结构在离源点一定距离的某处，由激励源产生的数值波将最终传播到达该结构，经过波与结构的相互作用，发生传输与反射，直到瞬态过程消失。这意味着传输波和反射波的时间过程均已完成，波已离开仿真区域，对该时间过程作傅里叶变换，可一次性得到传输波和反射波的宽频带振幅与相位信息。在稳态激励源作用下，电磁场达到稳定的时间步数与许多因素有关，电磁结构的复杂程度和性质起着非常重要的作用。对一般的结构，电磁波在其内至少传播两个来回，才能稳定。 任意两端口微带线结构 fft fft fft * * Nanjing University of Science Technology 计算电磁学 樊振宏 研究背景 近年来，微波和毫米波通信、航空航天、雷达、精确制导等民用军用系统朝着高频化、微型化、多功能、高可靠性以及低成本方向发展。 由于色散、不连续性和封装而产生的失真与时延，以及由于耦合而产生的串音噪声等问题变得十分严重，传统的准静设计方法已不能满足设计要求，必须采用精确的电磁场全波分析方法。 电磁分析的本质是求解Maxwell Equation 在特定初始条件特定边界条件下 边界的复杂，导致传统的解析分析方法无法胜任 计算电磁学的应用领域? 电磁场全波分析方法分类 基于微分方程模型的分析方法 时域有限差分 FDTD – Finite Difference Time Domain 频域有限差分 FDFD – Finite Difference Frequency Domain 有限元 FEM – Finite Element Method 基于积分方程模型的分析方法 矩量法 MoM – Method of Moment 基于矩量法的快速算法 3-D Maxwell’s Equations (Faraday’s law) (Ampere’s law) 各向同性介质中的本构关系为： 时域有限差分法 （Finite‐Dif