[nikesh]

I have been trying to implement a non equilibrium turbulent viscosity model (Akira Yoshizawa et al.) into OpenFOAM. The turbulence quantities (k and epsilon) are governed by the same equations of the standard k-epsilon model.

At the moment I am trying to validate my code for the case of a 2D flat plate turbulent boundary layer and compare the results I've obtained using LaunderSharmaKE model. I could obtain the results but they are not making sense to me. I have attached here the equations of the model and the code that I have used. I also have included the results.

I have as well made use of a low Reynolds-number corrections as proposed by Abe et al. I have included those equations too.

I would be really thankful if someone with insights in this field looked into it and shed me some light in this matter.

Thanks in advance!!
Nikesh.

[nikesh]

The result file is attached here.

What I have done here is that I have limited the blending functions.

1. min(fMu,1) and min(sqr(Rt), 50)
2. fMu-limited, Rt-not limited
3. fMu-not limited, Rt-limited