多核间接和并行迭代法及其在辐射流体力学中的应用.pdf

摘 要 大规模并行计算机的快速发展和应用，使得复杂物理系统的高分辨率数值模拟已成 为可能．在这些数值模拟中，系统隐式离散后，通常需要求解稀疏线性代数方程组，所耗 费的时问有的甚至达到模拟时间的80％．为了有效缩短数值模拟时间，要求相应的应用 程序必须具有高度的可扩展能力．因此．在大规模并行机上发展具有高可扩展能力的线 性解法器，是实现高性能数值模拟的重要组成部分，具有重要的理论意义和应用价值． 本文的工作主要有以下五个方面： 第一，通过减少(串行)GPBiCR方法中全局同步化点的个数，提出了并行GPBiCR方 稳定性，我们对两种方法在分布式存储并行机上实现时的并行性能，进行了分析，得 到新方法的同步开销减少为原米的三分之一．对比GPBiCR方法，从理论上证明 了PGPBiCR方法的可扩展性得到了提高，通讯性能提高比率趋向于66．7％．数值试验得 出了与理论分析相吻合的结果． 第二，通过消除(串jT)BiCGSTAB2方法中内积计算的数据相关性，提出了并 行BiCGs 法具有相I吲的数值稳定性和收敛性．为了保证性能分析的可靠性与公平性，把连续的偶 数步与奇数步两次迭代整体看成一个大迭代步．对L匕BiCGSTAB2方法，从理论上证明 了PBiCGS’IIAB2方法的可扩展性得到了提高，通讯性能提高比率趋向于66．7％．对两种 方法在分布式并行机上进行了性能测试，所得数值试验结果与理论分析结果相吻合．同 时，也通过数值例子验证了两种方法具有相同的收敛性和数值稳定性． 第三，基于目前占主流的多核体系结构的并行计算机，以及利用矩阵的低秩 structuredmultifrontal 性所产生的快速健壮的结构化多波前分解(robust factoriza- 对RSMF方法的嵌套分割排序、符号分解以及数值分解三部分进行了多核并行实现．其 中符号分解和数值分解是基于二叉消元树的分支结构来实现多核并行的．数值试验表 明MRSMF方法是有效的． 第四，基于-维三温辐射扩散方程及其离散系统，利用离散所得的三温矩阵所具有 的块对角结构，对RSMF方法中的重排序和符号分解部分进行了分块操作．在数值分解 部分利用了块矩阵操作技术．使得重排序和符号分解部分的时间效率得到了明显改善． 同时基于良好的排序和块矩阵操作技术，数值分解部分的时间效率也得到了提升．我们 万方数据 多核直接和并行迭代法及其在辐射流体力学·}·的心用 把这种分块形式的方法称作分块RSMF(记作BRSMF)方法．数值试验部分通过五点离 散模板所得的三温矩阵验证了BRSMF方法的稳定性与高效性，而且随着问题规模的增 大，改进效果越来越显著． 第五，基于前面所提出的MRSMF方法和BRSMF方法，进一步给出了适用于三温线 的特点与BRSMF方法分块操作的特点充分结合了起来，使得所产生的MBRSMF方法同 时具有了多核并行和分块操作的特性．数值试验部分基于二维三温辐射扩散方程的五 点离散模板，验证了MBRSMF方法L匕MRSMF和BRSMF方法更有效． 最后我们利用具有实际意义的三温模型所产生的线性方程组，对本文所给 测试结果进一步验证了所提方法的有效性． 关键词：稀疏线性代数方程组，K巧lov子空间迭代法，多核体系结构，直接法，二维三温 辐射扩散方程． 万方数据 Abstract The and of havemake rapiddevelopmentapplicationmassivelyparallelcomputers it forthe resolutionnumericalsimulationofcomplexphysicalsystems．