LOGAN K. KUIPERf(1987)

A COMPARISON OF ITERATIVE METHODS AS APPLIED TO THE SOLUTION OF THE NONLINEAR THREE-DIMENSIONAL GROUNDWATER FLOW EQUATION*

LOGAN K. KUIPERf

**应用于非线性三维地下水流方程求解的迭代方法比较**

**LOGAN K. KUIPERf**

**摘要：**
本文比较了十七种不同方法的效率，用于求解三维水位含水层中地下水流的同时非线性有限差分逼近方程。解决方法包括牛顿、皮卡和广义线性方法。通过利用几种不同的线性方法，包括逐次过松弛、强烈隐式过程和八种不同的预处理共轭梯度方法，实施了牛顿和皮卡迭代方法。发现最佳方法是使用实施预处理共轭梯度方法的皮卡迭代方法。

**关键词：** 非线性方程组的迭代求解，地下水流

**1. 引言。**
地下水流的逼近方程通常使用迭代方法求解。这不仅适用于线性的受限含水层流方程，也适用于非线性的水位含水层流方程。解决非线性流方程的最常用程序是使用皮卡迭代（Ortega和Rheinboldt [7]；Remson等[9]），结合线性迭代方法，如交替方向隐式过程（ADI），逐次过松弛（SOR）或强烈隐式过程（SIP）（Trescott和Larson [10]）。关于用于解决水位含水层模拟的同步非线性方程的许多其他方法的性能调查在文献中很少见。作者（Kuiper [6]）比较了与SIP结合的皮卡迭代（Picard-SIP）和与不完全Cholesky共轭梯度方法结合的皮卡迭代（Picard-ICCG），应用于二维水位含水层。Remson等[9]的评论称赞Ortega和Rheinboldt [7]的一种称为一步点SOR-牛顿的方法，其命名和分类程序也将在本文中使用。Concus [2]对二维非线性静磁场问题进行的调查比较了一步牛顿点SOR，一步牛顿线SOR和一步点SOR-牛顿。在他的比较中，一步点SOR-牛顿与其他两种方法相比较有利。然而，这些结果可能不直接适用于水位含水层，因为两个问题的非线性差异显著。Ortega和Rockoff [8]提出了几种广义线性方法的比较。

本文中进行的比较可能意在回答与通常的数值方法比较不同的一组问题。要处理的非线性是在地下水流的逼近方程中找到的特定类型。此外，系统的几何形状和边界条件的类型也是具体的。这些数值测试旨在比较各种方法在一系列测试问题上的表现，这些问题在地下水流方面是全面的，但对于同时求解一般非线性逼近方程而言非常严格。即使方法相当接近，结果仍然被认为是有用的。