---

Hallo und Grüße an alle FOAMer hier,

meine Frage bezieht sich auf folgendes

Löser: chtMultiRegionFoam

Tutorial: heatedDuct

Meine einzige Änderung an heatedDuct - in Allrun:

ideasUnvToFoam Mesh_1.unv
#runApplication blockMesh
#runApplication snappyHexMesh -overwrite

Mesh_1.unv  habe ich in Salome erstellt mit der Geometrie genau so wie in heatedDuct, nur die Längenmaße sind anders.
Alle Regionen, Regionen-Namen sind identisch mit Tutorial heatedDuct.

OpenFAOM 9 rechnet alle Zeitschritte ohne Fehlermeldung.

In paraView sehe ich dann fehlende Wärmeleitung zwischen "heater" und "metal".

Alle Region zu Region Kontktflächen werden in paraView gelistet, auch heater/heater_to_metal ist vorhanden.

Wo soll ich in heatedDuct korrigieren damit auch Wärmeleitung berechnet wird?

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

---

Grüß Dich und willkommen im Forum,

snappyHexMesh erstellt auch noch cellZones

Code:

---

fluidToMetal                                                           
        {                                                                     
            level (1 1);                                                       
            faceZone fluidToMetal;                                             
            cellZone metal;                                                   
            mode        insidePoint;                                           
            insidePoint (0.025 0.0025 0.075);                                 
        }                                                                     
                                                                               
        metalToHeater                                                         
        {                                                                     
            level (1 1);                                                       
            faceZone metalToHeater;                                           
            cellZone heater;                                                   
            mode        insidePoint;                                           
            insidePoint (0.025 -0.005 0.075);                                 
        }           

---

Die man anschließend Splittet:

Code:

---

runApplication splitMeshRegions -cellZones -overwrite 

---

Daher weiß ich nicht was Du wirklich machst und ob Dein Case auch wirklich so aufgebaut ist wie der originale Tutorial Case. Ohne Case keine Analyse, ohne Analyse kein Support. Ich könnte genauso gut Fragen wieso mein Kuchen nicht schmeckt.

------------------
Glück Auf,
Tobi

OpenFOAM® Community - Knowledge Base

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

---

heatedDuct-TEST.zip

Danke Tobi für schnelle Antwort,

mein Case - im Anhang.

Das ist original heatedDuct aus Tutorial mit Änderung nur in Allrun. Und mein Mesh-1.unv

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

---

Danke für den Case. Wärmeübertragung hast Du:

Code:

---

Solving for solid region metal
GAMG:  Solving for e, Initial residual = 0.472617, Final residual = 5.09155e-06, No Iterations 1000
Min/max T:300 300.003

---

Dein metal wird tatsächlich etwas wärmer. Aber 1000 Iterationen im Solid sagt uns, dass da was total falsch ist. Der GAMG ist dafür nicht wirklich geeignet:

Code:

---

system/metal/fvSolution
-----------------------

solvers                                                                       
{                                                                             
    e                                                                         
    {                                                                         
        solver          PBiCGStab;                                           
        preconditioner  FDIC;                                                 
        tolerance        1e-6;                                                 
        relTol          0.1;                                                 
    }           

---

Resultiert in:

Code:

---

Solving for solid region metal
FDICPBiCGStab:  Solving for e, Initial residual = 0.0194605, Final residual = 1.84419e-12, No Iterations 1
Min/max T:300 306.492

---

Allerdings ist das jetzt auch nicht das was das Problem ist. Es ist nämlich viel einfacher, dass auch aus dem GAMG (mit Erfahrung) ersichtlich ist. Ein Blick in Paraview und der Dimensionen ergibt, dass Du Deine Geometrie warhscheinlich um den Faktor 1000 zu groß hast.

Die Lösung:

Code:

---

transformPoints "scale=(0.001 0.001 0.001)" -allRegions

---

Dadurch wird auch der Zeitschritt wesentlich kleiner und man sieht in der Tat etwas.

Ich hab übrigens noch Deine ganzen Tets mit polyDualMesh 90 in ein anständigeres Netz konvertiert und die Numerik noch etwas angepasst, dass aber nicht notwendig ist.

------------------
Glück Auf,
Tobi

OpenFOAM® Community - Knowledge Base

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

---

Danke für ausführliche Analyse Tobi,

ich habe jetzt in Allrun die Skalierung nach Netzkonvertierung gesetzt

ideasUnvToFoam Mesh_1.unv
transformPoints "scale=(0.001 0.001 0.001)" -allRegions

und diese Meldung bekommen
------------------------------------
--> FOAM FATAL ERROR:
Cannot find file "points" in directory "fluid/polyMesh" in times "0" down to constant
------------------------------------

obwohl file "points" in "fluid/polyMesh" vorhanden ist.
Warum kann OF es nicht finden?

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

---

Das liegt daran, dass nach dem ideasUnv nur ein numerische Netz in constant/polyMesh liegen sollte. Anschließend sagst du aber mit -allRegions, dass er die Regionen skalieren soll.

Entweder du splittest Dein Netz vorher auf und machst -allRegions oder Du skaliert dasNetz um constant/polyMesh und splittest dann.

------------------
Glück Auf,
Tobi

OpenFOAM® Community - Knowledge Base

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

---

Nach dieser Korrektur funktioniert nun alles perfekt, nochmals danke!  

Als nächstes werde ich den case von laminar auf turbulent umrüsten.

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