Documentation in header of humidityTemperatureCoupledMixedFvPatchScalarField.H possibly wrong
In the header of the file src/thermophysicalModels/thermophysicalPropertiesFvPatchFields/liquidProperties/humidityTemperatureCoupledMixed/humidityTemperatureCoupledMixedFvPatchScalarField.H
the documentation might by incorrect:
Currently it states:
The mass transfer correlation used is:
\f[ h_m = D_{ab} \frac{Sc}{L} \f]
where:
\vartable
D_{ab} | mass vapour difussivity
L | characteristic length
Sc | Schmidt number
\endvartable
The Schmidt number is calculated using:
\f{eqnarray*}{
0.664 Re^\frac{1}{2} Sc^\frac{1}{3} & Re < 5.0E+05 \\
0.037 Re^\frac{4}{5} Sc^\frac{1}{3} & Re > 5.0E+05
\f}
But in the first equation, not the Schmidt number Sc should be used, but the Sherwood number Sh. And Sh is calculated by the correlations depending on the Schmidt number Sc and the Reynolds number Re. This is also the way, it is used in the code:
Calculation of Schmidt number and its use in the calculation of the Sherwood number:
const scalar Sc = nuf/Dab;
const scalar Sh = this->Sh(Re, Sc);
Calculation of Sherwood number:
Foam::scalar Foam::humidityTemperatureCoupledMixedFvPatchScalarField::Sh
(
const scalar Re,
const scalar Sc
) const
{
if (Re < 5.0E+05)
{
return 0.664*sqrt(Re)*cbrt(Sc);
}
else
{
return 0.037*pow(Re, 0.8)*cbrt(Sc);
}
}
Use of Sherwood number in the calculation of the mass transfer:
// Mass transfer coefficient [m/s]
const scalar hm = Dab*Sh/L_;
So the header should read:
The mass transfer correlation used is:
\f[ h_m = D_{ab} \frac{Sh}{L} \f]
where:
\vartable
D_{ab} | mass vapour difussivity
L | characteristic length
Sh | Sherwood number
\endvartable
The Sherwood number is calculated using:
\f{eqnarray*}{
0.664 Re^\frac{1}{2} Sc^\frac{1}{3} & Re < 5.0E+05 \\
0.037 Re^\frac{4}{5} Sc^\frac{1}{3} & Re > 5.0E+05
\f}
where:
\vartable
Re | Reynolds number
Sc | Schmidt number
\endvartable
(There is another typo: A space is missing between in \c p (currently \cp) in this line:
- requires \c rho, \c thickness and \cp
)