Thermal Degradation Kinetic Study of Rubber Seed Oil and Its Methyl

Aug 24, 2017 - The obtained overall average values of activation energy for RSO thermal decomposition are in good agreement with those of sunflower oi...
0 downloads 10 Views 2MB Size
Subscriber access provided by University of Texas Libraries

Article

Thermal degradation kinetic study of rubber seed oil and its methyl esters under inert atmosphere Ali Shemsedin Reshad, Pankaj Tiwari, and Vaibhav V. Goud Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b02249 • Publication Date (Web): 24 Aug 2017 Downloaded from http://pubs.acs.org on August 29, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

1

Thermal degradation kinetic study of rubber seed oil and its methyl esters under

2

inert atmosphere

3

Ali Shemsedin Reshad, Pankaj Tiwari*, Vaibhav V. Goud*

4 5

Department of Chemical Engineering, Indian Institute of Technology Guwahati, Guwahati,

6

Assam, 781039, India

7

∗Corresponding

8

E-mail address: [email protected] (Pankaj Tiwari), [email protected] (V. V. Goud)

9

ABSTRACT: Non-edible vegetable oil feedstocks are promising for sustainable production of

10

biodiesel. Thermal decomposition characteristics of the feedstocks and its biodiesel are crucial

11

for handling and quality control. Thermal degradation of rubber seed oil (RSO) and rubber seed

12

oil methyl esters (ROME) were investigated with the help of thermogravimetry. The samples

13

were pyrolysed from 30 °C to 800 °C for heating rates of 10 °C/min to 50 °C/min with 10

14

°C/min increment under nitrogen atmosphere. The temperature window for thermal degradation

15

of RSO and ROME were shifted towards higher range as the heating rate increased from 10

16

°C/min to 50 °C/min. Transesterification reaction leads to decrease the molecular weight of

17

triglycerides present in sample (RSO) and this causes to lower the thermal stability of the

18

produced product (ROME). Fourier transform infrared (FT-IR) analysis of evolved gaseous

19

products during pyrolysis revealed the formation of water, carbon dioxide, carbon monoxide,

20

saturated (alkanes), and unsaturated (alkenes) aliphatic hydrocarbons. Friedman (FRD), Flynn-

21

Wall-Ozawa (FWO), modified Coat-Redfern (MCR), Kissinger (KM) methods and Avrami

22

theory were applied to calculate the values of activation energy (E), order of reaction (n) and

author: Tel.: +91 361 2582263/2272; fax: +91 361 2582291

1 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 35

23

enthalpy (∆H). Furthermore, pre-exponential factor (A), entropy (∆S) and Gibbs free energy

24

(∆G) were also calculated.

25

INTRODUCTION

26

Biodiesel has gained worldwide attention to partially substitute fossil based diesel fuel.

27

Chemically, biodiesel is a mixture of long chain fatty acid alkyl esters derived from triglyceride

28

present in vegetable oil, animal fat and waste cooking oil through transesterification process.

29

Transesterification process improves the physico-chemical, thermal and flow properties of the

30

feedstock. Soybean oil, sunflower oil, palm oil, coconut oil, corn oil, rapeseed oil and olive oil

31

are the most widely used first generation feedstocks for the production of biodiesel. First

32

generation biodiesel feedstocks are usually categorized as a part of food chain and account

33

around 60%–80% of the total biodiesel production cost

34

biodiesel production can be overcome by using of non-edible feedstocks. Jatropha curcas L.

35

(Jatropha)

36

pinnata L. (Karanja), Calophyllum inophyllum L. (Polanga) 7, Croton megalocarpus (Musine) 7,

37

Cocos nucifera (coconut)

38

materials for biodiesel synthesis.

1, 2

, Ricinus communis (Castor)

7

and

3, 4

1, 2

. The challenge for higher cost of

, Hevea brasiliensis (Rubber tree)

Mesua ferrea (Nahor)

1, 4

5, 6

, Pongamia

are sustainable non-edible raw

39

Beside edibility of the feedstock, thermal degradation characteristics of both the feedstock

40

and product are of great concern for scientific applications 8-11. Thermal degradation analysis can

41

be used to estimate the thermal properties such as activation energy, enthalpy, Gibbs free energy,

42

entropy, heat capacity as well as the quality of the produced biodiesel

43

(TG) technique monitors the physical and chemical changes of sample happens with

44

temperature. Thermal degradation of vegetable oils and their methyl esters mainly involves 2 ACS Paragon Plus Environment

12, 13

. Thermogravimetry

Page 3 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

45

physico-chemical processes of volatilization and decomposition. The gaseous products evolved

46

during thermal decomposition of the sample can be identified in real time by coupling Fourier

47

transform infrared spectroscopy (FT-IR) with TG

48

such coupled technique has been used in various research fields to estimate sample structure and

49

composition. Li et al.

50

behavior of biodiesel samples derived from peanut oil, palm oil and waste cooking oil. TGA

51

results showed that the onset decomposition and peak temperatures for palm oil methyl esters are

52

higher than both peanut oil and waste cooking oil methyl esters due to lower content of

53

unsaturated fatty acids components in palm oil methyl esters. The real time analysis of evolved

54

products using attached FT-IR revealed the formation of alkanes, cyclic and aromatic

55

compounds along with CO, CO2 and H2O. Santos et al.13 reported thermal decomposition of

56

sunflower oil and its biodiesel using non-isothermal thermogravimetric analysis under inert

57

atmosphere. The values of kinetic parameters have been estimated as 155.62–200.12kJ/mol (E)

58

and 0.95–1.82 (n) for sunflower oil and 61.32kJ/mol–115.35kJ/mol (E) and 0.69–1.89 (n) for its

59

biodiesel. Souza et al.

60

under air and nitrogen atmospheres. The values of activation energy (E) for cotton oil have been

61

found higher than that of cotton oil based biodiesel for both air and nitrogen atmosphere.

10

15

10, 14

. The analysis of evolved gas by using

employed thermogravimetry analysis (TGA) to estimate thermal

evaluated thermal and kinetic behavior of cotton oil and its biodiesel

62

The investigations on thermal degradation of rubber seed oil and its biodiesel have been

63

rarely reported. In the present work, thermal degradation behavior of rubber seed oil and its

64

biodiesel under nitrogen atmosphere was studied. There are varieties of kinetic methods

65

available in open literature to deduce the kinetic parameters. However, Friedman model (FRD),

66

Flynn-Wall-Ozawa (FWO), and Coats-Redfern (CR) kinetic methods are consider more reliable

