The construction of projection vectors for a deflated ICCG method applied to problems with extreme contrasts in the coefficients

To predict the presence of oil and natural gas in a reservoir, it is import ant to know the fluid pressure in the rock formations. A mathematical model for the prediction of the fluid pressure history is given by a time-depend ent diffusion equation. Application of the finite-element method leads to s ystems of linear equations. A complication is that the underground consists of layers with very large contrasts in permeability. This implies that the symmetric and positive definite coefficient matrix has a very large condit ion number. Bad convergence behavior of the ICCG method has been observed, and a classical termination criterion is not valid in this problem. In [19] we have shown that the number of small eigenvalues of the diagonally scale d matrix is equal to the number of high-permeability domains, which are not connected to a Dirichlet boundary. In this paper the proof is extended to an Incomplete Cholesky decomposition. To annihilate the bad effect of these small eigenvalues on the convergence. the Deflated ICCG method is used. In [19] we have shown how to construct a deflation subspace for the case of a set of more or less parallel layers. That subspace proved to be a good app roximation of the span of the "small" eigenvectors. As a result of this, th e convergence Of DICCG is independent of the contrasts in the permeabilitie s. In this paper it is shown how to construct deflation vectors even in the case of very irregular shaped layers, and layers with so-called inclusions . A theoretical investigation and numerical experiments show that the DICCG method is not sensitive to small perturbations of the deflation vectors. T he efficiency of the DICCG method is illustrated by numerical experiments. (C) 2001 Academic Press.