zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

tidying headers + assignments to turb fields

Browse Source
This commit is contained in:
andy
2009-07-10 17:53:12 +01:00
parent 2158923839
commit 26e04aab19
14 changed files with 36 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/turbulenceModels/compressible/RAS/LRR/LRR.C
View File

@ -229,10 +229,10 @@ LRR::LRR
<< exit(FatalError); << exit(FatalError);
} }
mut_ == Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -418,7 +418,7 @@ void LRR::correct()
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/epsilon_; mut_ = rho_*Cmu_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

10
src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
View File

@ -27,8 +27,6 @@ License
#include "LaunderGibsonRSTM.H" #include "LaunderGibsonRSTM.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H" #include "wallFvPatch.H"
#include "wallDist.H"
#include "wallDistReflection.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -251,10 +249,10 @@ LaunderGibsonRSTM::LaunderGibsonRSTM
<< exit(FatalError); << exit(FatalError);
} }
mut_ == Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -349,7 +347,7 @@ void LaunderGibsonRSTM::correct()
if (!turbulence_) if (!turbulence_)
{ {
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity
@ -457,7 +455,7 @@ void LaunderGibsonRSTM::correct()
// Re-calculate turbulent viscosity // Re-calculate turbulent viscosity
mut_ == Cmu_*rho_*sqr(k_)/epsilon_; mut_ = Cmu_*rho_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

9
src/turbulenceModels/compressible/RAS/RNGkEpsilon/RNGkEpsilon.C
View File

@ -26,7 +26,6 @@ License
#include "RNGkEpsilon.H" #include "RNGkEpsilon.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -187,10 +186,10 @@ RNGkEpsilon::RNGkEpsilon
autoCreateAlphat("alphat", mesh_) autoCreateAlphat("alphat", mesh_)
) )
{ {
mut_ == Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -277,7 +276,7 @@ void RNGkEpsilon::correct()
if (!turbulence_) if (!turbulence_)
{ {
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity
@ -349,7 +348,7 @@ void RNGkEpsilon::correct()
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/epsilon_; mut_ = rho_*Cmu_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

9
src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C
View File

@ -26,7 +26,6 @@ License
#include "kEpsilon.H" #include "kEpsilon.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -169,10 +168,10 @@ kEpsilon::kEpsilon
autoCreateAlphat("alphat", mesh_) autoCreateAlphat("alphat", mesh_)
) )
{ {
mut_ == Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = Cmu_*rho_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -258,7 +257,7 @@ void kEpsilon::correct()
if (!turbulence_) if (!turbulence_)
{ {
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = rho_*Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity
@ -323,7 +322,7 @@ void kEpsilon::correct()
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rho_*Cmu_*sqr(k_)/epsilon_; mut_ = rho_*Cmu_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

9
src/turbulenceModels/compressible/RAS/kOmegaSST/kOmegaSST.C
View File

@ -26,7 +26,6 @@ License
#include "kOmegaSST.H" #include "kOmegaSST.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -259,10 +258,10 @@ kOmegaSST::kOmegaSST
autoCreateAlphat("alphat", mesh_) autoCreateAlphat("alphat", mesh_)
) )
{ {
mut_ == a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_))))); mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_)))));
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -352,7 +351,7 @@ void kOmegaSST::correct()
if (!turbulence_) if (!turbulence_)
{ {
// Re-calculate viscosity // Re-calculate viscosity
mut_ == mut_ =
a1_*rho_*k_ a1_*rho_*k_
/max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_))))); /max(a1_*omega_, F2()*sqrt(magSqr(symm(fvc::grad(U_)))));
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
@ -435,7 +434,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity // Re-calculate viscosity
mut_ == a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(S2)); mut_ = a1_*rho_*k_/max(a1_*omega_, F2()*sqrt(S2));
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

9
src/turbulenceModels/compressible/RAS/realizableKE/realizableKE.C
View File

@ -26,7 +26,6 @@ License
#include "realizableKE.H" #include "realizableKE.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -204,10 +203,10 @@ realizableKE::realizableKE
bound(k_, k0_); bound(k_, k0_);
bound(epsilon_, epsilon0_); bound(epsilon_, epsilon0_);
mut_ == rCmu(fvc::grad(U_))*rho_*sqr(k_)/(epsilon_ + epsilonSmall_); mut_ = rCmu(fvc::grad(U_))*rho_*sqr(k_)/(epsilon_ + epsilonSmall_);
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
alphat_ == mut_/Prt_; alphat_ = mut_/Prt_;
alphat_.correctBoundaryConditions(); alphat_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -291,7 +290,7 @@ void realizableKE::correct()
if (!turbulence_) if (!turbulence_)
{ {
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rCmu(fvc::grad(U_))*rho_*sqr(k_)/epsilon_; mut_ = rCmu(fvc::grad(U_))*rho_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity
@ -362,7 +361,7 @@ void realizableKE::correct()
bound(k_, k0_); bound(k_, k0_);
// Re-calculate viscosity // Re-calculate viscosity
mut_ == rCmu(gradU, S2, magS)*rho_*sqr(k_)/epsilon_; mut_ = rCmu(gradU, S2, magS)*rho_*sqr(k_)/epsilon_;
mut_.correctBoundaryConditions(); mut_.correctBoundaryConditions();
// Re-calculate thermal diffusivity // Re-calculate thermal diffusivity

4
src/turbulenceModels/incompressible/RAS/LRR/LRR.C
View File

@ -186,7 +186,7 @@ LRR::LRR
autoCreateNut("nut", mesh_) autoCreateNut("nut", mesh_)
) )
{ {
nut_ == Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_); nut_ = Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_);
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0) if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
@ -380,7 +380,7 @@ void LRR::correct()
// Re-calculate viscosity // Re-calculate viscosity
nut_ == Cmu_*sqr(k_)/epsilon_; nut_ = Cmu_*sqr(k_)/epsilon_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();

