exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Update on turbulence viscosity limiter

Browse source
This commit is contained in:
Hrvoje Jasak 2015-04-07 13:30:54 +01:00
parent d1f8b09ecb
commit 4a12f53a63
2 changed files with 6 additions and 16 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
src/turbulenceModels/compressible/RAS/RASModel/RASModel.C
View file

@ -173,14 +173,9 @@ tmp<volScalarField> RASModel::muEff() const
new volScalarField("muEff", mut() + mu())
);
if (gMax(tmuEff().internalField()) > muRatio_*gMax(mu().internalField()))
{
Info<< "Limiting turbulence viscosity" << endl;
// Apply mut limiter
tmuEff().internalField() =
Foam::min(tmuEff().internalField(), muRatio_*mu().internalField());
}
// Apply mut limiter
tmuEff().internalField() =
Foam::min(tmuEff().internalField(), muRatio_*mu().internalField());
return tmuEff;
}

11
src/turbulenceModels/incompressible/RAS/RASModel/RASModel.C
View file

@ -168,14 +168,9 @@ tmp<volScalarField> RASModel::nuEff() const
new volScalarField("nuEff", nut() + nu())
);
if (gMax(tnuEff().internalField()) > nuRatio_*gMax(nu().internalField()))
{
Info<< "Limiting turbulence viscosity" << endl;
// Apply nut limiter
tnuEff().internalField() =
Foam::min(tnuEff().internalField(), nuRatio_*nu().internalField());
}
// Apply nut limiter
tnuEff().internalField() =
Foam::min(tnuEff().internalField(), nuRatio_*nu().internalField());
return tnuEff;
}