38 template<
class BasicMomentumTransportModel>
44 0.5*(filter_(
magSqr(this->U_)) -
magSqr(filter_(this->U_))),
50 template<
class BasicMomentumTransportModel>
59 simpleFilter_(
dev(filter_(
sqr(this->U_)) - (
sqr(filter_(this->U_)))))
64 simpleFilter_(-2.0*this->
delta()*
sqrt(KK)*filter_(D))
69 simpleFilter_(0.5*(LL && MM))
81 template<
class BasicMomentumTransportModel>
90 simpleFilter_(this->nuEff()*(filter_(
magSqr(D)) -
magSqr(filter_(D))))
91 /simpleFilter_(
pow(KK, 1.5)/(2.0*this->
delta()))
99 template<
class BasicMomentumTransportModel>
107 template<
class BasicMomentumTransportModel>
114 this->nut_ = Ck(D, KK)*
sqrt(k_)*this->
delta();
115 this->nut_.correctBoundaryConditions();
120 template<
class BasicMomentumTransportModel>
127 template<
class BasicMomentumTransportModel>
135 dimVolume*this->rho_.dimensions()*k_.dimensions()
144 template<
class BasicMomentumTransportModel>
172 this->runTime_.timeName(),
180 simpleFilter_(this->mesh_),
182 filter_(filterPtr_())
184 bound(k_, this->kMin_);
186 if (type == typeName)
188 this->printCoeffs(type);
195 template<
class BasicMomentumTransportModel>
200 filter_.read(this->coeffDict());
211 template<
class BasicMomentumTransportModel>
222 template<
class BasicMomentumTransportModel>
225 if (!this->turbulence_)
256 -
fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_)
266 bound(k_, this->kMin_);
Eddy viscosity LES SGS model base class.
virtual bool read()
Read model coefficients if they have changed.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
fvMatrix< scalar > fvScalarMatrix
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
void clear() const
If object pointer points to valid object:
virtual void correct()
Correct Eddy-Viscosity and related properties.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
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)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
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)
label k
Boltzmann constant.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
One equation eddy-viscosity model.
BasicMomentumTransportModel::transportModel transportModel
virtual void correct()=0
Solve the turbulence equations and correct the turbulence viscosity.
const dimensionSet dimVolume(pow3(dimLength))
const dimensionSet & dimensions() const
Return dimensions.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
volScalarField Ck(const volSymmTensorField &D, const volScalarField &KK) const
Calculate Ck by filtering the velocity field U.
A class for handling words, derived from string.
static word groupName(Name name, const word &group)
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
virtual void correctNut()
volScalarField Ce() const
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)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
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.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
virtual tmp< volScalarField > epsilon() const
Return sub-grid disipation rate.
virtual tmp< fvScalarMatrix > kSource() const
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< volScalarField > KK() const
Return the KK field obtained by filtering the velocity field U.
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.
BasicMomentumTransportModel::alphaField alphaField
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
static autoPtr< LESfilter > New(const fvMesh &, const dictionary &, const word &filterDictName="filter")
Return a reference to the selected LES filter.
dynamicKEqn(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.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
const dimensionedScalar & G
Newtonian constant of gravitation.
BasicMomentumTransportModel::rhoField rhoField
const dimensionSet dimVelocity