Finite-Differenzen-Methode

Numerisches Verfahren zur Lösung von Differentialgleichungen durch Näherung mittels Differenzengleichungen
(Weitergeleitet von Differenzenverfahren)

Finite-Differenzen-Methoden (FDM), auch Finite-Differenzen-Verfahren, Methoden/Verfahren der endlichen Differenzen oder schlicht Differenzenverfahren, sind eine Klasse numerischer Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen (DGL).

Perspektivische Grafik der numerischen Lösung einer partiellen DGL in 2D nach der FDM in Zylinderkoordinaten für ein physikalisches Problem[1]

Die grundlegende Idee des Verfahrens ist es, die Orts- und/oder Zeitableitungen in der Differentialgleichung in einem vorgegebenen Intervall der unabhängigen Variablen an Gitterpunkten in diesem Intervall durch Differenzenquotienten zu approximieren. Diese Approximationen der Differenzialgleichung in den Gitterpunkten stellen dann ein Gleichungssystem dar, welches zu lösen ist.

Die erste Grafik illustriert die numerische Lösung einer partiellen DGL in 2D nach der FDM in einem krummlinigen -Koordinatensystem (Zylinderkoordinaten) mit 21 Gitterpunkten (20 Intervallen) in der r-Richtung und 11 Gitterpunkten (10 Intervallen) in der z-Richtung. Der diskretisierte Körper ist ein Zylinder mit einer Höhe von 100 cm und mit einem Radius von 200 cm. Es handelt sich eigentlich um einen (realen) 3D-Zylinder, der Besonderheiten besitzt. Die in die Berechnung eingehende Materialverteilung im Zylinder kann sich in Prinzip von Zelle zu Zelle[2] ändern. In dem vorliegenden Modell hängt jedoch die Materialverteilung nicht von der polaren Koordinate, der -Koordinate, ab. Diese Eigenschaft muss prinzipiell auch jede physikalische Lösungsfunktion besitzen, die eine Observable ist.

Somit wird das Problem zu einem mathematischen 2D-, einem -Problem. Da die Materialverteilung außerdem in z-Richtung eine Spiegelsymmetrieebene besitzt, wurde für als eine der Randbedingungen der Wert für den Gradienten vorgegeben und damit die Hälfte der Gitterpunkte eingespart. Der Zylinder ist also eigentlich 200 cm hoch. Als weitere Randbedingung wurde vorgegeben, dass die gesuchte Lösung an den Außenrändern verschwinden soll. Dass die Näherungslösung dem gehorcht, ist auf der Grafik zu erkennen. Es handelt sich bei diesem Zylindermodell um ein zwar einfaches, aber reales physikalisches Problem aus der Neutronendiffusionstheorie für ein spezielles Kernreaktormodell eines BN-Reaktors. Dargestellt wird der schnelle Neutronenfluss.

In den Jahren von 1950 bis 1980 dominierte die FDM in den numerischen Programmen der Reaktorphysik zur Berechnung des Neutronenflusses. Verfahren dieser Art finden verbreitete Anwendung unter anderem bei fluiddynamischen Simulationen, zum Beispiel in der Meteorologie und der Astrophysik. Eine gewisse Verbreitung findet das Differenzenverfahren in der Baustatik. Schon 1904 analysierte Friedrich Bleich den Durchlaufträger; 1909 untersuchte Lewis Fry Richardson elastische Scheiben und 1919 Henri Marcus elastische Platten mit dem Differenzenverfahren.

Zu den Pionieren des Finite-Differenzen-Verfahrens für partielle Differentialgleichungen zählen Lewis Fry Richardson, Richard Southwell, Richard Courant, Kurt Friedrichs, Hans Lewy, Peter Lax und John von Neumann.

Geschichte des Namens

Bearbeiten

