Newton-Iteration für nichtlineare PDE

8

Ich habe Schwierigkeiten zu verstehen, wie man die Newton-Iteration auf nichtlineare PDEs anwendet und dann ein vollständig implizites Schema für den Zeitschritt verwendet. Zum Beispiel möchte ich die Burgers-Gleichung lösen

ut+uuxuxx=0

Also Diskretisierung der Zeit mit einem Backward Euler

ut=un+1unh

wir glauben, dass

un+1- -unh+un+1(un+1)x- -(un+1)xx=0un+1- -h(un+1)xx+hun+1(un+1)x=un(ich- -hD.2)un+1+hun+1D.un+1=un(ich- -hD.2)un+1+N.(un+1)=un   (1)

Dabei steht für unseren nichtlinearen Term (beachten Sie, dass der nichtlineare Term implizit geschrieben ist). Jetzt wollen wir die Newton-Iteration auf diese nichtlineare ODE anwenden, aber hier stecke ich fest:N.

  1. Wenden wir einfach die Newton-Iteration auf die LHS von und ignorieren dabei den Term u n , dh lösen ( I - h D 2 ) u n + 1(1)un ? Oder sollen wir denBegriff u n aufnehmen ? (Nur zur Erinnerung, ich möchte nach Verwendung der Newton-Iteration einen Zeitschritt mit einem vollständig impliziten Schema ausführen, daher glaube ich, dass wir nur die LHS = 0 lösen möchten.)(ich- -hD.2)un+1+N.(un+1)=0un

  2. Was sollen wir dann mit den Informationen aus der ersten Vermutung und dem Ergebnis unserer Newton-Iteration tun? Wie verwenden wir diese Informationen in unserem Zeitschritt?

Da ich sicher bin, dass es schmerzlich offensichtlich ist, bin ich ziemlich verwirrt, wie ich dieses Problem angehen soll. Wenn jemand eine detaillierte Beschreibung der Anwendung der Newton-Iteration und des Zeitschritts auf nichtlineare PDEs (jedoch keine elliptischen PDEs) geben oder mir bei dem vorliegenden Problem helfen könnte, wäre ich sehr dankbar. Danke im Voraus.

Mattos
quelle

Antworten:

7

Es ist etwas einfacher zu sehen, ob Sie Ihre Gleichung in einem halbdiskretisierten System der Form und unter Anwendung der θ- Methode und unter Annäherung an u ' ( t ) schreiben. ( w n + 1 - w nu(t)=F(u(t))θ ergibt dies,u(t)(wn+1wn)/τ

wn+1- -wn- -(1- -θ)τF.(wn)- -θτF.(wn+1)=0

mit unbekanntem Vektor und Zeitschritt τ . Hier ist u ' ( t ) die nicht diskretisierte zeitliche partielle Ableitung und F ( u ( t ) ) repräsentiert Ihre diskretisierten räumlichen Ableitungen, die zum Zeitpunkt t zu u ( t ) ausgewertet werden . Mit dem θwn+1τu'(t)F.(u(t))u(t)tθ Methode bietet Ihnen ein wenig Flexibilität, da Sie die Methode zwischen vollständig implizit und vollständig explizit (und irgendwo dazwischen) ändern können.

Diese Gleichung kann mit der Newton-Iteration gelöst werden.

νk+1=νk- -(ich- -θτEINn)- -1(νk- -wn- -(1- -θ)τF.(wn)- -θτF.(wn+1))

Dabei ist der Iterationsindex ( k 0 ) und A n die Jacobi-Matrix von F ( w n)kk0EINn . Wir verwenden das Symbol ν k für Iterationsvariablenso dass sie ausLösung der Gleichung bei einer Echtzeit Punkt auszeichnen u n . Wie Sie sagten, benötigt die Iteration einen Anfangswert, es ist vollkommen gültig, ν 0 = w 0 zu wählen,oder für eine bessere Schätzung können wir eine Iteration ν 0 = w 0 + τ F ( w) vorberechnenF.(wn)νkunν0=w0. Genau genommen ist dies die sogenanntemodifizierteNewton-Iteration, da der Jacobian während der Iteration nicht aktualisiert wird, was bekanntermaßen für steife PDEs gut funktioniert.ν0=w0+τF.(w0)

Boyfarrell
quelle
Vielen Dank für Ihre Antwort. Ich hätte genauer sein sollen, ich bin tatsächlich schon am Punkt der Verwendung von Newton Iteration und habe das Problem als solches gelöst und habe Diagramme der Lösung. Aber mein Hauptproblem ist, wie man danach einen Zeitschritt macht? Ich dachte, wir wollten den nichtlinearen Term "vernichten" und danach eine lineare ODE zeitlich verschieben? Oder habe ich den Punkt der Newton-Iteration völlig verfehlt? Könnten Sie bitte auch eine schrittweise Methode angeben, wie Sie zu Ihrer Newton-Iterationsgleichung gekommen sind? Ich entschuldige mich dafür, dass ich so viele Fragen gestellt habe, dass ich im Moment nur ein bisschen Schwierigkeiten habe.
Mattos
Ich denke, Sie haben den Punkt verpasst. Sie diskretisieren Ihre nichtlinearen Gleichungen und iterieren dann (beginnend mit einem Schätzwert), um die Gleichungen zu lösen und zum Zeitpunkt t 1 zu ergebenw1t1νt2
θ=0
Weitere Informationen finden Sie auf S.127 books.google.co.uk/books?isbn=3540034404 . Ich habe damit numerische Algorithmen implementiert und es funktioniert.
Boyfarrell