songtianlun/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

define values for standardCmu, standardKappa, standardE in turbulenceModel.H

Browse Source
This commit is contained in:
Mark Olesen 2009-08-01 09:36:34 +02:00
parent 827e3659b9
commit 54a8b6a152
25 changed files with 156 additions and 187 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
applications/solvers/combustion/PDRFoam/PDRModels/turbulence/PDRkEpsilon/PDRkEpsilon.H
View File

@ -27,7 +27,7 @@ Class
Description
Standard k-epsilon turbulence model with additional source terms
corresponding to PDR basic drag model (\link basic.H \endlink)
corresponding to PDR basic drag model (@link basic.H @endlink)
The default model coefficients correspond to the following:
@verbatim
@ -47,10 +47,10 @@ Description
\f$ \kappa-\epsilon \f$ equation for the generation of turbulence due to
interaction with unresolved obstacles.
In the \f$ \epsilon \f$ equation \f$ C_{1} G_{R} \f$ is added as a source
In the \f$ \epsilon \f$ equation, \f$ C_{1} G_{R} \f$ is added as a source
term.
In the \f$ \kappa \f$ equation \f$ G_{R} \f$ is added as a source term.
In the \f$ \kappa \f$ equation, \f$ G_{R} \f$ is added as a source term.
SourceFiles
PDRkEpsilon.C

2
applications/solvers/multiphase/bubbleFoam/kEpsilon.H
View File

