www.wikidata.de-de.nina.az

Von Neumann Stabilitätsanalyse Die nach John von Neumann manchmal auch L2 Stabilitätsanalyse ist das Standardverfahren z

Von-Neumann-Stabilitätsanalyse

Die Von-Neumann-Stabilitätsanalyse (nach John von Neumann), manchmal auch L2-Stabilitätsanalyse, ist das Standardverfahren zur Untersuchung der Stabilität von numerischen Verfahren zur Lösung zeitabhängiger partieller Differentialgleichungen.

Das Verfahren wurde von John von Neumann in Los Alamos im Rahmen des Manhattan-Projekts entwickelt. Während des Krieges wurde die Methode unter Verschluss gehalten und erst 1947 von John Crank und Phyllis Nicolson publiziert. 1968 bewies Heinz-Otto Kreiss weitere zentrale Eigenschaften des Analyseverfahrens.

Der lineare, eindimensionale Fall Bearbeiten

Gegeben sei auf einem Intervall eine lineare, partielle Differentialgleichung mit konstanten Koeffizienten der Form

Anfangsdaten sowie ein numerisches Verfahren zur Lösung. Die Bedingung, dass das Verfahren in der L2-Norm stabil ist, besagt dann, dass der durch das numerische Verfahren produzierte Fehler bei gegebenen Schrittweiten und für beschränkt bleibt. Der erste Schritt der Von-Neumann-Stabilitätsanalyse besteht nun darin, die Lösung periodisch auf die kompletten reellen Zahlen fortzusetzen.

Der periodische Fehler zum Zeitpunkt im Diskretisierungspunkt kann nun in eine Fourier-Reihe

entwickelt werden. Hierbei bezeichnet die imaginäre Einheit. Das numerische Verfahren definiert dann eine Evolution der Koeffizienten der Fourierreihe in der Zeit mittels einer so genannten Amplifikationsmatrix . Die L2-Stabilitätsbedingung reduziert sich dann darauf, dass das numerische Verfahren genau dann stabil ist, wenn der Spektralradius der Amplifikationsmatrix betragsmäßig kleiner gleich eins ist.

Beispiel Bearbeiten

Der einfachste Fall ist die lineare Advektionsgleichung

wobei eine reelle Zahl ist. Eines der einfachsten vorstellbaren numerischen Verfahren zur Lösung solcher Gleichungen ist das explizite Euler-Verfahren zur Zeitintegration gekoppelt mit zentralen Differenzen auf einem äquidistanten Gitter im Raum. Der zweite Term wird also mittels

approximiert. Insgesamt ergibt sich das Verfahren

was auch die Entwicklung der Fehler definiert und auch jedes einzelnen Terms der Fourier-Reihenentwicklung. Betrachten wir den j-ten Summanden, so ergibt Einsetzen in die obige Formel und Division durch mit :

Die Amplifikationsmatrix ist nun gegeben durch

Das Verfahren ist -stabil, falls für alle , was hier nicht der Fall ist, da Damit ist das Verfahren unabhängig von der Wahl der Schrittweiten instabil. Dieses Verhalten beobachteten die Mitarbeiter des Manhattan Projekts, was von Neumann zur Entwicklung der Stabilitätsanalyse führte. Wird im Raum das Upwind-Verfahren

benutzt, so ergibt sich die Courant-Friedrichs-Lewy-Bedingung

also bedingte Stabilität.

Andere Gleichungen Bearbeiten

Im nichtlinearen Fall oder im Falle variabler Koeffizienten kann das Verfahren durch Linearisierung und Einfrieren der Koeffizienten angewandt werden, allerdings liefert die Analyse im allgemeinen Fall nur noch eine notwendige Bedingung für Stabilität, in Spezialfällen auch eine hinreichende. Ferner ist die Bedingung an den Spektralradius nun:

Ein allgemeines Verfahren zur vollständigen Stabilitätsanalyse nichtlinearer Gleichungen ist nicht bekannt.

Weblinks Bearbeiten

