Accelerating the solution of linear systems appearing in two-phase reservoir simulation by the use of POD-based deflation methods

Journal Article (2021)
Authors

Gabriela Berenice Diaz Cortes (Mexican Institute of Petroleum (IMP))

Kees Vuik (TU Delft - Numerical Analysis)

J. D. Jansen (TU Delft - Civil Engineering & Geosciences)

Research Group
Numerical Analysis
Copyright
© 2021 Gabriela Berenice Diaz Cortes, Cornelis Vuik, J.D. Jansen
To reference this document use:
https://doi.org/10.1007/s10596-021-10041-6
More Info
expand_more
Publication Year
2021
Language
English
Copyright
© 2021 Gabriela Berenice Diaz Cortes, Cornelis Vuik, J.D. Jansen
Research Group
Numerical Analysis
Issue number
5
Volume number
25
Pages (from-to)
1621-1645
DOI:
https://doi.org/10.1007/s10596-021-10041-6
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

We explore and develop a Proper Orthogonal Decomposition (POD)-based deflation method for the solution of ill-conditioned linear systems, appearing in simulations of two-phase flow through highly heterogeneous porous media. We accelerate the convergence of a Preconditioned Conjugate Gradient (PCG) method achieving speed-ups of factors up to five. The up-front extra computational cost of the proposed method depends on the number of deflation vectors. The POD-based deflation method is tested for a particular problem and linear solver; nevertheless, it can be applied to various transient problems, and combined with multiple solvers, e.g., Krylov subspace and multigrid methods.