A shrinking particle - shrinking core model for leaching of a zinc ore containing. Vida Safari, Gilnaz Arzpeyma, Fereshteh Rashchi, Navid Mostoufi

Size: px

Start display at page:

Download "A shrinking particle - shrinking core model for leaching of a zinc ore containing. Vida Safari, Gilnaz Arzpeyma, Fereshteh Rashchi, Navid Mostoufi"

Bertha Armstrong
6 years ago
Views:

Accepted Manuscript A shrinking particle - shrinking core model for leaching of a zinc ore containing silica Vida Safari, Gilnaz Arzpeyma, Fereshteh Rashchi, Navid Mostoufi PII: S0301-7516(09)00142-2

1 Accepted Manuscript A shrinking particle - shrinking core model for leaching of a zinc ore containing silica Vida Safari, Gilnaz Arzpeyma, Fereshteh Rashchi, Navid Mostoufi PII: S (09) DOI: doi: /j.minpro Reference: MINPRO 2179 To appear in: International Journal of Mineral Processing Received date: 10 January 2009 Revised date: 24 May 2009 Accepted date: 11 June 2009 Please cite this article as: Safari, Vida, Arzpeyma, Gilnaz, Rashchi, Fereshteh, Mostoufi, Navid, A shrinking particle - shrinking core model for leaching of a zinc ore containing silica, International Journal of Mineral Processing (2009), doi: /j.minpro This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

2 A shrinking particle - shrinking core model for leaching of a zinc ore containing silica Vida Safari 1, Gilnaz Arzpeyma 1, Fereshteh Rashchi 2, Navid Mostoufi 1,* 1 Department of Chemical Engineering, 2 Department of Metallurgy and Materials Engineering, University of Tehran, PO Box 11155/4563, Tehran, Iran Abstract A new mathematical model was developed for leaching of zinc ores containing silicates such as hemimorphite which produce a gel during leaching with sulfuric acid. This model is based on the shrinking core model in which the particle size and the reacting core shrink simultaneously. It was shown that the actual dissolution time of the ore particles is longer than the time corresponding to the dissolution of chemical zinc oxide itself. It was suggested that because of the existence of silicates in the ore, a gelatinous layer was formed around the reacting core. Since the gel product is soft, it breaks apart when the particles collide and as a result, the particles shrink. However, a thin gelatinous layer always covers the reacting core which increases the mass transfer resistance and increases the leaching time. This model was applied to leaching of a zinc-rich tailing containing hemimorphite and the thickness of the gelatinous layer as well as the diffusion coefficient in this layer was determined. Keywords: Leaching; Kinetics; Zinc silicate; Shrinking core model; gelatinous silica layer 1. Introduction The main source of zinc metal production is zinc sulfide ore. Currently, depletion of these sulfide ores has brought more emphasis on zinc extraction from oxides, silicates or even * Corresponding author, Phone: (98-21) ; Fax: (98-21) ; address: mostoufi@ut.ac.ir 1

3 secondary sources. Zinc extraction is performed mainly by hydrometallurgical methods. In the hydrometallurgical process, the ore is first leached by a solvent and then before electrowinning a purification process is used to prepare the solution for eletrolysis. During the leaching process of zinc oxidized ore, soluble zinc sulfate forms which stays in solution. In this process, the lead compounds form lead sulfate precipitates which transfer to the leaching filter cake during the solid/liquid separation. Leaching of the ore at ph ca. 2, transforms the silicate compounds of the ore to colloidal silica, i.e., a gel (Matthew and Elsner, 1997). Process kinetics and optimum operating conditions have been studied for the leaching of zinc silicate ore tailings. Monhemius and Terry (1983) investigated the influence of different parameters on the kinetic of acid dissolution of both natural and synthetic willemite (Zn 2 SiO 4 ) and hemimorphite (Zn 4 Si 2 O 7 (OH) 2.H 2 O). Specific rate constants were estimated for leaching of both willemite and hemimorphite in different acidic media. They found that the dissolution was mixed chemical/diffusion controlled in hemimorphite and chemically controlled in the case of willemite. Abdel-Aal (2000) investigated the kinetics of sulfuric acid leaching of low- grade zinc silicate ore. In their study, diffusion through the product layer was determined as the rate controlling step. Espiari et al. (2006) verified extraction of zinc from tailings of lead flotation plant. They studied the effect of different parameters on kinetics of zinc dissolution and found that the rate determining step is the physico-chemical desorption process. Souza et al. (2007) studied the effect of particle size, temperature and initial acid concentration on leaching of zinc silicate ores. They concluded that the grain model with porous diffusion control is the rate controlling step. 2

