00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 {
00034 labelList cellBoundaryFaceCount(epsilon_.size(), 0);
00035
00036 scalar Cmu25 = pow(Cmu.value(), 0.25);
00037 scalar Cmu75 = pow(Cmu.value(), 0.75);
00038
00039 const fvPatchList& patches = mesh_.boundary();
00040
00041
00042 forAll(patches, patchi)
00043 {
00044 const fvPatch& curPatch = patches[patchi];
00045
00046 if (isType<wallFvPatch>(curPatch))
00047 {
00048 forAll(curPatch, facei)
00049 {
00050 label faceCelli = curPatch.faceCells()[facei];
00051
00052 epsilon_[faceCelli] = 0.0;
00053 G[faceCelli] = 0.0;
00054 }
00055 }
00056 }
00057
00058
00059
00060 forAll(patches, patchi)
00061 {
00062 const fvPatch& curPatch = patches[patchi];
00063
00064 if (isType<wallFvPatch>(curPatch))
00065 {
00066 # include "checkPatchFieldTypes.H"
00067
00068 const scalarField& nuw = nu().boundaryField()[patchi];
00069 const scalarField& nutw = nut_.boundaryField()[patchi];
00070
00071 scalarField magFaceGradU = mag(U_.boundaryField()[patchi].snGrad());
00072
00073 forAll(curPatch, facei)
00074 {
00075 label faceCelli = curPatch.faceCells()[facei];
00076
00077 scalar yPlus =
00078 Cmu25*turbulenceModel::y_[patchi][facei]
00079 *sqrt(k_[faceCelli])
00080 /nuw[facei];
00081
00082
00083
00084
00085
00086 cellBoundaryFaceCount[faceCelli]++;
00087
00088 epsilon_[faceCelli] +=
00089 Cmu75*pow(k_[faceCelli], 1.5)
00090 /(kappa_*turbulenceModel::y_[patchi][facei]);
00091
00092 if (yPlus > yPlusLam_)
00093 {
00094 G[faceCelli] +=
00095 nutw[facei]
00096 *magFaceGradU[facei]
00097 *Cmu25*sqrt(k_[faceCelli])
00098 /(kappa_*turbulenceModel::y_[patchi][facei]);
00099 }
00100 }
00101 }
00102 }
00103
00104
00105
00106
00107 forAll(patches, patchi)
00108 {
00109 const fvPatch& curPatch = patches[patchi];
00110
00111 if (isType<wallFvPatch>(curPatch))
00112 {
00113 forAll(curPatch, facei)
00114 {
00115 label faceCelli = curPatch.faceCells()[facei];
00116
00117 epsilon_[faceCelli] /= cellBoundaryFaceCount[faceCelli];
00118 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];
00119 }
00120 }
00121 }
00122 }
00123
00124
00125