3.1-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
simpleh2o.hh
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 * See the file COPYING for full copying permissions. *
5 * *
6 * This program 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 3 of the License, or *
9 * (at your option) any later version. *
10 * *
11 * This program 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 this program. If not, see <http://www.gnu.org/licenses/>. *
18 *****************************************************************************/
25#ifndef DUMUX_SIMPLE_H2O_HH
26#define DUMUX_SIMPLE_H2O_HH
27
29
30#include <cmath>
31
35
36namespace Dumux {
37namespace Components {
38
46template <class Scalar>
48: public Components::Base<Scalar, SimpleH2O<Scalar> >
49, public Components::Liquid<Scalar, SimpleH2O<Scalar> >
50, public Components::Gas<Scalar, SimpleH2O<Scalar> >
51{
53
54 static const Scalar R; // specific gas constant of water
55
56public:
60 static std::string name()
61 { return "SimpleH2O"; }
62
66 static constexpr Scalar molarMass()
67 { return 18e-3; }
68
73 { return 647.096; /* [K] */ }
74
79 { return 22.064e6; /* [N/m^2] */ }
80
85 { return 273.16; /* [K] */ }
86
91 { return 611.657; /* [N/m^2] */ }
92
106 {
107 if (T > criticalTemperature())
108 return criticalPressure();
109 if (T < tripleTemperature())
110 return 0; // water is solid: We don't take sublimation into account
111
112 static const Scalar n[10] = {
113 0.11670521452767e4, -0.72421316703206e6, -0.17073846940092e2,
114 0.12020824702470e5, -0.32325550322333e7, 0.14915108613530e2,
115 -0.48232657361591e4, 0.40511340542057e6, -0.23855557567849,
116 0.65017534844798e3
117 };
118
119 Scalar sigma = T + n[8]/(T - n[9]);
120
121 Scalar A = (sigma + n[0])*sigma + n[1];
122 Scalar B = (n[2]*sigma + n[3])*sigma + n[4];
123 Scalar C = (n[5]*sigma + n[6])*sigma + n[7];
124
125 using std::sqrt;
126 Scalar tmp = Scalar(2.0)*C/(sqrt(B*B - Scalar(4.0)*A*C) - B);
127 tmp *= tmp;
128 tmp *= tmp;
129
130 return Scalar(1e6)*tmp;
131 }
132
141 { return 1976*(temperature - 293.15) + 2.45e6; }
142
151 {
152 return 4180*(temperature - 293.15);
153 }
154
170 {
171 return
173 1/molarMass()* // conversion from [J/(mol K)] to [J/(kg K)]
174 IdealGas::R*temperature; // = pressure *spec. volume for an ideal gas
175 }
176
185 { return
188
192 static constexpr bool gasIsCompressible()
193 { return true; }
194
198 static constexpr bool liquidIsCompressible()
199 { return false; }
200
204 static constexpr bool gasViscosityIsConstant()
205 { return true; }
206
210 static constexpr bool liquidViscosityIsConstant()
211 { return true; }
212
220 {
221 // Assume an ideal gas
223 }
224
234
238 static constexpr bool gasIsIdeal()
239 { return true; }
240
248 {
249 // Assume an ideal gas
251 }
252
260 {
261 return 1000.;
262 }
263
273
281 {
282 DUNE_THROW(Dune::InvalidStateException,
283 "The liquid pressure is undefined for incompressible fluids");
284 }
285
293 static Scalar gasViscosity(Scalar temperature, Scalar pressure, bool regularize=true)
294 {
295 return 1e-05;
296 }
297
305 {
306 return 1e-03;
307 }
308
317 {
318 return 4.2e3;
319 }
320
329 {
330 return 0.679;
331 }
332
341 {
342 return 0.025;
343 }
344
353 {
354 return 2.08e3;
355 }
356
357};
358
359template <class Scalar>
360struct IsAqueous<SimpleH2O<Scalar>> : public std::true_type {};
361
362template <class Scalar>
364
365} // end namespace Components
366
367} // end namespace Dumux
368
369#endif
Interface for components that have a gas state.
Interface for components that have a liquid state.
Relations valid for an ideal gas.
make the local view function available whenever we use the grid geometry
Definition: adapt.hh:29
std::string temperature() noexcept
I/O name of temperature for equilibrium models.
Definition: name.hh:51
std::string pressure(int phaseIdx) noexcept
I/O name of pressure for multiphase systems.
Definition: name.hh:34
std::string density(int phaseIdx) noexcept
I/O name of density for multiphase systems.
Definition: name.hh:65
Definition: components/base.hh:43
Definition: components/base.hh:47
Scalar Scalar
export the scalar type used by the component
Definition: components/base.hh:51
Interface for components that have a gas state.
Definition: gas.hh:41
Interface for components that have a liquid state.
Definition: liquid.hh:41
A much simpler (and thus potentially less buggy) version of pure water.
Definition: simpleh2o.hh:51
static Scalar criticalTemperature()
Returns the critical temperature of water.
Definition: simpleh2o.hh:72
static constexpr bool gasViscosityIsConstant()
Returns true if the gas phase viscostiy is constant.
Definition: simpleh2o.hh:204
static std::string name()
A human readable name for the water.
Definition: simpleh2o.hh:60
static Scalar criticalPressure()
Returns the critical pressure of water.
Definition: simpleh2o.hh:78
static Scalar tripleTemperature()
Returns the temperature at water's triple point.
Definition: simpleh2o.hh:84
static Scalar liquidMolarDensity(Scalar temperature, Scalar pressure)
The molar density of pure water in at a given pressure and temperature.
Definition: simpleh2o.hh:271
static constexpr bool liquidViscosityIsConstant()
Returns true if the liquid phase viscostiy is constant.
Definition: simpleh2o.hh:210
static constexpr Scalar molarMass()
The molar mass in of water.
Definition: simpleh2o.hh:66
static constexpr bool liquidIsCompressible()
Returns true if the liquid phase is assumed to be compressible.
Definition: simpleh2o.hh:198
static const Scalar gasInternalEnergy(Scalar temperature, Scalar pressure)
Specific internal energy of steam .
Definition: simpleh2o.hh:168
static Scalar liquidThermalConductivity(Scalar temperature, Scalar pressure)
Thermal conductivity of water. source: http://webbook.nist.gov/cgi/fluid.cgi?ID=C7732185&Action=Page...
Definition: simpleh2o.hh:328
static Scalar gasHeatCapacity(Scalar temperature, Scalar pressure)
Specific isobaric heat capacity of water steam . source: http://webbook.nist.gov/cgi/fluid....
Definition: simpleh2o.hh:352
static Scalar gasViscosity(Scalar temperature, Scalar pressure, bool regularize=true)
The dynamic viscosity of steam.
Definition: simpleh2o.hh:293
static const Scalar liquidEnthalpy(Scalar temperature, Scalar pressure)
Specific enthalpy of liquid water .
Definition: simpleh2o.hh:149
static Scalar liquidPressure(Scalar temperature, Scalar density)
The pressure of water in at a given density and temperature.
Definition: simpleh2o.hh:280
static Scalar liquidHeatCapacity(Scalar temperature, Scalar pressure)
Specific isobaric heat capacity of the component as a liquid. source: http://webbook....
Definition: simpleh2o.hh:316
static Scalar gasPressure(Scalar temperature, Scalar density)
The pressure of steam in at a given density and temperature.
Definition: simpleh2o.hh:247
static Scalar liquidViscosity(Scalar temperature, Scalar pressure)
The dynamic viscosity of pure water.
Definition: simpleh2o.hh:304
static const Scalar gasEnthalpy(Scalar temperature, Scalar pressure)
Specific enthalpy of water steam .
Definition: simpleh2o.hh:139
static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
Thermal conductivity of steam. source: http://webbook.nist.gov/cgi/fluid.cgi?ID=C7732185&Action=Page...
Definition: simpleh2o.hh:340
static Scalar gasMolarDensity(Scalar temperature, Scalar pressure)
The molar density of steam in at a given pressure and temperature.
Definition: simpleh2o.hh:232
static Scalar triplePressure()
Returns the pressure at water's triple point.
Definition: simpleh2o.hh:90
static constexpr bool gasIsCompressible()
Returns true if the gas phase is assumed to be compressible.
Definition: simpleh2o.hh:192
static Scalar vaporPressure(Scalar T)
The vapor pressure in of pure water at a given temperature.
Definition: simpleh2o.hh:105
static const Scalar liquidInternalEnergy(Scalar temperature, Scalar pressure)
Specific internal energy of liquid water .
Definition: simpleh2o.hh:183
static Scalar liquidDensity(Scalar temperature, Scalar pressure)
The density of pure water at a given pressure and temperature .
Definition: simpleh2o.hh:259
static constexpr bool gasIsIdeal()
Returns true if the gas phase is assumed to be ideal.
Definition: simpleh2o.hh:238
static Scalar gasDensity(Scalar temperature, Scalar pressure)
The density of steam at a given pressure and temperature.
Definition: simpleh2o.hh:219
A central place for various physical constants occuring in some equations.
Definition: constants.hh:39
Relations valid for an ideal gas.
Definition: idealgas.hh:37
static constexpr Scalar pressure(Scalar temperature, Scalar rhoMolar)
The pressure of the gas in , depending on the molar density and temperature.
Definition: idealgas.hh:60
static constexpr Scalar R
The ideal gas constant .
Definition: idealgas.hh:40
static constexpr Scalar density(Scalar avgMolarMass, Scalar temperature, Scalar pressure)
The density of the gas in , depending on pressure, temperature and average molar mass of the gas.
Definition: idealgas.hh:49
static constexpr Scalar molarDensity(Scalar temperature, Scalar pressure)
The molar density of the gas , depending on pressure and temperature.
Definition: idealgas.hh:70
Base class for all components Components provide the thermodynamic relations for the liquid,...