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
»
How to get -totals for SOLUTION_MODIFY?
« previous
next »
Print
Pages: [
1
]
Go Down
Author
Topic: How to get -totals for SOLUTION_MODIFY? (Read 992 times)
dlee05291
Contributor
Posts: 4
How to get -totals for SOLUTION_MODIFY?
«
on:
June 06, 2014, 05:53:30 PM »
Hi all,
I'm trying to couple COMSOL with PHREEQC. An issue that I have is getting -totals in SOLUTION_MODIFY. In my problem, transport is calculated using COMSOL, and chemical reaction and equilibrium is calculated using PHREEQC. COMSOL uses molarity so that I converted that into molality, and then the difference between current and previous molality is obtained. The mass_water calculated from the previous solution is multiplied to the difference in molality to get the change in the mass of solutes. Then the mass is subtracted from the current mass_water to get new mass_water. Finally, the current molality is multiplied by the new mass_water to obtain -totals for SOLUTION_MODIFY. This is o.k. as far as the flux is low enough. But when the flux become large, significant deviation is observed. Could anyone tell me whether my method is right? If not what is the correct way to get -totals for SOLUTION_MODIFY?
Thank you for your time!
Dooyoul Lee
Logged
dlparkhurst
Top Contributor
Posts: 1034
Re: How to get -totals for SOLUTION_MODIFY?
«
Reply #1 on:
June 07, 2014, 03:12:30 AM »
There will always be a discrepancy generated by the operator splitting. However, for the first step, I would take the molarity from COMSOL, and multiply by the volume of water (porosity * saturation) to get the moles of elements in the solution.
If there is no reaction, then you should be able to recreate the molarity by using the number of moles in solution (TOTMOL) divided by the solution volume. With some reaction, you will get some deviations that you have to decide how to deal with. Some strategies iterate to account for the differences, some simply ignore it.
Logged
dlee05291
Contributor
Posts: 4
Re: How to get -totals for SOLUTION_MODIFY?
«
Reply #2 on:
June 07, 2014, 07:59:39 PM »
Dear David,
Thank you for your quick reply. What you have mentioned is much simpler than I did. Thank you so much!
Dooyoul
Logged
Print
Pages: [
1
]
Go Up
« previous
next »
PhreeqcUsers Discussion Forum
»
Reactive Transport
»
Reactive Transport Modelling
»
How to get -totals for SOLUTION_MODIFY?