34 namespace phaseChangeTwoPhaseMixtures
39 phaseChangeTwoPhaseMixture,
53 const word& alpha1Name
56 phaseChangeTwoPhaseMixture(typeName, U, phi, alpha1Name),
58 n_(phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"n")),
59 dNuc_(phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"dNuc")),
60 Cc_(phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"Cc")),
61 Cv_(phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"Cv")),
63 p0_(
"0",
pSat().dimensions(), 0.0)
72 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::rRb
80 *limitedAlpha1/(1.0 + alphaNuc() - limitedAlpha1),
87 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::alphaNuc()
const
90 return Vnuc/(1 + Vnuc);
95 Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer::pCoeff
102 (limitedAlpha1*
rho1() + (scalar(1) - limitedAlpha1)*
rho2());
118 return Pair<tmp<volScalarField> >
120 Cc_*limitedAlpha1*pCoeff*
max(p -
pSat(), p0_),
122 Cv_*(1.0 + alphaNuc() - limitedAlpha1)*pCoeff*
min(p -
pSat(), p0_)
135 return Pair<tmp<volScalarField> >
137 Cc_*(1.0 - limitedAlpha1)*
pos(p -
pSat())*apCoeff,
139 (-Cv_)*(1.0 + alphaNuc() - limitedAlpha1)*
neg(p -
pSat())*apCoeff
152 phaseChangeTwoPhaseMixtureCoeffs_ = subDict(
type() +
"Coeffs");
154 phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"n") >> n_;
155 phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"dNuc") >> dNuc_;
156 phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"Cc") >> Cc_;
157 phaseChangeTwoPhaseMixtureCoeffs_.lookup(
"Cv") >> Cv_;