40 template<
class BasicTurbulenceModel>
41 const IDDESDelta& SpalartAllmarasIDDES<BasicTurbulenceModel>::setDelta()
const 43 if (!isA<IDDESDelta>(this->delta_()))
46 <<
"The delta function must be set to a " << IDDESDelta::typeName
50 return refCast<const IDDESDelta>(this->delta_());
54 template<
class BasicTurbulenceModel>
55 tmp<volScalarField> SpalartAllmarasIDDES<BasicTurbulenceModel>::alpha()
const 57 return max(0.25 - this->y_/IDDESDelta_.hmax(), scalar(-5));
61 template<
class BasicTurbulenceModel>
62 tmp<volScalarField> SpalartAllmarasIDDES<BasicTurbulenceModel>::ft
67 return tanh(
pow3(
sqr(Ct_)*this->r(this->nut_, magGradU, this->y_)));
71 template<
class BasicTurbulenceModel>
72 tmp<volScalarField> SpalartAllmarasIDDES<BasicTurbulenceModel>::fl
77 return tanh(
pow(
sqr(Cl_)*this->r(this->
nu(), magGradU, this->y_), 10));
81 template<
class BasicTurbulenceModel>
82 tmp<volScalarField> SpalartAllmarasIDDES<BasicTurbulenceModel>::fdt
87 return 1 -
tanh(
pow(Cdt1_*this->r(this->nut_, magGradU, this->y_), Cdt2_));
93 template<
class BasicTurbulenceModel>
123 fdTilda*(1 + fe)*lRAS + (1 - fdTilda)*lLES,
132 lerp(lLES, lRAS, fdTilda),
140 template<
class BasicTurbulenceModel>
143 const alphaField&
alpha,
149 const word& propertiesName,
211 IDDESDelta_(setDelta())
213 if (
type == typeName)
215 this->printCoeffs(
type);
222 template<
class BasicTurbulenceModel>
227 Cdt1_.readIfPresent(this->coeffDict());
228 Cdt2_.readIfPresent(this->coeffDict());
229 Cl_.readIfPresent(this->coeffDict());
230 Ct_.readIfPresent(this->coeffDict());
239 template<
class BasicTurbulenceModel>
dimensionedScalar tanh(const dimensionedScalar &ds)
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
errorManipArg< error, int > exit(error &err, const int errNo=1)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
error FatalError
Error stream (stdout output on all processes), with additional 'FOAM FATAL ERROR' header text and sta...
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual tmp< volScalarField > dTilda(const volScalarField &chi, const volScalarField &fv1, const volTensorField &gradU) const
Return the length scale.
A simple wrapper around bool so that it can be read as a word: true/false, on/off, yes/no, any/none. Also accepts 0/1 as a string and shortcuts t/f, y/n.
virtual bool read()
Re-read model coefficients if they have changed.
SpalartAllmarasDES DES turbulence model for incompressible and compressible flows.
Generic dimensioned Type class.
virtual bool read()
Re-read model coefficients if they have changed.
SpalartAllmaras IDDES turbulence model for incompressible and compressible flows. ...
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
dimensionedScalar exp(const dimensionedScalar &ds)
A class for handling words, derived from Foam::string.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedScalar pos0(const dimensionedScalar &ds)
virtual tmp< volScalarField > fd() const
Return the shielding function.
dimensioned< Type > lerp(const dimensioned< Type > &a, const dimensioned< Type > &b, const scalar t)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Base-class for all transport models used by the incompressible turbulence models. ...
const volScalarField & psi
A class for managing temporary objects.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].