热传导方程的一种新的具有并行本性的区域分解算法及其并行程序设计.pdf

热传导方程的一种新的具有并行本性的 区域分解算法及其并行程序设计1 左风丽*1 袁光伟*2 北京应用物理与计算数学研究所 应用科学计算部 高性能计算中心，北京 8009 信箱 10088 zfl_310@ ABSTRACT A domain decomposition method based on Du Fort-Frankel scheme and the fully implicit scheme for heat conduction is presented and some characters are proved by numerical experiments in this paper. Two kinds of parallel strategies are designed and some programming techniques are discussed, which is less memory storage than conventional programming method leading to obtaining good result. By performance analysis of parallel numerical experiments, the method with intrinsic parallelism has the same accurate as the fully implicit scheme, which suit large scale scientific computing and engineering computing. KEYWORDS Heat conduction; Domain decomposition method; High accuracy; High parallel efficiency 1. 引言 许多物理和力学问题的求解都可归结为热传导方程的求解，其数值模拟常采用有限差分 方法，在实际应用问题中通常采用隐式差分格式，这是因为隐式差分格式具有无条件稳定的 特征，对任意的网格比都能保证对计算误差的控制。随着超大规模并行计算机的发展，对所 要求解的问题的精度要求也越来越高，因而求解问题的规模变得越来越大。尽管单机的计算 速度和内存容量变得很大，但是仍难以完全满足需求，其有效的解决方法就是构造良好的具 有并行本性的区域分解方法。本文讨论求解一维热传导方程的区域分解方法，将求解的区域 剖分成多个子区域，在子区域的内边界点使用 Du Fort－Frankel 格式，在子区域的内点使 用全隐式格式。它不但具有并行本性的特征，而且具有与纯全隐式格式一样的精度和稳定性。 在构造并行算法的具体实现时，提出一种动态局部数据管理的方法，这样做的好处不仅可以 提高数据的局部性，而且可以任意改变内存的大小，而不改变通信的次数和通信量的大小， 因此可以获得良好的结果。 我们知道，对热传导方程的数值求解已有很多的方法，通常采用的格式为隐式格式。隐 式格式之所以被广泛的应用是因为它是无条件稳定和绝对收敛的，但是，求解隐式格式需要 1 本课题得到国家 973 项目（G1999032801），国家自然科学基金和中物院重点基金资助。 *1 左风丽，硕士，助理研究员，大规模科学与工程计算并行研究方向。 *2 袁光伟，博士，研究员，大规模科学与工程计算并行研究方向。 1 解一个代数方程组，并且它所求解的方程组的阶数与划分节点个数是一样的，这样，随着计 算节点的增加计算量变得越来越大，甚至难以计算下去。由于对求解问题的要求越来越高， 网格节点不断加密，即使有高性能大存储量的并行机，存储空间也难以完全满足需求。而基 于MPI编程的区域分解技术是减少单机上计算量和存储量的一条有效途径。例如，周毓鳞[1,2]、 沈隆钧、袁光伟[3]、张宝琳 [4