Unsaturated flow and solute transport in a porous column using spherical ore particles

Abstract

This paper dealt with the development of a two-dimensional (2D) mathematical model for column leaching and confirmed the important simulation parameters through experiment. The unsaturated state of the variably saturated flow column and the solute transport of copper ions were studied during leaching. The fluid flow problem was handled using the Richards equation on the premise of an ambient pressure column air, where the van Genuchten formulas were applied to define the nonlinear relationships of pressure head with the retention and permeability properties. The ore column permeability test gave a varied hydraulic conductivity, which was analyzed in the model. In the solute transport problem, the copper ion concentration was solved using the advection-diffusion-reaction equation whose reaction term was determined by the joint analysis of experimental copper leaching rate and the shrinking core model. Particle- and column-scale leaching tests were carried out to illustrate the difference and connection of copper extraction in both processes. This fluid flow and solute transport coupled model was determined through the finite element method using the numerical simulation software, COMSOL Multiphysics.

Keywords

mathematical models / leaching / unsaturated flow / solute transport

References

[1] Petersen J, Watling H, Dixon D, Franzmann P, Plumb J, Johnson J, Harrison S, Hansford G. Menacho J, Casas J. Progress on the development of comprehensive understanding and a model of copper heap bioleaching: The Amira P768 Project. III International Copper Hydrometallurgy Workshop, 2005, Chile, Santiago, 333

[2] Rodríguez Y, Ballester A, Blázquez ML, González F, Muñoz JA. Study of bacterial attachment during the bioleaching of pyrite, chalcopyrite, and sphalerite. Geomicrobiol. J., 2003, 20, 131.

[3] Rohwerder T, Gehrke T, Kinzler K, Sand W. Bioleaching review part A: Progress in bioleaching: fundamentals and mechanisms of bacterial metal sulfide oxidation. Appl. Microbiol. Biotechnol., 2003, 63, 239.

[4] Naderi H, Abdollahy M, Mostoufi N, Koleini MJ, Shojaosadati SA, Manafi Z. Kinetics of chemical leaching of chalcopyrite from low grade copper ore: behavior of different size fractions. Int. J. Miner. Metall. Mater., 2011, 18(6): 638.

[5] Halinen AK, Rahunen N, Kaksonen AH, Puhakka JA. Heap bioleaching of a complex sulfide ore: Part I: Effect of pH on metal extraction and microbial composition in pH controlled columns. Hydrometallurgy, 2009, 98, 92.

[6] Darezereshki E, Schaffie M, Lotfalian M, Seiedbaghery SA, Ranjbar M. Use of mesophilic and thermophilic bacteria for the improvement of copper extraction from a low-grade ore. Int. J. Miner. Metall. Mater., 2011, 18(2): 138.

[7] Cariaga E, Concha F, Sepúlveda M. Flow through porous media with applications to heap leaching of copper ores. Chem. Eng. J., 2005, 111, 151.

[8] McBride D, Cross M, Croft N, Bennett C, Gebhardt J. Computational modelling of variably saturated flow in porous media with complex three-dimensional geometries. Int. J. Numer. Methods Fluids, 2006, 50, 1085.

[9] Fick A. On liquid diffusion. Philos. Mag. J. Sci., 1855, 10, 30

[10] Levenspiel and Octave, Chemical Reaction Engineering, 2nd Ed., John Wiley & Sons, Inc., New York, 1972, p. 361.

[11] van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J., 1980, 44, 892.

[12] McBride D, Cross M, Gebhardt JE. Heap leach modeling employing CFD technology: a ‘process’ heap model. Miner. Eng., 2012, 33, 72.

[13] Bennett CR, McBride D, Cross M, Gebhardt JE. A comprehensive model for copper sulphide heap leaching: Part 1. Basic formulation and validation through column test simulation. Hydrometallurgy, 2012, 127–128, 150.

[14] Ogbonna N, Petersen J, Laurie H. An agglomerate scale model for the heap bioleaching of chalcocite. J. S. Afr. Inst. Min. Metall., 2006, 106, 433

[15] Nativ M, Goldstein S, Schmuckler G. Kinetics of ion-exchange processes accompanied by chemical reactions. J. Inorg. Nucl. Chem., 1975, 37, 1951.

[16] Miao XX, Liu JZ, Wu AX. Experimental study on copper leaching by square column and its numerical analysis. Min. Metall. Eng., 2012, 32, 74

[17] Sidborn M, Casas J, Martínez J, Moreno L. Two-dimensional dynamic model of a copper sulphide ore bed. Hydrometallurgy, 2003, 71, 67.

[18] Bear J. Hydraulics of Groundwater, 1979, New York, McGraw-Hill Inc., 231

[19] Leahy MJ, Davidson MR, Schwarz MP. A two-dimensional CFD model for heap bioleaching of chalcocite. ANZIAM J., 2005, 46, C439.

[20] Bischoff KB. Accuracy of the pseudo steady state approximation for moving boundary diffusion problems. Chem. Eng. Sci., 1963, 18, 711.

[21] Liddell KC. Shrinking core models in hydrometallurgy: what students are not being told about the pseudo-steady approximation. Hydrometallurgy, 2005, 79, 62.

[22] Lizama HM, Fairweather MJ, Dai Z, Allegretto TD. How does bioleaching start?. Hydrometallurgy, 2003, 69, 109.

[23] Decker DL, Tyler SW. Evaluation of flow and solute transport parameters for heap leach recovery materials. J. Environ. Qual., 1998, 28, 543.