Commit f6ec1cdc authored by Andrew Heather's avatar Andrew Heather
Browse files

updates to incompressible wall functions

parent 711ef825
......@@ -170,7 +170,7 @@ tmp<scalarField> nutRoughWallFunctionFvPatchScalarField::calcNut() const
const scalar Cmu25 = pow(Cmu_, 0.25);
tmp<scalarField> tnutw(scalarField(patch().size(), 0.0));
tmp<scalarField> tnutw(new scalarField(patch().size(), 0.0));
scalarField& nutw = tnutw();
forAll(nutw, faceI)
......
......@@ -130,9 +130,7 @@ nutSpalartAllmarasStandardRoughWallFunctionFvPatchScalarField::calcNut() const
// The flow velocity at the adjacent cell centre
scalarField magUp = mag(Uw.patchInternalField() - Uw);
scalarField magFaceGradU = mag(Uw.snGrad());
tmp<scalarField> tnutw(scalarField(patch().size(), 0.0));
tmp<scalarField> tnutw(new scalarField(patch().size(), 0.0));
scalarField& nutw = tnutw();
if (roughnessHeight_ > 0.0)
......@@ -262,22 +260,11 @@ nutSpalartAllmarasStandardRoughWallFunctionFvPatchScalarField::yPlus() const
{
notImplemented
(
"nutSpalartAllmarasWallFunctionFvPatchScalarField::yPlus() const"
"nutSpalartAllmarasStandardRoughWallFunctionFvPatchScalarField::yPlus()"
"const"
);
return tmp<scalarField>(NULL);
/*
const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
const scalarField& y = rasModel.y()[patch().index()];
const volScalarField& kw = db().lookupObject<volScalarField>(kName_);
const scalarField& kwc =
kw.boundaryField()[patch.index()].patchInternalField()
const scalarField& nuw =
patch().lookupPatchField<volScalarField, scalar>(nuName_);
return pow(Cmu_, 0.25)*y*sqrt(kwc)/nuw;
*/
}
......
......@@ -150,22 +150,11 @@ nutSpalartAllmarasStandardWallFunctionFvPatchScalarField::yPlus() const
{
notImplemented
(
"nutSpalartAllmarasWallFunctionFvPatchScalarField::yPlus() const"
"nutSpalartAllmarasStandardWallFunctionFvPatchScalarField::yPlus() "
"const"
);
return tmp<scalarField>(NULL);
/*
const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
const scalarField& y = rasModel.y()[patch().index()];
const volScalarField& kw = db().lookupObject<volScalarField>(kName_);
const scalarField& kwc =
kw.boundaryField()[patch.index()].patchInternalField()
const scalarField& nuw =
patch().lookupPatchField<volScalarField, scalar>(nuName_);
return pow(Cmu_, 0.25)*y*sqrt(kwc)/nuw;
*/
}
......
......@@ -82,7 +82,7 @@ tmp<scalarField> nutSpalartAllmarasWallFunctionFvPatchScalarField::calcUTau
+ 1/E_*(fkUu - 1.0/6.0*kUu*sqr(kUu));
scalar df =
1.0*y[facei]/nuw[facei]
y[facei]/nuw[facei]
+ magUpara/sqr(ut)
+ 1/E_*kUu*fkUu/ut;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment