---

Hallo,

ich möchte meine Rechnung etwas stabiler machen und für:

div((muEff*dev2(T(grad(U))))) Gauss linear;

ein anderes Schema nutzen (1.Ordnung dann). Hat jemand einen Vorschlag?

Grüße und Danke

Moritz

------------------
Moritz Franke

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Gibt es einen Grund wieso du hier ein Schema 1. Ordnung möchtest ?
Im Konvektionsterm macht die erst Ordnung Sinn aber das auch nur auf Kosten der Genauigkeit. Du erwähnst Stabilität, welchen Solver verwendest du?

------------------
Viele Grüße,
Tobias Holzmann

OpenFOAM Tutorials gibt es auf www.Holzmann-cfd.de
Diverse Veröffentlichungen und ein Buch gibt es auf www.Holzmann-cfd.de
Für Anfänger empfiehlt es sich die neue OpenFOAM Wiki Seite zu betrachten: wiki.openfoam.com

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Solver: rhoPorousSimpleFoam

Mir kracht ab ca. 6000 Iterationen für ein 20 Mio Netz alles zusammen: Druck steigt plötzlich stark an, Temperatur auch. Und zwar immer an der Stelle, wo mein poröses Medium liegt.

Das lustige: füge ich im Blockmesh 10 Zellen in jede Raumrichtung mehr ein, dann läuft die Rechnung stabil durch. Irgendwas haut da dann wahrscheinlich eher mit dem Netz und dem porösen Medium an dieser Stelle nicht hin

------------------
Moritz Franke

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP