30 #include "twoPhaseSystem.H" 34 Foam::RASModels::phasePressureModel::phasePressureModel
41 const transportModel& phase,
42 const word& propertiesName,
48 RASModel<EddyDiffusivity<ThermalDiffusivity
50 PhaseCompressibleTurbulenceModel<phaseModel>
64 alphaMax_(coeffDict_.
get<scalar>(
"alphaMax")),
65 preAlphaExp_(coeffDict_.
get<scalar>(
"preAlphaExp")),
66 expMax_(coeffDict_.
get<scalar>(
"expMax")),
86 RASModel<EddyDiffusivity<ThermalDiffusivity
88 PhaseCompressibleTurbulenceModel<phaseModel>
93 coeffDict().readEntry(
"alphaMax", alphaMax_);
94 coeffDict().readEntry(
"preAlphaExp", preAlphaExp_);
95 coeffDict().readEntry(
"expMax", expMax_);
96 g0_.readIfPresent(coeffDict());
137 dimensioned<symmTensor>(dimensionSet(0, 2, -2, 0, 0),
Zero)
145 tmp<volScalarField> tpPrime
150 exp(preAlphaExp_*(alpha_ - alphaMax_)),
156 tpPrime.ref().boundaryFieldRef();
160 if (!bpPrime[patchi].
coupled())
162 bpPrime[patchi] == 0;
173 tmp<surfaceScalarField> tpPrime
184 tpPrime.ref().boundaryFieldRef();
188 if (!bpPrime[patchi].
coupled())
190 bpPrime[patchi] == 0;
201 return devRhoReff(U_);
216 dimensioned<symmTensor>
218 rho_.dimensions()*dimensionSet(0, 2, -2, 0, 0),
Zero 233 rho_.dimensions()*dimensionSet(0, 4, -2, 0, 0)
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
GeometricBoundaryField< scalar, fvPatchField, volMesh > Boundary
Type of boundary fields.
GeometricField< vector, fvPatchField, volMesh > volVectorField
bool read(const char *buf, int32_t &val)
Same as readInt32.
#define forAll(list, i)
Loop across all elements in list.
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.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
dimensionedScalar exp(const dimensionedScalar &ds)
virtual tmp< volScalarField > k() const
Return the turbulence kinetic energy.
static tmp< T > New(Args &&... args)
Construct tmp with forwarding arguments.
virtual void correct()
Solve the kinetic theory equations and correct the viscosity.
virtual void printCoeffs(const word &type)
Print model coefficients.
virtual tmp< fvVectorMatrix > divDevRhoReff(volVectorField &U) const
Return the source term for the momentum equation.
const dimensionSet dimPressure
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
virtual tmp< volScalarField > pPrime() const
Return the phase-pressure'.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
static tmp< GeometricField< symmTensor, fvPatchField, volMesh > > New(const word &name, IOobjectOption::registerOption regOpt, const Mesh &mesh, const dimensionSet &dims, const word &patchFieldType=fvPatchField< symmTensor >::calculatedType())
Return tmp field (NO_READ, NO_WRITE) from name, mesh, dimensions and patch type. [Takes current timeN...
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
RASModel< EddyDiffusivity< turbulenceModel > > RASModel
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
virtual tmp< volSymmTensorField > devRhoReff() const
Return the effective stress tensor.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
A class for managing temporary objects.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
Do not request registration (bool: false)
const dimensionSet & dimensions() const noexcept
Return dimensions.
virtual tmp< surfaceScalarField > pPrimef() const
Return the face-phase-pressure'.
static constexpr const zero Zero
Global zero (0)
virtual bool read()
Re-read model coefficients if they have changed.