@ -1,4 +1,4 @@
if(turbulence)
if (turbulence)
{
if (mesh.changing())
{

6
applications/utilities/postProcessing/turbulence/createTurbulenceFields/createTurbulenceFields.C
View File

@ -41,6 +41,10 @@ Source files:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// turbulence constants - file-scope
static const scalar Cmu(0.09);
int main(int argc, char *argv[])
{
timeSelector::addOptions();
@ -104,8 +108,6 @@ int main(int argc, char *argv[])
if (!IOobject("omega", runTime.timeName(), mesh).headerOk())
{
const scalar Cmu = 0.09;
Info<< "creating omega" << endl;
volScalarField omega
(

24
src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/alphaSgsWallFunctions/alphaSgsJayatillekeWallFunction/alphaSgsJayatillekeWallFunctionFvPatchScalarField.C
View File

@ -46,12 +46,6 @@ scalar alphaSgsJayatillekeWallFunctionFvPatchScalarField::maxExp_ = 50.0;
scalar alphaSgsJayatillekeWallFunctionFvPatchScalarField::tolerance_ = 0.01;
label alphaSgsJayatillekeWallFunctionFvPatchScalarField::maxIters_ = 10;
//! @cond fileScope
static const scalar defaultPrt_(0.85);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
@ -121,9 +115,9 @@ alphaSgsJayatillekeWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF),
Prt_(defaultPrt_),
kappa_(defaultKappa_),
E_(defaultE_)
Prt_(standardPrt),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{
checkType();
}
@ -154,9 +148,9 @@ alphaSgsJayatillekeWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF, dict),
Prt_(dict.lookupOrDefault<scalar>("Prt", defaultPrt_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
Prt_(dict.lookupOrDefault<scalar>("Prt", standardPrt)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{
checkType();
}
@ -318,9 +312,9 @@ void alphaSgsJayatillekeWallFunctionFvPatchScalarField::write(Ostream& os) const
{
fvPatchField<scalar>::write(os);
writeEntryIfDifferent<scalar>(os, "Prt", defaultPrt_, Prt_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Prt", standardPrt, Prt_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

6
src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/alphaSgsWallFunctions/alphaSgsJayatillekeWallFunction/alphaSgsJayatillekeWallFunctionFvPatchScalarField.H
View File

@ -93,6 +93,12 @@ class alphaSgsJayatillekeWallFunctionFvPatchScalarField
public:
// Public data
//- Standard value for turbulent Prandtl number %Prt
static const scalar standardPrt = 0.85;
//- Runtime type information
TypeName("alphaSgsJayatillekeWallFunction");

4
src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/alphaSgsWallFunctions/alphaSgsWallFunction/alphaSgsWallFunctionFvPatchScalarField.C
View File

@ -67,7 +67,7 @@ alphaSgsWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF),
Prt_(0.85)
Prt_(standardPrt)
{
checkType();
}
@ -96,7 +96,7 @@ alphaSgsWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF, dict),
Prt_(dict.lookupOrDefault<scalar>("Prt", 0.85))
Prt_(dict.lookupOrDefault<scalar>("Prt", standardPrt))
{
checkType();
}

6
src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/alphaSgsWallFunctions/alphaSgsWallFunction/alphaSgsWallFunctionFvPatchScalarField.H
View File

@ -70,6 +70,12 @@ class alphaSgsWallFunctionFvPatchScalarField
public:
// Public data
//- Standard value for turbulent Prandtl number %Prt
static const scalar standardPrt = 0.85;
//- Runtime type information
TypeName("alphaSgsWallFunction");

20
src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/muSgsWallFunctions/muSgsWallFunction/muSgsWallFunctionFvPatchScalarField.C
View File

@ -39,14 +39,6 @@ namespace compressible
namespace LESModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
@ -59,8 +51,8 @@ muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
UName_("U"),
rhoName_("rho"),
muName_("mu"),
kappa_(defaultKappa_),
E_(defaultE_)
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{}
@ -92,8 +84,8 @@ muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
UName_(dict.lookupOrDefault<word>("U", "U")),
rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
muName_(dict.lookupOrDefault<word>("mu", "mu")),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{}
@ -205,8 +197,8 @@ void muSgsWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "rho", "rho", rhoName_);
writeEntryIfDifferent<word>(os, "mu", "mu", muName_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

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

@ -183,7 +183,7 @@ LaunderGibsonRSTM::LaunderGibsonRSTM
)
),
y_(mesh_),
yr_(mesh_),
R_
(
@ -371,7 +371,7 @@ void LaunderGibsonRSTM::correct()
if (mesh_.changing())
{
y_.correct();
yr_.correct();
}
volSymmTensorField P = -twoSymm(R_ & fvc::grad(U_));
@ -436,10 +436,10 @@ void LaunderGibsonRSTM::correct()
// wall reflection terms
+ symm
(
I*((y_.n() & reflect) & y_.n())
- 1.5*(y_.n()*(reflect & y_.n())
+ (y_.n() & reflect)*y_.n())
)*pow(Cmu_, 0.75)*rho_*pow(k_, 1.5)/(kappa_*y_*epsilon_)
I*((yr_.n() & reflect) & yr_.n())
- 1.5*(yr_.n()*(reflect & yr_.n())
+ (yr_.n() & reflect)*yr_.n())
)*pow(Cmu_, 0.75)*rho_*pow(k_, 1.5)/(kappa_*yr_*epsilon_)
);
REqn().relax();

4
src/turbulenceModels/compressible/RAS/LaunderGibsonRSTM/LaunderGibsonRSTM.H
View File

@ -42,10 +42,10 @@ Description
C2Ref 0.3;
Cs 0.25;
Ceps 0.15;
sigmah 1.0; // only for compressible
sigmaEps 1.3;
sigmaR 0.81967;
couplingFactor 0.0;
Prt 1.0; // only for compressible
}
@endverbatim
@ -105,7 +105,7 @@ class LaunderGibsonRSTM
// Fields
wallDistReflection y_;
wallDistReflection yr_;
volSymmTensorField R_;
volScalarField k_;

1
src/turbulenceModels/compressible/RAS/LaunderSharmaKE/LaunderSharmaKE.H
View File

@ -40,6 +40,7 @@ Description
alphah 1.0; // only for compressible
alphahk 1.0; // only for compressible
alphaEps 0.76923;
Prt 1.0; // only for compressible
}
@endverbatim

7
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentMixingLengthDissipationRateInlet/turbulentMixingLengthDissipationRateInletFvPatchScalarField.C
View File

@ -110,8 +110,11 @@ void turbulentMixingLengthDissipationRateInletFvPatchScalarField::updateCoeffs()
// Lookup Cmu corresponding to the turbulence model selected
const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
const scalar Cmu =
rasModel.coeffDict().lookupOrDefault<scalar>("Cmu", 0.09);
const scalar Cmu = rasModel.coeffDict().lookupOrDefault<scalar>
(
"Cmu",
turbulenceModel::standardCmu
);
const scalar Cmu75 = pow(Cmu, 0.75);

13
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/turbulentMixingLengthFrequencyInlet/turbulentMixingLengthFrequencyInletFvPatchScalarField.C
View File

@ -38,12 +38,6 @@ namespace Foam
namespace compressible
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
//! @endcond fileScope
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
turbulentMixingLengthFrequencyInletFvPatchScalarField::
@ -121,8 +115,11 @@ void turbulentMixingLengthFrequencyInletFvPatchScalarField::updateCoeffs()
// Lookup Cmu corresponding to the turbulence model selected
const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
const scalar Cmu =
rasModel.coeffDict().lookupOrDefault<scalar>("Cmu", defaultCmu_);
const scalar Cmu = rasModel.coeffDict().lookupOrDefault<scalar>
(
"Cmu",
turbulenceModel::standardCmu
);
const scalar Cmu25 = pow(Cmu, 0.25);

12
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatWallFunction/alphatWallFunctionFvPatchScalarField.C
View File

@ -38,12 +38,6 @@ namespace compressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultPrt_(0.85);
//! @endcond fileScope
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
alphatWallFunctionFvPatchScalarField::
@ -55,7 +49,7 @@ alphatWallFunctionFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF),
mutName_("mut"),
Prt_(defaultPrt_)
Prt_(standardPrt)
{}
@ -84,7 +78,7 @@ alphatWallFunctionFvPatchScalarField
:
fixedValueFvPatchScalarField(p, iF, dict),
mutName_(dict.lookupOrDefault<word>("mut", "mut")),
Prt_(dict.lookupOrDefault<scalar>("Prt", defaultPrt_))
Prt_(dict.lookupOrDefault<scalar>("Prt", standardPrt))
{}
@ -128,7 +122,7 @@ void alphatWallFunctionFvPatchScalarField::write(Ostream& os) const
{
fvPatchField<scalar>::write(os);
writeEntryIfDifferent<word>(os, "mut", "mut", mutName_);
writeEntryIfDifferent<scalar>(os, "Prt", defaultPrt_, Prt_);
writeEntryIfDifferent<scalar>(os, "Prt", standardPrt, Prt_);
writeEntry("value", os);
}