In keinem deutsch- oder englischsprachigen Lehrbuch bis etwa zum Jahr 1980 wird man die Termini Finite-Differenzen-Methode bzw. finite difference method finden. Zum Beispiel fehlt dieses Stichwort im Lexikon Mathematik von 1977, wo aber der Differenzenquotient enthalten ist.[3] Man sprach in der englischsprachigen Fachliteratur schon früher von finite-difference calculations, finite difference equations oder finite difference schemes, in der deutschsprachigen von Differenzenrechnungen, Differenzengleichungen und Diskretisierungsschemata. In der deutschsprachigen Fachliteratur tauchten Formen des englischen Adjektivs finite (begrenzt, endlich) in diesem Zusammenhang nicht auf. Streng genommen ist der Wortbestandteil finit überflüssig, denn Differenzen sind immer endlich. Sie sind das Gegenstück zum Differential, das einen unendlich kleinen Abschnitt auf der Achse eines Koordinatensystems beschreibt.[4] Auch handelt es sich nicht um eine spezielle Methode, sondern umfasst eine Vielzahl von Verfahren.

Finite difference method tritt erst in den 1980er Jahren auf, und zwar nach dem „Siegeszug“ der finite element method. Um 1960 wurde der Name finite element method für eine verwandte Methode eingeführt und mit Finite-Elemente-Methode eingedeutscht. Die wahrscheinlich erste Verwendung des Begriffs finite element findet sich 1960 in einem Zeitschriftenaufsatz von Ray W. Clough, einem Professor für Baustatik.[5] Diese Begriffe verbreiteten sich in den 1970er Jahren schnell.[6] Erst deutlich später muss den Mathematikern, Physikern und Ingenieuren aufgefallen sein, dass für die klassischen Diskretisierungsverfahren kein „griffiger“ Name existiert. Sie fanden, angelehnt an den Namen finite element method, dass der Name finite difference method recht passend sei. Ein Standardlehrbuch zu dieser Materie zeigt den Entwicklungsgang: Es trug 1965 noch den Titel Numerical solution of partial differential equations: With exercises and worked solutions.[7] Die Monographie wurde 1970 ins Deutsche übertragen.[8] Eine spätere, erweiterte Auflage der englischen Version hatte 1985 das finite difference methods dann schon im Untertitel.[9] Diese Auflage wurde wahrscheinlich nicht ins Deutsche übersetzt, so dass der Nachweis eines ersten Auftretens von Finite-Differenzen-Methode so nicht erbracht werden kann.

Es dauerte noch eine Weile, bis finite difference method in Finite-Differenzen-Methode eingedeutscht wurde. Man beachte einen kleinen Unterschied. Hätte man den englischen Namen direkt eingedeutscht, müsste es eigentlich Finite-Differenz-Methode heißen. In den 1970er, 1980er Jahren scheint die Beschäftigung mit numerischen Methoden, zumindest in den beiden deutschen Staaten, für Mathematiker kein attraktives Betätigungsfeld gewesen zu sein, wie man durch Recherchen in den Katalogen der deutschsprachigen Bibliotheksverbünde nach entsprechenden Fachbüchern nachweisen kann. Die Anwendung und Weiterentwicklung der Differenzenrechnung verlagerte sich auf Physiker und Ingenieure. Wer sich damit beschäftigte, bezog seine Fachliteratur ohnehin aus der englisch- und russischsprachigen Fachliteratur. In der russischsprachigen verwendet man heute Метод конечных разностей, also Methode der endlichen Differenzen. Eine Recherche in den Katalogen der deutschsprachigen Bibliotheksverbünde ergab, dass wahrscheinlich das erste Fachbuch in deutscher Sprache mit Finite-Differenzen-Methode im Titel im Jahr 1987 erschienen ist, das von Ingenieuren und nicht von Mathematikern verfasst wurde.[10] Es wird im Katalog der deutschen Nationalbibliothek nach den Sachgruppen Bergbau, Bautechnik, Umwelttechnik klassifiziert[11] und ist in vier Auflagen erschienen.[12]

Der Name klassisches Differenzenverfahren wäre sicher treffender gewesen.[13] Klassisch deshalb, weil die Methodik bereits seit Anfang des 19. Jahrhunderts bekannt ist, nachdem Augustin-Louis Cauchy der Differentialrechnung die heute übliche logische Strenge gegeben hatte. Er ging von den infinitesimalen Größen ab und definierte die Ableitung als Grenzwert von Sekantensteigungen.[14] Bereits Leibniz hatte das Wort Differentialquotient eingeführt, das den Grenzwert eines Differenzenquotienten benannte.[15] Richtig Fahrt nahm die Anwendung der Differenzenverfahren erst mit dem Einzug von Großrechnern in Forschungseinrichtungen auf. Das war Mitte der 1950er Jahre, als komplizierte Differentialgleichungen numerisch gelöst wurden. Motor der Differenzenrechnung und der Entwicklung iterativer Methoden zur Lösung der dabei entstehenden großen linearen Gleichungssysteme war die Reaktorphysik, insbesondere die Neutronendiffusionstheorie.[16]

