**Open Access**-
**Total Downloads**: 73 -
**Authors :**Olukunle E. Itabiyi , Pious O. Okekunle , Emmanuel O. Sangotayo , David O. Ogunwobi, Kehinde M. Adeleke, Olutosin O. Ilori -
**Paper ID :**IJERTV8IS080147 -
**Volume & Issue :**Volume 08, Issue 08 (August 2019) -
**Published (First Online):**26-08-2019 -
**ISSN (Online) :**2278-0181 -
**Publisher Name :**IJERT -
**License:**This work is licensed under a Creative Commons Attribution 4.0 International License

#### Effects of Pyrolysis Variables on the Product Yields During Pyrolysis of Palm Kernel Shells

1Olukunle E. Itabiyi*, 2Pious O. Okekunle, 3Emmanuel O. Sangotayo, 4David O. Ogunwobi 1,2,3,4 Department of Mechanical Engineering, Ladoke Akintola University of Technology,

P.M.B 4000, Ogbomoso Nigeria.

5Kehinde M. Adeleke, 6Olutosin O. Ilori 5,6 Department of Mechanical Engineering, Adeleke University,

Ede, Nigeria.

#### Abstract:- The effect of operating variables on the pyrolysis products when Palm Kernel Shell (PKS) was subjected to pyrolysis process was investigated in this work. 0.5 kg of dried sample of PKS was loaded into a steel retort, and the retort interior was rendered airtight. The retort was then placed into the furnace chamber and the PKS was pyrolysed at 300 oC between 10 – 30 minutes at 5 minutes interval. This was repeated for temperatures of 400, 500, 600 and 700 oC and in each case, the quantities of char, tar and pyro – gas produced were determined. Response surface methodology (RSM) was used to develop polynomial regression model and investigate the effect of changes in the level of pyrolysing temperature and duration on the product yields using Full Factorial Design (FFD). The contribution of pyrolysis temperature, duration and their squares (A, B, A2 and B2) to the model developed are significant. It was observed that the experimental data fitted better because of the Pred R-Squared of 0.9242 is in reasonable agreement with the Adj R-Squared of 0.9383. The agreement between the predicted and experimental values describe the accuracy of the model developed and can be used to navigate within the design space. The product yields minimum value of 0.89% was achieved at pyrolysis temperature and duration values of 300 oC and 10 min respectively. The optimum conversion product yields for dried PKS char, tar and gas at their respective pyrolysis conditions were 99 wt% char at 304 oC and 12 min., 35 wt% tar at 700 oC and 27 min., and 40 wt% gas at 700 oC and 30 min. The results obtained showed that PKS can be readily pyrolised to obtain optimum yield of bio fuels (gas, tar and char).

Keywords: Pyrolysis, Palm Kernel Shells, Bio fuels, Gas, Tar, Char and Response surface methodology.

1.0 INTRODUCTION

Periodic escalation in petroleum products prices and decreasing stock of crude oil deposits had brought non- conventional and renewable energy sources into greater focus (Itabiyi, 2017). Along with the solar and the wind energy, the long neglected but potentially rich biomass residues became the focus of intensive utilization for energy generation. Bio- fuel produced from biomass residues is an alternative to petroleum products and has received great attention during the last decades due to the environmental problems associated with the usage of fossil fuel. It already provides approximately 14% of the total worldwide energy needs (Funda, et al., 2006). Aside from its abundant availability, biomass residues have negligible contents of sulfur, nitrogen and ash which give lower emissions of SO2, NOx, soot and net emission of CO2 compared to conventional fossil fuels, thus keeping the environment and the public's health safe (Qi, et al., 2006; Tsai, et al., 2006a, Okekunle, et al., 2018). In

countries of excess production of agricultural residues such as Nigeria, eco-friendly system and more energy efficient utilization alternatives, that meets the needs of the present without compromising the ability of future generations must be developed (Itabiyi, et al., 2018).

