PhreeqcUsers Discussion Forum

Welcome, Guest. Please login or register.
Did you miss your activation email?

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?
 

  • SMF 2.0.3 | SMF © 2011, Simple Machines
  • XHTML
  • RSS
  • WAP2