快速多极边界元并行算法的研究与工程应用

论文摘要

边界元法作为有限元法、有限差分法等区域解法的重要补充,具有降维、精度高的特点,在各种工程领域有许多成功的应用。但传统边界元法形成的线性方程组的系数矩阵是满阵,因此一般只能应用于数千到数万自由度的小至中等规模问题的求解。然而,近20年以来,由于一些快速算法的出现并不断发展,情况正在发生改变。特别是由Greengard和Rokhlin首先提出的快速多极算法,可以将矩阵向量相乘操作的计算量级和存储量级同时降为O ( N ),其中N是未知量的个数,被美国工业与应用数学学会(SIAM)评为上个世纪十大算法之一。应用快速多极算法的边界元法被称为快速多极边界元法。而在此基础上对其可扩展性并行计算算法的研究则可以进一步扩大解题规模、提高计算速度、并通过网格加密和增加级数展开阶数两方面来提高求解的精度、从而进一步扩展边界元法的应用范围和优势领域。以此为目标,本文对适用于二维、三维弹性力学问题等的初始版本和新版本快速多极边界元法的并行算法作了研究,结合边界元二次等参单元的特点,提出了基于自适应树结构的快速多极边界元法的一种分布式并行计算格式,其中包括计算量预测、加权任务分配方式、通信关系的建立和数据通信过程等几部分的实现,并对几何形状不规则的结构在最多为64个处理器的并行机群上进行了测试并达到了满意的加速比,最大计算规模达到了二百多万自由度。在此基础上,本文通过数值算例对基于树结构的快速多极并行计算与常规矩阵并行计算进行了比较,指出了它们各自在计算规模上的合理应用范围。并给出了一种结合算法,从而在预处理方案、多子域问题等方面扩展了快速多极边界元并行算法的应用范围。最后,作为上述算法的具体应用,本文对短纤维复合材料进行了数值模拟,利用高性能计算的优势,对复杂形状纤维以及随机取向弯曲纤维复合材料中的应力分布规律进行了研究,给出了一些有参考价值的数值结果,并表明快速多极边界元法并行计算在复杂界面问题上的大规模计算与常用的有限元等区域方法相比具有明显的优势。

论文目录

• 摘要
• Abstract
• 第1章 绪论
• 1.1 边界元法及其快速多极算法概论
• 1.2 并行计算技术概论
• 1.3 并行计算在边界元法及快速多极算法中的实现概述
• 1.4 本文的主要内容
• 第2章 快速多极边界元法概述
• 2.1 引言
• 2.2 边界元法基本格式
• 2.3 边界元数值方法以及迭代求解算法
• 2.4 快速多极边界元法基本格式
• 2.5 快速多极边界元法数值实现步骤
• 2.6 数值算例
• 2.7 本章小结
• 第3章 快速多极边界元并行计算算法
• 3.1 引言
• 3.2 快速多极边界元并行计算流程图
• 3.3 不平衡树结构的一维映射
• 3.4 任务划分
• 3.5 通信关系的建立和数据通信过程
• 3.6 数值验证算例
• 3.7 本章小结
• 第4章 快速多极与矩阵分布式并行计算的比较与结合
• 4.1 引言
• 4.2 多子域问题中的矩阵向量相乘格式
• 4.3 快速多极算法与矩阵相结合的并行求解算法
• 4.4 数值算例
• 4.5 本章小结
• 第5章 并行快速多极边界元法在纤维增强复合材料中的应用
• 5.1 引言
• 5.2 快速多极边界元重复相似子域法并行实现基本格式
• 5.3 数值应用算例
• 5.4 本章小结
• 结论
• 参考文献
• 致谢
• 个人简历、在学期间发表的学术论文与研究成果