27 #ifndef OPM_IAPWS_COMMON_HPP
28 #define OPM_IAPWS_COMMON_HPP
53 template <
class Scalar>
61 static const Scalar
Rs;
98 template <
class Evaluation>
99 static Evaluation
viscosity(
const Evaluation& temperature,
const Evaluation& rho)
103 Evaluation rhoBar = rho/322.0;
107 const Scalar Hij[6][7] = {
108 { 5.20094e-1, 2.22531e-1,-2.81378e-1, 1.61913e-1,-3.25372e-2, 0, 0 },
109 { 8.50895e-2, 9.99115e-1,-9.06851e-1, 2.57399e-1, 0, 0, 0 },
110 { -1.08374, 1.88797 ,-7.72479e-1, 0, 0, 0, 0 },
111 { -2.89555e-1, 1.26613 ,-4.89837e-1, 0, 6.98452e-2, 0,-4.35673e-3 },
112 { 0, 0,-2.57040e-1, 0, 0, 8.72102e-3, 0 },
113 { 0, 1.20573e-1, 0, 0, 0, 0,-5.93264e-4 }
116 Evaluation tmp, tmp2, tmp3 = 1;
117 Evaluation muBar = 0;
118 for (
int i = 0; i <= 5; ++i) {
121 for (
int j = 0; j <= 6; ++j) {
122 tmp += Hij[i][j]*tmp2;
123 tmp2 *= (rhoBar - 1);
126 tmp3 *= 1.0/TBar - 1;
133 const Scalar H[4] = {
134 1.67752, 2.20462, 0.6366564, -0.241605
138 for (
int i = 0; i < 4; ++i) {
160 template <
class Evaluation>
165 static const Scalar thcond_tstar = 647.26 ;
166 static const Scalar thcond_rhostar = 317.7 ;
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 ;
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 ;
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] = {
194 Evaluation Tbar = T / thcond_tstar;
195 Evaluation rhobar = rho / thcond_rhostar;
199 Evaluation Tpow = Troot;
202 for(
int k = 0; k < thcond_a_count; ++k) {
203 lam += thcond_a[k] * Tpow;
208 thcond_b0 + thcond_b1
210 *
Toolbox::exp(thcond_B1 * ((rhobar + thcond_B2)*(rhobar + thcond_B2)));
214 Evaluation Q = 2. + thcond_c5 / DTbarpow;
220 S = thcond_c6 / DTbarpow;
226 (thcond_d1 /
Toolbox::pow(Tbar,10.0) + thcond_d2) * rhobar18 *
227 Toolbox::exp(thcond_c1 * (1 - rhobar * rhobar18))
228 + thcond_d3 * S * rhobarQ *
229 Toolbox::exp((Q/(1+Q))*(1 - rhobar*rhobarQ))
236 template <
class Scalar>
237 const Scalar Common<Scalar>::molarMass = 18.01518e-3;
238 template <
class Scalar>
240 template <
class Scalar>
241 const Scalar Common<Scalar>::criticalTemperature = 647.096;
242 template <
class Scalar>
243 const Scalar Common<Scalar>::criticalPressure = 22.064e6;
244 template <
class Scalar>
245 const Scalar Common<Scalar>::criticalDensity = 322.0;
246 template <
class Scalar>
247 const Scalar Common<Scalar>::criticalMolarVolume = molarMass/criticalDensity;
248 template <
class Scalar>
249 const Scalar Common<Scalar>::acentricFactor = 0.344;
250 template <
class Scalar>
251 const Scalar Common<Scalar>::tripleTemperature = 273.16;
252 template <
class Scalar>
253 const Scalar Common<Scalar>::triplePressure = 611.657;
Definition: Air_Mesitylene.hpp:31
static Evaluation thermalConductivityIAPWS(const Evaluation &T, const Evaluation &rho)
Thermal conductivity water (IAPWS) .
Definition: Common.hpp:161
Evaluation< Scalar, VarSetTag, numVars > sqrt(const Evaluation< Scalar, VarSetTag, numVars > &x)
Definition: Math.hpp:278
A central place for various physical constants occuring in some equations.
Evaluation< Scalar, VarSetTag, numVars > exp(const Evaluation< Scalar, VarSetTag, numVars > &x)
Definition: Math.hpp:295
static const Scalar criticalDensity
Density of water at the critical point .
Definition: Common.hpp:70
static const Scalar criticalMolarVolume
Critical molar volume of water .
Definition: Common.hpp:73
Evaluation< Scalar, VarSetTag, numVars > abs(const Evaluation< Scalar, VarSetTag, numVars > &)
Definition: Math.hpp:41
static const Scalar Rs
Specific gas constant of water .
Definition: Common.hpp:61
static const Scalar tripleTemperature
Triple temperature of water .
Definition: Common.hpp:79
Implements relations which are common for all regions of the IAPWS '97 formulation.
Definition: Common.hpp:54
Evaluation< Scalar, VarSetTag, numVars > pow(const Evaluation< Scalar, VarSetTag, numVars > &base, Scalar exp)
Definition: Math.hpp:312
static const Scalar molarMass
The molar mass of water .
Definition: Common.hpp:58
static const Scalar criticalPressure
Critical pressure of water .
Definition: Common.hpp:67
static const Scalar criticalTemperature
Critical temperature of water .
Definition: Common.hpp:64
static const Scalar triplePressure
Triple pressure of water .
Definition: Common.hpp:82
static Evaluation viscosity(const Evaluation &temperature, const Evaluation &rho)
The dynamic viscosity of pure water.
Definition: Common.hpp:99
A central place for various physical constants occuring in some equations.
Definition: Constants.hpp:39
static const Scalar acentricFactor
The acentric factor of water .
Definition: Common.hpp:76