MODFLOW的PCG4 Package引入了一个基于部分LU分解预处理的PCG求解器，旨在提高复杂场地模拟的效率和稳健性。该方法基于Meijerink和van der Vorst的两步程序，适用于解决大型稀疏对称矩阵方程。PCG4输入包括最大外部和内部迭代次数、回溯因子、残差减少标准等参数，适用于线性和非线性问题。
6.1 GENERAL
MODFLOW代码以有限差分形式生成描述地下水流系统的方程组。通常通过使用Slice Successive Over-Relaxation (SSOR)或Strongly Implicit Procedure (SIP) Package来求解这组代数方程。对于复杂的场地模拟，这些求解程序可能缺乏效率和稳健性。Hill（1994）提供了一个Preconditioned Conjugate Gradient (PCG2) Package，为MODFLOW提供了另一种求解器。
PCG2中的预处理方案有Saad（1985）提出的最小二乘多项式预处理器或Meyer等人（1989）提出的最优Chebyshev多项式预处理器。这些方案主要基于计算机存储的考虑选择，通常在大规模场地研究中性能较差。鉴于这一点，我们加入了PCG4 Package，其中包含一个基于部分LU分解作为预处理器的PCG求解器。与MODFLOW先前使用的求解器相比，计算机存储要求更多。然而，PCG4求解器简单、稳健且高效。该求解器已在许多不同规模和复杂性的地下水建模项目中使用。本章的其余部分致力于制定PCG4 Package的公式和其使用的输入说明。
6.2 FORMULATION OF STANDARD PCG SCHEME
预处理共轭梯度（PCG）方法（参见Kershaw，1978）已成为解决由有限元或有限差分逼近生成的大型稀疏对称矩阵方程的最有前途的迭代方法之一，这些逼近适用于多维场问题。许多最近的论文描述了PCG方案与有限元和有限差分方法相结合成功应用于地下水流问题。这些包括Gambolati（1980）、Kuiper（1981, 1987）、Gambolati等人（1986）、Kalf（1988）和Meyer等人（1
最低0.47元/天 解锁文章
扫一扫
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
