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

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
problem?
thank you very much.
Xiaomin
----------------
Dr. Xiaomin MAO
Institute for Infrastructure and Environment
School of Engineering and Electronics
The University of Edinburgh
Edinburgh, EH9 3JN
U.K.

Tel: +44 131 650 5810
Fax: +44 131 667 9238
http://myprofile.cos.com/maoxm


The input file is as follows,

Rates
Org_xcu
-start
10 Cu_conc=Tot("Cu")
20 H_conc=ACT("H+")
30 Org_xcue1=4.722*
(1.8197*Cu_conc)^0.53/((1.8197*Cu_conc)^0.53+(218.78*H_conc)^0.66)
*((1.8197*Cu_conc)^0.53+(218.78*H_conc)^0.66)^0.59/(1+((1.8197*Cu_conc)^0.53+(218.78*H_conc)^0.66)^0.59)

33 Org_xcue2=0.8811*
(181970086*Cu_conc)^0.36/((181970086*Cu_conc)^0.36+(398107171*H_conc)^0.76)
*((181970086*Cu_conc)^0.36+(398107171*H_conc)^0.76)^0.7/(1+((181970086*Cu_conc)^0.36+(398107171*H_conc)^0.76)^0.7)

#0.0001*Cu_conc/(Cu_conc+H_conc)*parm(1)
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
-end
Org_xh
-start
35 rate1=get(1)
40 moles=-rate1*time
50 save moles
-end
end


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
 Org_xcu
 -Formula  H -2.0 Cu +1
 -m0 0.00
 -m 0.00
 #-Parms 0.1
 Org_xh
 -Formula H 0
 -m0 0.1
 -Parms 0.1

TRANSPORT
  -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

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


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

END





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/msg00640.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.