一种孔隙介质中地下水流并行计算方法

计算机工程与应用 ›› 2011, Vol. 47 ›› Issue (20): 234-237.

一种孔隙介质中地下水流并行计算方法

程汤培1，王群1，季晓慧1，黄林显2，徐腾2

1.中国地质大学（北京）信息工程学院，北京 100083
2.中国地质大学（北京）水资源与环境学院北京 100083

收稿日期:1900-01-01 修回日期:1900-01-01 出版日期:2011-07-11 发布日期:2011-07-11

Parallel computing method for groundwater flow simulation in porous media

CHENG Tangpei1，WANG Qun1，JI Xiaohui1，HUANG Linxian2，XU Teng2

1.School of Information Engineering，China University of Geosciences（Beijing），Beijing 100083，China
2.School of Water Resources and Environmental Science，China University of Geosciences（Beijing），Beijing 100083，China

Received:1900-01-01 Revised:1900-01-01 Online:2011-07-11 Published:2011-07-11

摘要/Abstract

摘要： 针对孔隙介质中地下水流动问题提出了一种并行数值计算方法，并基于此设计了一套专用于求解大规模三维地下水流动方程的并行计算模块。计算模块基于区域分解的方法实现对模型区域的并行求解，采用了分布式内存和压缩矩阵技术解决大规模稀疏矩阵的存储及其计算，整合多种并行Krylov子空间方法和预条件子技术迭代求解大规模线性方程组。在Linux集群系统上进行了数值模拟实验，性能测试结果表明，程序具有良好的加速比和可扩展性。

关键词: 地下水数值模拟, 并行迭代算法, 有限差分法, Krylov子空间方法, 预条件子

Abstract: This paper presents an efficient parallel computing method for large-scale groundwater flow simulation in porous media，and it is implemented in a prototype tool.Domain decomposition method is adopted to partition the modeling area.Distributed memory and compressed matrices technology are introduced for both the storage and evaluation of the large-scale sparse matrices.All kinds of Krylov subspace methods and preconditioners are integrated to assemble and solve linear systems of equations.Experiments are carried out on Linux clusters，and the results demonstrate that this method can achieve desirable speedup and scalability.

Key words: groundwater simulation, parallel iterative algorithm, finite difference method, Krylov subspace methods, preconditioner

程汤培1，王群1，季晓慧1，黄林显2，徐腾2. 一种孔隙介质中地下水流并行计算方法[J]. 计算机工程与应用, 2011, 47(20): 234-237.

CHENG Tangpei1，WANG Qun1，JI Xiaohui1，HUANG Linxian2，XU Teng2. Parallel computing method for groundwater flow simulation in porous media[J]. Computer Engineering and Applications, 2011, 47(20): 234-237.

[1]	袁冬芳，曹富军. 三维对流扩散方程的稀疏存储及预条件迭代[J]. 计算机工程与应用, 2018, 54(4): 56-59.
[2]	孙蕾. 求解大型非对称稀疏线性方程组的FIMinpert算法[J]. 计算机工程与应用, 2016, 52(21): 63-67.
[3]	徐康丽1，杨志霞1，蒋耀林2. 基于Krylov-Schur重启技术的Arnoldi模型降阶方法[J]. 计算机工程与应用, 2016, 52(12): 251-255.
[4]	渐令1，张建松1，宋允全1，赵敏2. 应用有限差分法模拟高炉炉缸侵蚀[J]. 计算机工程与应用, 2012, 48(9): 27-29.
[5]	金涛斌，邹军. 基于CN-FDTD的新颖DGS带阻滤波器特性研究[J]. 计算机工程与应用, 2012, 48(8): 1-3.
[6]	戴宽平，吕全义. 改进的并行Arnoldi方法[J]. 计算机工程与应用, 2011, 47(34): 72-73.
[7]	金涛斌1，2，金杰1，李媛1，李可佳1. 可调谐矩形环DGS的滤波特性研究[J]. 计算机工程与应用, 2011, 47(25): 143-145.
[8]	周宏宇，王爱民. 非线性偏微分方程数值求解的自适应方法研究[J]. 计算机工程与应用, 2011, 47(20): 38-40.