Aus heutiger Sicht ist der Name Finite-Differenzen-Methode dennoch nicht so übel, denn es ist immer hilfreich, wenn sich Fachausdrücke in unterschiedlichen Sprachen ähneln. Zusammengefasst: Die Finite-Differenzen-Methode, die Finite-Elemente-Methode wie auch die Finite-Volumen-Verfahren sind Teilgebiete der Differenzenrechnung, der diskreten Entsprechung der Differentialrechnung.

FDM einer gewöhnlichen DGL 2. Ordnung

Bearbeiten

Wir suchen die Lösung einer Differentialgleichung 2. Ordnung[17] in einer Dimension (1D) in einem geradlinigen Koordinatensystem

 .

Zu berechnen sei die Funktion   im Intervall   der unabhängigen Variablen  ,  . Die sogenannte Inhomogenität   sei gegeben. Ohne weitere Bedingungen ist   nicht eindeutig bestimmt. Erst mit zwei Anfangs- oder Randbedingungen erhält man Eindeutigkeit. Nur im Fall von Randbedingungen favorisiert man zur numerischen Lösung finite Differenzen.

Das Vorgehen der numerischen Lösung nach der Differenzenmethode besteht darin, aus der Differentialgleichung ein Differenzengleichungssystem zu erzeugen und dieses zu lösen. Dazu hat der Anwender im ersten Schritt auf der  -Achse Stützpunkte (Gitterpunkte) zu wählen, über deren Lage er frei verfügen kann.

Äquidistante Stützstellen

Bearbeiten
 
Die FDM beruht auf der Diskretisierung einer Funktion auf einem Gitter, hier der Funktion einer unabhängigen Variablen mit äquidistanten Stützstellen

Wählt er ein äquidistantes Gitter, so muss er sich entscheiden, mit wie vielen Stützpunkten er starten möchte. Wählt er als Anzahl der „inneren“ Stützpunkte  , so sind   die Anzahl der Intervalle und   die Anzahl der Stützstellen insgesamt. Dann sind die Gitterschrittweite (Maschenweite, Intervallbreite)  

  und die Abszissenwerte   der Stützpunkte   für  .
 
Drei-Punkte-Stützstellenschema (englisch Three-Points stencil) für zentrale Differenzen in 1D mit äquidistanten Stützstellen

Im zweiten Schritt hat der Anwender eine Näherung für die zweite Ableitung zu wählen. Die entsprechenden Formeln kann er zum Beispiel aus dem Artikel Koeffizienten für Differenzenquotienten entnehmen. Die kleinst mögliche (und meist auch zu empfehlende) Anzahl von Gitterpunkten für eine Näherung der zweiten Ableitung ist drei. Für Stützstellen, die nicht auf dem Enden des Intervalls   oder   liegen, für die inneren Stützstellen, ist der zentrale Differenzenquotient die passende Wahl:

 

In der Grafik ist die rote die vorgegebene Stützstelle, die beiden Nachbarn sind die blauen Punkte.

Für die inneren Stützstellen folgen daraus die Differenzengleichungen

   für   

für die unbekannten numerischen Näherungswerte   von   und mit den vorgegebenen Werten  , den Werten der Inhomogenität   an den Stützstellen  .

Die Gesamtheit der Differenzengleichungen für innere Stützstellen bilden ein lineares Gleichungssystem der Ordnung   für die   unbekannten Funktionswerte von  .

Randwertproblem für äquidistante Stützstellen

Bearbeiten

Sind die Funktionswerte   und   vorgegeben, haben wir es mit einem Randwertproblem zu tun. Dann ist   eindeutig bestimmt.

