Advanced
Combustion Kinetics of Pulverized Indian Coal-Char in Different CO<sub>2</sub>-O<sub>2</sub> Mixture Isothermally
Combustion Kinetics of Pulverized Indian Coal-Char in Different CO2-O2 Mixture Isothermally
Journal of the Korean Chemical Society. 2009. Dec, 53(6): 635-639
Copyright © 2009, The Korean Chemical Society
  • Received : August 03, 2009
  • Accepted : October 26, 2009
  • Published : December 20, 2009
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Saravanan, V.
Central Power Research Institute, Bangalore, India
Shivakumar, R.
Niruguna babu, P.
, Ramakrishna
Central Power Research Institute, Bangalore, India

Abstract
Experimental work was carried out to study the combustion kinetics of the Indian coal-char in the varying mixture of carbon dioxide-oxygen (CO 2 -O 2 ). The coal sample was pulverized and sieved to less than 58 microns and charred using volatile furnace by passing the nitrogen gas. The experiments were carried out using the Thermo Gravimetric Analyzer (TGA-50) at CPRI, Bangalore, different proportions of (CO 2 -O 2 ) gas was allowed in to the TGA-50 (80-20, 60-40, 40-60, 20-80) mole basis were used to study the combustion kinetics of coal Isothermally, kinetic parameters like Activation energy (E) and the pre-exponential factors (A) are calculated using the unification approach and modified Arrhenius equation.
Keywords
INTRODUCTION
Many researchers in the past have studied the combustion kinetics of the coal-chars either by isothermal or non-isothermal analysis, 1 - 7 the reaction kinetics of coal-chars concerns correlation with char morphology, pore structure (porosity), surface area. Surface area and porosity of the prepared coal-char was determined using BET nitrogen gas adsorption. To estimate the combustion kinetics of coal-char, we need to consider the reactions occurring on surfaces (external and internal) of solid particles. There are different types of models farmed to study the kinetics of coal-chars, one such model is shrinking unreacted core model (SCM). 1 , 3 , 5 In this model chemical reaction occurring at the pore surface is the rate controlling step at low temperatures, here reaction starts on the outer surfaces and later ash layer is formed and then the reaction processed inside core of the particle leaving behind the ash or unreacted minerals, hence at any time during the reaction their exit a unreacted core which shrinks in size during the reaction hence the model named as shrinking unreacted core model. 1 , 3
The aim of this paper is to estimate the overall reaction rate, which is the product of conversion factor (f (X)) and intrinsic reaction rate (r s ). The activation energy (E) and pre-exponential factor (A) are determined by a unification approach using Random Pore Model (RPM), for experiments conducted isothermally using Thermogravimetric analyzer (TGA) for different CO 2 -O 2 proportion mixtures viz (80-20, 60-40, 40-60, 20-80), were the oxygen concentration was maintained constant for all the cases.
EXPERIMENTAL PROCEDURE
- Char preparation and Coal properties
Pulverized coal sample was heated in volatile furnace Leco (U.S.A), in flowing nitrogen gas environment at 950 ℃ at a rate of 20 liters/min and socked for 60 min, after 60 min the furnace was cooled to room temperature in flowing nitrogen gas environment and later the resulted coal-char was homogenized and stored for combustion studies in TGA. The char and coal samples, were analyzed for proximate, ultimate, particle size shown in 1 and BET nitrogen surface area, are shown in 2 .
Proximate and Ultimate analysis of coal and char samples
PPT Slide
Lager Image
Proximate and Ultimate analysis of coal and char samples
BET surface area analysis of char samples
PPT Slide
Lager Image
BET surface area analysis of char samples
Activation energy and lumped Pre exponential factor for different gas environment
PPT Slide
Lager Image
Activation energy and lumped Pre exponential factor for different gas environment
- Experimental work
Thermogravimetric analyzer (TGA) used here was TGA-50 model. Different proportions of (CO 2 -O 2 ) gas mixtures were fed to TGA at room temperatures and pressure more than atmospheric pressure (1.1 - 1.2 atm). The pre-mixed and calibrated (CO 2 -O 2 ) mixture cylinders were used for analysis. The oxygen concentration for all the case was maintained constant by mainly fixing the flow rate of the oxygen same for the cases i.e. (10.5 mL/min). Isothermal experiments were carried out with prepared coalchar (15 - 18 mg) in the presence of the (CO 2 -O 2 ) mixtures with constant heating rate 20 ℃/min for different temperatures 500, 525, 550, 575 and 600 ℃ respectively. Flow rate of the gas mixtures were varied by keeping oxygen flow rate constant at 10.5 mL/min.
ISOTHERMAL METHOD FOR EVALUATION OF MODEL AND ARRHENIUS PARAMETERS
A description of isothermal method used for our calculation, was earlier used by 1 , 9 , 10 ,11 with application to char combustion with mixture of (N 2 -CO 2 ) or pure CO 2 using Random pore model (RPM). The overall reaction rate or rate of the reaction incorporating the conversion factor f (X) is given by
PPT Slide
Lager Image
Were, rs is the intrinsic reaction rate (m/min). From RPM, conversion in terms of dimension less time is given by 9 , 10
PPT Slide
Lager Image
on differentiating equation (2), we get
PPT Slide
Lager Image
After integrating equation (3), the time constant
PPT Slide
Lager Image
Time factor,
PPT Slide
Lager Image
PPT Slide
Lager Image
But according to Arrhenius equation the reaction rate constant ‘K’ is given by
PPT Slide
Lager Image
Substituting equation (7) in to equation (6)
PPT Slide
Lager Image
Substituting equation (8) in to equation (5)
PPT Slide
Lager Image
In the above equation {(AS 0 )/(1-ε 0 )} = K 0 . Were K 0 is the lumped pre exponential factor.
PPT Slide
Lager Image
Taking logarithm of equation (9)
PPT Slide
Lager Image
Partial pressure of oxygen (P O2 ) is constant for a particular temperature. The order of the reaction can be obtained by plotting ln t f vs . ln P O2 . Further, plotting ln t f vs . 1/T, the E and K 0 will be obtained from the slope and intercept respectively.
Replace dτ with dt in equation (2)
PPT Slide
Lager Image
Where:
PPT Slide
Lager Image
for RPM model.
RESULTS AND DISCUSSION
- Property of coal char
Char prepared from volatile furnace was subjected to proximate (TGA-600, Leco), ultimate (CHNS-2000, Leco), particle size (Malvern particle size analyser) and BET nitrogen adsorption analysis (Micrometrics GEMINI 2360). The results are given in 1 .
The particle size of parental coal and its respective char was analysed using the Master sizer X (Malvern make), which works on the principle of laser ensemble light scattering technique.
A pore size distribution data were obtained from BET analysis using BJH approach this measurement was uniform and repeatable results were obtained, using this data we plotted a graph of pore radius against the pore volume is shown to get 2 nd order polynomial equation of pore volume in terms of pore radius which up on the integration will give the porosity ε o , External surface area per unit volume S o and pore length per unit volume L o using the following relation given by Bhatia. 9 , 10
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
- Conformation of model and determination of reactivity constants (isothermally)
The experiment was carried out with the char prepared in volatile furnace from the data obtained conversion on (ash dry free basis) as the ash is the minerals left out and will not take part in the combustion was calculated. Conversion X is calculated by
PPT Slide
Lager Image
The equation (11) was used for validating the shrinking core unreacted model (SCM) and the Random pore model (RPM), from BET results we can observe that the porosity values of the char sample was very high around 45%. Hence the SCM is not suitable model to calculate kinetics of the char because the SCM model does not take in to the account of the structural properties of the char, the reactivity of the char is confined to availability of the surface area. i.e. initially the pore surface area will grow therefore the rate of reaction will increase but latter due to the overlapping of growing surfaces will decrease the rate of reactions, this concept is very well explained in the RPM, therefore the RPM is used to calculate the kinetic parameters. The complete conversion was not able to be at temperatures 500 and 525 ℃, hence the Arrhenius plot was plotted for 50% conversion.We have got different activation energy (E), lumped pre-exponential factor (K 0 ) and pre exponential factor (A) for different (CO 2 -O 2 ) mixtures which are in range with earlier works.
PPT Slide
Lager Image
BJH adsorption pore distribution curve for char-A.
CONCLUSION
The combustion profiles of the char in CO 2 -O 2 mixtures for temperatures 600, 575 and 550 ℃ showed the influence of CO 2 concentration over the reactivity. The overall reactivity increases when the concentration of CO 2 decreases in the CO 2 -O 2 mixture for the constant O 2 concentration which indicates that the CO 2 have some cooling effect on the reactivity of char, because of which the differential weight loss with respect to time is slowed as concentration of CO 2 is increased even after maintaining same amount of O 2 concentration as shown in weight loss patterns at different temperatures in . 2 , 3 and 4 .
Here we have used a different approach to determine the Arrhenius parameters like activation energy (E) and lumped Pre exponential factor (K 0 ) experimentally which are very important inputs for numerical analysis of burning profiles of coal-chars in different (CO 2 -O 2 ) gas-environments using Computational Fluid Dynamics (CFD). The E and K 0 values were determined using the modified Arrhenius equation (10), were the partial pressure of the O 2 is constant through out the experiment for particular temperature, therefore lnP O2 term in the equation (10) will be constant for a particular temperature, if we plot the graph of lnt f versus 1/T the negative slope of the line give the Activation energy (E) in terms of KJ/mole and intercept will give the lumped pre-exponential factor (K 0 ) in terms of (1/min) pa -n . The E and K 0 values decreases as the concentration of CO 2 decreases which means as the peaking temperature comes down the activation energy required for the reaction to initiate decreases. The values experimentally determined are in range with the previous experimental results of other researchers.
PPT Slide
Lager Image
Rate of weight loss with time for char-A at 600 ℃.
PPT Slide
Lager Image
Rate of weight loss with time for char-A at 575 ℃.
PPT Slide
Lager Image
Rate of weight loss with time for char-A at 550 ℃.
PPT Slide
Lager Image
Determination of activation energies and lumped pre-exponential factors for char-A for 50% conversion.
PPT Slide
Lager Image
References
Everson R. C. , Neomagus H. W. J. P. , Kaitano R. , Falcon. R. , Du Cann V. M. Fuel article JFUE in press 3921 -
Orfao J. J. M. , Martins V. 2002 Therm. chim. act 390 195 - 211    DOI : 10.1016/S0040-6031(02)00133-8
Everson R. C. , Neomagus H. W. J. P. , Njapha D. 2006 Fuel 85 418 - 422    DOI : 10.1016/j.fuel.2005.07.003
Madhusudanan P. M. , Krishnan K. , Ninan K. N. 1986 Therm. chim. act 97 189 - 201    DOI : 10.1016/0040-6031(86)87019-8
Everson R. C. , Neomagus H. W. J. P. , Njapha D. , Kasaini H. 2006 Fuel 85 1067 - 1075    DOI : 10.1016/j.fuel.2005.10.025
Madhusudanan P. M. , Krishnan K. , Ninan K. N. 1993 Therm. chim. act 221 13 - 21    DOI : 10.1016/0040-6031(93)80519-G
Vyazovkin S. 2000 Therm. chim. act 355 155 - 163    DOI : 10.1016/S0040-6031(00)00445-7
Bhatia S. K. , Perlmutter D. D. 1980 AIChE Journal 26 (3) 379 - 386    DOI : 10.1002/aic.690260308
Bhatia S. K. , Perlmutter D. D. 1981 AIChE Journal 27 (2) 247 - 254    DOI : 10.1002/aic.690270211
Chodankar C. R. , Feng B. , Klimenko A. 2007 Proceed. Austr. Comb. Sym 82 - 85