From 86a0e9befcab40cda022a2b31301fb114a17aa4b Mon Sep 17 00:00:00 2001 From: Andreas Lauser Date: Thu, 4 Jan 2018 15:25:11 +0100 Subject: [PATCH] ECL thermal gas PVT: add the heat of vaporization to the gas enthalpy this does only matters for temperature changes due to dissolution and evaporization of gas in oil or oil in gas. The comprehensive solution would be to make the enthalpy of the phases composition dependent, but ECL does not seem to support this when using the black-oil model. --- opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp b/opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp index ded59a872..62b322078 100644 --- a/opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp +++ b/opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp @@ -131,7 +131,14 @@ public: std::vector hSamples(temperatureColumn.size()); - Scalar h = temperatureColumn[0]*cpGasColumn[0]; + // the specific enthalpy of vaporization. since ECL does not seem to + // feature a proper way to specify this quantity, we use the value for + // methane. A proper model would also need to consider the enthalpy + // change due to dissolution, i.e. the enthalpies of the gas and oil + // phases should depend on the phase composition + const Scalar hVap = 480.6e3; // [J / kg] + + Scalar h = temperatureColumn[0]*cpGasColumn[0] + hVap; for (size_t i = 0;; ++i) { hSamples[i] = h;