Commit 584f9970 authored by Katherina Baber's avatar Katherina Baber
Browse files

make calculation of gas mixture enthalpy consistent with the other

fluidsystems: call component enthalpies with phase pressure not with
partial pressure

this was discussed with Klaus, Holger, Philipp and Bernd. The
enthalpies can be called with the partial pressure and be summed up directly,
or they can be called  with the phase pressure and are then weighted by
the mass fractions. Klaus compared the two possibilities and found only
minor differences.


git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12402 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 1b072751
......@@ -630,12 +630,10 @@ public:
// "partial specific enthalpies" of the components.
Scalar hH2O =
fluidState.massFraction(nPhaseIdx, H2OIdx)
* H2O::gasEnthalpy(T,
p*fluidState.moleFraction(nPhaseIdx, H2OIdx));
* H2O::gasEnthalpy(T, p);
Scalar hN2 =
fluidState.massFraction(nPhaseIdx, N2Idx)
* N2::gasEnthalpy(T,
p*fluidState.moleFraction(nPhaseIdx, N2Idx));
* N2::gasEnthalpy(T, p);
return hH2O + hN2;
}
}
......
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