My Project
H2O_N2.hpp
Go to the documentation of this file.
1 // -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
2 // vi: set et ts=4 sw=4 sts=4:
3 /*
4  This file is part of the Open Porous Media project (OPM).
5 
6  OPM is free software: you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation, either version 2 of the License, or
9  (at your option) any later version.
10 
11  OPM is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with OPM. If not, see <http://www.gnu.org/licenses/>.
18 
19  Consult the COPYING file in the top-level source directory of this
20  module for the precise wording of the license and the list of
21  copyright holders.
22 */
27 #ifndef OPM_BINARY_COEFF_H2O_N2_HPP
28 #define OPM_BINARY_COEFF_H2O_N2_HPP
29 
30 #include "HenryIapws.hpp"
31 #include "FullerMethod.hpp"
32 
35 
36 namespace Opm {
37 namespace BinaryCoeff {
38 
43 class H2O_N2
44 {
45 public:
51  template <class Evaluation>
52  static Evaluation henry(const Evaluation& temperature)
53  {
54  const double E = 2388.8777;
55  const double F = -14.9593;
56  const double G = 42.0179;
57  const double H = -29.4396;
58 
59  return henryIAPWS(E, F, G, H, temperature);
60  }
61 
69  template <class Evaluation>
70  static Evaluation gasDiffCoeff(const Evaluation& temperature, const Evaluation& pressure)
71  {
72  typedef H2O<double> H2O;
73  typedef N2<double> N2;
74 
75  // atomic diffusion volumes
76  const double SigmaNu[2] = { 13.1 /* H2O */, 18.5 /* N2 */ };
77  // molar masses [g/mol]
78  const double M[2] = { H2O::molarMass()*1e3, N2::molarMass()*1e3 };
79 
80  return fullerMethod(M, SigmaNu, temperature, pressure);
81  }
82 
101  template <class Evaluation>
102  static Evaluation liquidDiffCoeff(const Evaluation& temperature, const Evaluation& /*pressure*/)
103  {
104  const double Texp = 273.15 + 25; // [K]
105  const double Dexp = 2.01e-9; // [m^2/s]
106 
107  return Dexp * temperature/Texp;
108  }
109 };
110 
111 } // namespace BinaryCoeff
112 } // namespace Opm
113 
114 #endif
Evaluation fullerMethod(const Scalar *M, const Scalar *SigmaNu, const Evaluation &temperature, const Evaluation &pressure)
Estimate binary diffusion coefficents in gases according to the method by Fuller.
Definition: FullerMethod.hpp:56
Material properties of pure water .
The IAPWS formulation of Henry coefficients in water.
Evaluation henryIAPWS(Scalar E, Scalar F, Scalar G, Scalar H, const Evaluation &temperature)
The Henry constants in liquid water using the IAPWS 2004 formulation.
Definition: HenryIapws.hpp:45
Properties of pure molecular nitrogen .
Binary coefficients for water and nitrogen.
Definition: H2O_N2.hpp:44
static Evaluation liquidDiffCoeff(const Evaluation &temperature, const Evaluation &)
Diffusion coefficent for molecular nitrogen in liquid water.
Definition: H2O_N2.hpp:102
static Evaluation gasDiffCoeff(const Evaluation &temperature, const Evaluation &pressure)
Binary diffusion coefficent for molecular water and nitrogen.
Definition: H2O_N2.hpp:70
static Evaluation henry(const Evaluation &temperature)
Henry coefficent for molecular nitrogen in liquid water.
Definition: H2O_N2.hpp:52
Material properties of pure water .
Definition: H2O.hpp:62
static const Scalar molarMass()
The molar mass in of water.
Definition: H2O.hpp:80
Properties of pure molecular nitrogen .
Definition: N2.hpp:49
static Scalar molarMass()
The molar mass in of molecular nitrogen.
Definition: N2.hpp:62