Die beiden „Ränder“ sind die erste und die letzte Stützstelle, die das vorgegebene  -Intervall links und rechts begrenzen. Randwertprobleme nennt man in der Mathematik die wichtige Klasse von Problemstellungen, bei denen zu einer vorgegebenen Differentialgleichung Lösungen gesucht werden, die auf dem Rand des Definitionsbereiches vorgegebene Funktionswerte, die Randbedingungen, annehmen sollen.

Die erste Differenzengleichung wird mit den Randwerten dann zu

  oder  ,

und die letzte

  oder  .

In Matrixform lautet das zu lösende System, wenn wir beide Seiten aller Differenzengleichungen mit dem Faktor   multiplizieren, dann:

 

Die Matrix heißt Koeffizientenmatrix und ist vom Rang  . Sie ist regulär und somit ist das erzeugte Gleichungssystem eindeutig lösbar und ihre Inverse ist sogar formelmäßig angebbar. Da in jeder Zeile maximal nur drei Unbekannte vorkommen, gehört sie zu den dünnbesetzten Matrizen. Sie ist eine Bandmatrix, bei der zusätzlich zur Hauptdiagonalen nur eine obere und eine untere Nebendiagonale Elemente ungleich Null aufweisen und ist somit eine Tridiagonalmatrix. Da die Elemente der Hauptdiagonalen und der beiden Nebendiagonalen konstant sind, zählt sie darüber hinaus zu den Tridiagonal-Toeplitz-Matrizen.

Ähnlich spezielle Koeffizientenmatrizen treten auch anderweitig bei der Diskretisierung nach der Finiten-Differenzen-Methode auf. Solche speziellen Matrizentypen verringert die Anzahl der Rechenoperationen bei der Lösung des Gleichungssystems. Insbesondere aber spart man Arbeitsspeicher, um die Elemente der Koeffizientenmatrix zu speichern. Bei einer Anzahl der Unbekannten   und einer voll besetzen Koeffizientenmatrix würde man   Speicherplätze benötigen, dagegen sind es im Fall einer Tridiagonalmatrix nur  . In der Frühzeit der Anwendung der FDM, in den 1950er Jahren, musste man mit dem Speicherplatz sehr haushalten. Deshalb konnte man vorrangig nur 1D-Probleme behandeln, wobei die Anzahl der Unbekannten in realen Problemstellungen durchaus schon die Millionengrenze erreichte.

Beispiel für äquidistante Stützstellen

Bearbeiten
 
Exakte Lösungsfunktion des Beispiels und Werte der Näherungslösungen mit 2, 3 und 4 äquidistanten Intervallen

Gegeben sei das Randwertproblem

   für  

mit den Randbedingungen

 .

Die Inhomogenität   ist also konstant.

Die Lösung   lässt sich in diesem Fall exakt angeben, und zwar durch die Festlegung  . Dass das tatsächlich eine Lösung ist, kann man durch Einsetzen in die Differentialgleichung leicht nachweisen. Sie erfüllt auch die Randbedingungen.

Wählen wir vier Intervalle, berechnen die Näherungswerte von   an drei inneren Stützstellen, ist die Gitterschrittweite  . Es ergibt sich das Gleichungssystem in Matrixform

 

mit der Lösung

 .

Die Grafik zeigt den Funktionsverlauf der exakten Lösung des Beispiels und Werte der Näherungslösungen mit 2, 3 und 4 Intervallen. Die diskreten Näherungslösungen an den Gitterpunkten sind durch farbige Kreise markiert, die durch Geraden in gleicher Farbe verbunden sind. Es ist anzumerken, dass alle Näherungslösungen dieses Beispiels Funktionswerte liefern, die genau gleich den Funktionswerten der exakten Lösung an den Stützstellen sind. Das liegt daran, dass für ein quadratisches Polynom als exakte Lösung der Konsistenzfehler bei dieser Diskretisierung gleich Null ist und damit auch der Fehler in den Gitterpunkten. Dies ist natürlich ein Ausnahmefall.

Nichtäquidistante Stützstellen

Bearbeiten
 
Drei-Punkte-Stützstellenschema für zentrale Differenzen in 1D mit nichtäquidistanten Stützstellen

