42template<
class BasicTurbulenceModel>
45 this->nut_ = this->Cmu_*
sqr(k_)/epsilon_;
46 this->nut_.correctBoundaryConditions();
49 BasicTurbulenceModel::correctNut();
55template<
class BasicTurbulenceModel>
64 const word& propertiesName,
215 if (
type == typeName)
228template<
class BasicTurbulenceModel>
233 Cmu_.readIfPresent(this->coeffDict());
234 C1_.readIfPresent(this->coeffDict());
235 C1s_.readIfPresent(this->coeffDict());
236 C2_.readIfPresent(this->coeffDict());
237 C3_.readIfPresent(this->coeffDict());
238 C3s_.readIfPresent(this->coeffDict());
239 C4_.readIfPresent(this->coeffDict());
240 C5_.readIfPresent(this->coeffDict());
242 Ceps1_.readIfPresent(this->coeffDict());
243 Ceps2_.readIfPresent(this->coeffDict());
244 Cs_.readIfPresent(this->coeffDict());
245 Ceps_.readIfPresent(this->coeffDict());
254template<
class BasicTurbulenceModel>
262 (Cs_*(this->k_/this->epsilon_))*this->R_ +
I*this->
nu()
268template<
class BasicTurbulenceModel>
276 (Ceps_*(this->k_/this->epsilon_))*this->R_ +
I*this->
nu()
282template<
class BasicTurbulenceModel>
285 if (!this->turbulence_)
307 epsilon_.boundaryFieldRef().updateCoeffs();
321 epsEqn.
ref().relax();
323 epsEqn.
ref().boundaryManipulate(epsilon_.boundaryFieldRef());
326 bound(epsilon_, this->epsilonMin_);
337 if (isA<wallFvPatch>(curPatch))
341 label celli = curPatch.
faceCells()[facei];
344 G[celli]/(0.5*
mag(
tr(P[celli])) + SMALL),
364 - ((1.0/3.0)*
I)*(((2.0 - C1_)*epsilon_ - C1s_*G)*
alpha*
rho)
380 this->boundNormalStress(
R);
383 bound(k_, this->kMin_);
388 this->correctWallShearStress(
R);
#define R(A, B, C, D, E, F, K, M)
Defines the attributes of an object for which implicit objectRegistry management is supported,...
A list of pointers to objects of type <T>, with allocation/deallocation management of the pointers....
Templated abstract base class for RAS turbulence models.
virtual void printCoeffs(const word &type)
Print model coefficients.
dimensionedScalar epsilonMin_
Lower limit of epsilon.
Speziale, Sarkar and Gatski Reynolds-stress turbulence model for incompressible and compressible flow...
BasicTurbulenceModel::alphaField alphaField
tmp< volSymmTensorField > DepsilonEff() const
Return the effective diffusivity for epsilon.
BasicTurbulenceModel::rhoField rhoField
virtual void correct()
Solve the turbulence equations and correct eddy-Viscosity and.
tmp< volSymmTensorField > DREff() const
Return the effective diffusivity for R.
virtual void correctNut()
Update the eddy-viscosity.
BasicTurbulenceModel::transportModel transportModel
virtual bool read()
Read model coefficients if they have changed.
Reynolds-stress turbulence model base class.
void boundNormalStress(volSymmTensorField &R) const
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
Generic dimensioned Type class.
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
virtual const labelUList & faceCells() const
Return faceCells.
A class for managing temporary objects.
A class for handling words, derived from Foam::string.
const polyBoundaryMesh & patches
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh > > grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
bool read(const char *buf, int32_t &val)
Same as readInt32.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
volScalarField & bound(volScalarField &, const dimensionedScalar &lowerBound)
Bound the given scalar field if it has gone unbounded.
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
static const Identity< scalar > I
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor innerSqr(const dimensionedSymmTensor &dt)
dimensionedTensor skew(const dimensionedTensor &dt)
#define forAll(list, i)
Loop across all elements in list.