Abstract
AbstractWe 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.
Funder
Instituto Mexicano del Petróleo
Secretaría de Energía
Publisher
Springer Science and Business Media LLC
Subject
Computational Mathematics,Computational Theory and Mathematics,Computers in Earth Sciences,Computer Science Applications
Reference50 articles.
1. Gaul, A.: Recycling Krylov subspace methods for sequences of linear systems. Ph.D. thesis, Technical University of Berlin (2014)
2. Gaul, A., Gutknecht, M., Liesen, J., Nabben, R.: A framework for deflated and augmented Krylov subspace methods. SIAM J. Matrix Anal. Appl. 34(2), 495–518 (2013)
3. Quarteroni, A., Rozza, G.: Reduced Order Methods for Modeling and Computational Reduction, vol. 9. Springer, New York (2014)
4. Van der Sluis, A., van der Vorst, H.A.: The rate of convergence of conjugate gradients. Numer. Math. 48(5), 543–560 (1986)
5. Van der Sluis, A., Van der Vorst, H.A.: The convergence behavior of Ritz values in the presence of close eigenvalues. Linear Algebra Appl. 88, 651–694 (1987)
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献