2 labelList cellBoundaryFaceCount(epsilon.size(), 0);
4 scalar Cmu25 =
::pow(
Cmu.value(), 0.25);
5 scalar Cmu75 =
::pow(
Cmu.value(), 0.75);
6 scalar kappa_ =
kappa.value();
7 scalar nub_ =
nub.value();
14 const fvPatch& currPatch = patches[
patchi];
16 if (isA<wallFvPatch>(currPatch))
20 label faceCelli = currPatch.faceCells()[facei];
22 epsilon[faceCelli] = 0.0;
32 const fvPatch& currPatch = patches[
patchi];
34 if (isA<wallFvPatch>(currPatch))
42 label faceCelli = currPatch.faceCells()[facei];
54 cellBoundaryFaceCount[faceCelli]++;
57 Cmu75*
::pow(
k[faceCelli], 1.5)
63 (nutbw[facei] + nub_)*magFaceGradU[facei]
64 *Cmu25*::
sqrt(
k[faceCelli])
76 const fvPatch& curPatch = patches[
patchi];
78 if (isA<wallFvPatch>(curPatch))
82 label faceCelli = curPatch.faceCells()[facei];
84 epsilon[faceCelli] /= cellBoundaryFaceCount[faceCelli];
85 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];