One of the technologies by which biomass could be converted to energy is through pyrolysis. Among the thermo- chemical conversion processes, pyrolysis is the viable process for biomass upgrading by cracking polymeric structure of lignocellulosic materials and converting them into fraction consisting of a solid carbon- rich residue in form of charcoal and a range of volatile products, which include condensable and non-condensable gases. Volatile fraction can be used as a fuel or as a chemical feedstock. The remaining solid fraction can find several applications, such as in the production of activated carbon or used directly as a solid fuel. Therefore, pyrolysis is an important process that can contribute much more, to solving energy problem in the world especially developing countries (Lucas and Itabiyi, 2012). Distribution of pyrolysis products depends on such operating conditions as type of feedstock, reaction time, pyrolysis temperature and sweep gas flow rate (Choi et al., 2010; Rodjeen et al., 2006). To obtain high liquid yield, the pyrolysis conditions require high heating rate, moderate temperature (450-550 oC) and short vapor residence time (Mohan et al., 2006). Higher proportion of gas product is obtained from applying high temperature, low heating rate and long residence time, while slow heating at lower temperatures and long vapor residence time favor the formation of char product (Mohan et al., 2006; Nugranad, 1997; Onay and KoÃ§kar, 2003). Also, cellulosic composition of biomass and process variables including pyrolysis temperature and heating rate can have a profound influence on the chemical compositions of both gas and liquid products (Yaman, 2004; Branca et al., 2003; Yang et al., 2007).

Weerachanchai, et al, 2011 studied the effects of pyrolysis temperature on the product yields of palm kernel cake and cassava pulp residue and noted that the char yield decreased sharply from 300 to 500 oC followed by a slow decrease at higher temperatures and approaching a constant value at 800 oC. Temperature is also found to have an influence on the composition of element in the products. For example, studies by few researchers (Zanzi, et. al. 2002) showed that carbon content increases with an increase in temperature, while the hydrogen and oxygen contents decrease.

Higher temperature, smaller particle size, and increased heating rate resulted in decreased char yield from pyrolysis of agricultural residues (Itabiyi and Lucas, 2013). The cracking of the hydrocarbons with an increase in the hydrogen content was favoured by a higher temperature and by using smaller particles (Zanzi et al., 2002). Wood gives more volatiles and less char than straw and olive waste. The effect of particle size in the range of 0.71 3.56 mm on product yields at 7000C was studied by Weerachanchai, et al. (2011) using palm kernel and cassava pulp residues. Their results showed that the pyrolysis of the two biomasses with the average of

2.03 mm gave the maximum in the liquid yield of 54.3 wt% and 42.4 wt% for palm kernel cake and cassava pulp residue, respectively. Their study further revealed that for particle sizes smaller than 2.03 mm, higher gas and char yields but a lowering in liquid yield were obtained as the particle size was decreased. Their study suggested that it is likely that the smaller size of biomass particles could affect greater heat transfer because of less temperature inside the particle, thus giving higher yields of released gases and volatiles. This work investigated the influence of pyrolysis temperature and time on the product yields during the pyrolysis of Palm kernel shells.

METHODOLOGY

Palm kernel shells preparation for the experiment

Palm kernel shells used for the pyrolysis experiments in this study was procured from an oil-palm industry in Iresa apa, Oyo State, South-Western Nigeria. The residues were cleaned in order to remove foreign particles such as stones, leaves, debris and other unwanted components. The weight of the sample (W1) was measured using Ohaus top loading digital weighing scale of sensitivity Â± 0.001 g (Model: PA4102,range: 0-4100 g, Ohaus company, Manufactured in Switzerland) and then oven-dried at a temperature of 105 oC until constant weight (W2) was obtained in accordance with official methods of the ASTM D5373-02 (2005).

Methods.

Pyrolysis experiments were carried out to determine the effect of operating variables on the product yields from PKS.

0.5 kg of dried PKS were fed into the retort. The retort was placed into the furnace and pyrolysed at around 300, 400, 500, 600 and 700 oC. The retort was connected through a pipe to the condensate receiver which was placed in an ice-cooling unit for the quick recovery of the condensable products (tar), and from the condensate receiver the uncondensed gases moved through a rubber hose into the gas collection unit.

The char in the retort and the tar in the condensate receiver were collected and weighed using Ohaus top loading digital weighing balance. The weight of gas was evaluated by subtraction. The percentage of product yields was calculated from equation 1.

Experimental Design

