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 »
  • Beginners »
  • PHREEQC Basics »
  • SOLUTION SPECIES problem
« previous next »
  • Print
Pages: [1]   Go Down

Author Topic: SOLUTION SPECIES problem  (Read 485 times)

EMH121

  • Frequent Contributor
  • Posts: 24
SOLUTION SPECIES problem
« on: January 31, 2017, 08:26:57 PM »
When I modified the LLNL DB by adding the following, and running for 25oC

SOLUTION_SPREAD
    -pH       6.0
  Number       pe     S
      1        2    28
      2        1    28
      3        0    28
      4       -1    28
      5       -2    28
      6       -3    28
      7    -3.66    28
      8       -4    28
      9       -5    28
     10       -6    28
     11       -8    28
     
   
 SOLUTION_SPECIES

SO4-2 + H+ = HS- + 2O2
   -log_k   -132.56

output was reasonable and different from an identical LLNL DB (log_k = -138.32) run that did not have the SOLUTION_SPECIES keyword block.

When I added the T-informing "analytical" line as well, as below, the output warned "max iterations exceeded 100" and the output did not include any real sulfate data. The analytical expression constants are of the prescribed form a + b*x + c/x +d*log10(x) + e/(x^2) + f*(x^2)

SO4-2 + H+ = HS- + 2O2
   -log_k   -132.56
   -analytical   1.71E+04   4.83E+00   -7.86E+05   -1.24E+03   3.53E+07   -1.49E-03

The analytical expression returns log_K = -132.54 when calculated in Excel. Can someone help please? Thanks in advance for your help!

Logged

dlparkhurst

  • Top Contributor
  • Posts: 1034
Re: SOLUTION SPECIES problem
« Reply #1 on: February 01, 2017, 05:15:33 AM »
The following produces the output:

HS- Log K:    1.3100e+04

So, I think your analytical expression is not correct.

SOLUTION_SPECIES
SO4-2 + H+ = HS- + 2O2
   -log_k   -132.56
   -analytical   1.71E+04   4.83E+00   -7.86E+05   -1.24E+03   3.53E+07   -1.49E-03
USER_PRINT
-start
10 PRINT "HS- Log K: ", LK_SPECIES("HS-")
-end
SOLUTION
END
Logged

EMH121

  • Frequent Contributor
  • Posts: 24
Re: SOLUTION SPECIES problem
« Reply #2 on: February 01, 2017, 06:59:29 PM »
Thank you Dr Parkhurst! My error was in the d*log10(T) term. The correct d constant in this case is -6583.929939. The keyword block below worked without error.

 SOLUTION_SPECIES

SO4-2 + H+ = HS- + 2O2
   -log_k   -132.56
   -analytical   1.71E+04   4.83E+00   -7.86E+05   -6583.929939   3.53E+07   -1.49E-03

Logged

MichaelZ20

  • Frequent Contributor
  • Posts: 16
Re: SOLUTION SPECIES problem
« Reply #3 on: February 02, 2017, 08:02:12 AM »
Calculated log_k equals -123.04
Logged

  • Print
Pages: [1]   Go Up
« previous next »
  • PhreeqcUsers Discussion Forum »
  • Beginners »
  • PHREEQC Basics »
  • SOLUTION SPECIES problem
 

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