Iraqi Journal of Chemical and Petroleum Engineering Vol.14 No.2 (June 2013) 13- 20 ISSN: 1997-4884 Kinetic Study and Simulation of Oleic Acid Esterification in Different Type of Reactors Ammar S. Abbas and Sura M. Abbas Chemical Engineering Department, College of Engineering – University of Baghdad Abstract Esterification reaction is most important reaction in biodiesel production. In this study, oleic acid was used as a suggested feedstock to study and simulate production of biodiesel. Batch esterification of oleic acid was carried out at operating conditions; temperature from 40 to 70 °C, ethanol to oleic acid molar ratio from 1/1 to 6/1, H2SO4 as the catalyst 1 and 5% wt of oleic acid, reaction time up to 180 min. The optimum conditions for the esterification reaction were molar ratio of ethanol/oleic acid 6/1, 5%wt H2SO4 relative to oleic acid, 70 °C, 90 min and conversion of oleic 0.92. The activation energy for the suggested model was 26625 J/mole for forward reaction and 42189 J/mole for equilibrium constant. The obtained results simulated to other types of reactors with different operating conditions using reactop cascade package. The conversion of oleic acid of simulation results at optimum operating conditions was 0.97 for isothermal batch and plug flow reactors, 0.67 for isothermal CSTR, while the conversions of oleic acid in the adiabatic mode were 0.82, 0.40, 0.74 for batch, CSTR, PFR reactors respectively. Key word: Oleic acid, Esterification, Kinetic study, Simulation. Introduction Energy sources can be classified into three groups fossil, fissile, and renewable. The term fossil refers to an earlier geological age. Fossil fuels were formed many years ago and are not renewable. The fossil energy sources are petroleum, coal, bitumen, natural gas, oil shale, and tar sands. The main fissile energy sources are uranium and thorium [1]. One of the main and most important non-renewable energy sources is oil, the production rate of oil is expected to peak in the next few years therefore, the world needs an alternative and renewable energy source that are able to meet world energy needing [2]. The renewable energy sources such as biomass, hydro, wind, solar (thermal and photovoltaic), geothermal, marine, and hydrogen will play an important role in the future. It is predicted that in 2025, approximately half of the global energy supply will come from renewable energy, and electricity generation from renewable will be more than 80% of the total global electricity supply [3, 4]. Biomass can be converted into liquid and gaseous fuels through thermo- chemical and biological methods. Biofuel is a non-polluting, locally available, accessible, sustainable, and reliable fuel obtained from renewable sources [5]. Liquid biofuels fall into Iraqi Journal of Chemical and Petroleum Engineering University of Baghdad College of Engineering Kinetic Study and Simulation of Oleic Acid Esterification in Different Type of Reactors 14 IJCPE Vol.14 No.2 (June 2013) -Available online at: www.iasj.net the following categories: (a) vegetable oils and biodiesels, (b) alcohols, and (c) bio-crude and synthetic oils [1]. Esterification is one of the most important reactions in chemical industry; and the relative reaction is usually conducted by using, in general, homogeneous acid catalysts [6]. Esterification reaction carried out between carboxylic acids (free fatty acids) and alcohols with the presence of an acid catalyst (such as sulfuric acid, orqanic sulfonic acid, and hydrochloric acid) to form ester and water (see equatin 1) [7]. ↔ …(1) FFA alcohol Ester Water The present work, intend to produce ethyl ester (biodiesel) from the oleic acid by esterification process. The effects of reaction temperature, amount of catalyst, and the molar ratio alcohol/oil were considered. Also, the kinetic of oleic acid esterification were studied and the results were simulated to different type of reactors and operating conditions. Experimental Work Materials 1. Oleic acid, obtained from local markets. The specific gravity of the oleic acid is 0.895, BDH Chemicals Ltd. 2. Ethyl alcohol obtained from local markets with specific gravity is 0.7692 (88 to 90 wt. %). 3. Sulfuric acid as an acid catalyst obtained from local markets, the purity of this acid is 98% (Sp.Gr. is 1.84) Aldrich. 4. Sodium hydroxide for titration (RIEDEL_DEHAEN AG SEELZE_HANNOVER Chem. rian, plozchen, DAB7, B.P.1968 M.Wt. 40). 5. Phenolphthalein (as indicator), Fluka. Apparatus The apparatus used in this study for esterification reaction is shown in Figure 1. The batch scale system consists of the followings: 1. Heat flat magnetic stirrer (Stuart (CB302)/USA. 2. Reflux Condenser (Germany). 3. Centrifuge (Griffin & George Loughborough/Britain) 4. Mercury thermometer from zero to 250 °C. 5. 3 Necks flask (500 ml). 6. Sensitive balance. Fig. 1, Schematic diagram of the reactor Esterification of Oleic Acid The esterification reaction carried out between acid and free fatty acid (FFA) to produce ester (biodiesel) and water. The system was maintained at atmospheric pressure and experiments were carried out at constant temperature. The agitation was kept constant at 300rpm. This process was studied at different percent of ethanol/oil mole ratio of 1/1, 3/1, 6/1 sulfuric acid as a catalyst of 1 and 5 wt% relative to oleic acid, reaction time up to 180 minutes and at different temperatures 40 to 70 ° C. The esterification reactor was loaded with 15 ml (13.43 g) of oleic acid, and the desired amount of ethanol. The mixer was agitated and preheated to Ammar S. Abbas and Sura M. Abbas -Available online at: www.iasj.net IJCPE Vol.14 No.2 (June 2013) 15 the desired temperature and then H2SO4 was added. At each period of time (15 or 30 minutes), 5 ml from the mixture reaction was taken and 2 drops of phenolphthalein was added as indicator and titrated with 1 molarity of NaOH in order to obtain oleic acid conversion[8]. Results and Discussion Effect of Operating Temperature and Time on the Oleic Acid Conversion Figures 2 and 3 show the oleic acid conversion with reaction time at various temperatures and ethanol/oleic acid molar ratio 6/1 with 1 and 5 wt.% of H2SO4. Time, min 0 20 40 60 80 100 120 140 160 180 200 O le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 T=40° C T=50° C T=60° C T=70° C EtOH/oil = 6/1 H2SO4 = 1% Fig. 2, effect of the reaction temperature on oleic acid conversion by esterification reaction, 6/1 ethanol/oil mol. ratio and 1 wt% H2SO4 Time, min 0 20 40 60 80 100 120 140 160 180 200 O le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 T=40° C T=50° C T=60° C T=70° C EtOH/oil = 6/1 H2SO4 = 5% Fig. 3, effect of the reaction temperature on oleic acid conversion by esterification reaction, 6/1 ethanol/oil mol. ratio and 5 wt% H2SO4 As shown in the fig 2 the conversion of oleic acid at 6/1 ethanol/oleic mol. ratio and 1% H2SO4 was 0.87 at 70 ° C after 150 minutes and was 0.92 at 70 ° C after 90 minutes with 5% H2SO4, which is the highest conversion achieved in this study (as shown in Fig.3). The oleic acid conversion increased with increasing temperature. Increasing the temperature may causes increase of molecule activity, which means that more molecules have more energy, thus, the possibility of molecule to react increased. These results are in a good agreement with the results of the production of fatty acid ethyl ester (FAEE) from oleic acid (FAA) with short-chain alcohols (ethanol, propanol and butanol) under ultrasonic irradiation reported by Hanh et al. [9] The higher conversion for the esterification process was 0.92 at the reaction temperature 60 ° C and 3/1 mol. ratio of ethanol to oleic acid with 5% H2SO4 after 2 hours irradiation. Effect of Ethyl Alcohol/Oil Mol. Ratio on the Oleic Acid Conversion Molar ratio of ethanol to oleic acid is one of the most important variables that affect the conversion of oleic acid. The experiments of esterification reaction of oleic acid with ethanol were carried out under various ethanol/oil mol. ratios (1/1 to 6/1). Figures 4 and 5 show the conversion of oleic acid with reaction time at various ethanol/oleic acid molar ratios using 1 and 5 wt. % H2SO4 (as catalyst) at 70 °C. It can be observed that the oleic acid conversion increased from 0.61 at a molar ratio of 1/1 to 0.87 at a molar ratio of 6/1 after 180 min using 1 wt% of H2SO4 at 70°C (Fig. 4). By using 5 wt% H2SO4 at 70°C, the oleic acid conversion increased from 0.87 at 1/1 of ethanol/oil mol. ratio after 180 min Kinetic Study and Simulation of Oleic Acid Esterification in Different Type of Reactors 16 IJCPE Vol.14 No.2 (June 2013) -Available online at: www.iasj.net to 0.92 at 6/1 ethanol/oil mol. ratio after 90 min, which is the highest conversion achieved in this study (in Fig.5). Time (min) 0 20 40 60 80 100 120 140 160 180 200 O le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 Etoh/oleic=1/1 Etoh/oleic=3/1 Etoh/oleic=6/1 Temp=70 o C H 2 SO 4 =1% Fig. 4, effect of the ethyl alcohol/oleic mol. ratio on oleic acid conversion by esterification reaction, at the temperature 70 ° C and 1 wt% H2SO4 Time (min) 0 20 40 60 80 100 120 140 160 180 200 O le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 Etoh/oleic=1/1 Etoh/oleic=3/1 Etoh/oleic=6/1 Temp=70 o C H 2 SO 4 =5% Fig. 5, effect of the ethyl alcohol/oleic mol. ratio on oleic acid conversion by esterification reaction, at the temperature 70 ° C and 5 wt% H2SO4 Increasing in ethanol/oleic mol. ratio causes increasing in the conversion in oleic acid, because the esterification of oleic acid with ethanol is an equilibrium-limited chemical reaction and the position of equilibrium controls the amount of ester formed [10]. These results are a good agreement with those obtained by Marchetti and Errazu [11], the highest conversion of FFA was 0.96 at the reaction conditions molar ratio of ethanol to acid oil was 6.126, reaction temperature 55° C, reaction time 240 min and 2.261 wt% H2SO4. Effect of Catalyst Amount on the Oleic Acid Conversion The amount of catalyst affect on the oleic acid conversion is given in Fig.6. It shows the relationships between the oleic acid conversion and time at various catalyst concentrations with ethanol/oleic acid molar ratio 6/1 and 70°C. The oleic acid conversion increased with increasing the catalyst concentration. Time (min) 0 20 40 60 80 100 120 140 160 o le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 1%wt H 2 SO 4 5% wt H 2 SO 4 Fig. 6, effect of the amount of catalyst concentration (H2SO4) on oleic acid conversion by esterification reaction, at the temperature 70 ° C and 6/1 ethanol /oleic molar ratio The reaction rate of esterification is directly proportional to the amount of a catalyst, so the catalyst is used to enhance the reaction rate and conversion. The amount of sulfuric acid employed as a catalyst is related to the formation of H + that catalyzes the reaction. Kinetic of Oleic Acid Esterification Kinetic obtained from laboratory unit are usually play an important role in modeling and scale up designs for new biodiesel production units. The data obtained by Differential method of analysis obtained data has Ammar S. Abbas and Sura M. Abbas -Available online at: www.iasj.net IJCPE Vol.14 No.2 (June 2013) 17 been used to find favorable kinetic model for oleic acid esterification. The differential method of analysis deals directly with the differential rate equation to be tested, evaluating all terms in the equation including the derivative -rA, and testing the goodness of fit of the equation with experiment. The conversion of oleic acid (x) fitted by the power degree formula according to the time at certain temperature, and then finding the rate coefficients (k1 and k2) as well as the orders of the reactants and products materials (n1, m1, n2 and m2) using least squares method for the suggested reaction kinetics equation (eqs. 2 to 6). ↔ ...(2) …(3) …(4) …(5) …(6) Where: k = rate coefficient of reaction. 1= forward reaction, 2= backward reaction, n and m are the reaction orders. Keq = equilibrium constant. CA, CB, CC, CD = reactants and products moles at any time. CAe, CBe, CCe, CDe = reactants and products moles at equilibrium time. According to Arrhenius law [12] a plot of Ln k vs. 1/T gives a straight line, with slope for E/R to find activation energy for forward reaction and equilibrium reaction (figs.7 and 8). The values of constants in the esterification reaction of oleic acid are summarized in Table 1. ln k 1 = -3202.3/T+5.28 1/T 0.0029 0.0030 0.0031 0.0032 ln k 1 -5.2 -5.0 -4.8 -4.6 -4.4 -4.2 -4.0 -3.8 Fig. 7, show the linear plot of ln k vs. 1/T resulting from the esterification of oleic acid for forward reaction ln k eq = -5074.5/T+16.87 1/T 0.0029 0.0030 0.0031 0.0032 ln k e q 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 Fig. 8, show the linear plot of ln keq vs. 1/T resulting from the esterification of oleic acid for equilibrium reaction Table 1, Constant values of the esterification reaction kinetic model Constant value Esterification reaction k°e 2.1*10 7 (J/mole) 42189 196.3 (J/mole) 26625 1.25 0.5 1 0.5 The obtained values of rate coefficients, equilibrium constant, orders of reactants and products materials, heat of reactions and activation energy were used to simulate the results. Predicted values calculated from empirical model and experimental data are shown in fig.9. Kinetic Study and Simulation of Oleic Acid Esterification in Different Type of Reactors 18 IJCPE Vol.14 No.2 (June 2013) -Available online at: www.iasj.net Real production rate, vol.%/min 0.000 0.002 0.004 0.006 0.008 0.010 P re d ic te d p ro d u c ti o n r a te , v o l. % /m in 0.000 0.002 0.004 0.006 0.008 0.010 Acid Fig. 9, Experimental and predicted values of apparent rate constant by using suggested model of esterification of oleic acid with 6/1 ethanol/oleic molar ratio, 5% H2SO4 as catalyst and temperature range from 40 to 70 ° C Statistical analysis of the experimental data shows that the solution of the model (eq. 6) is proportional with the experimental data. The distribution of the experimental data around the model solution of esterification (correlation coefficient (R)) is about 0.9061, standard deviation (S) is 6.9*10 -3 , and average relative error is 1.34% in 95 % confidence level, as summarized in Table 2. Table 2, statistical analysis of the model statistical analysis Esterification reaction Correlation coefficient (R) 0.9061 Standard deviation (S) 6.9 * 10 -3 Average relative error 1.34% Confidence level 95% Simulation Results Simulation for the behavior of batch, plug and mixed flow reactors in different operating modes (adiabatic and isothermal) was carried out at optimum conditions (70 °C, 6/1 ethanol/oleic mol ratio, 5% H2SO4 as catalyst) and the previously obtained results of kinetic models (Table 1) by using reactop cascade package to predict best reaction time. Fig. 10 shows the reactants and products moles with reaction time for simulated and experimental results. At the beginning of reaction the system has only liquid phase of reactant (oleic acid and ethanol with 12% water content). As time proceeds, reactant consumed, while water and ester formed. At the simulated results, the reactants consumed and products formed faster than the experimental results (as shown in Table 3). From this Table 3, simulated results are always higher than the experimental results which could be attributed to the effect of side reactions, which cause consumption of oleic acid to form undesirable products, such as water that reduced the formation of ester [8]. Time (min) 0 20 40 60 80 100 120 R e a c ta n ts a n d p ro d u c ts m o le s 0 1 2 3 4 5 Oleic (simulaltion) Ester (simulation) Ethanol (simulation) Water (simulation) 0 20 40 60 80 100 120 0 1 2 3 4 5 Oleic (calculated) Ester (calculated) Ethanol (calculated) Water (calculated) Fig. 10, show the reactants and products moles in experimental and simulation isothermal batch reactor with reaction time at best operating conditions Figs. 11 and 12 show the oleic acid conversion with reaction time in batch, PFR and CSTR operating in isothermal and adiabatic mode, respectively. The oleic acid conversion in isothermal batch and plug flow reactors were 0.97 after 90 min from the reaction, whereas the conversion in mixed flow reactor reached highest conversion 0.63 after 33 min and then the conversion dropped (fig. 11). Fig. 12 shows the oleic acid conversion with time in batch, PFR and CSTR in adiabatic mode. The Ammar S. Abbas and Sura M. Abbas -Available online at: www.iasj.net IJCPE Vol.14 No.2 (June 2013) 19 conversion of oleic acid were 0.84 and 0.75 after 90 min in a batch reactor and plug flow reactor, while the conversion in CSTR was 0.40 after 40 min and then dropped. Time (min) 0 20 40 60 80 100 120 o le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 isothermal batch isothermal CSTR isothermal PFR Batch and PFR responce Fig. 11, shows oleic acid conversion with time in batch, PFR and CSTR reactors in isothermal mode at the best operating conditions However, it is obvious that batch and PFR is substantially better than the CSTR for obtaining high conversions in both isothermal and adiabatic mode Time (min) 0 20 40 60 80 100 120 o le ic a c id c o n v e rs io n 0.0 0.2 0.4 0.6 0.8 1.0 adiabatic Batch adiabatic CSTR adiabatic PFR Fig. 12, shows the oleic acid conversion with time in batch, PFR and CSTR reactors in adiabatic mode at the best operating conditions Table 4 shows numerically the highest and best conversion in simulation and experimental results at different types of reactors (from figures 10 to 12). Table 3, higher results at the best conditions from experimental and simulated results Experimental results Simulation results Time Moles of reactants and products Time Moles of reactants and products oleic Ethanol Ester Water Oleic Ethanol Ester Water 0 1 4.44 0 1.15 0 1 4.44 0 1.15 15 0.3 3.74 0.7 1.85 15 0.22 3.64 0.79 1.94 30 0.23 3.67 0.77 1.92 30 0.08 3.51 0.92 2.07 60 0.13 3.57 0.87 2.02 60 0.03 3.47 0.96 2.12 90 0.08 3.52 0.92 2.07 90 0.02 3.46 0.97 2.12 Table 4, highest conversion in simulation and experimental at operating conditions (T=70 °C, ethanol/oleic molar ratio=6/1, H2SO4 =5%wt) Reactor type Operating mode Time (min) conversion Note Batch isothermal 90 0.92 Experimental Batch Isothermal 90 0.97 Simulation Batch Adiabatic 90 0.82 Simulation CSTR Isothermal 20 0.67 Simulation CSTR Adiabatic 40 0.40 Simulation PFR Isothermal 90 0.97 Simulation PFR Adiabatic 90 0.74 Simulation Conclusion 1. The maximum conversion of oleic acid was 0.92 at 6/1 of ethanol/oleic acid mol. ratio after 90 minutes, 70° C and 5% H2SO4. 2. Increasing the ethanol/oleic acid molar ratio from 1 to 6 increases the conversion of oleic acid. 3. The activation energies for the suggested kinetic model were 26625 Kinetic Study and Simulation of Oleic Acid Esterification in Different Type of Reactors 20 IJCPE Vol.14 No.2 (June 2013) -Available online at: www.iasj.net J/mole for the forward reaction and were 42189 J/mole for equilibrium constant. 4. The conversion of oleic acid in the simulated model in batch and plug flow reactors is higher than the conversion in CSTR reactor in each isothermal and adiabatic mode. Nomenclatures -rA = dx/dt reaction rate of oleic acid A = oleic acid , B= ethanol, C= ester , D= water k = rate coefficient of reaction. 1= forward reaction, 2= backward reaction, n and m are the reaction orders. n1 for oleic , m1 for ethanol, n2 for ester, m2 for water Keq = equilibrium constant. CA, CB, CC, CD = reactants and products moles at any time. CAe, CBe, CCe, CDe = reactants and products moles at equilibrium time. E1 = activation energy for forward reaction (J/mole) Eeq = activation energy for equilibrium reaction (J/mole) ko =frequency factor ((mol) -0.74 /minute) for forward reaction koe = frequency factor for equilibrium reaction References 1. Demirbas A., "Energy priorities and new energy strategies", Energy Education Science Technology, Vol.109, pp.16-53, 2006. 2. Hodgson P.E., a book of "Energy, The Environment and Climate Change", published by World Scientific Publishing Co. Pte. Ltd. ISBN-13 978-1-84816-415-4, pp.1- 9, 2010. 3. EWEA (European Wind Energy Association), "Large scale integration of wind energy in the European power", Issues and recommendations, Brussels, 2005. 4. EREC (European Renewable Energy Council), "Renewable energy scenario by 2040", EREC Statistics, Brussels, 2006. 5. Vasudevan P., Sharma S., Kumar A., "Liquid fuel from biomass: an overview", Journal of scientific & Industrial Research, Vol.64, pp.822- 831, 2005. 6. Yin P., Chen L., Wang Z., Qu R., Liu X., Ren S., "Production of biodiesel by esterification of oleic acid with ethanol over organophosphonic acid- functionalized silica", Bioresource Technology, Vol.110, pp.258-263, 2012. 7. Khan A. K., "Biodiesel Kinetics & Catalyst Development", Thesis, University of Queensland, 2002. 8. Abbas S. M., "Kinetic study and modeling of oleic acid esterification with ethanol in presence of water", Thesis, University of Baghdad, 2013. 9. Hanh H.D., Dong N., Okitsu K., Nishimur R., Maeda Y., ''Biodiesel production by esterification of oleic acid with short-chain alcohols under ultrasonic irradiation condition'', Renewable Energy, Vol.34, pp.780–783, 2008. 10. Zubir M.I., Chin S.Y., "Kinetic of Modified Zirconiz-catalyzed Heterogeneous Esterification Reaction for Biodiesel Production", Journal of Applied Sciences, Vol.10(21), pp.2584-2589, 2010. 11. Marchetti J.M., Errazu A.F., "Esterification of free fatty acids using sulfuric acid as catalyst in the presence of triglycerides" ,Biomass and Bioenergy, Vol.32, pp.892-895, 2008 12. Levenspiel O., a handbook of ''Chemical reaction engineering'', published by John Wiley and Sons, third edition, pp.14- 72, 1999.