Go to the documentation of this file.
41 template<
class BasicTurbulenceModel>
45 this->nut_ = Cmu_*phit_*k_*T_;
46 this->nut_.correctBoundaryConditions();
49 BasicTurbulenceModel::correctNut();
53 template<
class BasicTurbulenceModel>
73 template<
class BasicTurbulenceModel>
80 pow(k_, 1.5)/epsilon_,
98 template<
class BasicTurbulenceModel>
107 const word& propertiesName,
255 this->runTime_.timeName(),
267 this->runTime_.timeName(),
279 this->runTime_.timeName(),
291 this->runTime_.timeName(),
303 this->runTime_.timeName(),
318 bound(k_, this->kMin_);
319 bound(epsilon_, this->epsilonMin_);
320 bound(phit_, phitMin_);
323 if (
type == typeName)
325 this->printCoeffs(
type);
330 mag(sigmaK_.value()) < VSMALL
331 ||
mag(sigmaEps_.value()) < VSMALL
332 ||
mag(sigmaPhit_.value()) < VSMALL
336 <<
"Non-zero values are required for the model constants:" <<
nl
337 <<
"sigmaK = " << sigmaK_ <<
nl
338 <<
"sigmaEps = " << sigmaEps_ <<
nl
339 <<
"sigmaPhit = " << sigmaPhit_ <<
nl
347 template<
class BasicTurbulenceModel>
352 includeNu_.readIfPresent(
"includeNu", this->coeffDict());
353 Cmu_.readIfPresent(this->coeffDict());
354 Ceps1a_.readIfPresent(this->coeffDict());
355 Ceps1b_.readIfPresent(this->coeffDict());
356 Ceps1c_.readIfPresent(this->coeffDict());
357 Ceps2_.readIfPresent(this->coeffDict());
358 Cf1_.readIfPresent(this->coeffDict());
359 Cf2_.readIfPresent(this->coeffDict());
360 CL_.readIfPresent(this->coeffDict());
361 Ceta_.readIfPresent(this->coeffDict());
362 CT_.readIfPresent(this->coeffDict());
363 sigmaK_.readIfPresent(this->coeffDict());
364 sigmaEps_.readIfPresent(this->coeffDict());
365 sigmaPhit_.readIfPresent(this->coeffDict());
374 template<
class BasicTurbulenceModel>
377 if (!this->turbulence_)
409 Ceps1a_*(Ceps1b_ + Ceps1c_*
sqrt(1.0/phit_()))
414 epsilon_.boundaryFieldRef().updateCoeffs();
434 epsEqn.
ref().relax();
436 epsEqn.
ref().boundaryManipulate(epsilon_.boundaryFieldRef());
439 bound(epsilon_, this->epsilonMin_);
460 bound(k_, this->kMin_);
471 (Cf1_ - 1.0)*(phit_() - 2.0/3.0)/T_()
473 +(Cf2_*(2.0/3.0)*
divU)
500 /(k_()*sigmaPhit_*phit_())
507 phitEqn.
ref().relax();
511 bound(phit_, phitMin_);
Defines the attributes of an object for which implicit objectRegistry management is supported,...
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
virtual void correct()
Solve the transport equations and correct the turbulent viscosity.
A class for handling words, derived from Foam::string.
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
void clear() const noexcept
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
const dimensionedScalar G
Newtonian constant of gravitation.
static constexpr const zero Zero
Global zero (0)
Templated abstract base class for RAS turbulence models.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
BasicTurbulenceModel::alphaField alphaField
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
volScalarField & bound(volScalarField &, const dimensionedScalar &lowerBound)
Bound the given scalar field if it has gone unbounded.
bool read(const char *buf, int32_t &val)
Same as readInt32.
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
BasicTurbulenceModel::rhoField rhoField
dimensionedScalar pow025(const dimensionedScalar &ds)
tmp< volScalarField > Ls() const
Return the turbulent length scale, L.
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
BasicTurbulenceModel::transportModel transportModel
virtual void correctNut()
Update nut with the latest available k, phit, and T.
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensionedScalar pow3(const dimensionedScalar &ds)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Bound the given scalar field if it has gone unbounded.
SolverPerformance< Type > solve(faMatrix< Type > &, Istream &)
Solve returning the solution statistics given convergence tolerance.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
static dimensioned< Type > getOrAddToDict(const word &name, dictionary &dict, const dimensionSet &dims=dimless, const Type &deflt=Type(Zero))
Construct dimensioned from dictionary, with default value.
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
Eddy viscosity turbulence model base class.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
The k-epsilon-phit-f turbulence closure model for incompressible and compressible flows.
virtual bool read()
Re-read model coefficients if they have changed.
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
tmp< volScalarField > Ts() const
Return the turbulent time scale, T.
static Switch getOrAddToDict(const word &key, dictionary &dict, const Switch deflt=switchType::FALSE)
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)