Hi, Could anybody give an advise how to access pressure corrections, if it is possible? Thank you.

If you wish to print and plot the pressure correction, just add STORE(PCOR) to q1 (I hope I recall correctly the name is PCOR). If you need it in ground, I guess (but never tried) you can access it in the same way as other variables, assuming it is STOREd.

I hope it helps,

Rami

Dear Rami, thanks a lot, I'll try. I need PCOR in GROUND to modify them for the variable-density flow (I'm using single-fluid solver to simulate liquid-vapour flow with a phase transition). Unfortunately the solution doesn't converge for the steady-state and I look for methods to make iterations more stable. I have found that for shock-waves problems Malin and Sanchez have introduced some compressibility corrections (don't know details). Do you know, if this could be helpful and how to activate this option? Thank you for your help.

martin ,

I am familiar with the paper of Malin et al published in the PHOENICS J. I myself had a lot of difficulties with simulation of inviscid normal shock waves (underexpanded jet), but that was in an old version. I am not sure how far this problem was solved since then. Check POLIS - I vaguely recall there is a flag in group 8 or 9 dealing with something to do with variable density, but I had never tried it.

What exactly are you solving? Are you dealing with normal shocks?

Rami

I'm trying to simulate cavitation flows. For this purpose I have added an auxilary eqn. for the vapour mass fraction with a source term. I know there is SEM and HOL methods in PHOENICS, but not sure if SEM is compatible with the CCM solver and NONORT=F, which i use at present. Thanks a lot