Full-Factorial Design (FFD) of response surface methodology was used for the experimental design to optimise the pyrolysis product yields from PKS. FFD consisted of a two- factor, three-level design comprising the pyrolysis temperature and pyrolysis duration of the feedstock as the independent variables while pyrolysis product yields consisting of char, tar and gas as the dependent variables or the responses were used as shown in Table 1. A centre point for the design was selected with factors at a level of medium standards as shown in Table 2. With the centre point design selected, the actual values of each factor were calculated. The design was based upon the symmetrical selection of variation about the centre point and levels of variations were chosen to be within the boundary range of the variables. The coded and actual values of the variables at various levels and responses are given in the Table 2. Three replications were carried out for all experimental design conditions and the average recorded. Thirteen experimental runs were carried out and the order of the experiment was fully randomised to reduce the effect of the unexplained variability in the observed responses due to extraneous factor as recommended by Singh et al (2003).

Table1: Experimental Factors and Responses

Type

Variables

Symbols

Factors

Responses

Temperature Duration

Char yield Tar

yield Gas yield

A B

Yc

Yt

Yg

Table 2: Experimental Values of Coded Levels

Coded Levels

Factors

-1 0 +1

A (0C) 300 500 700

B (Min) 10 20 30

Analysis of Data and Response Equations.

Regression Models were developed for PKS product yield and each of the product yields as a function of the two factors. The Design Expert 6.0.8 software was used to analyse the data obtained from the pyrolysis of PKS for developing response equations, Analysis of Variance (ANOVA), to generate surface plots and determine optimum pyrolysis conditions and product yield using its optimization toolbox. In multiple regressions, as in the present case, R2, which is the square of the adjusted coefficient of determination and standard error are the indices. F statistics shows the significance of overall model while the t-statistics tests shows the significance of each of the variables of the model. The Functions was assumed to be approximated by a second degree polynomial equation as shown in equation 2.

Percentage product yields Y mass of

mass of

product 100

sample

(1)

Y b0

m

m

bi xi

i1

m

m

bii

i1

m

x i b

x i b

2

ij

i j

xi x j

(2)

where Y is the predicted response, b0 is the value of the fitted response at the centre point, and bi, bii, bij are linear, quadratic

Y

Y

CPKS

64.77 22.07 A 17.06B

and cross product regression terms respectively. m is the

number of factors considered in the study which is equal to 2.

R2 = 0.8709 (4)

Optimization of the Product Yields.

A nonlinear programming problem of the form of equation 2 was formed from the vector of equation 2 as shown in equation 3. The optimization problem statement to maximize the product yields was formulated as shown in equation 3.

Y

t

t

PKS

g

g

Y

PKS

25.83 9.77 A 9.13B 7.52B2

R2 = 0.8711 (5)

15.2311.35A8.09B 5.14AB

R2 = 0.9383 (6)

Maximize Y = Subject to

f AB

(3)

where:

C

C

Y

PKS

t

t

Y

PKS

= Yield of char from PKS (wt%)

= tar yield from PKS (wt%)

LA

A U A

Y

Y

g PKS

= gas yield from PKS (wt%)

LB B UB

Where Y is the product yields, Li is the lower limit of the factors and Ui is the upper boundary of the factors. The line search problem stated in equation 2 was embedded and solved in the optimization routine of design expert 6.0.8 version to obtain the optimal yields and the corresponding optimal process variables.

RESULTS AND DISCUSSIONS

Based on t-test, the regression coefficient that are not significant at 95% confidence level were discarded while only those ones that are significant were used to develop the final model.

Response Equations for PKS Product Yields.

The effect of FFD on the PKS pyrolysis product yields (char, tar and gas yields) is as shown on Table 3 that was subsequently used to fit the response equations for product yields. Multiple regression analysis was used as tools of assessment of the effects of two or more independent factors on the dependent variables (Boomee et al, 2010). The coefficients of determination (R2) is a measure of the total variation of the observe values of the product yields about the mean explained by the fitted model (Shridhar et al, 2010). The factors of the models, their parameters estimates and the statistics of the estimates for the best functions adopted, taking into consideration all main effects, linear, quadratic, and interaction for each model are as shown on Table 4. The coefficients of determination (R2) for the responses (char, tar and gas) were 0.8709, 0.8711 and 0.9383 respectively. The coefficient of determination (R2) were high for response surfaces, and indicated that the fitted quadratic models accounted for more than 89% of the variance in the experimental data. Base on the p values, the regression coefficient that were significant at p<95% were selected for the models that resulted in equations 4 – 6. Analyses of variance (ANOVA) were conducted to evaluate the adequacy and consistency of the models using F-statistic. The analysis of variance of the models is presented in Table 5. The results presented on Table 5 showed the F- values for char, tar and gas as 41.49, 17.22 and 61.86 respectively. These values were significant at p< 0.05 indicating good model fit.

