39 template<
class BasicTurbulenceModel>
42 return max(k_/epsilon_, 6.0*
sqrt(this->
nu()/epsilon_));
46 template<
class BasicTurbulenceModel>
55 template<
class BasicTurbulenceModel>
58 this->nut_ =
min(CmuKEps_*
sqr(k_)/epsilon_, this->Cmu_*v2_*Ts());
59 this->nut_.correctBoundaryConditions();
62 BasicTurbulenceModel::correctNut();
68 template<
class BasicTurbulenceModel>
77 const word& propertiesName,
190 this->runTime_.timeName(),
202 this->runTime_.timeName(),
214 this->runTime_.timeName(),
226 this->runTime_.timeName(),
236 bound(k_, this->kMin_);
237 bound(epsilon_, this->epsilonMin_);
241 if (type == typeName)
243 this->printCoeffs(type);
250 template<
class BasicTurbulenceModel>
255 Cmu_.readIfPresent(this->coeffDict());
256 CmuKEps_.readIfPresent(this->coeffDict());
257 C1_.readIfPresent(this->coeffDict());
258 C2_.readIfPresent(this->coeffDict());
259 CL_.readIfPresent(this->coeffDict());
260 Ceta_.readIfPresent(this->coeffDict());
261 Ceps2_.readIfPresent(this->coeffDict());
262 Ceps3_.readIfPresent(this->coeffDict());
263 sigmaK_.readIfPresent(this->coeffDict());
264 sigmaEps_.readIfPresent(this->coeffDict());
275 template<
class BasicTurbulenceModel>
278 if (!this->turbulence_)
307 1.0/Ts*((C1_ - N)*v2_ - 2.0/3.0*k_*(C1_ - 1.0))
313 1.4*(1.0 + 0.05*
min(
sqrt(k_/v2_), scalar(100.0)))
317 epsilon_.boundaryFieldRef().updateCoeffs();
327 -
fvm::SuSp(((2.0/3.0)*Ceps1 + Ceps3_)*alpha*rho*divU, epsilon_)
328 -
fvm::Sp(Ceps2_*alpha*rho/Ts, epsilon_)
332 epsEqn.
ref().relax();
334 epsEqn.
ref().boundaryManipulate(epsilon_.boundaryFieldRef());
337 bound(epsilon_, this->epsilonMin_);
348 -
fvm::SuSp((2.0/3.0)*alpha*rho*divU, k_)
349 -
fvm::Sp(alpha*rho*epsilon_/k_, k_)
357 bound(k_, this->kMin_);
366 - 1.0/L2/k_*(v2fAlpha - C2_*G)
383 alpha*rho*
min(k_*f_, C2_*G - v2fAlpha)
384 -
fvm::Sp(N*alpha*rho*epsilon_/k_, v2_)
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
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.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
BasicTurbulenceModel::rhoField rhoField
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensionedScalar pow025(const dimensionedScalar &ds)
Generic dimensioned Type class.
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
Eddy viscosity turbulence model base class.
BasicTurbulenceModel::transportModel transportModel
tmp< volScalarField > Ls() const
Return length scale, Ls.
Templated abstract base class for RAS turbulence models.
Abstract base-class for v2-f models to provide BCs access to the v2 and f fields. ...
bool read(const char *, int32_t &)
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
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.
v2f(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.
virtual void correctNut()
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)
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()
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
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)
virtual bool read()
Read RASProperties dictionary.
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.
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...
tmp< volScalarField > Ts() const
Return time scale, Ts.
BasicTurbulenceModel::alphaField alphaField