2
src/turbulenceModels/incompressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.C
View File

@ -421,7 +421,7 @@ void LaunderGibsonRSTM::correct()
// Re-calculate turbulent viscosity // Re-calculate turbulent viscosity
nut_ == Cmu_*sqr(k_)/epsilon_; nut_ = Cmu_*sqr(k_)/epsilon_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();

5
src/turbulenceModels/incompressible/RAS/LienCubicKE/LienCubicKE.C
View File

@ -26,7 +26,6 @@ License
#include "LienCubicKE.H" #include "LienCubicKE.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -227,7 +226,7 @@ LienCubicKE::LienCubicKE
) )
) )
{ {
nut_ == Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_) + C5viscosity_; nut_ = Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_) + C5viscosity_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -382,7 +381,7 @@ void LienCubicKE::correct()
- 2.0*pow(Cmu_, 3.0)*pow(k_, 4.0)/pow(epsilon_, 3.0) - 2.0*pow(Cmu_, 3.0)*pow(k_, 4.0)/pow(epsilon_, 3.0)
*(magSqr(gradU_ + gradU_.T()) - magSqr(gradU_ - gradU_.T())); *(magSqr(gradU_ + gradU_.T()) - magSqr(gradU_ - gradU_.T()));
nut_ == Cmu_*sqr(k_)/epsilon_ + C5viscosity_; nut_ = Cmu_*sqr(k_)/epsilon_ + C5viscosity_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
nonlinearStress_ = symm nonlinearStress_ = symm

5
src/turbulenceModels/incompressible/RAS/RNGkEpsilon/RNGkEpsilon.C
View File

@ -26,7 +26,6 @@ License
#include "RNGkEpsilon.H" #include "RNGkEpsilon.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -156,7 +155,7 @@ RNGkEpsilon::RNGkEpsilon
autoCreateNut("nut", mesh_) autoCreateNut("nut", mesh_)
) )
{ {
nut_ == Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_); nut_ = Cmu_*sqr(k_)/(epsilon_ + epsilonSmall_);
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -295,7 +294,7 @@ void RNGkEpsilon::correct()
// Re-calculate viscosity // Re-calculate viscosity
nut_ == Cmu_*sqr(k_)/epsilon_; nut_ = Cmu_*sqr(k_)/epsilon_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
} }

1
src/turbulenceModels/incompressible/RAS/kEpsilon/kEpsilon.C
View File

@ -26,7 +26,6 @@ License
#include "kEpsilon.H" #include "kEpsilon.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"

5
src/turbulenceModels/incompressible/RAS/kOmega/kOmega.C
View File

@ -26,7 +26,6 @@ License
#include "kOmega.H" #include "kOmega.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -138,7 +137,7 @@ kOmega::kOmega
autoCreateNut("nut", mesh_) autoCreateNut("nut", mesh_)
) )
{ {
nut_ == k_/(omega_ + omegaSmall_); nut_ = k_/(omega_ + omegaSmall_);
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -268,7 +267,7 @@ void kOmega::correct()
// Re-calculate viscosity // Re-calculate viscosity
nut_ == k_/(omega_ + omegaSmall_); nut_ = k_/(omega_ + omegaSmall_);
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
} }

5
src/turbulenceModels/incompressible/RAS/kOmegaSST/kOmegaSST.C
View File

@ -26,7 +26,6 @@ License
#include "kOmegaSST.H" #include "kOmegaSST.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -237,7 +236,7 @@ kOmegaSST::kOmegaSST
autoCreateNut("nut", mesh_) autoCreateNut("nut", mesh_)
) )
{ {
nut_ == nut_ =
a1_*k_ a1_*k_
/max /max
( (
@ -395,7 +394,7 @@ void kOmegaSST::correct()
// Re-calculate viscosity // Re-calculate viscosity
nut_ == a1_*k_/max(a1_*omega_, F2()*sqrt(S2)); nut_ = a1_*k_/max(a1_*omega_, F2()*sqrt(S2));
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
} }

5
src/turbulenceModels/incompressible/RAS/realizableKE/realizableKE.C
View File

@ -26,7 +26,6 @@ License
#include "realizableKE.H" #include "realizableKE.H"
#include "addToRunTimeSelectionTable.H" #include "addToRunTimeSelectionTable.H"
#include "wallFvPatch.H"
#include "backwardsCompatibilityWallFunctions.H" #include "backwardsCompatibilityWallFunctions.H"
@ -182,7 +181,7 @@ realizableKE::realizableKE
bound(k_, k0_); bound(k_, k0_);
bound(epsilon_, epsilon0_); bound(epsilon_, epsilon0_);
nut_ == rCmu(fvc::grad(U_))*sqr(k_)/(epsilon_ + epsilonSmall_); nut_ = rCmu(fvc::grad(U_))*sqr(k_)/(epsilon_ + epsilonSmall_);
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
printCoeffs(); printCoeffs();
@ -324,7 +323,7 @@ void realizableKE::correct()
// Re-calculate viscosity // Re-calculate viscosity
nut_ == rCmu(gradU, S2, magS)*sqr(k_)/epsilon_; nut_ = rCmu(gradU, S2, magS)*sqr(k_)/epsilon_;
nut_.correctBoundaryConditions(); nut_.correctBoundaryConditions();
} }