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

Journal Article (2021)
Author(s)

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

Cornelis Vuik (TU Delft - Numerical Analysis)

Jan Dirk Jansen (TU Delft - Civil Engineering & Geosciences)

DOI related publication
https://doi.org/10.1007/s10596-021-10041-6 Final published version
More Info
expand_more
Publication Year
2021
Language
English
Issue number
5
Volume number
25
Pages (from-to)
1621-1645
Downloads counter
199
Collections
Institutional Repository
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.