38 template<
class BasicMomentumTransportModel>
49 template<
class BasicMomentumTransportModel>
59 chi3/(chi3 +
pow3(Cv1_))
64 template<
class BasicMomentumTransportModel>
75 1.0 - chi/(1.0 + chi*fv1)
80 template<
class BasicMomentumTransportModel>
95 template<
class BasicMomentumTransportModel>
111 + fv2(chi, fv1)*nuTilda_()/
sqr(kappa_*dTilda),
118 template<
class BasicMomentumTransportModel>
146 template<
class BasicMomentumTransportModel>
165 template<
class BasicMomentumTransportModel>
182 template<
class BasicMomentumTransportModel>
188 if (this->mesh_.cacheTemporaryObject(modelName(
"LESRegion")))
192 modelName(
"LESRegion"),
199 template<
class BasicMomentumTransportModel>
205 this->nut_ = nuTilda_*fv1;
211 template<
class BasicMomentumTransportModel>
214 correctNut(fv1(this->chi()));
220 template<
class BasicMomentumTransportModel>
279 Cw1_(Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_),
340 this->runTime_.timeName(),
350 if (type == typeName)
352 this->printCoeffs(type);
359 template<
class BasicMomentumTransportModel>
364 sigmaNut_.readIfPresent(this->coeffDict());
365 kappa_.readIfPresent(*
this);
367 Cb1_.readIfPresent(this->coeffDict());
368 Cb2_.readIfPresent(this->coeffDict());
369 Cw1_ = Cb1_/
sqr(kappa_) + (1.0 + Cb2_)/sigmaNut_;
371 Cw3_.readIfPresent(this->coeffDict());
372 Cv1_.readIfPresent(this->coeffDict());
373 Cs_.readIfPresent(this->coeffDict());
375 CDES_.readIfPresent(this->coeffDict());
376 ck_.readIfPresent(this->coeffDict());
387 template<
class BasicMomentumTransportModel>
394 (nuTilda_ + this->nu())/sigmaNut_
399 template<
class BasicMomentumTransportModel>
409 "dTildaExtrapolated",
410 this->mesh_.time().timeName(),
415 zeroGradientFvPatchScalarField::typeName
417 dTildaExtrapolated.
ref() = dTilda(chi, fv1,
fvc::grad(this->U_));
425 sqr(this->
nut()/ck_/dTildaExtrapolated)
430 volScalarField::Boundary& kBf = tk.
ref().boundaryFieldRef();
434 if (isA<wallFvPatch>(patches[
patchi]))
444 template<
class BasicMomentumTransportModel>
447 if (!this->turbulence_)
469 this->Stilda(chi, fv1, Omega, dTilda)
480 Cb1_*
alpha()*
rho()*Stilda*nuTilda_()
483 Cw1_*
alpha()*
rho()*fw(Stilda, dTilda)*nuTilda_()/
sqr(dTilda),
489 nuTildaEqn.
ref().relax();
494 nuTilda_.correctBoundaryConditions();
499 cacheLESRegion(dTilda);
Eddy viscosity LES SGS model base class.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
BasicMomentumTransportModel::rhoField rhoField
#define forAll(list, i)
Loop across all elements in list.
void clear() const
If object pointer points to valid object:
static tmp< DimensionedField< Type, GeoMesh > > New(const word &name, const Mesh &mesh, const dimensionSet &)
Return a temporary field constructed from name, mesh.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedTensor skew(const dimensionedTensor &dt)
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)
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Internal &, const PtrList< fvPatchField< scalar >> &)
Return a temporary field constructed from name,.
dimensionedScalar sqrt(const dimensionedScalar &ds)
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 > &)
DimensionedField< scalar, volMesh > Internal
Type of the internal field from which this GeometricField is derived.
dimensionedScalar neg(const dimensionedScalar &ds)
static const wallDist & New(const fvMesh &mesh)
virtual void correct()=0
Solve the turbulence equations and correct the turbulence viscosity.
A class for handling words, derived from string.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
virtual tmp< volScalarField::Internal > dTilda(const volScalarField::Internal &chi, const volScalarField::Internal &fv1, const volTensorField::Internal &gradU) const
Length scale.
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::Internal > r(const volScalarField::Internal &nur, const volScalarField::Internal &Omega, const volScalarField::Internal &dTilda) const
tmp< volScalarField > fv1(const volScalarField &chi) const
tmp< volScalarField > DnuTildaEff() const
Return the effective diffusivity for nuTilda.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
Internal & ref()
Return a reference to the dimensioned internal field.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
BasicMomentumTransportModel::alphaField alphaField
volScalarField & bound(volScalarField &, const dimensionedScalar &lowerBound)
Bound the given scalar field if it has gone unbounded.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
SpalartAllmarasDES(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &type=typeName)
Construct from components.
tmp< volScalarField::Internal > Omega(const volTensorField::Internal &gradU) const
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
BasicMomentumTransportModel::transportModel transportModel
virtual void cacheLESRegion(const volScalarField::Internal &dTilda) const
Cache the LES region indicator field.
void correctBoundaryConditions()
Correct boundary field.
dimensionedScalar pow6(const dimensionedScalar &ds)
tmp< volScalarField::Internal > Stilda(const volScalarField::Internal &chi, const volScalarField::Internal &fv1, const volScalarField::Internal &Omega, const volScalarField::Internal &dTilda) const
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
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 void correct()
Correct nuTilda and related properties.
tmp< volScalarField::Internal > fw(const volScalarField::Internal &Omega, const volScalarField::Internal &dTilda) const
tmp< volScalarField::Internal > fv2(const volScalarField::Internal &chi, const volScalarField::Internal &fv1) const
virtual void correctNut()
virtual tmp< volScalarField > k() const
Return SGS kinetic energy.