4 Some researchers have reported the possible existence of internal diffusion resistance. For instance, Pecina et al. (2007) have reported formation of a sulfur layer in zinc sulfide leaching with high concentration of sulfuric acid solution containing hydrogen peroxide and indicated that this layer reduces leaching efficiency. Mulak et al. (2005) found that in leaching of spent nickel oxide catalyst with sulfuric acid, an aluminium-rich layer surrounds the unreacted core of the particle and grows inward as the particle reacts. In leaching process of sphalerite containing a lower concentration of iron, the decrease in zinc dissolution rate (as compared to a greater iron-containing sphalerite mineral) was attributed to formation and growth of a polysulfide surface layer during the initial rapid leach period (Weisener et al., 2004). Based on the above evidences, it can be concluded that when a core shrinks, an internal resistance layer, either a reaction product layer or a gel film forms around the core and results in a decrease in the extraction yield in the leaching process. However, in many cases, the size of the particle (including unreacted core and the layer) decreases with time. In other words, although the product layer forms around the core, it shrinks as the core shrinks. In all the previous studies, existence of this gelatinous layer around the particle was neglected during the kinetics calculations and its effect on the kinetics of leaching was not taken into account. In the present study, a mathematical model has been developed based on the shrinking core model in which the resistance of the gelatinous product film is also considered. This kinetic model consists of three steps: external diffusion in the liquid, internal diffusion in the gelatinous product film and chemical reaction on the surface of the core. 2. Model Development 3

5 Zinc silicate ore, as hemimorphite, reacts with sulfuric acid according to the following reaction: Zn Si O ( OH ). H O + 4H SO 4 ZnSO + Si O( OH ) + 3H O (1) Zinc sulfate is soluble in water. However, disilicic acid, Si 2 O(OH) 6, which apparently polymerizes to produce polysilicic acid, forms a gelatinous phase at specific acidic ph which remains around the particles surface and leads to a decrease in the extraction yield. The most common models considered in leaching are illustrated in Fig. 1. Fig. 1a corresponds to the case when the reaction takes place on the exposed surface of the particle and the product completely dissolves in the liquid. This shrinking particle model (SPM) has been used by researchers such as Espiari et al. (2006), Aydoğan et al. (2006) and Velardo et al. (2002). If the product does not dissolve in the liquid, the particle size would not change but the reacting core shrinks inside the particle. This situation is shown in Fig. 1b and the model is called Shrinking Core-Constant Particle Size. It has been used by researchers such as Liu et al. (2006), Liddell (2005) and Szubert et al. (2006). Fig 1c demonstrates schematic of a model called Shrinking Core-Shrinking Particle. In this case as the reaction proceeds, the unreacted core of particle shrinks while a gelatinous silica layer forms around the core. However, since this layer is soft, it breaks apart when the particles collide. Nevertheless, a thin layer of silica remains around the core. The silica layer creates a resistance during acid transfer from the solution to the surface of the core. In the present work, the last model was considered as the base kinetic model since the silica product does not dissolve completely in the acid solution. Therefore, the internal diffusion through the gel film should be taken into account Assumptions 4

6 The assumptions of the model are as follows: The particles are spherical. During the process, the particle shrinks uniformly, thus, it maintains its spherical shape. In the absence of adequate information, the thickness of the silica layer around the core was assumed to be constant during the leaching process. Although there exist many reactions in the leaching of zinc by sulfuric acid from the ore, for the sake of simplicity, the main reaction considered in this work is dissolution of zinc oxide in acid. In other words, the main source of zinc in the leaching was assumed to be zinc oxide. Hemimorphite particles were considered as the source of silica responsible for the gel formation. Other substances present in the ore do not have any significant effect on the kinetics. The temperature remains constant during the process. The particle and the gelatinous layer are both non-porous. Thus, mass transfer occurs through molecular/ion diffusion in these phases Kinetic Modeling Considering all the above assumptions, the first step for developing the model is to define a criterion indicative of the advancement of reaction versus time. The rate of reaction per unit surface of the core can be related to the dissolution rate of zinc oxide as follows: R r M S = ZnO ZnO (2) c dn dt The rate of zinc oxide disappearance can be expressed as: 5

