I haven't done it, but Tony Appelo describes how to set up radial problems
for PHREEQM in his book, Geochemistry groundwater and pollution (Balkema,
he is working on a second edition). The same approach applies to PHREEQC.
I think all you really have to do is define cell lengths that decrease with
distance by the appropriate amount.
From p. 423:
length(1) = R/ncell^.5, where length(1) is the length of the first cell,
ncell is the number of cells you want to use, and R is the maximum radial
distance.
length(n) = length(1)(n^.5 - (n-1)^.5), for cell n.
He gives an example on pager 437. If you have problems, he can probably
help you better than I (appt@xxxxxxxxx).
David
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
J Sam Arey
<arey@xxxxxxx> To: dlpark@xxxxxxxx
cc:
05/19/03 10:40 PM Subject: PHREEQ
Dear David,
I'm very interested to know if any version of PHREEQ is currently set
up to perform 1-d radial transport calculations, for either a
spherically or cylindrically symmetric system.
Thanks -
Sam Arey
--------------
48-336 ph 617 253 3897
Dept of Civil and Environmental Engineering fax 617 253 7475
Massachusetts Institute of Technology
Cambridge, MA 02139 USA
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/msg00637.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.