Hallo Foamer's
ich schriebe der Zeit meine Maserarbeit. In diesem Zusammenhang ich bin grade dabei eine Randbedingung zu programmieren, mit welcher ich die thermisch träge Masse der Wand abbilden kann (ohne diese mit dem chtMultiRegion zuvernetzen)[Nutzung von bouyantPimpleFoam]. Da die Geometrie extrem aufwendig ist und es in erster Linie um den Fluidraum geht. Somit soll die Wand nur "fiktiv" in Form einer BC abgebildet werden.
Dabei habe ich ein Problem, da ich noch nicht so recht weiß, wie ich auf eine Variable zugreife (und ggf. speichere als Feld für jedes Patch. Quasi als zweites Temperaturfeld) welche innerhalb der BC berechnet wird. Ich würde mich sehr freuen, wenn ihr mir helfen könntet.
Die Randbedingung bau ich wie folgt auf (siehe Skizze):
gegeben:
- Wandaußentemperatur im Fernfeld (T_U)
- Abgeschätztes Alpha(außen) (alpha_a)
- Materialparameter der Wand (cp, s, lamda, ....)
es ergeben sich die beiden DGL's für die beiden Massenpunkte:
Code:
=>Für die Wandinnentemperatur (Wandtemperatur für CFD-Berechnung):
dTi= [(Q°i-Q°l)/(m*cp*0.5)]*dt
Q°l= A*(lamda/s)*(Ti-Ta)Ti = TiOld + dTi
(Q°i -> Konvektiver Wärmestrom kommt aus der CFD)
=>Für die Wandaußentemperatur:
dTa= [(Q°l-Q°a)/(m*cp*0.5)]*dt
Q°a= alpha_a*A*(Ti-Ta)
Ta = TaOld + dTa
Mein Problem besteht darin, dass ich für das Euler-Verfahren den Wert des vorherigen Zeitschrittes benötige. Für die Wandinnentemperatur über:
Code:
Tiold=Ti.oldTime().boundaryFiled();
Ein Problem stellt für mich die Abfrage von Ta dar. Der Zeit weiß ich nicht wie ich diesen Wert am besten abfrage.
-> Kann ich dies über den selben Weg erfolgen wie bei Ti?
-> ist es möglich innerhalb einer BC ein IOobject zu definieren? (So wäre es mir möglich diesen Wert auch in Paraview anzuschauen) oder muss ich dazu ein IOobject im Solver (createFields.) definieren?
------------------
------------------------------------------------------
Vielen Dank und viele Grüße
Christopher
[Diese Nachricht wurde von CF-D am 30. Aug. 2017 editiert.]
[Diese Nachricht wurde von CF-D am 30. Aug. 2017 editiert.]
[Diese Nachricht wurde von CF-D am 30. Aug. 2017 editiert.]
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP