omegaWallFunction kinetic energy production term inconsistent with cited articles
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 .
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 .
- OpenFOAM version : v2112
Change the G0 term to match the reference or add the reference from which the current relation was taken from.