7 dnzno ρoregs c drc = (3) dt M dt ZnO It has been shown in several references that the rate of reaction for dissolution of zinc oxide, R r, is first order with respect to the concentration of the solvent (e.g., Espiari et al., 2006; Monhemius and Terry, 1983): R kc A r = (4) M ZnO Therefore, the rate of shrinkage can be expressed as: drc dt kc ρ G A = (5) Ore The acid concentration, C A, used in Eq. (5) should be evaluated at the surface of the core. To find the acid concentration at reaction surface, mass balances should be written for both the liquid film layer and the gelatinous layer. In both cases, mass transfer occurs only in radial direction and the mass transfer equation becomes: 1 d dc A 2 r 2 = r dr dr The boundary conditions are: 0 dc = + δ = (7) ( ) A r rc hd C A C 0 A D g dr (6) dc A r = rc D g = kc A (8) dr Considering that the process is in quasi-steady-state conditions, the rate of reaction at the surface of the core would be equal to the rate of mass transfer to and through the gel film: ( 0 r +δ ) h C C = k C (9) D A A A c rc 6

8 Solving the differential Eq. (6) with boundary conditions (7) and (8) and then solving Eq. (9) for finding C A (r c ) would result in obtaining C A : C A = C A0 Dg krc rc r D 1 1 D + + kr r r h r g g 2 2 c c c + δ D( c + δ ) Therefore, the dissolution rate of zinc oxide can be expressed as: R r k C A r = r C 1 c A0 = = 2 M ZnO M ZnO 1 r c r c rc k Dg rc + δ hd ( rc + δ) or in form of the shrinkage rate: dr C c A 1 0 = 2 dt ρoreg 1 r c r c rc k D g rc + δ hd ( rc + δ ) 2 2 (10) (11) (12) Eq. (12) can be expressed in terms of advancement of the reaction. The conversion can be determined based on the residual volume of particle as follows: rc ( ) 1 ( t X t ) = r0 3 (13) Therefore, the final differential equation, from which the extent of dissolution of particles as a function of time can be obtained, can be then achieved by combining Eqs. (12) and (13): dx dt 3C A 1 = 3 ρ 2/ 3 OreGr0 2 2/ ( 1 X ) r0 1 1 (1 X ) r / 3 1/ 3 k Dg r0 (1 X ) r0 (1 X ) + δ hd r (1 X ) + 1/3 ( 0 δ ) 2 (14) 7

9 2.3. Chemical and Physical Properties The only chemical reaction considered in this work on the surface of the core is dissolution of zinc oxide by sulfuric acid. Reaction rate constant for dissolution of pure zinc oxide in sulfuric acid is calculated from the following equation (Kristovnikov and Davydovskaya, 1936): k = exp (m/s) (15) RT It is worth mentioning that in all previous investigations on kinetics of leaching of zinc with sulfuric acid, different rate constants were reported (e.g., Espiari et al., 2006; Abdel-Aal., 2000). However, since in all these investigations the effect of gelatinous layer was not taken into account, the reported rate constant is in fact a combination of reaction rate constant and dispersion coefficient of acid in the gelatinous layer. In the present study, the effect of gelatinous layer is separated from the reaction. Therefore, only dissolution of pure zinc oxide in sulfuric acid was considered on the surface of the ore. Of course, based on the assumptions listed above, effect of other substances in the ore, on the dissolution kinetic of zinc oxide was neglected in this work. The mass transfer coefficient was calculated from (Ranz and Marshall, 1952): Sh 1/3 1/2 1/3 1/2 µ d pu p = Sc Re = ρ D µ (16) Evaluation of mass transfer coefficient from Eq. (16) requires estimation of viscosity, density and diffusion coefficient of the acid solution. These properties are not simple functions of their compositions. There can be found several equations for viscosity, density and diffusion coefficient of the mixture in this process in various literatures. In the present study, the following correlations of Guerra et al. (2006) were used: 8

