A Reactor Model For Gold Elution From
A Reactor Model For Gold Elution From
A Reactor Model For Gold Elution From
303-310,1995
Copyright 1995 Canadian Institute of Mining and Metallurgy
Printed in Great Britain. All rights reserved
0008-4433/95 $9.50+0.00
~ Pergamon
0008-4433(95)00020-8
Abstract-A model for the Zadra elution process was developed based on pore diffusion inside the carbon
particles, the dispersed plug flow in the elution column, and the gold desorption from carbon as a first
order reaction. A set of process description equations was obtained. Numerical solutions for the system
equations are given using the finite difference method, and the stability for the method was also analyzed.
Fitting of experimental elution concentration profiles using this model was presented. The model can fit
both the concentration profile and the recovery curve, and can simulate the elution process up to the final
stage. It was found that the elution process is controlled by a first order rate constant and the effective
intraparticle diffusivity in the carbon particle. The parameter analysis showed that the rate constant was
dependent on the elution temperature and the chemical composition of the eluant. The effective intraparticle
diffusivity was relatively constant under the experimental conditions.
NOMENCLATURE
BN,PN
C
C(t,y)
Ccal
Cexp
Cp
Cp(t,y, r)
dp
DI
De
D1B
k;
km
L
n
INTRODUCTION
N
Pe
Q
Qo
e,
r
rAu
R
Re
S
t
u
W
Wt
y
IX, {3, W
(J
r
8p
8
f1
P
Pc
>
S-1
constants
constriction factor for carbon pores
tortuosity of carbon pores
porosity of carbon particle
void fraction of the packed column
viscosity of solution (g cm- I S-I)
density of solution (g em -3)
density of carbon particles (g em:")
sum of the squares of the deviations (g em -3)
304
(rate of\
flow
in-out
)axialdispersion
(1)
rate of gold desorptiOn)
from particles
+(
rate of
= accumulation
The mass flow of gold cyanide per unit time may be defined as
follows:
entering by bulk flow = eliSCy;
-(EDIS~~)I ;
C)
.I'
-(eDIS~~)1
CJ
;
y+!'.y
Loaded Carbon
rate of desorption
rate of accumulation
Substituting
gives
of gold cyanide
ac
eS~Yat'
Elution
Column
C,H,-C,
-u ..
~)'
. +D1
~~I)HY-~~I,
+
~)'
l-e
--rAu
ac
= ---;::-.(2)
ct
Heat
~F(L)
Exchanger
,L----------:"<:.-.--~---------Eluted
Carbon
305
(8)
2
where W is the mass flux of gold cyanide (mol S-l cm- ) and
km is the mass transfer coefficient (cm S-l). Defining rAu as the
gold desorption rate per unit volume of carbon in mol S-l cmc-3,
i.e. rAu = W(4nR2)/(~nR3), we have
3km
rAu =:: -'R [Cplr=R - C].
Initialand
boundary conditions
C(t, y)
As ~y .~ 0, the above equation
oC
at =
becomes:
02C
oC
D1 oy2 -u oy
1-
(3)
+ -G-rAu'
(4)
where kc is the rate constant (S-l) and -oQ/ot
is the gold
desorption rate per unit weight of carbon.
Then, a mass balance (molls) for gold cyanide between rand
r + ~r in the spherical particle is
2)1
GP.
( -Dey
oCp4nr 2)1
2
The first
the normal
desorption.
pendent of
0;
(10)
at
(11)
Qo(r)
at
0,
(12)
-=
oy
at
y =L;
(13)
at
(14)
D10C
oQ
2
oCp
~rpCat = 4nr ~rGp7it.
OCp)_
or
we obtain
~ oQ.
Gp ot
C = - - u oy
(5)
the
(6)
Q(t, y, r)
r+~r
Gp
-4nr
at
ey
Cp(t,y,r)
( -D oCp4nr
r
(9)
+ ~ OCp) _
r or
~ oQ .
Gp ot
y = O.
Boundary conditions (13) and (14) are the Danckwerts boundary conditions [10]. The spherical symmetry condition is
applied to the carbon particles:
oCp
or
-=
at
=0
(15)
(16)
(7)
PARAMETER
ESTIMATION
306
P/.1ok
..
+(1-2w)C
P\I-l).j.(k-l)
(19)
w=D-e
Ci,j
(24)
1-8
)0.67
33
UO.
rt.=D1--,
II
I1t
(21)
where /1b and /12' are the liquid viscosities at temperatures T1
and T2, respectively.
(25)
(26)
(11.1')2
~t
211y'
(27)
(20)
1'Au,
where
13 =
8;p
I1t
(111')2 .
= (rt.-!3)Cu-1),u+I)+(1-2rt.)Cu-I),j
DO
PI/-I}.).,'.()
km = 1.09(
k+l
+w--Ck
where
..
P\I-I).j);
(18)
k-l
C .. =w--Ck
..
l.j.P:-;
D1
-~-
~
II
l1y ~ 2'
(29)
2'
(30)
307
(33)
Qo (kg t-I)
Bp
B
rc (g cm-3)
Sample A
Sample B
Sample C
9.63
0.55
0.4
0.45
8.03
0.55
0.4
0.45
4.61
0.55
0.4
0.45
I
2'
Equations (3), (4) and (7) were used to calculate the gold
concentration profile and recovery as a function of time. Figures
4(a), 5(a) and 6(a) demonstrate the gold concentration profile
of the eluant at the exit of the column. The corresponding
recovery curves for the elution process are shown in Figs 4(b),
5(b) and 6(b). The results demonstrate that the calculated curves
agree well with the experimental data, both in concentration
profiles and recovery curves. The recovery curve is the integrated desorption curve of the gold concentration against the
bed volume of the elution. It is apparent that both De and kc
(31)
1500
1350
....J
........
0'1
1200
.2 1050
"
C
Q)
k =0.70
900
-1
-5
2
em
/s
750
::J
De=1.11x10
experimental
fitted curve
600
Q)
::J
W
300
150
0
<I> =
450
[Cexp -
Ccalf,
(32)
where Cexp is the experimental result and Cca1 is the gold concentration predicted by the model at the exit of the elution
column.
To minimize the objective function <I> for the given search
parameters, the combination of two series of search values, De
and kc' are calculated on an IBM mainframe with a FORTRAN
program.
Experimental data from a laboratory elution column (2.5 cm
dia., 20 cm high) were used to examine the model. A mass of
50 g of loaded carbon was packed in the column, which was
then filled with water and heated to the desired temperature
with an electrical heating pad. The temperature was maintained
at 2C with a Parr temperature controller. The eluant was
prepared by dissolving analytical grade reagents in de-ionized
water. The elution flow rate of two bed volumes per hour was
controlled by compressed nitrogen gas.
Three different loaded carbon samples were obtained from
the production CIP circuits of Canadian gold mines. The specification of the carbons is shown in Table 1. Other parameters
used in the model, such as particle porosity, void fraction of
packed column and density were assumed to be the same for
all samples at no measurements.
The average equivalent particle radii, R, were determined by
measuring the weight (Wt) of about 200 randomly selected
particles (n) and calculated by the following equation:
4
Time
k=1
(Hour)
Fig. 4(a). Elution of gold from the carbon with NaCN-NaOH : NaCN,
1 g 1-1; NaOH, 10 g 1-1; 148C; Q = 9.63 kg t-1 Au.
100
90
~
-0.
80
experimental
Calculated
70
60
k=0.70
()
'0
>
-1
-5
De=1.11x10
(])
em
/s
50
(])
40
30
20
10
4
Time
(Hour)
Fig. 4(b). Elution of gold from the carbon with NaCN-NaOH: NaCN,
1 g 1-1; NaOH, 10 g 1-1; 148C; Q = 9.63 kg t-I Au.
308
400
1500
1350
....J
....J
"-
"-
01
E
c-
.2
1200
experimental
fitted
c-
curve
1050
~
~
"
C
Q)
900
k =0.78 h
c
De=1.11x10
750
-1
C
Q)
-5
/s
cm
250
k =0.40
c
-0
W
100
/s
0
2
345
Time
(Hour)
NaCN,
(Hour)
Fig. 6(a). Elution of gold from the carbon with NaCN-NaOH : NaCN,
2 g 1-1; NaOH, 10 g 1-1; 140:lC; 2 BV h-I.
100
100
0000
90
90
80
experimental
70
Calculated
"0
""6
<.:>
""6
<.:>
60
k =0.78
De=1.11x10
50
'0
h-
-5
2
cm
/s
Q)
>
o
--
80
experimental
Calculated
70
60
k =0.40
-1
-5
De=1.11x10
50
2
cm
/s
u
Q)
cm
o
Time
lY
-5
50
>
-1
150
Q)
:J
450
150
curve
De=1.11x10
200
300
'0
fitted
:J
Q)
-0
-6
experimental
<{
600
Q)
:J
W
o
-
300
:J
<{
350
01
Q)
40
0:::
40
30
30
20
20
10
10
345
Time
(Hour)
4
Time
NaCN,
of gold elution in a
DISCUSSION
The numerical solution to the pore diffusion-axial dispersion
model, using the finite difference method, fits the experimental
elution profiles and the recovery curves of the pressurized Zadra
process. The finite difference approximation
of the model is
accurate as long as the convergence criteria [eqns (29}-(31)] are
satisfied. The agreement between the experimental and fitted
concentration profiles also indicates that the difference equation
(Hour)
NaCN,
system [eqns (3), (4) and (7)] is accurate for the gold elution
column.
A long tail in the elution profile, as shown in Figs 4(a), 5(a)
and 6(a), indicates that the elution of a small portion of the
gold cyanide is delayed by the pores inside the core of the
carbon particles. The demonstrated ability to fit the long tail is
seen as an indication of the importance of the porous diffusion
in the elution of gold cyanide from loaded carbon. The suggested model can simulate the complete elution process. Other
carbon elution models, such as that seen in Ref. [4], are hardly
applicable to the final stage of elution, owing to neglect of the
porous diffusion inside carbon particles.
D1B(J
=--
309
Acknowledgements-The research was partially supported by the Natural Science and Engineering Research Council of Canada (NSERC).
The authors wish to thank Dr R. M. Erdahl in the Department of
Mathematics and Statistics in Queen's University for many helpful
discussions.
(34)
,
!
where (J is the constriction factor, which accounts for the variation in the cross-sectional area of carbon particle pores and is
a function of the ratio of maximum to minimum pore areas,
and! is the tortuosity.
Typical values of the system are: ! = 3-5, (J = 0.8 [9], so
De/D1B = 0.26-0.16. For the experiments, the ratio of De/D1B
is 0.22, which indicates that the pores inside the carbon particles
reduce the diffusivity of gold cyanide complex within the carbon
particles to about 22 % of that in a normal dilute solution. This
further confirms the importance of pore diffusion within the
carbon particles for the gold elution process.
CONCLUSIONS
A model for the Zadra carbon elution process was developed
and the system description equations and the initial and boundary conditions have been defined. In this model, the gold
* BV h-1,
9.63
148
1.0
10
2.0
0.70
1.11
Sample B
8.03
148
4.0
10
2.0
0.78
1.11
Sample C
4.61
140
2.0
10
2.0
0.40
1.11
REFERENCES
1. J. B. Zadra, A. L. Engel and H. J. Heinen, Rep. Invest. Bur. Mines.
4843 (1952).
2. R. J. Davidson and D. Duncanson, J. S. Afr. Inst. Min. Metall. 77,
254 (1977).
3. D. M. Muir, Hydrometallurgy 14(1), 47 (1985).
4. M. D. Adams and M. J. Nicol, Gold 100 Proc. Int. Con! on
Gold, Vol. 2, Extractive Metallurgy of Gold. SAIMM, p. 111,
1986.
5. M. D. Adams, Trans. Inst. Min. Metall., Sect. C. 99, C71 (1990).
6. N. M. Vegter and R. F. Sandenbergh, Hydrometallurgy 28(2), 205
(1992).
7. W. Stange, Miner. Engng 4(12), 1279 (1991).
8. J. S. J. Van Deventer and V. E. Ross, Miner. Engng 4(7-11), 667
(1991).
9. H. Scott Fogler, Elements of Chemical Reaction Engineering, 2nd
Edn. Prentice-Hall, pp. 543-92, 759-94, 1992.
10. P. V. Danckerts, Chem. Engng Sci. 2, 1 (1953).
11. N. M. Vegter, Hydrometallurgy 30, 229 (1992).
12. C. Y. Wen and L. T. Fan, Modelsfor Flow Systems and Chemical
Reactors. Marcel Dekker, p. 570, 1975.
13. E. J. Wilson and C. J. Geankoplis, Ind. Engng Chem. Fundam. 1,9
(1966).
14. G. M. Schmid and M. E. Curley-Fiorino, Encyclopedia of Electrochemistry of the Elements, Vol. IV. Marcel Dekker, pp. 87-178,
1975.
15. B. Carnahan, H. A. Luther and J. O. Wilkes, Applied Numerical
Methods. John Wiley & Sons, NY, pp. 429-519, 1969.
16. A. Constantinides, Applied Numerical Methods with Personal Computers. McGraw-Hill, pp. 458-519, 1987.
17. B. A. Finlayson, Nonlinear Analysis in Chemical Engineering.
McGraw-Hill, p. 216, 1980.
18. L. Lapidus and G. F. Pinder, Numerical Solution of Partial Differential Equations in Science and Engineering. John Wiley & Sons,
NY, pp. 149-350, 1982.
19. G. J. McDougall and C. A. Fleming, in Ion Exchange and Sorption
Processes in Hydrometallurgy (edited by M. Streat and D. Naden).
The Society of Chemical Industry, pp. 56-126, 1987.
310
GOLD
ELUTION
FROM
ACTIVATED
CARBON
Interface
APPENDIX
Finite-difference
solution
approximations
When the gold cyanide diffuses from the core of the carbon particle
to its surface, the mass flow at the interface between carbon and solution
behaves as shown in Fig. AI. We wish to derive the relevant finitedifference approximation
for gold concentration
Cp at point PN on the
interface between carbon and solution. The following procedure is
based on the continuity of molar flux at the interface.
In carbon (from Taylor's expansion) we have, approximately,
C
Pi,j.(P:-i-I)
The subscript
i.e. :
_
Pi.j.P:-i
dr(aCp)
ar
PN.c
l'
2
(A)2
~--
PN.c
[C
(ac)ar
+dr
Pi,j.P:-;
is approximated
p)
p)
PN . .11'
cr
Pc Qu
PN,c
...I),j.PN -
Qi,j.PN
(A4)
dt
l;p
(A2)
CPU.I),j.P:-; =
PN,c
[rc. +?wC
~
I,J
P(i-I).j.(P:-;-1l
- ~(Q('.
p
1-
-Q. t,},'PN)]I(1
I) .p .
,J,.~
+C+2w)
-
by
Cp(i+ll.j.P:-i
PN,c -
dt
CPi.j.P:-;
D(_2_[c
-_Cp-'----ij_.P:-;
=
dt
e (dr)2
_C_P(i+_!l-'----.j.P_:-i
_2_ (8:
(A3)
where
v
Substitution
PN
(AS)
ac
( at
Fig. A I. Interface
(AI)
PN,c
_P
IJ.r
'" 2
-C
Pi.j.(P:-i-ll
p
C )
or
-_P
PN-2 PN-l
+ (dr)2
2
(aaC)2
(a
solution
Carbon
-C
P(i+l).j.(P:-;-I)
P(i+I).j.P:-;
+.1r(_aC_p)]
ar PN.c
~-
2(1 +PN)dt
PN.1r
Dedt
w=--.
(dr):!
k
m
(A6)
(A7)