42 template<
class BasicTurbulenceModel>
45 this->nut_ = Cmu_*
sqr(k_)/epsilon_;
46 this->nut_.correctBoundaryConditions();
49 BasicTurbulenceModel::correctNut();
53 template<
class BasicTurbulenceModel>
61 dimVolume*this->rho_.dimensions()*k_.dimensions()
68 template<
class BasicTurbulenceModel>
76 dimVolume*this->rho_.dimensions()*epsilon_.dimensions()
85 template<
class BasicTurbulenceModel>
88 const alphaField&
alpha,
94 const word& propertiesName,
190 bound(k_, this->kMin_);
191 bound(epsilon_, this->epsilonMin_);
193 if (
type == typeName)
195 this->printCoeffs(
type);
202 template<
class BasicTurbulenceModel>
207 Cmu_.readIfPresent(this->coeffDict());
208 C1_.readIfPresent(this->coeffDict());
209 C2_.readIfPresent(this->coeffDict());
210 C3_.readIfPresent(this->coeffDict());
211 sigmak_.readIfPresent(this->coeffDict());
212 sigmaEps_.readIfPresent(this->coeffDict());
221 template<
class BasicTurbulenceModel>
224 if (!this->turbulence_)
230 const alphaField&
alpha = this->alpha_;
231 const rhoField&
rho = this->rho_;
238 eddyViscosity<RASModel<BasicTurbulenceModel>>
::correct();
240 const volScalarField::Internal
divU 246 const volScalarField::Internal GbyNu
248 IOobject::scopedName(this->
type(),
"GbyNu"),
251 const volScalarField::Internal
G(this->GName(),
nut()*GbyNu);
255 epsilon_.boundaryFieldRef().updateCoeffs();
257 epsilon_.boundaryFieldRef().template evaluateCoupled<coupledFvPatch>();
260 tmp<fvScalarMatrix> epsEqn
273 epsEqn.ref().relax();
275 epsEqn.ref().boundaryManipulate(epsilon_.boundaryFieldRef());
278 bound(epsilon_, this->epsilonMin_);
281 tmp<fvScalarMatrix> kEqn
298 bound(k_, this->kMin_);
virtual tmp< fvScalarMatrix > epsilonSource() const
void correct(GeometricField< Type, PatchField, GeoMesh > &field)
Apply correction to field.
fvMatrix< scalar > fvScalarMatrix
const dimensionedScalar G
Newtonian constant of gravitation.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
virtual tmp< fvScalarMatrix > kSource() const
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
Generic dimensioned Type class.
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
zeroField SuSp(const Foam::zero, const GeometricField< Type, faPatchField, areaMesh > &)
A no-op source.
Eddy viscosity turbulence model base class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
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.
const dimensionSet dimVolume(pow3(dimLength))
Standard k-epsilon turbulence model for incompressible and compressible flows including rapid distort...
A class for handling words, derived from Foam::string.
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &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/v2312/OpenFOAM-v2312/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H" 1{ alphav=clamp((rho - rholSat)/(rhovSat - rholSat), zero_one{});alphal=1.0 - alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
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.
Base-class for all transport models used by the incompressible turbulence models. ...
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
virtual void correctNut()
A class for managing temporary objects.
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...
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
SymmTensor< Cmpt > devTwoSymm(const SymmTensor< Cmpt > &st)
Return the deviatoric part of twice the symmetric part of a SymmTensor.