Modelling unsaturated flow in porous media using an improved iterative scheme

Richards' equation is often used in unsaturated flow problems and has a wide range of applications. In the numerical solution, the Richards' equation is linearized first, and then the finite difference method is used for numerical discretization and iterative calculation. The traditional iterative methods such as Jacobi, Gauss-Seidel (GS) and SOR iterative methods have a slower convergence rate, especially when the discrete space step is small and the time step is large. Therefore, we adopt the integral correction method and the multistep preconditioner to improve the traditional iterative methods and propose an improved Gauss-Seidel iterative method (ICMP(m)-GS) with multistep preconditioner based on the integral correction method to solve the linear equations derived from linearized Richards' equation. Through examples of unsaturated flow, convergence rate and computational accuracy of the proposed algorithm are validated by comparing the traditional methods and analytical solutions. The results show that the proposed ICMP(m)-GS can greatly improve the ill-condition of linear equations. Compared with the conventional methods GS, SOR and a single improvement method, ICMP(m)-GS has a faster convergence rate, higher calculation efficiency and calculation accuracy. The application example shows that the proposed method can also well simulate the time-varying law of pressure head in the rainfall infiltration of unsaturated soil slopes, and has a nice application effect.

Automated summary

Richards' equation is widely used in unsaturated flow problems and its numerical solution can be improved by linearization and the use of integral correction method and multistep preconditioner. The proposed improved Gauss-Seidel iterative method (ICMP(m)-GS) demonstrates faster convergence rate, higher calculation efficiency and accuracy compared to traditional methods. The method is validated through examples of unsaturated flow and shows promising results in simulating rainfall infiltration of unsaturated soil slopes.