Nichtäquidistante Stützstellen sind angebracht, wenn die Lösung Besonderheiten aufweist, z. B. Singularitäten oder Grenzschichten. Wir betrachten wieder die Differentialgleichung

 

mit der vorgegebenen Inhomogenität   und mit gegebenenfalls vorgegebenen Randwerten   und  

Die unabhängige Variable   werde mit folgenden beliebigen, nichtäquidistanten Stützstellen diskretisiert:

 

mit   und  

Da die Nummerierungen   der Gitterschrittweiten   von Autor zu Autor unterschiedlich sind, könnte folgende anschauliche „Merkregel“ nützlich sein: Betrachten wir den Gitterpunkt  , so bedeuten bei der Nummerierung in diesem Abschnitt   die Gitterschrittweite vor und   die Gitterschrittweite nach dem Gitterpunkt.

Es ist folgende zentrale Drei-Punkte-Diskretisierung zweckmäßig:

    für   

Ordnen wir die Formel wieder in der Reihenfolge aufsteigender Stützstellen und multiplizieren beide Seiten mit  , so ergibt sich

    für    

oder

    für    

Die sich daraus ergebende Koeffizientenmatrix ist wie im Fall äquidistanter Stützstellen eine Tridiagonalmatrix. Da die Elemente auf der Hauptdiagonalen und der beiden Nebendiagonalen nicht konstant sind, ist es keine Matrix vom Tridiagonal-Toeplitz-Typ.

Setzen wir   für alle  , dann ergibt sich wieder die Formel für äquidistante Stützstellen

    für    

Zur Untersuchung des Konsistenzfehlers wird die exakte Lösung in die Diskretisierung eingesetzt und man hat folgenden Term abzuschätzen:

 

Taylorentwicklung liefert

 

Auf äquidistanten Gittern verschwindet der erste Anteil und der Konsistenzfehler ist von der Ordnung 2. Auf beliebigen Gittern jedoch reduziert sich der Konsistenzfehler auf die Ordnung Eins! Trotzdem ist der Fehler selbst in der Maximumnorm wieder von der Ordnung Zwei. Das kann man mit Hilfe der diskreten Greenschen Funktion oder verbesserten Stabilitätsaussagen verifizieren. Es gilt also

 

Randwertproblem für nichtäquidistante Stützstellen

Bearbeiten

Sind wie oben wieder die Randwerte   und

  gegeben, wird die erste Differenzengleichung zu

 

und die letzte

 

Ansonsten gilt auch im Fall nichtäquidistanter Stützstellen das im Abschnitt Randwertproblem für äquidistante Stützstellen gesagte.

Beispiel für nichtäquidistante Stützstellen

Bearbeiten
 
Exakte Lösungsfunktion des Beispiels und Werte der Näherungslösungen mit 2, 3 und 4 nichtäquidistanten Intervallen

Wie im Fall äquidistanter Intervalle soll wieder folgendes Randwertproblem numerisch gelöst werden

   für  

mit den Randbedingungen

 .

Die Inhomogenität   ist also wieder konstant.

Die Lösungsfunktion   kennen wir schon. Sie ist gegeben durch  .

Bevor wir die Intervalle wählen, eine Vorbemerkung zu physikalischen Größen. Das ist wichtig für Beispiele, wo Zahlenwerte für Größen eingesetzt werden. In mathematischen Lehrbüchern findet man zum Beispiel solche Texte wie: Die Gitterschrittweite   beträgt  . Das ist, vom physikalischen Standpunkt (und mit Differentialgleichungen werden bevorzugt physikalische Probleme gelöst), nicht korrekt. Die Gitterschrittweite ist eine physikalische Größe und besteht immer aus dem Produkt einer Zahl und einer Einheit, also zum Beispiel  . Darauf sollten auch Mathematiker hinweisen: Die Gitterschrittweite beträgt   (mit welcher Längeneinheit auch immer) oder Der Zahlenwert der Gitterschrittweite beträgt   Die nachfolgend verwendeten Symbole   stehen nicht für physikalische Größen, sondern für Zahlenwerte der physikalischen Größen.

Behandeln wir für dieses Beispiel wieder den Fall mit vier Intervallen, berechnen folglich die Näherungswerte der Funktion   an drei inneren Stützstellen. Als Zahlenwerte der nicht gleich breiten Gitterschrittweiten sollen vorgegeben werden

 

