36namespace phaseChangeTwoPhaseMixtures
47 const volVectorField&
U,
48 const surfaceScalarField&
phi
51 phaseChangeTwoPhaseMixture(typeName,
U,
phi),
53 UInf_(
"UInf",
dimVelocity, phaseChangeTwoPhaseMixtureCoeffs_),
54 tInf_(
"tInf",
dimTime, phaseChangeTwoPhaseMixtureCoeffs_),
55 Cc_(
"Cc",
dimless, phaseChangeTwoPhaseMixtureCoeffs_),
56 Cv_(
"Cv",
dimless, phaseChangeTwoPhaseMixtureCoeffs_),
58 p0_(pSat().dimensions(), Zero),
60 mcCoeff_(Cc_*
rho2()/tInf_),
75 return Pair<tmp<volScalarField>>
77 mcCoeff_*
sqr(limitedAlpha1)
90 return Pair<tmp<volScalarField>>
92 mcCoeff_*
sqr(limitedAlpha1)*(1.0 - limitedAlpha1)
93 *
pos0(
p - pSat())/
max(
p - pSat(), 0.01*pSat()),
95 (-mvCoeff_)*limitedAlpha1*
neg(
p - pSat())
108 phaseChangeTwoPhaseMixtureCoeffs_ = optionalSubDict(
type() +
"Coeffs");
110 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"UInf", UInf_);
111 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"tInf", tInf_);
112 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cc", Cc_);
113 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cv", Cv_);
115 mcCoeff_ = Cc_*
rho2()/tInf_;
116 mvCoeff_ = Cv_*
rho2()/(0.5*
rho1()*
sqr(UInf_)*tInf_);
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
An ordered pair of two objects of type <T> with first() and second() elements.
const Type & lookupObject(const word &name, const bool recursive=false) const
virtual bool read()=0
Read the transportProperties dictionary and update.
const dimensionedScalar & pSat() const
Return const-access to the saturation vapour pressure.
Kunz cavitation model slightly modified so that the condensation term is switched off when the pressu...
virtual Pair< tmp< volScalarField > > mDotP() const
Return the mass condensation and vaporisation rates as coefficients.
virtual Pair< tmp< volScalarField > > mDotAlphal() const
Return the mass condensation and vaporisation rates as a.
virtual void correct()
Correct the Kunz phaseChange model.
virtual bool read()
Read the transportProperties dictionary and update.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
const dimensionSet dimless
Dimensionless.
dimensionedScalar pos0(const dimensionedScalar &ds)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const dimensionSet dimVelocity
fileName::Type type(const fileName &name, const bool followLink=true)
Return the file type: DIRECTORY or FILE, normally following symbolic links.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedScalar neg(const dimensionedScalar &ds)