电磁散射问题退化核快速解法的研究

摘要

采用数值方法可以灵活地对任意目标建模仿真求解，已成为获得目标散射特征的重要方法，是计算电磁学的重要研究领域。在各种数值方法中，矩量法由精度高，未知量规模小而得到了广泛应用，但是求解矩量法得到的线性方程组需要大量的存储空间和计算时间。为了克服上述问题发展出了一系列快速方法，其中远区作用矩阵块的高效表示是构建快速算法的重点和难点。退化核技术可以有效用于矩阵的低秩近似表示。
　　本文以层次型矩阵（H矩阵）框架为基础，深入研究了格林函数的快速退化核构造方法，用来构建快速算法求解电磁散射问题。主要研究内容和贡献主要分为以下几个方面：
　　1.提出一种线性计算复杂度的公共边生成方法，用于采用RWG基函数矩量法通用求解方法的前处理中。该算法基本思路为只通过比较邻近三角形面元的顶点编号来得到公共边信息，从而将计算复杂度由O(N2)减少为O(N)，这里N为三角形面元数目。该算法可以集成到通用的仿真软件和算法测试平台中。此外，该算法已成功推广应用于快速生成四面体体元剖分的公共面。
　　2.研究了采用Lagrange多项式和径向基函数（RBF）构造的插值退化核。格林函数为一个震荡函数，基于距离的径向基函数插值的效率更高，只需要更少的插值点就能得到较高的插值精度。针对导体散射问题，提出一种近表面插值网格来减少插值点的数目。相比传统的均匀插值网格，近表面插值网格只选用靠近目标体表面的点，从而可以将插值插值点的数目由O(k3)下降为O(k2)，其中k为每一维上的插值点数目。计算效率大幅提高，同时保持精度基本不变。
　　3.基于分层基层次型矩阵（H2矩阵）框架，提出了一种用于求解导体电磁散射问题的混合插值网格的快速算法。该方法采用高斯径向基函数构造插值退化核函数，根据不同层划分得到的立方体块大小：对于较小的块，插值网格选用均匀网格，其它块选用近表面插值网格。该方法结合了两种网格的优点。相比于只采用均匀网格的方法，所提方法的计算量和存储量都有所下降。
　　4.提出了一种新颖的指数退化核函数。指数退化核函数可以将格林函数源点与场点作用分离，得到的转移矩阵是对角矩阵，相比于多极子展开，达到相同精度时所需展开项数减少。从而可以有效地加快迭代求解算法中的矩阵向量乘法运算。在此基础上，提出一种新的方向选择方案，可以进一步减少展开项数至多极子展开的一半，且随着块尺寸变大比多极子展开增长更慢。数值算例验证了算法的有效性。

目录

封面

声明

中文摘要

英文摘要

目录

专用术语注释表

第一章 绪论

1.1课题研究背景与意义

1.2矩量法快速算法的现状

1.3本文的主要内容与创新点

第二章 矩量法及其前处理

2.1引言

2.2矩量法

2.3公共边的快速生成方法

2.4公共面和非公共面的快速生成方法

2.5本章小结

第三章 层次型矩阵及方法

3.1引言

3.2层次型矩阵

3.3分层基层次型矩阵

3.4核函数的压缩近似

3.5本章小结

第四章 插值退化核函数

4.1引言

4.2 Lagrange插值退化核

4.3径向基函数插值退化核

4.4基于层次型矩阵的快速算法

4.5基于混合插值网格的快速算法

4.6本章小结

第五章 指数退化核函数

5.1引言

5.2快速多极子方法回顾

5.3指数退化核函数

5.4数值算例

5.5退化核函数比较

5.6本章小结

第六章 总结与展望

6.1本文总结

6.2展望

参考文献

附录1 攻读博士学位期间撰写的论文

附录2 攻读博士学位期间申请的专利

附录3 攻读博士学位期间参加的科研项目

致谢