Es ergibt sich das Gleichungssystem in Matrixform

 

mit der Lösung

 

Die Grafik zeigt den Funktionsverlauf der exakten Lösung des Beispiels und Werte der Näherungslösungen mit zwei, drei und vier Intervallen. Bei zwei Intervallen wurden die Schrittweiten

 

gewählt, bei drei Intervallen

 

Die diskreten Näherungslösungen an den Gitterpunkten der Grafik sind wieder durch farbige Kreise markiert, die durch Geraden in gleicher Farbe verbunden sind. Ergebnis: Für diese Differentialgleichung liefern alle Näherungslösungen in Fall nichtäquisistanter Stützstellen ebenfalls Werte, die genau gleich den Funktionswerten der exakten Lösung an den Stützstellen sind, wie auch bei äquidistanten Gitterschrittweiten. Das war nach der oben gegebenen Formel für den Konsistenzfehler   auch zu erwarten.

Upwind finite Differenzen für ein Konvektion-Diffusionsproblem

Bearbeiten

Betrachtet wird das Konvektions-Diffusionsproblen

 

mit einem kleinen Parameter  . Durch den kleinen Parameter dominiert in diesem Problem der Konvektionsterm   und nicht wie oft üblich der Diffusionsterm  . Das normale Differenzenverfahren mit der äquidistanten Schrittweite   lautet

 

Mit   und   kann man das diskrete Problem exakt lösen und erhält mit  

 

Wenn   extrem klein ist, verbietet es sich, mit Schrittweiten zu arbeiten, die noch kleiner als   sind, insbesondere für mehrdimensionale Probleme. Ist aber  , so ist   negativ, und die numerische Lösung oszilliert und approximiert die exakte Lösung extrem schlecht!

Deshalb ist es üblich, eine upwind Finite-Differenzen-Methode einzusetzen. Dabei wird   statt durch den obigen zentralen Differenzenquotienten durch einen einseitigen ersetzt, also durch   oder  . Dies geschieht in Abhängigkeit von der Stromrichtung (deshalb upwind), hier in Abhängigkeit vom Vorzeichen von  . Im eindimensionalen Fall wird so erreicht, dass die Koeffizientenmatrix des diskreten Problems positive Diagonalelemente und negative Außerdiagonalelemente besitzt (dies garantiert Stabilität des diskreten Problems). Das upwind-Verfahren lautet also

 

es liefert eine robuste Approximation der Lösung auch im Fall  .

Beispiel zur numerischen Lösung der Wärmeleitungsgleichung

Bearbeiten

Im Folgenden wird die numerische Lösung der Wärmeleitungsgleichung auf einem beschränkten Gebiet   betrachtet:

 
 
 

Numerische Lösung im 1D

Bearbeiten

Im 1D-Fall ist   ein beschränktes Intervall. Da in diesem Fall nur eine Ortsableitung betrachtet wird, kann die Wärmeleitungsgleichung folgendermaßen geschrieben werden:

 

Diskretisierung

Bearbeiten

Um die Finite-Differenzen-Methode anwenden zu können, muss das Intervall   zunächst in endlich viele Teilintervalle unterteilt werden. Hierfür werden   äquidistante Stützstellen verwendet:

 , für  .

Die Gitterweite dieser Diskretisierung ist also  . Nach Voraussetzung verschwindet die gesuchte Funktion   an den Randwerten, d. h.  , sodass diese Werte nicht weiter betrachtet werden müssen. Damit lassen sich die Funktionsauswertungen von   an Stützstellen als Vektor im   darstellen:

 

Approximation der Ableitung

Bearbeiten

Die zweite Ableitung von   bzgl. des Orts kann nun an den Stützstellen durch Differenzenquotienten zweiter Ordnung approximiert werden:

 

Wird die Wärmeleitungsgleichung nach   umgestellt, ergibt sich damit folgendes System gewöhnlicher Differentialgleichungen erster Ordnung:

 

wobei   und  .

Dieses System kann nun durch beliebige Verfahren zur Lösung von gewöhnlichen Differentialgleichungen, wie z. B. das Runge-Kutta-Verfahren oder das Euler-Verfahren, gelöst werden.

