Commit 1b94025a authored by sergio's avatar sergio
Browse files

BUG: adding limiter for KK to dynamicKEqn. Fixing access to tmp of nuEff.

parent 63ecd073
......@@ -48,7 +48,13 @@ volScalarField dynamicKEqn<BasicTurbulenceModel>::Ck
const volSymmTensorField MM
(
simpleFilter_(-2.0*this->delta()*sqrt(KK)*filter_(D))
simpleFilter_
(
-2.0*this->delta()*sqrt
(
max(KK, dimensionedScalar("zero", KK.dimensions(), 0.0))
)*filter_(D)
)
);
const volScalarField Ck
......
......@@ -109,7 +109,7 @@ void Foam::fixedShearStressFvPatchVectorField::updateCoeffs()
)
);
scalarField nuEff(turbModel.nuEff()()[patch().index()]);
tmp<scalarField> nuEff(turbModel.nuEff(patch().index()));
const vectorField Uc(patchInternalField());
......@@ -117,7 +117,7 @@ void Foam::fixedShearStressFvPatchVectorField::updateCoeffs()
const scalarField& ry = patch().deltaCoeffs();
operator==(tauHat*(tauHat & (tau0_*(1.0/(ry*nuEff)) + Uc)));
operator==(tauHat*(tauHat & (tau0_*(1.0/(ry*nuEff())) + Uc)));
fixedValueFvPatchVectorField::updateCoeffs();
}
......
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