42 template<
class BasicTurbulenceModel>
45 this->nut_ = Ck_*
sqrt(this->
k())*this->
delta();
46 this->nut_.correctBoundaryConditions();
49 BasicTurbulenceModel::correctNut();
55 template<
class BasicTurbulenceModel>
58 const alphaField&
alpha,
64 const word& propertiesName,
117 if (
type == typeName)
127 template<
class BasicTurbulenceModel>
132 Ck_.readIfPresent(this->coeffDict());
133 Cm_.readIfPresent(this->coeffDict());
134 Ce_.readIfPresent(this->coeffDict());
135 Cs_.readIfPresent(this->coeffDict());
144 template<
class BasicTurbulenceModel>
147 if (!this->turbulence_)
153 const alphaField&
alpha = this->alpha_;
154 const rhoField&
rho = this->rho_;
160 ReynoldsStress<LESModel<BasicTurbulenceModel>>
::correct();
171 tmp<fvSymmTensorMatrix> REqn
188 this->boundNormalStress(
R);
void correct(GeometricField< Type, PatchField, GeoMesh > &field)
Apply correction to field.
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
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. ...
GeometricField< tensor, fvPatchField, volMesh > volTensorField
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual void printCoeffs(const word &type)
Print model coefficients.
label k
Boltzmann constant.
Generic dimensioned Type class.
GeometricField< vector, fvPatchField, volMesh > volVectorField
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
bool read(const char *buf, int32_t &val)
Same as readInt32.
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.
virtual void correct()
Correct sub-grid stress, eddy-Viscosity and related properties.
static const Identity< scalar > I
A class for handling words, derived from Foam::string.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
void boundNormalStress(volSymmTensorField &R) const
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/chef2/andy/OpenFOAM/release/v2406/OpenFOAM-v2406/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()
#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.
Base-class for all transport models used by the incompressible turbulence models. ...
virtual bool read()
Read model coefficients if they have changed.
const dimensionedScalar & D
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: [].
Reynolds-stress turbulence model base class.