Galerkin-Methode

Die Galerkin-Methode (auch Galerkin-Verfahren bzw. Galerkin-Ansatz, nach Boris Galerkin, 1915) ist ein numerisches Verfahren zur näherungsweisen Lösung von Operatorgleichungen, wie beispielsweise 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.

Kurzfassung

John William Strutt und Walter Ritz haben die in Variationsproblemen gesuchte Funktion als Linearkombination von Basisfunktionen angesetzt und damit das Variationsproblem auf ein gewöhnliches Problem der Optimierung einer Funktion von endlich vielen Parametern zurückgeführt.[1]

Für eine Operatorgleichung

kann die gesuchte Funktion ebenso angesetzt werden, etwa als

was substituiert in die Operatorgleichung auf der linken Seite des Gleichheitszeichens eine von den Koeffizienten abhängige Funktion ergibt. Nach der Methode der gewichteten Residuen wählt man die freien Koeffizienten so, dass diese Funktion im Testraum, der von gewissen Basisfunktionen aufgespannt wird, verschwindet, d. h. orthogonal zu diesen Basisfunktionen wird. Damit erhält man folgende Gleichungen für alle Kombinationen

zur Bestimmung der . Diese Bedingung wird auch Galerkin-Orthogonalität genannt. Falls der Operator linear ist, lassen sich diese Gleichungen als ein lineares Gleichungssystem darstellen mit

Für , wobei eine feste Funktion ist (häufig eine Radiale Basisfunktion), erhält man ein Punkt-Kollokationsverfahren. Für erhält man das Galerkin-Verfahren, das vor allem in russischen Büchern auch Iwan Grigorjewitsch Bubnow (1911, 1913) zugeschrieben wird, dort also Bubnov-Galerkin-Verfahren heißt.[2]

Herleitung

Ausgangspunkt für die Galerkin-Methode ist eine sogenannte „variationelle“ Formulierung der Anfangswertaufgabe.[3]

Sei also die Anfangswertaufgabe (AWA) gegeben mit auf einem Intervall

Außerdem nehmen wir die AWA als (global) Lipschitz-stetig an. D. h., dass eindeutige Lösungen existieren. (Satz von Picard-Lindelöf)

Dann wird die Differentialgleichung der Anfangswertaufgabe nun zunächst mit einer „Testfunktion“ multipliziert und über das Lösungsintervall integriert. Wir erhalten aus der AWA

Diese Beziehung ist sinnvoll für jede stetige und stückweise stetig differenziertere Funktion . Der Vektorraum all dieser Funktionen sei ab hier mit benannt. „Stückweise“ bedeutet hier, dass die Differenzierbarkeit nur bis auf endlich viele Ausnahmestellen in gefordert wird. Das linke Integral ist stückweise als Summe von Teilintegralen zu verstehen.[3]

Jede Funktion , die die Anfangsbedingungen des Startwertes erfüllt und der integralen Beziehung für jede Testfunktion genügt, ist auch Lösung der Anfangswertaufgabe.

Die Galerkin-Methode bestimmt eine Näherungslösung in einem endlich dimensionalen Teilraum durch die Vorschriften des Startpunktes und der Integralgleichung im Teilraum: für ein beliebiges .

Der diskrete Testraum ist in der Regel anders als zu wählen. Wähle bspw.

Man kann die integrale Bestimmungsgleichung auf jedes einzelne Teilintervall einschränken, da die Testfunktionen nur stückweise stetig sein müssen.[3]

Das bedeutet: Die Galerkin-Methode ist ein „Zeitschrittverfahren“.[3] Wertet man bspw. das Integral auf der rechten Seite mit der Trapezregel aus, dann erhalten wir für die Werte

Zur „variationellen“ Formulierung der Anfangswertaufgabe

