[hiep.nguyentrong]

Hi all,

Sorry for stupid question. Im new to OpenFOAM, when i read source code of epsilonWallfunction. But i dont know what is the G field.
Can anyone explain it for me.
Thank you.

231 // Set epsilon and G
232 forAll(nutw, facei)
233 {
234 label celli = patch.faceCells()[facei];
235
236 scalar w = cornerWeights[facei];
237
238 epsilon[celli] += w*Cmu75*pow(k[celli], 1.5)/(kappa_*y[facei]);
239
240 G[celli] +=
241 w
242 *(nutw[facei] + nuw[facei])
243 *magGradUw[facei]
244 *Cmu25*sqrt(k[celli])
245 /(kappa_*y[facei]);
246 }
247 }

[HPE]

afaik: turbulent kinetic energy production rate [m2/s3] (the same units of epsilon).

[hiep.nguyentrong]

Originally Posted by HPE

afaik: turbulent kinetic energy production rate [m2/s3] (the same units of epsilon).

Thank you for reply.
Do you have any document or any key word for this. I still don't know why need to have this field

[Tobermory]

Perhaps read up a little on RAS closure models - start with the k-epsilon model, https://cpp.openfoam.org/v8/kEpsilon_8C_source.html. You will quickly see why this term is useful.