| | | 3D-Druck: 7 Gründe für den Einsatz in der Medizin, ein Fachartikel
|
Autor
|
Thema: divergenz durch omega (1433 mal gelesen)
|
mump0 Mitglied Student
Beiträge: 20 Registriert: 24.10.2011
|
erstellt am: 16. Dez. 2011 09:36 <-- editieren / zitieren --> Unities abgeben:
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 = 58smoothSolver: 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 |
TTB Mitglied CFD Engineer
Beiträge: 353 Registriert: 02.10.2008 BIM HVACTool für Windows OpenFOAM-2.2.x
|
erstellt am: 20. Dez. 2011 12:33 <-- editieren / zitieren --> Unities abgeben: Nur für mump0
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 |
| Anzeige.:
Anzeige: (Infos zum Werbeplatz >>)
|