My Project
Loading...
Searching...
No Matches
Common.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_IAPWS_COMMON_HPP
28#define OPM_IAPWS_COMMON_HPP
29
32#include <opm/common/utility/gpuDecorators.hpp>
33
34#include <cmath>
35
36namespace Opm {
37namespace IAPWS {
38
54template <class Scalar>
55class Common
56{
57public:
59 static const Scalar molarMass;
60
62 static const Scalar Rs;
63
65 static const Scalar criticalTemperature;
66
68 static const Scalar criticalPressure;
69
71 static const Scalar criticalDensity;
72
74 static const Scalar criticalVolume;
75
77 static const Scalar criticalMolarVolume;
78
80 static const Scalar acentricFactor;
81
83 static const Scalar tripleTemperature;
84
86 static const Scalar triplePressure;
87
102 template <class Evaluation>
103 OPM_HOST_DEVICE static Evaluation viscosity(const Evaluation& temperature, const Evaluation& rho)
104 {
105 Evaluation rhoBar = rho/322.0;
106 Evaluation TBar = temperature/criticalTemperature;
107
108 // muBar = muBar_1
109 const Scalar Hij[6][7] = {
110 { 5.20094e-1, 2.22531e-1,-2.81378e-1, 1.61913e-1,-3.25372e-2, 0, 0 },
111 { 8.50895e-2, 9.99115e-1,-9.06851e-1, 2.57399e-1, 0, 0, 0 },
112 { -1.08374, 1.88797 ,-7.72479e-1, 0, 0, 0, 0 },
113 { -2.89555e-1, 1.26613 ,-4.89837e-1, 0, 6.98452e-2, 0,-4.35673e-3 },
114 { 0, 0,-2.57040e-1, 0, 0, 8.72102e-3, 0 },
115 { 0, 1.20573e-1, 0, 0, 0, 0,-5.93264e-4 }
116 };
117
118 Evaluation tmp, tmp2, tmp3 = 1;
119 Evaluation muBar = 0;
120 for (int i = 0; i <= 5; ++i) {
121 tmp = 0;
122 tmp2 = 1;
123 for (int j = 0; j <= 6; ++j) {
124 tmp += Hij[i][j]*tmp2;
125 tmp2 *= (rhoBar - 1);
126 };
127 muBar += tmp3 * tmp;
128 tmp3 *= 1.0/TBar - 1;
129 };
130 muBar *= rhoBar;
131 muBar = exp(muBar);
132
133 // muBar *= muBar_0
134 muBar *= 100*sqrt(TBar);
135 const Scalar H[4] = {
136 1.67752, 2.20462, 0.6366564, -0.241605
137 };
138
139 tmp = 0, tmp2 = 1;
140 for (int i = 0; i < 4; ++i) {
141 tmp += H[i]/tmp2;
142 tmp2 *= TBar;
143 };
144 muBar /= tmp;
145
146 return 1e-6*muBar;
147 }
148
162 template <class Evaluation>
163 OPM_HOST_DEVICE static Evaluation thermalConductivityIAPWS(const Evaluation& T, const Evaluation& rho)
164 {
165 static const Scalar thcond_tstar = 647.26 ;
166 static const Scalar thcond_rhostar = 317.7 ;
167 /*static const Scalar thcond_kstar = 1.0 ;*/
168
169 static const Scalar thcond_b0 = -0.397070 ;
170 static const Scalar thcond_b1 = 0.400302 ;
171 static const Scalar thcond_b2 = 1.060000 ;
172 static const Scalar thcond_B1 = -0.171587 ;
173 static const Scalar thcond_B2 = 2.392190 ;
174
175 static const Scalar thcond_c1 = 0.642857 ;
176 static const Scalar thcond_c2 = -4.11717 ;
177 static const Scalar thcond_c3 = -6.17937 ;
178 static const Scalar thcond_c4 = 0.00308976 ;
179 static const Scalar thcond_c5 = 0.0822994 ;
180 static const Scalar thcond_c6 = 10.0932 ;
181
182 static const Scalar thcond_d1 = 0.0701309 ;
183 static const Scalar thcond_d2 = 0.0118520 ;
184 static const Scalar thcond_d3 = 0.00169937 ;
185 static const Scalar thcond_d4 = -1.0200 ;
186 static const int thcond_a_count = 4;
187 static const Scalar thcond_a[thcond_a_count] = {
188 0.0102811
189 ,0.0299621
190 ,0.0156146
191 ,-0.00422464
192 };
193
194 Evaluation Tbar = T / thcond_tstar;
195 Evaluation rhobar = rho / thcond_rhostar;
196
197 /* fast implementation... minimised calls to 'pow' routine... */
198 Evaluation Troot = sqrt(Tbar);
199 Evaluation Tpow = Troot;
200 Evaluation lam = 0;
201
202 for(int k = 0; k < thcond_a_count; ++k) {
203 lam += thcond_a[k] * Tpow;
204 Tpow *= Tbar;
205 }
206
207 lam +=
208 thcond_b0 + thcond_b1
209 * rhobar + thcond_b2
210 * exp(thcond_B1 * ((rhobar + thcond_B2)*(rhobar + thcond_B2)));
211
212 Evaluation DTbar = abs(Tbar - 1) + thcond_c4;
213 Evaluation DTbarpow = pow(DTbar, 3./5);
214 Evaluation Q = 2. + thcond_c5 / DTbarpow;
215
216 Evaluation S;
217 if(Tbar >= 1)
218 S = 1. / DTbar;
219 else
220 S = thcond_c6 / DTbarpow;
221
222 Evaluation rhobar18 = pow(rhobar, 1.8);
223 Evaluation rhobarQ = pow(rhobar, Q);
224
225 lam +=
226 (thcond_d1 / pow(Tbar,10.0) + thcond_d2) * rhobar18 *
227 exp(thcond_c1 * (1 - rhobar * rhobar18))
228 + thcond_d3 * S * rhobarQ *
229 exp((Q/(1+Q))*(1 - rhobar*rhobarQ))
230 + thcond_d4 *
231 exp(thcond_c2 * pow(Troot,3.0) + thcond_c3 / pow(rhobar,5.0));
232 return /*thcond_kstar * */ lam;
233 }
234};
235
236template <class Scalar>
237const Scalar Common<Scalar>::molarMass = 18.01518e-3;
238template <class Scalar>
239const Scalar Common<Scalar>::Rs = Constants<Scalar>::R/molarMass;
240template <class Scalar>
241const Scalar Common<Scalar>::criticalTemperature = 647.096;
242template <class Scalar>
243const Scalar Common<Scalar>::criticalPressure = 22.064e6;
244template <class Scalar>
245const Scalar Common<Scalar>::criticalDensity = 322.0;
246template <class Scalar>
247const Scalar Common<Scalar>::criticalMolarVolume = molarMass/criticalDensity;
248template <class Scalar>
249const Scalar Common<Scalar>::criticalVolume = 5.595e-2;
250template <class Scalar>
251const Scalar Common<Scalar>::acentricFactor = 0.344;
252template <class Scalar>
253const Scalar Common<Scalar>::tripleTemperature = 273.16;
254template <class Scalar>
255const Scalar Common<Scalar>::triplePressure = 611.657;
256
257} // namespace IAPWS
258} // namespace Opm
259
260#endif
A traits class which provides basic mathematical functions for arbitrary scalar floating point values...
Definition Constants.hpp:40
Implements relations which are common for all regions of the IAPWS '97 formulation.
Definition Common.hpp:56
static OPM_HOST_DEVICE Evaluation viscosity(const Evaluation &temperature, const Evaluation &rho)
The dynamic viscosity of pure water.
Definition Common.hpp:103
static const Scalar criticalVolume
Critical volume of water .
Definition Common.hpp:74
static const Scalar criticalPressure
Critical pressure of water .
Definition Common.hpp:68
static const Scalar criticalDensity
Density of water at the critical point .
Definition Common.hpp:71
static const Scalar criticalMolarVolume
Critical molar volume of water .
Definition Common.hpp:77
static const Scalar criticalTemperature
Critical temperature of water .
Definition Common.hpp:65
static OPM_HOST_DEVICE Evaluation thermalConductivityIAPWS(const Evaluation &T, const Evaluation &rho)
Thermal conductivity water (IAPWS) .
Definition Common.hpp:163
static const Scalar tripleTemperature
Triple temperature of water .
Definition Common.hpp:83
static const Scalar triplePressure
Triple pressure of water .
Definition Common.hpp:86
static const Scalar molarMass
The molar mass of water .
Definition Common.hpp:59
static const Scalar acentricFactor
The acentric factor of water .
Definition Common.hpp:80
static const Scalar Rs
Specific gas constant of water .
Definition Common.hpp:62
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30