[pod]

Hello everybody,
I am trying to simulate a phase transition (liquid vapor) without involving in chtFoam neither eulerFoam method nor interFoam method. So basically just by accounting for the density change according to temperature. The simulation seems to be really unstable because of the density correction in the pressure calculation. So I was wondering:
- I tried to provide the solver with interpolating tables of density and calculating the compressibility just by dividing by the pressure; I also tried to provide tables of compressibility as well; but none of these options seems to work, what should it be the right way to calculate rho and psi while having a stable code?
- is it even possible to run a simulation in chtMultiRegionFoam with such big density jumps without implementing additional parts like in multiphase solvers? Could it be an option to change the pressure algorithm of this solver to the psi-based one of compressible solvers like rhoPimpleFoam to make it stable and absorb the density correction term?
thank you for any advice,
cheers