37 namespace compressible
44 void epsilonWallFunctionFvPatchScalarField::checkType()
46 if (!isA<wallFvPatch>(patch()))
48 FatalErrorIn(
"epsilonWallFunctionFvPatchScalarField::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 epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
69 GName_(
"RASModel::G"),
80 epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
93 mutName_(ptf.mutName_),
102 epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
123 epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
129 UName_(ewfpsf.UName_),
130 kName_(ewfpsf.kName_),
131 GName_(ewfpsf.GName_),
132 muName_(ewfpsf.muName_),
133 mutName_(ewfpsf.mutName_),
135 kappa_(ewfpsf.kappa_),
142 epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
149 UName_(ewfpsf.UName_),
150 kName_(ewfpsf.kName_),
151 GName_(ewfpsf.GName_),
152 muName_(ewfpsf.muName_),
153 mutName_(ewfpsf.mutName_),
155 kappa_(ewfpsf.kappa_),
164 void epsilonWallFunctionFvPatchScalarField::updateCoeffs()
173 const scalar Cmu25 =
pow(Cmu_, 0.25);
174 const scalar Cmu75 =
pow(Cmu_, 0.75);
182 (db().lookupObject<
volScalarField>(dimensionedInternalField().name()));
200 label faceCellI = patch().faceCells()[faceI];
202 epsilon[faceCellI] = Cmu75*
pow(k[faceCellI], 1.5)/(kappa_*y[faceI]);
205 (mutw[faceI] + muw[faceI])
207 *Cmu25*
sqrt(k[faceCellI])
217 void epsilonWallFunctionFvPatchScalarField::evaluate
226 void epsilonWallFunctionFvPatchScalarField::write(
Ostream& os)
const
229 writeEntryIfDifferent<word>(os,
"U",
"U", UName_);
230 writeEntryIfDifferent<word>(os,
"k",
"k", kName_);
231 writeEntryIfDifferent<word>(os,
"G",
"RASModel::G", GName_);
232 writeEntryIfDifferent<word>(os,
"mu",
"mu", muName_);
233 writeEntryIfDifferent<word>(os,
"mut",
"mut", mutName_);
237 writeEntry(
"value", os);