Simulation of a closed system with interThermalPhaseChangeFoam
 

Hello to everyone, I’m trying to simulate a closed system whit a fluid that change phase inside of it using the interThermalPhaseChangeFoam (Solver developed by Nabil et al.), running in openFoam 2.4.
More in particular I want to situate a cylinder with some fluid on the bottom, this cylinder is heated at the bottom and cooled on the top, theoretically, the fluid must evaporate, flow through the domain, and then it will condense in the cold region.
My problem is the fact that having only Neuman BC for the pressure I need to set pRefValue in the pRefCell but when the fluid change phase the pressure in the whole domain change everywhere except that in the pRefCell, this generate unfisically flow in the direction of that cell.
Someone can help me?
Thank you so much