Hi,

Does anyone have a DEFINE_SOURCE UDF that implements this correction?

I have given it a go (see below) but the solution either diverges or becomes essentially laminar (no mixing at all...). Can anyone see any errors?

DEFINE_SOURCE(pope_correction,c,t,dS,eqn) { double C_eps_3; double chi, radius, source; double x[ND_ND];

C_eps_3 = 0.79;

C_CENTROID(x,c,t); radius = x[1];

chi = (C_DUDY(c,t) - C_DVDX(c,t)) * (C_DUDY(c,t) - C_DVDX(c,t)); chi = (C_U(c,t) / radius) * chi; chi = C_K(c,t) * C_K(c,t) / C_D(c,t) * chi; chi = 0.25 * C_eps_3 * C_R(c,t) * chi;

source = chi; dS[eqn] = 0.0;

return source; }

[xhaza]

Been looking for Pope correction using UDF aswell