Research Article Numerical Analysis on Flow and Solute Transmission during Heap Leaching Processes

Transcription

1 Mathematical Problems in Engineering Volume 25, Article ID , 5 pages Research Article Numerical Analysis on Flow and Solute Transmission during Heap Leaching Processes J. Z. Liu, A. X. Wu, 2 and L. W. Zhang College of Information Technology, Shanghai Ocean University, Shanghai 236, China 2 School of Civil and Environment Engineering, University of Science and Technology Beijing, Beijing 83, China Correspondence should be addressed to J. Z. Liu; jzliu@shou.edu.cn and A. X. Wu; liujinzhi2cn@sina.com Received 26 August 24; Accepted 3 October 24 Academic Editor: Kim M. Liew Copyright 25 J. Z. Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Based on fluid flow and rock skeleton elastic deformation during heap leaching process, a deformation-flow coupling model is developed. Regarding a leaching column with m height, solution concentration unit, and the leaching time being days, numerical simulations and indoors experiment are conducted, respectively. Numerical results indicate that volumetric strain and concentration of solvent decrease with bed s depth increasing; while the concentration of dissolved mineral increases firstly and decreases from a certain position, the peak values of concentration curves move leftward with time. The comparison between experimental results and numerical solutions is given, which shows these two are in agreement on the whole trend.. Introduction Solution mining is conceptualized as the removal of dissolved metals from original solid matrix [ 3]. In general, in situ leaching and heap leaching are adopted, and the latter is more often used. During heap leaching processes, factors, such as fluid flow, pore pressure, chemical or biochemical reaction between target metals and leaching solution, target metals dissolution, and reaction byproduct deposition, all result in deformation of the heap, affecting the leaching rate [4]. Of all these factors, elastic deformation caused by pore pressure is the main skeleton deformation. In recent years, some mathematical models have been developed to describe the processes of heap leaching. Bouffard and Dixon studied the hydrodynamics of heap leaching processes deeply. They derived three mathematical models in dimensionless form to simulate the transport of solutes through the flowing channels and the stagnant pores of an unsaturated heap [5]. Lasaga investigated the chemical kinetics of water-rock interactions and gave the description of rock deformation regularity [6]. Solute transport and flow through porous media with applications to heap leaching of copper were studied deeply [7 9]. Sheikhzadeh et al. developed an unsteady and two-dimensional model based on the mass conservation equations of liquid phase in the ore bed and in the ore particle, respectively. The model equations were solved using a fully implicit finite difference method, and the results gave the distributions of the degree of saturation and the vertical flowing velocity in the bed []. Wu et al. built the basic equations describing the mass transmission in heap leaching. They gave the analytic solution omitting convection with small application rate and determine the hydrodiffusion coefficient []. The models discussed above concentrated on the steady flowing conditions without considering the effect of elastic deformation. The purpose of this work is to apply an elastic deformation model for simulating the column leaching processes and develop the governing equations of coupled flow and deformation behavior with mass transfer. These equations are solved numerically by Comsol Multiphysics. The changeable regularity of volumetric strain and concentration distributionsofthesolventandthesoluteisgiven.thevalidationof the mathematical model and numerical analysis is concerned through experiment.

