Der Effekt der Entkopplung eines gekoppelten Systems von PDEs

8

Ich habe vorher eine etwas ähnliche Frage gestellt, aber vielleicht war sie zu spezifisch, als dass jemand sie wirklich beantworten könnte. Hier ist eine etwas allgemeinere Frage, mit der ich zu kämpfen habe. Betrachten Sie das folgende System:

u 2

(D1(u2)u1)=f1(u2)
u2t+f2(u1,u2)(D2(u2)u2)=0

unter der Annahme einer allgemeinen Menge von BCs: D iu in = u i , N ,

ui=ui,D,onΓD
Diuin=ui,N,onΓN

Verwendung von DGFEM für räumliche Diskretisierungen und Rückwärts-Euler für die Zeitableitung. Wir entkoppeln uns so:

  1. Löse nach mit u k 2 : - ( D 1 ( u k 2 ) u k + 1 1 ) = f 1 ( u k 2 )u1k+1u2k

    (D1(u2k)u1k+1)=f1(u2k)
  2. Löse nach mit u k + 1 1 : u 2u2k+1u1k+1

    u2t+f2(u1k+1,u2k,u2k+1)(D2(u2k+1)u2k+1)=0

Die Newtonsche Methode wird verwendet, um die Nichtlinearität zu behandeln.

f2(u1,u2)=u2(u1+u2u2)

damit

f2(u1k+1,u2k,u2k+1)=u2k+1(u1k+1+u2ku2k)

f2u1k+1u2ku2k+1f2u2ku1k+1+u2ku2k

f2=f2(u1,exactk+1,u2k+1,u2,exactk+1)=u2k+1(u1,exactk+1+u2,exactk+1u2,exactk+1)

Wird diese suboptimale Leistung im ersten Szenario angesichts der Entkopplung erwartet? Ich weiß, dass die Entkopplung den Zeitschritt begrenzt, aber ich bin mir ziemlich sicher, dass ich den Zeitschritt ausreichend klein nehme.

f2u2u2u1k+1

Ich hatte nicht viel Literatur darüber finden können, welche Auswirkungen die Entkopplung auf die Konvergenz haben würde. Wenn mich also jemand in die richtige Richtung weisen oder Ratschläge geben könnte, wäre das großartig.

Justin Dong
quelle
3
Sie entkoppeln sie nicht - Sie wenden eine Aufteilungsmethode an.
David Ketcheson
Danke für die Korrektur. Ich habe diese Frage nicht vergessen und arbeite nur noch daran. Ich werde aktualisieren, wenn es funktioniert ...
Justin Dong

Antworten:

3

Dies ist definitiv das, was als Aufteilungsschema bekannt ist und sowohl in der nichtlinearen Optik als auch in Quantensimulationen sehr beliebt ist. Ein Beispiel ist die Dirac-Gleichung, aber Sie können dies auch auf nichtlinearen PDES tun.

u1k+1u2ku1k+1u2k+1

u1u2u1u2Lu1f=f2(u2)G(x;u2)LG=δ(xs)u1

u1=ΩG((xs);u2)fds

u2k+1u1k+1u2k+1

Es gibt definitiv Flexibilität bei dieser Entkopplung oder Aufteilung, obwohl dies zu einem schönen Schema in Ihren Gleichungen führen könnte. Eine schöne Einführung in die Aufteilung in dem Kontext, an dem Sie interessiert sind, finden Sie in Kapitel 13 'Aufteilung und ihre Verwandten' von Boyd 'Chebyshev- und Fourier-Spektralmethoden'.

Hoffe das hilft

Mike MacNeil
quelle
1

Es ist ein bisschen schwer zu sagen, ohne genau zu wissen, welche PDE Sie lösen. Das heißt, ein paar Orte, um vielleicht zu schauen, sind

  1. Für gekoppelte lineare Systeme sind Uzawa-Methoden und erweiterte Lagrange-Methoden Beispiele für allgemeine Aufteilungsschemata. Es sind Ergebnisse bezüglich der Konvergenz dieser Methoden bekannt, die auf den Spektren der Blöcke und dem Schur-Komplement des gekoppelten Blocksystems basieren (siehe Bramble et al. Zur Lösung von Sattelpunktproblemen).
  2. Es gibt eine Fülle von Materialien zu verschiedenen Aufteilungsmethoden für PDEs. Obwohl diese weitgehend für lineare Probleme formuliert sind, gibt es reichlich Verwendung und Theorie für ihre Erweiterung auf nichtlineare Gleichungen - zumindest für die nichtlinearen PDEs, die Arten von Konvektionsdiffusion und inkompressiblem Fluss regeln, haben Aufteilungs- / Projektionsschemata eine ziemlich reiche Geschichte .

u1u2u1u2

Jesse Chan
quelle