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

TurbulenceModels: Correct nut during construction for single-phase solvers only.

Browse source 
For multiphase solvers the phase construction is not complete at this point.
This commit is contained in:
Henry Weller 2015-11-23 13:49:03 +00:00
parent 56913f6307
commit d5f8d691b0
23 changed files with 183 additions and 27 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LamBremhorstKE/LamBremhorstKE.C
View file

@ -180,8 +180,15 @@ LamBremhorstKE::LamBremhorstKE
if (type == typeName)
{
correctNut();
printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LienCubicKE/LienCubicKE.C
View file

@ -330,8 +330,15 @@ LienCubicKE::LienCubicKE
if (type == typeName)
{
correctNut();
printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/LienLeschziner/LienLeschziner.C
View file

@ -224,8 +224,15 @@ LienLeschziner::LienLeschziner
if (type == typeName)
{
correctNut();
printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/ShihQuadraticKE/ShihQuadraticKE.C
View file

@ -222,8 +222,15 @@ ShihQuadraticKE::ShihQuadraticKE
if (type == typeName)
{
correctNut();
printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C
View file

@ -207,8 +207,15 @@ qZeta::qZeta
if (type == typeName)
{
correctNut();
printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

4
src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
View file

@ -106,10 +106,10 @@ LaheyKEpsilon<BasicTurbulenceModel>::LaheyKEpsilon
{
if (type == typeName)
{
this->printCoeffs(type);
// Cannot correct nut yet: construction of the phases is not complete
// correctNut();
this->printCoeffs(type);
}
}

11
src/TurbulenceModels/turbulenceModels/LES/DeardorffDiffStress/DeardorffDiffStress.C
View file

@ -110,9 +110,16 @@ DeardorffDiffStress<BasicTurbulenceModel>::DeardorffDiffStress
{
if (type == typeName)
{
this->boundNormalStress(this->R_);
correctNut();
this->printCoeffs(type);
this->boundNormalStress(this->R_);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/Smagorinsky/Smagorinsky.C
View file

@ -113,8 +113,15 @@ Smagorinsky<BasicTurbulenceModel>::Smagorinsky
{
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/SpalartAllmarasDES/SpalartAllmarasDES.C
View file

@ -318,8 +318,15 @@ SpalartAllmarasDES<BasicTurbulenceModel>::SpalartAllmarasDES
{
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/WALE/WALE.C
View file

@ -143,8 +143,15 @@ WALE<BasicTurbulenceModel>::WALE
{
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/dynamicKEqn/dynamicKEqn.C
View file

@ -185,8 +185,15 @@ dynamicKEqn<BasicTurbulenceModel>::dynamicKEqn
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C
View file

@ -122,8 +122,15 @@ dynamicLagrangian<BasicTurbulenceModel>::dynamicLagrangian
{
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/LES/kEqn/kEqn.C
View file

@ -113,8 +113,15 @@ kEqn<BasicTurbulenceModel>::kEqn
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

12
src/TurbulenceModels/turbulenceModels/RAS/LRR/LRR.C
View file

@ -200,11 +200,19 @@ LRR<BasicTurbulenceModel>::LRR
{
if (type == typeName)
{
this->printCoeffs(type);
this->boundNormalStress(this->R_);
bound(epsilon_, this->epsilonMin_);
k_ = 0.5*tr(this->R_);
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/LaunderSharmaKE/LaunderSharmaKE.C
View file

@ -204,8 +204,15 @@ LaunderSharmaKE<BasicTurbulenceModel>::LaunderSharmaKE
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/RNGkEpsilon/RNGkEpsilon.C
View file

@ -205,8 +205,15 @@ RNGkEpsilon<BasicTurbulenceModel>::RNGkEpsilon
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

12
src/TurbulenceModels/turbulenceModels/RAS/SSG/SSG.C
View file

@ -209,11 +209,19 @@ SSG<BasicTurbulenceModel>::SSG
{
if (type == typeName)
{
this->printCoeffs(type);
this->boundNormalStress(this->R_);
bound(epsilon_, this->epsilonMin_);
k_ = 0.5*tr(this->R_);
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/SpalartAllmaras/SpalartAllmaras.C
View file

@ -255,8 +255,15 @@ SpalartAllmaras<BasicTurbulenceModel>::SpalartAllmaras
{
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C
View file

@ -188,7 +188,14 @@ kEpsilon<BasicTurbulenceModel>::kEpsilon
if (type == typeName)
{
this->printCoeffs(type);
correctNut();
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/kOmega/kOmega.C
View file

@ -147,8 +147,15 @@ kOmega<BasicTurbulenceModel>::kOmega
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/kOmegaSST/kOmegaSST.C
View file

@ -350,8 +350,15 @@ kOmegaSST<BasicTurbulenceModel>::kOmegaSST
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/realizableKE/realizableKE.C
View file

@ -226,8 +226,15 @@ realizableKE<BasicTurbulenceModel>::realizableKE
if (type == typeName)
{
correctNut();
this->printCoeffs(type);
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}

9
src/TurbulenceModels/turbulenceModels/RAS/v2f/v2f.C
View file

@ -239,7 +239,14 @@ v2f<BasicTurbulenceModel>::v2f
if (type == typeName)
{
this->printCoeffs(type);
correctNut();
// Correct nut for single-phase solvers only.
// For multiphase solvers the phase construction is not complete
// at this point.
if (isType<geometricOneField>(alpha))
{
correctNut();
}
}
}