2 2 Mathematical Problems in Engineering 2. Model Development 2.. Flow and Solid Elastic Deformation Model. Supposing the solution flows in a deformational and homogeneous porous medium, the basic seepage equation for column leaching is [2] [χ p ( n) +χ f n] p + [ κ η ( p + ρ fg e)]+ ε V =Q s, () where χ p, χ f are pore deformation coefficient and fluid deformation coefficient, p is liquid pressure, e is elevation, κ is permeability, η is viscosity, g is acceleration of gravity, ρ f istheliquiddensity,ε V is volumetric strain, Q s is the source term, and n is the porosity. Solid elastic deformation equations describing the plain strain deformation state are [3] where y is the axis along ore column; t is leaching time; C, C 2 are the concentrations of reagent and dissolved metal; s, which can be written as s(y, t), is the absorbed solute mass on unit pore area; u is the flowing velocity; D is the dispersion coefficient; b is opening width of pore; J d is diffusion flux; R i is chemical reaction rate; β is the stoichiometric coefficient. The chemical reaction rate R i can be expressed as follows [6]: R i = C 2 = ρ s n ρ l n = C max n G, (7) where C max is the maximum concentration of dissolved metal in solution and n and n are instant and initial porosity. Assuming that the absorption on pore surface is linear, balanced, and thermal, the relationship between dissolved term and absorption term is σ + p=, σ = Dε, (2) s= ds dc C =k f C. (8) where σ is stress matrix; ε is strain matrix; D, theelasticity matrix, is a functionof Young s moduluse and Poisson s ratio ]. Consider σ =[σ xx σ yy σ xy ] T, ε =[ε xx ε yy ε xy ] T, D = E ( υ) (+υ)( 2υ) υ υ υ [ υ. 2υ ] [ 2 ( υ)] With S being the displacement vector, strain matrix ε and volumetric strain ε V canbeexpressedasfollows: ε xy = 2 ( S x y + S y x ), S =[S x S y ] T, ε xx = S x x, ε V =ε xx +ε yy. ε yy = S y y, 2.2. Mass Transfer. Both H + of solvent and Cu 2+ of solute are transported by the leaching solution. The couple mass relationship is deduced based on the continuous reaction rates between them. The equations describing mass transfer in pore during leaching processes are C C 2 (3) (4) + 2 s b +u C C y D 2 y 2 + 2J d b = βr i, (5) + 2 s b +u C 2 C 2 y D 2 y 2 + 2J d b =R i, (6) That is, s =k C f, (9) where k f is distributed coefficient [7]. Considering the diffusion flux J d, according to the first Fick theorem, J d = nd C x. () Substituting (6) and (7) into (3) and (4) introduces the retardation coefficient R. Consider R=+ 2 b k f. () The solute transmission equations can be written as follows: C C 2 + u C R y D 2 C R y 2 + u C 2 R y D 2 C 2 R y 2 3. Numerical Analysis 2nD C br x = βc max n G, 2nD C 2 br x = C max n G. (2) Regarding a leaching column with m height and solution concentration unit being continually supplied from the top of the column for days, application rate is w =.25 6 m 3 /(m 2 s).thecalculatedmodelisillustratedin Figure.

3 Mathematical Problems in Engineering 3 Table : The chemical analysis of main element contained in ore sample. Component Cu Fe S Mo SiO 2 Al 2 O 3 CaO MgO As Percentage (%) y Barren solution C n = C 2 n = C n = C 2 n = Volumetric strain Pregnant solution Figure : Schematic of numerical calculation. The initial conditions, top boundary conditions, and bottom boundary conditions for the flow, deformation, and mass transfer coupled equations are p(y,)=, n [ κ η ( p + ρ fg e)] (, t) =.25 6, p (, t) =, S (y, ) =, S (, t) free, S (, t) =, C (y, ) =, C 2 (y, ) =, C (, t) =, C 2 (, t) =, n [θd C (, t)] =, n [θd C 2 (, t)] =. (3) During calculation process, initial porosity n and final porosity n f areassumedtobe.3and.35,respectively; the stoichiometric coefficient β is. Equations (), (2), and (2) are solved by Comsol Multiphysics Software for the given problem. Figure 2 shows the variations of the volumetric strain in leaching column with respect to the bed s depth at different time intervals. It indicates the volumetric strain decreases with the bed s depth increasing. This is because reagent reacts with valuable metal and consumes gradually. Figure 3 shows the spatial and temporal distributions of dissolved mineral and reagent at different time. (a) indicates Time (d) Heap height y (m) Figure 2: Distribution of volumetric strain in leaching column. the solute concentration increases rapidly at the first stage and reaches the peak value and decreases gradually towards the heapbottom.thepeakvaluesmoverightwardswithleaching duration. The reason is that, at the beginning of leaching, solvent concentration is higher and chemical reaction speed isquicker.moreover,thecontentoftargetmetalsisalso higher. (b) indicates the concentration of solvent decreases with the depth increasing which is because chemical reaction consumes reagent. 4. Experiment and Discussion To verify the numerical simulations, indoors physical experiment is done according to dump leaching in Dexing copper mine, Jiangxi province. The chemical content analysis of ore sample is.2% sulphide copper,.7% sulphide copper,.2% free oxide copper, and.72% combined oxide copper. Ore component analysis is conducted by X-Ray Diffractometer M2X and is shown in Table. The maximum diameter of ore particle in dump leaching field in Dexing mine is 8 mm. It is very difficult to carry on experiment according to field situation. What is more, the general apparatus is not large enough to hold such large ore sample, so most theoretical research works are conducted

