26 #include "phaseModel.H" 28 #include "diameterModel.H" 43 const phaseSystem& fluid,
44 const word& phaseName,
52 IOobject::groupName(
"alpha", phaseName),
55 IOobject::READ_IF_PRESENT,
69 fluid.subDict(phaseName).
lookup(
"residualAlpha")
71 alphaMax_(fluid.subDict(phaseName).lookupOrDefault(
"alphaMax", 1.0))
80 return autoPtr<phaseModel>(NULL);
118 return residualAlpha_;
130 return diameterModel_().d();
136 diameterModel_->correct();
158 return diameterModel_->read(fluid_.subDict(name_));
171 static tmp<Foam::volScalarField> divU_(NULL);
179 <<
"Attempt to set the dilatation rate of an incompressible phase" 200 <<
"Attempt to set the dilatation rate of an incompressible phase" intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
virtual void correctKinematics()
Correct the kinematics.
Ostream & endl(Ostream &os)
Add newline and flush stream.
Generic dimensioned Type class.
phaseModel(const word &phaseName, const volScalarField &p, const volScalarField &T)
Construct from components.
const dimensionedScalar & residualAlpha() const
Return the residual phase-fraction for given phase.
static autoPtr< diameterModel > New(const dictionary &dict, const phaseModel &phase)
autoPtr< phaseModel > clone() const
Return clone.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual ~phaseModel()
Destructor.
stressControl lookup("compactNormalStress") >> compactNormalStress
virtual bool read()
Read phase properties dictionary.
Class to represent a system of phases and model interfacial transfers between them.
A class for handling words, derived from string.
const word & keyword() const
void correct()
Correct the laminar viscosity.
defineRunTimeSelectionTable(reactionRateFlameArea, dictionary)
defineTypeNameAndDebug(combustionModel, 0)
tmp< volScalarField > d() const
const word & name() const
static const GeometricField< scalar, fvPatchField, volMesh > & null()
Return a null geometric field.
virtual const tmp< volScalarField > & divU() const
Return the phase dilatation rate (d(alpha)/dt + div(alpha*phi))
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
#define WarningInFunction
Report a warning using Foam::Warning.
label index() const
Return the index of the phase.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual const surfaceScalarField & DbyA() const
Return the phase diffusivity divided by the momentum coefficient.
virtual void correctThermo()
Correct the thermodynamics.
volScalarField divU(fvc::div(fvc::absolute(phi, U)))
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
A class for managing temporary objects.
virtual void correctEnergyTransport()
Correct the energy transport e.g. alphat.
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
virtual void correctTurbulence()
Correct the turbulence.
const phaseSystem & fluid() const
Return the system to which this phase belongs.
virtual bool compressible() const
Return true if the phase is compressible otherwise false.
virtual const volScalarField & K() const
Return the phase kinetic energy.
scalar alphaMax() const
Return the maximum phase-fraction (e.g. packing limit)