Go to the documentation of this file.
33 template<
class BasicTurbulenceModel>
38 Info<< coeffDict_.dictName() << coeffDict_ <<
endl;
45 template<
class BasicTurbulenceModel>
55 const word& propertiesName
69 LESDict_(this->subOrEmptyDict(
"LES")),
70 turbulence_(LESDict_.getOrDefault<
Switch>(
"turbulence",
true)),
71 printCoeffs_(LESDict_.getOrDefault<
Switch>(
"printCoeffs",
false)),
72 coeffDict_(LESDict_.optionalSubDict(
type +
"Coeffs")),
116 IOobject::groupName(
"delta", alphaRhoPhi.group()),
124 this->mesh_.deltaCoeffs();
130 template<
class BasicTurbulenceModel>
140 const word& propertiesName
147 IOobject::groupName(propertiesName, alphaRhoPhi.group()),
150 IOobject::MUST_READ_IF_MODIFIED,
161 dict.getCompat<
word>(
"model", {{
"LESModel", -2006}})
164 Info<<
"Selecting LES turbulence model " << modelType <<
endl;
166 auto* ctorPtr = dictionaryConstructorTable(modelType);
175 *dictionaryConstructorTablePtr_
179 return autoPtr<LESModel>
181 ctorPtr(
alpha,
rho,
U, alphaRhoPhi,
phi, transport, propertiesName)
188 template<
class BasicTurbulenceModel>
193 LESDict_ <<= this->subDict(
"LES");
194 LESDict_.readEntry(
"turbulence", turbulence_);
196 coeffDict_ <<= LESDict_.optionalSubDict(
type() +
"Coeffs");
198 delta_().read(LESDict_);
200 Ce_.readIfPresent(LESDict_);
202 kMin_.readIfPresent(LESDict_);
211 template<
class BasicTurbulenceModel>
219 IOobject::groupName(
"epsilon", this->alphaRhoPhi_.group()),
220 this->mesh_.time().timeName(),
223 this->Ce()*
pow(this->
k(), 1.5)/this->
delta()
228 template<
class BasicTurbulenceModel>
232 const scalar betaStar = 0.09;
239 IOobject::groupName(
"omega", this->alphaRhoPhi_.group()),
240 this->mesh_.time().timeName(),
243 this->
epsilon()/(betaStar*(this->
k() + k0))
248 template<
class BasicTurbulenceModel>
BasicTurbulenceModel::rhoField rhoField
IOdictionary is derived from dictionary and IOobject to give the dictionary automatic IO functionalit...
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
A simple wrapper around bool so that it can be read as a word: true/false, on/off,...
A class for handling words, derived from Foam::string.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
A class for managing temporary objects.
const dimensionSet dimVelocity
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
static autoPtr< LESModel > New(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName)
Return a reference to the selected LES model.
Ostream & endl(Ostream &os)
Add newline and flush stream.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
BasicTurbulenceModel::transportModel transportModel
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
#define FatalIOErrorInLookup(ios, lookupTag, lookupName, lookupTable)
Report an error message using Foam::FatalIOError.
virtual bool read()
Read model coefficients if they have changed.
messageStream Info
Information stream (stdout output on master, null elsewhere)
virtual void printCoeffs(const word &type)
Print model coefficients.
const dictionary & subDict(const word &keyword, enum keyType::option matchOpt=keyType::REGEX) const
Find and return a sub-dictionary.
LESModel(const LESModel &)=delete
No copy construct.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Generic dimensioned Type class.
errorManipArg< error, int > exit(error &err, const int errNo=1)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
label k
Boltzmann constant.
BasicTurbulenceModel::alphaField alphaField
const dimensionSet dimless
Dimensionless.