---

Screenshot.jpg

Hallo Foamers,

Ich sitze schon seit mehreren Tagen an dem Problem, dass meine Rechnung immer wegen omega divergiert.Rechne mit rhoPorousMRFSimpleFoam. Habe mit den Relaxationsfaktoren rumgespielt, habe den solver gewechselt(PCG zu GAMG für p) und habe auch schon ein wenig mit den schemes rumgespielt. Bei jeder Variante schießen irgendwann die Residuen von omega in die Höhe...Bei meiner letzten Variante hatte ich zudem noch starke Schwingungen von p. Am Netz kann es nicht liegen, da ich eine Vergleichsrechnung mit CFX gemacht habe, welche sehr gute Werte geliefert hat. Hat jemand ne Idee was ich falsch mache oder wie ich es ändern kann?
Wäre für jede Hilfe sehr dankbar!

Anbei Residuenverlauf und letzten steps im log-file.

Code:

---

Time = 58

smoothSolver:  Solving for Ux, Initial residual = 0.120016, Final residual = 0.0107659, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.123988, Final residual = 0.0112834, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.135026, Final residual = 0.0121009, No Iterations 2
smoothSolver:  Solving for h, Initial residual = 0.0933729, Final residual = 0.00750188, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.491722, Final residual = 0.0164022, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0880264, Final residual = 0.00438395, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0370632, Final residual = 0.00169125, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.012555, Final residual = 0.000583616, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0041581, Final residual = 0.000194771, No Iterations 13
GAMG:  Solving for p, Initial residual = 0.00151314, Final residual = 7.48414e-005, No Iterations 14
GAMG:  Solving for p, Initial residual = 0.000596276, Final residual = 2.8951e-005, No Iterations 19
GAMG:  Solving for p, Initial residual = 0.000272889, Final residual = 1.35552e-005, No Iterations 38
GAMG:  Solving for p, Initial residual = 0.000140946, Final residual = 1.21018e-005, No Iterations 1000
time step continuity errors : sum local = 0.885654, global = -0.00538908, cumulative = 0.192475
rho max/min : 4.93471 0.0951114
smoothSolver:  Solving for omega, Initial residual = 0.0466452, Final residual = 0.00220318, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0928436, Final residual = 0.00278674, No Iterations 2
ExecutionTime = 32971.1 s  ClockTime = 32971 s

Time = 59

smoothSolver:  Solving for Ux, Initial residual = 0.119863, Final residual = 0.0108241, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.123468, Final residual = 0.0113648, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.134143, Final residual = 0.0122472, No Iterations 2
smoothSolver:  Solving for h, Initial residual = 0.098266, Final residual = 0.0077809, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.497361, Final residual = 0.0170733, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0898226, Final residual = 0.00436005, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0328101, Final residual = 0.00150678, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.0134029, Final residual = 0.000647348, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.00378356, Final residual = 0.00018327, No Iterations 15
GAMG:  Solving for p, Initial residual = 0.00171656, Final residual = 8.16549e-005, No Iterations 12
GAMG:  Solving for p, Initial residual = 0.000587276, Final residual = 2.83591e-005, No Iterations 21
GAMG:  Solving for p, Initial residual = 0.000301085, Final residual = 1.48817e-005, No Iterations 26
GAMG:  Solving for p, Initial residual = 0.000146648, Final residual = 1.24209e-005, No Iterations 1000
time step continuity errors : sum local = 0.948546, global = -0.0035274, cumulative = 0.188948
rho max/min : 4.94124 0.0856002
smoothSolver:  Solving for omega, Initial residual = 1.#QNAN, Final residual = 1.#QNAN, No Iterations 1000

---

Solltet ihr zusätzlich nach Daten zur Fehleranalyse brauchen, kann ich die gerne noch reinstellen.

Irgendwie kommt es auch immer zur Divergenz, wenn rho min/max fast die von mir eingestellten Werte erreichen. Was beschreibt rho min/max eigentlich? bin davon ausgegangen, dass es die Grenzen sind, welche mein rho nicht überschreitet bzw. überschreiten darf. Liege ich da richtig?

Gruß, Eric

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

---

Hallo,

hört sich immer so einfach an, aber prüfe doch mal, ob du dir sicher bist, dass die Randbedingungen so stimmen...

"GAMG:  Solving for p, Initial residual = 0.000140946, Final residual = 1.21018e-005, No Iterations 1000"

Das der GAMG 1000 Iterationen und mehr bräuchte macht einen nachdenklich...

Gruß Thomas

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