7
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/alphatWallFunctions/alphatWallFunction/alphatWallFunctionFvPatchScalarField.H
View File

@ -28,6 +28,7 @@ Class
Description
Boundary condition for turbulent thermal diffusivity when using wall
functions
- replicates OpenFOAM v1.5 (and earlier) behaviour
- Turbulent Prandtl number defaults to 0.85 if unspecified
@ -69,6 +70,12 @@ class alphatWallFunctionFvPatchScalarField
public:
// Public data
//- Standard value for turbulent Prandtl number %Prt
static const scalar standardPrt = 0.85;
//- Runtime type information
TypeName("compressible::alphatWallFunction");

26
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace compressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void epsilonWallFunctionFvPatchScalarField::checkType()
@ -79,9 +71,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
rhoName_("rho"),
muName_("mu"),
mutName_("mut"),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_)
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{
checkType();
}
@ -124,9 +116,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
muName_(dict.lookupOrDefault<word>("mu", "mu")),
mutName_(dict.lookupOrDefault<word>("mut", "mut")),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{
checkType();
}
@ -255,9 +247,9 @@ void epsilonWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "mu", "mu", muName_);
writeEntryIfDifferent<word>(os, "mut", "mut", mutName_);
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

26
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/mutWallFunctions/mutWallFunction/mutWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace compressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
void mutWallFunctionFvPatchScalarField::checkType()
@ -115,9 +107,9 @@ tmp<scalarField> mutWallFunctionFvPatchScalarField::calcMut() const
void mutWallFunctionFvPatchScalarField::writeLocalEntries(Ostream& os) const
{
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
}
@ -130,9 +122,9 @@ mutWallFunctionFvPatchScalarField::mutWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_),
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE),
yPlusLam_(calcYPlusLam(kappa_, E_))
{}
@ -161,9 +153,9 @@ mutWallFunctionFvPatchScalarField::mutWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF, dict),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_)),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE)),
yPlusLam_(calcYPlusLam(kappa_, E_))
{}

