Welcome,
Guest
. Please
login
or
register
.
Did you miss your
activation email
?
1 Hour
1 Day
1 Week
1 Month
Forever
Login with username, password and session length
Forum Home
Login
Register
PhreeqcUsers Discussion Forum
»
Reactive Transport
»
Reactive Transport Modelling
»
Attempting Two Column Separate Transfers with additional compexity
« previous
next »
Print
Pages: [
1
]
Go Down
Author
Topic: Attempting Two Column Separate Transfers with additional compexity (Read 900 times)
laurafie
Contributor
Posts: 2
Attempting Two Column Separate Transfers with additional compexity
«
on:
June 14, 2015, 11:41:18 PM »
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!
Logged
dlparkhurst
Top Contributor
Posts: 1033
Re: Attempting Two Column Separate Transfers with additional compexity
«
Reply #1 on:
June 16, 2015, 02:32:24 PM »
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.
Logged
laurafie
Contributor
Posts: 2
Re: Attempting Two Column Separate Transfers with additional compexity
«
Reply #2 on:
June 23, 2015, 07:10:27 PM »
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.
Logged
Print
Pages: [
1
]
Go Up
« previous
next »
PhreeqcUsers Discussion Forum
»
Reactive Transport
»
Reactive Transport Modelling
»
Attempting Two Column Separate Transfers with additional compexity