A = Temperature (oC) B = Time (Minutes).

Optimization of Pyrolysis Process.

Response surface methodology was used for the optimization of the pyrolysis process of the PKS and for understanding the factors affecting the pyrolysis process. The models were useful for indicating the direction in which to change the variable in order to maximise the yields of char, tar and gas. The multiple regression equations were solved using Design Expert 6.0.8. The regression equation was optimized for maximum value, to obtain the optimum conditions. The optimum actual values obtained for PKS pyrolysis roduct yields and their respective pyrolysis conditions are: 99.27% char at A = 304.38 oC and B = 12.43 minutes, 35.22% tar at A = 700.00 oC and B = 26.94 minutes and 39.81% gas at A = 700 oC and B = 30 minutes.

The linear effects of temperature and time are the primary determining factors of the responses as shown in Table 4. Pyrolysing temperature as a single factor was the most influential factor, because of its higher F-value. The temperature at which pyrolysis process was conducted is highly significant (p<0.05) with an F-value of 112.85 as shown in Table 4.

Figures 1-3 show three-dimensional (3D) surface plots and accompany contour plot for the relationship between the independent and dependent variables for chosen model. The cubic response surface plot shown in Figure 1(a) depicts the effect of the pyrolysing temperature and time on the PKS char yield. From the contour plot in Fig 1(b) it is observed that the surface area decreases as the pyrolysing temperature and time increase. Fig 1(b) shows that, char yield of PKS decreases as the pyrolysing temperature and time increase. Itabiyi and Lucas (2016) ; Okekunle et al., (2016a) ; Okekunle et al., (2016b) and Adeleke et al., (2018) observed similar trend when they conducted pyrolysis experiment on oil palm trunk, cassava chaff and cassava peel respectively in a fixed bed pyrolysis reactor. Mohamad (2008) reported that, the decrease in char yield with an increase in pyrolysing temperature could either be due to secondary decomposition

of the char or through the greater primary decomposition of the PKS at higher temperatures.

The cubic response surface plot shown in Figure 2(a) depicts the effect of pyrolysing temperature and time on the PKS tar yield. It was observe from the contour plot in Figure 2(b) that the surface area increases as the pyrolysing temperature and time increases. Figure 2(a) cubic response surface indicates that the tar yield increases as the pyrolysing temperature and time increase to optimum condition while further increase in pyrolysing temperature and time led to decrease in tar yield. This shows that, there was a mutual interaction between the pyrolysing temperature and time on tar yield. Itabiyi and Lucas (2016) ; Okekunle et al., (2016a) ; Okekunle et al., (2016b) and Adeleke et al., (2018) observed similar trend when they conducted pyrolysis experiment on oil palm trunk,

cassava chaff and cassava peel respectively in a fixed bed pyrolysis reactor. Pyrolysis process at higher temperature might have led to more tar cracking resulting into higher gas yield and lower tar yield.

The cubic response surface plot shown in figure 3(a) depicts the effect of pyrolysing temperature and time on the PKS gas yield. It is observed from the contour plot in figure 3(b) that the surface area increases as the pyrolysing temperature and time increased. Figure 3(a) the cubic response surface indicates that the gas yield increases as the pyrolysing temperature and time increase. The increase in gaseous products as the reaction temperature increases might be due to the secondary cracking of the pyrolysis vapours at higher temperatures, or secondary decomposition of the char at the higher temperatures (Mohamad, 2008).

Table 3: Full Factorial Design Arrangement and Responses for PKS

Coded Level

Actual Values

Responses

Exp. No.

A(oC)

B(min)

Temp. ( oC)

Time (min)

YC

PKS

YT

PK S

YG

PKS

1

1

0

700

20

34.98