Ausgangspunkt sind, wie im Abschnitt der Herleitung, die AWA mit Lipschitz-Bedingung und . Auftretende Funktionen können auch vektorwertig sein und bezeichnet das euklidische Skalarprodukt. Für eine Funktion , also eine einmal komplexwertig-differenziertere Funktion mit Dimension d (da jeder Fall höhere Ordnung auf den Fall erster Ordnung zurückgeführt werden kann) mit dem Anfangswert ist die AWA und die äquivalente Formulierung in der Herleitung äquivalent zu:

Da die Funktionen beliebig variieren dürfen, wird diese Formulierung der AWA „variationell“ genannt.[4]

Geometrisch ausgedrückt besagt die variationelle Formulierung der AWA, dass das Residuum der Lösungsfunktion :

bzgl. des Skalarproduktes von orthogonal zu allen Testfunktionen [4]

Detailliertere Darstellung

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 ein Differentialoperator. Gesucht ist die Lösung der Differentialgleichung:

(Gleichung 1)

mit einer vorgegebenen Funktion und zusätzlich Randbedingungen für . Dazu wird eine Näherungslösung für angesetzt als Linearkombination von Basisfunktionen aus einem Funktionenraum :

mit noch zu bestimmenden Koeffizienten . Für die Wahl der Basis eignen sich zum Beispiel B-Splines. Die Funktion erfüllt im Allgemeinen noch nicht die Differentialgleichung (1), es bleibt ein Residuum

In dem Raum sei ein inneres Produkt definiert, für das gilt, dass ist, falls für alle Funktionen aus ist. Das innere Produkt ist häufig definiert als

Häufig kann man nicht die exakte Lösung bestimmen, für die für jede Testfunktion 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“ verschwindet:

Beim Galerkin-Verfahren werden als Gewichtsfunktionen gerade die Basisfunktionen , von gewählt, so dass sich ein Gleichungssystem für die Koeffizienten ergibt:

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 Olgierd Cecil 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.[5]

Weiterführende Literatur

  • H. R. Schwarz: Methode der Finiten Elemente (= Leitfäden der angewandten Mathematik und Mechanik. Band 47). Teubner, Stuttgart 1984, ISBN 3-519-12349-5.
  • Olgierd Cecil Zienkiewicz, Robert Lee Taylor: The Finite Element Method. 4. Auflage, Band 1: Basic Formulation and Linear Problems. McGraw-Hill Book Company, London 1989, ISBN 0-07-084174-8.
  • Junuthula Narasimha Reddy: Energy Principles And Variational Methods In Applied Mechanics. 2. Auflage, John Wiley & Sons, New York 2002, ISBN 0-471-17985-X (Leseprobe, books.google.de).
  • Daniel S. Weile, Raymond A. Wildman, Greeshma Pisharody, Anuraag Mohan: Galerkin Method (Rayleigh–Ritz Method). In: Encyclopedia of RF and Microwave Engineering. John Wiley, Hoboken, N.J. 15. April 2005, ISBN 0-471-65450-7, doi:10.1002/0471654507.eme142.
  • Rolf Rannacher: Numerik 1. Numerik gewöhnlicher Differentialgleichungen (= Lecture Notes Mathematik). Universitätsverlag, Heidelberg 2017, ISBN 978-3-946054-32-0, S. 9–10, doi:10.17885/heiup.258.342.

Einzelnachweise

  1. Zienkiewicz: The Finite Element Method. 4. Auflage, Band 1, S. 35.
  2. Zienkiewicz: The Finite Element Method. 4. Auflage, Band. 1, S. 215.
  3. a b c d Rolf Rannacher: Numerik 1. Numerik gewöhnlicher Differentialgleichungen. Heidelberg 2017, S. 9–10.
  4. a b Rolf Rannacher: Numerik 1. Numerik gewöhnlicher Differentialgleichungen. Heidelberg 2017, S. 151 f.
  5. Finite-Elemente-Methode. und Galerkin-Methode. In: Guido Walz (Hrsg.): Lexikon der Mathematik. Band 2: Eig bis Inn. Springer Spektrum, Mannheim 2017, ISBN 978-3-662-53503-5, S. 157 S. 227–228, doi:10.1007/978-3-662-53504-2 (archive.org).