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& rhow = rho_.boundaryField()[patchi];
00069
00070 const scalarField muw = mu().boundaryField()[patchi];
00071 const scalarField& mutw = mut_.boundaryField()[patchi];
00072
00073 scalarField magFaceGradU = mag(U_.boundaryField()[patchi].snGrad());
00074
00075 forAll(curPatch, facei)
00076 {
00077 label faceCelli = curPatch.faceCells()[facei];
00078
00079 scalar yPlus =
00080 Cmu25*turbulenceModel::y_[patchi][facei]*sqrt(k_[faceCelli])
00081 /(muw[facei]/rhow[facei]);
00082
00083
00084
00085
00086
00087 cellBoundaryFaceCount[faceCelli]++;
00088
00089 epsilon_[faceCelli] +=
00090 Cmu75*pow(k_[faceCelli], 1.5)
00091 /(kappa_*turbulenceModel::y_[patchi][facei]);
00092
00093 if (yPlus > yPlusLam_)
00094 {
00095 G[faceCelli] +=
00096 mutw[facei]*magFaceGradU[facei]
00097
00098 *Cmu25*sqrt(k_[faceCelli])
00099 /(kappa_*turbulenceModel::y_[patchi][facei]);
00100 }
00101 }
00102 }
00103 }
00104
00105
00106
00107
00108 forAll(patches, patchi)
00109 {
00110 const fvPatch& curPatch = patches[patchi];
00111
00112 if (isType<wallFvPatch>(curPatch))
00113 {
00114 forAll(curPatch, facei)
00115 {
00116 label faceCelli = curPatch.faceCells()[facei];
00117
00118 epsilon_[faceCelli] /= cellBoundaryFaceCount[faceCelli];
00119 G[faceCelli] /= cellBoundaryFaceCount[faceCelli];
00120 }
00121 }
00122 }
00123 }
00124
00125
00126