[su_junwei]

Hi Foamers
The following thread discussed the implementation of compressibleInterFoam.
http://www.cfd-online.com/Forums/ope...interfoam.html

I am confused with the implementation of dgdt in the file pEqn.H
In this thread, we can get.
dgdt = { ( psi2 / rho2) - (psi1 / rho1) } DpDt
However, in the implementation, it gives

dgdt = (pos(alpha2)*psi2/rho2 - pos(alpha1)*(psi1/rho1))*(p_rghEqnComp & p_rgh);

Is this right? shoud it be
dgdt = (psi2/rho2 - psi1/rho1))*(p_rghEqnComp & p_rgh);

where are alpha1 and alpha2 from?

would you please give me some tips?

Best Regards, Junwei