Galerkin-Methode

Die Galerkin-Methode (auch Galerkin-Ansatz, nach Boris Galerkin) ist ein numerisches Verfahren zur näherungsweisen Lösung von partiellen Differentialgleichungen. Sie stellt die gebräuchlichste Variante der "Methode der gewichteten Residuen" dar, bei der das resultierende Residuum einer Näherungslösung minimiert wird. Der Ausdruck "Methode der gewichteten Residuen" wird häufig auch synonym mit "Galerkin-Methode" benutzt.

Vorgehensweise

Das Residuum ist in dem betrachteten Gebiet verteilt. Es wird mit geeigneten Wichtungsfunktionen gewichtet, daher der Ausdruck "gewichtete Residuen". Das Integral des über dem Gebiet gewichteten Residuums soll möglichst klein sein oder besser noch ganz verschwinden. Die Wichtungsfunktionen haben Parameter, deren Anzahl der Zahl der Freiheitsgrade des Systems entspricht. Diese führen zu genauso vielen Gleichungen und damit zu dem gleichen großen Gleichungssystem, das aus der Finite-Elemente-Methode bekannt ist. Bei der Galerkin-Methode sind die Wichtungsfunktionen identisch mit den Ansatzfunktionen in den Elementen.

Beispiel

Sei D ein Differentialoperator. Gesucht ist die Lösung u(x) der Differentialgleichung

D(u(x)) + f(x) = 0 (Gleichung 1)

mit einer vorgegebenen Funktion f(x) und zusätzlich Randbedingungen für u. Dazu wird eine Näherungslösung v(x) für u angesetzt als Linearkombination von Basisfunktionen Φ_i(x) aus einem Funktionen-Vektorraum V:

$v (x)= \sum_{i=1}^N c_i \Phi_i (x)$

mit noch zu bestimmenden Koeffizienten c_i. v(x) erfüllt im Allgemeinen noch nicht die Differentialgleichung (1), es bleibt ein Residuum:

r(x) = D(v(x)) + f(x)

In dem Raum V von Funktionen ist ein inneres Produkt < h,g > definiert, für das gilt, dass g=0 ist falls < h,g > = 0 für beliebige Funktionen h aus V. Das innere Produkt ist häufig definiert als

$\langle h,g \rangle = \int h(x) g(x) dx$

Häufig kann man nicht die exakte Lösung bestimmen, für die < w,r > für jede Testfunktion w(x) verschwindet (und damit das Residuum auch), sondern nur eine Näherungslösung, für die das innere Produkt des Residuums mit einer Menge ausgewählter linear unabhängiger „Gewichtsfunktionen“ w(x) verschwindet:

$\langle w(x), r(x) \rangle =0.$

Beim Galerkin Verfahren werden als Gewichtsfunktionen gerade die Basisfunktionen von V Φ_j (j=1, .., N) gewählt, so dass sich ein Gleichungssystem für die Koeffizienten c_i ergibt.

$\langle \Phi_j, D ( \sum_{i=1}^N c_i \Phi_i (x)) + f(x) \rangle = 0.$

Anwendungsgebiet

Die Galerkin-Methode ist anwendbar, wenn kein natürliches Extremalprinzip für die Lösung der Differentialgleichung existiert. Sie ist somit eine Grundlage der Finite-Elemente-Methode und dehnt deren Anwendbarkeit auf weitere physikalische Problemstellungen (Kontinuumsprobleme) aus, die ein solches natürliches Extremalprinzip nicht besitzen. Beispiele dafür sind stationäre oder instationäre Strömungen. Ein natürliches Extremalprinzip (natürliches Variationsprinzip) existiert dagegen z. B. bei mechanischen Problemen der Festkörpermechanik, bei denen der Energieinhalt ein Minimum haben muss.

Nach Zienkiewicz ist die Galerkin-Lösung identisch mit einer natürlichen Variationslösung oder lässt sich zumindest so interpretieren. Die Finite-Elemente-Methode (FEM) ist ein spezielles Ritz-Galerkin-Verfahren.

Weiterführende Literatur

• O. C. Zienkiewicz: Methode der Finiten Elemente. (The Finite Element Method) 1977.
• J.N. Reddy: Energy Principles And Variational Methods In Applied Mechanics. Second Edition, John Wiley & Sons 2002, ISBN 978-0-471-17985-6.
• H. R. Schwarz: Methode der Finiten Elemente. Stuttgart 1984.

Weblinks

• wikibooks: Ritz-Galerkin-Verfahren
• Methode der gewichteten Residuen
• Galerkin's Method
• Galerkin Method, Springer Online Reference