[Horus]

Hello,

I'm studying the implementation of the incompressible omegaWallFunction.

http://foam.sourceforge.net/docs/cpp/a09165_source.html

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

This clearly becomes oo when y approaches 0. Does this account for the turbulence damping near the wall?

Next question: What is G?

00173     volScalarField& G =
00174         const_cast<volScalarField&>
00175         (
00176             db().lookupObject<volScalarField>
00177             (
00178                 turbulence.type() + ".G"
00179             )
00180         );

And set here:

00213         G[faceCellI] =
00214             (nutw[faceI] + nuw[faceI])
00215            *magGradUw[faceI]
00216            *Cmu25*sqrt(k[faceCellI])
00217            /(kappa_*y[faceI]);
00218     }

The paper mentioned in the source code is available nowhere...

I'm especially clueless that the last part with the G is about?

Thanks!
Florian