26
src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace compressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void omegaWallFunctionFvPatchScalarField::checkType()
@ -79,9 +71,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
GName_("RASModel::G"),
muName_("mu"),
mutName_("mut"),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_)
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{
checkType();
}
@ -124,9 +116,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
GName_(dict.lookupOrDefault<word>("G", "RASModel::G")),
muName_(dict.lookupOrDefault<word>("mu", "mu")),
mutName_(dict.lookupOrDefault<word>("mut", "mut")),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{
checkType();
}
@ -244,9 +236,9 @@ void omegaWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "mu", "mu", muName_);
writeEntryIfDifferent<word>(os, "mut", "mut", mutName_);
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

12
src/turbulenceModels/compressible/turbulenceModel/turbulenceModel.H
View File

@ -98,6 +98,18 @@ private:
public:
// Public data
//- Standard value for turbulence constant %Cmu
static const scalar standardCmu = 0.09;
//- Standard value for turbulence constant %kappa
static const scalar standardKappa = 0.41;
//- Standard value for turbulence constant %E
static const scalar standardE = 9.8;
//- Runtime type information
TypeName("turbulenceModel");

20
src/turbulenceModels/incompressible/LES/derivedFvPatchFields/wallFunctions/nuSgsWallFunctions/nuSgsWallFunction/nuSgsWallFunctionFvPatchScalarField.C
View File

@ -25,6 +25,7 @@ License
\*---------------------------------------------------------------------------*/
#include "nuSgsWallFunctionFvPatchScalarField.H"
#include "LESModel.H"
#include "fvPatchFieldMapper.H"
#include "volFields.H"
#include "addToRunTimeSelectionTable.H"
@ -38,13 +39,6 @@ namespace incompressible
namespace LESModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
nuSgsWallFunctionFvPatchScalarField::
@ -57,8 +51,8 @@ nuSgsWallFunctionFvPatchScalarField
fixedValueFvPatchScalarField(p, iF),
UName_("U"),
nuName_("nu"),
kappa_(defaultKappa_),
E_(defaultE_)
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{}
@ -90,8 +84,8 @@ nuSgsWallFunctionFvPatchScalarField
fixedValueFvPatchScalarField(p, iF, dict),
UName_(dict.lookupOrDefault<word>("U", "U")),
nuName_(dict.lookupOrDefault<word>("nu", "nu")),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{}
@ -196,8 +190,8 @@ void nuSgsWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "U", "U", UName_);
writeEntryIfDifferent<word>(os, "nu", "nu", nuName_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

13
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/turbulentMixingLengthFrequencyInlet/turbulentMixingLengthFrequencyInletFvPatchScalarField.C
View File

@ -38,12 +38,6 @@ namespace Foam
namespace incompressible
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
//! @endcond fileScope
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
turbulentMixingLengthFrequencyInletFvPatchScalarField::
@ -121,8 +115,11 @@ void turbulentMixingLengthFrequencyInletFvPatchScalarField::updateCoeffs()
// Lookup Cmu corresponding to the turbulence model selected
const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
const scalar Cmu =
rasModel.coeffDict().lookupOrDefault<scalar>("Cmu", defaultCmu_);
const scalar Cmu = rasModel.coeffDict().lookupOrDefault<scalar>
(
"Cmu",
turbulenceModel::standardCmu
);
const scalar Cmu25 = pow(Cmu, 0.25);

26
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace incompressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void epsilonWallFunctionFvPatchScalarField::checkType()
@ -78,9 +70,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
GName_("RASModel::G"),
nuName_("nu"),
nutName_("nut"),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_)
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{
checkType();
}
@ -121,9 +113,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
GName_(dict.lookupOrDefault<word>("G", "RASModel::G")),
nuName_(dict.lookupOrDefault<word>("nu", "nu")),
nutName_(dict.lookupOrDefault<word>("nut", "nut")),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{
checkType();
}
@ -243,9 +235,9 @@ void epsilonWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "nu", "nu", nuName_);
writeEntryIfDifferent<word>(os, "nut", "nut", nutName_);
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

