37 namespace compressible
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
60 omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
70 GName_(
"RASModel::G"),
82 omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
92 rhoName_(ptf.rhoName_),
96 mutName_(ptf.mutName_),
106 omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
129 omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
135 UName_(owfpsf.UName_),
136 rhoName_(owfpsf.rhoName_),
137 kName_(owfpsf.kName_),
138 GName_(owfpsf.GName_),
139 muName_(owfpsf.muName_),
140 mutName_(owfpsf.mutName_),
142 kappa_(owfpsf.kappa_),
144 beta1_(owfpsf.beta1_)
150 omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
157 UName_(owfpsf.UName_),
158 rhoName_(owfpsf.rhoName_),
159 kName_(owfpsf.kName_),
160 GName_(owfpsf.GName_),
161 muName_(owfpsf.muName_),
162 mutName_(owfpsf.mutName_),
164 kappa_(owfpsf.kappa_),
166 beta1_(owfpsf.beta1_)
174 void omegaWallFunctionFvPatchScalarField::updateCoeffs()
184 const scalar Cmu25 =
pow(Cmu_, 0.25);
190 (db().lookupObject<
volScalarField>(dimensionedInternalField().name()));
211 label faceCellI = patch().faceCells()[faceI];
213 scalar omegaVis = 6.0*(muw[faceI]/rhow[faceI])/(beta1_*
sqr(y[faceI]));
214 scalar omegaLog =
sqrt(k[faceCellI])/(Cmu25*kappa_*y[faceI]);
215 omega[faceCellI] =
sqrt(
sqr(omegaVis) +
sqr(omegaLog));
218 (mutw[faceI] + muw[faceI])
220 *Cmu25*
sqrt(k[faceCellI])
230 void omegaWallFunctionFvPatchScalarField::write(
Ostream& os)
const
233 writeEntryIfDifferent<word>(os,
"U",
"U", UName_);
234 writeEntryIfDifferent<word>(os,
"rho",
"rho", rhoName_);
235 writeEntryIfDifferent<word>(os,
"k",
"k", kName_);
236 writeEntryIfDifferent<word>(os,
"G",
"RASModel::G", GName_);
237 writeEntryIfDifferent<word>(os,
"mu",
"mu", muName_);
238 writeEntryIfDifferent<word>(os,
"mut",
"mut", mutName_);
243 writeEntry(
"value", os);