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 .

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 .

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.

assigned to @kuti

added bug community labels

Hello @macakpav - many thanks for the ticket.

G0 term was not taken from that reference. Also, I'm not sure there is any resemblance between the code you quoted and the governing equation you posted. Could you elaborate it a bit?

Thank you for the reply.

Looking at the quoted code, the resemblance is the following (excluding w corner weighting coefficient):

(nutw[facei] + nuw[facei]) *magGradUw[facei] is , since

, where

, hence

and therefore , which is the quoted part of code.

Cmu25*sqrt(k[celli]) is the definition of .

That is why I think the code for computing the kinetic energy production (denoted as P in the quoted article) can be interpreted as .

If this relation for kinetic energy production term is not taken from this reference (and the other two articles cited in the source code do not elaborate on this), where is it from?