Using REACTION to simulate water evaporation from a group (SOLUTION_SPREAD)

[peterwadeuk]

Greetings, fellow GeoModellers.

I have a problem with simulating evaporation on a group of Solutions, which I define with SOLUTION_SPREAD.

I attempt to remove water using REACTION, but I only receive information for REACTION steps involving the first SOLUTION.

I believe I need advice.

Is what I am attempting within PhreeqC feasible? If so, what is a better (functional) way of achieving my aim? Else I will simulate evaporation using a spreadsheet (and charge balance on pH).

Bon Voyage,
Peter

Here is my script:

TITLE Can-Wol All WQ constituents per L_effluent
TITLE RAW DATA

PHASES
Moganite
    SiO2 = SiO2
    log_k     -3.025

SOLUTION_SPREAD
    -units    mg/l
         Description    Temperature    Number       C       Ca      Cl        K      Mg    N(5)      Na       P      S(6)      Si      pH
                                                                                                                  charge                
  GCL0229-002-HCT-S2          20.55         2    10.3     55.7    17.3     4.06    4.66    39.8    15.0    0.13      38.0    69.6    7.72
  GCL0229-002-HCT-S3          20.68         3    14.9     29.6    13.5     8.05    5.13    9.32    27.1    0.36      15.2    32.2    7.58
  GCL0229-002-HCT-S4          20.61         4    12.9    24.15    0.50     1.55    3.56    0.50    2.51    0.14      4.20    62.4    8.13
  GCL0229-002-HCT-S5          20.58         5    12.8     24.6    0.50     1.32    3.81    0.50    1.89    0.12      2.92    58.7    8.62
  GCL0229-002-HCT-S6          20.60         6    13.8     28.2    0.50     1.27    4.15    0.50    1.86    0.15      7.90    66.5    8.44
  GCL0229-002-HCT-S7          20.66         7    15.0     30.5    0.50     1.20    3.93    0.50    1.79    0.13      8.63    67.5    8.55
  GCL0229-002-HCT-S8          23.72         8    13.3     32.0    0.50     1.17    3.69    0.50    1.71    0.09      7.38    69.0    8.33
  GCL0229-002-HCT-S9          19.81         9    13.6     32.0    0.50     1.04    3.10    0.50    1.54    0.09      6.82    69.9    8.22
 GCL0229-002-HCT-S10          21.39        10    13.2     32.3    0.50    0.919    2.73    0.50    1.44    0.10      8.27    70.1    8.13
 GCL0229-002-HCT-S11          21.27        11    13.7     33.8    0.50    0.863    2.51    0.50    1.44    0.08      8.30    79.2    8.19
 GCL0229-002-HCT-S12          18.43        12    19.4     31.8    0.50    0.444    1.17    0.50    1.02    0.13      4.44    81.0    8.42
 GCL0229-002-HCT-S13          19.95        13    15.0     33.4    0.50    0.755    1.96    0.50    1.33    0.12      9.18    79.5    8.23
 GCL0229-002-HCT-S14          20.55        14    14.5     32.6    0.50    0.683    1.74    0.50    1.24    0.09      7.80    76.9    8.30
 GCL0229-002-HCT-S15          20.88        15    14.5     32.7    0.50    0.652    1.55    0.50    1.21    0.10      8.95    75.5    8.21
 GCL0229-002-HCT-S16          19.80        16    13.7     34.8    0.50    0.643    1.54    0.50    1.24    0.09      12.8    83.1    7.85
 GCL0229-002-HCT-S17          20.49        17    15.0     40.7    0.50    0.670    1.57    0.50    1.29    0.13      20.9    92.3    8.30
 GCL0229-002-HCT-S18          21.51        18    14.4     39.7    0.50    0.591    1.32    0.50    1.12    0.11      32.0    92.0    8.50
 GCL0229-002-HCT-S19          20.19        19    14.8     31.0    0.50    0.897    3.01    0.50    3.26    0.55      38.7    62.6    8.46

