SciELO - Scientific Electronic Library Online

 
vol.45Bienvenido el nuevo año 2025Modelación de la transferencia de calor convectiva en tubos verticales e inclinados bajo efectos viscoso-gravitacionales índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Articulo

Indicadores

  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Tecnología Química

versión On-line ISSN 2224-6185

RTQ vol.45  Santiago de Cuba  2025  Epub 25-Feb-2025

 

Artículo original

Lixiviación carbonato amoniacal: Simulación del perfil de amoniaco y concentración de sólidos a escala industrial por balance de masa

Ammoniacal Carbonate Leaching: Simulation of ammonia profile and solid concentration on a full-scale by mass balance

0000-0002-8927-2023Armando Rojas-Vargas1  * 

1Empresa de Servicios Técnicos de Computación, Comunicaciones y Electrónica del Níquel (SERCONI), Cuba

RESUMEN

La modelación y simulación son útiles para evaluar, analizar o entrenar un sistema real o conceptos representados por un modelo. Esta investigación se realizó en la tecnología de lixiviación amoniacal para la recuperación de níquel del mineral laterita, con el propósito de desarrollar un algoritmo de simulación del perfil de amoniaco y la concentración de sólidos por balance de masa en régimen estacionario. El diagrama de flujo del sistema de lixiviación asumido como modelo a escala industrial, consistió en reactores continuos aireados con agitación y sedimentadores dispuestos en ocho etapas de lixiviación y lavado, en configuración contracorriente. Como resultado, se determinaron los flujos y puntos de alimentación del solvente para alcanzar la mayor aproximación de los parámetros de operación a los valores normados. En posteriores trabajos se va a incorporar al algoritmo la función de distribución del tiempo de residencia y la modelación cinética con vistas a simular la eficiencia de lixiviación de níquel.

Palabras-clave: lixiviación; balance de masa; níquel; simulación

ABSTRACT

Modeling and simulation are useful for testing, analysis or training a real system or concept represented by a model. This research was carried out in the ammoniacal leaching technology for nickel recovery from laterite ore, with the aim of performer a simulation algorithm of ammonia profile and solid concentration by mass balance in steady state. Flowsheet of leaching system assumed as model on a full-scale, consisted of aerated continues stirring tank reactors and thickeners tanks disposed in eight steps of leaching and washing, in countercurrent operation. As a result, the solvent flows and feed points were determined to achieve a closer approximation of operating parameters to standardized values. In future work, the residence time distribution function and kinetic modeling will be incorporated into the algorithm to simulate the efficiency of nickel leaching.

Key words: leaching; mass balance; nickel; simulation

Introduction

Modeling and simulation (M&S) provides virtual duplication of a system entity, process or phenomenon, and represents those in readily available and operationally valid environments. A model is a physical, mathematical or logical representation, while a simulation is the implementation of a model over time.1,2

M&S are useful for testing, analysis or training a real system or concept represented by a model.

Systems are modeled as having a number of inputs x1, x2…, xn; outputs or responses y1, y2... yn; and system parameters P1, P2…, Pn.

An essential stage in the development of any model is the formulation of mass and energy balance equations. In addition, kinetic reaction, rate of transfer of momentum, heat and mass, and property changes physic-chemistry and thermodynamic are considered.

Models can be static, representing a system at a particular point in time, or dynamic, representing how a system changes with time.

Several numerical or experimental works deal with M&S in the ammoniacal carbonate leaching technology have been presented, in the process steps: laterite ore beneficiation, drying 3,4,5, grinding 6, reduction-roasting 7,8, leaching of roaster product in ammoniacal solvent, precipitation of nickel as basic nickel carbonate (BNC) intermediate product in steam-stripped columns 9,10, and thermal decomposition to NiO. 11,12

In addition, M&S of thickener tanks have been made by non-linear partial differential equations, from both mass and momentum balance equations, to attain a better knowledge of the decantation process.13,14 Likewise, combination of Residence Time Distribution (RTD) function and kinetic modelling is used to predict the course of chemical reactions in Continuous Stirring Tank Reactors (CSTR).15

In this work, M&S was performed in the leaching system on a full-scale by total mass balances and component mass balances at steady-state, in order to simulate the ammonia profile and solids concentration in the system, and analyze deviations from operating parameters. Flowsheet assumed as model consisted of aerated CSTR and thickeners tanks disposed in eight steps of leaching and washing, in countercurrent operation.

