The parallel subdomain-levelset deflation method in reservoir simulation

DOI10.1016/J.JCP.2015.10.016MaRDI QIDQ2374901zbMATH OpenOpenAlexFDO

Authors J. H. van der Linden, T. B. Jönsthövel, A. A. Luk'yanov, C. Vuik

Publication date 5 December 2016

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1510.02148

zbMATH Keywords

preconditionersdeflationreservoir simulationextreme eigenvaluesharmonic Ritz deflationphysics-based deflation

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Numerical solution of discretized equations for boundary value problems involving PDEs (65N22) Flows in porous media; filtration; seepage (76S05)

Abstract: Extreme and isolated eigenvalues are known to be harmful to the convergence of an iterative solver. These eigenvalues can be produced by strong heterogeneity in the underlying physics. We can improve the quality of the spectrum by "deflating" the harmful eigenvalues. In this work, deflation is applied to linear systems in reservoir simulation. In particular, large, sudden differences in the permeability produce extreme eigenvalues. The number and magnitude of these eigenvalues is linked to the number and magnitude of the permeability jumps. Two deflation methods are discussed. Firstly, we state that harmonic Ritz eigenvector deflation, which computes the deflation vectors from the information produced by the linear solver, is unfeasible in modern reservoir simulation due to high costs and lack of parallelism. Secondly, we test a physics-based subdomain-levelset deflation algorithm that constructs the deflation vectors a priori. Numerical experiments show that both methods can improve the performance of the linear solver. We highlight the fact that subdomain-levelset deflation is particularly suitable for a parallel implementation. For cases with well-defined permeability jumps of a factor 104 or higher, parallel physics-based deflation has potential in commercial applications. In particular, the good scalability of parallel subdomain-levelset deflation combined with the robust parallel preconditioner for deflated system suggests the use of this method as an alternative for AMG.

Recommendations

• On POD-based deflation vectors for DPCG applied to porous media problems
• Accelerating the solution of linear systems appearing in two-phase reservoir simulation by the use of POD-based deflation methods
• A comparison of various deflation vectors applied to elliptic problems with discontinuous coefficients
• On the construction of deflation-based preconditioners
• The construction of projection vectors for a deflated ICCG method applied to problems with extreme contrasts in the coefficients

Cites work

• scientific article; zbMATH DE number 1561761 (Why is no real title available?)
• scientific article; zbMATH DE number 2164717 (Why is no real title available?)
• scientific article; zbMATH DE number 1406963 (Why is no real title available?)
• A comparison of various deflation vectors applied to elliptic problems with discontinuous coefficients
• A multiscale finite element method for elliptic problems in composite materials and porous media
• A parallel linear solver exploiting the physical properties of the underlying mechanical problem
• Adaptive multiscale finite-volume method for nonlinear multiphase transport in heterogeneous formations
• Algebraic multiscale solver for flow in heterogeneous porous media
• An efficient preconditioned CG method for the solution of a class of layered problems with extreme contrasts in the coefficients
• Conjugate gradient method with preconditioning by projector
• Damped Jacobi Preconditioning and Coarse Grid Deflation for Conjugate Gradient Iteration on Parallel Computers
• Deflated and Augmented Krylov Subspace Techniques
• Deflation of Conjugate Gradients with Applications to Boundary Value Problems
• Efficient deflation methods applied to 3-D bubbly flow problems
• Fast and robust solvers for pressure-correction in bubbly flow problems
• GMRES with Deflated Restarting
• GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
• Iterative Galerkin-enriched multiscale finite-volume method
• Matrix algorithms. Vol. 2: Eigensystems
• Methods of conjugate gradients for solving linear systems
• Multi-scale finite-volume method for elliptic problems in subsurface flow simulation.
• Multiscale Finite Element Methods
• New variants of deflation techniques for pressure correction in bubbly flow problems
• On the Conjugate Gradient Solution of the Schur Complement System Obtained from Domain Decomposition
• On the construction of deflation-based preconditioners
• On the use of rigid body modes in the deflated preconditioned conjugate gradient method
• Restarted GMRES preconditioned by deflation
• The construction of projection vectors for a deflated ICCG method applied to problems with extreme contrasts in the coefficients
• The superlinear convergence behaviour of GMRES

Cited in

(4)

• Zonal multiscale finite-volume framework
• Exploiting spatial symmetries for solving Poisson's equation
• Accelerating the solution of linear systems appearing in two-phase reservoir simulation by the use of POD-based deflation methods
• Subdomain deflation combined with local AMG: a case study using AMGCL library

This page was built for publication: The parallel subdomain-levelset deflation method in reservoir simulation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2374901)