Problem 1 - Brine Density Calculation (dnh)

This problem is from the sample problem No. 1 in TOUGH3/EWASG user manual.

The problem demonstrates the brine correlations. Several element subproblems are simulated, which are entirely independent of each other (no flow connections between subproblems), except that being run together they all must go through the same sequence of time steps. Only a single time step (1E-8 sec) is computed to calculate the brine density at different initial salt mass fraction. The salt content approximately varies from 1 to 6 molar. The initial condition for these elements are presented in Figure 10-10.

Figure 10-10 Initial conditions representing different salt content for the 6 elements.

Brine properties can be calculated with different approaches through the options specified by IE array. Discussion of these IE options are present in Section 7/EWASG. Data block OUTPU is used to select the liquid density as a single print-out variable. Two simulation cases are run:

(1) IE(4)=6, using Driesner (2007) for brine density calculation.

(2) IE(4)=1 and IE(15)=4, using Andersen et al. (1992) for density calculation, using Lorenz et al. (2000) for brine enthalpy calculation.

Note that if only either the brine density (IE(4)) or enthalpy (IE(15)) is specified to use the Driesner’s correlation, TOUGH4 internally assigns the other brine property uses the Driesner’s correlation as well. Table 10-2 shows the results of calculated brine density for the both cases. The brine densities are slightly overestimated with IE(4) = 1, compared to those by Driesner (2007).

Table 10-2 Simulated brine densities for case 1 and case 2

Element	case 1 (kg/ m^3)	case 2 (kg/m^3)
F 1	0.85027E+03	0.86241E+03
F 2	0.89252E+03	0.89902E+03
F 3	0.93441E+03	0.94535E+03
F 4	0.97668E+03	0.99803E+03
F 5	0.10198E+04	0.10571E+04
F 6	0.10642E+04	0.11286E+04

Input Files: INFILE_case1, INFILE_case2

Output Files: dnh_results_case1.zip, dnh_results_case2.zip