Materials and methods

Flowsheet of leaching system on a full-scale

Roaster product is mixed with ammonia/ammonium carbonate solution in the quench tank at a liquid-to-solid ratio 5.4-6.3 kg/kg. The resulting slurry is leached in aerated CSTR. The Pregnant Liquor is separated from the undissolved solid in thickener tanks. The ammoniacal liquors (Fresh Liquor and Weak Liquor) are feed in countercurrent configuration from washing to leaching steps (figure 1).

Fig. 1  Flowsheet of leaching system on a full-scale Product Liquor CNi: 8.0 g/L 

Input parameters of simulation algorithm

Leaching conditions are based in the following independent variables: (10,16

Roaster product flow (t/h) > 100 t/h

Product Liquor: ammonia > 60 g/L, nickel > 6.5 g/L

Feed sludge density in first leaching stage: 1142-1160 kg/m3

Thickened sludge density: 1700-1750 kg/m3 for leaching-steps and 1800-1850 kg/m3 for washing-steps

Nickel extraction for each step (%)

Ammonia loses in each step (g/L)

Ammonia concentration: 120-130 g/L for Fresh Liquor, and 10-15 g/L Weak Liquor

Density real of roaster product: 4200 kg/m3 for leaching-steps and 4000 kg/m3 for washing-steps

Product Liquor density: 1024 kg/m3 for leaching and 1000 kg/m3 for washing-steps

Output parameters of simulation algorithm

The output parameters of simulation algorithm are:

Fresh Liquor flow, second washing stage

Weak Liquor flow, third washing stage

Density profile in the feed sludge

Ammonia concentration profile

Flow of thickened sludge (underflow) and clear liquor (overflow)

Ammonia losses, by the “Black box model”

Air flow and theoretical retention time in the leaching reactors

Makeup Fresh Liquor flow, in the second and third leaching stage and first washing stage

Equations for the mass balance

Flowsheet of thickener tanks are shown in figure 2 and the parameters: Sludge Feed (F, m3/h), overflow or clear liquor (O, m3/h), underflow or thickened sludge (U, m3/h), density (ρ, kg/m3) and solid concentration (C, kg/m3).

Fig. 2 Flowsheet of thickener tanks 

Mass balance for liquor and solid are represented by (1-2)

Concentration of solid is determined by (4), liquid-to-solid ratio (L/S, kg/kg) (5), where sub-indexes “p, s, l” represent particle, sludge and liquor, respectively.

Feed sludge density profile in the thickener tanks was determined fallowing the algorithm from figure 3, were roasted product flow is RP (t/h). Sludge density is supposed (II) to check there is no solid accumulation (V). Then, a component mass balance was performed to determine ammonia concentration for each stage.6 Finally, both ammonia concentration and feed density are supposed for the last step of washing (T8A and T8B); using figure 3, the solid accumulation is checked, and then ammonia concentration is tested (7).

Fig. 3 Algorithm for density profile in the thickeners feed 

Restrictions and validation of simulation algorithm

Mass balance was performed in steady state. The ammonia concentration in the third washing step is approximately 30 g/L.16 Sludge density in the thickeners feed should as following: First-step of leaching 1142-1160 kg/m3, second-step of leaching 1200-1250 kg/m3, from third to fifth step of leaching 1300-1350 kg/m3, and from first to third step of washing: 1350-1450 kg/m3.

The simulation algorithm was validated by the “Black-box model” for the flows of sludge 8, solid and the component ammonia, the interpretation of simulation results and industrial operation.

Results and discussion

Mass balance was performed at the operating parameters showed in table 1, simulating different ammonia gradient (∇NH3) in the reactors of first leaching step.

Table 1 Input/output parameter of simulation algorithm (CNi 8.0 g/L)  

From table 1, the higher the ammonia gradient in first-step leach reactors is, the lower Weak Liquor flow and the higher the makeup Fresh Liquor flow are. No Fresh Liquor was fed to thickener tanks T7A and T7B because an ammonia concentration is obtained in the third washing step (T8A, T8B) much higher than the recommended value (30 g/L). This avoids greater energy consumption in the steam stripping column for the ammonia recovery from washed solid residue (figure 4).

Furthermore, when ammonia gradient increases the slurry density feed to thickeners tanks also tend to increase. For ∇NH3 = 5.5 g/L the density is higher than the target values. Therefore, operation of the leaching system at these flow and concentration conditions could result in a higher concentration of suspended solids in the overflows and a lower settling velocity, which is determined experimentally. Dirty thickener overflows affect the efficiency of the ammoniacal liquor cooling system, by increasing deposits accumulation onto the internal surfaces of the heat exchangers. (17 The higher the temperature in the leaching system the greater the evaporation of ammonia, thus generating a chain of disturbances that tend to reduce leaching efficiency (figure 5).

Fig. 4 Simulation of density profile of leaching process on a full-scale at different ammonia gradients in the reactors of leaching first-step, CNi 8 g/L (table 1

Fig. 5 Simulation of ammonia profile of leaching process on a full-scale at different ammonia gradients in the reactors of leaching first-step, CNi 8 g/L (table 1

In addition, a mass balance was performed at a nickel concentration of 6.5 g/L in the Product Liquor (PL), at the same input parameters (table 2) (figure 6).

Table 2 Input/output parameter of simulation algorithm (CNi 6.5 g/L) 

As shown in table 2, the lower the concentration of nickel in the Product Liquor, the flow of Product Liquor, Weak Liquor, Fresh Liquor, and the second stage overflow liquor flow (313 m3/h) tend to increase. Then, an increase in the liquor velocity can be expected, with a greater possibility of entraining fine particles.

Density profile of fed sludge to the thickener tanks decreased, thus achieving a density closer to operating standards (figure 7).

Ammonia profile was approximately the same compared to the simulation at CNi 8.0 g/L, for the same NH3 concentration gradient (figure 8).

The algorithm was tested using the “Black Box Model”, and as a result the flows of sludge, solids and ammonia in the inlet and outlet streams of the system presented an absolute error in the estimate of less than 0,5 %.

The results contribute to decision-making regarding the feeding point of ammoniacal solvent (Fresh Liquor and Weak Liquor) to achieve a lower deviation of operating parameters with respect to standardized values in the leaching process. Future work will incorporate the residence time distribution function and kinetic modeling into the algorithm to simulate nickel leaching efficiency.18 Furthermore, it is necessary to perform an experimental evaluation of the simulation results.

Fig. 6  Flowsheet of leaching system on a full-scale Product Liquor CNi: 6,5 g/L 

Fig. 7 Simulation of density profile of leaching process on a full-scale at different ammonia gradients in the reactors of leaching first-step, CNi 6,5 g/L (table 2

Fig. 8 Simulation of ammonia profile of leaching process on a full-scale at different ammonia gradients in the reactors of leaching first-step, CNi 6,5 g/L (table 2

Conclusions

  1. A numeric algorithm was performed by mass balance in steady state to model ammonia profile and solid concentration in the nickel leaching process on a full-scale, which was validated by the “Black Box model”, the interpretation of simulation results and the industrial operation.

  2. Simulation of operating parameters was carried out at different ammonia gradients in the leaching reactors, then the nickel concentration in the Product Liquor must be considered to determine solvent flows and feed points to achieve a closer approximation to standardized values.

Nomenclature

References

1.  MENNER, W. A. Introduction to Modeling and Simulation. Johns Hopkins APL Technical Digest1.  . 1995. 16(1). [ Links ]

2.  Maria, A. Introduction to modeling and simulation. In Proceedings of the 29th conference on Winter simulation. 1997. 7-13. https://doi.org/10.1145/268437.268440 [ Links ]

3.  BELETTE, O-F. Comportamiento del secado solar de lateritas a la intemperie en condiciones de clima húmedo tropical. HOLOS3.  . [online]. 2020, 7. 1-14. https://doi.org/10.15628/holos.2020.9767 [ Links ]

4.  ROJAS, V.A., PÉREZ, G.L, SÁNCHEZ, G.C., et al. Performance evaluation of a flighted rotary dryer for lateritic ore in concurrent configuration. Heliyon 9.4.  [online]. 2023. e21345. https://doi.org/10.1016/j.heliyon.2023.e21345 [ Links ]

5.  ZALAZAR, C.O, GÓNGORA, E.L., LEGRÁ, A.A.L, et al. Modelación matemática del secado de menas lateríticas de Ni en un horno cilíndrico rotatorio. Minería y Geología5.  . [online]. 2023, 38(4). 346-362. ISSN 0258 5979 [ Links ]

6.  ANGULO, H.J.P, MANUEL, R.N.G., LEGRÁ, A.L.L., et al. Modelación de la distribución de tamaño de menas saprolíticas alimentadas al proceso Caron.Tecnología Química6.  [online]. 2024, 44(1). 91-106. e-ISSN: 2224-6185. [ Links ]

7.  CHEN, S., GUO, S., JIANG, L., XU, Y., DING, W. Thermodynamic of selective reduction of laterite ore. Trans. Nonferrous Met. Soc. China7.  . [online]. 2015. 25. 3133−3138. https://doi.org/10.1016/S1003-6326(15)63943-7 [ Links ]

8.  RHAMDHANI, M.A., HAYES, P.C., JAK, E. Ni-laterites Part II - Thermodynamic analysis of phase transformations. Trans. Inst. Min. Metall. C.8.  [online]. 2009. 118(3).146-155. https://doi.org/10.1179/174328509X431409 [ Links ]

9.  ROJAS, V.A., MAGAÑA, H.M.E., RICARDO, A.R. Lixiviación carbonato amoniacal: Estimación del niquel disuelto en el efluente de destilación. Revista de Metalurgia9.  [online]. 2019. 55(3). 2-13. https://doi.org/10.3989/revmetalm.1499.  [ Links ]

10.  ROJAS, V.A., NIEVES, M.E.T., DÍAZ, Y.G. Ammoniacal Carbonate Leaching: effect of dissolved sulfur in the distillation operation, Acta Chim. Slov10.  . [online]. 2020. 67. 1239-1249. https://doi.org/10.17344/acsi.2020.6147 [ Links ]

11.  GAINZA, Y.D., MAGAÑA, M.E.M., ROJAS, A.R. Determinación del índice de consumo de petróleo del horno de calcinación de la empresa ECG. Tecnología Química11.  . [online]. 2006. 26(3). 407-416. e-ISSN: 2224-6185. [ Links ]

12.  COLUMBIE, A.O.N, RODRÍGUEZ, J.G., GUZMÁN, D.R.G., et al12.  . Modelo matemático del proceso de calcinación del carbonato básico de níquel en el horno tubular. Minería y Geología12.  . [online]. 2000. 17(2). ISSN 0258 5979 [ Links ]

13.  VALENTIN, C., DOCHAIN, D., JALLUT, C., DOS SANTOS MARTINS, V. Representation of a Continuous Settling Tank by Hybrid Partial Differential Non Linear Equations for Control Design. World congress IFAC. 2020, Berlin, Germany. https://hal.science/hal-0237264313.  [ Links ]

14.  BAKIRI, Z., CHEBLI, D., NACEF, S. Dynamic modelling of the secondary settler of a wastewater treatment via activated sludge to low-load. Energy Procedia14.  [online]. 2012. 18, p. 1-9. ISSN: 1876 6102 [ Links ]

15.  Levenspiel, O. Chemical reaction engineering (3a ed.), Limusa Wiley, Mexico, 2004, p. 678. ISBN 0-471-25424-X [ Links ]

16.  CHANG, C.A.R., ROJAS, V.A. La lixiviación del proceso Caron: síntesis del conocimiento para su perfeccionamiento industrial. Parte I. Tecnología Química16.  . [online]. 2009. 29(1). 98-107. e-ISSN: 2224-6185. [ Links ]

17.  TORRES, E.T., RETIRADO, M.Y., GONGORA, L.E. Coeficientes de transferencia de calor experimental para el enfriamiento de licor. Ingeniería Mecánica17.  . [online]. 2014. 17(1). 68-77. ISSN 1815-5944. [ Links ]

18.  ROJAS, V.A., GARCÍA, A.P Análisis de la curva de distribución del tiempo de residencia en un sistema de lixiviación. Tecnología Química18.  . [online]. 2010. 30(1). 61-68. e-ISSN: 2224-6185. [ Links ]

Recibido: 06 de Septiembre de 2024; Aprobado: 12 de Diciembre de 2024

*Autor para la correspondencia. Correo electrónico: arojas@eros.moa.minem.cu

The author declares that he has not competing interest that could have appeared to influence the work reported in this paper.

Armando Rojas Vargas: conceptualization, ethodology, Investigation, Writing - Reviewing and Editing.

Creative Commons License