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

Re: PHREEQC problem




I have been working on a method for integrating the rate equations that
works better with stiff equations. I have implemented a package called
CVODE from Lawrence Livermore National Lab as an alternative to the
Runge-Kutta method that should speed up the calculations. This version
should work better when you are trying to simulate equilibrium by using
fast rates, which is the situation in the problem you sent me earlier. You
should be able to use the enclosed version without having to increase the
rate constants to be on the order of the time step. I have done some
preliminary testing and the program seems to work, but you should do some
additional testing to make sure it is working correctly on your problems.
The enclosed program has identical functionality to PHREEQC, except this
version uses CVODE to integrate rate expressions, which could be slower
than RK in some (many?) cases. Please let me know if you find any bugs and,
if you use it, whether it works satisfactorily.

David

(See attached file: phreeqc.cvode.exe)


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

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

Attachment: phreeqc.cvode.exe
Description: Binary data



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

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: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/mail/msg00427.html
Email:dlpark@usgs.gov
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number [an error occurred while processing this directive] since Jan 22, 1998.