42 template<
class BasicTurbulenceModel>
45 this->nut_ = k_/omega_;
46 this->nut_.correctBoundaryConditions();
49 BasicTurbulenceModel::correctNut();
55 template<
class BasicTurbulenceModel>
58 const alphaField&
alpha,
64 const word& propertiesName,
154 if (type == typeName)
163 template<
class BasicTurbulenceModel>
168 Cmu_.readIfPresent(this->coeffDict());
169 beta_.readIfPresent(this->coeffDict());
170 gamma_.readIfPresent(this->coeffDict());
171 alphaK_.readIfPresent(this->coeffDict());
172 alphaOmega_.readIfPresent(this->coeffDict());
181 template<
class BasicTurbulenceModel>
184 if (!this->turbulence_)
190 const alphaField&
alpha = this->alpha_;
191 const rhoField&
rho = this->rho_;
197 eddyViscosity<RASModel<BasicTurbulenceModel>>
::correct();
213 omega_.boundaryFieldRef().updateCoeffs();
216 tmp<fvScalarMatrix> omegaEqn
228 omegaEqn.ref().relax();
230 omegaEqn.ref().boundaryManipulate(omega_.boundaryFieldRef());
233 bound(omega_, this->omegaMin_);
237 tmp<fvScalarMatrix> kEqn
253 bound(k_, this->kMin_);
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
SolverPerformance< Type > solve(faMatrix< Type > &, const dictionary &)
Solve returning the solution statistics given convergence tolerance.
dimensionedScalar kMin_
Lower limit of k.
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.
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Generic dimensioned Type class.
Eddy viscosity turbulence model base class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
DimensionedField< scalar, volMesh > Internal
The internal field type from which this GeometricField is derived.
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
bool read(const char *buf, int32_t &val)
Same as readInt32.
Templated abstract base class for RAS turbulence models.
constexpr const char *const group
Group name for atomic constants.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
A class for handling words, derived from Foam::string.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
zeroField SuSp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
virtual void printCoeffs(const word &type)
Print model coefficients.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
virtual void correct()
Solve the turbulence equations and correct the turbulence viscosity.
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(psi *p+alphal *rhol0+((alphav *psiv+alphal *psil) - psi) *pSat, rhoMin);# 1 "/home/chef2/andy/OpenFOAM/release/v2212/OpenFOAM-v2212/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()
dimensionedScalar omegaMin_
Lower limit for omega.
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)
Base-class for all transport models used by the incompressible turbulence models. ...
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
static options & New(const fvMesh &mesh)
Construct fvOptions and register to database if not present.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
Defines the attributes of an object for which implicit objectRegistry management is supported...