Gegeben ist die Poisson-Gleichung
Es soll ein Gitter mit den Punkten verwendet werden. Für die Neumann-Randbedingung soll die Ableitung durch einen zentralen Differenzenquotienten approximiert werden. Dazu sind zusätzliche Punkte
einzuführen.
Stellen Sie das Gleichungssystem auf, das bei der Diskretisierung mittels finiter Differenzen entsteht!
Lösung
Wie in der Aufgabenstellung gefordert, zerlegen wir das Einheitsquadrat in neun gleiche Felder. Die gesuchte Funktion ist auf den Randpunkten bis
bereits bekannt. Auf den anderen beiden Randpunkten kennen wir nur die Ableitung:
Den Laplace-Operator modellieren wir mit Hilfe der Finite-Differenzen-Methode, indem wir in beiden Raumdimensionen durch Differenzenquotienten approximieren und dabei die oben, unten, links und rechts an den jeweiligen Knoten ankreuzenden Nachbarn nutzen:
Es ergibt sich:
Beispielsweise lautet der Fünf-Punkte-Stern für den Knoten 5:
Die 5-Punkte-Sterne für die Knoten 2, 3 und 6 ergeben sich analog. Problematisch wird es bei den Knoten 1 und 4, da hier keine linken Nachbarn vorhanden sind. An diesen Punkten ist aber die Ableitung gegeben. Wir können daher mit der zentralen Differenz die Hilfspunkte 0 und -1 erschaffen:
Die Differenzen sind hier so herum, da die Normalenableitung am linken Rand nach außen (links) zeigt, statt in Koordinatenrichtung (rechts). Mit den Hilfspunkten ergeben sich die fehlenden Differenzensterne:
Knoten 1:
Knoten 4:
Es ergibt sich das folgende Gleichungssystem: