An improved matrix split-iteration method for analyzing underground water flow

Shuairun Zhu

doi:10.1007/s00366-021-01551-z

2022-01-30

PRE

朱帅润 (Shuairun Zhu)

武立州 (Lizhou Wu) *

Xiaolin Song *

DOI：10.1007/s00366-021-01551-z

原文链接

原文求助

摘要

En 中文

The Hermitian and skew-Hermitian splitting iteration method (HSS) is commonly an effective linear iterative method for solving sparse non-Hermite positive definite equations. However, it is time-consuming to solve linear equations. Hence, inexact Hermitian and skew-Hermitian splitting iteration approaches with multistep preconditioner (PIHSS(m)) are proposed for analyzing underground water flow. For unsaturated porous media, an exponential model is adopted to linearize the Richards equation. The governing equations are discretized using the finite element method to produce a system of linear equations. Furthermore, the inexact Hermitian and skew-Hermitian splitting iteration methods (IHSS) and PIHSS(m) are used to solve the linear equations. The results show that PIHSS(m) can effectively solve the 1D unsaturated flow problem and 2D transient drainage problem in partially and completely saturated soils. The IHSS has higher numerical accuracy than the classical methods such as Picard method and Gauss-Seidel iterative method. Compared with IHSS, PIHSS(m) achieves faster convergence rate and higher computational efficiency, particularly for solving groundwater flow problems with high grid density. Additionally, the numerical results reveal that PIHSS(m) has excellent acceleration, that is, at least 50% acceleration compared with the IHSS.

Keyword:

Richards equation

Splitting iteration

Multistep preconditioner

Convergence rate

Acceleration