PhreeqcUsers Discussion Forum
Click here to donate to keep PhreeqcUsers open

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 »
  • Beginners »
  • PHREEQC basics »
  • Evaporate a solution using PHREEQC
« previous next »
  • Print
Pages: [1]   Go Down

Author Topic: Evaporate a solution using PHREEQC  (Read 6323 times)

Tom

  • Admin
  • Posts: 123
Evaporate a solution using PHREEQC
« on: April 26, 2014, 10:56:41 AM »
Q. How can I simulate evaporation of a solution in PHREEQC? I am looking to investigate the minerals which are likely to precipitate.

A. Evaporation can be simulated using the REACTION keyword:

In the example below, 1 kg of a hypothetical solution (100 %) is evaporated.

----------------
SOLUTION 1
pH 7 charge #adjust pH in accordance with charge balance
units ppm
S(6) 1700
Mg 50
Ca 200
Fe 10
Si 10
water 1.0 #1 kg of water

REACTION 1
H2O -1.0 # -1.0 is the relative moles of the reaction.
55.51 moles in 10 steps #1 kg of water is 55.51 moles. 5.551 moles H2O removed per step (=100 ml H2O per step)
END

--------------

From the output you can assess the SIs of phases at each step.
« Last Edit: April 27, 2014, 09:11:22 AM by Charlie »
Logged

Charlie

  • Admin
  • Posts: 159
Re: Evaporate a solution using PHREEQC
« Reply #1 on: April 26, 2014, 05:08:45 PM »

Adding this block will help assess the reaction progress:

SELECTED_OUTPUT
-reset                false
-step         true
-water               true
Logged

dlparkhurst

  • Top Contributor
  • Posts: 2736
Re: Evaporate a solution using PHREEQC
« Reply #2 on: April 26, 2014, 06:04:01 PM »
For high evaporation factors, PHREEQC can fail to find the amount of water left in solution, and then the numerical method fails. Another approach that can be useful is to fix the activity of water with a pseudo phase as in the following script.

SOLUTION
   Na   1
   Cl   1
END
PHASES
water
    H2O = H2O
    log_k     0
END
USE solution 1
EQUILIBRIUM_PHASES 1
    water     -0.1 10
END
Logged

liamabrahamsen

  • Contributor
  • Posts: 9
Re: Evaporate a solution using PHREEQC
« Reply #3 on: March 10, 2015, 11:36:49 AM »
That's an interesting approach. What effect does introducing the pseudo water phase have? Doesn't it simply replace the water that was removed in the evaporation step?
Logged

dlparkhurst

  • Top Contributor
  • Posts: 2736
Re: Evaporate a solution using PHREEQC
« Reply #4 on: March 10, 2015, 02:26:16 PM »
You could have the same effect using the H2O(g) phase, but the target saturation index would be a bit harder to set up. With H2O(g), the target saturation index is equal to the log of the partial pressure of water that is desired. Whereas, the "water" phase simply removes (or adds) water until the log activity of water is equal to the target saturation index.
« Last Edit: February 04, 2021, 03:03:05 PM by dlparkhurst »
Logged

MichaelZ20

  • Top Contributor
  • Posts: 101
Re: Evaporate a solution using PHREEQC
« Reply #5 on: February 04, 2021, 10:27:34 AM »
Should we consider reliable the results of modeling strong evaporation when the ionic strength exceeds 6-10?
Logged

dlparkhurst

  • Top Contributor
  • Posts: 2736
Re: Evaporate a solution using PHREEQC
« Reply #6 on: February 04, 2021, 02:54:22 PM »
The caveats related to high ionic strength always apply. Pitzer.dat would be best, perhaps followed by sit.dat, followed by the other ion association models, in order of reliability.

The ionic composition of the solution makes a difference. Pitzer parameters are fit over a range of concentration, and some elements may have been fit only to a lower ionic strength. Ion association models often rely on activity coefficients fit from chloride and sodium salts, so they will work best in a NaCl medium.
Logged

  • Print
Pages: [1]   Go Up
« previous next »
  • PhreeqcUsers Discussion Forum »
  • Beginners »
  • PHREEQC basics »
  • Evaporate a solution using PHREEQC
 

  • SMF 2.0.17 | SMF © 2019, Simple Machines | Terms and Policies
  • XHTML
  • RSS
  • WAP2