00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00027 #ifndef OPM_BINARY_COEFF_H2O_AIR_HPP
00028 #define OPM_BINARY_COEFF_H2O_AIR_HPP
00029
00030 #include <opm/material/common/MathToolbox.hpp>
00031
00032 #include <cmath>
00033
00034 namespace Opm {
00035 namespace BinaryCoeff {
00036
00041 class H2O_Air
00042 {
00043 public:
00054 template <class Evaluation>
00055 static Evaluation henry(const Evaluation& temperature)
00056 { return 1.0/((0.8942+1.47*Opm::exp(-0.04394*(temperature-273.15)))*1e-10); }
00057
00069 template <class Evaluation>
00070 static Evaluation gasDiffCoeff(const Evaluation& temperature, const Evaluation& pressure)
00071 {
00072 double Theta=1.8;
00073 double Daw=2.13e-5;
00074 double pg0=1.e5;
00075 double T0=273.15;
00076
00077 return Daw*(pg0/pressure)*Opm::pow((temperature/T0),Theta);
00078 }
00079
00100 template <class Evaluation>
00101 static Evaluation liquidDiffCoeff(const Evaluation& temperature, const Evaluation& )
00102 {
00103 const double Texp = 273.15 + 25;
00104 const double Dexp = 2.01e-9;
00105 return Dexp/Texp*temperature;
00106 }
00107 };
00108
00109 }
00110 }
00111
00112 #endif