天津大学金杰电磁场与电磁波电磁算法及仿真.ppt

第四章 电磁算法及仿真 主要内容 电磁场数值算法 MOM、FEM、FDTD、MRTD 电磁仿真软件 Maxwell、CST、HFSS 电磁场数值算法 加权余量法 边界积分法 边界元法 矩量法 小波基 快速算法 时域多分辨分析法 内域积分法（伽略金法） 有限元法 麦克斯韦方程组 时域积分方程法 有限差分法 时域有限差分法 矩量法 MOM (Method of Moment) 概况 R.F.Harrington在20世纪60年代对矩量法求解电磁问题做了全面深入分析。 核心思想：根据线性空间理论，N个线性方程的联立方程组、微分方程、差分方程及积分方程均属于希尔伯特空间中的算子方程，它们可化作矩阵方程予以求解，在求解过程中需计算广义矩量。 基于MOM的电磁场计算软件：ADS的momentum，Sonnet等 特点 频域矩量法比较成熟，时域矩量法有待发展； 矩阵规模的大小涉及到占用内存的多少，在很大程度上影响了计算的速度。 理论 在静电学中，在由点的电荷分布在点产生的电位分布可以表示为 设 的一个解是 将解带入电位函数中可得 即 其中 亦即 或 有限元法 FEM (Finite Element Method ) 概况 起源于土木工程和航空工程中的弹性和结构分析问题的研究，它的发展可以追溯到Alexander Hrennikoff(1941)和Richard Courant(1942)的工作。 核心思想：由解给定的泊松方程化为求解泛函的极值问题。 基于MOM的电磁场计算软件：HFSS、ANSYS。 原理 将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。 特点 近似性仅限于相对小的子域； 将函数定义在简单几何形状的单元域上，不考虑整个定义域的复杂边界条件。 时域有限差分法 FDTD (Finite-Difference Time Domain) 概况 1966年K.S.Yee提出的； 核心思想：把带时间变量的Maxwell旋度方程转化为差分形式，模拟出电子脉冲和理想导体作用的时域响应； 是目前计算电磁学界最受关注，最时髦的算法，但还在发展完善之中 基于FDTD的电磁场计算软件：XFDTD等 关键要素 差分格式 解的稳定性 吸收边界条件 特点 广泛的应用性 节约运算和存储空间 适合并行计算 计算程序的通用性 简单直观，容易掌握 计算步骤 采用一定的网格划分方式离散化场域； 对场内的偏微分方程及各种边界条件进行差分离散化处理，建立差分格式，得到差分方程组； 结合选定的代数方程组的解法，编制程序，求边值问题的数值解。 原理 麦克斯韦方程组 本构关系式 Yee元胞 Yee采用矩形网格进行空间离散，将每个节点进行编号，节点的编号和其空间坐标位置按照下面的方式对应起来。 任意函数F(x, y, z, t)在时刻nΔt的值可以表示为 计算方法 若已知t1=t0=nt时刻空间各处E的值 计算t2=t1+t/2时刻空间各处H的值 计算t1=t2+nt/2时刻空间各处E的值 差分方程 其他节点 可同样求得。 关键技术 数值稳定性（解的收敛性） Courant稳定性条件 数值色散性 相速与网格尺寸和传播角的关系： 通常要求： 吸收边界条件 目前吸收边界条件研究方向主要有： 1. 通过波动方程的因子分解获得单行波方程并取近似来建立吸收边界条件，Mur吸收边界条件和Liao吸收边界条件是两个典型形式； 2. 在边界上引入吸收材料，电磁波在无反射地进入吸收材料后被迅速衰减掉，如Berenger提出的完全匹配层吸收边界条件。 激励源 分析电磁问题的一个重要任务是对激励源的模拟，即选择合适的入射波形式以及用适当方法将入射波加入到FDTD迭代中，也就是说，应将被研究媒质在真实源激励下这一完整条件在数值计算中尽可能近似地复现出来。 远—近场外推 吸收边界 散射场区 散射场区 总场区 目标 输出边界 连接边界 辐射场区 辐射场区 激励源 吸收边界 输出边界 散射计算区域划分 辐射计算区域划分 概况 1996年美国学者Krumpholz和Katehi提出的； 核心思想：将麦克斯韦微分方程按照电场与磁场在时间和空间分别正交的 方式用小波展开，应用伽略金法，得到了按照时间和空间迭代的差分方程。 时域多分辨分析法 MRTD (Multi-Resolution Time Domain) 发展趋势 提高计算精度的同时，保持算法的稳定性。 加强关于边界问题的处理和介质的处理。 增强计算电大尺寸对象的能力