My Project
SimpleCO2.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 */
30 #ifndef OPM_SIMPLE_CO2_HPP
31 #define OPM_SIMPLE_CO2_HPP
32 
35 
37 
38 #include <cmath>
39 
40 namespace Opm {
41 
49 template <class Scalar>
50 class SimpleCO2 : public Component<Scalar, SimpleCO2<Scalar> >
51 {
52  typedef ::Opm::IdealGas<Scalar> IdealGas;
53 
54 public:
58  static const char* name()
59  { return "CO2"; }
60 
64  static Scalar molarMass()
65  { return 44e-3; }
66 
70  static Scalar criticalTemperature()
71  { return 273.15 + 30.95; /* [K] */ }
72 
76  static Scalar criticalPressure()
77  { return 73.8e5; /* [N/m^2] */ }
78 
82  static Scalar tripleTemperature()
83  { return 273.15 - 56.35; /* [K] */ }
84 
88  static Scalar acentricFactor() { return 0.224; }
89 
93  // Critical volume [m3/kmol]
94  static Scalar criticalVolume() {return 9.412e-5; }
95 
99  static Scalar triplePressure()
100  { return 5.11e5; /* [N/m^2] */ }
101 
105  static bool gasIsCompressible()
106  { return true; }
107 
111  static bool gasIsIdeal()
112  { return true; }
113 
117  template <class Evaluation>
118  static Evaluation gasEnthalpy(const Evaluation& temperature,
119  const Evaluation&)
120  { return 350.0e3 + temperature*0.85e3; }
121 
125  template <class Evaluation>
126  static Evaluation gasHeatCapacity(const Evaluation&,
127  const Evaluation&)
128  { return 0.85e3; }
129 
133  template <class Evaluation>
134  static Evaluation liquidEnthalpy(const Evaluation& temperature,
135  const Evaluation&)
136  { return temperature*2e3; }
137 
141  template <class Evaluation>
142  static Evaluation liquidHeatCapacity(const Evaluation&,
143  const Evaluation&)
144  { return 2e3; /* TODO: UNKNOWN! */ }
145 
149  template <class Evaluation>
150  static Evaluation gasInternalEnergy(const Evaluation& temperature,
151  const Evaluation& pressure)
152  {
153  return
154  gasEnthalpy(temperature, pressure) -
155  1/molarMass()* // conversion from [J/(mol K)] to [J/(kg K)]
156  IdealGas::R*temperature; // = pressure * spec. volume for an ideal gas
157  }
158 
162  template <class Evaluation>
163  static Evaluation gasDensity(const Evaluation& temperature, const Evaluation& pressure)
164  {
165  // Assume an ideal gas
166  return IdealGas::density(Evaluation(molarMass()), temperature, pressure);
167  }
168 
177  template <class Evaluation>
178  static Evaluation gasViscosity(const Evaluation& temperature, const Evaluation& /*pressure*/)
179  {
180  const Scalar Tc = criticalTemperature();
181  const Scalar Vc = 93.9; // critical specific volume [cm^3/mol]
182  const Scalar omega = 0.239; // accentric factor
183  const Scalar M = molarMass() * 1e3; // molar mas [g/mol]
184  const Scalar dipole = 0.0; // dipole moment [debye]
185 
186  Scalar mu_r4 = 131.3 * dipole / std::sqrt(Vc * Tc);
187  mu_r4 *= mu_r4;
188  mu_r4 *= mu_r4;
189 
190  Scalar Fc = 1 - 0.2756*omega + 0.059035*mu_r4;
191  Evaluation Tstar = 1.2593 * temperature/Tc;
192  Evaluation Omega_v =
193  1.16145*pow(Tstar, -0.14874) +
194  0.52487*exp(- 0.77320*Tstar) +
195  2.16178*exp(- 2.43787*Tstar);
196  Evaluation mu = 40.785*Fc*sqrt(M*temperature)/(std::pow(Vc, 2./3)*Omega_v);
197 
198  // convertion from micro poise to Pa s
199  return mu/1e6 / 10;
200  }
201 };
202 
203 } // namespace Opm
204 
205 #endif
Abstract base class of a pure chemical species.
Relations valid for an ideal gas.
A number of commonly used algebraic functions for the localized OPM automatic differentiation (AD) fr...
Abstract base class of a pure chemical species.
Definition: Component.hpp:42
Relations valid for an ideal gas.
Definition: IdealGas.hpp:38
static const Scalar R
The ideal gas constant .
Definition: IdealGas.hpp:41
static Evaluation density(const Evaluation &avgMolarMass, const Evaluation &temperature, const Evaluation &pressure)
The density of the gas in , depending on pressure, temperature and average molar mass of the gas.
Definition: IdealGas.hpp:48
A simplistic class representing the fluid properties.
Definition: SimpleCO2.hpp:51
static Scalar criticalVolume()
Critical volume of [m2/kmol].
Definition: SimpleCO2.hpp:94
static Evaluation liquidHeatCapacity(const Evaluation &, const Evaluation &)
Specific isobaric heat capacity of the component [J/kg] as a liquid.
Definition: SimpleCO2.hpp:142
static Evaluation gasHeatCapacity(const Evaluation &, const Evaluation &)
Specific isobaric heat capacity of the component [J/kg] as a gas.
Definition: SimpleCO2.hpp:126
static const char * name()
A human readable name for the component.
Definition: SimpleCO2.hpp:58
static Scalar criticalTemperature()
Returns the critical temperature of .
Definition: SimpleCO2.hpp:70
static Scalar triplePressure()
Returns the pressure at the triple point of .
Definition: SimpleCO2.hpp:99
static Scalar acentricFactor()
Acentric factor of .
Definition: SimpleCO2.hpp:88
static Evaluation liquidEnthalpy(const Evaluation &temperature, const Evaluation &)
Specific enthalpy of the pure component in liquid.
Definition: SimpleCO2.hpp:134
static bool gasIsCompressible()
Returns true iff the gas phase is assumed to be compressible.
Definition: SimpleCO2.hpp:105
static bool gasIsIdeal()
Returns true iff the gas phase is assumed to be ideal.
Definition: SimpleCO2.hpp:111
static Scalar tripleTemperature()
Returns the temperature at the triple point of .
Definition: SimpleCO2.hpp:82
static Scalar criticalPressure()
Returns the critical pressure of .
Definition: SimpleCO2.hpp:76
static Evaluation gasViscosity(const Evaluation &temperature, const Evaluation &)
The dynamic viscosity of the pure component at a given pressure in and temperature in .
Definition: SimpleCO2.hpp:178
static Evaluation gasEnthalpy(const Evaluation &temperature, const Evaluation &)
Specific enthalpy of the pure component in gas.
Definition: SimpleCO2.hpp:118
static Evaluation gasInternalEnergy(const Evaluation &temperature, const Evaluation &pressure)
Specific internal energy of the pure component in gas.
Definition: SimpleCO2.hpp:150
static Evaluation gasDensity(const Evaluation &temperature, const Evaluation &pressure)
The density in of the component at a given pressure in and temperature in .
Definition: SimpleCO2.hpp:163
static Scalar molarMass()
The molar mass in of the component.
Definition: SimpleCO2.hpp:64