|
OpenFOAM : Flux Limiters
Shor-ty am 14.10.2020 um 09:10 Uhr (1)
Hey Friendly,ich habe mit dem rhoCentralFoam noch nicht gearbeitet. Allerdings sollte hier eben eine Gleichung für die Dichte gelöst werden und über die EOS dann der Druck. Bei den druck-basierten Lösern ist es ja genau andersrum.Ob hier eine Limiter Verwendung findet kann ich Dir nicht sagen. Allerdings ist ein Limiter nicht automatisch eine Reduktion auf ein Upwind-Verfahren. Die Limiter-Funktionen haben ja eine spezielle Form.Übrigens Reconstructing ist das Verfahren, um von den Face Values auf den Cell ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Layers hinzufügen in Multi-Region Case
Shor-ty am 01.07.2021 um 20:32 Uhr (1)
Hi,also auf meiner Seite gibt es verschiedene Tutorials. Weiß nicht welches Du jetzt meinst. Einfach mal etwas stöbern und die Beschreibungen durchlesen (ggf. Suchfunktion nutzen). Prinzipiell ist Dein Vorgehen korrekt:Netz erstellen ohne LayerSplitten der RegionenEinzeln ein constant/DeineRegion/polyMesh in constant/polyMesh kopierensnappy mit layer Erstellung ausführen (-overwrite nicht vergessen)Das Netz dann wieder von constant/polyMesh in constant/DeineRegion/polyMesh kopierenFür alle Netze/Regionen w ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Wärmetauscher mit Druckeinlass
Shor-ty am 29.11.2021 um 15:35 Uhr (1)
Hey,das p Feld ist belanglos und die Randbedingungen darin sollten alle auf calculated stehen. Wir rechnen mit p_rgh. Damit ist wohl klar wieso nichts passiert:p_rgh braucht den DruckgradientFixiere InletFixiere OutletU mit pressureInletVelocityp_rgh wird am Inlet prghTotalPressure sein bzw. prghTotalHydrostaticPressure in Kombination mit der hydrostaticPressure FO für die Initialisierung von ph_rghp_rgh wird am Outlet identisch zum Inlet sein (fixedValue ist ungeschickt weil das kein Fixed value ist aufgr ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : druck graient zwischen refinement zones
Shor-ty am 17.10.2021 um 09:04 Uhr (1)
Ich glaub wir reden aneinander vorbei:Code:simulationType LES;LES{ LESModel SpalartAllmarasIDDES; turbulence on; printCoeffs on; delta IDDESDelta; IDDESDeltaCoeffs { hmax maxDeltaxyzCubeRoot; maxDeltaxyzCubeRootCoeffs { } }} Die Filterfunktion braucht das delta. Normal sollte für LES das Netz einheitlich sein und Co ~ 1. Praktiabel ist das natürlich nicht. Man verwendet in der Industrie hier gerne DES (oder DDES) (wird nur über Tubulenzmodell gesteuert), daher weiß ich ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : AMI bewegt sich nicht
Shor-ty am 15.03.2021 um 18:17 Uhr (1)
Wissen tut man das schon nur kannst Du nicht erwarten, dass sich zu jederzeit jemand um Deine Themen kümmert pRef muss man angeben, wenn man keine Druckrandbedingung definiert und sich somit unendlich viele Lösungen einstellen könnten. Es ist nur der Druckgradient entscheidend und ob wir jetzt 1000 - 300 als Gradient berechnen oder 25735 - 25035 ist vollkommen egal.Wo bekommst Du denn deine Extremwerte? Ist dein Case noch aktuell oben?Du erwähnst immer ein Tutorial - welches denn?Und ja movingWallVelocit ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : AMI bewegt sich nicht
Shor-ty am 15.04.2021 um 14:03 Uhr (1)
Eine Zeitschrittweite von 1e-9 habe ich nur bei Strömungen bei Ma 1. Ich gehe davon aus, dass Du das nicht hast. Daher ist schonmal irgendwas falsch oder macht Probleme. Wohl möglich eine Zelle die Dir das versaut. Schon mal deine Geschwindigkeiten angeschaut? Bei 1e-9 wäre ich sehr kritisch. Zum anderen Thema. Ja, mit PIMPLE kann man die Courant-Zahl deutlich 1 halten. Bspw. auch 200. Aber dann verliert man Zeitinformationen. Sofern das nicht von großer Bedeutung ist, kann man das natürlich durchaus mac ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Decomposepar
Shor-ty am 17.12.2021 um 13:33 Uhr (1)
Frage 1: Ja Du liegst falsch. Decomposing auf 2 Kernen gibt ein anderes Ergebnis wie mit 4 oder 2000 (überspitzt ausgedrückt).Vor allem wenn nur kleine Kräfte die Strömung bewegen, kannst Du durch Processor-Boundaries das etwas verfälschen.Du solltest herausfinden was am Case mit 4 Kernen anders ist (schau Dir die Zeitschritte an).Zudem:- Vernetzung niemals mit scotch -- hierarchical oder simple- Simulation dann mit scotchFrage 2:- Verwende fvOptions: velocityDamping + limitTemperature (sofern Du eine Ener ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Bestimmung des Strömungswiderstandkoeffizienten für einen porösen Zylinder in OF
Shor-ty am 07.10.2020 um 08:55 Uhr (1)
Hallo zusammen,schon mal die Funktion "forces" angeschaut? Ich verweise mal auf das hier:Code: dict.readIfPresent("porosity", porosity_); if (porosity_) { Info(log_) " Including porosity effects" endl; } else { Info(log_) " Not including porosity effects" endl; }Gegenenfalls etwas, dass Ihr euch anschauen könnt. Ist übrigens auch bei »forceCoeffs« integriert.Ach übrigens, man kann in FOAM schon interne Faces auswählen und intern ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Unterschied grad(U) vs. grad(phi), paraview
Shor-ty am 24.02.2023 um 12:01 Uhr (1)
Heieieieie...meine Antwort ist hier wohl nie angekommen. Dann nochmals aber kürzer.Der Grund ist einfach:phi ist ein surfaceScalarField - also eine skalare Größe an jedem Face. Ein grad(phi) macht daraus eine Vektorgröße (super für Glyps) (Tensor 1ster Ordnung)U ist ein volVectorField - also ein Vektor im Zellmittelpunkt. Ein grad(U) macht daraus dann einen Tensor (wie soll ein Glyp Filter eine Matrix analysieren?) (Tensor 2ter Ordnung)Was Du suchst ist grad(mag(U)). Damit sollte alles dann geklärt sein ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Boundary-Condition-Typ für Inlet gesucht: bewegtes Netz plus Polynomprofil
Shor-ty am 16.01.2023 um 22:32 Uhr (1)
Hallo,das gibt es so nicht. Ich denke da musst Du über codedFixedValue oder exprFixedValue arbeiten (letzteres bin ich mir nicht sicher ob es geht). Mit dem ersten gehts. Das gleiche (beide oben genannten) kannst du mit der uniformFixedValue BC erreichen, wobei bei value dann entweder coded oder expression steht.Die Rotation des Profils musst Du herleiten und dann eben draufgeben. Da gibts wohl viele Möglichkeiten. Man könnte sicher auch mittels der coded BC auf die Rotations-Klassen zurückgreifen. Hab ich ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Temperatursteigerung durch Zufügung einer Leistung
Shor-ty am 13.10.2020 um 10:38 Uhr (1)
Hallo und willkommen im Forum Becoeu,ich entnehme Deiner Schilderung, dass Du beide Regionen einzeln vernetzt hast. Deshalb wirst Du die Kugeloberfläche und die entsprechende Fläche in Deiner Fluiddomain wohl als wall definiert haben. Hier muss eine mappedWall Typdefinition stehen (ich bezieh mich auf Deine boundary Datei die Du unter constant/region/polyMesh findest).Ansonsten kann ich nur sagen, dass ohne Informationen auch keine Hilfestellung möglich ist. Wäre genauso, als würde ich sagen, mein Kuchen s ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Error: Fehlendes interFaceCompositionModel, bitte um Tipps zur Vorgehensweise
Shor-ty am 17.12.2020 um 09:14 Uhr (1)
Hi Teresa,Deine Fehlermeldung sagt dir, dass Du diese Kombination für die Thermodynamik nicht vorliegen hast. Die Boussinesq Approximation gibt es in diesem Zusammenhang nicht (nur incompressiblePerfectGas und rhoConst). Wenn Du die Code-Zeilen dazufügst kann es sein, dass Du es kompilieren kannst, gehe aber davon asu, dass Du damit nichts anfangen kannst, sofern die Modellkombination nicht vorliegt. Das findest Du aber schnell selber raus. Deine Aufgabe ist ziemlich komplex, dass möchte ich Dir nur mal sa ...
|
| In das Form OpenFOAM wechseln |
|
OpenFOAM : Temperatursteigerung durch Zufügung einer Leistung
Shor-ty am 14.10.2020 um 13:31 Uhr (1)
Hey Beckoeu,mit dem Case kann man arbeiten. Es hätte auch gereicht, wenn Du sagst, du machst das Netz mit snappyHexMesh. Wie man damit arbeitet ist mir klar (siehe meine Fußzeile). Aber ich hab deine Randbedingung nicht richtig gelesen:Code: luft_to_gross { type externalWallHeatFluxTemperature; value uniform 300; mode flux; q uniform 0; ...
|
| In das Form OpenFOAM wechseln |