67

and widely used. Hence, in the present study, Friedman (FRD), Flynn-Wall-Ozawa (FWO), 3 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 35

68

Modified Coats-Redfern (MCR) and Kissinger (KM) kinetics methods were applied to estimate

69

the value of activation energy (E) and enthalpy (∆H) for thermal degradation reactions.

70

MATERIALS AND METHODS

71

Materials. Rubber seeds, collected from Assam India were de-shelled manually and the

72

obtained kernels were subjected for oil extraction. Rubber seed oil (RSO) and rubber seed oil

73

methyl esters (ROME) were obtained through Soxhlet extractor and ultrasonic- assisted

74

transesterification, respectively. The detail can be found in our previous studies

75

extracted oil and produced biodiesel under optimum conditions were used for thermal analysis.

76

The physico-chemical properties of obtained RSO and ROME are presented in Table 1

77

13

78

(biodiesel) are depicted in Fig. S1 and Fig. S2, respectively. The signal at 69 ppm and 62 ppm in

79

13

80

O–) and 62 ppm (CH2–C–O–)) (Fi. S1) while the signal is absent in the biodiesel (ROME) (Fig.

81

S2). It can be clearly seen that the glyceride backbone of triglyceride is totally absent in the

82

ROME sample. The methoxy carbon of methyl esters of ROME illustrates the signal at 51.49

83

ppm. The unsaturation signal (–C=C–) obtained between 133–120 ppm in

84

the

85

(monounsaturated) fatty acids and ester in RSO (Fig. S1) and ROME (Fig. S2), respectively.

16, 17

. The

16

. The

C Nuclear magnetic resonance (NMR) spectra of rubber seed oil and its methyl esters