36.83

28.2

2

0

0

500

20

61.1

25.77

13.12

3

-1

-1

300

10

98.13

0.89

0.98

4

0

0

500

20

61.1

25.77

13.12

5

0

0

500

20

61.1

25.77

13.12

6

0

-1

500

10

82.13

5.89

11.98

7

0

1

500

30

43.16

31

26.84

8

-1

0

300

20

94.38

7.12

4.2

9

0

0

500

20

61.1

25.77

13.12

10

0

0

500

20

61.1

25.77

13.12

11

1

1

700

30

29.92

30.19

39.89

12

-1

1

300

30

79.29

13.15

7.51

13

1

-1

700

10

74.46

12.77

12.77

A = Temperature (oC) B = Time (min)

Y = Yield of tar from PKS (wt%)

t

t

PKS

Y

Y

CPKS

= Yield of char from PKS (wt%)

Y = Yield of gas from PKS (wt%)

g

g

PKS

Table 4: Parameter Estimation from Regression Analysis of PKS

Estimated Coefficient of the fitted model for properties based on t-statistics

Responses

Model Factors

Coefficients

F-Values

p-Values

Model

64.77

41.49

0.0001*

Yield of Char YC

PKS

A

-22.07

51.95

0.0001*

B

-17.06

31.03

0.0002*

R2

0.8709

Model

25.83

17.22

0.0008*

Yield of Tar Yt

PKS

A

9.77

35.56

0.0006*

B

9.13

31.05

0.0008*

A2

-3.99

2.73

0.1427

B2

-7.52

9.69

0.0170*

AB

1.29

0.41

0.5409

R2

0.811

Model

15.23

61.86

0.0001*

Yield of Gas Yg

PKS

A

11.35

112.85

0.0001*

B

8.09

57.35

0.0001*

AB

5.14

15.39

0.0035*

R2

0.9383

* Significant at p <0.05 level

MODEL EQUATION OF OPH PRODUCT YIELDS

C

C

Y

OPH

t

t

Y

OPH

g

g

Y

OPH

= 38.78 1698A 25.32B + 6.29A2 + 14.13B2 R2 = 0.9690

= 47.42 + 7.53A + 13.63B2 17.08A2 13.12B2 R2 = 0.8934

= 13.74 + 9.45A + 11.69B + 10.83A2 R2 = 0.8937

Responses

Source of

Degree of

Sum of

Mean

Adjusted

Variance

Freedom

Squares

Square

F

R2

Regression

5

6572.58

1314.66

76.02

0.969

Residual

7

121.05

17.29

Total

12

6693.63

Lack of fit

3

121.05

40.35

Regression

5

3512.49

702.5

21.1

0.8934

Residual

7

233

33.29

Total

12

3745.49

Lack of fit

3

233

77.67

Regression

5

1736.29

347.26

21.19

0.8937

Residual

7

114.74

16.39

Total

12

1851.03

Lack of fit

3

114.74

38.25

Responses

Source of

Degree of

Sum of

Mean

Adjusted

Variance

Freedom

Squares

Square

F

R2

Regression

5

6572.58

1314.66

76.02

0.969

Residual

7

121.05

17.29

Total

12

6693.63

Lack of fit

3

121.05

40.35

Regression

5

3512.49

702.5

21.1

0.8934

Residual

7

233

33.29

Total

12

3745.49

Lack of fit

3

233

77.67

Regression

5

1736.29

347.26

21.19

0.8937

Residual

7

114.74

16.39

Total

12

1851.03

Lack of fit

3

114.74

38.25

Table 5: Analysis of Variance (ANOVA) for the Responses

C

C

Y

OPH

t

t

Y

OPH

g

g

Y

OPH

Yield of Char

30.00

100.399

79.2394

25.00

58.0795

29.8662

Yield of Char

Yield of Char

36.9195

15.7596

B: Time

B: Time

20.00

43.9728 5

30.00

25.00

20.00

500.00

600.00

700.00

15.00

10.00

86.2927

72.1861

58.0795

B: Time

15.00

10.00

300.00

400.00

A: Temperature

300.00 400.00 500.00 600.00 700.00

A: T emperature

(b)

Figure 1: (a) Response Surface Cubic Plot showing the 3D Effects of Temperature, Time and their Interaction on the Optimum Char yield from

