Statistical approach was employed to optimize biodiesel production from sesame oil in this work. Precisely, response surface methodology was applied, and the effects of four variables, viz. reaction temperature, catalyst amount, reaction time, and methanol/oil molar ratio, and their reciprocal interactions were determined. Central composite rotatable design was used to generate a total of 30 individual experiments, which were designed to study the effects of these variables during alkali-catalyzed methanolysis of sesame oil. A statistical model predicted the highest conversion yield of sesame biodiesel to be 99.71% at the following optimized variable conditions: reaction temperature of 63°C, catalyst amount of 1.04 wt.%, and methanol/oil molar ratio of 6.24, with a reaction time of 51.09 min. Using these variables under experimental condition in four independent replicates, an actual biodiesel yield of 98.36% was accomplished. The fuel properties of biodiesel produced were found to be within the ASTM D6751 and DIN EN 14214 biodiesel specifications.
Keywords:Biodiesel; Sesame oil; Methanolysis; Optimization; Response surface methodology
The need for alternative fuel, which can solve the problems of finite nature of fossil fuels and their environmental concerns, is imperative. About 98% of carbon emissions result from fossil fuel combustion . The continued and increasing use of diesel derived from fossil oil intensifies air pollution and magnifies the global warming problem caused by CO2 and other particulate matters . Biodiesel, which is renewable in nature, non-toxic, and environmentally friendly with higher octane number than diesel, is the most promising candidate. It can be defined as the alkyl esters of long-chain fatty acids derived from renewable feedstock like vegetable oils or animal fats, obtained by transesterifying oil or fat with an alcohol .
Despite the food vs. fuel challenge facing the world, both edible and non-edible vegetable oil have been successfully used to produce biodiesel with superior qualities to diesel. These feedstocks include soybean oil , moringa oil , sunflower oil , rapeseed oil , palm oil , jatropha oil , camelina oil , neem oil [11,12], karanja oil , and sesame oil . The ideal vegetable oil for biodiesel must be readily available, its plant should be easy to cultivate, and its composition must include a high percentage of monounsaturated fatty acids (C16:1, C18:1), a low proportion of polyunsaturated acids (C18:2, C18:3) and a controlled amount of saturated fatty acids (C16:0, C18:0) .
Transesterification of oil or fat to produce a high yield of biodiesel is typically investigated by optimizing the following reaction variables: alcohol/oil molar ratio, catalyst concentration, reaction temperature, and reaction time involved in the process. Both the effects of the variables and their reciprocal interactions have been evaluated using response surface methodology (RSM). RSM is a useful statistical tool, which has been applied in research for optimizing various processes including pretreatment of the Jatropha curcas fruit hull for enzymatic hydrolysis , bioremediation of contaminated soil , and transesterification reaction of vegetable oils. Transesterification of Moringa oleifera oil to biodiesel was optimized via RSM by Rashid et al. , whereas Tiwari et al.  optimized biodiesel production from jatropha oil with high free fatty acids using RSM. Jeong et al.  optimized the transesterification process of animal fat ester, while Yuan et al.  used RSM for optimization of waste rapeseed oil with high FFA to biodiesel. Also, Zhang et al.  applied RSM for the optimization of Zanthoxylum bungeanum seed oil transesterification to biodiesel using CaO as catalyst. Optimization of transesterification variables for biodiesel production from cottonseed oil using RSM has also been reported .
Sesamum indicum, commonly called beniseed in Africa but globally known as sesame, is a fast-growing crop whose oilseed can be utilized for biodiesel production through transesterification process. It is an important oilseed crop cultivated in many parts of the world. Its seed is composed of 50% to 52% oil, 17% to 19% protein, and 16% to 18% carbohydrate . Fatty acid compositions of the oil are mostly oleic (32.7% to 53.9%), linoleic (39.3% to 59%), palmitic (8.3% to 10.9%), and stearic (3.4% to 6.0%) acids . However, the fatty acid composition of this oil varies considerably among the different varieties worldwide [22,23]. Sesame oil has emerged as a potential feedstock for biodiesel production. Although the transesterification of this seed oil from Turkey and Pakistan cultivars to biodiesel has been reported [14,24,25], the process has not yet been optimized using a statistical approach such as response surface methodology, and the maximum yield reported was 92% .
In this present study, an effort was made to optimize the process conditions for transesterification of sesame oil from a Nigeria cultivar to biodiesel. RSM with central composite rotatable design (CCRD) was applied, and the effect of four variables (reaction temperature, catalyst amount, reaction time, and methanol/oil molar ratio) at five levels and their reciprocal interactions were evaluated.
S. indicum samples were collected from Lokoja, Nigeria. Chaffs were separated from the oilseeds by winnowing. The cleaned oilseeds were milled to powder by grinding with a plate machine. All chemicals and reagents used for this work were of analytical grades.
Extraction of sesame oil
Soxhlet apparatus and n-hexane as solvent were used for the oil extraction. Initially, the apparatus was charged with a known weight of the oilseed powder in a muslin cloth placed in a thimble of the apparatus. A round bottom flask containing a known volume of n-hexane was fixed to the end of the apparatus, and a condenser was tightly fixed at the bottom end of the extractor. The whole setup was heated up in a heating mantle at a temperature of 70°C. The excess solvent in the oil was recovered by distillation.
Experimental design of biodiesel production
In this work, the central composite experimental design was employed to optimize the transesterification of sesame oil to biodiesel. A five-level-four-factor design was applied, which generated 30 experimental runs. These included 16 factorial points, 8 axial points, and 6 central points to provide information regarding the interior of the experimental region, making it possible to evaluate the curvature effect. Selected variables for the transesterification process from the sesame oil were reaction temperature (X1), catalyst amount (X2), reaction time (X3), and methanol/oil molar ratio (X4). The coded levels of the independent variables are given in Table 1. The experiments were randomized to minimize the effects of unexplained variability in the observed response due to extraneous factors.
Table 1. Factors and their levels for composite central design
Sesame biodiesel production procedure
The base-catalyzed reaction was carried out according to the work of Schinas et al.. The transesterification reaction of sesame oil was carried out in a glass reactor (1-L capacity). The reactor was charged with a known weight of the oil sample, which was preheated to the desired temperature. A known weight of NaOH pellet (catalyst) was dissolved in a specific volume of methanol and was transferred into the preheated oil on the hot plate with a magnetic stirrer. The resulting mixture was transferred into a separating funnel for glycerol and methyl ester separation. The glycerol was allowed to separate by gravity for 24 h. Two phases separated clearly, the less dense methyl ester at the top and the denser phase of glycerol at the bottom. The glycerol was tapped off the bottom of the separating funnel, leaving behind methyl ester in the separating funnel. The methyl ester was washed with distilled water to remove residual catalyst, glycerol, methanol, and soap. The methyl ester was further dried over a heated CaCl2 powder. The final product was the sesame biodiesel, and the yield was calculated gravimetrically as described in Equation 1:
Statistical data analysis
Sesame biodiesel production data were analyzed statistically using response surface methodology, in order to fit the quadratic polynomial equation generated by the Design Expert software version 126.96.36.199 (Stat-Ease Inc., Minneapolis, MN, USA). To correlate the response variable to the independent variables, multiple regressions were used to fit the coefficient of the polynomial model of the response. The quality of the fit of the model was evaluated using a test of significance and analysis of variance (ANOVA). The fitted quadratic response model is described by
where Y is the response variable (sesame biodiesel yield), bo is the intercept value, bi (i = 1, 2,…, k) is the first-order model coefficient, bij is the interaction effect, and bii represents the quadratic coefficients of Xi, and e is the random error.
Oil and fuel properties
Fuel properties, namely, moisture content, density at 15°C, kinematic viscosity at 40°C, iodine value, acid value, higher heating value, flash point, pour point, cloud point, cetane number, diesel index, API, and aniline point of both sesame oil and biodiesel, were determined following standard methods and compared with American and European standards (ASTM D 6751 and DIN EN 14214).
Analysis of fatty acid compositions of sesame biodiesel
Fatty acid profile of the sesame biodiesel was determined using gas chromatography (HP 6890 powered with HP ChemStation Rev. A 09.01 (1206) software, Agilent Technologies Inc., Santa Clara, CA, USA). The sesame biodiesel sample (50 mg) was esterified for 5 min at 95°C with 3.4 mL of the 0.5 M KOH in dry methanol. The mixture was neutralized using 0.7 M HCl, and 3 mL of 14% boron triflouride in methanol was added. The mixture was heated for 5 min at a temperature of 90°C to achieve complete methylation process. The fatty acids were thrice extracted from the mixture with redistilled n-hexane. The content was concentrated to 1 μL for gas chromatography analysis, and 1 μL was injected into the injection port of gas chromatography (GC).
Results and discussion
Properties of the sesame seed oil used as feedstock
The analysis of the seed oil showed that it has a moisture content of 0.09%, specific gravity of 0.88, and kinematic viscosity of 31.39 mm2/s. The acid value of the oil was 0.50 mg KOH/g oil, while the iodine value was 108 g I2/100 g oil. Whereas the saponification value and the peroxide value of the oil were 190 mg KOH/g oil and 7.80 meq O2/kg oil, respectively, its higher heating value and cetane number were 40.02 MJ/kg and 50.73, respectively. These results are within earlier reported values in the literature [14,27].
Optimization of sesame biodiesel production
The coded factors considered in this work with experimental results, predicted values, as well as the residual values obtained, are shown in Table 2. The parity plot showed a satisfactory correlation between the actual and predicted values of sesame biodiesel yield (Figure 1). Table 3 showed the results of the test of significance for every regression coefficient. The results revealed that the p value of the model terms was significant, i.e., p < 0.05. In this particular case, the four linear terms (X1, X2, X3,X4), six cross products (X1X2, X1X3, X1X4, X2X3, X2X4, X3X4), and the four quadratic terms (X12, X22, X32, and X42) were all remarkably significant model terms at 95% confidence level. Based on the large F values and low corresponding p values, all the linear terms had very strong effects on the biodiesel yield. However, methanol/oil molar ratio X4 (F value = 9831.27) was the most significant variable followed by the catalyst amount X2 (F value = 3364.87), followed by reaction time X3 (F value = 1835.80), and lastly, reaction temperature X1 (F value = 1192.78). In order to minimize error, all the coefficients were considered in the design. The model coefficients and probability values, i.e., coded values, are shown in Table 4. The model also proved suitable for adequate representation of the real relationship among the selected independent variables (X1, X2, X3,X4). Results of the ANOVA of the regression equation are presented in Table 5. The model F value of 3,234.83 with a low probability value (p < 0.0001) implied a high significance for the regression model . The goodness of fit of the model was checked by the coefficient of determination (R2), which should be at least 0.80 for the good fit of a model . The R2 value of 0.9995 indicated that the sample variation of 99.95% for biodiesel was attributed to the independent variables, and only 0.05% of the total variations are not explained by the model. The value of the adjusted determination coefficient (adjusted R2 of 0.9989) was also very high, supporting a high significance of the model [29,30], and all p values were less than 0.0001 except for X3X4 with a p value of 0.0012, which indicated the model proved suitable for the adequate representation of the actual relationship among the selected variables. The lack-of-fit term of 0.9589 was not significant relative to the pure error. The final equation in terms of coded factors for the central composite response surface quadratic model is expressed in Equation 3.
Table 2. CCRD and experimental, predicted, and residual values for five-level four-factor response surface analysis
Figure 1. Parity plot showing the distribution of actual vs. predicted values of sesame biodiesel yield.
Table 3. Test of significance for every regression coefficient of CCRD
Table 4. ANOVA of regression equation
Table 5. ANOVA for response surface quadratic model for intercept
All the linear terms and cross-product term X1X3 had a positive effect on the sesame biodiesel yield, while the rest of the surface quadratic model terms had negative effects on the yield. The low values of standard error observed in the intercept and all the model terms demonstrated that the regression model fits the data well, and the prediction was good (Table 5). The variance inflation factor (VIF) obtained in this study showed that the center points are orthogonal to all other factors in the model.
Generally, the 3D response surface plots are graphical representations of the regression equation for the optimization of the reaction variables, and they are presented in Figure 2. The curvatures’ nature of the three-dimensional surfaces in Figure 2a,c,f suggested mutual interaction of catalyst amount with reaction temperature, methanol/oil molar ratio with reaction temperature, and methanol/oil molar ratio with reaction time, respectively. On the other hand, the curvatures’ nature of the three-dimensional surfaces in Figure 2b,d,e indicated moderate interactions of reaction time with reaction temperature, reaction time with catalyst amount, and methanol/oil molar ratio with catalyst amount, respectively.
Figure 2. Surface plots for conversion of sesame seed oil to biodiesel. (a) catalyst amount and reaction temperature, (b) reaction time and reaction temperature, (c) methanol/oil molar ratio and reaction temperature, (d) reaction time with catalyst amount, (e) methanol/oil molar ratio and catalyst amount, and (f) methanol/oil molar ratio and reaction time.
The optimal values of the independent variables selected for the transesterification of sesame biodiesel were obtained by solving the regression equation (Equation 2) using the Design Expert software. The optimal conditions for this process were statistically predicted as X1 = 63°C, X2 = 1.04 wt.%, X3 = 51.09 min, and X4 = 6.24. The predicted sesame biodiesel yield under the above set of conditions was Y = 99.71%. In order to verify the prediction of the model, the optimal conditions were applied to three independent replicates, and the average biodiesel yield obtained was 98.63%, which is well within the predicted value for the model equation.
Sesame biodiesel quality characterization and fuel properties
The results of GC analysis revealed that sesame biodiesel is mainly composed of methyl esters of linoleic acid (45.17%) and oleic acid (39.02%), palmitic acid (8.47%) followed by 5.15% stearic acid, 0.93% behenic acid, 0.77% arachidic acid, 0.15% linolenic acid, 0.09% myristic acid, 0.04% palmitoleic acid, 0.03% arachidonic acid, and 0.01% erucic acid. These observations further showed that the sesame biodiesel produced under the optimal condition was within ASTM D 6751 specifications for free and total glycerol contents. In addition, Table 6 shows the results of the evaluation of sesame seed oil and fuel properties of the sesame biodiesel. Most of the fuel properties of sesame biodiesel from this study compared well with ASTM D6751 and EN 14214 biodiesel specifications.
Table 6. Properties of sesame biodiesel in comparison to biodiesel standards
The results obtained in this study using response surface methodology to determine the effects of four reaction variables, namely, oil/methanol molar ratio, reaction temperature, catalyst concentration, and reaction time on sesame biodiesel yield via methanolysis, indicate that the biodiesel produced had fuel properties that satisfied both the ASTM D 6751 and EN 1424 standards. The statistical model predicted the optimal conditions for the selected transesterification variables as reaction temperature of 63°C, catalyst amount of 1.04 wt.%, reaction time of 51.09 min, and methanol/oil molar ratio of 6.24 with an actual sesame biodiesel yield of 98.63%. Hence, the work established the usefulness of RSM for optimum conversion of sesame (S. indicum) seed oil to biodiesel and also suggests that sesame oil could be effectively used as feedstock for biodiesel production.
Authors declare that they have no competing interests.
EB conceived and designed the research work. He is the major contributor in the preparation of the manuscript. TFA was involved in the laboratory work, and he drafted the manuscript. Both authors read and approved the final manuscript.
EB is a senior lecturer in the Department of Chemical Engineering, Obafemi Awolowo University, Ile-Ife, Nigeria. He holds B.Sc. and M.Sc. degrees in Chemical Engineering from Obafemi Awolowo University. He also received a Ph.D. degree in Biotechnology from Technical University Carolo-Wilhelmina, Braunschweig, Germany in 2005. He was a DAAD scholar and a guest scientist at the National Research Centre for Biotechnology, Germany between April 2001 and December 2004, respectively. His research focus is in the area of biofuel development and bioconversion of agricultural crops and wastes to useful products. TFA is a lecturer in the Department of Chemical Engineering, Landmark University, Nigeria. He holds a M.Sc. degree in Chemical Engineering from Obafemi Awolowo University, Ile-Ife, Nigeria.
E. Betiku gratefully acknowledged the equipment donation by the World University Service, Germany and provision of relevant literature by the DAAD. The authors also thank A.K. Omole and O.E. Aluko for their technical support.
Shay, EG: Diesel fuel from vegetable oils: status and opportunities. Biomass Bioenerg.. 4(4), 227–242 (1993). Publisher Full Text
Knothe, G: Biodiesel and renewable diesel: a comparison. Prog. Energ. Combust.. 36, 364–373 (2010). Publisher Full Text
Rashid, U, Anwar, F, Ashraf, M, Saleem, M, Yusup, S: Application of response surface methodology for optimizing transesterification of Moringa oleifera oil: biodiesel production. Energ. Convers. Manage.. 52, 3034–3042 (2011). Publisher Full Text
Berrios, M, Siles, J, Martin, M, Martin, A: A kinetic study of the esterification of free fatty acids (FFA) in sunflower oil. Fuel. 86, 2383–2388 (2007). Publisher Full Text
Yuan, X, Liu, J, Zeng, G, Shi, J, Tong, J, Huang, G: Optimization of conversion of waste rapeseed oil with high FFA to biodiesel using response surface methodology. Renew. Energ.. 33, 1678–1684 (2008). Publisher Full Text
Margaretha, YY, Prastyo, HS, Ayucitra, A, Ismadji, S: Calcium oxide from Pomacea sp. shell as a catalyst for biodiesel production. Int. J. Energy Env. Eng.. 3, 33 (2012). BioMed Central Full Text
Tiwari, AK, Kumar, A, Raheman, H: Biodiesel production from jatropha oil (Jatropha curcas) with high free fatty acids: an optimized process. Biomass Bioenerg.. 31, 569–575 (2007). Publisher Full Text
Wu, X, Leung, DYC: Optimization of biodiesel production from camelina oil using orthogonal experiment. Appl. Energy. 88(11), 3615–3624 (2011). Publisher Full Text
Muthu, H, SathyaSelvabala, V, Varathachary, T, Kirupha Selvaraj, D, Nandagopal, J, Subramanian, S: Synthesis of biodiesel from neem oil using sulfated zirconia via transesterification. Braz J Chem Eng. 27(4), 601–608 (2010)
Aransiola, EF, Betiku, E, Ikhuomoregbe, DIO, Ojumu, TV: Production of biodiesel from crude neem oil feedstock and its emissions from internal combustion engines. Afr. J. Biotech.. 11(22), 6178–6186 (2012)
Meher, LC, Dharmagadda, VSS, Naik, SN: Optimization of alkali-catalyzed transesterification of Pongamia pinnata oil for production of biodiesel. Bioresour Technol. 97(12), 1392–1397 (2006). PubMed Abstract | Publisher Full Text
Saydut, A, Duz, MZ, Kaya, C, Kafadar, AB, Hamamci, C: Transesterified sesame (Sesamum indicum L.) seed oil as a biodiesel fuel. Bioresour. Technol. 99(14), 6656–6660 (2008). PubMed Abstract | Publisher Full Text
Wang, R, Hanna, MA, Zhou, W-W, Bhadury, PS, Chen, Q, Song, B-A, Yang, S: Production and selected fuel properties of biodiesel from promising non-edible oils: Euphorbia lathyris L., Sapium sebiferum L. and Jatropha curcas L. Bioresour. Technol.. 102(2), 1194–1199 (2011). PubMed Abstract | Publisher Full Text
Marasabessy, A, Kootstra, AMJ, Sanders, JPM, Weusthuis, RA: Dilute H2SO4-catalyzed hydrothermal pretreatment to enhance enzymatic digestibility of Jatropha curcas fruit hull for ethanol fermentation. Int. J. Energy Env. Eng.. 3, 15 (2012). BioMed Central Full Text
Agarry, SE, Ogunleye, OO: Box-Behnken design application to study enhanced bioremediation of soil artificially contaminated with spent engine oil using biostimulation strategy. Int. J. Energy. Env. Eng.. 3, 31 (2012). BioMed Central Full Text
Zhang, J, Chen, S, Yang, R, Yan, Y: Biodiesel production from vegetable oil using heterogenous acid and alkali catalyst. Fuel. 89(10), 2939–2944 (2010). Publisher Full Text
Fan, X, Wang, X, Chen, F: Biodiesel production from crude cottonseed oil: an optimization process using response surface methodology. The Open Fuels Energy Sci. J.. 4, 1–8 (2011). Publisher Full Text
Yeramanos, DM, Hemstreet, S, Saleeb, W, Huszar, CK: Oil content and composition of the seed in the world collection of sesame introductions. J. Am. Oil Chem. Soc.. 49, 20–23 (1972). Publisher Full Text
Schinas, P, Karavalakis, G, Davaris, C, Anastopoulos, G, Karonis, D, Zannikos, F, Stournas, S, Lois, E: Pumpkin (Cucurbita pepo L.) seed oil as an alternative feedstock for the production of biodiesel in Greece. Biomass Bioenerg. 33, 44–49 (2009). Publisher Full Text
Elleuch, M, Besbes, S, Roiseux, O, Blecker, C, Attia, H: Quality characteristics of sesame seeds and by-products. Food Chem. 103(2), 641–650 (2007). Publisher Full Text
Guan, X, Yao, H: Optimization of viscozyme L-assisted extraction of oat bran protein using response surface methodology. Food Chem. 106, 345–351 (2008). Publisher Full Text