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

*To*: hui Tan <hut2+@xxxxxxxx>*Subject*: Re: kinetic problem*From*: "David L Parkhurst" <dlpark@xxxxxxxx>*Date*: Wed, 28 May 2003 09:02:12 -0600*Bcc*: "David L Parkhurst" <dlpark@xxxxxxxx>*In-reply-to*: <439551406.1054059751@Yucca.geology.pitt.edu>

Long time no talk. Hope everything is going right for you. > My student Hui Tan and I have encountered the problem when using kinetics with phreeqc. After reaching eqbm after about 0.1 year, it oscillates between zero and negative rates --- meaning feldspar is precipitating. I suspect that this is a numerical problem. Do you have any ideas how we can fix or trick it? Specify -tol 1e-12 in the KINETICS definition. This tolerance reduces the oscillations to the 1e-15 range. By default PHREEQC integrates the equations to a precision of 1e-8 moles. In this case, you need a smaller tolerance. > We have math profs here who are expert. If you describe clear enough what is the problem, they may be able to help devise some numerical solutions. It would be advantageous to use an implicit method for integrating the rate equations. I worked on this once with an implicit Runge-Kutta method, but did not finish. This approach would help with stiff equations where the rates differ by orders of magnitude, like when you want one reaction to be nearly equilibrium and others slow. You may run into problems with stiff equations with your set of rates, I'm not sure. 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

**Follow-Ups**:**Re: kinetic problem***From:*hui Tan

**References**:**kinetic problem***From:*hui Tan

- Prev by Date:
**Re: PHREEQC problem** - Next by Date:
**Question about PHAST** - Previous by thread:
**kinetic problem** - Next by thread:
**Re: kinetic problem** - Index(es):

Project Home Page

Complete Water Resources Division Software

USGS Home Page

Water Resources Division Home Page

NRP Home Page

Help Page

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/msg00646.html

Email:dlpark@usgs.gov

Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $

Visitor number 1502 since Jan 22, 1998.