PKS. (b) Contour Plot of Figure 1a.

51.6394

30.00

Yield

37.3962

25.00

23.1531

8.90991

-5.33325

20.00

5

51.6394

30.00

Yield

37.3962

25.00

23.1531

8.90991

-5.33325

20.00

5

of Tar

42.1439

42.1439

32.6485

32.6485

Yield of Tar

Yield of Tar

B: Time

B: Time

15.00

30.00

25.00

20.00

500.00

600.00

700.00

10.00

13.6576

4.16219

23.1531

B: Time

15.00

10.00

300.00

(a)

400.00

A: Temperature

300.00 400.00 500.00 600.00 700.00

A: T emperature

Figure 2 (a) Response Surface Cubic Plot showing the Effects of Temperature, Time and their interaction on the Optimum Tar

yield from PKS. (b) Contour Plot of Figure 2a.

Yield of Gas

30.00

39.3463

47.2242

35.4073

25.00

31.4683

23.5904

23.5904

B: Time

B: Time

23.5904

15.7125

15.7125

Yield of Gas

Yield of Gas

11.7735

-0.0434328

20.00 5

7.83451

7.83451

15.00

30.00

25.00

20.00

500.00

600.00

700.00

10.00

300.00 400.00 500.00 600.00 700.00

B: Time

15.00

10.00

300.00

400.00

A: Temperature

A: T emperature

(b)

Figure 3: (a) Response Surface Cubic Plot showing the Effects of Temperature, Time and their Interaction on the Optimum Gas yield from PKS.

Contour Plot of Figure 3a.

4.0 CONCLUSION

This study has shown that PKS lend itself to pyrolysis process under different pyrolysis conditions and these have produced different amounts of pyrolytic products (Char, tar and gas). In general, as the pyrolysis temperature increases, the char or solid production decreases and vice versa. The study has also demonstrated the applicability of response surface methodology in selecting pyrolysis variables that maximises product yields from PKS. Pyrolysis of PKS gave the optimum char yield of 99.27 w% at 304.28 oC, optimum tar yield of 35.22 wt% at 700 oC and optimum gas yield of 39.81wt% at 700 oC.

5.0 REFERENCES

Adeleke K. M., Itabiyi, O. E., Ilori O. O. (2018). Temperature Effect on the Products Yield from Pyrolysis of Cassava Peel. International Journal of Scientific and Engineering Research, Vol. 9 (4) Pp 953- 958.

ASTM, Standard D5373-02, (2005). Standard Test Methods for Instrumental Determination of Carbon, Hydrogen and Nitrogen in Laboratory Samples of Coal and Coke, ASTM International, Philadelphia.

Boonmee, K., Chuntranuluck, S., Punsuvon, V. and Silayoi P. (2010). Optimization of biodiesel Production from Jatropha (Jatropha Curcas L.), using surface response methodology, Kasetsart J. Nat Sci. 44:290-299.

Branca, C. Giudicianni, P. and Blasi, C.D. (2003). Ind. Eng. Chem. Res., Vol. 42, Pp 3190.

Choi, H. S., Choi, Y. S., and Park, H.C. (2010). Korean J. Chem. Eng., Vol. 27, Pp 1164.

Funda A., Ayse E., PÃ¼tÃ¼n and Ersan P. (2006). Pyrolysis of two different biomass sample in a fixed-bed reactor combined with two different catalysts, Fuel. Vol. 85 Pp 1851-1859.

Itabiyi, O. E., Okekunle, P. O. and Ogunwobi, D. O. (2018). Influence of Operating Parameters on the Product Yields During Pyrolysis of Oil Palm Husk, Journal of Energy Technologies and Policy. Vol. 8, No. 1, Pp 6-13.

Itabiyi, O. E. (2017). Pyrolysis of Oil-palm Husk and Optimization of Pyrolytic Oil Yield, Humboldt-Kolleg International Conference on Harvesting Research outcomes: a Practical Plan to Confirm Achievement of the Millennium Development Goals, Ogbomoso Nigeria. Pp 199-205.

Itabiyi, O. E. and Lucas E. B. (2016). Determining the Optimal Product Yields from the Pyrolytic of Oil Palm Trunck, Journal of Energy Technologies and Policy. Vol. 6, No. 6, Pp 19-25.

