38 template<
class BasicTurbulenceModel>
41 return nuTilda_/this->
nu();
45 template<
class BasicTurbulenceModel>
52 return chi3/(chi3 +
pow3(Cv1_));
56 template<
class BasicTurbulenceModel>
63 return 1.0 - chi/(1.0 + chi*fv1);
67 template<
class BasicTurbulenceModel>
77 template<
class BasicTurbulenceModel>
87 template<
class BasicTurbulenceModel>
101 + fv2(chi, fv1)*nuTilda_/
sqr(kappa_*dTilda),
108 template<
class BasicTurbulenceModel>
132 tr.
ref().boundaryFieldRef() == 0.0;
138 template<
class BasicTurbulenceModel>
152 template<
class BasicTurbulenceModel>
161 min(tdTilda.
ref().ref(), tdTilda(), y_);
166 template<
class BasicTurbulenceModel>
172 this->nut_ = nuTilda_*fv1;
176 BasicTurbulenceModel::correctNut();
180 template<
class BasicTurbulenceModel>
183 correctNut(fv1(this->chi()));
189 template<
class BasicTurbulenceModel>
198 const word& propertiesName,
250 Cw1_(Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_),
311 this->runTime_.timeName(),
321 if (type == typeName)
323 this->printCoeffs(type);
330 template<
class BasicTurbulenceModel>
335 sigmaNut_.readIfPresent(this->coeffDict());
336 kappa_.readIfPresent(*
this);
338 Cb1_.readIfPresent(this->coeffDict());
339 Cb2_.readIfPresent(this->coeffDict());
340 Cw1_ = Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_;
342 Cw3_.readIfPresent(this->coeffDict());
343 Cv1_.readIfPresent(this->coeffDict());
344 Cs_.readIfPresent(this->coeffDict());
346 CDES_.readIfPresent(this->coeffDict());
347 ck_.readIfPresent(this->coeffDict());
358 template<
class BasicTurbulenceModel>
365 (nuTilda_ + this->
nu())/sigmaNut_
370 template<
class BasicTurbulenceModel>
379 template<
class BasicTurbulenceModel>
398 template<
class BasicTurbulenceModel>
401 if (!this->turbulence_)
421 const volScalarField Stilda(this->Stilda(chi, fv1, Omega, dTilda));
430 Cb1_*alpha*rho*Stilda*nuTilda_
433 Cw1_*alpha*rho*fw(Stilda, dTilda)*nuTilda_/
sqr(dTilda),
439 nuTildaEqn.
ref().relax();
444 nuTilda_.correctBoundaryConditions();
Eddy viscosity LES SGS model base class.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
SpalartAllmarasDES(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName, const word &type=typeName)
Construct from components.
BasicTurbulenceModel::transportModel transportModel
tmp< volScalarField > Stilda(const volScalarField &chi, const volScalarField &fv1, const volScalarField &Omega, const volScalarField &dTilda) const
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedTensor skew(const dimensionedTensor &dt)
tmp< volScalarField > fv2(const volScalarField &chi, const volScalarField &fv1) const
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
T & ref() const
Return non-const reference or generate a fatal error.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
tmp< volScalarField > r(const volScalarField &nur, const volScalarField &Omega, const volScalarField &dTilda) const
bool readIfPresent(const dictionary &)
Update the value of dimensioned<Type> if found in the dictionary.
virtual bool read()
Read model coefficients if they have changed.
Generic dimensioned Type class.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensionedScalar neg(const dimensionedScalar &ds)
BasicTurbulenceModel::alphaField alphaField
static const wallDist & New(const fvMesh &mesh)
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Mesh &, const dimensionSet &, const word &patchFieldType=fvPatchField< scalar >::calculatedType())
Return a temporary field constructed from name, mesh, dimensionSet.
virtual tmp< volScalarField > LESRegion() const
Return the LES field indicator.
virtual void correct()=0
Solve the turbulence equations and correct the turbulence viscosity.
const dimensionSet & dimensions() const
Return dimensions.
A class for handling words, derived from string.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
void constrain(fvMatrix< Type > &eqn)
Apply constraints to equation.
dimensioned< scalar > magSqr(const dimensioned< Type > &)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< volScalarField > fv1(const volScalarField &chi) const
tmp< volScalarField > DnuTildaEff() const
Return the effective diffusivity for nuTilda.
tmp< volScalarField > S(const volTensorField &gradU) const
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
volScalarField & bound(volScalarField &, const dimensionedScalar &lowerBound)
Bound the given scalar field if it has gone unbounded.
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
void correctBoundaryConditions()
Correct boundary field.
dimensionedScalar pow6(const dimensionedScalar &ds)
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
BasicTurbulenceModel::rhoField rhoField
dimensioned< scalar > mag(const dimensioned< Type > &)
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
A class for managing temporary objects.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const dimensionedVector & g
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
tmp< volScalarField > chi() const
virtual tmp< volScalarField > dTilda(const volScalarField &chi, const volScalarField &fv1, const volTensorField &gradU) const
Length scale.
virtual void correct()
Correct nuTilda and related properties.
virtual void correctNut()
tmp< volScalarField > fw(const volScalarField &Omega, const volScalarField &dTilda) const
tmp< volScalarField > Omega(const volTensorField &gradU) const
virtual tmp< volScalarField > k() const
Return SGS kinetic energy.