4 4 Mathematical Problems in Engineering Table 2: The distribution of ore particle diameter after crashing. Particle diameter (mm) < Content (%) Solute concentration (kg/m 3 ) Reagent concentration (kg/m 3 ) Heap height y (m) Heap height y (m) Time (d) Time (d) (a) Solute concentration distribution (b) Reagent concentration distribution Figure 3: Spatial and temporal distributions of reagent (a) and dissolved mineral (b). indoors. The inner diameter of the column leaching cylinder used in experiment is 5 mm; it is necessary to crash ore sample to let the diameter be less than mm according to the research conclusions obtained by Bear [3]. The distribution of ore particle diameter after crashing is shown in Table 2. The samples were bioleached in PVC (5 cm in diameter and cm in height) for days. Solution with a concentration of unit is continually supplied from the top of the column; the application rate is w =.25 6 m 3 /(m 2 s). As shown in Figure 4, numerical results and experimental values of copper ion concentration at a certain point (nearly themiddlepartofthetrunk)areconsistentonthewhole trend, which indicates that the mathematical model, the numerical method, and parameters can describe the transmission process in leaching ore column. 5. Conclusions (i) With respect to the mineral skeleton deformation, a flow and solid elastic model is developed to describe Concentration (mol/l) Calculated values Experimental values Time (min) Figure 4: Comparison between calculated result and experimental result of the concentration of copper ions in ore heap.

5 Mathematical Problems in Engineering 5 theflowreactionandmasstransferprocessesinheap leaching. (ii) The model equations are solved by Comsol Multiphysics Software. The distributions of volumetric strain and concentrations of reagent and dissolved mineral are given based on numerical results. (iii) The numerical simulation results show that the straight strain decreases with the bed s depth increasing; the concentration of the solvent decreases with the bed s depth increasing; the concentration of dissolved mineral increases firstly and decreases from a certain position: the peak values of the curves move leftward with time. (iv) The numerical results are compared with the experimental results; these two are in agreement on the wholetrend,whichindicatesthatthemathematical model, the numerical method, and parameters can describe the multifactor coupled processes in heap leaching. [9] E.Cariaga,F.Concha,andM.Sepúlveda, Flow through porous media with applications to heap leaching of copper ores, Chemical Engineering Journal, vol., no. 2-3, pp. 5 65, 25. [] G. A. Sheikhzadeh, M. A. Mehrabian, S. H. Mansouri, and A. Sarrafi, Computational modelling of unsaturated flow of liquid in heap leaching using the results of column tests to calibrate the model, International Heat and Mass Transfer,vol. 48,no.2,pp ,25. [] A. X. Wu, J. Z. Liu, and S. H. Yin, Mathematical model and analytic solution of mass transfer in heap leaching process, Mining and Metallurgical Engineering,vol.5,pp.7,25. [2] F. A. L. Dullien, Porous Media: Fluid Transport and Pore Structure, Academic Press, New York, NY, USA, 979. [3] J. Bear, Dynamics of Fluids in Porous Media,translatedbyJ.S.Li & Z. X. Chen, China Architecture & Building Press, 983. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgment This project is supported by the Natural Science Fund of China (54, 53476, and 5743). References [] Z.-K. Li, B.-W. Gui, and X.-X. Duan, The production practice and technical study in heap leaching mill of Dexing Copper Mine, Mining and Metallurgical Engineering, vol. 22, no., p. 46, 22 (Chinese). [2] H. Z. Liu and Y. K. Zhang, The status and prospect of intractable low-grade gold mineral bacterium heap leaching, Sichuan Geology, vol. 8, no. 3, pp , 998 (Chinese). [3] T. Yuan, G. Zibin, and G. Renxi, The advance of the uranium andgoldoreheapleaching, Uranium Metallurgy,vol.7,no.2, pp. 2 26, 998 (Chinese). [4] R. W. Bartlett, Solution Mining: Leaching and Fluid Recovery of Minerals, Gordon and Breach Science Publisher, Singapore, 992. [5] S. C. Bouffard and D. G. Dixon, Investigative study into the hydrodynamics of heap leaching processes, Metallurgical and Materials Transactions B,vol.3,pp ,2. [6] A. C. Lasaga, Chemical kinetics of water-rock interations, Geophysical Research, vol.89,no.6,pp , 982. [7] G. E. Grisak and J. A. Cherry, Solute transport through fractured media: 2. Column study of fractured till, Water Resources Research, vol. 6, no. 4, pp. 4 45, 998. [8] M. Sidborn, J. Casas, J. Martínez, and L. Moreno, Twodimensional dynamic model of a copper sulphide ore bed, Hydrometallurgy,vol.7,no.-2,pp.67 74,23.

6 Advances in Operations Research Advances in Decision Sciences Applied Mathematics Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Discrete Mathematics Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis International Stochastic Analysis Optimization