34 labelList cellBoundaryFaceCount(epsilon_.size(), 0);
36 scalar yPlusLam = this->yPlusLam(kappa_.value(), E_.value());
43 const fvPatch& curPatch = patches[
patchi];
45 if (isA<wallFvPatch>(curPatch))
49 label faceCelli = curPatch.faceCells()[facei];
51 epsilon_[faceCelli] = 0.0;
61 const fvPatch& curPatch = patches[
patchi];
63 if (isA<wallFvPatch>(curPatch))
74 label faceCelli = curPatch.faceCells()[facei];
77 scalar Cmu25 =
pow(Cmu_[faceCelli], 0.25);
78 scalar Cmu75 =
pow(Cmu_[faceCelli], 0.75);
89 cellBoundaryFaceCount[faceCelli]++;
91 epsilon_[faceCelli] +=
92 Cmu75*
pow(k_[faceCelli], 1.5)
93 /(kappa_.value()*y_[
patchi][facei]);
98 (nutw[facei] + nuw[facei])
100 *Cmu25*
sqrt(k_[faceCelli])
101 /(kappa_.value()*y_[
patchi][facei])
102 - (nonlinearStress_[faceCelli] && gradU_[faceCelli]);
112 const fvPatch& curPatch = patches[
patchi];
114 if (isA<wallFvPatch>(curPatch))
118 label faceCelli = curPatch.faceCells()[facei];
120 epsilon_[faceCelli] /= cellBoundaryFaceCount[faceCelli];
121 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];