10 µ = ( T T T [ ] [ ] + T ZnSO + H SO ) exp( (17) [ Fe ] [ H SO ] Zn Fe + + ρ = [ ] ZnSO4 + T 0.5 ([ 4 ] [ 2 4 ]) [ 4 ] [ H 2SO4 ] ([ ZnSO 4 ] [ H 2SO4 ]) D = [ ZnSO H SO ZnSO ] 10 (18) (19) 3. Results and discussion Performance of the proposed model was examined using the experimental data reported by Espiari et al. (2006) for leaching of zinc from a zinc-rich oxide silicate tailing with sulfuric acid. Their XRF analysis results showed that the sample contained 37% zinc oxide and 23.7% silica. The leaching data reported by Espiari et al. (2006) are available as zinc recovery vs. time at various temperatures. Eq. (14) was solved for different operating conditions reported by Espiari et al. (2006). In this equation, thickness and diffusion coefficient of the gel layer were considered as fitting parameters and their values at different temperatures were determined by fitting the equation to the experimental data. Fig. 2 illustrates a sample solution of the model as well as the corresponding experimental data. This figure shows the conversion of the ore as a function of time in the batch system. It can be seen in Fig. 2 that the model fits satisfactorily to the experimental data (Espiari et al., 2006). Prediction of the conversion when the effect of the gel film was neglected (only zinc oxide reaction with acid and liquid film resistance was considered) is also shown in the same figure. As it can be seen, the particles dissolved very fast if the effect of gel film formation is neglected. However, the trend of the experimental data suggests that the leaching process is not as fast as dissolution of zinc oxide in sulfuric acid. Slower reaction can be justified by 9

11 adding the resistance of the gel film which is modeled in the present study. Therefore, it is not possible to neglect the effect of the gel film because its effect on the leaching is substantial. Espiari et al. (2006) reported conversion vs. time at various temperatures. Using these experimental data it was found that the thickness of the layer is almost constant in the range of temperature considered in this work and its average was determined to be 1.2 µm with standard deviation of 0.8 µm. However, diffusion coefficient is a strong function of temperature. The Arrhenius plot of diffusion coefficient of the gel against temperature is shown in Fig. 3 from which temperature dependency of diffusion coefficient was found to be: ln D g = (20) RT There are three terms in the denominator of Eq. (14) which correspond to resistances due to chemical reaction, diffusion through gel film and mass transfer in the liquid film, respectively. Resistance of the liquid film is negligible as compared to the other two resistances for the operating conditions considered in this work. Therefore, in order to investigate the effect of temperature on the leaching rate, only chemical reaction and diffusion through the gel were considered. Fig. 4 demonstrates the reaction and internal diffusion resistances as a function of temperature. It can be seen in this figure that both of the resistances decrease by increasing the temperature. At low temperatures, reaction resistance is negligible and kinetic mechanism would be reduced to gel diffusion control. Increasing the temperature has a significant effect on both reaction and mass transfer rates. At high temperatures these two rates are of the same order of magnitude as shown in figure 4. Therefore, at high temperature, the process is controlled by both chemical reaction and diffusion though the gelatinous layer. 10

12 4. Conclusion A mathematical model was developed for leaching of zinc from a zinc ore containing silica. It was shown that considering only the chemical reaction of zinc dissolution is not sufficient for estimating the leaching time of the ore. A thin gelatinous layer was considered to cover the reacting core of the ore in order to correct the model prediction. Adding the mass transfer resistance of this layer to the model considerably improved the predictions. It was shown that the thickness of the gel is almost constant but the diffusion coefficient of the gel decreases with temperature. At low temperatures, the reaction rate is significantly lower than the rate of mass transfer through this layer which alters the mechanism to diffusion controlled. At high temperatures, both chemical reaction and mass transfer control the dissolution rate of the zinc ore particles containing silica. Acknowledgement The authors would like to thank Professor Fathi Habashi from University of Laval, Canada, for his valuable comments during the work. Nomenclature C A acid concentration at the surface of the core (kg.m -3 ) C A0 acid concentration in the bulk (kg.m -3 ) d p particle diameter (m) D g diffusivity of gel (m 2.s -1 ) D ZnSO4 diffusion coefficient of zinc sulfate in solution (m 2.s -1 ) G zinc oxide grade h D mass transfer coefficient (m.s -1 ) 11

13 k reaction rate constant (m.s -1 ) M ZnO zinc oxide molecular weight (kg.kmol -1 ) n ZnO r c r 0 moles of zinc oxide (kmol) core radius (m) initial radius of the particle (m) R gas constant (8.314 J.mol -1.K -1 ) R r rate of reaction for dissolution of zinc oxide (kg.m -2.s -1 ) Re Reynolds number; Re=ρu p d p /µ S c surface of the core (m 2 ) Sc Schmidt number; Sc=µ/ρD ZnSO4 Sh Sherwood number; Sh=h D d p /µ t time (sec.) T temperature (K) u p terminal velocity (m.s -1 ) X volumetric conversion Greek letters δ gelatinous layer thickness (m) µ viscosity of solution (kg.m -1.s -1 ) ρ Ore density of zinc ore (kg.m -3 ) ρ density of solution (kg.m -3 ) References Abdel-Aal, E. A., Kinetics of sulfuric acid leaching of low-grade zinc silicate ore. Hydrometallurgy 55,

14 Aydoğan, S., Erdemoğlu, M., Aras, A., Uçar, G., Özkan, A., Dissolution kinetics of celestite (SrSO4) in HCl solution with BaCl2, Hydrometallurgy 84, Espiari, S., Rashchi, F., Sadrnezhaad, S. K., Hydrometallurgical treatment of tailings with high zinc content. Hydrometallurgy 82, Guerra, E., Bestetti, M., Physicochemical properties of ZnSO 4 -H 2 SO 4 -H 2 O electrolytes of relevance to zinc electrowinning, J. Chem. Eng. Data 51, Kristovnikov, A. H., Davydovskaya, E. A., Zh. Fiz. Khim (In Russian). 8, Liddell, K. C., Shrinking core models in hydrometallurgy: What students are not being told about the pseudo-steady approximation, Hydrometallurgy 79, Liu, Y., Qi, T., Chu, J., Tong, Q., Zhang, Y., Decomposition of ilmenite by concentrated KOH solution under atmospheric pressure, Int. J. Miner. Process. 81, Matthew, G., Elsner, D., Hydrometallurgical treatment of zinc silicate ores, Metallurgical Transactions 8B, Monhemius, A. J., Terry., Acid dissolution of willemite and hemimorphite, Metallurgical Transactions 14B, Mulak, W., Miazga, B., Szymczycha, A., Kinetics of nickel leaching from spent catalyst in sulphuric acid solution, Int. J. Miner Process 77, Pecina, T., Franco, T., Castillo, P., Orrantia, E., Leaching of zinc concentrate in H 2 SO 4 solutions containing H 2 O 2 and complexing agents, Minerals Eng. 21(1), Ranz, W. E., Marshall, W. R. 1952, Evaporation from drops: part 1, Chem. Eng. Prog. 48, Souza, A. D., Peina, P. S., Lima, E. V. O., dasilva, C. A., Leão, V. A., Kinetics of sulphuric acid leaching of a zinc silicate calcine, Hydrometallurgy 89,

15 Szubert, A., Łupiński, M., Sadowski, Z., Application of shrinking core model to bioleaching of black shale particles, Physicochemical Problems of Mineral Processing, 40, Velardo, A., Giona, M., Adrover, A., Pagnanelli, F., Toro, L., Two-layer shrinkingcore model: parameter estimation for the reaction order in leaching processes, Chem. Eng. J. 90, Weisener, C. G., Smart, R. St. C., Gerson, A. R., A comparison of the kinetics and mechanism of acid leaching of sphalerite containing low and high concentrations of iron, Int. J. Miner. Process 74,

16 Figure Captions Figure 1. Schematics of different mechanisms of leaching Figure 2. Model verification with and without gelatinous layer Figure 3. Gel diffusivity versus temperature Figure 4. Comparing reaction and gel diffusion resistances. 15

17 (a) Shrinking Particle (b) Shrinking Core Constant Particle Size (c) Shrinking Core Shrinking Particle Leaching Time Figure 1. Schematics of different mechanisms of leaching 16

18 X Espiari et al. (2006) Present model Neglecting gel time (min) Figure 2. Model verification with and without gelatinous layer 17

19 ln D g /T (K -1 ) Figure 3. Gel diffusivity versus temperature 18

20 Resistance (min./m) Diffusion Reaction Temperature (K) Figure 4. Comparing reaction and gel diffusion resistances. 19

DISSOLUTION KINETICS OF SMITHSONITE IN BORIC ACID SOLUTIONS

DISSOLUTION KINETICS OF SMITHSONITE IN BORIC ACID SOLUTIONS Physicochem. Probl. Miner. Process. 53(1), 2017, 161 172 www.minproc.pwr.wroc.pl/journal/ Physicochemical Problems of Mineral Processing ISSN 1643-1049 (print) ISSN 2084-4735 (online) Received February