CFL-Zustand in diskontinuierlichen Galerkin-Schemata

9

Ich habe ein ADER-diskontinuierliches Galerkin-Schema zur Auflösung linearer Systeme von Erhaltungsgesetzen vom Typ implementiert und festgestellt, dass die CFL-Bedingung sehr restriktiv ist. Im Literaturverzeichnis, eine obere für den Zeitschritt gebunden Δ t htU.+EINxU.+B.yU.=0 kann gefunden werden, wobeihdie Zellengröße ist,ddie Anzahl der Dimensionen ist undNder maximale Grad der Polynome ist.Δthd(2N.+1)λmeinxhdN.

Gibt es eine Möglichkeit, dieses Problem zu umgehen? Ich hatte mit WENO-ADER Finite-Volumenschemata gearbeitet und die CFL-Beschränkungen waren viel lockerer. Beispielsweise muss für ein Schema 5. Ordnung eine CFL von weniger als 0,04 auferlegt werden, wenn DG verwendet wird, während CFL = 0,4 weiterhin in einem WENO-ADER FV-Schema verwendet werden kann.

Warum sollten DG-Schemata anstelle von ADER-FV beispielsweise in der rechnergestützten Aeroakustik (linearisierte Euler-Gleichungen) oder ähnlichen Anwendungen (Gasdynamik, Flachwasser, Magnetohydrodynamik) verwendet werden? Sind die Gesamtrechnungskosten des Schemas trotz des viel geringeren Zeitschritts ähnlich wie die des ADER-FV?

Gedanken und Vorschläge dazu sind willkommen.

Adr
quelle

Antworten:

6

L.2hN.Ö(h/.N.2)

N.

N.N.

DG-Methoden sind teurer, können jedoch problemlos mit unstrukturierten Netzen umgehen und effizient implementiert werden. Es gibt Versionen höherer Ordnung von WENO (oder ähnlichen Rekonstruktionen) für unstrukturierte Gitter, obwohl diese zusätzliche mathematische oder implementierende Komplikationen mit sich bringen können.

Jesse Chan
quelle
Vielen Dank für Ihre ausführliche Antwort Jesse, sie hat mir eine breitere Sicht auf dieses Thema gegeben. Bei meinen numerischen Versuchen mit dem DG-ADER habe ich festgestellt, dass bei Verwendung strukturierter Vierecksnetze (mit beliebiger Viereckform, z. B. Quadraten, Trapezoiden oder Parallelogrammen ...) die numerische Lösung nicht oszillierend und konvergent zur exakten Lösung ist Wenn Sie sich jedoch zu unstrukturierten Netzen bewegen, treten selbst bei quasi strukturierten Netzen Schwingungen auf, die durch zufälliges Verschieben der Knoten eines strukturierten Netzes um einen kleinen Abstand erzeugt werden. Ist das ein erwartetes Verhalten?
Adr
1
h
h