omegaWallFunctionFvPatchScalarField in OpenFoam 3.0.0
 

Hi guys, I'm using a k omega SST with a mesh with y+<1 and so that is able to resolve the viscous sub-layer. I've founded a discussion http://www.cfd-online.com/Forums/ope...tml#post438465 in which is reported a modification of omegawallfunction's:

scalar omegaVis = 6.0*nuw[faceI]/(beta1_*sqr(y[faceI]));
scalar omegaLog = sqrt(k[faceCellI])/(Cmu25*kappa_*y[faceI]);
omega[faceCellI] = sqrt(sqr(omegaVis) + sqr(omegaLog));

we have found cases for which this causes a sudden change in the viscosity near the wall if the mesh is sufficiently fine and that just using the logarithmic part give more continuous behavior:

omega[faceCellI] = omegaLog;

So, according you, i have to modify the code like above, to allow to model k omega SST to work correctly in viscous sub-layer?