**标题:共轭梯度法在三维随机异质多孔介质中的地下水流问题中的数值研究**
**PHILIP D. MEYER 和 ALBERT J. VALOCCHI**
*伊利诺伊大学香槟分校土木工程系*
**STEVEN F. ASHBY**
*劳伦斯利弗莫尔国家实验室,利弗莫尔,加利福尼亚*
**PAUL E. TAYLOR**
*伊利诺伊大学香槟分校计算机科学系*
**摘要:**
本文考虑在三维随机异质多孔介质中对地下水流进行场地规模建模。由该问题的有限差分离散化得到的线性方程组可能涉及超过一百万个未知数。这种规模的问题只能在超级计算机上实际求解,并需要一种对所使用的特定计算机体系结构非常适用的高效迭代求解方法。经过预处理的共轭梯度法对于这个地下水流问题非常高效,通过使用适当的预处理矩阵,该方法可以适应不同的计算机。数值软件包CgCode增强了这种适应性;CgCode的用户只需提供用于预处理和矩阵-向量乘法的子例程。这些子例程应根据具体问题和机器体系结构进行定制。文章呈现了在矢量机(Cray X-MP/48)和矢量-并行机(Alliant FX/8)上使用多项式预处理来提高共轭梯度法效率的数值结果。
**引言:**
使用随机过程理论对含水层的异质性进行建模在地下水科学家和工程师中越来越被接受,作为研究场地尺度溶质运移期间显现的宏观弥散现象的手段。在场地尺度上进行建模不可避免地涉及大量的计算工作,其规模的大小来自两个条件。首先,宏观弥散是一种大尺度现象。在随机过程的术语中,宏观弥散发生在许多过程的相关长度上。因此,建模域必然很大。引起大量计算工作的第二个条件是必须对随机过程的相关结构进行建模的细尺度。域的离散化可能涉及在相关长度中的若干网格点。
本研究涉及使用七点有限差分逼近模型的三维稳态饱和地下水流建模。水力导数场是呈对数正态分布的随机场的单个实现,具有指数协方差。上述两个条件,加上问题的三个维度,可能导致一个具有一百万个或更多网格点的域。求解由此产生的线性方程组需要一个高效的迭代方法;由于问题规模庞大,直接方法是不适用的。本文描述了共轭梯度法,这是一种对于这里考虑的问题类型非常有优势的迭代技术。虽然共轭梯度法不需要像其他方法那样对迭代参数进行先验估计,但它的有效性取决于预处理器和所采用的计算机体系结构。讨论了几种预处理器,并探讨了预处理器与其在具有串行、矢量和矢量-并行体系结构的计算机上的性能之间的关系。给出了地下水流问题的数值结果。
这些数值实验的目标不是确定地下水流问题的最佳可能迭代方法。共轭梯度法已被证明与其他迭代方法(如逐步超松弛法(SOR)、强制隐式程序(SIP)和高斯-赛德尔法)相比,在地下水流问题上表现得较好[Kuiper, 1981, 1987; Jackson and Robinson, 1985; Gambolati and Perdon, 1984]。我们的目标也不是确定共轭梯度法的最佳预处理器。我们需要一种易于编程、高效的大规模问题求解方法。我们的实验记录了共轭梯度法在矢量(Cray X-MP/48)和矢量-并行(Alliant FX/8)计算机上的效率,并且展示了用于地下水流问题的多项式预处理器的有效性。