On Solving Groundwater Flow and Transport Models with Algebraic Multigrid Preconditioning

by   M. A. Sbai, et al.

Sparse iterative solvers preconditioned with the algebraic multigrid has been devised as an optimal technology to speed up the response of large linear systems. In this work, this technique was introduced into the framework of the dual delineation approach. This involves a single groundwater flow solution and a scalar advective transport solution with different right-hand sides. The method was compared with traditional preconditioned iterative methods and a direct sparse solver on several two- and three-dimensional benchmarks spanning homogeneous and heterogeneous formations. The algebraic multigrid preconditioning enabled speedups lying between one and two orders of magnitude for the groundwater flow problems. However, the sparse direct solver was the fastest for the pure advective transport processes such as the forward travel time simulations. This leads to conclude that the best sparse solver for the general advection-dispersion transport equation is Péclet number dependent. When equipped with the best solvers, the dual delineation technique was run on a took only a few seconds to solve multi-million grid cell problems paving the way for comprehensive sensitivity analysis. The paper gives practical hints on the strategies and conditions under which the algebraic multigrid preconditioning for nonlinear and/or transient problems would remain competitive.



page 22

page 24


Effective Preconditioners for Mixed-Dimensional Scalar Elliptic Problems

Discretization of flow in fractured porous media commonly lead to large ...

Efficient Reordered Nonlinear Gauss-Seidel Solvers With Higher Order For Black-Oil Models

The fully implicit method is the most commonly used approach to solve bl...

Learning Algebraic Multigrid Using Graph Neural Networks

Efficient numerical solvers for sparse linear systems are crucial in sci...

Acceleration of Radiation Transport Solves Using Artificial Neural Networks

Discontinuous Finite Element Methods (DFEM) have been widely used for so...

AMPS: An Augmented Matrix Formulation for Principal Submatrix Updates with Application to Power Grids

We present AMPS, an augmented matrix approach to update the solution to ...

Multigrid-augmented deep learning preconditioners for the Helmholtz equation

In this paper, we present a data-driven approach to iteratively solve th...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.