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();
13 const fvPatch& curPatch = patches[
patchi];
15 if (isA<wallFvPatch>(curPatch))
19 label faceCelli = curPatch.faceCells()[facei];
21 epsilon[faceCelli] = 0.0;
31 const fvPatch& curPatch = patches[
patchi];
33 if (isA<wallFvPatch>(curPatch))
45 label faceCelli = curPatch.faceCells()[facei];
49 /(muw[facei]/rhow[facei]);
55 cellBoundaryFaceCount[faceCelli]++;
58 Cmu75*
rho[faceCelli]*
::pow(
k[faceCelli], 1.5)
64 mutw[facei]*magFaceGradU[facei]
77 const fvPatch& curPatch = patches[
patchi];
79 if (isA<wallFvPatch>(curPatch))
83 label faceCelli = curPatch.faceCells()[facei];
85 epsilon[faceCelli] /= cellBoundaryFaceCount[faceCelli];
86 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];