Constructing projection vectors for DICCG
Gjalt Wijma

Plaats van afstuderen:
Shell Research
Volmerlaan 8
2280 AB Rijswijk
start van afstuderen: oktober 1997

De afstudeeropdracht is augustus 1998 afgerond met het schrijven van het afstudeerverslag. Huidige adresgegevens etc. zijn te vinden op onze alumnipagina.

Korte omschrijving van de afstudeeropdracht:
Shell wants to use CG for solving pressure problems in the upper three to six kilometers of the earth. Due to large contrasts in permeability the usual scaling methods like Incomplete Choleski are not powerful enough to ensure a smooth convergence behaviour. The mathematical reason for this is that Incomplete Choleski does not remove all of the small eigenvalues in the original system. There are a few eigenvalues left of order of the low permeability.

A projection method is introduced to remove these small eigenvalues. The projection method transforms the preconditioned system into a new one which kernel contains the eigenvectors belonging to the small eigenvalues. The new system can be solved and has a much better convergence bahviour than the preconditioned system. The projection method is practical in this case because the projection vectors have a clear link to the physical problem. Therefore it is possible to estimate them upfront. My research is focussed on the prediction and construction of these eigenvectors. The resulting method is implemented in SEPRAN.
Back to home page or Afstudeer pagina of Kees Vuik