[linch]

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