An Efficient Preconditioned CG Method for the Solution of a Class of Layered Problems with Extreme Contrasts in the Coefficients

Knowledge of fluid pressure is important to predict the presence of oil and gas in reservoirs. A mathematical model for the prediction of fluid pressures is given by a time-dependent diffusion equation. Application of the finite element method leads to a system of linear equations. A complication is that the underground consists of layers with very large differences in permeability. This implies that the symmetric and positive definite coefficient matrix has a very large condition number. Bad convergence behavior of the CG method has been observed; moreover, a classical termination criterion is not valid in this problem. After diagonal scaling of the matrix the number of extreme eigenvalues is reduced and it is proved to be equal to the number of layers with a high permeability. For the IC preconditioner the same behavior is observed. To annihilate the effect of the extreme eigenvalues a deflated CG method is used. The convergence rate improves considerably and the termination criterion becomes again reliable. Finally a cheap approximation of the eigenvectors is proposed.

[1]  Frédéric Guyomarc'h,et al.  A Deflated Version of the Conjugate Gradient Algorithm , 1999, SIAM J. Sci. Comput..

[2]  Ian G. Graham,et al.  Unstructured Additive Schwarz-Conjugate Gradient Method for Elliptic Problems with Highly Discontinuous Coefficients , 1999, SIAM J. Sci. Comput..

[3]  Alan B. Williams,et al.  A Deflation Technique for Linear Systems of Equations , 1998, SIAM J. Sci. Comput..

[4]  Melvyn R. Giles,et al.  Diagenesis: A Quantitative Perspective: Implications for Basin Modelling and Rock Property Prediction , 1997 .

[5]  Cornelis Vuik,et al.  How fast the Laplace equation was solved in 1995 , 1997 .

[6]  Peter Deuflhard,et al.  Domain decomposition with subdomain CCG for material jump elliptic problems , 1997 .

[7]  Y. Saad.  Analysis of Augmented Krylov Subspace Methods , 1997, SIAM J. Matrix Anal. Appl..

[8]  Yousef Saad,et al.  ILUM: A Multi-Elimination ILU Preconditioner for General Sparse Matrices , 1996, SIAM J. Sci. Comput..

[9]  K. Burrage,et al.  Restarted GMRES preconditioned by deflation , 1996 .

[10]  F. Wubs,et al.  Nested grids ILU-decomposition (NGILU) , 1996 .

[11]  Ronald B. Morgan,et al.  A Restarted GMRES Method Augmented with Eigenvectors , 1995, SIAM J. Matrix Anal. Appl..

[12]  O. Axelsson.  Iterative solution methods , 1995 .

[13]  Yvan Notay,et al.  DRIC: A dynamic version of the RIC method , 1994, Numer. Linear Algebra Appl..

[14]  A. V. Popov,et al.  Efficient Computational Schemes of the Conjugate Gradient Method for Solving Linear Systems , 1994 .

[15]  Y. Notay,et al.  Incomplete factorizations of singular linear systems , 1989 .

[16]  E. F. Kaasschieter,et al.  Preconditioned conjugate gradients for solving singular systems , 1988 .

[17]  E. F. Kaasschieter,et al.  A practical termination criterion for the conjugate gradient method , 1988 .

[18]  R. Nicolaides.  Deflation of conjugate gradients with applications to boundary value problems , 1987 .

[19]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[20]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[21]  O. Axelsson,et al.  On the eigenvalue distribution of a class of preconditioning methods , 1986 .

[22]  A. Brandt,et al.  The Multi-Grid Method for the Diffusion Equation with Strongly Discontinuous Coefficients , 1981 .

[23]  W. Greub.  Linear Algebra , 1981 .

[24]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[25]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[26]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[27]  H. Walker,et al.  GMRES On (Nearly) Singular Systems , 1997, SIAM J. Matrix Anal. Appl..

[28]  Yousef Saad,et al.  Deflated and Augmented Krylov Subspace Techniques , 1997, Numer. Linear Algebra Appl..

[29]  A. van der Ploeg,et al.  Nested grids ILU-decomposition , 1996 .

[30]  S. A. Kharchenko,et al.  Eigenvalue translation based preconditioners for the GMRES(k) method , 1995, Numer. Linear Algebra Appl..

[31]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[32]  E. Polak.  Introduction to linear and nonlinear programming , 1973 .