Turbulence viscosity UDF
 

Hi,
I am using k epsilon model for cavitation model. In order to modify the turbulence viscosity with the compressibility effec, I use the Reboud function as: cmut= f(rho) * cmu* k*k/epsilon
In that f(rho)=rho_liquid+rho_vapor*(1-void)^n
I tryed to make the UDF as in Fluent UDF manual: DEFINE_TURBULENT_VISCOSITY(name, c,t)
But I dont know how to call the density of phases: liquid and vapor in the code in order to calculate the density function.
Could you help me how to do this? Thank you very much!