From 6862817f439b2d9a8a78f7445809b9bc9f73bb4b Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Fri, 12 May 2017 12:17:07 +0200
Subject: [PATCH] [h2oheavyoil] Fix return values of the Henry coefficients
(was the fugacitycoefficient before)
However, the results are unaffected.
---
dumux/material/fluidsystems/h2oheavyoil.hh | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/dumux/material/fluidsystems/h2oheavyoil.hh b/dumux/material/fluidsystems/h2oheavyoil.hh
index fc271d61f8..bd811f3b7e 100644
--- a/dumux/material/fluidsystems/h2oheavyoil.hh
+++ b/dumux/material/fluidsystems/h2oheavyoil.hh
@@ -326,13 +326,12 @@ public:
assert(0 <= compIdx && compIdx < numComponents);
const Scalar T = fluidState.temperature(phaseIdx);
- const Scalar p = fluidState.pressure(phaseIdx);
if (compIdx == NAPLIdx && phaseIdx == wPhaseIdx)
- return Dumux::BinaryCoeff::H2O_HeavyOil::henryOilInWater(T)/p;
+ return Dumux::BinaryCoeff::H2O_HeavyOil::henryOilInWater(T);
else if (phaseIdx == nPhaseIdx && compIdx == H2OIdx)
- return Dumux::BinaryCoeff::H2O_HeavyOil::henryWaterInOil(T)/p;
+ return Dumux::BinaryCoeff::H2O_HeavyOil::henryWaterInOil(T);
else
DUNE_THROW(Dune::InvalidStateException, "non-existent henry coefficient for phase index " << phaseIdx
--
GitLab