EOS5#

The EOS5 fluid property module is developed to study the behavior of groundwater systems in which hydrogen releases take place. For instance, in a number of waste disposal projects, corrosive metals are to be emplaced in geologic formations beneath the water table. These will evolve a mixture of gases, with hydrogen being the chief constituent. EOS5 is a close “cousin” of EOS3, the main difference being that the air component is replaced by hydrogen, with considerably different thermophysical properties (see Table 71 to Table 73). The assignment and handling of primary thermodynamic variables in EOS5 is identical to EOS3 (see Table 74). The main differences in the assignment of secondary parameters are as follows. Density of gaseous hydrogen is computed from the ideal gas law. Viscosity and water solubility of hydrogen are interpolated from the data given in Table 1. For temperatures in excess of 25˚C, the solubility at 25˚C is used.

Table 71 Density of hydrogen at \(P\) = 1 bar.#

Temperature

Experimental*

Ideal gas law†

280 K

0.086546 kg/m3

0.08660 kg/m3

300 K

0.080776 kg/m3

0.08082 kg/m3

Table 72 Viscosity* of hydrogen.#

Pressure

\(T\) = 0˚C

\(T\) = 100˚C

1 bar

8.40 x 10-6 Pa·s

10.33 x 10-6 Pa·s

100 bar

8.57 x 10-6 Pa·s

10.44 x 10-6 Pa·s

Table 73 Solubility in water at \(P\) = 1 bar§.#

Temperature

Solubility

\(T\) = 0˚C

1.92 x 10-6 g H2/g H2O

\(T\) = 25˚C

1.54 x 10-6 g H2/g H2O

Note

* from Vargaftik, 1975, p. 39.
† universal gas constant \(R\) = 8314.56 J/mol/˚C; molecular weight of hydrogen 2.0160 kg/mol.
§ after Dean, 1999; solubility at different pressures is computed from Henry’s law.

Specifications#

A summary of EOS5 specifications and parameters is given in Table Table 74. The default parameter settings are (NK, NEQ, NPH, NB) = (2, 3, 2, 6). The option NEQ = 2 is available for constant temperature conditions. The choice of primary thermodynamic variables is (\(P\), \(X\), \(T\)) for single-phase, (\(P_g\), \(S_g\) + 10, \(T\)) for two-phase conditions. As a convenience to users, we retain the capability to optionally initialize flow problems with TOUGH-style primary variables by setting MOP(19) = 1. In TOUGH we have (\(P\), \(T\), \(X\)) for single-phase conditions, (\(P_g\), \(S_g\), \(T\)) for two-phase conditions.

Table 74 Summary of EOS5.#

Components

#1: water
#2: hydrogen

Parameter choices

(NK, NEQ, NPH, NB) =
(2, 3, 2, 6) water and hydrogen, nonisothermal (default)
(2, 2, 2, 6) water and hydrogen, isothermal
Molecular diffusion can be modeled by setting NB = 8

Primary variables*

Single-phase conditions:
(\(P\), \(X\), \(T\)): (pressure, hydrogen mass fraction, temperature)
Two-phase conditions:
(\(P_g\), \(S_g\) + 10, \(T\)): (gas phase pressure, gas saturation plus 10, temperature)

Note

* By setting MOP(19) = 1, initialization can be made with TOUGH-style variables (\(P\), \(T\), \(X\)) for single-phase, (\(P_g\), \(S_g\), \(T\)) for two-phase.