Hi Nima,

Update: I've added 2 rhoEqns, for rho1 and rho2 at the end of my time step. Now I force the continuity of each phase. The drawback is the the EOS is not exactly fulfilled in some interface cells, that caused problems earlier, but I can live with it, since the continuity is more important.

Thanks a lot for the discussion, it helped me a lot. If if you're interested, you can still give me your solver some day, and I could test it with my case. I have a high thermal expansivity of both phases and it may be interesting to know, if your solver still performs well. It must not be now, some day you wish.

Best wishes and a nice weeked!
Illya

Any references on your work?
 

Hi,
first need to thank both of you for the valuable discussion and I'm really sorry as I'm trying to refer your work 6 years ago. I'm trying to simulate high pressure gas atomization problem and molten metal is atomizing due to the high pressure gas. The gas inlet is a pressure boundary condition and pressure is in the range of 15-30 atm. However, I tried compressibleInterFoam in OF40 and it crashes with negative temperatures. I know it's due to the sudden expansion of the gas. I tried different solver conditions but still getting the same error. I assume it's a result of taking density as a function of pressure only.

I trying to workout the math to obtain the effect of rho(P,T), following the method described in "A pressure based compressible, two phase flow finite volume method for underwater explotions" by Miller et al. I'm still trying to understand the PEqn.H. However, seems like you to have made some progress on this.

Can you please give me some references of your work so that I can follow.

Thank you.
Kalpana