### “Simulation of a Hydrometallurgical Leaching Reactor modeled as a DAE system”

**Authors:**Marta Duenas Diez, G. Ausland, Magne Fjeld and Bernt Lie,

**Affiliation:**Telemark University College and Elkem

**Reference:**2002, Vol 23, No 2, pp. 93-115.

**Keywords:**Population Balance Equation, Differential and Algebraic Equations, leaching, dynamic modeling, ferrosilicon

**Abstract:**An existing dynamic model of the main reactor in the Silgrain process for the production of Si from FeSi has been extended here in order to resemble more closely the behavior of the real reactor. The previous model was based on the application of macroscopic mass conservation law, the population balance equation and the assumptions of complete mixing and isothermic conditions. The major modifications are the inclusion of the condition governing the entrainment of particles in the outflow, and the formulation of the energy balance. The extended model consists of 1 integrodifferential equation, 4 implicit ordinary differential equations, 7 algebraic equations and 3 integral equations. After discretization in the particle size space, a system of differential and algebraic equations (DAE) is obtained. DAEs are not ODEs and they require analysis and characterization and may require reformulation. After such analysis, it was concluded that the system is implicit index-one for the usual range of operation and that a method based on the Backward Differentiation formulas (BDF) can be used for its solution. The model was implemented in Matlab and the ode15s code was used for solving the system of equations. The simulation results are satisfactory and seem to match qualitatively with the known operation of the reactor. The model is suitable for further use in designing a model-based control scheme.

PDF (2548 Kb) DOI: 10.4173/mic.2002.2.2

**DOI forward links to this article:**

[1] Marta Dueñas Díez, B. Erik Ydstie, Magne Fjeld and Bernt Lie (2008), doi:10.1016/j.compchemeng.2007.01.007 |

[2] Marta Dueñas Díez, Magne Fjeld, Einar Andersen and Bernt Lie (2006), doi:10.1016/j.ces.2005.01.047 |

[3] Christy M. White, Paul Ege and B. Erik Ydstie (2006), doi:10.1016/j.powtec.2006.01.005 |

[4] V. M Zarochentsev and A. L Rutkovskiy (2024), doi:10.31857/S0869573324021828 |

**References:**

[1] Aas, H. (1971). The silgrain process: Silicon metal from 90, .

[2] Brenan, K. E., Campbell, S. L. Petzold, L. R. (1996). Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations, SIAM Society for Industrial and Applied Mathematics.

[3] Chen, Z., Prüss, J. Warbecke, H. (1998). A population balance model for disperse systems: Drop size distribution in emulsion, Chemical Engineering Science 53(5), 1059-1066 doi:10.1016/S0009-2509(97)00328-X

[4] Coulson, J. Richardson, J. (1978). Chemical Engineering: Fluid Flow, Heat Transfer and Mass Transfer, Vol 1, Pergamon Press, Inc.

[5] Dueñas, M. Lie, B. (2000). Modelling and simulation of a hydrometallurgical leaching reactor, In B. Elmegaard, N. Houbak, A. Jakobsen and F. J. Wagner, eds, 'Proceedings SIMS 2000', SIMS- Scandinavian Simulation Society and Technical University of Denmark, Lyngby.Denmark, pp. 199-225.

[6] Herbst, J. (1979). Rate processes in multiparticle metallurgical systems, In H. Sohn and M.Wadsworth, eds, 'Rate Processes in Extractive Metallurgy', Plenum Press, pp. 53-112.

[7] Herbst, J. Asihene, S. (1993). Modelling and simulation of hydrometallurgical processes using the population balance approach, In V. Papangelakis and G. Demopoulos, eds, 'Modelling, Simulation and Control of Hydrometallurgical Processes. Proceedings of the 32nd Annual Conference of Metallurgists of CIM', Canadian Institute of Mining, metallurgy and petroleum, Quebec, pp. 3-44.

[8] Hulburt, H. M. Katz, S. (1964). Some problems in particle technology - A statistical mechanical formulation, Chemical Engineering Science 19, 555-574 doi:10.1016/0009-2509(64)85047-8

[9] Kapur, P. (1995). Particle population balance in granulation of iron ores by an autolayering mode, In S. Mehrotra and R. Shekhar, eds, 'Mineral Processing. Recent Advances and Future Trends.', Allied Publishers ltd, pp. 703-717.

[10] Kumar, S. Ramkrishna, D. (1996). On the solution of population balance equations by discretization - I: A fixed pivot technique, Chemical Engineering Science 51(8), 1311-1332 doi:10.1016/0009-2509(96)88489-2

[11] Lefkopoulos, A. Stadherr, M. A. (1993). Index analysis of unsteady-state chemical process systems - I: an algorithm for problem formulation, Computers Chem. Engng. 17(4), 399-413 doi:10.1016/0098-1354(93)80031-H

[12] Margarido, F., Martins, J. P., Figueiredo, M. O. Bastos, M. H. (1993). Kinetics of acid leaching refining of an industrial fe-si alloy, Hydrometallurgy 34, 1-11 doi:10.1016/0304-386X(93)90077-Q

[13] Petzold, L. R. (1985). Differential-algebraic equations are not ODEs, SIAM J. Sci. Stat. Compt. .3, 49-95.

[14] Ramkrishna, D. (1985). The status of population balances, Reviews in Chemical Engineering 3, 49-95.

[15] Randolph, A. (1964). A population balance for countable entities, The Canadian Journal of Chemical engineering 42, 280-281.

[16] Randolph, A. D. Larson, M. A. (1988). Theory of Particulate Processes: Analysis and Techniques of Continuous Crystallization, Second edn, Academic Press, Inc.

[17] Rubisov, D. H. Papangelakis, V. G. (1996). Mathematical modelling of the transient behaviour of CSTRs with reactive particulates: Part 1 - The population balance framework, The Canadian Journal of Chemical Engineering 74, 353-362 doi:10.1002/cjce.5450740306

[18] Rubisov, D. H. Papangelakis, V. G. (1996). Mathematical modelling of the transient behaviour of CSTRs with reactive particulates: Part 2 - Application to pyrite pressure oxidation, The Canadian Journal of Chemical Engineering 74, 363-371 doi:10.1002/cjce.5450740307

[19] Rubisov, D. Papangelakis, V. (1995). Model-based analysis of pressure oxidation autoclave behaviour during process upsets, Hydrometallurgy 39, 377-389 doi:10.1016/0304-386X(95)00042-F

[20] Shampine, L. F., Reichelt, M. W. Kierzenka, J. A. (1999). Solving index-1 DAEs in Matlab and Simulink, SIAM Review 4.3, 538-552 doi:10.1137/S003614459933425X

**BibTeX:**

@article{MIC-2002-2-2,

title={{Simulation of a Hydrometallurgical Leaching Reactor modeled as a DAE system}},

author={Duenas Diez, Marta and Ausland, G. and Fjeld, Magne and Lie, Bernt},

journal={Modeling, Identification and Control},

volume={23},

number={2},

pages={93--115},

year={2002},

doi={10.4173/mic.2002.2.2},

publisher={Norwegian Society of Automatic Control}

};