38namespace temperaturePhaseChangeTwoPhaseMixtures
43 temperaturePhaseChangeTwoPhaseMixture,
54 const thermoIncompressibleTwoPhaseMixture&
mixture,
58 temperaturePhaseChangeTwoPhaseMixture(
mixture,
mesh),
63 optionalSubDict(
type() +
"Coeffs")
69 optionalSubDict(
type() +
"Coeffs")
81 const twoPhaseMixtureEThermo&
thermo =
82 refCast<const twoPhaseMixtureEThermo>
91 return Pair<tmp<volScalarField>>
105 min(
max(mixture_.alpha1(), scalar(0)), scalar(1))
110 min(
max(mixture_.alpha2(), scalar(0)), scalar(1))
115 const twoPhaseMixtureEThermo&
thermo =
116 refCast<const twoPhaseMixtureEThermo>
127 "mDotE", coeffE_*mixture_.rho1()*limitedAlpha1*
max(
T - TSat,
T0)
131 "mDotC", coeffC_*mixture_.rho2()*limitedAlpha2*
max(TSat -
T,
T0)
134 if (mesh_.time().outputTime())
140 return Pair<tmp<volScalarField>>
153 min(
max(mixture_.alpha1(), scalar(0)), scalar(1))
158 min(
max(mixture_.alpha2(), scalar(0)), scalar(1))
163 const twoPhaseMixtureEThermo&
thermo =
164 refCast<const twoPhaseMixtureEThermo>
171 return Pair<tmp<volScalarField>>
173 coeffC_*mixture_.rho2()*limitedAlpha2*
pos(TSat -
T),
174 coeffE_*mixture_.rho1()*limitedAlpha1*
pos(
T - TSat)
185 tmp<fvScalarMatrix> tTSource
196 const twoPhaseMixtureEThermo&
thermo =
197 refCast<const twoPhaseMixtureEThermo>
208 min(
max(mixture_.alpha1(), scalar(0)), scalar(1))
213 min(
max(mixture_.alpha2(), scalar(0)), scalar(1))
218 coeffE_*mixture_.rho1()*limitedAlpha1*
L*
pos(
T - TSat)
222 coeffC_*mixture_.rho2()*limitedAlpha2*
L*
pos(TSat -
T)
242 subDict(
type() +
"Coeffs").readEntry(
"coeffC", coeffC_);
243 subDict(
type() +
"Coeffs").readEntry(
"coeffE", coeffE_);
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
An ordered pair of two objects of type <T> with first() and second() elements.
static const word dictName
const dimensionedScalar & rho1() const
Return const-access to phase1 density.
const dimensionedScalar & rho2() const
Return const-access to phase2 density.
Constant dispersed-phase particle diameter model.
const Type & lookupObject(const word &name, const bool recursive=false) const
constant condensation/saturation model.
const fvMesh & mesh_
Reference to fvMesh.
const thermoIncompressibleTwoPhaseMixture & mixture_
Reference to the thermoIncompressibleTwoPhaseMixture.
virtual bool read()
Read the transportProperties dictionary and update.
virtual Pair< tmp< volScalarField > > mDotAlphal() const
Return the mass condensation and vaporisation rates as a.
virtual Pair< tmp< volScalarField > > mDotDeltaT() const
Return the mass condensation and vaporisation rates as a.
virtual tmp< fvScalarMatrix > TSource() const
Source for T equarion.
virtual void correct()
Correct the constant phaseChange model.
virtual bool read()
Read the transportProperties dictionary and update.
virtual Pair< tmp< volScalarField > > mDot() const
Return the mass condensation and vaporisation rates as coefficients.
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
A class for managing temporary objects.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Calculate the gradient of the given field.
Calculate the finiteVolume matrix for implicit and explicit sources.
zeroField Sp(const Foam::zero, const GeometricField< Type, fvPatchField, volMesh > &)
A no-op source.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedScalar pos(const dimensionedScalar &ds)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
const dimensionSet dimless
Dimensionless.
const dimensionSet dimEnergy
fvMatrix< scalar > fvScalarMatrix
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
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 dimTemperature(0, 0, 0, 1, 0, 0, 0)
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)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
const vector L(dict.get< vector >("L"))