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

Bearbeiten

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

Bearbeiten

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

Bearbeiten

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

Bearbeiten

Vorgehensweise

Bearbeiten

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

Bearbeiten

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

Bearbeiten

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

Bearbeiten
  • 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.
Bearbeiten

Einzelnachweise

Bearbeiten
  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).