40 template<
class BasicMomentumTransportModel>
43 this->nut_ = Cmu_*
sqr(k_)/epsilon_;
44 this->nut_.correctBoundaryConditions();
49 template<
class BasicMomentumTransportModel>
57 dimVolume*this->rho_.dimensions()*k_.dimensions()
64 template<
class BasicMomentumTransportModel>
73 dimVolume*this->rho_.dimensions()*epsilon_.dimensions()
82 template<
class BasicMomentumTransportModel>
182 this->groupName(
"k"),
183 this->runTime_.
name(),
194 this->groupName(
"epsilon"),
195 this->runTime_.
name(),
206 if (
type == typeName)
208 this->printCoeffs(
type);
215 template<
class BasicMomentumTransportModel>
220 Cmu_.readIfPresent(this->coeffDict());
221 C1_.readIfPresent(this->coeffDict());
222 C2_.readIfPresent(this->coeffDict());
223 C3_.readIfPresent(this->coeffDict());
224 sigmak_.readIfPresent(this->coeffDict());
225 sigmaEps_.readIfPresent(this->coeffDict());
226 eta0_.readIfPresent(this->coeffDict());
227 beta_.readIfPresent(this->coeffDict());
238 template<
class BasicMomentumTransportModel>
241 if (!this->turbulence_)
287 ((eta*(-eta/eta0_ + scalar(1)))/(beta_*eta3 + scalar(1)))
291 epsilon_.boundaryFieldRef().updateCoeffs();
307 epsEqn.
ref().relax();
309 epsEqn.
ref().boundaryManipulate(epsilon_.boundaryFieldRef());
312 bound(epsilon_, this->epsilonMin_);
334 bound(k_, this->kMin_);
Bound the given scalar field where it is below the specified minimum.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
Generic GeometricField class.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Templated abstract base class for RAS turbulence models.
virtual tmp< fvScalarMatrix > epsilonSource() const
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
virtual void correctNut()
virtual tmp< fvScalarMatrix > kSource() const
RNGkEpsilon(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const viscosity &viscosity, const word &type=typeName)
Construct from components.
virtual bool read()
Re-read model coefficients if they have changed.
static autoPtr< dictionary > New(Istream &)
Construct top-level dictionary on freestore from Istream.
Generic dimensioned Type class.
Eddy viscosity turbulence model base class.
Finite volume constraints.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
A special matrix type and solver, designed for finite volume solutions of scalar equations....
BasicMomentumTransportModel::alphaField alphaField
BasicMomentumTransportModel::rhoField rhoField
A class for managing temporary objects.
void clear() const
If object pointer points to valid object:
T & ref() const
Return non-const reference or generate a fatal error.
Abstract base class for all fluid physical properties.
A class for handling words, derived from string.
Foam::fvConstraints & fvConstraints(Foam::fvConstraints::New(mesh))
Foam::fvModels & fvModels(Foam::fvModels::New(mesh))
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
void correct(const RdeltaTType &rDeltaT, const RhoType &rho, volScalarField &psi, const surfaceScalarField &phiCorr, const SpType &Sp, const SuType &Su)
const dimensionedScalar G
Newtonian constant of gravitation.
tmp< VolField< typename outerProduct< vector, Type >::type > > grad(const SurfaceField< Type > &ssf)
tmp< VolField< Type > > div(const SurfaceField< Type > &ssf)
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const VolField< Type > &)
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const VolField< Type > &)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const VolField< Type > &vf, const word &name)
tmp< fvMatrix< Type > > ddt(const VolField< Type > &vf)
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
bool read(const char *, int32_t &)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
const dimensionSet dimTime
dimensionedScalar sqrt(const dimensionedScalar &ds)
const dimensionSet dimVolume
dimensioned< scalar > mag(const dimensioned< Type > &)
word typedName(Name name)
Return the name of the object within the given type.
bool bound(volScalarField &, const dimensionedScalar &min)
Bound the given scalar field where it is below the specified min value.
static scalar R(const scalar a, const scalar x)
word name(const complex &)
Return a string representation of a complex.
fileType type(const fileName &, const bool checkVariants=true, const bool followLink=true)
Return the file type: directory or file.
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.