experimental condition range T=25-40C, pH=7.3-9.4
1 rem unit should be mol,kgw-1 and second-1
2 rem parm(1) is surface area in the unit of m2/kgw
3 rem calculation of surface area can be found in the note
4 rem M is current moles of minerals. M0 is the initial moles of minerals
5 rem parm(2) is a correction factor
10 rem acid solution parameters
11 a1=0
12 E1=0
13 n1=0
20 rem neutral solution parameters
21 a2=4.95E-09
22 E2=8700
23 n3=0.180
30 rem base solution parameters
31 a3=1.36E-16
32 E3=8700
33 n2=-1.48
36 rem rate=0 if no minerals and undersaturated
40 SR_mineral=SR("ORPIMENT")
41 if (M<0) then goto 200
42 if (M=0 and SR_mineral<1) then goto 200
43 if (M0<=0) then SA=PARM(1) else SA=PARM(1)*(M/M0)ˆ0.67
50 if (SA<=0) then SA=1
60 R=8.31451
75 Rate1=a1*EXP(-E1/R/TK)*ACT("H+")ˆn1 #acid rate expression
80 Rate2=a2*EXP(-E2/R/TK)*ACT("O2")ˆn3 #neutral rate expression
85 Rate3=a3*EXP(-E3/R/TK)*ACT("H+")ˆn2 #base rate expression
90 Rate=(Rate1+Rate3)*(1-Sr_mineral)*SA*parm(2)
100 moles= rate*Time
200 save moles
##############
#pyrite
############
pyrite