30#include "twoPhaseSystem.H"
36 const volScalarField&
alpha,
37 const volScalarField&
rho,
38 const volVectorField&
U,
39 const surfaceScalarField& alphaRhoPhi,
40 const surfaceScalarField&
phi,
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_);
143 dimensioned<symmTensor>(dimensionSet(0, 2, -2, 0, 0),
Zero)
151 tmp<volScalarField> tpPrime
156 exp(preAlphaExp_*(alpha_ - alphaMax_)),
162 tpPrime.ref().boundaryFieldRef();
166 if (!bpPrime[patchi].
coupled())
168 bpPrime[patchi] == 0;
179 tmp<surfaceScalarField> tpPrime
190 tpPrime.ref().boundaryFieldRef();
194 if (!bpPrime[patchi].
coupled())
196 bpPrime[patchi] == 0;
207 return devRhoReff(U_);
228 dimensioned<symmTensor>
230 rho_.dimensions()*dimensionSet(0, 2, -2, 0, 0),
Zero
242 return tmp<fvVectorMatrix>
247 rho_.dimensions()*dimensionSet(0, 4, -2, 0, 0)
GeometricBoundaryField< scalar, fvPatchField, volMesh > Boundary
Type of boundary fields.
static word groupName(StringType base, const word &group)
Create dot-delimited name.group string.
Particle-particle phase-pressure RAS model.
virtual tmp< volScalarField > epsilon() const
Return the turbulence kinetic energy dissipation rate.
virtual tmp< volSymmTensorField > R() const
Return the Reynolds stress tensor.
virtual void correct()
Solve the kinetic theory equations and correct the viscosity.
virtual tmp< volScalarField > k() const
Return the turbulence kinetic energy.
virtual tmp< volScalarField > pPrime() const
Return the phase-pressure'.
virtual tmp< surfaceScalarField > pPrimef() const
Return the face-phase-pressure'.
virtual tmp< volScalarField > omega() const
Return the specific dissipation rate.
virtual tmp< volSymmTensorField > devRhoReff() const
Return the effective stress tensor.
virtual bool read()
Re-read model coefficients if they have changed.
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
bool readIfPresent(const dictionary &dict)
tmp< volVectorField > divDevRhoReff()
Return the effective viscous stress (laminar + turbulent).
A class for managing temporary objects.
bool coupled(solutionDict.getOrDefault("coupledEnergyField", false))
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
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.
const dimensionSet dimPressure
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedScalar exp(const dimensionedScalar &ds)
GeometricField< vector, fvPatchField, volMesh > volVectorField
fvMatrix< vector > fvVectorMatrix
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
static constexpr const zero Zero
Global zero (0)
#define forAll(list, i)
Loop across all elements in list.