Finite-Differenzen-Koordinatentransformation für sphärische Polarkoordinaten

8

Ich habe einen Teil eines Problems, das durch die Impulserhaltungsgleichung beschrieben wird:

ρt+1sinθθ(ρusinθ)=0

Wobei und (konstante Geschwindigkeit).u=f(θ)ρ=f(θ,t)

Naiv könnte man eine der hier aufgeführten Lösungen anwenden . Das vorliegende Problem lässt sich am besten mit sphärischen Polarkoordinaten (dünne Kugelschale) und diesen Lösungen auf Kartesisch beschreiben. Muss ich vor der Diskretisierung dieser Gleichung eine Art Koordinatentransformation durchführen oder kann ich sie direkt diskretisieren?

Zweitens gibt es einen Grund, warum man zuerst die Ableitung inθ und dann versuchen sollte, zu diskretisieren?

Als Anmerkung: Ich habe mehrere der oben genannten Schritte ausgeführt und Lösungen erhalten, die nicht konsistent zu sein scheinen (physisch scheint ein Paar sinnvoll zu sein). Ich bin daran interessiert, ob eine ordnungsgemäße Koordinatentransformation durchgeführt werden sollte oder ob eine der zuvor genannten Methoden ausreicht.

BEARBEITEN:

Ich definiere Fluss als:

Φi+1/2=ui+1/2+|ui+1/2|2ρisinθi+ui+1/2|ui+1/2|2ρi+1sinθi+1

Φi1/2=ui1/2+|ui1/2|2ρi1sinθi1+ui1/2|ui1/2|2ρisinθi

Ich vermute, der 'richtige' Weg, den Fluss zu definieren, besteht darin, an den 'Zellgrenzen' und nicht im auszuwerten . Dies würde eher der Definition des Flusses entsprechen.sinθ±12

Eine letzte Frage - was soll ich an den Grenzen tun ( ist besonders ein Problem und ich vermeide diesen Punkt einfach ganz).θ=0

Marm0t
quelle
2
Ich sehe in deiner Gleichung. Bitte korrigieren Sie sin \ theta ( ) in \ sin \ theta ( ): Die richtige mathematische Formatierung ist wichtig! hsinθsinθ
Stefano M

Antworten:

3

Es gibt keinen Grund, die Koordinaten zu transformieren, bevor Sie die Gleichung diskretisieren. Da Sie jedoch sphärische Koordinaten verwenden, erhalten Sie eher ein nichtlineares System als das lineare System, das in kartesischen Koordinaten erscheinen würde. Wenn Sie mit einem linearen System arbeiten möchten, müssen Sie variable Transformationen durchführen (ähnlich wie bei der Erstellung von , die in der ersten Antwort Ihrer verknüpften Frage angezeigt wird).Φ

In Bezug auf die Erweiterung der Ableitung in würde ich empfehlen, dies nicht zu tun, es sei denn, Sie müssen für die obigen Transformationen. Die Gleichung in ihrem aktuellen Zustand liegt in der konservativen Form vor, wobei alle mit gekoppelt sind . Wenn Sie die Ableitung erweitern, verlieren Sie diese Kopplung. Die beiden Lösungen konvergieren mit immer feineren Maschen zueinander, es besteht jedoch ein Unterschied ungleich Null zwischen ihnen. Jemand anderes mit mehr Erfahrung kann möglicherweise auf die Vor- und Nachteile der beiden Formen eingehen, da Finite-Differenzen-Methoden, obwohl sie als konservative Form der Gleichung bezeichnet werden, von Natur aus auf einer bestimmten Ebene nicht konservativ sind .θuρ

Godric Seer
quelle
Vielen Dank für die Antwort. Ich habe mögliche Lösungen für dieses Problem untersucht, aber noch keine vernünftige Lösung gefunden. Ich habe die BCs ausgefallen - aber keine Benchmarks, an denen ich mein Schema testen kann.
März
Was war Ihre Meinung hinter der Formulierung von ? Es scheint, dass Sie Ihren Index für basierend auf dem Vorzeichen von auswählen . Ist das für die Stabilität? Φρsin(θ)u
Godric Seer