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

Re: PHREEQC problem

Your rate expression is pretty complicated, so I have not completely
analyzed your input file. I'll just make a couple comments.

(1) The formulas for the kinetic reactions must be charge balanced. From
what I see, you have
H -2 Cu +1 OH -1. This does not look right to me. Should it be CuOH 2? or
just Cu+1 H -2? There does not appear to be a need for the second kinetic
reaction because it is the same rate as the first, couldn't you do the same
thing with an additional OH -1 in the first rate?

(2) I think another  problem may be the fraction of a kilogram of water in
your solutions. TOT returns molality. However, the SAVEd value is moles.
You probably need to multiply all TOT functions by TOT("water")  or the
final moles by TOT("water")  to keep the units consistent. TOT("water") is
the mass of water in kilograms.

I'd try a first order rate expression first to make sure everything is
working as you expect before you go to the full expression that you have.


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

                      xiaomin Mao                                                                          
                      <xmao@xxxxxxxx>          To:       dlpark@xxxxxxxx                                   
                      05/25/03 06:21 AM        Subject:  PHREEQC problem                                   

Dear Parkhurst,

I am a research staff from Edinburgh Univ. in U.K. PHREEQC is a
programme I like, and recently we even coupled Modflow, Mt3dms with
PHREEQC for reactive transport modelling with dencity dependent flow.
Now, I have a trouble with PHREEQC modelling and need you help.
I am using PHREEQC for a transport modelling with kinetic sorption of Cu
by organic matter. The kinetic sorption was written based on Nica-Donna
model, and the parameter used is  from Christopher J. Milne (2003).
I am wondering why the simulation is so slow(or stopped) if the organic
matter get a little larger. Is it because the numerical convergence
problem? or, could you give me any suggestion on how to avoid such
thank you very much.
Dr. Xiaomin MAO
Institute for Infrastructure and Environment
School of Engineering and Electronics
The University of Edinburgh
Edinburgh, EH9 3JN

Tel: +44 131 650 5810
Fax: +44 131 667 9238

The input file is as follows,

10 Cu_conc=Tot("Cu")
20 H_conc=ACT("H+")
30 Org_xcue1=4.722*

33 Org_xcue2=0.8811*

35 rate1=-(0.0001*(Org_xcue1+Org_xcue2)-Kin("Org_xcu"))/10 #
0.0001--based on the quantity of organic matter
37 put(rate1,1)
40 moles=rate1 *time
45 if Tot("Cu")+moles < 0 then moles = -TOT("Cu")/3
50 save moles
35 rate1=get(1)
40 moles=-rate1*time
50 save moles

SOLUTION 0  CCA concentration
  units mol/kgw
  pH 7.0 charge
  pe 4
  -water 0.03375
   Na 1
   Cl 1

SOLUTION 1  CCA concentration
  units mol/kgw
  pH 4.0
  pe 4
  -water 0.03375

   Cu  0.123562
   Na 1
   Cl 1.3 charge

SOLUTION 2-16  CCA concentration
  units mol/kgw
  pH 4.0 charge
  pe 4
  -water 0.03375
  Cu  0.00
Na 1
Cl 1

kinetics 1-16
 -Formula  H -2.0 Cu +1
 -m0 0.00
 -m 0.00
 #-Parms 0.1
 -Formula H 0
 -m0 0.1
 -Parms 0.1

  -cells 16
  -shifts  33
  -time_step 30857  #10.5d/28*86400s/d
  -flow_direction forward
  -lengths 16*5e-3
  -boundary_conditions flux flux
  -dispersivities 16*0.012
  -correct_disp true
  -diffusion_coefficient 0.0000
  -punch_cells 16
  -punch_frequency 1
  -warnings false

  -file C:\cca_tmp\cca1
  -selected_out true
  -user_punch true
  -reset false

10  PUNCH ToT("Cu"), -LM("H+"),Kin("Org_xcu"),Kin("Org_xh")


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/msg00641.html
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 3068 since Jan 22, 1998.