06_2,FEKO电磁仿真数值方法.ppt

二维多层快速多极子方法中的分层级结构 基于上图的分层级结构，多层快速多极子方法由上行过程、下行过程两部分组成。上行过程分为最高层的多极展开(Mexp)、子层到父层的多极聚合(MM)。上行过程在多极聚合到第二层后，经远亲转移计算转向下行过程。下行过程则分为父层到子层的多极配置(LL)、同层间远亲组的转移(ML)和最高层的部分场展开(Lexp)。 最初层 次细层 最细层 多层快速多极子 对于N体互耦，多层快速多极子方法采用多层分区计算。即对于附近区强耦合量直接计算；对于非附近区耦合量则用多层快速多极子方法实现。多层快速多极子方法基于树形结构计算，其特点是：逐层聚合、逐层转移，逐层配置、逐层递推。由于每层数值计算量均为 量级，共有 层，所以多层快速多极子方法计算矩阵与矢量相乘的工作量为 量级。内存需求也为 量级。 次相邻中心的概念 多层快速多极子的复杂度 和快速多级子一样多层快速多级子技术也是将矩阵矢量乘分解为聚集，转移发散三个步骤，然而聚集和发散的过程是从最高层进行的，后通过移置中心，以及插值来完成低层的聚集和发散，至于转移过程，则和FMM方法是完全一样的，只是多层快速多级子只在同一层的次相邻中心进行。次相邻中心：是指在本层不属于相邻中心，但在第一层隶属的组是相邻组的那些中心。 计算流程 网格剖分与几何信息读入 最细层的近场信息 聚集 (内插值) 近场矩阵 近场作用 转移 (次相邻组) 矩量法 迭代法求解 向量运算 矩阵-向量乘积 模型分组 发散 (外插值) 各层的远场信息 远场作用 预条件器 线性方程组 电磁场积分方程 多种算法的比较 精度： 矩量法精度最高，有限元次之，时域有限差分最差 矩量法没有数值色散误差，其它两种都有 时域有限差分不仅有数值色散误差，且模拟复杂几何形状的误差一般也要大于其它两种算法 实施性 矩量法最难，有限元次之，时域有限差分最易 MOM要注意基函数的选择，处理奇异点，构思快速求解技术 FEM只需注意基函数选取及稀疏矩阵存储方式 －－摘自　盛新庆　《计算电磁学要论》 FEKO计算流程及对应算法步骤说明 计算流程—FEKO案例 案例：飞机RCS分析 几何信息读入检查 模型分组 近场作用 向量运算 预条件器 远场作用 矩阵矢量乘与迭代求解 计算RCS * 说明 在这三部分内容上，相信咱们已有了较深的体会，因为咱们已经使用过试用的版本 说明在后面的例子中，也都通过实例向大家进行展示 作为一款成熟的商业软件，FEKO的后处理能力不容质疑。 各种参数的直接显示、坐标及图形的直接修改编辑等能满足绝大多数用户需要。 此外，所有的计算数据完全开放，以文本的形式存放，每个参数存放的位置固定。这一特征让用户可以轻松从Matlab、VC等程序中读写调用计算数据，便于二次开发。 * 强调 cem软件 非常重要的是内核，是求解的算法，而feko在这一点上也是非常突出的 * 有限元法必须在分析模型的外部添加虚拟截断边界，对于多个物体还必须分析目标之间的传播空间，导致计算量很大。 矩量法求解的是积分麦克斯韦方程，有限元求解的是麦克斯韦微分方程，虽然方程形式不同，但是利用MoM/FEM交界区域的切向场连续性，可以把两种算法结合起来，我们知道，矩量法自动满足辐射条件，无需设定外部截断边界就可以进行场的计算。如果把二者结合起来，就可以避免分析无效传播空间和在分析模型的FEM区域外部添加截断边界造成的大量未知量，从而大大节约计算资源，提高计算效率。 * 快速多级子技术始于20世纪80年代，在90年代得到大力发展，目前学术界研究最热的就是快速多级子技术。 FEKO融合了最新的已成熟的理论成果。 快速多级子技术不会导致精度下降。 对这个例子，矩量法所需的内存在300G左右 * 对于电大尺寸的问题，局部效应起主要作用，从理论上来说，高频方法能够得到满意的精度。 物理光学法极大地降低了计算资源需求，使得可求解地未知量达到1000万量级。 一致性几何绕射理论的计算资源需求不再受结构尺寸的限制，大大提高了可计算的问题规模。 * 对于电大尺寸问题，有时即使采用最高效的算法，计算资源需求和计算时间也是相当大的。因此并行处理就势在必行。 FEKO从创建初就非常重视并行处理的实现，如今FEKO的并行技术已经非常成熟，在整个计算过程（包括矩阵运算、线性方程组求解等）中实现了真正意义上的并行（不是其它软件的“伪”并行）。 左边的图表统计了矩量法的并行处理效率，可以看出，并行处理的确可以提高计算的效率，并行效率约为60%。 右图是以基站天线为例来分析Feko软件UTD算法的并行效率问题，右下图给出了实际的并行加速比与理想的并行加速比关系。 可见，Feko软件的并行效率很