Modeling the Nonisothermal Devolatilization Kinetics of Typical South

Dec 26, 2013 - South African Coals. Burgert B. Hattingh,* Raymond C. Everson, Hein W. J. P. Neomagus, John R. Bunt, Daniel van Niekerk, and Ben P. Ash...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/EF

Modeling the Nonisothermal Devolatilization Kinetics of Typical South African Coals Burgert B. Hattingh,* Raymond C. Everson, Hein W. J. P. Neomagus, John R. Bunt, Daniel van Niekerk, and Ben P. Ashton Research Focus Area for Chemical Resource Beneficiation, North-West University, Potchefstroom Campus, Private Bag X6001, Potchefstroom, 2520 North West, South Africa Energy Systems, School of Chemical and Minerals Engineering, North-West University, Potchefstroom Campus, Private Bag X6001, Potchefstroom, 2520 North West, South Africa Sasol Technology (Pty) Ltd., Research & Development, Coal & Gas Processing Technology, Box 1, Sasolburg, 1947 Free State, South Africa S Supporting Information *

ABSTRACT: Multicomponent model fitting was conducted in order to evaluate the devolatilization rate behavior of four typical South African coals, with the aid of nonisothermal thermogravimetry. Rate evaluation was conducted at four different heating rates (5, 10, 25, and 40 K/min) by heating the samples under an inert N2 atmosphere to 950 °C. Evaluation of the kinetic parameters of each coal involved the numerical regression of nonisothermal rate data in MATLAB 7.1.1 according to a pseudocomponent modeling philosophy. The number of pseudocomponents used ranged between three and eight, as larger values induced the risk of over fitting. Quality of fit (QOF) was found to decrease with decreasing heating rate as a result of improved separation of the individual component reactions at the lower heating rates. All four coals showed the occurrence of similar pseudocomponent reactions, although significant differences were observed in the fractional contributions of the different pseudocomponents to the overall reaction rates. Modeling results indicated that the assumption of eight pseudocomponents produced the lowest QOF values and subsequently the best fit to the devolatilization profiles of each coal. For the vitrnite-rich coals (G#5 and TSH), no remarkable decrease in QOF could be observed after 6 pseudocomponent reactions, suggesting that even 6 or 7 pseudocomponent reactions would have provided accurate experimental predictions. Activation energies determined from the selected number of pseudocomponents (between 3 and 8) were found to range between 20 and 250 kJ/mol.

1. INTRODUCTION Coal devolatilization plays an important role in not only the metallurgical industry for producing coke but also during coal gasification where it normally constitutes the initial step of the process.1,2 It is therefore important to evaluate the devolatilization behavior of a coal feedstock in order to assess and optimize the production of valuable products such as char/coke, tar, and gas. Extensive research during the past few decades has advanced our knowledge of the kinetics and mechanisms of the devolatilization process. Furthermore, it has also provided us with valuable techniques for predicting, to a reasonable extent, the behavior of coals.3−5 Devolatilization modeling is quite straightforward if the chemical reaction step is rate controlling and the fuels are of a simple characteristic nature. Kinetic models for describing thermal decomposition range in different levels of complexity from free radical mechanistic models for simple hydrocarbon species such as propane6 to more complex reaction schemes. The latter involves a number of individual reactions, incorporating extra transport steps such as in the case of naphtha devolatilization.7 The kinetic description of more complex polyaromatic substances such as coal presents a challenging task, due to a vast amount of reactions involved. Kinetic evaluation of these substances is therefore normally conducted using pseudomechanistic models, which attribute the overall measured reaction rate to the cumulative effect of a number of separate reactions.8,9 © 2013 American Chemical Society

A large number of possible modeling strategies are currently available, of which the simplest are empirical in nature and employ global kinetics.8−11 Available models can be divided into either general weight-loss models or structural models. Typical weight-loss models include models employing a (1) single rate, (2) two rates, (3) multiple rates, and (4) distributed rates.12−21 Although simplistic in nature, the validity of a single reaction rate mechanism is quite limited. Kinetic parameters derived at a single heating rate has been generally shown not to be appropriate to other heating rates.22 Currently, the Distributed Activation Energy Model (DAEM) (Anthony-Howard model) has been shown to be the most powerful model for predicting devolatilization behavior.10,18,19,22−25 This complicated model was first proposed by Pitt26 and assumes the devolatilization process to consist of an infinite series of independent parallel first-order reactions. Accordingly, coal devolatilization can be explained by a distribution of activation energies about some mean activation energy value (Ea,0). The function within the model f(Ea) accounts for a distribution of activation energies and is assumed to be of Gaussian form. A simplistic approach for solving the DAEM considers a common, Received: October 25, 2013 Revised: December 25, 2013 Published: December 26, 2013 920

dx.doi.org/10.1021/ef402124f | Energy Fuels 2014, 28, 920−933

Energy & Fuels

Article

first-order reaction model impractical in determination of the kinetic parameters. In a response to this, authors such as Alonso et al.28 formulated a lumped first-order model allowing for the fractional contribution (ξ) made by the different individual peaks or parallel first-order reactions. Extensive research has been conducted in an attempt to understand the devolatilization kinetics of northern hemisphere (Carboniferous) coals, which are normally characterized by their high proportions of vitrinite, large amounts of reactive inertinite, and very low ash contents.40−45 These coals have been shown to be very different from South African coals, due to their differences in depositional environment.46−48 In contrast, Permian aged Gondwanaland coals are mainly rich in inertinite with a high abundance of minerals, while some coal deposits (Grootegeluk and Venda-Pafuri districts) also display higher abundances of vitrinite-rich coals. Although very limited, investigations such as those conducted by Aboyade et al.49 have provided a means of describing the devolatilization behavior of some South African coals. However, this work did not include the evaluation of any decomposition reactions below 200 °C, therefore neglecting the release of inherent moisture, crystal water, and amorphous phases from the coal particles. An investigation was therefore undertaken to gain a better understanding of the kinetic behavior of four typical South African coals (including for the release of moisture). Evaluation of the devolatilization behavior of the selected coals entailed the determination of the total mass loss behavior of the amount of volatile matter and not of individual species. Mass loss behavior during devolatilization was therefore described at the hand of the pseudocomponent modeling approach as proposed by Alonso et al.28 Kinetic parameters derived using this approach can be ultimately used in the evaluation of large particle devolatilization models (including for heat and mass transport effects) for effectively describing commercial coal conversion processes. This paper therefore aims to provide a comprehensive kinetic description of a variety of South African coal types, ranging from differences in maceral composition to differences in swelling behavior.

constant frequency factor applicable for all activation energies (or reactions) in the distribution.27 This simplification has, however, been criticized by Alonso et al.28 who demonstrated that the isokinetic effect, which involves the mutual relationship between pre-exponential constant and activation energy, cannot be neglected in the calculation procedure. Furthermore, the validity of using a constant pre-exponential factor becomes questionable when the function f(Ea) spreads over a wide range of Ea values.29 The use of continuous distribution curves (i.e., Gaussian), for the description of overall activation energy by a constant average activation energy and a standard deviation of energies (σ), has also raised concerns.29,30 From this perspective, the use of a Gaussian distribution does not necessarily guarantee the correct prediction of the characteristic f(Ea) curve of the substance under investigation. Although a very powerful model for estimating devolatilization behavior, the DAEM requires an extensive numerical discretization procedure for solving the required parameters. Recent advances in the understanding of the coal molecular structure has led to the development of some structural/network models of which the FG-DVC (Functional Group-Depolymerization Vaporization Cross-linking), the FLASHCHAIN, and the percolation lattice theory are the most commonly known.31−33 More comprehensive chemical models have also been formulated for describing char and tar formation.31 These models address the influence of coal molecular structure on volatile evolution rate by assuming that aromatic groups are linked by bridges and peripheral groups.34,35 However, evaluation of the kinetic parameters relevant to the different bridge and peripheral functional groups requires extensive knowledge of the coal molecular structure as determined using advanced characterization techniques such as NMR (Nuclear Magnetic Resonance spectroscopy), etc. The choice of a suitable kinetic model should therefore bear relevance to the process under investigation. The norm for kinetic studies is to conduct reactions under isothermal conditions, especially for fast reactions such as devolatilization at high temperatures. Time-resolved measurements of coal devolatilization are therefore very difficult and present uncertainty due to the fact that the devolatilization process normally completes within a few seconds before the isothermal state is reached.9 Currently, nonisothermal techniques have proven to be more useful than isothermal techniques for deriving the kinetic triplets (activation energy, preexponential factor, and reaction order).9,28 In addition, the process of deconvolution of differential thermogravimetric (DTG) curves into pseudocomponent curves has been found to be much easier for model- and less complex carbon-containing compounds (such as biomass and oil shales) than coals.36−38 This challenge necessitates the need for further elaboration on the behavior of typical DTG curves of coals, in order to formulate and evaluate an appropriate model. In general, a typical coal DTG curve is characterized by a mass loss peak in the low temperature region with a maximum rate occurring between 40 and 100 °C. This corresponds to the initial release of absorbed moisture. In some cases an adjacent peak to the absorbed moisture peak is observed, which has been attributed to either the release of crystal water associated with inherent minerals or chemically bonded moisture.28,39 The main devolatilization zone (>300 °C) constitutes the release of tar, primary gases, secondary gases, and the subsequent formation of char.14,28 The extensive asymmetric nature of a typical DTG curve of coal devolatilization therefore makes the use of a single

2. EXPERIMENTAL SECTION 2.1. Coal Selection. Four coals were selected from South African collieries: three noncaking coals originating respectively from the Witbank no.2 (INY), 4 (UMZ), and 5 (G#5) seams and one coking coal (TSH) originating from the Venda-Pafuri sector of the Soutpansberg coalfields.50 The three noncoking coals were beneficiated samples from the respective mines, while coal TSH was density separated elsewhere.51 The choice of the four coals was based on their (1) similarity in bituminous rank, (2) varying vitrinite content, (3) relatively low ash content ( UMZ > INY > TSH. Further investigation of the respective TG results of each coal revealed the division of the overall TG profile into clear regions of mass loss i.e., the release of moisture, occluded gases, and secondary pyrolytic water at temperatures below 350 °C and a main devolatilization regime where the parent coal structure decomposes to form tars, gases, and char at temperatures exceeding 350 °C. From DTG results it was evident that all four coals exhibited a clear, primary peak of devolatilization confined to the temperature range between 350 and 600 °C, while a subsequent shift in peak temperature occurred for the maximum evolution rate of the main peak of coal TSH. This subsequent shift in peak temperature for coal TSH could be attributed to both its higher aromaticity (rank) as well as its extensive thermoplasticity, which could have resulted in repolymerization reactions becoming more favorable in the first

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors would like to acknowledge and thank the following parties: Sasol for financial support, Anglo Coal and Exxaro for providing the necessary coal samples, and Dr. Marion Carrier (University of Stellenbosch) for her valuable inputs regarding coal devolatilization kinetics. The work presented in this paper is based on the research supported by the South African Research Chairs Initiative of the Department of Science and Technology and National Research Foundation of South Africa (NRF). Any opinion, finding or conclusion or recommendation expressed in this material is that of the author(s), and the NRF does not accept any liability in this regard.



NOMENCLATURE

Abbreviations

d.a.f. = dry, ash-free basis d.b. = dry basis FSI = free swelling index m.m.f.b. = mineral matter free basis OBF = objective function QOF = quality of fit 930

dx.doi.org/10.1021/ef402124f | Energy Fuels 2014, 28, 920−933

Energy & Fuels

Article

Greek symbols

(18) Braun, R. L.; Burnham, A. K. Analysis of chemical reaction kinetics using a distribution of activation energies and simpler models. Energy Fuels 1987, 1 (2), 153−161. (19) Cai, J.; Liu, R. New distributed activation energy model: numerical solution and application to pyrolysis kinetics of some types of biomass. Bioresour. Technol. 2008, 99, 2795−2799. (20) Navarro, M. V.; Aranda, A.; Garcia, T.; Murillo, R.; Mastral, A. M. Application of the distributed activation energy model to blends devolatilization. Chem. Eng. J. 2008, 142, 87−94. (21) Solomon, P. R.; Hamblen, D. G.; Carangelo, R. M.; Serio, M. A.; Deshpande, G. V. General model of coal devolatilization. Energy Fuels 1988, 2 (4), 405−422. (22) Donskoi, E.; McElwain, D. L. S. Approximate modelling of coal pyrolysis. Fuel 1999, 78, 825−835. (23) Donskoi, E.; McElwain, D. L. S. Optimization of coal pyrolysis modelling. Combust. Flame 2000, 122, 359−367. (24) Heidenreich, C. A.; Yan, H. M.; Zhang, D. K. Mathematical modelling of pyrolysis of large coal particles − estimation of kinetic parameters for methane evolution. Fuel 1999, 78, 557−566. (25) Mani, T.; Murugan, P.; Mahinpey, N. Determination of distributed activation energy model kinetic parameters using simulated annealing optimization method for non-isothermal pyrolysis of lignin. Ind. Eng. Chem. Res. 2009, 48, 1464−1467. (26) Pitt, G. J. The kinetics of the evolution of volatile products from coal. Fuel 1962, 41, 267−274. (27) Burnham, A. K.; Schmidt, B. J.; Braun, R. L. A test of the parallel reaction model using kinetic measurements on hydrous pyrolysis residues. Org. Geochem. 1995, 23 (10), 931−939. (28) Alonso, M. J. G.; Alvarez, D.; Borrego, A. G.; Menéndez, R.; Marbán, G. Systematic effects of coal rank and type on the kinetics of coal pyrolysis. Energy Fuels 2001, 15 (2), 413−428. (29) Miura, K. A new and simple method to estimate f(E) and k0(E) in the distributed activation energy model from three sets of experimental data. Energy Fuels 1995, 9 (2), 302−307. (30) Mianowski, A.; Radko, T. Isokinetic effect in coal pyrolysis. Fuel 1993, 72 (11), 1537−1539. (31) Solomon, P. R.; Fletcher, T. H.; Pugmire, R. J. Progress in coal pyrolysis. Fuel 1993, 72, 587−597. (32) Niksa, S.; Kerstein, A. R. FLASHCHAIN theory for rapid coal devolatilization kinetics. 1. Formulation. Energy Fuels 1991, 5, 647−665. (33) Grant, D. M.; Pugmire, R. J.; Fletcher, T. H.; Kerstein, A. R. Chemical model of coal devolatilization using percolation lattice statistics. Energy Fuels 1989, 3, 175−186. (34) De Souza-Santos, M. L. Solid fuels, combustion and gasification: Modelling, simulation and equipment operation; Marcel Dekker: New York, 2005; 462p. (35) Niksa, S.; Kerstein, A. R. On the role of macromolecular configuration in rapid coal devolatilization. Fuel 1987, 66, 1389−1399. (36) Aboulkas, A.; El Harfi, K.; Nadifiyine, M.; Benchanaa, M. Pyrolysis behaviour and kinetics of Moroccan oil shale with polystyrene. J. Pet. Gas Eng. 2011, 2 (6), 108−117. (37) Aboyade, A. O.; Hugo, T. J.; Carrier, M.; Meyer, E. L.; Stahl, R.; Knoetze, J. H.; Görgens, J. F. Non-isothermal kinetic analysis of the devolatilization of corn cobs and sugar cane bagasse in an inert atmosphere. Thermochim. Acta 2011, 517, 81−89. (38) Grønli, M. G.; Várhegyi, G.; Di Blasi, C. Thermogravimetric analysis and devolatilization kinetics of wood. Ind. Eng. Chem. Res. 2002, 41, 4201−4208. (39) Boiko, E. A. Research on kinetics of the thermal processing of brown coals of various oxidative aging degree using the non-isothermal method. Thermochim. Acta 2000, 348, 97−104. (40) Ahmaruzzaman, M.; Sharma, D. K. Non-isothermal kinetic studies on co-processing of vacuum residue, plastics, coal and petrocrop. J. Anal. Appl. Pyrolysis 2005, 73, 263−275. (41) Cai, J.; Wang, Y.; Zhou, L.; Huang, Q. Thermogravimetric analysis and kinetics of coal/plastic blends during co-pyrolysis in nitrogen atmosphere. Fuel Process. Technol. 2008, 89, 21−27.

β = heating rate σ = standard deviation of DAEM distribution ξi = fractional contribution of component, i

Roman symbols

Ea = activation energy Ea,i = activation energy of component i k0 = pre-exponential constant k0,i = pre-exponential constant of the ith pseudocomponent m0 = initial sample mass mf = final sample mass mt = logged mass at a certain time t Nk = number of heating rates Nm = number of DTG experimental points R = molar gas constant t = time T = temperature X = fractional conversion Xi = fractional conversion of the ith pseudocomponent Xt = overall fractional conversion of volatiles



REFERENCES

(1) Bunt, J. R.; Waanders, F. B. Identification of the reaction zones occurring in a commercial-scale Sasol-Lurgi FBDB gasifier. Fuel 2008, 87 (10−11), 1814−1823. (2) Miura, K. Mild conversion of coal for producing valuable chemicals. Fuel Process. Technol. 2000, 62, 119−135. (3) Smoot, L. D.; Pratt, D. T. Pulverised coal combustion and gasification; Plenum Press: New York, USA, 1979. (4) Smoot, L. D. Pulverised coal diffusion flames: A perspective through modelling. Eighteenth Symposium (International) on Combustion, The Combustion Institute 1981, 18 (1), 1185−1202. (5) Yue, C. Experimental and modelling study of pitch pyrolysis kinetics. Ph.D. thesis, University of British Columbia, Canada, 1995, 242p. (6) Trimm, D. L.; Turner, C. J. The pyrolysis of propane. I. Production of gases, liquids and carbon. J. Chem. Technol. Biotechnol. 1981, 31, 195− 204. (7) Kumar, P.; Kunzru, D. Modelling of naphtha pyrolysis. Ind. Eng. Chem. Process Des. Dev. 1985, 24, 774−782. (8) Conesa, J. A.; Marcilla, A.; Caballero, J. A.; Font, R. Comments on the validity and utility of the different methods for kinetic analysis of thermogravimetric data. J. Anal. Appl. Pyrolysis 2001, 58−59, 617−633. (9) Lázaro, M.-J.; Moliner, R.; Suelves, I. Non-isothermal versus isothermal technique to evaluate kinetic parameters of coal pyrolysis. J. Anal. Appl. Pyrolysis 1998, 47, 111−125. (10) Arenillas, A.; Rubiera, F.; Pevida, C.; Pis, J. J. A comparison of different methods of predicting coal devolatilization kinetics. J. Anal. Appl. Pyrolysis 2001, 58−59, 685−701. (11) Lu, G. Q.; Do, D. D. A kinetic model for coal reject pyrolysis at low heating rates. Fuel Process. Technol. 1991, 28, 35−48. (12) Badzioch, S.; Hawksley, P. G. W. Kinetics of thermal decomposition of pulverised coal particles. Ind. Eng. Chem. Process Des. Dev. 1970, 9, 521−527. (13) Kristiansen, A. Understanding coal gasification; IEA Coal Research: London, UK, 1996; 70 p. (14) Serio, M. A.; Hamblen, D. G.; Markham, J. R.; Solomon, P. R. Kinetics of volatile product evolution in coal pyrolysis: Experiment and theory. Energy Fuels 1987, 1 (2), 138−152. (15) Gürüz, G. A.; Ü çtepe, Ü .; Durusoy, T. Mathematical modelling of thermal decomposition of coal. J. Anal. Appl. Pyrolysis 2004, 71, 537− 551. (16) Lu, G. Q. Preparation and characterization of an effective sorbent for SO3 and NOx removal using coal washery reject. Ph.D. thesis, University of Queensland, Australia, 1991; 230p. (17) Sommariva, S.; Maffei, T.; Migliavacca, G.; Faravelli, T.; Ranzi, E. A predictive multi-step kinetic model of coal devolatilization. Fuel 2010, 89, 318−328. 931

dx.doi.org/10.1021/ef402124f | Energy Fuels 2014, 28, 920−933

Energy & Fuels

Article

(42) Dhumal, S. S.; Saha, R. K. Application of genetic algorithm for evaluation of kinetic parameters of coal pyrolysis. J. Energy Environ. 2006, 5, 112−124. (43) Sutcu, H. Pyrolysis by thermogravimetric analysis of blends of peat with coals of different characteristics and biomass. J. Chin. Inst. Chem. Eng. 2007, 38, 245−249. (44) Wiktorrson, L.-P.; Wanzl, W. Kinetic parameters for coal pyrolysis at low and high heating rates- a comparison of data from different laboratory equipment. Fuel 2000, 79, 701−716. (45) Zhang, C.; Jiang, X.; Wei, L.; Wang, H. Research on pyrolysis characteristics and kinetics of super fine and conventional pulverized coal. Energy Convers. Manage. 2007, 48, 797−802. (46) Cadle, A. B.; Cairncross, B.; Christie, A. D. M.; Roberts, D. L. The Karoo basin of South Africa: type basin for the coal-bearing deposits of southern Africa. Int. J. Coal Geol. 1993, 23, 117−157. (47) Cai, H.-Y.; Kandiyoti, R. Effect of changing inertinite concentration on pyrolysis yields and char reactivities of two South African coals. Energy Fuels 1995, 9 (6), 956−961. (48) Cairncross, B. An overview of the Permian (Karoo) coal deposits of southern Africa. African Earth Sci. 2001, 33, 529−562. (49) Aboyade, A. O.; Carrier, M.; Meyer, E. L.; Knoetze, J. H.; Görgens, J. F. Model fitting kinetic analysis and characterisation of the devolatilization of coal blends with corn and sugarcane residues. Thermochim. Acta 2012, 530, 95−106. (50) Pinheiro, H. J. A techno-economic and historical review of the South African coal industry in the 19th and 20th centuries and analyses of coal product samples of South African collieries 1998−1999. In Bulletin 113. SABS: Pretoria, 1999; 97p. (51) Hattingh, B. B.; Everson, R. C.; Neomagus, H. W. J. P.; Bunt, J. R.; Van Niekerk, D.; Jordaan, J. H. L.; Mathews, J. P. Elucidation of the structural and molecular properties of typical South African coals. Energy Fuels 2013, 27 (6), 3161−3172. (52) ISO Standard 7404-3, 1994. Methods for the petrographic analysis of coals - Part 3. Method of determining maceral group composition, ISO Standards, Geneva. (53) ISO Standard 7404-5, 1994. Methods for the petrographic analysis of coals - Part 3. Method of determining microscopically the reflectance of vitrinite, ISO Standards, Geneva. (54) ISO Standard 12902, 2007. Solid mineral fuels - Determination of total carbon, hydrogen and nitrogen - Instrumental methods, ISO Standards, Geneva. (55) SABS ISO Standard 1171, 1997. Standards South Africa, Pretoria. ISBN 978-0-626-24518-4, 13p. (56) SABS ISO Standard 562, 2011:2. Standards South Africa, Pretoria. ISBN 978-0-626-24517-7, 17p. (57) Antal, M. J.; Várhegyi, G. Cellulose pyrolysis kinetics: The current state of knowledge. Ind. Eng. Chem. Res. 1995, 34, 703−717. (58) Yang, H.; Chen, H.; Ju, F.; Yan, R.; Zhang, S. Influence of pressure on coal pyrolysis and char gasification. Energy Fuels 2007, 21, 3165− 3170. (59) Abhari, R.; Isaacs, L. L. Drying kinetics of lignite, subbituminous coals and high-volatile bituminous coals. Energy Fuels 1990, 4, 448−452. (60) Biagini, E.; Lippi, F.; Petarca, L.; Tognotti, L. Devolatilization rate of biomasses and coal-biomass blends: an experimental investigation. Fuel 2002, 81, 1041−1050. (61) Brown, M. E. Introduction to thermal analysis: techniques and applications, 2nd ed.; Kluwer Academic Publishers: London, UK, 2001; 280p. (62) Várhegyi, G. Aims and methods in non-isothermal reaction kinetics. J. Anal. Appl. Pyrolysis 2007, 79, 278−288. (63) Vyazovkin, S. Computational aspects of kinetic analysis. Part C. The ICTAC kinetics project − the light at the end of the tunnel? Thermochim. Acta 2000, 355, 155−163. (64) White, J. E.; Catallo, W. J.; Legendre, B. L. Biomass pyrolysis kinetics: a comparative critical review with relevant agricultural residue case studies. J. Anal. Appl. Pyrolysis 2011, 91, 1−33. (65) Bellais, M. Modelling of the pyrolysis of large wood particles. Ph.D. Thesis, The Royal Institute of Technology (KTH), Stockholm, Sweden, 2007; 105p.

(66) Li, X.; Song, H.; Wang, Q.; Meesri, C.; Wall, T.; Yu, J. Experimental study on drying and moisture re-adsorption kinetics of an Indonesian low rank coal. J. Environ. Sci. 2009, 21 (Suppl. 1), S127− S130. (67) Yip, K.; Wu, H.; Zhang, D-.K. Mathematical modelling of Collie coal pyrolysis considering the effect of steam produced in situ from coal inherent moisture and pyrolytic water. Proc. Combust. Inst. 2009, 32, 2675−2683. (68) Caballero, J. A.; Conesa, J. A. Mathematical considerations for non-isothermal kinetics in thermal decomposition. J. Anal. Appl. Pyrolysis 2005, 73, 85−100. (69) Brazier, D. W.; Nickel, G. H. Thermo-analytical methods in vulcanizate analysis II. Derivative thermogravimetric analysis. Rubber Chem. Technol. 1975, 48 (1), 661−677. (70) Brazier, D. W.; Schwartz, N. V. The effect of heating rate on the thermal degradation of polybutadiene. J. Appl. Polym. Sci. 1978, 22 (1), 113−124. (71) Sircar, A. K.; Lamond, T. G. Identification of elastomers in tire sections by total thermal analysis. I. Tread and black sidewall. Rubber Chem. Technol. 1975, 48 (1), 301−309. (72) Yang, J.; Kaliaguine, S.; Roy, C. Improved quantitative determination of elastomers in tire rubber by kinetic simulation of DTG curves. Rubber Chem. Technol. 1993, 66, 213−229. (73) Conesa, J. A.; Caballero, J. A.; Marcilla, A.; Font, R. Analysis of different kinetic models in the dynamic pyrolysis of cellulose. Thermochim. Acta 1995, 254, 175−192. (74) Roduit, B. Computational aspects of kinetic analysis. Part E: The ICTAC Kinetics Project-numerical techniques and kinetics of solid state processes. Thermochim. Acta 2000, 355, 171−180. (75) Seo, D. K.; Park, S. S.; Kim, Y. T.; Hwang, J.; Yu, T-.U. Study of coal pyrolysis by thermo-gravimetric analysis (TGA) and concentration measurements of the evolved species. J. Anal. Appl. Pyrolysis 2011, 92, 209−216. (76) Solomon, P. R.; Serio, M. A.; Suuberg, E. M. Coal pyrolysis experiments: Kinetic rates and mechanisms. Prog. Energy Combust. Sci. 1992, 18, 133−220. (77) Wang, S.; Tang, Y.; Schobert, H. H.; Mitchell, G. D.; Liao, F.; Liu, Z. A thermal behaviour study of Chinese coals with high hydrogen content. Int. J. Coal Geol. 2010, 81, 37−44. (78) Elder, J. P.; Harris, M. B. Thermogravimetry and differential scanning calorimetry of Kentucky Bituminous coals. Fuel 1984, 63 (2), 262−267. (79) Howard, H. C. Fundamentals of coal pyrolysis and hydro pyrolysis. In The chemistry of coal utilization. Second supplementary volume; Elliott, M. A., Ed.; Wiley: New York, USA, 1981; 2374p. (80) Kabe, T.; Ishihara, A.; Qian, E. W.; Sutrisna, I. P.; Kabe, Y. Coal and coal related compounds: Structures, reactivity and catalytic reactions; Kodansha Ltd.: Amsterdam, Elsevier: Tokyo, 2004; 350p. (81) Ladner, W. R. The products of coal pyrolysis: properties, conversion and reactivity. Fuel Process. Technol. 1988, 20, 207−222. (82) Singh, S.; Wu, C.; Williams, P. T. Pyrolysis of waste materials using TGA-MS and GA-FTIR as complementary characterisation techniques. J. Anal. Appl. Pyrolysis 2012, 94, 99−107. (83) Vuthaluru, H. B. Thermal behaviour of coal/biomass blends during co-pyrolysis. Fuel Process. Technol. 2004, 85 (2−3), 141−155. (84) Seo, D. K.; Park, S. S.; Hwang, J.; Yu, T-.U. Study of the pyrolysis of biomass using thermo-gravimetric analysis (TGA) and concentration measurements of the evolved species. J. Anal. Appl. Pyrolysis 2010, 89, 66−73. (85) Williams, P. T.; Besler, S. The influence of temperature and heating rate on the slow pyrolysis of biomass. Renewable Energy 1996, 7 (3), 233−250. (86) Caballero, J. A.; Conesa, J. A. New approach to thermal analysis kinetics by considering several first order reactions. Thermochim. Acta 2011, 525, 40−49. (87) Van Heek, K. H.; Hodek, W. Structure and pyrolysis behaviour of different coals and relevant model substances. Fuel 1994, 73 (6), 886− 896. 932

dx.doi.org/10.1021/ef402124f | Energy Fuels 2014, 28, 920−933

Energy & Fuels

Article

(88) Branca, C.; Albano, A.; Di Blasi, C. Critical evaluation of global mechanisms of wood devolatilization. Thermochim. Acta 2005, 429, 133−141. (89) Lázaro, M.-J.; Moliner, R.; Suelves, I. Non-isothermal versus isothermal technique to evaluate kinetic parameters of coal pyrolysis. J. Anal. Appl. Pyrolysis 1998, 47, 111−125. (90) Li, X.; Song, H.; Wang, Q.; Meesri, C.; Wall, T.; Yu, J. Experimental study on drying and moisture re-adsorption kinetics of an Indonesian low rank coal. J. Environ. Sci. 2009, 21 (Suppl. 1), S127− S130. (91) Anthony, D. B.; Howard, J. B. Coal devolatilization and hydrogasification. J. Am. Inst. Chem. Eng. 1976, 22, 625−656. (92) Jüntgen, H.; Van Heek, K. H. An update of German nonisothermal coal pyrolysis work. Fuel Process. Technol. 1979, 2 (4), 261− 293. (93) Suuberg, E. Rapid pyrolysis and hydropyrolysis of coal. D.Sc. Dissertation, Massachusetts Institute of Technology, USA, 1977, 396p. (94) Fitzer, E.; Mueller, K.; Schaefer, W. Conversion of organic compounds to carbon. In Chemistry and physics of carbon; Walker, P., Ed.; Marcel Dekker: New York, USA, 1971; Vol. 7, pp 238−283. (95) Biswas, S.; Sharma, D. K. Synergistic co-processing/co-cracking of Jatropha oil, petroleum vacuum residue, and high density polyethylene. J. Renewable Sustainable Energy 2012, 4, 043112. (96) Biswas, S.; Sharma, D. K. Studies on cracking of Jatropha oil. J. Anal. Appl. Pyrolysis 2013, 99, 122−129. (97) Sadhukhan, A. K.; Gupta, P.; Saha, R. K. Modelling and experimental investigations on the pyrolysis of large coal particles. Energy Fuels 2011, 25, 5573−5583. (98) Sharma, S.; Ghoshal, A. K. Study of kinetics of co-pyrolysis of coal and waste LDPE blends under argon atmosphere. Fuel 2010, 89, 3943− 3951. (99) Zhang, K.; You, C.; Li, Y. Experimental and numerical investigation on the pyrolysis of single course lignite particles. Korean J. Chem. Eng. 2012, 29 (4), 540−548. (100) Joseph, J. T.; Fisher, R. B.; Masin, C. A.; Dyrkacz, G. R.; Bloomquist, C. A.; Winans, R. E. Coal maceral chemistry. 1. Liquefaction behaviour. Energy Fuels 1991, 5, 724−729. (101) Zhao, Y.; Hu, H.; Jin, L.; He, X.; Wu, B. Pyrolysis behaviour of vitrinite and inertinite from Chinese Pingshuo coal by TG-MS and in a fixed bed reactor. Fuel Process. Technol. 2011, 92, 780−786.

933

dx.doi.org/10.1021/ef402124f | Energy Fuels 2014, 28, 920−933