Ein hierarchischer iterativer Löseransatz, der die Konvergenz für diskretisierte lineare (und einige nichtlineare) PDEs beschleunigt, indem er auf einer Folge von Diskretisierungsstufen arbeitet: Glättung hochfrequenter Fehler auf feinen Gittern, Übertragung der Residuen auf gröbere Gitter zur Eliminierung niederfrequenter Fehler, Korrektur der Approximationen und Interpolation der Korrekturen zurück auf feinere Stufen.

Prinzip

Kombiniere komplementäre Relaxation (Glättung) auf feinen Gittern mit Grobgitterkorrektur, um Fehlerkomponenten in ihren natürlichen Skalen zu behandeln: Hochfrequente Fehler werden durch lokale Relaxationen gedämpft, während niederfrequente Fehler effizient auf gröberen Repräsentationen adressiert werden, was zu schneller, oft problemgrößenunabhängiger Konvergenz führt.

Demonstration

Lösung einer Poisson-Gleichung auf einem uniformen Gitter: Wende einige Gauss–Seidel-Relaxationen (Smoother) auf dem feinen Gitter an, um oszillatorische Fehler zu dämpfen, restriktiere das Residuum auf ein gröberes Gitter und löse es dort (ungefähr), prolongiere die Grobgitterkorrektur auf das feine Gitter und iteriere in V- oder W-Zyklen bis die Residuen die Toleranz erfüllen.

Fehlanwendung

Multigrid mit einem ungeeigneten Smoother verwenden (z. B. einer Relaxation, die hochfrequente Moden für die gegebene Diskretisierung nicht dämpft), oder zwischen inkompatiblen Diskretisierungen ohne konsistente Restriktions-/Prolongationsoperatoren transferieren, was zu Stagnation oder Divergenz statt Beschleunigung führt.

Konsequenz

Richtig entworfen erreichen Multigrid-Verfahren konvergenzraten, die netzunabhängig sind, und oft optimale O(N)-Rechenkomplexität für elliptische Operatoren, sodass sehr große diskretisierte PDE-Systeme effizient in Zeit und Speicher gelöst werden können.

Umkehrung

Einstufige iterative Methoden ohne Grobgitterkorrektur: sich nur auf Relaxation (z. B. Jacobi, Gauss–Seidel) verlassen, die typischerweise hochfrequente Fehler reduziert, aber für glatte (niederfrequente) Komponenten langsam konvergiert, was bei wachsender Problemgröße netzabhängige Verlangsamung verursacht.

Abgrenzung

Geltungsbereich: diskretisierte elliptische und viele parabolische PDEs, für die Smoother und Transferoperatoren definiert werden können; ausgenommen sind Probleme mit stark nichtlokaler Kopplung, die keine wirksame Grobgitterkorrektur erlauben, schlecht konditionierte Diskretisierungen ohne geeignete Smoother und Konfigurationen, in denen Grobgittermodelle inkonsistent oder nicht verfügbar sind.

Semantische Spannung

Spannung zwischen Komplexität der Intergitteroperatoren und Robustheit: aufwändige Transfer- und Grobgitteroperatoren verbessern Robustheit und operatorabhängige Konvergenz, erhöhen aber Implementierungsaufwand und Aufwand pro Zyklus, während einfachere Wahlmöglichkeiten billigere Zyklen liefern, jedoch bei schwierigen Problemen versagen können.

Synthese

Das Mehrgitterverfahren ist ein mehrstufiger Algorithmus, der lokales Glätten auf feinen Gittern mit Grobgitterkorrektur durch Restriktion und Prolongation kombiniert, um alle Fehlerfrequenzen effizient zu reduzieren und schnelle, oft netzunabhängige Konvergenz für eine breite Klasse diskretisierter PDEs zu erzielen.