多层快速多极子加速区域分解迭代算法.pdfVIP

第22卷增刊 电波科学学报 、，01．22，Sup CHDmSE oFRADIoSCⅡNCE 20cr7年9月 JoURNAL Septl锄b％2007 多层快速多极子加速的区域分解迭代算法 陈 明 叶宗民陈如山 丁大志 樊振宏 (南京理工大学通信工程系，南京210094) 摘要本文提出的基于多层快速多极子的区域分解迭代算法，将目标结构表面分 成若干子区域，每个子区用其附近的缓冲区扩展成迭代区，用积分方程对扩展后的 区域求解，求解后舍弃缓冲区的电流，保留原子区域的电流并更新到整个区域的电 流，只需要两三次外迭代即可获得准确的电流和散射场，大大节省了内存的消耗， 大幅提高了对电大尺寸目标的分析能力。并将算法进行了并行，更进一步的极大提 高了分析超电大尺寸目标的能力。 关键字多层快速多极子Q缸n嗄M)，区域分解(DDM)，并行算法 法相比较，重叠型区域分解迭代算法能够更有效地 1．引言 抑制子区域边界的电流奇异性。 如何准确快速的对电大尺寸复杂结构的电磁 特性进行仿真分析一直是计算电磁学研究领域持 ； 。卜——呻’ 1． ．i 1 3 4{ M 久不衰的热点和难点问题。矩量法[1】具有求解区 i2 {5 — l一 ；一 一i 一 - 一 I一 域仅仅限制在源和等效源区的优点，精度较高，对 !一 。； !’ 电大问题而言，存储量和计算量过大，因此众多学 圈l区域分解示意图 、 者致力于扩大矩量法适用范围的研究。九十年代以 来，伴随着各种快速算法的出现和计算机技术的飞 现给出子区f的缓冲区易(f)的概念：在整个去 速发展，计算电磁学取得了长足的进步。其中的快 内不属于f子区而与f子区边界相邻的带形区域。 更一般的结构如图4-l所示，f子区的边界和虚线 速多极子方法(FMM：蠡lstmul卸olememod)【3】和多 fast 所围的带形区域为f子区的缓冲区。我们用QM表 层快速多极子方法(M【n压M：multilevel 示子区珀勺缓冲区，Q，‘=Q；+Q珊表示扩展后的求 multipolemethod)[2】就是当今积分方程研究成果 的杰出代表。目前快速多极子方法和多层快速多极 解子区域，Q；=Q—Q；’是扩展子区域Q；’的补余区 域。 子方法已经广泛应用于各种复杂目标的电磁辐射 与散射分析。 根据上述原理，可建立子区域之间积分迭代公 式，定义两个关于电流．，的线性算子【5】，对于磁 BreImaIl等在处理三维电场积分方程的前向后 向迭代时，为了抑制电流的奇异性，提出了前向后 场方程 向延伸迭代算法(FBBS)。最近，又有重叠型的区