Beispiel zur numerischen Lösung der zweidimensionalen Poissongleichung

Bearbeiten

Betrachtet wird exemplarisch

 

Auf dem Gitter mit den Gitterpunkten   mit   mit   und   wird eine Näherung   für   gesucht.

Analog wie im eindimensionalen Fall diskretisiert man die zweiten Ableitungen in jeder Koordinatenrichtung in den inneren Gitterpunkten und erhält für   das diskrete Problem

 

Dabei setzt man natürlich die Werte auf dem Rand   gleich Null. Diese Diskretisierung heißt Fünf-Punkte-Stern für die Poissongleichung. Man kann zeigen, dass das erzeugte Gleichungssystem eine eindeutige Lösung besitzt. Zudem gilt für den Fehler

 

falls die exakte Lösung des Problems ausreichend glatt ist. Dies ist im zweidimensionalen Fall insbesondere in Gebieten mit Ecken keinesfalls trivial.

Konsistenz, Stabilität und Konvergenz

Bearbeiten

Eine Finite-Differenzen-Methode erzeugt ein lineares Gleichungssystem (analog Gleichung im Kapitel Beispiel)

 

wobei   die numerische Approximation der Lösung ist und   die Abhängigkeit vom Gitter explizit darstellen soll. Sei   die exakte Lösung und   die endliche Darstellung mittels  . Sowohl   als auch   sind sogenannte Gitterfunktionen, sie sind nur definiert in allen Gitterpunkten des verwendeten Gitters.

Eine FDM heißt konsistent von Ordnung  , falls es ein   gibt, so dass man den Konsistenzfehler folgendermaßen abschätzen kann:

 

Dabei ist   eine Norm für Gitterfunktionen. Oft verwendet wird die Maximumnorm  , sie erlaubt punktweise Fehlerabschätzungen.

Eine FDM heißt stabil, falls es ein   gibt, sodass für alle Gitterfunktionen   gilt

 

Aus Konsistenz der Ordnung   und Stabilität folgt Konvergenz der Ordnung  

 

Stabilität vorausgesetzt, ist Konsistenz der Ordnung   für Konvergenz der Ordnung   hinreichend, aber nicht notwendig.

Klassen von partiellen Differentialgleichungen zweiter Ordnung

Bearbeiten

Betrachtet werden folgende partielle Differentialgleichung zweiter Ordnung:

 

Gesucht ist die Funktion  , von der oft noch Anfangs- oder Randbedingungen vorgeschrieben werden.

Die Funktionen  ,  ,  ,  ,   und   sind gegebene Koeffizienten und sollten daher beschränkt und im gegebenen Integrationsbereich stetig sein. Bei vielen Gleichungen, die spezielle technische Probleme beschreiben, nehmen diese Funktionen konstante Werte an. Diese Gleichungen vom werden oft nach dem Wert der Diskriminante   in eine der folgenden Gruppen eingeteilt: hyperbolisch für  , parabolisch für   und elliptisch für  .

Typische Zusatzbedingungen für elliptische Probleme sind Randbedingungen, für parabolische und hyperbolische Probleme jedoch Anfangsbedingungen und Randbedingungen.

Beispiele für partielle Differentialgleichungen, die für physikalische Probleme formuliert werden, sind:

 
 
 
 
 

Literatur

Bearbeiten
  • Clarence Hudson Richardson: An introduction to the calculus of finite differences. Van Nostrand, New York 1954 (142 S.).
  • Christian Großmann, Hans-Görg Roos: Numerische Behandlung partieller Differentialgleichungen. 3., völlig überarb. und erw. Auflage. B. G. Teubner, Wiesbaden 2005, ISBN 978-3-519-22089-3 (572 S., eingeschränkte Vorschau in der Google-Buchsuche).}
  • Stig Larsson, Vidar Thomée: Partielle Differentialgleichungen und numerische Methoden. Springer, Berlin, Heidelberg, New York 2005, ISBN 978-3-540-20823-5 (272 S., eingeschränkte Vorschau in der Google-Buchsuche).
  • Claus-Dieter Munz, Thomas Westermann: Numerische Behandlung gewöhnlicher und partieller Differenzialgleichungen: Ein anwendungsorientiertes Lehrbuch für Ingenieure. 4., verbesserte und überarbeitete Auflage. Springer Vieweg, Berlin 2019, ISBN 978-3-662-55886-7 (394 S., eingeschränkte Vorschau in der Google-Buchsuche).