SELECTED_OUTPUT 1
    -file                 modl-evap.tsv
    -water                true
    -percent_error        true
    -totals               C(4)  Ca  Cl  Fe  K  Mg  N(5)
                          Na  P  S(6)  Si
    -saturation_indices   Calcite  Moganite  Fe(OH)3  SiO2(am)
                          CO2(g)
#END
REACTION 2-19
H2O -1.0
54.951 moles in 5 steps # Evaporate 99% of solutions in 5 steps
END

[dlparkhurst]

You can do it with RUN_CELLS.

Code: [Select]

TITLE RAW DATA
SOLUTION_SPREAD
    -units    mg/l
         Description    Temperature    Number       C       Ca      Cl        K      Mg    N(5)      Na       P      S(6)      Si      pH
                                                                                                                  charge               
  GCL0229-002-HCT-S2          20.55         2    10.3     55.7    17.3     4.06    4.66    39.8    15.0    0.13      38.0    69.6    7.72
  GCL0229-002-HCT-S3          20.68         3    14.9     29.6    13.5     8.05    5.13    9.32    27.1    0.36      15.2    32.2    7.58
  GCL0229-002-HCT-S4          20.61         4    12.9    24.15    0.50     1.55    3.56    0.50    2.51    0.14      4.20    62.4    8.13
  GCL0229-002-HCT-S5          20.58         5    12.8     24.6    0.50     1.32    3.81    0.50    1.89    0.12      2.92    58.7    8.62
  GCL0229-002-HCT-S6          20.60         6    13.8     28.2    0.50     1.27    4.15    0.50    1.86    0.15      7.90    66.5    8.44
  GCL0229-002-HCT-S7          20.66         7    15.0     30.5    0.50     1.20    3.93    0.50    1.79    0.13      8.63    67.5    8.55
  GCL0229-002-HCT-S8          23.72         8    13.3     32.0    0.50     1.17    3.69    0.50    1.71    0.09      7.38    69.0    8.33
  GCL0229-002-HCT-S9          19.81         9    13.6     32.0    0.50     1.04    3.10    0.50    1.54    0.09      6.82    69.9    8.22
 GCL0229-002-HCT-S10          21.39        10    13.2     32.3    0.50    0.919    2.73    0.50    1.44    0.10      8.27    70.1    8.13
 GCL0229-002-HCT-S11          21.27        11    13.7     33.8    0.50    0.863    2.51    0.50    1.44    0.08      8.30    79.2    8.19
 GCL0229-002-HCT-S12          18.43        12    19.4     31.8    0.50    0.444    1.17    0.50    1.02    0.13      4.44    81.0    8.42
 GCL0229-002-HCT-S13          19.95        13    15.0     33.4    0.50    0.755    1.96    0.50    1.33    0.12      9.18    79.5    8.23
 GCL0229-002-HCT-S14          20.55        14    14.5     32.6    0.50    0.683    1.74    0.50    1.24    0.09      7.80    76.9    8.30
 GCL0229-002-HCT-S15          20.88        15    14.5     32.7    0.50    0.652    1.55    0.50    1.21    0.10      8.95    75.5    8.21
 GCL0229-002-HCT-S16          19.80        16    13.7     34.8    0.50    0.643    1.54    0.50    1.24    0.09      12.8    83.1    7.85
 GCL0229-002-HCT-S17          20.49        17    15.0     40.7    0.50    0.670    1.57    0.50    1.29    0.13      20.9    92.3    8.30
 GCL0229-002-HCT-S18          21.51        18    14.4     39.7    0.50    0.591    1.32    0.50    1.12    0.11      32.0    92.0    8.50
 GCL0229-002-HCT-S19          20.19        19    14.8     31.0    0.50    0.897    3.01    0.50    3.26    0.55      38.7    62.6    8.46
END

REACTION 2-19
H2O -1.0
54.951 moles in 5 steps # Evaporate 99% of solutions in 5 steps
END

RUN_CELLS
-cells 2-19
END


[peterwadeuk]

Thank you very much, kind Sir :)