Hi,
ich bin kein Fluent Experte aber die Gleichungen sollte ich recht gut verstehen. Es kann natürlich auch sein, dass Fluent irgendwelche anderen Gleichungen löst aber im Allgmeinen ist es so, dass die Referenzdichte nur dann benötigt wird, wenn man die Boussinesq Approximation verwendet (-> inkompressible). Sollte man kompressible Fluide betrachten, hat man die Dichte in den Gleichungen, daraus resultiert die Notwendigkeit der Verwendung des idealen Gasgesetzes (oder einem Realgasgesetz), wodurch man die Dichte berechnet, die dann wieder in der Impulserhaltung auftaucht und als Quelle für den Auftrieb wirkt.
Anmerkung: Normalerweise wird in FOAM dieser Term anders behandelt, also nicht mit rho*g in der Impulserhaltung sonder, als extra Term in der Druckgleichung; man rekonstruiert die Bouyancy, dass bessere Ergebnisse erzielt - geht hier aber zu weit.
Kurzum: OpenFOAM bietet dir folgendes an (du weißt das sicher schon)
- Inkompressible mit Buossinesq Approximation für die Auftriebskraft
- Kompressible ohne Buossinesq Approximation; Auftriebskraft bzw die Dichte Änderung wird mittels einer Gasgleichung berechnet
Da wir beim inkompressiblen Löser alles durch die Dichte teilen, spielt der absolute Druck auch keine Rolle mehr, da die Dichte nicht mit dem idealen Gasgesetz berechnet wird. Es interessiert quasi nur noch die Druckdifferenz (dp). Daher setzt man einfach p auf 0 [Pa / kg / m^3]. Da nur dp interessiert kann man auch negative Druckwerte erhalten, basierend auf deinem pseudo Druckwert den du einstellst. Wenn du aber aus einem Grund auf deine Atmosphäre Bezug nehmen möchtest, dann wäre das für mich wie folgt:
- pRef = 1e5 Pa
- p für den Patch = pRef / Rho
Anstelle von p = 0 (fixedValue) setzt du einfach p = 1e5/rho (fixedValue). An der Berechnung ändert sich nichts, da nur dp zählt und ob du jetzt (0 - 10) oder (800 - 810) hast ist egal. Den Druck in Paraview musst du natürlich mit der Dichte multiplizieren, damit du deine korrekten Werte bekommst. Andere Möglichkeit wäre, die Erstellung eines neue Feldes.
Vielleicht war das ja die Antwort auf deine Frage.
------------------
Viele Grüße,
Tobias Holzmann
OpenFOAM Tutorials | Publikationen | Für Anfänger wiki.openfoam.com
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP