research/OpenFOAM-4.x-lab
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

buoyantKEpsilon: filter sources depending on mag(g) so that it operates as kEpsilon if mag(g) = 0

Browse source
This commit is contained in:
Henry 2015-01-28 16:34:18 +00:00
parent 3ecb02c4f0
commit 870d71db1f
1 changed files with 27 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

35
src/TurbulenceModels/compressible/RAS/buoyantKEpsilon/buoyantKEpsilon.C
View file

@ -116,7 +116,18 @@ template<class BasicTurbulenceModel>
tmp<fvScalarMatrix>
buoyantKEpsilon<BasicTurbulenceModel>::kSource() const
{
return -fvm::SuSp(Gcoef(), this->k_);
const uniformDimensionedVectorField& g =
this->mesh_.objectRegistry::template
lookupObject<uniformDimensionedVectorField>("g");
if (mag(g.value()) > SMALL)
{
return -fvm::SuSp(Gcoef(), this->k_);
}
else
{
return kEpsilon<BasicTurbulenceModel>::kSource();
}
}
@ -128,15 +139,23 @@ buoyantKEpsilon<BasicTurbulenceModel>::epsilonSource() const
this->mesh_.objectRegistry::template
lookupObject<uniformDimensionedVectorField>("g");
vector gHat(g.value()/mag(g.value()));
if (mag(g.value()) > SMALL)
{
vector gHat(g.value()/mag(g.value()));
volScalarField v(gHat & this->U_);
volScalarField u
(
mag(this->U_ - gHat*v) + dimensionedScalar("SMALL", dimVelocity, SMALL)
);
volScalarField v(gHat & this->U_);
volScalarField u
(
mag(this->U_ - gHat*v)
+ dimensionedScalar("SMALL", dimVelocity, SMALL)
);
return -fvm::SuSp(this->C1_*tanh(mag(v)/u)*Gcoef(), this->epsilon_);
return -fvm::SuSp(this->C1_*tanh(mag(v)/u)*Gcoef(), this->epsilon_);
}
else
{
return kEpsilon<BasicTurbulenceModel>::epsilonSource();
}
}