Bearbeiten
Commons: Finite-Differenzen-Methode – Sammlung von Bildern, Videos und Audiodateien

Siehe auch

Bearbeiten

Einzelnachweise

Bearbeiten
  1. Modell eines Kernreaktors, berechnet mit dem Programm MCT (Christian Reiche: Das zweidimensionale Mehrgruppenprogramm MCT und einige Untersuchungen über den Einfluss des Axialblankets bei schnellen Reaktoren. In: Kernenergie. Band 9, 1966, S. 232–240 (Artikeldaten online).), Ergebnisse als Grafik dargestellt mit dem Tabellenkalkulationsprogramm Microsoft Excel 2013, und zwar mit der Diagramm-Option 3D-Säulen.
  2. Eine Zelle ist hier in 2D eine von vier benachbarten Gitterpunkten festgelegte Fläche, in 3D ein von acht benachbarten Gitterpunkten festgelegter Raumbereich.
  3. Walter Gellert, Herbert Kästner, Siegfried Neubert (Hrsg.): Lexikon der Mathematik. 1. Auflage. Bibliographisches Institut, Leipzig 1977 (624 S.).
  4. Jürgen Schmidt: Basiswissen Mathematik. Springer-Verlag, 2014, ISBN 978-3-662-43545-8, S. 307 (eingeschränkte Vorschau in der Google-Buchsuche).
  5. Ray W. Clough: The Finite Element Method in Plane Stress Analysis. American Society of Civil Engineers, 1960 (35 S.).
  6. Finite element method. In: Books Ngram Viewer. Abgerufen am 15. November 2022.
  7. Gordon D. Smith: Numerical solution of partial differential equations: With exercises and worked solutions. Oxford University Press, London 1965 (179 S.).
  8. Gordon D. Smith: Numerische Lösung von partiellen Differentialgleichungen: Mit Aufgaben und Lösungen. Studienausgabe Auflage. Vieweg; C. F. Winter, Braunschweig, Basel 1970, ISBN 3-528-08296-8 (246 S.).
  9. Gordon D. Smith: Numerical solution of partial differential equations: Finite difference methods. Clarendon, Oxford 1985, ISBN 0-19-859650-2 (337 S.).
  10. Herbert Fritsch, Rolf Wohlfahrt, Manfred Teige (Hrsg.): Finite-Differenzen-Methode. 1. Auflage. IRB-Verlag, Stuttgart 1987, ISBN 3-8167-1279-7 (75 S.).
  11. Finite-Differenzen-Methode / [Hrsg.: Informationszentrum Raum u. Bau. d. Fraunhofer-Ges. (IRB). Red. Bearb.: Manfred Teige], Belegexemplar DNB 942077040 bei der Deutschen Nationalbibliothek.
  12. Herbert Fritsch, Rolf Wohlfahrt, Manfred Teige (Hrsg.): Finite-Differenzen-Methode. 4., erweiterte Auflage. IRB-Verlag, Stuttgart 1994, ISBN 3-8167-1279-7 (121 S.).
  13. Lothar Collatz: Eigenwertaufgaben mit technischen Anwendungen. 2., durchges. Auflage. Akad. Verl.-Ges. Geest & Portig, Leipzig 1963, S. 354 (XIV, 500).
  14. Thomas Sonar: 3000 Jahre Analysis, Springer, S. 506–514.
  15. Richard Courant: Vorlesungen über Differential- und Integralrechnung: Erster Band, Funktionen einer Veränderlichen. Zweite, verbesserte Auflage. Julius Springer, Berlin 1930, S. 79 (410 S.).
  16. Eugene Leon Wachspress: Iterative solution of elliptic systems and applications to the neutron diffusion equations of reactor physics. Prentice-Hall, Englewood Cliffs, N. J. 1966 (ix, 299 S.).
  17. Die höchste vorkommende Ableitungsordnung   wird Ordnung der Differentialgleichung genannt.