wikipedia, wiki, deutsches, deutschland, buch, bücher, bibliothek artikel lesen, herunterladen kostenlos kostenloser herunterladen, MP3, Video, MP4, 3GP, JPG, JPEG, GIF, PNG, Bild, Musik, Lied, Film, Buch, Spiel, Spiele
Veröffentlichungsdatum:
Die Von Neumann Stabilitatsanalyse nach John von Neumann manchmal auch L2 Stabilitatsanalyse ist das Standardverfahren zur Untersuchung der Stabilitat von numerischen Verfahren zur Losung zeitabhangiger partieller Differentialgleichungen Das Verfahren wurde von John von Neumann in Los Alamos im Rahmen des Manhattan Projekts entwickelt Wahrend des Krieges wurde die Methode unter Verschluss gehalten und erst 1947 von John Crank und Phyllis Nicolson publiziert 1968 bewies Heinz Otto Kreiss weitere zentrale Eigenschaften des Analyseverfahrens Inhaltsverzeichnis 1 Der lineare eindimensionale Fall 1 1 Beispiel 2 Andere Gleichungen 3 WeblinksDer lineare eindimensionale Fall BearbeitenGegeben sei auf einem Intervall 0 L displaystyle 0 L nbsp eine lineare partielle Differentialgleichung mit konstanten Koeffizienten der Form u t A u x 0 displaystyle u t Au x 0 nbsp Anfangsdaten u 0 u 0 displaystyle u 0 u 0 nbsp sowie ein numerisches Verfahren zur Losung Die Bedingung dass das Verfahren in der L2 Norm stabil ist besagt dann dass der durch das numerische Verfahren produzierte Fehler bei gegebenen Schrittweiten D x displaystyle Delta x nbsp und D t displaystyle Delta t nbsp fur t displaystyle t to infty nbsp beschrankt bleibt Der erste Schritt der Von Neumann Stabilitatsanalyse besteht nun darin die Losung periodisch auf die kompletten reellen Zahlen fortzusetzen Der periodische Fehler E i n displaystyle E i n nbsp zum Zeitpunkt t n displaystyle t n nbsp im Diskretisierungspunkt x i displaystyle x i nbsp kann nun in eine Fourier Reihe E i n N N c j n e i j D x I displaystyle E i n sum N N c j n e ij Delta xI nbsp entwickelt werden Hierbei bezeichnet I displaystyle I nbsp die imaginare Einheit Das numerische Verfahren definiert dann eine Evolution der Koeffizienten der Fourierreihe in der Zeit mittels einer so genannten Amplifikationsmatrix G displaystyle G nbsp Die L2 Stabilitatsbedingung reduziert sich dann darauf dass das numerische Verfahren genau dann stabil ist wenn der Spektralradius der Amplifikationsmatrix r G displaystyle rho G nbsp betragsmassig kleiner gleich eins ist Beispiel Bearbeiten Der einfachste Fall ist die lineare Advektionsgleichung u t a u x 0 displaystyle u t au x 0 nbsp wobei a displaystyle a nbsp eine reelle Zahl ist Eines der einfachsten vorstellbaren numerischen Verfahren zur Losung solcher Gleichungen ist das explizite Euler Verfahren u n 1 u n D t a u x displaystyle u n 1 u n Delta t au x nbsp zur Zeitintegration gekoppelt mit zentralen Differenzen auf einem aquidistanten Gitter im Raum Der zweite Term wird also mittels a u x a 2 D x u i 1 n u i 1 n displaystyle au x approx frac a 2 Delta x u i 1 n u i 1 n nbsp approximiert Insgesamt ergibt sich das Verfahren u i n 1 u i n D t a 2 D x u i 1 n u i 1 n 0 displaystyle frac u i n 1 u i n Delta t frac a 2 Delta x u i 1 n u i 1 n 0 nbsp was auch die Entwicklung der Fehler definiert und auch jedes einzelnen Terms der Fourier Reihenentwicklung Betrachten wir den j ten Summanden so ergibt Einsetzen in die obige Formel und Division durch e i j D x I displaystyle e ij Delta xI nbsp mit ϕ j D x displaystyle phi j Delta x nbsp c j n 1 c j n a D t 2 D x c j n e I ϕ e I ϕ 0 displaystyle c j n 1 c j n frac a Delta t 2 Delta x c j n e I phi e I phi 0 nbsp Die Amplifikationsmatrix ist nun gegeben durch G ϕ c j n 1 c j n 1 I a D t D x sin ϕ displaystyle G phi frac c j n 1 c j n 1 I frac a Delta t Delta x sin phi nbsp Das Verfahren ist L 2 displaystyle L 2 nbsp stabil falls G 1 displaystyle G leq 1 nbsp fur alle ϕ displaystyle phi nbsp was hier nicht der Fall ist da G 2 1 a D t D x 2 sin 2 ϕ displaystyle G 2 1 left frac a Delta t Delta x right 2 sin 2 phi nbsp Damit ist das Verfahren unabhangig von der Wahl der Schrittweiten instabil Dieses Verhalten beobachteten die Mitarbeiter des Manhattan Projekts was von Neumann zur Entwicklung der Stabilitatsanalyse fuhrte Wird im Raum das Upwind Verfahren a u x a D x u i n u i 1 n displaystyle au x approx frac a Delta x u i n u i 1 n nbsp benutzt so ergibt sich die Courant Friedrichs Lewy Bedingung a D t D x lt 1 displaystyle frac a Delta t Delta x lt 1 nbsp also bedingte Stabilitat Andere Gleichungen BearbeitenIm nichtlinearen Fall oder im Falle variabler Koeffizienten kann das Verfahren durch Linearisierung und Einfrieren der Koeffizienten angewandt werden allerdings liefert die Analyse im allgemeinen Fall nur noch eine notwendige Bedingung fur Stabilitat in Spezialfallen auch eine hinreichende Ferner ist die Bedingung an den Spektralradius nun r G 1 c D t displaystyle rho G leq 1 c Delta t nbsp Ein allgemeines Verfahren zur vollstandigen Stabilitatsanalyse nichtlinearer Gleichungen ist nicht bekannt Weblinks BearbeitenMartin Neumann Erklarung der Analyse in einem Skript Computational Physics I Grundlagen der Uni Wien Abgerufen von https de wikipedia org w index php title Von Neumann Stabilitatsanalyse amp oldid 180572590