Itabiyi, O. E. and Lucas E. B. (2013). Development of Pyrolysis Unit for the Conversion of Biomass into Fuels and Chemial Products, LAUTECH Journal of Engineering and Technology. Vol. 7, No. 2, Pp 136-142.

Lucas, E. B. and Itabiyi, O. E., (2012). Pyrolysis of Palm Kernel Shells to Fuels and Chemicals Feedstocks, LAUTECH Journal of Engineering and Technology. Vol. 7, No. 1, Pp 124-128.

Mohamad, A.S. (2008). Pyrolysis of Empty Oil Palm Fruit Bunches Using the Quartz Fluidised Fixed Bed Reactor, Available online at: http://dspace.fsktm.um.edu.mybitstream 18125821 M.Sc%20 Thesis-mohamad%20 Azri%20 Sukiran-SGR%20060006.

Mohan, D. Pittman, C. U. and Steele, P. H. (2006). Pyrolysis of Wood/Biomass for Bio-oil: A Critical Review, Journal of Energy and Fuel, Vol. 20 (3) Pp 848-889.

Nugranad, N. (1997). Pyrolysis of biomass, An Unpublished Ph.D. Dissertation, University of Leeds, Leeds.

Okekunle, P. O., Itabiyi, O. E., Bello M., Adeleke A.G., Olayanju A. and Olapade O. (2018). Biofuel Yields from Sawdust Pyrolysis of different Woods in a Fixed Bed Reactor, International Conference of Mechanical Engineering, Energy Technology and Management, IMEETMCON 2018 September 4-7, 2018, International Conference Centre, university of Ibadan, Nigeria. Pp 217-223.

Okekunle, P.O., Itabiyi, O.E., Adetola, S.O., Alayande, I.O., Ogundiran, H.O. and Odeh, K.G. (2016a). Biofuel production by pyrolysis of cassava peel in a fixed bed reactor, International Journal of Energy for a Clean Environment, Vol. 17 (1) Pp 57-65.

Okekunle, P.O., Adetola, S.O., Itabiyi, O.E., Alayande, I.O., Ogundiran, H.O., and Odeh, K.G. (2016b). Experimental investigation of the effect of microwave drying and reactor temperature on product yields from pyrolysis of cassava chaff, Nigerian Journal of Technology, Vol. 34 (4) Pp 874-879.

Onay, O. and KoÃ§kar, O. M. (2003). Energy, Vol. 28, Pp 2417.

Rodjeen, S., Mekasut, L., Kuchontara, P. and Piumsomboon, P. (2006), Korean J. Chem. Eng., Vol. 32 (2) Pp 216.

Qi Z., Chang J., Wang T. and Xu Y. (2007). Review of Biomass Pyrolysis Oil Properties and Upgrading Research, Energy Conversion and Management. Vol. 48 Pp 87-92.

Shridhar, B.S. Beena, K.V. Anita, M.V. and Paramjeet K.B. (2010). Optimization and Characterization of castor seed oil, Leonardo Journal of Science 17:59-70.

Tsai, W.T., Lee, M.K. and Chang, Y.M., (2006a). Fast pyrolysis of rice straw, sugarcane bagasse and coconut shell in an induction- heating reactor, J. Anal. Appl. Pyrolysis Vol. 76 Pp 230237.

Weerechanchai, P., Tanggsathitkulchai, C., and Tanggsathitkulchai, M., (2011). Characterization of products of pyrolysis of palm kernel cake and cassava pulp residue, Korean Journal of Chemical Engineering, Vol. 28 (2), Pp. 2262 2272.

Yaman, S., (2004). Pyrolysis of Biomass to produce Fuels and Chemical Feed Stocks, Energy Conversion and Management, Vol. 45 Pp. 651 671.

Yang, H., Yan, R., Chen, H. Lee, D. H. and Zheng, C. (2007). Fuel, Vol. 86, Pp 1781.

Zanzi, R., Sjotrom, K. and Bjornbom E. (2002). Rapid Pyrolysis of Agricultural Residues at High Temperature, Biomass and Bioenergy, Vol. 23 (5), Pp. 357 366.