26
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutWallFunction/nutWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace incompressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
void nutWallFunctionFvPatchScalarField::checkType()
@ -114,9 +106,9 @@ tmp<scalarField> nutWallFunctionFvPatchScalarField::calcNut() const
void nutWallFunctionFvPatchScalarField::writeLocalEntries(Ostream& os) const
{
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
}
@ -129,9 +121,9 @@ nutWallFunctionFvPatchScalarField::nutWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_),
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE),
yPlusLam_(calcYPlusLam(kappa_, E_))
{
checkType();
@ -164,9 +156,9 @@ nutWallFunctionFvPatchScalarField::nutWallFunctionFvPatchScalarField
)
:
fixedValueFvPatchScalarField(p, iF, dict),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_)),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE)),
yPlusLam_(calcYPlusLam(kappa_, E_))
{
checkType();

26
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
View File

@ -40,14 +40,6 @@ namespace incompressible
namespace RASModels
{
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
//! @cond fileScope
static const scalar defaultCmu_(0.09);
static const scalar defaultKappa_(0.41);
static const scalar defaultE_(9.8);
//! @endcond fileScope
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void omegaWallFunctionFvPatchScalarField::checkType()
@ -78,9 +70,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
GName_("RASModel::G"),
nuName_("nu"),
nutName_("nut"),
Cmu_(defaultCmu_),
kappa_(defaultKappa_),
E_(defaultE_)
Cmu_(turbulenceModel::standardCmu),
kappa_(turbulenceModel::standardKappa),
E_(turbulenceModel::standardE)
{
checkType();
}
@ -121,9 +113,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
GName_(dict.lookupOrDefault<word>("G", "RASModel::G")),
nuName_(dict.lookupOrDefault<word>("nu", "nu")),
nutName_(dict.lookupOrDefault<word>("nut", "nut")),
Cmu_(dict.lookupOrDefault<scalar>("Cmu", defaultCmu_)),
kappa_(dict.lookupOrDefault<scalar>("kappa", defaultKappa_)),
E_(dict.lookupOrDefault<scalar>("E", defaultE_))
Cmu_(dict.lookupOrDefault<scalar>("Cmu", turbulenceModel::standardCmu)),
kappa_(dict.lookupOrDefault<scalar>("kappa", turbulenceModel::standardKappa)),
E_(dict.lookupOrDefault<scalar>("E", turbulenceModel::standardE))
{
checkType();
}
@ -233,9 +225,9 @@ void omegaWallFunctionFvPatchScalarField::write(Ostream& os) const
writeEntryIfDifferent<word>(os, "nu", "nu", nuName_);
writeEntryIfDifferent<word>(os, "nut", "nut", nutName_);
writeEntryIfDifferent<scalar>(os, "Cmu", defaultCmu_, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", defaultKappa_, kappa_);
writeEntryIfDifferent<scalar>(os, "E", defaultE_, E_);
writeEntryIfDifferent<scalar>(os, "Cmu", turbulenceModel::standardCmu, Cmu_);
writeEntryIfDifferent<scalar>(os, "kappa", turbulenceModel::standardKappa, kappa_);
writeEntryIfDifferent<scalar>(os, "E", turbulenceModel::standardE, E_);
writeEntry("value", os);
}

12
src/turbulenceModels/incompressible/turbulenceModel/turbulenceModel.H
View File

@ -96,6 +96,18 @@ private:
public:
// Public data
//- Standard value for turbulence constant %Cmu
static const scalar standardCmu = 0.09;
//- Standard value for turbulence constant %kappa
static const scalar standardKappa = 0.41;
//- Standard value for turbulence constant %E
static const scalar standardE = 9.8;
//- Runtime type information
TypeName("turbulenceModel");