[immortality]

in mutkWallFunctionFvPatchScalarField the equation uses for y-plus is related to y-star not y-plus

Code:

y-star=pow(0.09,0.25)*sqr(k)*nearDist()/mu

Code:

y-plus=rho*u_star*y/mu which u*=aqrt(taw_w/rho)

y-star is related to turbulent kinetic energy while y-plus relates to friction velocity(that relates to wall skin friction in its own)

in the code of mut
:
Code:

forAll(mutw, faceI)     {         label faceCellI = patch().faceCells()[faceI];          scalar yPlus =             Cmu25*y[faceI]*sqrt(k[faceCellI])/(muw[faceI]/rhow[faceI]);          if (yPlus > yPlusLam_)         {             mutw[faceI] = muw[faceI]*(yPlus*kappa_/log(E_*yPlus) - 1);         }     }

thanks.