39 template<
class BasicTurbulenceModel>
42 this->nut_ = k_/omega_;
43 this->nut_.correctBoundaryConditions();
46 BasicTurbulenceModel::correctNut();
50 template<
class BasicTurbulenceModel>
58 dimVolume*this->rho_.dimensions()*k_.dimensions()
65 template<
class BasicTurbulenceModel>
81 template<
class BasicTurbulenceModel>
90 const word& propertiesName,
157 this->runTime_.timeName(),
169 this->runTime_.timeName(),
177 bound(k_, this->kMin_);
178 bound(omega_, this->omegaMin_);
180 if (type == typeName)
182 this->printCoeffs(type);
189 template<
class BasicTurbulenceModel>
194 Cmu_.readIfPresent(this->coeffDict());
195 beta_.readIfPresent(this->coeffDict());
196 gamma_.readIfPresent(this->coeffDict());
197 alphaK_.readIfPresent(this->coeffDict());
198 alphaOmega_.readIfPresent(this->coeffDict());
209 template<
class BasicTurbulenceModel>
212 if (!this->turbulence_)
236 nut.
v()*(
dev(
twoSymm(tgradU().v())) && tgradU().v())
241 omega_.boundaryFieldRef().updateCoeffs();
250 gamma_*
alpha()*
rho()*G*omega_()/k_()
257 omegaEqn.
ref().relax();
259 omegaEqn.
ref().boundaryManipulate(omega_.boundaryFieldRef());
262 bound(omega_, this->omegaMin_);
283 bound(k_, this->kMin_);
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:
BasicTurbulenceModel::transportModel transportModel
kOmega(const alphaField &alpha, const rhoField &rho, const volVectorField &U, const surfaceScalarField &alphaRhoPhi, const surfaceScalarField &phi, const transportModel &transport, const word &propertiesName=turbulenceModel::propertiesName, const word &type=typeName)
Construct from components.
void correct(GeometricField< Type, fvPatchField, volMesh > &field)
Apply correction to field.
const dimensionedScalar G
Newtonian constant of gravitation.
T & ref() const
Return non-const reference or generate a fatal error.
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Generic dimensioned Type class.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
BasicTurbulenceModel::alphaField alphaField
Eddy viscosity turbulence model base class.
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
Templated abstract base class for RAS turbulence models.
const dimensionSet dimVolume(pow3(dimLength))
bool read(const char *, int32_t &)
const Internal & v() const
Return a const-reference to the dimensioned internal field.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
BasicTurbulenceModel::rhoField rhoField
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 correct()
Solve the turbulence equations and correct the turbulence viscosity.
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)
Bound the given scalar field if it has gone unbounded.
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(rho0+psi *p, rhoMin);# 1 "/home/ubuntu/OpenFOAM-7/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()
virtual void correctNut()
virtual bool read()
Read RASProperties dictionary.
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.
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.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
virtual tmp< fvScalarMatrix > kSource() const
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.
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.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
virtual tmp< fvScalarMatrix > omegaSource() const