论文部分内容阅读
随着科技的不断发展,电磁波的应用越来越广泛,对于超电大目标的电磁散射特性分析一直是计算电磁学研究的热点。特别是在军事领域,对于飞机、航母这类超电大复杂的目标的电磁散射分析,是将来信息化战争的研究重点。积分方程是分析电磁散射问题的强力手段,但用传统积分方程法求解问题时,会遇到迭代速度慢、迭代不收敛的问题。多层快速多极子是一种快速算法,能够将矩阵矢量相乘的时间复杂度降低到O(NlogN)。区域分解算法能够有效地解决多尺度问题。但仅通过这两种手段,对于超电大目标的电磁散射分析问题,也是非常棘手困难的。对于这一问题,本文开展了对并行的区域分解算法的研究,研究了基于并行多层快速多极子算法的并行区域分解算法。首先,本文基于等效原理推导了表面积分方程,并介绍了矩量法求解积分方程的重要步骤:网格离散、基函数选取、测试方法和矩阵求解。重点介绍了基于Krylov子空间的GMRES迭代方法,详细阐述了多层快速多极子算法的基本原理。简单介绍了基于共享内存式并行的OpenMP和基于分布式内存并行的MPI。接着,详细地推导了不连续伽辽金积分方程,先介绍了对于非共形网格常使用的mono-polar RWG基函数,然后基于对偶配对测试准则,建立了两个内罚项,得到了使用mono-polar RWG基函数的不连续伽辽金积分方程,并用数值算例验证了其精度。最后,针对超电大目标的电磁散射特性分析,介绍了区域分解的基本原理和基于OpenMP和MPI混合并行多层快速多极子的并行策略——高层采用盒子并行,低层采用角谱并行。区域分解算法的内迭代和外迭代都使用多层快速多极子算法进行加速,但由于两个迭代处理的数据不相同,所以不能简单地实现并行。本文提出了一种沿用多层快速多极子并行算法数据结构的区域分解并行算法,并用该方法完成了对未知量为2700万飞机RCS的计算。