[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: convergence problem




> As the problem occurs during a transport
simulation it seems the best if I used
the last option you proposed.
When writing the input for the
reaction step I will check if the O(0) conc
is 0 and then set it to a small value ... or
do you think it would be a problem if I generally
avoided O(0) to be 0 ?


> As far as I remember I previously fixed a problem
by replacing very small concentrations (eg 1e-35)
with zeros in the input file(s).

PHREEQC ignores concentrations less than about 1e-25, at least in some of
the time. Concentrations less than these could cause problems, so what you
did (setting to zero) probably makes sense for elements other than H and O.
For the O2, I think you are best to set it, either to the correct
concentration if greater than 0, or to a value between 1e-15 and 1e-20 if
equal to zero. I think you will get some erroneous results if you are
generating dissolved O2 by the specification of pe and adjusting pH.
PHREEQC will calculate an O2 that corresponds (maybe millimoles to moles
given your example pH 9, pe 12) and this O2 will then be in the system,
which would then probably oxidize pyrite. However, the added O2 would be an
artifact of the calculation (spontaneous generation of O2) causing
mass-balance and other misleading results.

David


David Parkhurst (dlpark@xxxxxxxx)
U.S. Geological Survey
Box 25046, MS 413
Denver Federal Center
Denver, CO 80225

Project web page: http://wwwbrr.cr.usgs.gov/projects/GWC_coupled




Project Home Page
Complete Water Resources Division Software
USGS Home Page
Water Resources Division Home Page
NRP Home Page
Help Page
USGS Privacy Statement        Disclaimer

Please note that some U.S. Geological Survey (USGS) information accessed through this page may be preliminary in nature and presented prior to final review and approval by the Director of the USGS. This information is provided with the understanding that it is not guaranteed to be correct or complete and conclusions drawn from such information are the sole responsibility of the user.

Any use of trade, product, or firm names in this publication is for descriptive purposes only and does not imply endorsement by the U.S. Government. 


The URL of this page is: http://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/mail/msg00050.html
Email:dlpark@usgs.gov
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 1042 since Jan 22, 1998.