TY - JOUR
T1 - Kinetics of spodumene calcination (α-LiAlSi2O6)
AU - Abdullah, Arif A.
AU - Dlugogorski, Bogdan Z.
AU - Oskierski, Hans C.
AU - Senanayake, Gamini
N1 - Publisher Copyright:
© 2024 The Authors
PY - 2024/9/15
Y1 - 2024/9/15
N2 - Kilns consume about half of the energy necessary to operate lithium refineries and their decarbonisation requires accurate modelling of the calcination of spodumene concentrates fed to the process. This contribution applies the isoconversional methodology to investigate the kinetic parameters of the transition of α-spodumene to its high-temperature polymorph of β-spodumene, using the heat flux measurements from the differential scanning calorimetry (DSC). The normalised energy demand (αH), presented as a function of temperature, characterises these measurements. The activation energy Eα and the product of the reaction model fα and the frequency factor Aα, (fA)α, depend on αH. As the process involves multi-step reactions, we deploy the Friedman differential method and the accurate flexible integral method of Vyazovkin to obtain the kinetic parameters. We also modify the method of Ortega to acquire additional estimates of Eα and (fA)α and apply the rigid integral method of Starink for comparison. The Friedman, Vyazovkin and modified methods deliver the same estimates of the kinetic parameters within their error bands. The Starink method works surprisingly well for predicting the conversion time despite the inaccuracies in the derived values of Eα and (fA)α. This comes to pass because of the compensation effect between these parameters. The activation energy declines rapidly from around 1000 kJ mol−1 at the commencement of the heat treatment to 668 kJ mol−1 at αH = 0.22, then, decreases gently to 577 kJ mol−1 at αH = 0.98, during successive recrystallisation events. Average uncertainties in these results amount to 13 kJ mol−1. The frequency factors fall between 58.5 (±1.0) min−1 and 51.0 (±3.2) min−1, as computed at αH = 0.23 and 0.98, respectively. The so-called false compensation analysis reveals that the first-order reaction model (in αH) governs the energy demand for calcination for αH ≥ 0.23, but, initially, the transformation proceeds through the dissociation-diffusion regime that is not part of the established reaction models. This regime must not be ignored in modelling the calcination of α-spodumene, as it consumes around 20 % of energy required for the transformation reactions. The results reveal significant differences in the predictions of the treatment time, by more than two orders of magnitude, from the existing kinetic models, and explain the differences by the experimental conditions to collect the data for the models. The dissociation of Si-O bonds and diffusion of Si4+ ions out of their tetrahedral cages govern the onset of the thermal treatment of α-spodumene and account for the elevated values of the activation energy in the dissociation-diffusion regime. The two recrystallisation events are limited by the multicomponent diffusion, especially Si4+, in partly crystallised structure. The recrystallisation of γ- to β-spodumene defines the required retention time of concentrate particles in the kiln, to maximise the effectiveness of the subsequent recovery of lithium from the treated material. Fitting Eα and ln(fA)α to polynomials allows a convenient integration of the model equation for making predictions under any heating program. The model forecasts well the transformation of particles of α-spodumene characterised by d80 = 315 µm, studied in additional experiments, using the X-ray powder diffraction to quantify the conversion of α-spodumene. The predictions from the isoconversional model also concur well with other conversion measurements available in the literature, within the expected variability of different spodumene concentrates.
AB - Kilns consume about half of the energy necessary to operate lithium refineries and their decarbonisation requires accurate modelling of the calcination of spodumene concentrates fed to the process. This contribution applies the isoconversional methodology to investigate the kinetic parameters of the transition of α-spodumene to its high-temperature polymorph of β-spodumene, using the heat flux measurements from the differential scanning calorimetry (DSC). The normalised energy demand (αH), presented as a function of temperature, characterises these measurements. The activation energy Eα and the product of the reaction model fα and the frequency factor Aα, (fA)α, depend on αH. As the process involves multi-step reactions, we deploy the Friedman differential method and the accurate flexible integral method of Vyazovkin to obtain the kinetic parameters. We also modify the method of Ortega to acquire additional estimates of Eα and (fA)α and apply the rigid integral method of Starink for comparison. The Friedman, Vyazovkin and modified methods deliver the same estimates of the kinetic parameters within their error bands. The Starink method works surprisingly well for predicting the conversion time despite the inaccuracies in the derived values of Eα and (fA)α. This comes to pass because of the compensation effect between these parameters. The activation energy declines rapidly from around 1000 kJ mol−1 at the commencement of the heat treatment to 668 kJ mol−1 at αH = 0.22, then, decreases gently to 577 kJ mol−1 at αH = 0.98, during successive recrystallisation events. Average uncertainties in these results amount to 13 kJ mol−1. The frequency factors fall between 58.5 (±1.0) min−1 and 51.0 (±3.2) min−1, as computed at αH = 0.23 and 0.98, respectively. The so-called false compensation analysis reveals that the first-order reaction model (in αH) governs the energy demand for calcination for αH ≥ 0.23, but, initially, the transformation proceeds through the dissociation-diffusion regime that is not part of the established reaction models. This regime must not be ignored in modelling the calcination of α-spodumene, as it consumes around 20 % of energy required for the transformation reactions. The results reveal significant differences in the predictions of the treatment time, by more than two orders of magnitude, from the existing kinetic models, and explain the differences by the experimental conditions to collect the data for the models. The dissociation of Si-O bonds and diffusion of Si4+ ions out of their tetrahedral cages govern the onset of the thermal treatment of α-spodumene and account for the elevated values of the activation energy in the dissociation-diffusion regime. The two recrystallisation events are limited by the multicomponent diffusion, especially Si4+, in partly crystallised structure. The recrystallisation of γ- to β-spodumene defines the required retention time of concentrate particles in the kiln, to maximise the effectiveness of the subsequent recovery of lithium from the treated material. Fitting Eα and ln(fA)α to polynomials allows a convenient integration of the model equation for making predictions under any heating program. The model forecasts well the transformation of particles of α-spodumene characterised by d80 = 315 µm, studied in additional experiments, using the X-ray powder diffraction to quantify the conversion of α-spodumene. The predictions from the isoconversional model also concur well with other conversion measurements available in the literature, within the expected variability of different spodumene concentrates.
KW - Decrepitation of spodumene
KW - Kinetics of spodumene activation
KW - Lithium extraction and refining
KW - Lithium hydroxide and carbonate
KW - Phase change
KW - Reconstructive transformation
UR - http://www.scopus.com/inward/record.url?scp=85200959306&partnerID=8YFLogxK
U2 - 10.1016/j.mineng.2024.108902
DO - 10.1016/j.mineng.2024.108902
M3 - Article
AN - SCOPUS:85200959306
SN - 0892-6875
VL - 216
SP - 1
EP - 24
JO - Minerals Engineering
JF - Minerals Engineering
M1 - 108902
ER -