C NMR spectrum of rubber seed oil are due to the carbonyl methylene groups (69 ppm (H–C–

presence

of

linoleic

(polyunsaturated),

linolenic

13

C NMR are due to

(polyunsaturated)

and

oleic

86

Methods. Thermogravimetric analysis. Thermal degradation of rubber seed oil (RSO) and

87

rubber seed oil methyl esters (ROME) were evaluated using TG analyzer (Netzch STA449F300)

88

at various heating rates (10, 20, 30 40 and 50 °C/min) under nitrogen atmosphere. The samples

89

were weighed ~10 mg in alumina crucible. TGA experiments were conducted from 30 °C to 800 4 ACS Paragon Plus Environment

Page 5 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

90

°C for respective heating rates, where the nitrogen gas (99.999% purity) flow rate was set at 60

91

ml/min. The evolved gas products during thermal decomposition at heating rate of 40 °C/min

92

were continuously monitored and measured using Perkin-Elmer TGA and FT-IR coupled system.

93

Gas transfer tube and gas cell were heated up to 250 °C to prevent condensation of evolved

94

products.

95

Thermal degradation study. The values of activation energy (E) for thermal decomposition of 18

RSO and ROME were calculated using Friedman (FRD)

97

Kissinger (KM)

98

rate, considering Arrhenius temperature dependency for constant heating rate can be expressed

99

by the following Eq.1.

20

, Flynn-Wall-Ozawa (FWO)

19

96

and Modified Coats-Redfern (MCR) methods

(

11, 21

,

. Thermal decomposition

)

100

β ⋅ dα

101

Where β, f(α), A, α, E, T and R refer to heating rate (°C/min), reaction mechanism model, pre-

102

exponential factor (1/min), degree of decomposition of the samples (conversion), activation

103

energy (J/mol), temperature (K) and gas constant (8.314 J/mol K), respectively. The value of α at

104

appropriate time/temperature can be calculated using TG data (Eq.2)

105

α=

106

Where: wo, wt and wf are initial weight, weight at time t, and final weight, respectively.

107

dT

= f (α ) ⋅ A exp − E

( w0 − wt )

RT

(Eq.1)

(Eq.2)

( w0 − w f )

Friedman model (FRD) (Eq.3) is the first and general isoconversional method on the basis of

108

model free differential technique, and obtained by taking natural logarithm both side of Eq. 1

109

and becomes (Eq. 3).

5 ACS Paragon Plus Environment

18

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(

)

Page 6 of 35

110

ln dα

111

The aspects of Friedman 18 and Coats-Redfern 21 methods can be combined for estimation of

112

kinetic parameters for multiple heating rates. The general expression modified Coats-Redfern for

113

nth order is as follows (Eq. 4) 22-24;

114

  −E ln  β 2  = ln A ⋅ R E ⋅ g ( ) R ⋅ Ti, j α  T i , j i, j

dt i, j

= ln ( f (α ) ⋅ A ) − E

(

(Eq.3)

RTi, j

)

(Eq.4)

115

Where: g(α) is the integral form of the reaction model.

116

Beside the differential approaches, the fundamental rate expression (Eq.1) can also be used

117

by integral method of Flynn-Wall-Ozawa model (Eq. 5) using Doyle approximation to estimate

118

the kinetic parameters 19.

119

ln βi, j = C − 1.052 E

(Eq.5)

RTi, j

120

Kissinger proposed that the maximum rate occurs when d(dα/dt)/dt is zero. Therefore, the

121

differentiation of fundamental Arrhenius expression (Eq. 1) for constant heating rate at which

122

maximum rate occurs (at peak temperature Tmax) is equal to zero. The simplified Kissinger model

123

for first order thermal decomposition at peak temperature (Tmax) expression is as follow (Eq. 6);

124

 β   AR   E   1  ln  = ln     −   T2  E    R   Tmax  j  max  j

(Eq.6)

125

Based on the same degree of thermal degradation (i) at different heating rates (j), linear plots

126

of ln(dα/dt)i,j versus 1/Ti,j (FRD, Eq. 3), ln(βj/T2i,j) versus 1/Ti,j (MCR, Eq. 4), ln(βj) versus 1/Ti,j

6 ACS Paragon Plus Environment

Page 7 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

127

(FWO, Eq. 5) and ln(β/Tmax2) versus 1/Tmaxj (KM, Eq. 6) were constructed. The slope of the

128

straight line was used to calculate the value of activation energy.

129

The thermodynamic parameters such as enthalpy (∆H) (Eq.7), Gibbs free energy (∆G) and

130

entropy (∆S) of RSO and ROME samples were calculated at the maximum peak temperatures 10,

131

13, 25

132

calculated using Eq. 8 10, 25. The obtained pre-exponential factors were used to calculate ∆G (Eq.

133

9) and ∆S (Eq.10) 26-28.

134

∆ H = E − RT

135

A=

136

∆G = E + RT ln 

137

∆S =

. For this purpose , the value of pre-exponential factor (A) for RSO and ROME were

(Eq.7)

β ⋅E

 E  ⋅ exp    RT  R ⋅T 2

(Eq.8)

 KB ⋅T    h⋅ A 

(Eq.9)

∆ H − ∆G

(Eq.10)

T

138

Where KB is Boltzmann constant (1.3806×10-23 m2 kg s-2 K-1) and h is Planck constant

139

(6.6261×10-34 m2 kg s-1).

140

Beside the values of activation energy and pre-exponential factor of thermal degradation,

141

reaction order is also an important index. Avrami theory was applied to calculate the order of

142

thermal degradation of both RSO and ROME at various temperatures using (Eq.11) 10, 13, 29:

143

α = 1 − exp 

144

Eq.11 can be simplified to the following expression (Eq.12);

 A ⋅ exp ( − E / RT )    βn  

(Eq.11)

7 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

E − n ln β RT

Page 8 of 35

145

ln(− ln(1 − α )) = ln A −

146

At a particular degradation temperature, T, the points of ln(-ln(1-α)) versus lnβ at various heating

147

rates could be fitted to a linear line. The reaction order (n) can be calculated from the slope of the

148

linear line.

(Eq.12)

149

RESULTS AND DISCUSSION

150

TGA analysis. Thermal decomposition behavior of both, RSO and ROME samples was

151

investigated. TGA and DTG profiles for RSO and ROME at various heating rates are shown Fig.

152

1. Volatilization occurs during early stage, when the lighter components evolved. During

153

decomposition the heavier components break to low molecular weight components. Further, the

154

evolved products, low molecular weight components during decomposition go through the

155

volatilization process. Thus, in active thermal degradation stage both the phenomena take place.

156

A change in the slope of TG profile was considered as the beginning of new stage. However,

157

only active thermal degradation stage was subjected for kinetic analysis. Therefore, overall 3

158

stages were found for both the samples, RSO and ROME (Fig. 1). Active pyrolysis stages for

159

RSO (stage I, IIa and IIb) and for ROME (stage I and II) were selected for kinetic analysis. The

160

two split peaks of RSO (stage IIa and IIb) and ROME (stage Ia and Ib) were considered as single

161

stage (II for RSO and stage I for ROME) for kinetic analysis. The main (active) thermal

162

decomposition of ROME sample at all the heating rates considered (10, 20, 30, 40 and 50

163

°C/min) occurred in single stage, that describes the decomposition and volatilization

164

phenomenas. However, the two-stage thermal decomposition was observed for RSO sample

165

(Table 2). The first stage (RSO-I) showed the presence of higher free fatty acids in the sample. In

166

addition to free fatty acids decomposition, mass loss due to moisture removal and degradation of 8 ACS Paragon Plus Environment

Page 9 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

167

light volatile compounds also takes place within this stage. The active decomposition of RSO in

168

second stage (RSO-II) is due to the degradation and volatilization of triglyceride. At a heating

169

rate of 20 °C/min, the thermal decomposition of RSO and its methyl ester (ROME) occurred in

170

the temperature range of 242–480 °C (95.7 wt%) and 165–456 °C (97.2 wt%), respectively. The

171

active thermal decomposition and volatilization of RSO sample were started relatively at higher

172

temperature values, started around 242 °C and completed around 480 °C compared to its methyl

173

ester (ROME) (i.e. started around 165 °C and completed around 456 °C). This is due to the fact

174

that RSO has higher molecular weight compounds and stronger intermolecular force (higher

175

viscosity) as compared to ROME. The value of onset temperature for ROME (165 °C) obtained

176

in the present study is proximity similar with palm oil methyl ester (164.5 °C), peanut oil methyl

177

ester (155.8 °C), waste cooking oil methyl ester (142.2 °C) as reported by Li et al.10. The onset

178

temperature values for palm oil methyl ester (164.5 °C) 10, peanut oil methyl ester (155.8 °C) 10,

179

waste cooking oil methyl ester (142.2 °C) 10 were approximately similar with ROME (165 °C).

180

It can be seen from Fig. 1(A–D) that the onset temperature and the temperature at which the

181

rate of mass loss is maximum (Tmax) were shifted towards higher temperatures with increasing

182

heating rate. This is due to low heat distribution (heat transfer limitation) 10, 22, 30. With respect to

183

heat transfer phenomenon, the initial thermal degradation temperature values for rubber seed oil

184

in first stage (RSO-I) and second stage (RSO-II), and for rubber seed oil methyl esters (ROME)

185

were shifted from 223 °C to 264 °C, 253 °C to 295 °C and 120 °C to 193 °C as the heating rate

186

increased from 10 °C/min to 50 °C/min, respectively (Table 2). Furthermore, the values of peak

187

temperature, Tmax were also changed from 271 °C to 328 °C, 416 °C to 446 °C and 250 °C to

188

287.9 °C for thermal degradation of RSO-I, RSO-II and ROME, respectively. Similarly, some

189

other parameters such as Tf and wmax values were also increased (Table 3). From Fig. 1(B and D), 9 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

190

it can also be clearly observed that the heating rate has a significant effect on the rate of thermal

191

degradation of the samples.

Page 10 of 35

192

Fourier transform infrared (FT-IR) analysis of evolved products. FT-IR spectra obtained

193

for RSO and ROME were found with similar characteristics due to similar nature of chemical

194

structure of functional group present in the samples (Fig. S3, supplementary data). However, the

195

signals specific to hydroxyl group of free fatty acid can be observed only in the spectrum of RSO

196

at 3480 cm-1. In addition, single peak at 1456 cm-1 was observed for bending vibration of CH2

197

and CH3 group only in RSO sample. Furthermore, a signal peak specific to ester functional group

198

for triglycerides for RSO and methyl esters of ROME was clearly observed at wave number

199

1735–1740 cm-1. Fig. 2 shows FT-IR signatures of evolved products at various mass loss

200

temperatures during thermal volatilization and decomposition of RSO and its methyl esters for

201

heating rate of 40 °C/min. Absorbance peaks corresponding to gaseous and liquid water

202

molecules can be seen only for RSO thermal degradation at 3500–3950 cm−1 and 3400–

203

3500 cm−1, respectively

204

oxygen-containing group in ROME (i.e. R1COOR2) mainly decomposed into C=O– and C–O 10,

205

12

206

clearly shown in Fig. 2. Symmetric and asymmetric stretching vibrations of –CH– and, –CH3

207

asymmetric deformation vibration in the range of wave number, 3000–2700 cm-1 and 1475–1000

208

cm-1, respectively revealed the presence of alkanes in the evolved products during thermal

209

degradation of RSO and ROME. Carbonyl groups of aldehydes and ketones were also observed

210

as H–C=O– and –C=O– in plane bending vibrations at 1720–1740 cm-1 and 1735–1750 cm-1

211

(Fig. 2), respectively. Furthermore, C–O–C (stretching vibration at 1000–1300 cm-1), C=O

212

(bending vibration at 2250–2400 cm-1 and 580–730 cm-1), C=O (stretching vibration at 2000–

31

. This shows the moisture content for RSO higher than ROME and

. The characteristics infrared absorption peaks for volatile components functional group were

10 ACS Paragon Plus Environment

Page 11 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

213

2250 cm-1) and C–O (stretching vibration at 2200–2100 cm-1) were observed in the evolved

214

gaseous product (Fig. 2) 14.

215

Taking the absorbance of identified volatile compounds such as alkanes, alkenes,

216

aldehydes, ketones, ethers and CO2, the intensity of the evolved compounds with increasing

217

thermal decomposition temperature is presented in Fig. 3. As the thermal decomposition

218

temperature increased, ether and aldehydes were formed due to de-oxygenation of ester12.

219

Presence of CO2 in the evolved product revealed that the de-carboxylation of ester and

220

triglyceride occurred during the thermal decomposition of ROME and RSO.10. It can be seen

221

from Fig. 3 that the maximum absorbance characteristic peaks for alkenes, alkanes, aldehydes

222

and ether occurred at same temperature (448 °C) and, for ketones and CO2 at temperature of 536

223

°C during thermal degradation of RSO. Similarly, during the thermal degradation of ROME, the

224

maximum rate of alkenes, alkanes, aldehydes and ether production occurred at temperature of

225

307 °C while for ketones and CO2, the maximum rate values were found at 571 °C. As can be

226

observed in TGA profiles (Fig. 1) and FT-IR spectra (Fig. 2), beyond 500°C, the mass losses

227

during thermal decomposition of ROME (1.44 wt%) and RSO (1.19 wt%) are mainly due to the

228

formation of CO2 and ketone (Fig. 3). Taking the Lambert-Beer law into consideration, the

229

concentrations of alkanes in evolved products were found maximum during the thermal

230

decomposition of RSO and ROME

231

ROME with temperature (Fig. 3) were found similar with that of DTG curves (Fig. 1).

12, 14

. The appearance of absorbance profiles of RSO and

232

Kinetic parameter calculation. It can be seen from Fig. 4 that the degree of conversion of

233

both the samples were greatly varied with temperature and heating rate. To estimate the

234

dependency of activation energy on temperature and degree of conversion during the active

235

decomposition (pyrolysis) process, nine conversion fractions from 0.1 to 0.9 with increment of 11 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

236

0.1 were selected at all heating rates. Based on FRD, FWO and MCR isoconversional methods,

237

the activation energies (E) for selected conversions were calculated from the slopes of linear

238

regression (Fig. 5–7).

Page 12 of 35

239

The deduced regression lines for RSO-I, RSO-II and ROME are presented in Fig. 5–7 and

240

the values of activation energies (E) are summarizes in Table 3. Parallel lines shown in Fig. 5–7

241

indicate that the values of activation energy (E) for thermal degradation of the samples (RSO and

242

ROME) for free fatty acids (RSO-I), triglycerides (RSO-II) and fatty acid methyl esters (ROME)

243

in the respective sample follow the same reaction rate or intensity. In other words, the values of

244

activation energy (E) obtained at different degree of conversions and temperatures were

245

proximity similar for RSO-I and ROME. This suggests that a single mechanism or unification of

246

multiple reactions mechanism was followed in the active thermal degradation of RSO-I and

247

ROME

248

to each other’s and this shows the change in activation energy at different degree of conversions

249

due to multiple and parallel reactions during thermal decomposition of triglycerides of RSO.

250

Change in the slope of the lines at different conversion show that the rates of thermal

251

decomposition differ due to multiple reactions occurring. The higher values of activation

252

energies were found at later stage of conversion (α=0.9). The values of R2 were found to be more

253

than 0.99 for the selected conversions (0.1 to 0.9) which show the fitness of the methods

254

considered, FRD, FWO and MCR.

255

The values of activation energy (E) for RSO-I and RSO-II vary from 69.8–92.7 kJ/mol and 144–

256

433.3 kJ/mol, respectively. The difference in the values of activation energy at lower (α=0.1) and

257

higher (α=0.9) degree of conversion for RSO-II reveals that incomplete decomposition of free

258

fatty acids occurred in first stage. Further, it also suggests that the thermal decomposition of

10, 11

. However, for RSO-II, fitted lines at different conversion were slightly not parallel

12 ACS Paragon Plus Environment

Page 13 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

259

RSO samples is a complex reaction which involves several parallel, competitive and consecutive

260

types of reactions. Activation energy (E) is the critical energy barrier to be overcome to generate

261

a chemical reaction and also it represent the minimum energy required to break the chemical

262

bond between atoms 29. Hence, the higher values of activation energy for RSO-II as compared to

263

RSO-I indicate that more difficult reactions have taken place during secondary stage (RSO-II).

264

This is maybe due to higher free fatty acids content in RSO sample. Considering all three

265

methods (FRD, FWO and MCR), the overall average values of activation energy (E) for thermal

266

decomposition of RSO were found to be 167.8 kJ/mol (FRD), 144.85 kJ/mol (FWO) and 142.05

267

kJ/mol (MCR). The obtained overall average values of activation energy for RSO thermal

268

decomposition are in good agreement with sunflower oil (170–210 kJ/mol)

269

160.2 kJ/mol)

270

kJ/mol)

