Attempting Two Column Separate Transfers with additional compexity

[laurafie]

Hello,

I am attempting to create a two column model with transfers occurring separately along each column and mixing at the end into the final solution.  I have attached a depiction of my concept and will describe below my explanation of my issues. Following my description below is the script I'm working off of.

Issues:
1.   Key issue-Difficulties with separate transfers
                       Is it even possible to do transferes on separate columns in the same model?
                       If not is it possible to connect the results so that I can mix them simply?
                       Is there a way to define the solution set used for transfer?
When I’ve been working separately with solution 80-81 column,
the transfer keyword came up with errors saying it couldn’t find solution zero, or solutions 1-3.
Is there any way to tell it to look for other solution names?

2.   Side note- Column 81-83 is messy and it is difficult to pinpoint where problems are occurring
                         I am looking at using a combination of SAVE, END, and use to break it up.
                         Is there a better way of getting outputs specifically for the separate parts.

Description:

There is a well that is getting a mixture of river water inflow and aquifer inflow:
final water… mix of solution 40 and 83.

The river water is coming in along around a two year path:
River Water Column Solutions 0-40
(These specifics have been worked out satisfactorily in previous modeling work).

The aquifer water column has evidence of salt rich water leaking in:
Artificial input of Solution 60 mixing with solution 81.

The aquifer water  column also has evidence of cation exchange occurring somewhere else:
exchange occurring in solution 82. 

I need transfer to occur along the path combining the leaking salt solution and cation exchange before mixing with incoming river water.

#Flow from Fox River to Well RL255
SOLUTION 0 infilling Isco& Fox2(2014 April-Dec., pH from Fox 2 4/07-10/14)
    temp      10
    pH        7.86
    pe        8.4
    redox     pe
    units     mmol/kgw
    density   1
    Alkalinity 4.59
    Ca        2.03
    Cl        7.43
    K         0.12
    Mg        1.58
    Na        6.22
    S         0.63
    -water    1 # kg
SOLUTION 1-40 Inital RL255(1April2005)
    temp      10
    pH        7.8
    pe        8.4
    redox     pe
    units     mmol/kgw
    density   1
    Alkalinity 6.8
    Ca        2.99
    Cl        3.38
    K         0.01
    Mg        2.3
    Na        1.35
    S         0.8
    -water    1 # kg
TRANSPORT
    -cells                 40
    -shifts                160
    -lengths               14*7.62 7.6 25*7.62
    -dispersivities        40*2
    -stagnant              1  0  0  0
    -print_cells           40
    -print_frequency       20
    -punch_cells           40
    -punch_frequency       20
#Flow through aquifer with CEC to Well RL 255
SOLUTION 81-83 Pre-aquifer Water-WK947(1April2005)Equalized Na/Cl
    temp      10
    pH        7.11
    pe        8.4
    redox     pe
    units     umol/kgw
    density   1
    Alkalinity 6.46
    Ca        2.4
    Cl        2.32
    Mg        2.11
    Na        2.32
    S         1.07
    -water    1 # kg
SOLUTION 60 Artifical NaCl input
    temp      10
    pH        7.11
    pe        8.4
    redox     pe
    units     mmol/kgw
    density   1
    Alkalinity 6.46
    Ca        2.4
    Cl        2.57
    K         0.7
    Mg        2.11
    Na        2.57
    -water    1 # kg
MIX 1
    60    0.5
    81    0.5
TRANSPORT
    -cells                 3
    -lengths               3*7.62
    -dispersivities        3*2
    -print_cells           1-3
    -punch_cells           1-3
EXCHANGE 82 pre-aquifer cation exchange
    X     1
#Mixing in Well RL255
MIX 2 Mixing in RL255: River and Aquifer Water
    40    0.44
    83    0.56
EQUILIBRIUM_PHASES 1 Model A (6.8.15)
    Dolomite  2.793e-04 10
    Calcite   -5.639e-04 10
    CO2(g)    2.934e-03 10
SELECTED_OUTPUT 1
    -file                 C:\~.sel
    -solution             true
    -time                 true
    -step                 true
    -totals               Alkalinity  Ca  Cl  K  Mg  Na  S(6)
    -molalities           CaCO3  CO2  MgCO3
    -equilibrium_phases   Calcite  CO2(g)  Dolomite
    -saturation_indices   Calcite  CO2(g)  Dolomite
END

Thank you very much for your time and energy!

[dlparkhurst]

1.   Key issue-Difficulties with separate transfers
                       Is it even possible to do transferes on separate columns in the same model?

It is probably possible if you were Tony Appelo. You can use TRANSPORT and MIXes to define the flow configuration. I can't say that I know how to do it.


                       If not is it possible to connect the results so that I can mix them simply?

Not sure.

                       Is there a way to define the solution set used for transfer?

MIX 80 will define the solutions that are mixed to produce the solution in cell 80.

When I’ve been working separately with solution 80-81 column,
the transfer keyword came up with errors saying it couldn’t find solution zero, or solutions 1-3.

Columns for TRANSPORT are defined for cells 1-n. It is not possible to define a column that starts at 80. You would have to define a column 1-82, and then define MIXes that produce the connections that you want.

Is there any way to tell it to look for other solution names?

No. TRANSPORT always uses 1-n, with infilling either from solution 0 or solution n+1, depending on the flow direction.

2.   Side note- Column 81-83 is messy and it is difficult to pinpoint where problems are occurring
                         I am looking at using a combination of SAVE, END, and use to break it up.

You can use SAVE and END to define the initial conditions for the column and then use TRANSPORT. It is also possible to use MIX and RUN_CELLS, but it would require input for each time step.

                         Is there a better way of getting outputs specifically for the separate parts.

You can use multiple SELECTED_OUTPUT n and USER_PUNCH n to define different output files for different types of output. You can use Basic IF statements with cell numbers (cell_no) to subset the cells into different files, for example.

[laurafie]

Thank you very much for your help.

It doesn't seem like any of the methods you mentioned would allow for the Aquifer column to be modeled over time in sync with the river column.  I think you are correct, I might just need to be Tony Appelo to model this situation.