Skip to content

omegaWallFunction kinetic energy production term inconsistent with cited articles

Summary

This reference

Exponential/Max blending of the viscous and inertial sublayers (tag:PH): Popovac, M., & Hanjalić, K. (2007). Compound wall treatment for RANS computation of complex turbulent flows and heat transfer. Flow, turbulence and combustion, 78(2), 177-202. DOI:10.1007/s10494-006-9067-x

from omegaWallFunctionFvPatchScalarField.H recommends on page 193 that the production term G0 should be image.

However in code omegaWallFunctionFvPatchScalarField.C

G0[celli] += w *(nutw[facei] + nuw[facei]) *magGradUw[facei] *Cmu25*sqrt(k[celli]) /(nutw.kappa()*y[facei]);

the term for G0 corresponds rather to image.

Environment information

  • OpenFOAM version : v2112

Possible fixes

Change the G0 term to match the reference or add the reference from which the current relation was taken from.