Gegeben sei ein lineares Gleichungssystem $Au = f$ mit Koeffizientenmatrix $A$, welche eine spezielle block-tridiagonale Struktur besitzt. Solche lineare Gleichungssysteme entstehen bei der Diskretisierung dreidimensionaler elliptischer Randwertprobleme mit 7- oder 27-Punkte-Stern. In geophysikalischen Anwedungen, insbesondere bei Aufgaben aus der Geoelektrik, haben die Randwertprobleme unstetige Koeffizienten und sind meistens auf nicht-uniformen Gittern diskretisiert. Klassische geometrische Mehrgitterverfahren konvergieren um so langsamer, je stärker die Koeffizientensprünge ausfallen. Außerdem kann die Konvergenz durch die Variation der Gitterabstände beeinträchtigt werden. Zur Lösung wird ein matrix-abhängiges Mehrgitterverfahren vorgestellt. Als Glätter wird eine unvollständige Block LU-Zerlegung verwendet. Die Gittertransferoperationen werden anhand der Einträge der Matrix $A$ ermittelt. Das resultierende Verfahren erweist sich als sehr robust, insbesondere wenn es als Vorkonditionierung für das Verfahren der konjugierten Gradienten eingesetzt wird.
Identifer | oai:union.ndltd.org:DRESDEN/oai:qucosa.de:bsz:105-682918 |
Date | 16 December 2009 |
Creators | Zhebel, Elena |
Contributors | TU Bergakademie Freiberg, Mathematik und Informatik, Prof. Dr. Michael Eiermann, Prof. Dr. Michael Eiermann, Prof. Dr. Klaus Spitzer, Prof. Dr. Yuri Farkov |
Publisher | Technische Universitaet Bergakademie Freiberg Universitaetsbibliothek "Georgius Agricola" |
Source Sets | Hochschulschriftenserver (HSSS) der SLUB Dresden |
Language | English |
Detected Language | German |
Type | doc-type:doctoralThesis |
Format | application/pdf |
Page generated in 0.0027 seconds