5 labelList cellBoundaryFaceCount(epsilon_.size(), 0);
8 scalar Cmu75 =
pow(CmuWall_.value(), 0.75);
15 const fvPatch& curPatch = patches[
patchi];
17 if (isA<wallFvPatch>(curPatch))
21 label faceCelli = curPatch.faceCells()[facei];
23 epsilon_[faceCelli] = 0.0;
30 const fvPatch& curPatch = patches[
patchi];
32 if (isA<wallFvPatch>(curPatch))
36 label faceCelli = curPatch.faceCells()[facei];
41 cellBoundaryFaceCount[faceCelli]++;
43 epsilon_[faceCelli] +=
44 Cmu75*
pow(k_[faceCelli], 1.5)
46 kappa_.value()*y_[faceCelli]
47 *(1.0 -
exp(-Aepsilon_.value()*yStar_[faceCelli]))
49 *
exp(-Amu_.value()*
sqr(yStar_[faceCelli]));
59 const fvPatch& curPatch = patches[
patchi];
61 if (isA<wallFvPatch>(curPatch))
65 label faceCelli = curPatch.faceCells()[facei];
67 epsilon_[faceCelli] /= cellBoundaryFaceCount[faceCelli];