36 namespace phaseChangeTwoPhaseMixtures
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_/(0.5*
sqr(UInf_)*tInf_)),
74 return Pair<tmp<volScalarField>>
76 mcCoeff_*
max(
p - pSat(), p0_),
77 mvCoeff_*
min(
p - pSat(), p0_)
87 return Pair<tmp<volScalarField>>
89 mcCoeff_*(1.0 - limitedAlpha1)*
pos0(
p - pSat()),
90 (-mvCoeff_)*limitedAlpha1*
neg(
p - pSat())
103 phaseChangeTwoPhaseMixtureCoeffs_ = optionalSubDict(
type() +
"Coeffs");
105 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"UInf", UInf_);
106 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"tInf", tInf_);
107 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cc", Cc_);
108 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cv", Cv_);
110 mcCoeff_ = Cc_/(0.5*
sqr(UInf_)*tInf_);
111 mvCoeff_ = Cv_*
rho1()/(0.5*
sqr(UInf_)*tInf_*
rho2());
virtual bool read()=0
Read the transportProperties dictionary and update.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
virtual Pair< tmp< volScalarField > > mDotP() const
Return the mass condensation and vaporisation rates as coefficients.
const dimensionSet dimless
Dimensionless.
GeometricField< vector, fvPatchField, volMesh > volVectorField
Merkle(const volVectorField &U, const surfaceScalarField &phi)
Construct from components.
dimensionedScalar neg(const dimensionedScalar &ds)
Macros for easy insertion into run-time selection tables.
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.
virtual void correct()
Correct the Merkle phaseChange model.
An ordered pair of two objects of type <T> with first() and second() elements.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedScalar pos0(const dimensionedScalar &ds)
defineTypeNameAndDebug(combustionModel, 0)
virtual bool read()
Read the transportProperties dictionary and update.
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
virtual void correct()
Correct the Kunz phaseChange model.
virtual Pair< tmp< volScalarField > > mDotAlphal() const
Return the mass condensation and vaporisation rates as a.
dimensionSet clamp(const dimensionSet &a, const dimensionSet &range)
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity