[dinhanh]

Dear All!

I am trying to make UDF code with the following terms in the Realizable k-epsilon model: W=(Sij*Sjk*Ski)/S^3 and U* = sqrt(Sij*Sij+Oij*Oij).
Here, Sij and Oij are the mean strain and rotation-rate tensors.

For the 2-dimensional problem, these above terms can be expressed as:
W=(Sii+Sij+2Sij)/S^3
U* = sqrt(Sii^2+Sjj^2+2*Sij^2+Oij*Oij)

with Sii=0.5*(C_DUDX(c,t)+C_DUDX(c,t));
Sij=0.5*(C_DUDY(c,t)+C_DVDX(c,t));
Sjj=0.5*(C_DVDY(c,t)+C_DVDY(c,t));
Oij=0.5*(C_DUDY(c,t)-C_DVDX(c,t));
S = sqrt(Sij*Sij)=sqrt(Sii*Sii+2*Sij*Sij+Sjj*Sjj)

Are these expressions correct? Please help.
Thank.