基于Popescu法的高铁三水铝土矿非等温热分解动力学
来源期刊:中国有色金属学报(英文版)2015年第7期
论文作者:柳政根 王 峥 唐 珏 王宏涛 龙红明
文章页码:2415 - 2421
关键词:高铁三水铝土矿;热分解动力学;活化能;指前因子;Popescu法
Key words:high iron gibbsite ore; thermal decomposition kinetics; activation energy; pre-exponential factor; Popescu method
摘 要:在非等温条件下研究高铁三水铝土矿的热分解动力学。利用Pcpescu法分析高铁三水铝土矿的热分解反应机理。结果表明,Jander方程的三维扩散模型是高铁三水铝土矿热分解反应的最概然机理函数,函数机理代码为D3。由Popescu法得到的高铁三水铝土矿热分解活化能和指前因子分别为75.36 kJ/mol和1.51×105 s-1。由等转化率方法求出的活化能验证了所求机理函数的正确性。Popescu法是分析高铁三水铝土矿热分解反应机理较为合理和可靠的方法。
Abstract: The thermal decomposition kinetics of high iron gibbsite ore was investigated under non-isothermal conditions. Popescu method was applied to analyzing the thermal decomposition mechanism. The results show that the most probable thermal decomposition mechanism is the three-dimensional diffusion model of Jander equation, and the mechanism code is D3. The activation energy and pre-exponential factor for thermal decomposition of high iron gibbsite ore calculated by the Popescu method are 75.36 kJ/mol and 1.51×105 s-1, respectively. The correctness of the obtained mechanism function is validated by the activation energy acquired by the iso-conversional method. Popescu method is a rational and reliable method for the analysis of the thermal decomposition mechanism of high iron gibbsite ore.
Trans. Nonferrous Met. Soc. China 25(2015) 2415-2421
Zheng-gen LIU1, Zheng WANG1, Jue TANG1, Hong-tao WANG1, Hong-ming LONG2
1. School of Materials and Metallurgy, Northeastern University, Shenyang 110819, China;
2. School of Metallurgical Engineering, Anhui University of Technology, Ma’anshan 243002, China
Received 16 September 2014; accepted 15 April 2015
Abstract: The thermal decomposition kinetics of high iron gibbsite ore was investigated under non-isothermal conditions. Popescu method was applied to analyzing the thermal decomposition mechanism. The results show that the most probable thermal decomposition mechanism is the three-dimensional diffusion model of Jander equation, and the mechanism code is D3. The activation energy and pre-exponential factor for thermal decomposition of high iron gibbsite ore calculated by the Popescu method are 75.36 kJ/mol and 1.51×105 s-1, respectively. The correctness of the obtained mechanism function is validated by the activation energy acquired by the iso-conversional method. Popescu method is a rational and reliable method for the analysis of the thermal decomposition mechanism of high iron gibbsite ore.
Key words: high iron gibbsite ore; thermal decomposition kinetics; activation energy; pre-exponential factor; Popescu method
1 Introduction
High iron gibbsite ore consists of goethite (FeOOH), gibbsite (Al(OH)3), diaspore (AlOOH), and little other minerals. The high iron gibbsite ore spreads over some countries, such as China, Australia, Greece, Brazil and America. The national reserve of high iron gibbsite ore in China attains to more than 1.5 billion tons [1]. With the aim to achieve comprehensive utilization of this ore, many researches have been carried out, and several treatment processes have been presented successively, such as direct reduction-magnetic separation of the iron-alumina leaching [2-5], sintering agglomeration-blast furnace ironmaking-alumina leaching [6,7] and direct reduction-smelting separation of the iron-alumina leaching [8]. These processes have a common characteristic of ironmaking first and then alumina leaching. If the problems of high energy consumption and high cost in these processes can be solved, these processes would be able to develop well.
During the ironmaking step in these processes, the thermal decomposition of high iron gibbsite ore would be generated. The thermal decomposition kinetics of high iron gibbsite ore contributes to the development and improvement of these processes. However, the thermal decomposition mechanism of high iron gibbsite ore has not been reported. Most studies of thermal decomposition kinetics of bauxite focused on the gibbsite. ZHU et al [9] researched dehydration reaction kinetic of gibbsite with Kissinger equation. MACKENZIE et al [10] summarized the thermal decomposition of mechanically activated gibbsite. AMIRPIRAN et al [11] studied the non-isothermal dynamic model for the thermal decomposition of a gibbsite particle. And also AMIRPIRAN et al [12] developed a multi-stage, multi-reaction shrinking core model for self-inhibiting gas-solid reactions of gibbsite. However, the researches on the non-isothermal thermal decomposition kinetics of high iron gibbsite ore are quite limited.
This work focused on the non-isothermal thermal decomposition kinetics of high iron gibbsite ore. And its corresponding kinetics parameters were researched by Popescu method. This work could provide a deep insight into the thermal decomposition mechanism of high iron gibbsite ore. And also it could ensure optimal and proper control for the thermal decomposition conditions of high iron gibbsite ore, which are helpful to the efficient energy saving and prospective utilization of high iron gibbsite ore in various processes.
2 Experimental
2.1 Samples preparation
The high iron gibbsite ore samples used in this investigation were obtained from Guangxi region in China. Before experiments, the samples were firstly ground to particle size less than75 μm, and then they are dried at 105 °C for 5 h.
Table 1 shows the major chemical composition of high iron gibbsite ore. It can be seen that the alumina and ferric oxide contents are 23.85% and 49.21%, respectively. The high iron gibbsite ore is a run-of-mine ore with characteristics of high iron, high silica and low A/S ratio.
Table 1 Major chemical compositions of high iron gibbsite ore (mass fraction, %)
The mineral composition of high iron gibbsite ore was investigated by means of X-ray diffraction. Figure 1 shows the XRD pattern of high iron gibbsite ore analyzed by the search match software. It can be concluded that the main crystalline phases are goethite (FeOOH), gibbsite (Al(OH)3), diaspore (AlOOH), hematite (Fe2O3), and silica (SiO2).
Fig. 1 XRD pattern of high iron gibbsite ore
2.2 Thermogravimetric experiments
Non-isothermal decomposition experiments were carried out using a Netzsch STA 409 C/CD in Ar atmosphere (99.99%, volume fraction). A sample of about 10 mg was taken for each run and heated from room temperature to 800 °C with a steady Ar flow (40 mL/min) maintained through the reactor tube. The samples were used at different heating rates, 5, 10, 15, 20 °C/min, respectively.
2.3 Popescu analysis method
Popescu method, as well as a variant on the Flynn-Wall-Ozawa method, was put forward by POPESCU [13] in 1996. The Popescu method is a method of multiple scanning rate. The activation energy E, the pre-exponential factor A and the mechanism functions f(α) were obtained by dealing with the degree of conversion at the same temperature on several TG curves at different heating rates. The main advantages of this method are that it neither includes any assumption concerning the temperature integral, nor takes the form of k(T) into account, and the results have a higher precision [14]. And the method allows the activation energy to be determined as a function of the extent of conversion and/or temperature without making any assumption about the reaction model.
In general, the governing equation for kinetics analysis of solid state decomposition is as follows:
(1)
where t is the time, T is the temperature, α is the extent of conversion and k(T) is the temperature dependence of the rate constant.
By the usual change of the variable time into temperature, Eq. (1) becomes
(2)
where β=dT/dt is the heating rate.
Because in most experiments, the heating rate is kept constant, in Popescu method, the integral form of Eq. (2) is
(3)
where αm and αn are two different degrees of conversion and Tm and Tn are their corresponding temperatures. By using the notations
(4)
(5)
(6)
Equation (6) is actually the basis of the method for finding the most probable kinetic mechanism of the studied reaction.
It has to be assumed that over certain ranges of α and β values, the kinetics of the reaction does not change. For example, the analytical forms of f(α) and k(T), consequently G(α)mn and I(T)mn, are not changed when either α or β is varied.
On the basis of the experimental data, a pairs of α, such as (αm1, αn1), (αm2, αn2), …, (αmi, αni) will be determined, when T=Tm and Tn, under the different heating rates β. With the help of these pairs and various functions, the values of G(α)mn1, G(α)mn2, …, G(α)mni could be computed based on Eq. (4). And for each conversation function f(α), the corresponding values were acquired. As the temperatures Tm and Tn are the same for all the experiments, according to Eq. (5), it follows that I(T)mn is constant, and from Eq. (6), a plot of the values of G(α)mn versus 1/β will lead to a straight line with an intercept of zero if the analytical form of f(α) is properly selected. With the data analysis, a series of straight lines will indicate the best kinetics model by the best linear correlation coefficient.
The equations for calculating activation energy E and pre-exponential factor A are as follows:
(7)
where . By substituting the values of (Tm1, Tn1), (Tm2, Tn2), …, (Tm5, Tn5) and mechanism functions G(α)mn into Eq. (7), the kinetic parameters of E, A and linear correlation r were obtained by the linear least squares method with ln[β/(Tn-Tm)] versus 1/Tξ.
3 Results and discussion
3.1 Analysis on TG-DSC curves
Figures 2-4 give the TG curves, DSC curves and conversion degree curves of high iron gibbsite ore heated at different heating rates, respectively. From Fig. 2, the TG curves of high iron gibbsite ore could be divided into three stages roughly. The segment points of each stages have some changes and hysteretic phenomenon with the increase of the heating rate. From Fig. 3, there is an obvious endothermic peak of each curve with the temperature between 250 and 290 °C. However, the endothermic peaks with the temperature between 50 and 150 °C are not remarkable. As shown in Fig. 1, the high iron gibbsite ore mainly consists of goethite (FeOOH), gibbsite(Al(OH)3) and diaspore(AlOOH). These hydrates would decompose when they are heated to a certain temperature. Combining the characteristic of goethite, gibbsite and diaspore with the TG-DSC curves, it could be concluded that the thermal decomposition could divide into three stages [15,16]. The first is the stage of losing absorption water (0-200 °C), and the next is dehydroxylation of gibbsite and goethite (200-300 °C), and the last is dehydroxylation of diaspore (300-500 °C). Figure 5 shows the XRD pattern of high iron gibbsite ore heated at the temperatures of 600 and 800 °C, respectively. It can be seen that there is almost no change of the mineral compositions of high iron gibbsite ore heated at 600 and 800 °C. Based on Figs. 2 and 5, it could be drawn that the decomposition of minerals in the samples is completed when the temperature reaches 800 °C. As the characteristic of thermal decomposition of high iron gibbsite ore is the dehydration of goethite (FeOOH), gibbsite (Al(OH)3) and diaspore (AlOOH), it is possible that the thermal decomposition kinetics of high iron gibbsite ore is calculated by the mass loss or heat flow.
Fig. 2 TG curves of high iron gibbsite ore heated at different heating rates
Fig. 3 DSC curves of high iron gibbsite ore heated at different heating rates
Fig. 4 Conversion degree curves of high iron gibbsite ore heated at different heating rates
Fig. 5 XRD patterns of high iron gibbsite ore heated at 600 and 800 °C
3.2 Determination on thermal decomposition mechanism function
Forty-one kinds mechanism functions in reference were analyzed [17]. Among them, 6 kinds of diffusion model, 2 kinds of random nucleation-subsequent growth model and 2 kinds of reaction order model present a good linearity. The integral form G(α) and the differential form f(α) of these models are shown in Table 2 [18-21].
In order to determinate the thermal decomposition mechanism function of high iron gibbsite ore, as shown in Table 3, some temperatures Tm and Tn and the corresponding conversion degrees for each of them are selected. The selected temperatures are around 300 °C (573.15 K), because these temperatures could represent the relationship among conversion degrees, temperatures and heating rates at the greatest extent, so as to reduce the random error resulted from selecting thermal decomposition temperature.
By substituting the conversion degrees for each of them in Table 3 and kinetics mechanism functions in Table 2 into Eq. (4), respectively, the values of G(α)mn at different heating rates for every conversion function are obtained. Plotting G(α)mn versus 1/β and applying a linear fitting least squares method with an intercept of zero, the inferring schematic diagram of the most probability mechanism function is shown in Fig. 6.
The calculated linear correlation coefficients r and standard deviations rSD for each kinetic mechanism functions are listed in Table 4. Liner correlation coefficients r reflect the goodness of fitting for linear logarithmic stronghold. The standard deviations rSD stand for the variation size between actual value and the return value of linear, which are expected as smaller as possible. Thus, taking linear correlation coefficients r and standard deviations rSD into the account, the mechanism functions would be ensured.
From Table 4, it can be seen that the data of No.3 are better, the correlation coefficients r are greater and the standard deviations rSD are relatively smaller. Therefore, three-dimensional diffusion model of Jander equation is the most probable mechanism function for the thermal decomposition of high iron gibbsite ore, and the integral form and differential form are G(α)= [1-(1-α)1/3]2 and f(α)=(3/2)(1-α)2/3[1-(1-α)1/3]-1, respectively. And the mechanism code is D3, which is different from the mechanism function of gibbsite [9,10].
3.3 Calculation of kinetics parameters
In Propescu method, according to Eq. (7), the plot of ln[β/(Tn-Tm)] versus 1/Tξ leads to a straight line, whose slope is the value of the activation energy E. Values of the activation energy are thus computed, as well as the corresponding correlation coefficients. The values of the pre-exponential factor A are calculated from the values of the intercept of the plot. These results are listed in Tables 5 and 6. From Table 6, it can be seen that the largest correlation coefficient r is 0.95855. Thus, the most probable activation energy E and pre- exponential factor A are 75.36 kJ/mol and 1.51×105 s-1, respectively.
Table 2 Typical kinetics mechanism functions
Table 3 Conversion degrees measured for given temperatures with different heating rates
Fig. 6 Inferring schematic diagram of most probability mechanism function
3.4 Activation energy verification
Only the correct mechanism function corresponds to the right activation energy. Thus, the right activation energy also can test whether the mechanism function is correct [18,19].
In order to verify the thermal decomposition activation energy, Flynn-Wall-Ozawa (FWO) method was used to calculate the activation energy. FWO method was a typical method of iso-conversional method [22]. Compared with other methods, the method avoids the possibility of error, which resulted from the different kinetic mechanism function hypotheses, so the activation energy obtained by FWO method is more reasonable [23].
The equation of FWO method is as follows:
(8)
where TP is the peak temperatures of corresponding DSC curves. In this work, the TP is obtained from Fig. 3. Table 7 gives the activation energy E and correlation coefficients r of the endothermic peaks at different conversion degrees. As shown in Table 7, the average activation energy E is the average value of the activation energy of the endothermic peaks at different conversion degrees, and the value is of 79.44 kJ/mol, which relatively matches the value (75.36 kJ/mol) obtained by Popescu method. And the values satisfy the condition of |(Ep-EF)/Ep|=0.05≤0.1 [17]. Thus, the correctness of the thermal decomposition mechanism function is proved from the view of the activation energy.
Table 4 Linear fitting results of kinetic mechanism functions
Table 5 Temperatures measured for various conversion degrees on curves from Fig. 4
Table 6 Values of kinetic parameters for thermal decomposition of high iron gibbsite ore computed by Eq. (7)
Table 7 Activation energy and correlation coefficients of endothermic peaks at different conversion degrees
4 Conclusions
1) The thermal decomposition mechanism of high iron gibbsite ore follows the three dimensional diffusion model of Jander equation, whose integral form and differential form are G(α)=[1-(1-α)1/3]2 and f(α)=(3/2)(1-α)2/3[1-(1-α)1/3]-1, respectively.
2) The activation energy and pre-exponential factor for thermal decomposition kinetics of high iron gibbsite ore are 75.36 kJ/mol and 1.51×105 s-1, respectively.
3) The correctness of the thermal decomposition mechanism function is proved from the view of the activation energy by the FWO method, and the Popescu method is a rational and reliable method for the research on the thermal decomposition mechanism of high iron gibbsite ore.
References
[1] MU Xin-he. Discussion on reasonable exploitation of bauxite mineral resources in China [J]. Mineral Resources and Geology, 2002, 16(5): 313-315. (in Chinese)
[2] LI Guang-hui, LIU Mu-dan, JIANG Tao, ZHOU Tai-hua, FAN Xiao-hui. Mineralogy characteristics and separation of aluminum and iron of high-aluminum iron ores [J]. Journal of Central South University: Science and Technology, 2009, 40(5): 1165-1171. (in Chinese)
[3] HU Wen-tao, WANG Hua-jun, LIU Xin-wei, SUN Chuan-yao. Effect of nonmetallic additives on iron grain grindability [J]. International Journal of Mineral Processing, 2014, 130(7): 108-113.
[4] LIU Zheng-gen, CHU Man-sheng, TANG Jue, HAN Yuan-ting, WU Xiang-long. Appropriate reduction and Fe-Al separation of high iron gibbsite [C]//SADLER B A. Light Metals 2013. San Antonio, TX: TMS, 2013: 223-227.
[5] PICKLES C A, CHAMBERS T L B, FORSTER J. A study of reduction and magnetic separation of iron from high iron bauxite ore [J]. Canadian Metallurgical Quarterly, 2012, 51(4): 424-433.
[6] SHI Guo-song, FANG Jue, YANG Gai-yan, GAO Yan-jia. Feasibility study of Guangxi high Fe bauxite sintering [J]. Journal of Hebei Polytechnic University: Natural Science Edition, 2011, 33(2): 11-13. (in Chinese)
[7] LI Yin-tai, BI Shi-wen, DUAN Zheng-yin, YANG Yi-hong, ZHANG Jing-dong. Discussion on the comprehensive utilization technology of high iron gibbsite ore in Guigang city Guangxi region [J]. Light Metals, 1992(9): 6-14. (in Chinese)
[8] ZHANG Zuo-liang, LI Qiang, ZOU Zong-shu. Reduction properties of high alumina iron ore cold bonded pellet with CO-H2 mixture [J]. Ironmaking and Steelmaking, 2014, 41(8): 561-567.
[9] ZHU Bo-quan, FANG Bin-xiang, LI Xiang-cheng. Dehydration reaction and kinetic parameters of gibbsite [J]. Ceramics International, 2010, 36(8): 2493-2498.
[10] MACKENZIE K J D, TEMUUJIN J, OKADA K. Thermal decomposition of mechanically activated gibbsite [J]. Thermochimica Acta, 1999, 327(1-2): 103-108.
[11] AMIRPIRAN A, ANDREY V B, GORDON D I, IZTOK L, NICOLETA E M. A 1-D non-isothermal dynamic model for the thermal decomposition of a gibbsite particle [J]. Chemical Engineering Research and Design, 2013, 91(3): 485-496.
[12] AMIRPIRAN A, GORDON D I, ANDREY V B, IZTOK L, NICOLETA E M. A multi-stage, multi-reaction shrinking core model for self-inhibiting gas–solid reactions [J]. Advanced Powder Technology, 2013, 24(4): 728-736.
[13] POPESCU C. Integral method to analyze the kinetics of heterogeneous reactions under non-isothermal conditions a variant on the Ozawa-Flynn-Wall method [J]. Thermochimica Acta, 1996, 285(2): 309-323.
[14] ZHANG Jian-jun, REN Ning, BAI Ji-hai. Non-isothermal decomposition reaction kinetics of the magnesium oxalate dihydrate [J]. Chinese Journal of Chemistry, 2006, 24(3): 360-364.
[15] GAO Zhen-xin, HE Zhong-yang, ZHENG Xiao-ping, FU Qing-hua, SHI Xu-bo. Phase transformation and morphology of Bayer-gibbsite during heating [J]. Journal of the Chinese Ceramic Society, 2008, 36(S1): 117-123. (in Chinese)
[16] LIU Xue-fei, WANG Qing-fei, ZHANG Qi-zuan, ZHOU Fang, GAO Bang-fei, XU Hao. Thermal analysis of core of the ore body X in the Xinxu bauxite deposit, Jinxi county, Guangxi region [J]. Journal of Mineralogy and Petrology. 2008, 28(4): 53-59. (in Chinese)
[17] HU Rong-zu. Thermal analysis kinetics [M]. 2nd ed. Beijing: Science Press, 2008. (in Chinese)
[18] DICKINSON C F, HEAL G R. A review of the ICATC kinetics project, 2000 part 1: Isothermal results [J]. Thermochimica Acta, 2009, 494(1-2): 1-14.
[19] ZHANG Kun, LIN Shao-kun, LIN Mu-liang. Introduction to a new method for thermal analysis kinetics multivariate non-linear regression and its application [J]. Modern Scientific Instruments, 2002(5): 15-18. (in Chinese)
[20] CRADO J M, MDEK J, ORTEG A. Applicability of the master plots in kinetic analysis of non-isothermal data [J]. Thermochimica Acta, 1989, 147(2): 377-385.
[21] LIU Zhong-suo, WANG Qi, ZOU Zong-shu, TAN Guang-lei. Popescu method used for study on reaction mechanism of carbon gasification with CO2 [J]. Research on Iron & Steel. 2011, 39(6): 12-17. (in Chinese)
[22] OZAWA T. A new method of analyzing thermogravimetric data [J]. Bulletin of the Chemical Society of Japan, 1965, 38(11): 1881-1886.
[23] ZHANG Xiang-hui, HE Chuan, WANG Ling, LIU Jing, DENG Miao, FENG Qian. Non-isothermal kinetic analysis of thermal dehydration of La2(CO3)3·3.4H2O in air [J]. Transactions of Nonferrous Metals Society of China, 2014, 24(10): 3378-3385.
柳政根1,王 峥1,唐 珏1,王宏涛1,龙红明2
1. 东北大学 材料与冶金学院,沈阳 110819;
2. 安徽工业大学 冶金工程学院,马鞍山 243002
摘 要:在非等温条件下研究高铁三水铝土矿的热分解动力学。利用Pcpescu法分析高铁三水铝土矿的热分解反应机理。结果表明,Jander方程的三维扩散模型是高铁三水铝土矿热分解反应的最概然机理函数,函数机理代码为D3。由Popescu法得到的高铁三水铝土矿热分解活化能和指前因子分别为75.36 kJ/mol和1.51×105 s-1。由等转化率方法求出的活化能验证了所求机理函数的正确性。Popescu法是分析高铁三水铝土矿热分解反应机理较为合理和可靠的方法。
关键词:高铁三水铝土矿;热分解动力学;活化能;指前因子;Popescu法
(Edited by Yun-bin HE)
Foundation item: Project (51374058) supported by the National Natural Science Foundation of China
Corresponding author: Hong-ming LONG; Tel: +86-555-2311571; E-mail: yaflhm@126.com
DOI: 10.1016/S1003-6326(15)63857-2