Add a co2brine pvt model

The CO2BRINE model is activated by setting CO2STOR in the RUNSPEC section

The CO2 and brine pvt properties are computed based on pvt models in opm-material
- CO2 density is from Span and Wager (1996) as given in co2table.inc
- CO2 viscosity is from Fenhour et al (1998)
- Brine density and viscosity is based on H20 + correction based on Batzle and Wang (1992)
- H20 density is from Hu et al (2007)
- H20 viscosity is computed from density based on correlation given in IAPWS 97

At the current stage the oil-phase is used to model the brine. If a proper gas-water
simulator is made, the Brine PVT should be moved to the water phase.

Known limitations:
- Currently the viscosity of the Brine does not depend on the composition
- Salinity is assumed to be constant and given by SALINITY [mol/kg].

CMakeLists.txt

@@ -103,6 +103,7 @@ opm_add_test(test_ConditionalStorage)
 opm_add_test(test_eclblackoilfluidsystem CONDITION HAVE_ECL_INPUT)
 opm_add_test(test_eclblackoilpvt CONDITION HAVE_ECL_INPUT)
 opm_add_test(test_eclmateriallawmanager CONDITION HAVE_ECL_INPUT)
+opm_add_test(test_co2brinepvt CONDITION HAVE_ECL_INPUT)
 opm_add_test(test_fluidmatrixinteractions)
 opm_add_test(test_pengrobinson)
 opm_add_test(test_densead)

opm/material/binarycoefficients/Brine_CO2.hpp

@@ -28,9 +28,6 @@
 #ifndef OPM_BINARY_COEFF_BRINE_CO2_HPP
 #define OPM_BINARY_COEFF_BRINE_CO2_HPP
 
-#include <opm/material/components/Brine.hpp>
-#include <opm/material/components/H2O.hpp>
-#include <opm/material/components/CO2.hpp>
 #include <opm/material/IdealGas.hpp>
 
 namespace Opm {
@@ -40,10 +37,8 @@ namespace BinaryCoeff {
  * \ingroup Binarycoefficients
  * \brief Binary coefficients for brine and CO2.
  */
-template<class Scalar, class CO2Tables, bool verbose = true>
+template<class Scalar, class H2O, class CO2, bool verbose = true>
 class Brine_CO2 {
-    typedef Opm::H2O<Scalar> H2O;
-    typedef Opm::CO2<Scalar, CO2Tables> CO2;
     typedef Opm::IdealGas<Scalar> IdealGas;
     static const int liquidPhaseIdx = 0; // index of the liquid phase
     static const int gasPhaseIdx = 1; // index of the gas phase