38 template<
class BasicTurbulenceModel>
41 this->nut_ = Ck_*
sqrt(this->
k())*this->
delta();
42 this->nut_.correctBoundaryConditions();
45 BasicTurbulenceModel::correctNut();
51 template<
class BasicTurbulenceModel>
60 const word& propertiesName,
113 if (type == typeName)
115 this->printCoeffs(type);
116 this->boundNormalStress(this->R_);
123 template<
class BasicTurbulenceModel>
128 Ck_.readIfPresent(this->coeffDict());
129 Cm_.readIfPresent(this->coeffDict());
130 Ce_.readIfPresent(this->coeffDict());
131 Cs_.readIfPresent(this->coeffDict());
142 template<
class BasicTurbulenceModel>
154 this->runTime_.timeName(),
165 template<
class BasicTurbulenceModel>
168 if (!this->turbulence_)
200 + (4.0/5.0)*alpha*rho*k*D
201 - ((2.0/3.0)*(1.0 - Cm_/this->Ce_)*
I)*(alpha*rho*this->
epsilon())
209 this->boundNormalStress(R);
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Templated abstract base class for LES SGS models.
Differential SGS Stress Equation Model for incompressible and compressible flows. ...
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
T & ref() const
Return non-const reference or generate a fatal error.
BasicTurbulenceModel::rhoField rhoField
dimensionedScalar sqrt(const dimensionedScalar &ds)
BasicTurbulenceModel::alphaField alphaField
label k
Boltzmann constant.
Generic dimensioned Type class.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
BasicTurbulenceModel::transportModel transportModel
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual void correct()
Correct sub-grid stress, eddy-Viscosity and related properties.
bool read(const char *, int32_t &)
static const Identity< scalar > I
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)
void constrain(fvMatrix< Type > &eqn)
Apply constraints to equation.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(psi *p+alphal *rhol0+((alphav *psiv+alphal *psil) - psi) *pSat, rhoMin);# 1 "/home/ubuntu/OpenFOAM-5.0/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H" 1{ alphav=max(min((rho - rholSat)/(rhovSat - rholSat), scalar(1)), scalar(0));alphal=1.0 - alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
fileName::Type type(const fileName &, const bool followLink=true)
Return the file type: DIRECTORY or FILE.
#define R(A, B, C, D, E, F, K, M)
dimensionedSymmTensor symm(const dimensionedSymmTensor &dt)
tmp< fvMatrix< Type > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
virtual void correctNut()
Update the eddy-viscosity.
virtual bool read()
Read model coefficients if they have changed.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
A class for managing temporary objects.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Reynolds-stress turbulence model base class.