37 namespace phaseChangeTwoPhaseMixtures
42 phaseChangeTwoPhaseMixture,
58 phaseChangeTwoPhaseMixture(typeName,
U,
phi),
61 dNuc_(
"dNuc",
dimLength, phaseChangeTwoPhaseMixtureCoeffs_),
62 Cc_(
"Cc",
dimless, phaseChangeTwoPhaseMixtureCoeffs_),
63 Cv_(
"Cv",
dimless, phaseChangeTwoPhaseMixtureCoeffs_),
65 p0_(pSat().dimensions(),
Zero)
74 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::rRb
82 *limitedAlpha1/(1.0 + alphaNuc() - limitedAlpha1),
89 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::alphaNuc()
const 92 return Vnuc/(1 + Vnuc);
97 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::pCoeff
105 limitedAlpha1*
rho1() + (scalar(1) - limitedAlpha1)*
rho2()
110 *rRb(limitedAlpha1)/(
rho*
sqrt(
mag(
p - pSat()) + 0.01*pSat()));
122 return Pair<tmp<volScalarField>>
124 Cc_*limitedAlpha1*pCoeff*
max(
p - pSat(), p0_),
126 Cv_*(1.0 + alphaNuc() - limitedAlpha1)*pCoeff*
min(
p - pSat(), p0_)
140 return Pair<tmp<volScalarField>>
142 Cc_*(1.0 - limitedAlpha1)*
pos0(
p - pSat())*apCoeff,
144 (-Cv_)*(1.0 + alphaNuc() - limitedAlpha1)*
neg(
p - pSat())*apCoeff
157 phaseChangeTwoPhaseMixtureCoeffs_ = optionalSubDict(
type() +
"Coeffs");
159 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"n", n_);
160 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"dNuc", dNuc_);
161 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cc", Cc_);
162 phaseChangeTwoPhaseMixtureCoeffs_.readEntry(
"Cv", Cv_);
const Type & lookupObject(const word &name, const bool recursive=false) const
Lookup and return const reference to the object of the given Type. Fatal if not found or the wrong ty...
virtual bool read()=0
Read the transportProperties dictionary and update.
virtual void correct()
Correct the SchnerrSauer phaseChange model.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
virtual void correct()=0
Correct the phaseChange model.
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual Pair< tmp< volScalarField > > mDotP() const
Return the mass condensation and vaporisation rates as coefficients.
const dimensionSet dimless
Dimensionless.
GeometricField< vector, fvPatchField, volMesh > volVectorField
virtual bool read()
Read the transportProperties dictionary and update.
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.
const dimensionSet dimVolume(pow3(dimLength))
An ordered pair of two objects of type <T> with first() and second() elements.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
constexpr scalar pi(M_PI)
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)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
virtual Pair< tmp< volScalarField > > mDotAlphal() const
Return the mass condensation and vaporisation rates as a.
A class for managing temporary objects.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
SchnerrSauer(const volVectorField &U, const surfaceScalarField &phi)
Construct from components.
dimensionSet clamp(const dimensionSet &a, const dimensionSet &range)
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)
static constexpr const zero Zero
Global zero (0)