[vkrastev]

Hi all, actually I'm not sure if this is a real bug or not, but I have some questions about the Launder-Sharma model implementation in OpenFOAM (1.6). The literature about this model tells me that, in order to avoid a singular behavior of the sink term in the epsilon equation, the epsilon itself could be replaced with a quantity called epsilonTilda, which is equal to:

epsilonTilda=epsilon-2*nu*magSqr(grad(sqrt(k))

Well, looking at the source code, it seems like the epsilon equation is solved for epsilonTilda, and after that the "corrected" epsilon (epsilonTilda + 2*nu*...) is used as a sink term for the k-equation. However, in all other formulas or equations where epsilon is needed (in particular the turbulent eddy viscosity formula and the calculation of Ret for the damping functions), epsilonTilda is used instead...Maybe I'm missing something, but this sounds not so consistent to me as, in general, epsilon and epsilonTilda are not the same. Can anybody correct me or give some explanation for this?

Thank you all

V.