Solver preconditioning using the combinatorial multilevel method on reservoir simulation

Yuhe Wang, John E. Killough

Research output: Contribution to journalArticle

Abstract

The purpose of this paper is to report the first preliminary study of the recently introduced combinatorial multilevel (CML) method for solver preconditioning in large-scale reservoir simulation with coupled geomechanics. The CML method is a variant of the popular algebraic multigrid (AMG) method yet with essential differences. The basic idea of this new approach is to construct a hierarchy of matrices using the discrete geometry of the graph, based on support theory for preconditioners. In this way, the CML method combines the merits of both geometric and algebraic multigrid methods. The resulting hybrid approach not only provides a simpler and faster setup phase compared to AMG, but the method can be proven to exhibit strong convergence guarantees for arbitrary symmetric diagonally dominant matrices. In addition, the underlying theoretical soundness of the CML method contrasts to the heuristic AMG approach, which often can show slow convergence for difficult problems. This new approach is implemented in a reservoir simulator for both pressure and poroelastic displacement preconditioners in the multistage preconditioning technique. We present results based on several known benchmark problems and provide a comparison of performance and complexity with the widespread preconditioning schemes used in large-scale reservoir simulation. An adaptation of CML for non-symmetric matrices is shown to exhibit excellent convergence properties for realistic cases.

Original languageEnglish
Pages (from-to)695-708
Number of pages14
JournalComputational Geosciences
Volume19
Issue number4
DOIs
Publication statusPublished - 22 Aug 2015

Fingerprint

Reservoir Simulation
Multilevel Methods
Preconditioning
Algebraic multigrid Method
Algebraic multigrid
Geomechanics
simulation
Preconditioner
Diagonally Dominant Matrix
Discrete Geometry
Nonsymmetric Matrix
matrix
Preconditioning Techniques
Simulators
Hybrid Approach
Soundness
Strong Convergence
Symmetric matrix
Geometry
Convergence Properties

Keywords

  • Combinatorial multilevel method
  • Linear solver
  • Preconditioner
  • Reservoir simulation

ASJC Scopus subject areas

  • Computer Science Applications
  • Computers in Earth Sciences
  • Computational Theory and Mathematics
  • Computational Mathematics

Cite this

Solver preconditioning using the combinatorial multilevel method on reservoir simulation. / Wang, Yuhe; Killough, John E.

In: Computational Geosciences, Vol. 19, No. 4, 22.08.2015, p. 695-708.

Research output: Contribution to journalArticle

@article{d592d7847d544ab9958f9bdb05c578f6,
title = "Solver preconditioning using the combinatorial multilevel method on reservoir simulation",
abstract = "The purpose of this paper is to report the first preliminary study of the recently introduced combinatorial multilevel (CML) method for solver preconditioning in large-scale reservoir simulation with coupled geomechanics. The CML method is a variant of the popular algebraic multigrid (AMG) method yet with essential differences. The basic idea of this new approach is to construct a hierarchy of matrices using the discrete geometry of the graph, based on support theory for preconditioners. In this way, the CML method combines the merits of both geometric and algebraic multigrid methods. The resulting hybrid approach not only provides a simpler and faster setup phase compared to AMG, but the method can be proven to exhibit strong convergence guarantees for arbitrary symmetric diagonally dominant matrices. In addition, the underlying theoretical soundness of the CML method contrasts to the heuristic AMG approach, which often can show slow convergence for difficult problems. This new approach is implemented in a reservoir simulator for both pressure and poroelastic displacement preconditioners in the multistage preconditioning technique. We present results based on several known benchmark problems and provide a comparison of performance and complexity with the widespread preconditioning schemes used in large-scale reservoir simulation. An adaptation of CML for non-symmetric matrices is shown to exhibit excellent convergence properties for realistic cases.",
keywords = "Combinatorial multilevel method, Linear solver, Preconditioner, Reservoir simulation",
author = "Yuhe Wang and Killough, {John E.}",
year = "2015",
month = "8",
day = "22",
doi = "10.1007/s10596-015-9485-8",
language = "English",
volume = "19",
pages = "695--708",
journal = "Computational Geosciences",
issn = "1420-0597",
publisher = "Springer Netherlands",
number = "4",

}

TY - JOUR

T1 - Solver preconditioning using the combinatorial multilevel method on reservoir simulation

AU - Wang, Yuhe

AU - Killough, John E.

PY - 2015/8/22

Y1 - 2015/8/22

N2 - The purpose of this paper is to report the first preliminary study of the recently introduced combinatorial multilevel (CML) method for solver preconditioning in large-scale reservoir simulation with coupled geomechanics. The CML method is a variant of the popular algebraic multigrid (AMG) method yet with essential differences. The basic idea of this new approach is to construct a hierarchy of matrices using the discrete geometry of the graph, based on support theory for preconditioners. In this way, the CML method combines the merits of both geometric and algebraic multigrid methods. The resulting hybrid approach not only provides a simpler and faster setup phase compared to AMG, but the method can be proven to exhibit strong convergence guarantees for arbitrary symmetric diagonally dominant matrices. In addition, the underlying theoretical soundness of the CML method contrasts to the heuristic AMG approach, which often can show slow convergence for difficult problems. This new approach is implemented in a reservoir simulator for both pressure and poroelastic displacement preconditioners in the multistage preconditioning technique. We present results based on several known benchmark problems and provide a comparison of performance and complexity with the widespread preconditioning schemes used in large-scale reservoir simulation. An adaptation of CML for non-symmetric matrices is shown to exhibit excellent convergence properties for realistic cases.

AB - The purpose of this paper is to report the first preliminary study of the recently introduced combinatorial multilevel (CML) method for solver preconditioning in large-scale reservoir simulation with coupled geomechanics. The CML method is a variant of the popular algebraic multigrid (AMG) method yet with essential differences. The basic idea of this new approach is to construct a hierarchy of matrices using the discrete geometry of the graph, based on support theory for preconditioners. In this way, the CML method combines the merits of both geometric and algebraic multigrid methods. The resulting hybrid approach not only provides a simpler and faster setup phase compared to AMG, but the method can be proven to exhibit strong convergence guarantees for arbitrary symmetric diagonally dominant matrices. In addition, the underlying theoretical soundness of the CML method contrasts to the heuristic AMG approach, which often can show slow convergence for difficult problems. This new approach is implemented in a reservoir simulator for both pressure and poroelastic displacement preconditioners in the multistage preconditioning technique. We present results based on several known benchmark problems and provide a comparison of performance and complexity with the widespread preconditioning schemes used in large-scale reservoir simulation. An adaptation of CML for non-symmetric matrices is shown to exhibit excellent convergence properties for realistic cases.

KW - Combinatorial multilevel method

KW - Linear solver

KW - Preconditioner

KW - Reservoir simulation

UR - http://www.scopus.com/inward/record.url?scp=84941879456&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84941879456&partnerID=8YFLogxK

U2 - 10.1007/s10596-015-9485-8

DO - 10.1007/s10596-015-9485-8

M3 - Article

VL - 19

SP - 695

EP - 708

JO - Computational Geosciences

JF - Computational Geosciences

SN - 1420-0597

IS - 4

ER -