[helena97]

Hello,
I am simulating multiphase problem of cavitation with k-omega SST turbulence model and I am trying to make UDF to reduce the turbulent visocosity, especially in areas with low void ratios. I have already done it for k-epsilon RNG model (and it worked well) according to:
COUTIER-DELGOSHA, O., R. FORTES-PATELLA a J. L. REBOUD. Evaluation of the Turbulence Model Influence on the Numerical Simulations of Unsteady Cavitation. Journal of Fluids Engineering. 2003, 125(1), 38-45. ISSN 0098-2202. Available at: doi:10.1115/1.1524584

However, I am struggling with the UDF function for k-omega SST. When I try to interpret the function in Fluent, it says: line 31: unknown data op (ld.pointer) in block data.
Nonetheless, there is no line 31. I have already checked for semicolons and parantheses...
The equations for k-omega SST model, I have used, are described e.g. in Appendix of this article: https://pdfs.semanticscholar.org/767...640.1634903874

And here is my function:

#include "udf.h"
#include "mem.h"
#define mojeMax(a,b) (((a)>(b))?(a)b))

DEFINE_TURBULENT_VISCOSITY(cavitation_mu_t, c, t)
{
int n = 10;
Thread** pt = THREAD_SUB_THREADS(t);

real a1 = 0.31;
real k = C_K(c,t);
real y = C_WALL_DIST(c,t);
real omega = C_O(c,t);
float OMEGA[3] = {C_DWDY(c,t) - C_DVDZ(c,t), C_DUDZ(c,t) - C_DWDX(c,t), C_DVDX(c,t) - C_DUDY(c,t)};
real OMEGA_V = pow(SQR(OMEGA[0]) + SQR(OMEGA[1]) + SQR(OMEGA[2]), 0.5);

real arg2 = mojeMax((2-pow(k,0.5)/(0.09*omega*y)), (500*0.0000010034/(SQR(y)*omega)));

real F2 = tanh(SQR(arg2));

real nut = (a1 * k)/mojeMax(a1*omega, OMEGA_V*F2);


real rhoM = C_R(c, t);
real rhoL = C_R(c, pt[P_PHASE]);
real rhoV = C_R(c, pt[S_PHASE]);
real f_of_rho = rhoV + pow((rhoM - rhoV)/(rhoL - rhoV), n)*(rhoL - rhoV);

return f_of_rho*nut;
}


Do you know, where the problem could possibly be?
Thank you.