Commit d3d32d32 authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[material]

- enthalpy function in the h2oairfs now use the given phaseIdx and not
automatically the nPhaseIdx, this could affect local-non-equilibrium
results
- fixed one description and one unit given

reviewed by nicolas


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@15048 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent e2f8db60
......@@ -167,7 +167,7 @@ public:
};
/*!
* \brief Specific enthalpy of liquid water \f$\mathrm{[J/kg]}\f$
* \brief Specific enthalpy of air \f$\mathrm{[J/kg]}\f$
* with 273.15 K as basis.
* See:
* W. Kays, M. Crawford, B. Weigand
......
......@@ -296,7 +296,7 @@ public:
{ density_[phaseIdx] = value; }
/*!
* \brief Set the specific enthalpy of a phase [J/m^3]
* \brief Set the specific enthalpy of a phase [J/kg]
*/
void setEnthalpy(int phaseIdx, Scalar value)
{ enthalpy_[phaseIdx] = value; }
......@@ -327,7 +327,7 @@ public:
Valgrind::CheckDefined(pressure_[phaseIdx]);
Valgrind::CheckDefined(saturation_[phaseIdx]);
Valgrind::CheckDefined(density_[phaseIdx]);
//Valgrind::CheckDefined(enthalpy_[phaseIdx]);
Valgrind::CheckDefined(enthalpy_[phaseIdx]);
Valgrind::CheckDefined(viscosity_[phaseIdx]);
}
......
......@@ -660,8 +660,8 @@ public:
int phaseIdx,
int componentIdx)
{
Scalar T = fluidState.temperature(nPhaseIdx);
Scalar p = fluidState.pressure(nPhaseIdx);
Scalar T = fluidState.temperature(phaseIdx);
Scalar p = fluidState.pressure(phaseIdx);
Valgrind::CheckDefined(T);
Valgrind::CheckDefined(p);
......@@ -671,7 +671,7 @@ public:
}
else if (phaseIdx == nPhaseIdx)
{
if (componentIdx == H2OIdx)
if (componentIdx == H2OIdx)
{
return H2O::gasEnthalpy(T, p);
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment