Autor
|
Thema: simpleFoam - Einzelne Elemente am Inlet mit zu hohen Druckwerten (1906 mal gelesen)
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 05. Sep. 2019 08:54 <-- editieren / zitieren --> Unities abgeben:
Hallo zusammen, ich bin schon einige Zeit mit OpenFOAM beschäftigt. In meiner aktuellen Simulation habe ich ein Problem an welchem ich einfach nicht weiter komme. Die Simulation läuft zügig und schön durch. Sie konvergiert recht schnell. Die Druckwerte weichen leicht von der Realität ab. Grund scheinen einige Zellen am inlet zu sein. Am Rand sind stets keine Zellen mit hohen Werten. (siehe Bild anbei). Ich habe das schon einiges getan. Netz verbessert am inlet, etc. Es verändert sich jedoch nichts. Die Werte bleiben. Vielleicht kennt hier jemand das Problem und hat einen Tipp für mich??? Da paste-copy nicht funktioniert werde ich einen Screenshot von U und p beifügen. Vielen Dank schon mal. ------------------ Viele Grüße an alle
[Diese Nachricht wurde von Tempolini am 05. Sep. 2019 editiert.] [Diese Nachricht wurde von Tempolini am 27. Sep. 2019 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 05. Sep. 2019 09:47 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 05. Sep. 2019 10:03 <-- editieren / zitieren --> Unities abgeben:
|
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 05. Sep. 2019 16:51 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 06. Sep. 2019 08:09 <-- editieren / zitieren --> Unities abgeben:
Hallo Shorty, nicht die komplette Seite. Nur ab Downloads. Das ist normal hier. Ich denke du meinst das Thema "Layer". Die Geometrie ist sehr kompliziert. Nur Ein- und Auslauf sind "Rohrabschnitte". Ich bekomme den Layer nicht bis zum inlet vorgezogen. Von daher hast du sicher recht, es ist das Netz. Ich muss hier wohl noch etwas Arbeit reinstecken. Eventuell weiß ich wo der Fehler liegen könnte. An der zu feinen Vernetzung der Kante. Ich schalte das mal aus. Ich habe gestern gelesen man sollte das FeatureEdge refinement ausschalten, wenn man Layer haben möchte. ------------------ Gruß an alle Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 06. Sep. 2019 08:21 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 16. Sep. 2019 10:32 <-- editieren / zitieren --> Unities abgeben:
Hallo Shorty, ich habe inzwischen so einiges mit snappyHexMesh durch und mir auch dein Tutorial angeschaut (vielen Dank hierfür. Ich habe meinen Case auf ein einfaches Rohr runtergebrochen. Selbst hier habe ich das Phänomen. Wandfunktionen sind hierbei auch berücksichtigt. Hast du noch eine Idee? Oder hat jemand anders noch eine Idee? (Lege ich einen Slice knapp hinter den inlet und lese über den calculator in paraview die Druckwerte aus, stimmen sie mit Messergebnissen überein.). ------------------ Gruß an alle Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 16. Sep. 2019 10:53 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 16. Sep. 2019 11:13 <-- editieren / zitieren --> Unities abgeben:
|
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 16. Sep. 2019 18:22 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 17. Sep. 2019 09:51 <-- editieren / zitieren --> Unities abgeben:
Das war keine Absicht. Das Problem ist noch nicht gelöst. Ich habe irgendetwas versehentlich angeklickt. Ich teste gerade die Einstellungen. k läuft im residual-Plot etwas unruhiger. mit "Gauss linear limited corrected 0.5" für laplacianSchemes und snGradSchemes inkl. die Einstellungen für k und Epsilon ist das Ergebnis ein wenig besser. Ich schau mal was mit 0.33 passiert. Gebe ich 1 ein tut sich nichts. ------------------ Gruß an alle Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 17. Sep. 2019 10:17 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
Hi, hab dein "Gelöst" nochmals herausgenommen (kann man mit dem Grünen Pfeil im 'Grafik für den Beitrag' anwählen). Ich habe: Code:
divSchemes { default none; div(phi,U) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div(phi,k) bounded Gauss upwind; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited corrected 0.33; } interpolationSchemes { default linear; } snGradSchemes { default limited corrected 0.33; } wallDist { method meshWave; }
Hab den PCG und PBiCG verwendet und das hier: Code:
SIMPLE { nNonOrthogonalCorrectors 1; residualControl { p 1e-3; U 1e-4; // possibly check turbulence fields "(k|epsilon|omega)" 1e-3; } } relaxationFactors { fields { p 0.8; rho 0.01; } equations { U 0.4; "(k|omega)" 0.3; epsilon 0.3; } }
Zudem eben noch k und epsilon angepasst. Mir ist auch aufgefallen, dass du eine sehr kleine Geometrie hast. Das SST ist hier besser geeignet. k-Epsilon ist schlecht für Wandnahe Bereiche, und das hast du eben genau hier am Inlet. Deine Layer könnte man noch enger ziehen um dem ggf. entgegenzuwirken.
------------------ Viele Grüße, Tobi OpenFOAM® Tutorials | Training | Publikationen | Für Anfänger wiki.openfoam.com Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 17. Sep. 2019 16:18 <-- editieren / zitieren --> Unities abgeben:
Hi, ich habe am Netz der Originalgeometrie noch etwas "rumgefummelt". Aktuell sieht es gut aus. Die Rechnung läuft aber noch. Und das sehr, sehr, sehr viel langsamer als zuvor mit feinerem Netz. Bin mal gespannt wie es am Ende aussieht...
------------------ Gruß an alle Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 18. Sep. 2019 09:49 <-- editieren / zitieren --> Unities abgeben:
Guten Morgen, die Ergebnisse sehen schon super aus. Die Rechnung ist auch schön konvergiert. Vielen Dank für deine Tipps. Ist es normal mit den Solvereinstellungen, dass die Rechnung sehr viel länger dauert obwohl das Netz relativ grob ist? ------------------ Gruß an alle
[Diese Nachricht wurde von Tempolini am 18. Sep. 2019 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 27. Sep. 2019 10:39 <-- editieren / zitieren --> Unities abgeben:
Problem gelöst. Vielen Dank für die Hilfe ;-) ------------------ Gruß an alle [Diese Nachricht wurde von Tempolini am 27. Sep. 2019 editiert.] [Diese Nachricht wurde von Tempolini am 27. Sep. 2019 editiert.] Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |
Shor-ty Moderator
Beiträge: 2463 Registriert: 27.08.2010 OpenFOAM-dev (Foundation) OpenFOAM-xxxx (ESI)
|
erstellt am: 27. Sep. 2019 16:21 <-- editieren / zitieren --> Unities abgeben: Nur für Tempolini
|
Tempolini Mitglied
Beiträge: 16 Registriert: 05.09.2019 Salome Meca 9.3.0 OpenFOAM 7
|
erstellt am: 30. Sep. 2019 09:03 <-- editieren / zitieren --> Unities abgeben:
Hallo Shorty, ich bin mir ziemlich sicher, dass es am Solver lag. Die von dir empfohlenen Solver sind laut meinem Buch, in das nochmal geschaut habe, deutlich besser für mein Problem gedacht. (Danke für den Impuls deinerseits nochmal.) Auch das kOmegaSST hat noch dazu beigetragen letztenendes etwas schneller durchzulaufen. Denn mit kEpsilon war es nicht so stabil und lief damit auch etwas schleppender. Das Netz musste auch gar nicht mehr so superfein sein (habe ich getestet. Ergebnisse unterschieden sich kaum).
------------------ Gruß an alle Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP |