Kinetics of AlOOH dissolving in caustic solution studied by high-pressure DSC
BAO Li, ZHANG Ting-an, DOU Zhi-he, L? Guo-zhi,
GUO Yong-nan, NI Pei-yuan, WU Xu-jian, MA Jia
School of Materials and Metallurgy, Northeastern University, Shenyang 110004, China
Received 17 March 2010; accepted 13 June 2010
Abstract:
Mechanism functions and kinetic parameters of AlOOH (boehmite or diaspore) dissolving in sodium hydroxide solution were researched. The mixture of boehmite or diaspore and caustic solution was scanned by high-pressure differential scanning calorimetry (DSC) instrument with heating rate of 10 °C/min, and differential equation method was used to analyse the DSC curves, combining with iterative method and linear least square method. The most probable mechanism functions for both boehmite or diaspore and caustic solution reactions were logically selected from 30 types of non-isothermal kinetics differential equations, according to the calculated results obtained by Matlab program. The most probable differential mechanism function of boehmite dissolving in caustic solution is f(α)=1-α, which reveals the first-order reaction with apparent activation energy of 79.178 kJ/mol and the preexponential constant 1.031×108 s-1. The function, f(α)=2(1-α)3/2, can describe the dissolution of diaspore sample in sodium hydroxide solution. The calculated results of kinetic parameters are apparent activation energy of 73.858 kJ/mol, preexponential constant of 5.752×107 s-1 and reaction order of 1.5.
Key words:
differential mechanism function; boehmite; diaspore; apparent activation energy; preexponential constant;
1 Introduction
Boehmite, α-AlOOH, and diaspore, γ-AlOOH, which are natural heterogeneous ores and primarily consist of one aluminium hydroxide mineral, are the bauxite used to produce alumina. Bayer process is the most popular method for alumina producing all over the world. Digestion as one of the key steps of the Bayer process aims at dissolving the utmost aluminium available in the ore into caustic solution. The chemical reaction between boehmite or diaspore and sodium hydroxide solution is heterogeneous reaction, the main form of metallurgical reaction can be presented as[1]:
(1)
Despite of many open reports on kinetic modelling of bauxite digestion, there is not a consensus. The basic models, -ln(1-α)=kt and 1-(1-α)1/3=kt were directly employed to investigate the kinetics of leaching gibbsite and low-quality boehmite bauxite with hydrochloric acid by REDDY et al[2] and ZIVKOVSC et al[3]. Based on the shrinking core model, PEREIRA et al[4] presented four models for dissolution of gibbsite in NaOH solution controlled by chemical reaction. Some semi empirical kinetic models were developed by ZAFAR[5] for dissolution of low-grade bauxite ore with sulphuric acid. BRITTAN[6] presented a variable activation energy model postulating that the rate-limiting factors can be lumped together as an overall resistance to reaction equivalent to an activation energy barrier. GU et al[7-8] developed a mathematical model of leaching diaspore on the basis of kinetic steady-state principle, while BI et al[9] presented that the model for leaching diaspore should not be the shrinking unreacted-core model and proposed the models for three leaching temperature stages theoretically. BORNTRAG et al[10] measured the change regulation of solution conductivity to get the kinetic model of gibbsite leaching in alkali solution, while LI et al[11] employed isothermal calorimetric method and YIN et al[12] used non-isothermal analysis with differential scanning calorimetry (DSC) technology. No matter what method the authors above used, all the developed models were on the base of hypothesis simple models.
DSC is a thermoanalytical technique in which the difference in the amount of heat required to increase the temperature of a sample is measured as a function of temperature. The whole reaction process at a certain temperature range could be detected with the help of DSC, and thus it is much convenient to assess the kinetics of reacting process compared with the conventional method. In this work, in order to study the kinetics of dissolving in caustic solution, the high-pressure DSC technique is employed to obtain the non-isothermal kinetics for boehmite and diaspore. The differential equation method is used to analyse the sets of curve data. Finally, the most probable mechanism functions for boehmite and diaspore are selected logically.
2 Experimental
2.1 Preparation of AlOOH
The boehmite sample used was synthesized by hydrothermal synthesis method[13] and the diaspore sample was the diaspore flotation concentrate[14-16].
The material for preparing boehmite was pure gibbsite with the diaspore flotation concentrate as crystal seed, and the reaction media was sodium hydroxide solution with the concentration of 1mol/L. The mixture of materials in proportion reacted in autoclave under the conditions of temperature of 400 °C, stirring speed of 200 r/min and pressure of 17 MPa. Keeping reaction for 25 d, the new sample was washed to be neutral by distilled water, and then dried in oven at 60 °C. Figure 1 shows the XRD pattern of the dried sample, indicating that the sample is boehmite.
Fig.1 XRD pattern of sample prepared by hydrothermal synthesis method and JCPDS pattern of boehmite
The diaspore from Tiandong, Guangxi Province of China, concentrated for 5 times by direct flotation method, was used as the sample. Sodium oleate was taken as the flotation agent, sodium hexamet aphosphate as the regulator, pine oil as the foaming agent and sodium carbonate as the pH regulator. The results of chemical analysis showed that Al2O3 content in the diaspore concentrate was 74%, while the theoretical content in pure diaspore is 85%.
2.2 DSC experiment
The DSC experiment was carried out on 204HP differential scanning analyzer made by NETZSCH Instrument Manufacture Co. Ltd., Germany. A gold crucible with the stainless steel cap was used as the reactor to enduring high pressure. The mixture of solid sample and sodium hydroxide solution was sealed in the gold crucible. Since the leaching condition of boehmite was different from that of diaspore, the caustic concentration of solution for boehmite sample and diaspore sample were 220 and 240 g/L, respectively. The setting parameters of DSC experiment were different as well. For boehmite, the range of scanning temperature was 20-250 °C, while the range for diaspore sample was 20-300 °C. The heating rate was 10 °C/min for both samples. Proteus software developed by NETZSCH was used to collect and process the differential thermal signals.
2.3 Differential equation method
Differential equation method was employed to analyse the DSC curve to research the reaction kinetics. The model of differential equation for non-isothermal kinetics is[17]
where f(α) is the reaction differential mechanism function, α is the reaction fraction, T is the scanning temperature, T0 is the initial reaction temperature on non-isothermal kinetic curve at which DSC curve deviates from the baseline, E is the apparent activation energy of reaction, A is the preexponential constant, β is the linear heating rate and R is the gas constant. Thus, T=T0+βt, where t is the scanning time.
There are 30 mechanism functions of differential equation[17], shown as Table 1.
3 Results and discussion
3.1 Mechanism function of boehmite dissolution
The DSC curve of boehmite and caustic solution scanning from 20 °C to 250 °C is shown in Fig.2. It is seen that an endothermic peak appears at 118 °C. The reaction between boehmite and solution happens around 90 °C, and ends at 180 °C. Besides, the curve between 170 °C and 250 °C tends to a straight line, because the
Table 1 30 mechanism functions of differential equation
reaction rate, α, at 170 °C reaches 99.998%. Therefore, the temperature of starting reaction, T0, is taken as 363.5 K, and the range of T was from 364.5 K to 443.5 K. According to thermal analysis kinetics, the reaction fraction, α, equals the ratio of the heat absorbed at a transient time, t, to the total absorbed heat of reaction. The total absorbed heat is equivalent to the whole peak area above the baseline of DSC curve, H, while the heat absorbed is the part peak area before a transient time, Ht. Proteus software gives a set of values of reaction fraction and the rate of endotherm at time t, dH/dt, at the temperature from 364.5 K to 443.5 K. From α=Ht/H and T=T0+βt, it is easy to express dα/ dT as
Thus, a set of corresponding value of dα/dT can be got.
Fig.2 DSC curve of boehmite and caustic solution scanning from 20 to 250 °C with heating rate of 10 °C/min
Since there is a linear relationship between the left side of Eq.(2) and 1/T, Eq.(2) can be solved by iterative method. Giving E any positive value, the linear least square method could be carried out to obtain new values of E from the slope of straight line and A from the intercept. Considering the new value of E as initial value and performing iterative calculation again, another new value of E would be gained. After several times of iteration, reasonable values of E and A can be got. Using Matlab software[18] to process data, the difference value of two adjacent E less than 10-6 was set as iterating termination condition. If iteration times was more than 500, the corresponding E was set as 0. Each differential mechanism function f(α) in Table 1 is taken into Eq.(2) and processed by the procedures above and the one with the highest correlation coefficient, r, and the least residual variance Q may be chosen as the proper function. The expressions of r and Q are as follows:
where
, L is the number of the data, and Alast and Elast are the values of A and E got by the last iteration, respectively.
Using 30 points on the DSC curve of boehmite sample dissolving in sodium hydroxide solution to analyse the 30 differential functions, the analysis results are shown in Table 2. It is noted that the boehmite or diaspore dissolves in caustic solution is a type of reaction between ions and neutral molecules, shown as Eq.(1). According to transition state theory[19] and numbers of empirical values, the value of preexponential constant A in Arrhenius equation for the reaction between ions and neutral molecules is in the range of 106 to 1013. In this case, only four functions, namely No.1, No.2, No.9 and No.22 may be suitable for the reaction process. However, the function No.9, f(α)=1-α, shows the highest correlation (0.985 6) and the lowest variance (0.237 5) among the four, which thus is considered the most proper expression to describe the kinetics of boehmite dissolving in sodium hydroxide solution. After 8 times of iteration, the obtained apparent activation energy, E, is
Table 2 Kinetic parameters obtained by iterative method of 30 differential mechanism functions using data of DSC curve of boehmite sample dissolving in caustic solution
79.178 kJ/mol and the preexponential constant, A, is 1.031×108 s-1.
3.2 Mechanism function of diaspore dissolution
Figure 3 displays the DSC curve of diaspore sample dissolving in sodium hydroxide solution with scanning temperature from 20 °C to 300 °C. On the start of the curve, there is a fluctuation before 70 °C, which is led by the error of equipment and the unstable voltage. There is a small peak appearing around 235 °C resulting from the impurities in the diaspore sample gained by flotation. Around 110 °C, a main endothermic peak turns up, and climbs up with the increase of temperature. The endothermic peak reaches the highest at 121 °C, and ends around 195 °C. The initial value of parameter in
Fig.3 DSC curve of diaspore and caustic solution scanning from 20 °C to 300 °C with heating rate of 10 °C/min
Eq.(2), T0, is set as 384.5 K. The range of T in accordance with the start and end of peaks is considered to be from 385.5 K to 465.5 K.
Iterative method and linear least square method were carried out using the data of 30 points in the temperature range of 112 °C to 192 °C on the DSC curve obtained by scanning mixture of diaspore sample and caustic solution. Following the same procedure described in section 3.1, the calculated results run by Matlab program are shown in Table 3. It is indicated that eight functions show the non-existence results because the iteration is more than 500 times. And the values of apparent activation energy of three functions are negative and the preexponential constant values of 13 functions are beyond the value range of 106 to 1013. Among the left six functions, Nos.28-30 show the correlation coefficient over 0.99, and the low residual variance. However, GU et al[7-8] revealed the apparent activation energy of diaspore ore digesting in sodium aluminate solution was around 44 kJ/mol and 85 kJ/mol for two reaction phases, respectively. Therefore, the last function (30) is the most suitable differential mechanism function to describe the process of diaspore dissolving in caustic solution. The kinetic parameter, apparent activation energy, with respect to the No.30 function, f(α)=2(1-α)3/2, is 73.858 kJ/mol, and the preexponential constant A is 5.752×107 s-1.
Table 3 Kinetic parameters obtained by iterative method of 30 differential mechanism functions using data of DSC curve of diaspore sample dissolving in caustic solution
4 Conclusions
1) A single non-isothermal DSC curve is the simple method to investigate the kinetics of AlOOH dissolving in sodium hydroxide solution, by carrying out high- pressure DSC technology. Though boehmite and diaspore have the same chemical formula, AlOOH, their differential mechanism functions, which are logically selected from 30 types of functions by carrying out iterative method and linear least square method, are different.
2) The most probable differential mechanism function of boehmite dissolving in caustic solution is f(α)=1-α, which shows the first-order reaction with apparent activation energy of 79.178 kJ/mol and the preexponential constant of 1.031×108 s-1.
3) The mechenism function, f(α)=2(1-α)3/2, can probably describe the dissolution of diaspore sample in alkali solution. The calculated results of kinetic parameters are apparent activation energy of 73.858 kJ/mol, preexponential constant of 5.752×107 s-1 and reaction order of 1.5.
References
[1] HABASHI F. Kinetics of metallurgical processes [M]. Canada: Metallurgie Extractive Quebec, 1999.
[2] REDDY B R, MISHRA S K, BANERJEE G N. Kinetics of leaching of a gibbsitic bauxite with hydrochloric acid [J]. Hydrometallurgy, 1999, 51(1): 131-138.
[3] ZIVKOVIC Z D, STRBAC N D. Kinetics and mechanism of leaching of low-quality boehmite bauxite with hydrochloric-acid [J]. Hydrometallurgy, 1994, 36(2): 247-258.
[4] PEREIRA J A M, SCHWAAB M, DELLORO E, PINTO J C, MONTEIRO J L F, HENRIQUES C A. The kinetics of gibbsite dissolution in NaOH [J]. Hydrometallurgy, 2009, 96(1/2): 6-13.
[5] ZAFAR Z I. Determination of semi empirical kinetic model for dissolution of bauxite ore with sulfuric acid: parametric cumulative effect on the Arrhenius parameters [J]. Chemical Engineering Journal, 2008, 141(1/3): 233-241.
[6] BRITTAN M I. Variable activation energy model for leaching kinetics [J]. International Journal of Mineral Processing, 1975, 2: 321-331.
[7] GU Song-qing, CAO Rong-jiang, CHEN Xin-min. Steady state mathematical model of diasporic bauxite digestion [J]. Nonferrous Metals, 1986, 38(2): 66-78. (in Chinese)
[8] GU Song-qing, CAO Rong-jiang, CHEN Xin-min. Study on digestion kinetics of diasporic bauxite [J]. Acta Metallurgica Sinica, 1987, 23(6): B269-B276. (in Chinese)
[9] BI Shi-wen, LI Dian-feng, YANG Yi-hong, MASANELI D. Kinetic model for diaspore digestion [J]. Journal of Northeastern University (Natural Science), 1995, 19(3): 302-306. (in Chinese)
[10] BORNTRAG A, KNIRSCH D, FETTING F. Rate of dissolution of hydrargillite-containing and boehmite-containing bauxites in caustic soda [J]. Chemie Ingenieur Technik, 1974, 46(15): 659-665.
[11] LI Chao-qun, CHEN Qi-yuan, CHEN Xin-min. Investigation of digestion process kinetics of gibbsite [J]. Nonferrous Metals, 1991, 43(4): 52-55. (in Chinese)
[12] YIN Ai-jun, CHEN Qi-yuan, ZHANG Ping-min. Studies on the kinetics of the digestion process of synthetic gibbsite by DSC [J]. Chemical Journal of Chinese University, 1991, 12(11): 1507-1509. (in Chinese)
[13] ZHOU Qiu-sheng, LI Xiao-bin, XIONG Xiang, LIU Ye-xiang, ZHAI Yu-chun. Hydrothermal synthesis of nano-meter sized diaspore powder [J]. Mining and Metallurgical Engineering, 2001, 21(1): 51-53. (in Chinese)
[14] LIU Shui-hong, FANG Qi-xue. Commentary on present situation of research into bauxite desilication technology [J]. Mining and Metallurgy, 2004, 13(4): 24-29. (in Chinese)
[15] ZHAO Shi-min, WANG Dian-zuo, HU Yue-hua, XU Jing, ZHAO Xiao-ling. Current study status of the pre-desilicification for bauxite [J]. Mining R&D, 2004, 24(5): 37-44. (in Chinese)
[16] CHEN Xiang-qing, BAI Wang-quan, YAN Wei-zhe. Status and research progress of bauxite floatation desilication [J]. Light Metals, 2006(2): 8-12. (in Chinese)
[17] HU Rong-zu, GAO Sheng-li, ZHAO Feng. Thermal analysis kinetics [M]. Beijing: Science Press, 2008. (in Chinese)
[18] LIU Wei-guo. Design and application of MATLAB program [M]. Beijing: Higher Education Press, 2006. (in Chinese)
[19] MOORE J W, PEARSON R G. Kinetics and mechanism [M]. New York: John Wiley & Sons, Inc., 1981.
利用高压DSC研究一水铝石在碱液中的溶出动力学
鲍 丽，张廷安，豆志河，吕国志，郭永楠，倪培远，吴许建，马佳
东北大学 材料与冶金学院, 沈阳 110004
摘 要：研究一水铝石(包括一水软铝石和一水硬铝石)在氢氧化钠溶液中溶出过程的机理函数并确定相应的动力学参数。将一水铝石矿和氢氧化钠溶液在高压DSC设备中以10 °C/min的速度加热，用微分方程法分析得到DSC曲线，同时借助Matlab程序利用迭代法和线性最小二乘法计算微分方程，从30种机理函数的微分形式中分别逻辑地选择了两种一水铝石矿和苛性碱溶液的反应过程的最可能的机理函数。计算结果表明：一水软铝石溶解于苛性碱溶液的最可能的机理函数为f(α)=1-α，表观活化能为79.178 kJ/mol，指前因子为1.031×108 mol/s，反应级数为1；一水硬铝石在苛性碱溶液中溶解的最可能的机理函数为f(α)=2(1-α)3/2，表观活化能为73.858 kJ/mol，指前因子为5.752×107 mol/s，反应级数为1.5。
关键词：微分机理函数；一水软铝石；一水硬铝石；表观活化能；指前因子
(Edited by LAI Hai-hui)
Foundation item: Project (2007BC13504) supported by the National Basic Research Program of China; Project (20050145029) supported by Research Fund for the Doctoral Program of Higher Education; Project (2005221012) supported by the Science and Technology Talents Fund for Excellent Youth of Liaoning Province, China
Corresponding author: ZHANG Ting-an; Tel: +86-24-83687715; E-mail: zta2000@163.net
DOI: 10.1016/S1003-6326(11)60695-X
