37 namespace incompressible
44 void omegaWallFunctionFvPatchScalarField::checkType()
46 if (!isA<wallFvPatch>(
patch()))
48 FatalErrorIn(
"omegaWallFunctionFvPatchScalarField::checkType()")
49 <<
"Invalid wall function specification" <<
nl
50 <<
" Patch type for patch " <<
patch().
name()
51 <<
" must be wall" <<
nl
52 <<
" Current patch type is " <<
patch().type() <<
nl <<
endl
69 GName_(
"RASModel::G"),
94 nutName_(ptf.nutName_),
132 UName_(owfpsf.UName_),
133 kName_(owfpsf.kName_),
134 GName_(owfpsf.GName_),
135 nuName_(owfpsf.nuName_),
136 nutName_(owfpsf.nutName_),
138 kappa_(owfpsf.kappa_),
140 beta1_(owfpsf.beta1_)
153 UName_(owfpsf.UName_),
154 kName_(owfpsf.kName_),
155 GName_(owfpsf.GName_),
156 nuName_(owfpsf.nuName_),
157 nutName_(owfpsf.nutName_),
159 kappa_(owfpsf.kappa_),
161 beta1_(owfpsf.beta1_)
179 const scalar Cmu25 =
pow(Cmu_, 0.25);
207 scalar omegaVis = 6.0*nuw[faceI]/(beta1_*
sqr(y[faceI]));
208 scalar omegaLog =
sqrt(k[faceCellI])/(Cmu25*kappa_*y[faceI]);
209 omega[faceCellI] =
sqrt(
sqr(omegaVis) +
sqr(omegaLog));
212 (nutw[faceI] + nuw[faceI])
214 *Cmu25*
sqrt(k[faceCellI])
227 writeEntryIfDifferent<word>(os,
"U",
"U", UName_);
228 writeEntryIfDifferent<word>(os,
"k",
"k", kName_);
229 writeEntryIfDifferent<word>(os,
"G",
"RASModel::G", GName_);
230 writeEntryIfDifferent<word>(os,
"nu",
"nu", nuName_);
231 writeEntryIfDifferent<word>(os,
"nut",
"nut", nutName_);