271

RSO-I and RSO-II. The values of activation energy obtained by Kissinger method (E) (67.89

272

kJ/mol, RSO-I; 183.85 kJ/mol, RSO-II) are lower than that of obtained by FRD (76.7 kJ/mol,

273

RSO-I; 258.9 kJ/mol, RSO-II), FWO (80.9 kJ/mol, RSO-I; 208.8 kJ/mol, RSO-II) and MCR

274

(75.8 kJ/mol, RSO-I; 208.3 kJ/mol, RSO-II) methods.

11

11

, karanja seed oil (156.5–160.7 kJ/mol)

11

13

, soybean (146.6–

and mustard seed oil (142.4–148.1

. Fig. 6C shows the linear plots for Kissinger method for thermal decomposition of

275

Table 3 demonstrates that the average values of activation energy (E) of ROME were

276

estimated as 87.7kJ/mol, 85.5 kJ/mol and 80.9 kJ/mol calculated using FRD, FWO and MCR

277

approaches, respectively. The values are lower than that of RSO-II. HIgher molecular weight of

278

triglycerides of RSO require high energy for thermal decomposition and volatilization

279

differences in the values of activation energy (E) of RSO and ROME revealed that

280

decomposition or/and volatilization mechanism of the samples occur in different manner, and

281

RSO was chemically modified through transesterification process (Table 3). However, the

13 ACS Paragon Plus Environment

13

. The

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 35

282

activation energy range for first stage decomposition of RSO was found to be similar with

283

ROME thermal degradation with 5–12% error. Fig.7D shows the fitness of Kissinger method for

284

ROME thermal decomposition and the activation energy was found to be 92.5 kJ/mol.

285

Considering all the above methods, the average activation energy for ROME was obtained as

286

86.65 kJ/mol. The calculated average activation energy is in good agreement with for sunflower

287

oil methyl ester

288

cooking oil methyl ester (50.07 kJ/mol) and palm oil methyl ester (54.09 kJ/mol) 10. This may be

289

due to the fact that the physico-chemical-thermal properties of the parent feedstocks and

290

produced biodiesels differ.

291

13

. However, it is higher than peanut oil methyl ester (49.71 kJ/mol), waste

Order of reactions for RSO and its methyl ester. Most of the investigations assumed zero11

292

order or first-order reaction for thermal decomposition of oil and biodiesel samples

293

present study, the dependency of order of reaction (n) on temperature for thermal degradation of

294

RSO and ROME were evaluated through the Avrami theory. The regression plots of RSO (stage

295

I and II) and ROME are shown in Fig. 8(A–C) and the calculated values of order of reaction for

296

thermal degradation of samples are presented in Table 4. Based to the R2 values (Table 4), the

297

Avrami theory is suitable and well fitted to estimate the values of n for thermal degradation of

298

RSO and ROME. As the decomposition temperature increased from 267 °C to 307 °C (within

299

stage I), the reaction order of RSO first decreased from 1.59 to 1.26 and then increased to 1.42.

300

Further increasing the decomposition temperature for RSO, 402 °C to 457 °C (within stage II),

301

the value of n decreased from 1.41 to 0.38. The average order of reaction for RSO-I (1.37) was

302

found higher than that of RSO-II (0.92). The signatures of different compounds in FTIR spectra

303

of RSO-I, RSO-II and ROME also suggest that the values of overall order of the reactions may

304

differ significantly. The values of n for sunflower oil reported Santos et al. 13 vary in the range of 14 ACS Paragon Plus Environment

. In the

Page 15 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

305

0.95 to 1.82. Similarly, for overall decomposition of RSO, the n value is within the range 0.39 to

306

1.59. Varying decomposition temperature for ROME, the reaction order (n) initially increased

307

from 0.92 (227 °C) to 1.10 (257 °C) and then decreased to 0.19 (317 °C). Li et al. 10 and Santos

308

et al.

309

1.77 to 1.95 (257 °C–277 °C), 0.71 to 1.13 (257 °C–277 °C) and 0.69 to 1.89 for peanut oil,

310

palm oil, waste cooking oil and sunflower methyl esters, respectively 10, 13. Similarly, the order of

311

reaction (n) was found to be varied from 0.96 to 1.1 (237 °C–267 °C) for ROME thermal

312

degradation. Furthermore, the average value of order of reaction (n) for ROME (0.69) is in good

313

agreement with soybean (0.5)

314

esters thermal decomposition.

13

reported that the values of order reactions (n) varied from 1.6 to 1.68 (269 °C–277 °),

32

, higuereta (0.7)

32

, babassu (1.4)

33

and palm (0.4)

33

oil ethyl

315

Thermodynamic parameter calculation for RSO and its methyl esters. In addition to the

316

values of activation energy and order of reaction, important thermodynamic parameters (∆H, ∆G

317

and ∆S) for thermal decomposition of RSO and it biodiesel were calculated using Eq. 7–10. It

318

can be seen from Fig. 8D that all the calculated ∆H values are positives; thermal decomposition

319

of RSO and ROME within active degradation stages is endothermic processes. Due to high

320

molecular weight of triglycerides of vegetable oil as compared to that of their fatty acid esters,

321

∆H for vegetable oil is higher than that of biodiesel

322

calculated by FWO and MCR methods for RSO within conversion interval of 0.1 to 0.9 fractions

323

are higher than that of ROME. The values of ∆H by FWO and FRD were found within the range

324

of 69.27–79.34 kJ/mol (RSO-I), 138.39–318.01 kJ/mol (RSO-II) and 61.76–84.52 kJ/mol

325

(ROME). The average values of ∆H, considering the active thermal degradation for above two

326

methods, were determined as 146.26 kJ/mol and 78.69 kJ/mol for RSO and ROME, respectively.

327

The values of ∆H are positive which indicate endothermic reactions. High values of ∆H for RSO

12

, similar to E values. The values of ∆H

15 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

34, 35

Page 16 of 35

. Li et al.12 evaluated the average values of ∆H for

328

infer the high degree of endothermicity

329

thermal decompositions of peanut oil (118.54 kJ/mol) and its biodiesel (48.08 kJ/mol) using

330

applying different kinetics methods. Oliveira et al.

331

kJ/mol) and babassu (80.38 kJ/mol) oil biodiesel. Additionally, ∆S and ∆G at To and Tmax for

332

RSO and ROME thermal degradation were evaluated and, the obtained data are presented in

333

Table 5. The negative values of ∆S and positive values of ∆G confirm that the thermal

334

decomposition of both the samples are non-spontaneous process10,

335

since the samples were subjected to forced thermal decomposition by non-isothermal conditions.

336

The values of ∆S obtained at To and Tmax for RSO and ROME thermal degradation are negative,

337

which indicates that the activated complex has a more ordered structure than the reactants, and

338

that the reactions are slower

339

required to reduce the degree of disorder at To as compared to disorder degree at Tmax for RSO-I,

340

RSO-II and ROME (Table 5). The higher value of ∆G reveals lower favorability of a reaction10,

341

12, 25, 37, 38

342

to be 145.01–201.37 kJ/mol and 115.73–152.19 kJ/mol, respectively. Higher value of ∆G for

343

RSO indicates that larger amount of heat is required for thermal decomposition as compared its

344

biodiesel, which is similar with the average activation energy values (E).

34-36

33

reported ∆H values of for palm (90.53

12, 25, 33

. This was expected

. The absolute values of ∆S indicate that higher energy is

. The favorability order for thermal degradation process of RSO and ROME were found

345

CONCLUSIONS

346

Thermal decomposition of RSO and its methyl esters (ROME) under inert atmosphere was

347

taken place in two and one stage, respectively. The thermal stability of RSO was found greater

348

than that of its biodiesel. Rate of maximum weight loss (wmax) was increased from 12.1%/min to

349

76.7%/min (RSO) and 15.9 %/min to 76.2%/min (ROME) as the heating rate increased from 10

350

°C/min to 50 °C/min. Similarly, the active thermal degradation temperature range were also 16 ACS Paragon Plus Environment

Page 17 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

351

shifted from 223 °C–471 °C to 264 °C–488 °C for RSO and 120 °C–407 °C to 193 °C–480 °C

352

for ROME. The values of activation energy (E) and enthalpy (∆H) of RSO, calculated by several

353

methods (FRD, FWO, MCR and KM) were found greater than that of ROME. The values of

354

average order of reaction obtained using Avrami theory were found to be 1.14 and 0.69 for RSO

355

and ROME, respectively. The values of order of favorability (∆G) for thermal decomposition of

356

RSO and ROME were 145.01–201.37 kJ/mol and 115.73–152.19 kJ/mol, respectively.

357

Furthermore, the positive value of ∆G and negative value of ∆S at initial (To) and maximum

358

(Tmax) thermal degradation temperatures indicate that thermal decomposition for RSO and

359

ROME are non-spontaneous process. The absorbance peaks for alkanes, alkenes, aldehydes,

360

ketones, ethers, water, carbon dioxide and carbon monoxide were detected in the evolved

361

products. The absence of absorbance peak for water in ROME degradation shows the quality of

362

produced ester from RSO through transesterification. Thermal decomposition of RSO and

363

ROME after 500 °C (~1.19 wt% loss for RSO and 1.44 wt% loss for ROME) was mainly due to

364

the formation of ketone and CO2. From TGA-FTIR, it can be concluded that, absorbance

365

characteristics for the formation of evolved product such as alkenes, alkane ethers and aldehydes

366

were maximum within temperature range of 400 °C–500 °C and 200 °C–350 °C for RSO and

367

ROME, respectively.

368

ACKNOWLEDGMENTS

369

Authors would like to acknowledge the Central Instruments Facility (CIF) and Center for

370

Energy at Indian Institute of Technology Guwahati (IITG) for providing the characterization

371

facility to conduct the sample analyses. Authors also would like to acknowledge the Center of

372

Excellence on Sustainable Polymers at IITG for providing access to TGA-FTIR instrument.

17 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

373 374 375 376 377 378 379

REFERENCES (1) Chhetri, A. B.; Tango, M. S.; Budge, S. M.; Watts, K. C.; Islam, M. R. International Journal of Molecular Sciences 2008, 9, 169–180. (2) Atabani, A. E.; Silitonga, A. S.; Ong, H. C.; Mahlia, T. M. I.; Masjuki, H. H.; Badruddin, I. A.; Fayaz, H. Renewable and Sustainable Energy Reviews 2013, 18, 211–245. (3) Saez-Bastante, J.; Pinzi, S.; Jimenez-Romero, F. J.; de Castro, M. D. L. Energy Conversion and Management 2015, 96, 561–567.

380

(4) No, S. Y. Renewable and Sustainable Energy Reviews 2011, 15, 131–149.

381

(5) Ahmad, J.; Yusup, S.; Bokhari, A.; Kamil, R. N. M. Energy Conversion and Management

382 383 384 385 386

2014, 78, 266–275. (6) Reshad, A. S.; Barman, P.; Chaudhari, A. J.; Tiwari, P.; Kulkarni, V.; Goud, V. V.; Sahoo, N. Energy and Fuels 2015, 29, 5136–5144. (7) Atabani, A. E.; Badruddin, I. A.; Mahlia, T. M. I.; Masjuki, H. H.; Mofijur, M.; Lee, K. T.; Chong, W. T. Energy Technology 2013, 1, 685–694.

387

(8) Math, M. C. Energy for Sustainable Development 2007, 1, 100–104.

388

(9) Zhao, H.; Cao, Y.; Orndorff, W.; Cheng, Y.; Pan, W. Journal of Thermal Analysis and

389

Calorimetry 2012, 109, 1145–1150.

390

(10) Li, H.; Niu, S.; Lu, C.; Cheng, S. Energy Conversion and Management 2015, 98, 81–88.

391

(11) Volli, V.; Purkait, M. K. Fuel 2014, 117, 1010–1019.

392

(12) Li, H.; Niu, S.; Lu, C.; Wang, Y. Energy and Fuels 2015, 29, 5145–5153.

393

(13) Santos, A. G. D.; Caldeira, V. P. S.; Farias, M. F.; Araujo, A. S.; Souza, L. D.; Barros, A.

394

K. Journal of Thermal Analysis and Calorimetry 2011, 106, 747–751.

18 ACS Paragon Plus Environment

Page 18 of 35

Page 19 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

395 396 397

Energy & Fuels

(14) Ma, Z.; Chen, D.; Gu, J.; Bao, B.; Qisheng, Z. Energy Conversion and Management 2015, 89, 251–259. (15) Sauza, A. G.; Danta, H. J.; Silva, M. C. D.; Santos, I. M. G.; Fernandes, V. J.; Sinfronio, F.

398

S. M.; Teixeira, L. S. G.; Novak, C. Journal of Thermal Analysis and Calorimetry 2007, 90,

399

945–949.

400 401

(16) Reshad, A. S.; Panjiara, D.; Tiwari, P.; Goud, V. V. Journal of Cleaner Production 2017, 142, 3490–3499.

402

(17) Reshad, A. S.; Tiwari, P.; Goud, V. V. Fuel 2015, 150, 636–644.

403

(18) Friedman, H. L. Journal of Polymer Science: Part C 1964, 6, 183–195.

404

(19) Flynn, J. H.; Wall, L. A. Journal of Research of the National Bureau of Standards —A.

405

Physics and Chemistry 1966, 6, 487–523.

406

(20) Kissinger, H. E. Analytical Chemistry 1957, 29, 1702–1706.

407

(21) Ebrahim-Kahrizsangi, R.; Abbasi, M. H. Transactions of Nonferrous Metals Society of

408 409 410

China 2008, 18, 217–221. (22) Niu, S.; Liu, M.; Lu, C.; Li, H.; Huo, M. Journal of Thermal Analysis and Calorimetry

2014, 115, 73–79.

411

(23) Niu, S.; Han, K.; Lu, C.; Sun, R. Applied Energy 2010, 87, 2237–2242.

412

(24) Vayazovkin, S.; Sbirrazzuoli, N. Analytica Chimica Acta 1997, 355, 175–180.

413

(25) Kim, Y. S.; Kim, Y. S.; Kim, S. H. Environmental Science and Technology 2010, 44,

414

5313–5317.

415

(26) Sunitha, M.; Nair, R. C. P.; Krishnan, K. Thermochimica Acta 2001, 374, 159–169.

416

(27) Straszko, J.; Olszak-Humienik, M.; Mozejko, J. Thermochimica Acta 1997, 292, 145–150.

19 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

417 418

(28) Pourmortazavi, S. M.; Kohsari, I.; Teimouri, M. B.; Hajimirsadeghi, S. S. Materials Letter

2007, 61, 4670–4673.

419

(29) Gai, C.; Dong, Y.; Zhang, T. Bioresource Technology 2013, 127, 298–305.

420

(30) Ko, K.; Rawal, A.; Sahajwalla, V. Energy Conversion and Management 2014, 86, 154–

421 422 423 424 425 426 427

164. (31) Liu, Z.; Qi, N.; Luan, Y.; Sun, X. Advances in Materials Science and Engineering 2015, 2015, 1–8. (32) Raodriguez, R. P.; Sierens, R.; Verhelst, S. Journal of Thermal Analysis and Calorimetry

2009, 96, 897–901. (33) Oliveira, L. E.; Giordani, D. S.; Paiva, E. M.; De Castro, H. F.; Da silva, M. L. C. P. Journal of Thermal Analysis and Calorimetry 2013, 111, 155–160.

428

(34) Sokoto, M. A.; Singh, R.; Krishna, B. B.; Kumar, J.; Bhaskar, T. Heliyon 2016, 2.

429

(35) Aravindakshan, K. K.; Muraleedharan, K. Thermochim Acta 1990, 159, 101–107.

430

(36) Markovska, I. G.; Bogdanov, B.; Nedelchev, N. M.; Gurova, K. M.; Zagorcheva, M. H.;

431 432 433 434 435

Lyubchev, L. A. Journal of the Chinese Chemical Society 2010, 57, 411–416. (37) Mo, Y.; Zhao, L.; Chen, C.; Tan, G. Y. A. T.; Wang, J. W. Journal of Thermal Analysis and Calorimetry 2013, 111, 781–788. (38) Vecchio, S.; Cerretani, L.; Bendini, A.; Chiavaro, E. Journal of Agricultural and Food Chemistry 2009, 57, 4793–4800.

436

437 438 439 440 20 ACS Paragon Plus Environment

Page 20 of 35

Page 21 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

441

LIST OF TABLE CAPTIONS

442

Table 1: Physico-chemical properties of rubber oil (RSO) and its methyl esters (ROME) samples

443

Table 2: TG characteristic properties for active pyrolysis stages of RSO and ROME samples

444

Table 3: The value of activation energy deduced from FRD, FWO and MCR methods

445

Table 4: The value of reaction order for thermal decomposition of RSO and ROME calculated

446

using Avrami theory

447

Table 5: Thermodynamic parameters of RSO and ROME samples at TO and Tmax for heating rate

448

of 10 °C/min

449 450 451 452 453 454 455 456 457 458 459 460 461 462

21 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

463

Table 1: Physico-chemical properties of rubber oil (RSO) and its methyl esters (ROME) samples Properties

Unit –

Specific Gravity @ 24 °C

464 465

Page 22 of 35

RSO

ROMEα Biodiesel#

Diesel§

0.91

0.883

0.86–0.9*

0.846*

Kinematic Viscosity @24°C

mm2/s

30

5.82

1.9–6*

1.9–4.1*

Kinematic Viscosity @40°C

mm2/s

13.13

3.81

1.9–6

1.9–4.1

Calorific value

MJ/kg

39.34

39.53

35min

45.62–46.48

Iodine Value

g I2/100g

113

114

120max

N/A

Acid Value

mg KOH/g oil

24

0.4

0.5max

0.35

Saponification value

mg KOH/g oil

235.28

190

N/S

N/A

Refractive index @ 24°C



1.47

1.45

N/S

N/S

Cloud point

°C

3

2.5

(-3) –12

(-15) –5

Pour point

°C

2**

-3 **

(-15) –10

-20

Cetane No.



44

49.9

47min

40min

Moisture content

%

0.27

0.33

0.05max

0.05max

Flash point

o

273

131

130min

52 – 96

Fire point

o

282

146

N/S

N/S

C C

Rubber seed oil methyl ester, *measured at 40 oC, #Standards ASTM test methods, §Fossil fuel, ** DSC method, N/S: not specified, N/A: not applicable α

466 467 468 469 470 471 472 22 ACS Paragon Plus Environment

Page 23 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

473 474

Table 2: TG characteristic properties for active pyrolysis stages of RSO and ROME samples Sample

Properties

Heating rate (°C/min) 10

20

30

40

50

Stage

I

II

I

II

I

II

I

II

I

II

To, °C

223

353

242

373

249

382

253

384

264

395

Tf, °C

309

471

334

480.0

345.6

483

350

487

365

488

∆wt, %

20

69.5

20

69.5

20

69.5

20

69.5

20

69.5

Tmax, °C

271

416

296

424

310

438

318

443

328

446

wmax,%/min 3.26

12.1

5.9

35.2

8.8

44.6

11.7

60

13.6

76.7

RSO

To, °C

120



165



172



181



193



Tf, °C

407



456



470



474



480



∆wt., %

97.2



97.2



97.2



97.2



97.2



Tmax, °C

250



268.5



277.5



283.5



287.9



wmax,%/min 15.9



29.1



41.4



55.2



76.2



ROME

475

To: Initial temperature for the main mass loss (°C), Tf: Final temperature for the main mass loss

476

(°C), Tmax: Temperature for maximum rate of mass loss (°C), ∆wt: change of mass loss and wmax:

477

maximum mass loss rate (%/min)

23 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

478

Page 24 of 35

Table 3: The value of activation energy deduced from FRD, FWO and MCR methods Activation Energy (kJ/mol) at Conversion (%)* Average Sample

Stage

I RSO

II

ROME

479 480

I

Method #

10

20

30

40

50

60

70

80

90

FRD

72.8

70.7

69.8

73.5

74.1

75.2

77.72

83.7

92.7

76.7

FWO

83.8

81.6

80.2

79.6

79.4

79.5

79.9

81.1

83.7

80.9

MCR

79.3

76.8

75.2

74.4

74.1

74.1

74.3

75.5

78.1

75.8

FRD

164.2

196.9

220.7

212.1

231.2

261.1

283.1

328.6

433.3

258.9

FWO

147.8

162.8

177.9

188.9

197.2

209.4

225.4

250.9

319.2

208.8

MCR

144

160

176.3

187.3

196

209.5

225.3

252.1

324.1

208.3

FRD

71.8

89.6

93.9

93.3

90.9

87.1

83.9

78.2

100.9

87.7

FWO

70.5

78.4

84.6

88.1

89.9

90.4

90.1

88.1

89.7

85.5

MCR

65.9

74.0

80.3

83.8

85.6

86.1

85.5

83.3

84.2

80.9

*

Stage by stage conversion (α) of active pyrolysis of the samples, # average R2 was found to be ~0.994 and overall average activation energy (E) for RSO decomposition was found to be 167.8 kJ/mol (FRD), 144.9 kJ/mol (FWO) and 152.5 kJ/mol (MCR)

481 482 483 484 24 ACS Paragon Plus Environment

Page 25 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

485 486

Energy & Fuels

Table 4: The value of reaction order for thermal decomposition of RSO and ROME calculated using Avrami theory RSO-I

RSO-II

ROME

Temp (°C)

n

R2

Temp (°C)

n

R2

Temp (°C)

n

R2

267

1.59

0.99

402

1.41

0.99

227

0.92

0.99

272

1.41

0.99

407

1.24

0.99

237

0.96

0.99

282

1.31

0.99

417

1.16

0.99

247

1.04

0.99

292

1.26

0.99

427

1.06

0.99

257

1.10

0.99

297

1.26

0.98

432

0.88

0.99

267

1.05

0.99

307

1.42

0.98

442

0.69

0.98

277

0.86

0.99

452

0.50

0.98

287

0.61

0.98

457

0.38

0.99

297

0.39

0.98

307

0.24

0.98

317

0.19

0.99

327

0.19

0.99

Average (n)

1.37

0.92

487 488 489 490 491 492 493 494

25 ACS Paragon Plus Environment

0.69

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

495 496

Page 26 of 35

Table 5: Thermodynamic parameters of RSO and ROME samples at To and Tmax for heating rate of 10 °C/min Sample

Temp (°C)

Log A (s-1)

∆H (kJ/mol)

∆G (kJ/mol)

∆S (kJ/mol K)

=223

4.89

+63.75

+145.01

-163.82

Tmax =271

5.67

+78.02

+159.41

-149.62

To

=353

7.46

+122.33

+195.35

-116.64

Tmax =416

11.83

+178.11

+201.37

-33.76

To

=120

1.01

+22.89

+115.73

-236.22

Tmax =250

7.2

+89.44

+152.19

-119.97

To RSO-I

RSO-II

ROME 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523

26 ACS Paragon Plus Environment

Page 27 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

524

LIST OF FIGURE CAPTIONS

525

Fig. 1: TGA and DTG profiles (A, B) for RSO and (C, D) for ROME samples

526

Fig. 2: FT-IR spectra for evolved products during thermal decomposition of (A) RSO and (B)

527

ROME at heating rate of 40 °C/min

528

Fig. 3: Formation of evolved products during thermal decomposition of (A) RSO and (B) ROME

529

at heating rate of 40 °C/min

530

Fig. 4: TGA mass conversion for (A) RSO-I, (B) RSO-II and (C) ROME samples

531

Fig. 5: Fig.5: Regression plots based on (A) FRD for RSO-I, (B) FRD for RSO-II, (C) FWO for

532

RSO-I and (D) FWO for RSO-II

533

Fig. 6: Regression plots based on (A) MCR for RSO-I, (B) MCR for RSO-II and (C) KM for

534

RSO-I and RSO-II

535

Fig. 7: Regression plots based on (A) FRD, (B) FWO, (C) MCR and (D) KM for ROME thermal

536

degradation

537

Fig. 8: Regression plots to calculate order of reaction proposed by Avrami theory for (A) RSO-I,

538

(B) RSO-II and (C) ROME; (D) Enthalpy change for RSO and ROME based on FRD and MCR

539

methods

540 541 542 543 544 27 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

545

546 547

Fig.1: TGA and DTG profiles (A, B) for RSO and (C, D) for ROME samples

548 549 550 551 552 553 554 555

28 ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

556 557 558

559 560 561

Fig. 2: FT-IR spectra for evolved products during thermal decomposition of (A) RSO and (B) ROME at heating rate of 40 °C/min

562 563 564 565 566 29 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

567 568 569 570 571 572

573 574

Fig.3: Formation of evolved products during thermal decomposition of (A) RSO and (B) ROME

575

at heating rate of 40 °C/min

576 577 578 579 580 581

30 ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

582 583 584 585 586 587

588 589

Fig.4: TGA mass conversion for (A) RSO-I, (B) RSO-II and (C) ROME samples

31 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

590 591

Fig.5: Regression plots based on (A) FRD for RSO-I, (B) FRD for RSO-II, (C) FWO for RSO-I and (D) FWO for RSO-II

32 ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Energy & Fuels

592

593 594

Fig.6: Regression plots based on (A) MCR for RSO-I, (B) MCR for RSO-II and (C) KM for RSO-I and RSO-II 33 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

595

596 597

Fig.7: Regression plots based on (A) FRD, (B) FWO, (C) MCR and (D) KM for ROME thermal degradation 34 ACS Paragon Plus Environment

Page 34 of 35

Page 35 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

598 599 600 601

602 603

Fig.8: Regression plots to calculate order of reaction proposed by Avrami theory for (A) RSO-I,

604

(B) RSO-II and (C) ROME; (D) Enthalpy change for RSO and ROME based on FWO and MCR

605

methods

606

35 ACS Paragon Plus Environment