The kinetic of the transesterification reaction over the CsXH3−XPW12O40/Fe-SiO2 catalyst prepared using sol-gel and impregnation procedures was investigated in different operational conditions. Experimental conditions were varied as follows: reaction temperature 323–333 K, methanol/oil molar ratio = 12/1, and the reaction time 0–240 min. The H3PW12O40 heteropolyacid has recently attracted significant attention due to its potential for application in the production of biodiesel, in either homogeneous or heterogeneous catalytic conditions. Although fatty acids esterification reaction has been known for some time, data is still scarce regarding kinetic and thermodynamic parameters, especially when catalyzed by nonconventional compounds such as H3PW12O40. Herein, a kinetic study utilizing Gc-Mas in situ allows for evaluating the effects of operation conditions on reaction rate and determining the activation energy along with thermodynamic constants including ΔG, ΔS, and ΔH. It indicated that the CsXH3−XPW12O40/Fe-SiO2 magnetic nanocatalyst can be easily recycled with a little loss by magnetic field and can maintain higher catalytic activity and higher recovery even after being used 5 times. Characterization of catalyst was carried out by using scanning electron microscopy (SEM), X-ray diffraction (XRD), Fourier transform-infrared spectroscopy (FT-IR), N2 adsorption-desorption measurements methods, thermal gravimetric analysis (TGA), and differential scanning calorimetry (DSC).
1. Introduction
Biodiesel is a fuel composed of monoalkyl esters of long chain fatty acids derived from renewable sources, such as vegetable oils and animal fats [1]. Due to oxygen content, biodiesel is a clean, nontoxic, and biodegradable fuel with low exhaust emissions, without the sulphur and carcinogen content [2]. Biodiesel is prepared via reaction between triglycerides and alcohol in the presence of a catalyst [3]. Transesterification reaction of oil and alcohol with homogeneous catalyst is the most common method for the preparation of biodiesel [4, 5]. However, the homogeneous catalyst has many drawbacks, such as the difficulty in product separation, and equipment corrosion requirement of large quantity of water, environmental pollution [6]. The use of heterogeneous catalysts to replace homogeneous ones is easily regenerated and has a less corrosive nature, leading to safer, cheaper, and more environment-friendly operations, [7]. Currently, the research is focused on sustainable solid acid catalysts for transesterification reaction. In addition, it is believed that solid acid catalysts have the strong potential than liquid acid catalyst [8]. HPW catalyst is a strong Brønsted acid, with high thermal stability and high solubility in polar solvents, and it has shown to be more active in FFA esterification reactions than mineral acid catalysts [8–12]. However, heteropolyacids have low surface areas and high solubility in polar solvents. These features make some difficulties in catalyst recovery and catalyst lifetime [13]. Therefore, salts of HPAs with large single valence ions, such as Cs+, NH4+, and Ag+, have attracted much fondness because these salts will increase in surface area and profound changes in solubility over the parent heteropolyacid [14, 15]. HPW can be supported on several kinds of support such as SiO2, Al, ZrO2, activated carbon, SiO2-Al, and MCM-41 that SiO2 is cheap, easily available, and easily surface modifiable [16–18].
In the present kinetic study, pseudo-first-order model was applied to correlate the experimental kinetic data of catalytic performance of the CsXH3-XPW12O40/Fe-SiO2 for sunflower oil transesterification with methanol and the main thermodynamic parameters such as the effects of temperature on reaction rate, activation energy, entropy variation (ΔS), and enthalpy variation (ΔH) were determined; moreover, the effects of temperature on reaction rate and the order of reaction were assessed. Characterization of catalyst was carried out by using scanning electron SEM, XRD, FT-IR, N2 adsorption-desorption measurements methods, TGA, and DSC methods.
2. Experimental2.1. Fe-SiO2 Support Preparation
All materials with analytical purity were purchased from Merck and used without further purification. Ferric nitrate (Fe(NO3)3·9H2O) and tetraethyl orthosilicate (TEOS) were selected as the source of ferric and silica, respectively. A typical procedure for the preparation of ferric-silica mixed oxide containing 60 wt% ferric was followed. Firstly, 30 mL TEOS was mixed with certain amount anhydrous ethyl alcohol (C2H5OH). Secondly, 35.234 gr Fe(NO3)3·9H2O was dissolved with certain amount of anhydrous ethyl alcohol under stirring; also 30 gr of oxalic acid was dissolved in certain amount of anhydrous alcohol under stirring. In the final step, Fe and Si sols and oxalic acid were added simultaneously into the beaker under constant stirring to obtain a gel form. After the end of the above operations, the samples were aged for 90 min at 50°C. The obtained gel was dried in an oven (120°C, 12 h) to give a material denoted as the catalyst precursor. Finally, the catalyst precursor was calcined at 600°C for 6 h to produce magnetic solid catalyst. The Fe-SiO2 supported 12-tungstophosphoric acid catalyst with 4 wt.% of aqueous solution HPW (based on the Fe-SiO2 weight) were prepared by incipient wetness impregnation. The impregnated precursor was dried at 120°C for overnight and calcined at 600°C for 6 h.
Finally, the promoted catalyst by Cs was synthesized with Cs/H3PW12O40 = 2 wt.%. The salt obtained by this procedure will be designated hereinafter as CsXH3-XPW12O40/Fe-SiO2, where X is the amount of protons per [PW12O40]3− anion in the salt.
2.2. Characterization of Catalyst2.2.1. N2 Physisorption Measurements
The specific surface area, total pore volume, and the mean pore diameter were measured using a N2 adsorption-desorption isotherm at liquid nitrogen temperature (−196°C), using a NOVA 2200 instrument (Quantachrome, USA). Prior to the adsorption-desorption measurements, all the samples were degassed at 110°C in a N2 flow for 2 h to remove the moisture and other adsorbates.
2.2.2. Scanning Electron Microscopy (SEM)
The morphology of catalyst and precursor was observed by means of an S-360 Oxford Eng scanning electron microscopy.
Fourier transform infrared (FT-IR) spectra of the samples were obtained using a Bruker Vector 22 spectrometer in the region of 400–4000 cm−1.
2.2.4. X-Ray Diffraction (XRD)
X-ray diffraction (XRD) patterns of the catalysts were recorded on a diffractometer using CuKα radiation. The intensity data were collected over a 2θ range of 15–75.
2.2.5. Thermal Gravimetric Analysis (TGA) and Differential Scanning Calorimetry (DSC)
The TGA and DSC were carried out using simultaneous thermal analyzer (PerkinElmer) under a flow of dry air with a flow rate of 50 mL min−1. The temperature was raised from 20 to 700°C using a linear programmer at a heating rate of 3°C min−1.
2.3. Catalytic Tests
The type and quantity of methyl esters in the biodiesel samples were determined using gas chromatography-mass spectrometry (GC Agilent 6890N model and Mass Agilent 5973N model) equipped with a flame ionized detector (FID). A capillary column (HP-5) with column length (60 m), inner diameter (0.25 mm), and 0.25 μm film thickness was used with helium as the carrier gas. The temperature program for the biodiesel samples started at 50°C and ramped to 150°C at 10°C min−1. The temperature was held at 150°C for 15 min and ramped to 280°C at 5°C min−1. The holding time at the final temperature (250°C) was 5 min. Also, the injector was used from kind split/splitless.
2.4. Kinetic Studies
The transesterification of sunflower oil was carried out in a round bottomed flask fitted with a condenser and magnetic stirring system. The reaction system was heated to selected temperature in 50, 55, and 60°C. When the oil reached selected temperature, methanol/oil molar ratio = 12/1 and the catalyst amount (3 wt% related to oil weight) were added with continuous stirring (500 rpm). After completion of the reaction time (0–240 min), the sample concentration is calculated according to mole fraction at any time. The results can be seen in Table 1.
Concentration of methyl ester (based on mole fraction) at temperatures of 60 (A), 55 (B), and 50°C (C), respectively.
Reaction time (s)
XME
−ln (1-XME)
A
B
C
A
B
C
0
0.00
0.00
0.00
0.000000
0.000000
0.000000
30
0.32
0.24
0.14
0.385660
0.274437
0.150823
60
0.40
0.28
0.15
0.510826
0.328504
0.162519
90
0.42
0.34
0.24
0.590000
0.415515
0.274437
120
0.5819
0.43
0.30
0.872035
0.562119
0.356675
150
0.63
0.51
0.34
0.994252
0.713350
0.415515
180
0.70
0.57
0.38
1.203973
0.843970
0.478036
210
0.75
0.63
0.45
1.386294
0.994252
0.597837
240
0.81
0.65
0.50
1.660731
1.049822
0.693147
XME: concentration of methyl ester.
−ln (1-XME): triglyceride concentration.
3. Kinetic Model
The mole fraction for the transesterification reaction was established. For the reaction stoichiometry requires 3 Mol of methanol (M) and 1 mol of triglyceride (TG) to give 3 mol of methyl ester (ME) and 1 Mol of glycerol (GL). Transesterification reaction comprises three consecutive reversible reactions, where 1 mol of ME is produced in each step and monoglycerides (MG) and diglycerides (DG) are intermediate products [19]. The kinetic model used in this work is based on the following assumptions.
Firstly, keq should be considered not to depend on methanol concentration (due to its excess) (the reaction is considered pseudo-first order [20, 21]).
Production of intermediate species is negligible (the result reaction is a one step).
The chemical reaction occurred in the oil phase.
Based on assumption (1),
(1)-r=-d[TG]dt=k[TG][ROH]3.
And based on assumption (2),
(2)k′=k[ROH]3-r=-d[TG]dt=k′[TG]⟹lnTG0-lnTG=k′·t.
According to the mass balance,
(3)XME=1-[TG][TG0][TG]=[TG0][1-XME]dXMEdt=k′[1-XME]⟹-ln(1-XME)=k′·t.
Based on this model and the experimental data, at first we calculated the concentration of methyl ester at different times (based on the moles fraction). Second, the rate constants at each temperature were obtained. Third, the preexponential factors and activation energies are obtained by plotting the logarithm of the rate constants (k) versus 1/T of absolute temperature using the Arrhenius equation and in the final stage thermodynamic parameters were obtained such as ΔS and ΔH.
4. Results and Discussion4.1. Characterization of the Catalyst
The FT-IR spectra of H3PW12O40, CsXPW12O40 and CsXH3-XPW12O40/Fe-SiO2 are shown in Figure 1. The Keggin anion of HPW consists of a central phosphorous atom tetrahedrally coordinated by four oxygen atoms and surrounded by twelve octahedral WO6 units that share edges and corners in the structure.bands related to a Keggin structure (i(O–P–O) = 550 cm−1, indicative of the bending of the central oxygen of P–O–P; νas(W–Oe–W) at 798 cm−1, related to asymmetric stretching of tungsten with edge oxygen in W–O–W; νas(W–OcW) = 893 cm−1, related to the asymmetric stretching of corner oxygen in W–O–W; νas(W=O) = 983 cm−1, indicative of the asymmetric stretching of the terminal oxygen; and ν(P–O) at 1080 cm−1, assigned to asymmetric stretching of oxygen with a central phosphorous atom) [22, 23]. HPA salts maintain their corner Keggin structure with the addition of different amount of metallic Cs. When CsXH3-XPW12O40 is supported on Fe-SiO2, these bands have somewhat changed. The bands at 1080 and 890 cm−1 are overlapped by the characteristic band of SiO2, while these bands at 985 and 794 cm−1 shift to 966 and 805 cm−1, respectively.
FT-IR spectrums of the CsXH3-XPW12O40/Fe-SiO2 (a), H3PW12O40 (b), and CsXH3-XPW12O40 (c) nanocatalysts.
The XRD pattern of CsXH3-XPW12O40/Fe-SiO2 catalyst was presented in Figure 2. Supported HPA samples do not show diffraction patterns probably due to the following reasons: (i) after treatment at 500°C, the HPA practically loses its crystalline structure, (ii) HPA species were highly dispersed, and (iii) the deposited amount of HPA was not big enough to be detected by this technique [24].
XRD patterns of the Cs1H2PW12O40/Fe-SiO2 (a) and Cs1H2PW12O40 (b) nanocatalysts.
The TGA-DSC experiment on the catalyst precursor has shown four steps of mass loss (Figure 3). The first step at the temperature of 70–110°C was attributed to the evaporation of residual moistures in the catalyst precursor and loss of physisorbed waters. The second stage (190–280°C) is accompanying weight loss of the crystallization water expelling which is most likely the hydrated proton. The peak around 390–460°C is due to the decomposition of iron and Si oxalates to oxide phases. Most of weight loss happened from 580°C to 640°C due to the phase transition and formation of Fe2SiO4 (cubic). The weight loss curve is involved with a total overall weight loss of ca. 69 wt%. DSC measurement was performed in order to provide further evidence for the presence of the various species and evaluate their thermal behavior. As shown in Figure 3, the endothermic curve represents the removal of the physically adsorbed water from the material (70–110°C). Two exothermic peaks at around 190–280°C and 390–480°C are due to the crystallization water expelling which is most likely the hydrated proton and the decomposition of iron and Si oxalates to oxide phases, respectively. The exothermic peak at around 580–640°C is due to formation of iron silicate phase [25].
TGA and DSC curves for CsXH3-XPW12O40/Fe-SiO2 precursor.
Figure 4 shows SEM pictures of the precursor (a) and calcined catalyst (b). After calcination catalyst particle aggregated together and formed a spherical shape and more uniform particles which are beneficial to the activity and augmenting the surface of the catalyst that exhibited a large amount of aggregates than the precursor. The measured BET surface areas are 237.5 m2 cm3 g−1 for CsXH3-XPW12O40/Fe-SiO2 catalyst and corresponded pore volume is 0.7672 cm3 g−1 obtained from analysis of the desorption using the BJH (Barrett-Joyner-Halenda) method. The particle size could be calculated by Scherer-equation form XRD pattern (Figure 2). It is clear that the catalyst particle size was in nanodimension (45 nm). The CsXH3-XPW12O40/Fe-SiO2 catalyst was characterized with SEM (Figure 4). It is obvious in this figure that the crystal sizes were from 38–47 nm. This result confirmed the obtained results studied by using the Scherrer equation.
The SEM images of CsXH3-XPW12O40/Fe-SiO2 nanocatalyst, (a) precursor and (b) calcined catalyst.
4.2. Calculation of Rate Constant, Activation Energy, and Preexponential Factor
Plots of -ln(1-XME) versus (T) are given in Figure 5. Rate constant has been calculated using Figure 5. And activation energy is calculated through the Arrhenius equation:
(4)k=Ae-Ea/RT(5)⟹lnk=lnA-EaRT,
where (k) is the reaction constant, A is the frequency or preexponential factor, Ea is the activation energy of the reaction, R is the gas constant, and T is the absolute temperature.
Plots of -ln(1-X) versus time (min) at temperatures 60, 55, and 50°C for reaction of sunflower oil with methanol.
Therefore, plots of lnk versus 1/T are given in Figure 6. Activation energy (Ea) and preexponential factor have been calculated using Figure 6 to be 79.805 kJ/mol and 8.9 × 108 kJ/mol, respectively. Based on the proposed kinetic model, the kinetic parameters for this catalyst were determined. The experiments demonstrate that the reactions follow first-order kinetics. The proposed kinetic model describes the experimental results well and the rate constants follow the Arrhenius equation.
Arrhenius plot of lnk versus 1/T for reaction of sunflower oil with methanol.
4.3. Thermodynamic Parameters (ΔS and ΔH)
Based on the definition of Gibbs energy free using (7) and using linear plot of lnkeqversus ln1/T which is given in Figure 7 and by using (7) the respective values of ΔS and ΔH were calculated which are 0.0197 kJ/mol, 79.784 kJ/Kmol respectively:
(6)ΔG=-RTlnkeq(7)lnkeq=ΔSR-ΔHRT.
From the results (Table 2) of the thermodynamic parameters it can be found that the transesterification reaction is an endothermic reaction and, with increasing temperature, reaction rate increases. Moreover, enthalpy and entropy change are not affected by methanol concentration due to its excess [26].
Calculated values of thermodynamic parameters.
T (K)
ΔG (kJ/K·mol)
ΔH (kJ/K·mol)
ΔS (KJ/mol)
323
26.83272
328
25.98118
79.784
0.0197
333
25.19303
Plot of ΔG (J/K·mol) versus T (K) for reaction of sunflower oil with methanol.
5. Conclusions
The magnetic CsXH3-XPW12O40/Fe-SiO2 nanocatalyst was prepared for biodiesel production. Experimental conditions were varied as follows: reaction temperature 328–338 K, methanol/oil molar ratio = 12/1, and the reaction time 0–240 min. Thermodynamic properties such as ΔS and ΔH were successfully determined from equilibrium constants measured at different temperature. Activation energy (Ea) and preexponential factor have been calculated to be 79.805 kJ/mol and 8.9 × 108 kJ/mol, respectively. The experiments demonstrate that the reactions follow first-order kinetics. Also notably, recovery of the catalyst can be achieved easily with the help of an external magnet in a very short time (<20 seconds) with no need for expensive ultracentrifugation.
Acknowledgment
The authors are grateful to the Iran Nanotechnology Initiative Council (INIC) for its partial support of this project.
MathM. C.KumarS. P.ChettyS. V.Technologies for biodiesel production from used cooking oil: a review20101443393452-s2.0-7864963542610.1016/j.esd.2010.08.001BoeyP.-L.ManiamG. P.HamidS. A.Performance of calcium oxide as a heterogeneous catalyst in biodiesel production: a review2011168115222-s2.0-7995243995210.1016/j.cej.2011.01.009PrydeE. H.Vegetable oils as fuel alternatives: symposium overview19846110160916102-s2.0-5124918387710.1007/BF02541644DiasJ. M.Alvim-FerrazM. C. M.AlmeidaM. F.Comparison of the performance of different homogeneous alkali catalysts during transesterification of waste and virgin oils and evaluation of biodiesel quality20088717-18357235782-s2.0-5064912060010.1016/j.fuel.2008.06.014HelwaniZ.OthmanM. R.AzizN.KimJ.FernandoW. J. N.Solid heterogeneous catalysts for transesterification of triglycerides with methanol: a review20093631-21102-s2.0-6764938360110.1016/j.apcata.2009.05.021MeherL. C.Vidya SagarD.NaikS. N.Technical aspects of biodiesel production by transesterification: a review20061032482682-s2.0-2914450845410.1016/j.rser.2004.09.002QiuF.LiY.YangD.LiX.SunP.The use of heterogeneous catalysts to replace homogeneous ones can be expected to eliminate the problems associated with homogeneous catalysts201110241504156JacobsonK.GopinathR.MeherL. C.DalaiA. K.Solid acid catalyzed biodiesel production from waste cooking oil2008851-286912-s2.0-5654911678610.1016/j.apcatb.2008.07.005ZabetiM.Wan DaudW. M. A.ArouaM. K.Activity of solid catalysts for biodiesel production: a review20099067707772-s2.0-6564911296710.1016/j.fuproc.2009.03.010ShanmugamS.ViswanathanB.VaradarajanT. K.Esterification by solid acid catalysts: a comparison20042231-21431472-s2.0-744423966710.1016/j.molcata.2004.02.030LiuY.LoteroE.GoodwinJ. G.Jr.A comparison of the esterification of acetic acid with methanol using heterogeneous versus homogeneous acid catalysis200624222782862-s2.0-3374682772810.1016/j.jcat.2006.05.026HelwaniZ.OthmanM. R.AzizN.FernandoW. J. N.KimJ.Technologies for production of biodiesel focusing on green catalytic techniques: a review20099012150215142-s2.0-7164910761810.1016/j.fuproc.2009.07.016KimH.-J.ShulY.-G.HanH.Synthesis of heteropolyacid (H3PW12O40)/SiO2 nanoparticles and their catalytic properties20062991-246512-s2.0-2914447555210.1016/j.apcata.2005.10.010MizunoN.MisonoM.Pore structure and surface area of Cs (x = 0–3, M = W, Mo)198716967970OkuharaT.IchikiT.LeeK. Y.MisonoM.Dehydration mechanism of ethanol in the pseudoliquid phase of H3-XCsXPW12O40198955293301PizzioL. R.VázquezP. G.CáceresC. V.BlancoM. N.AlessoE. N.TorvisoM. R.LantañoB.MoltrasioG. Y.AguirreJ. M.C-alkylation reactions catalyzed by silica-supported Keggin heteropolyacids20052871182-s2.0-2034436859110.1016/j.apcata.2004.11.023KimH.-J.ShulY.-G.HanH.Synthesis of heteropolyacid (H3PW12O40)/SiO2 nanoparticles and their catalytic properties20062991-246512-s2.0-2914447555210.1016/j.apcata.2005.10.010ZhaoX.HanY.SunX.WangY.Structure and catalytic performance of H3PW12O40/SiO2 prepared by several methods200728191962-s2.0-3424867763510.1016/S1872-2067(07)60011-4ShuQ.GaoJ.LiaoY.WangJ.Reaction kinetics of biodiesel synthesis from waste oil using a carbon-based solid acid catalyst20111911631682-s2.0-7995252222110.1016/S1004-9541(09)60193-2SakaS.KusdianaD.Biodiesel fuel from rapeseed oil as prepared in supercritical methanol1998802225231KusdianaD.SakaS.Kinetics of transesterification in rapeseed oil to biodiesel fuel as treated in supercritical methanol20018056936982-s2.0-003586493410.1016/S0016-2361(00)00140-XSrilathaK.SreeR.Prabhavathi DeviB. L. A.PrasadP. S. S.PrasadR. B. N.LingaiahN.Preparation of biodiesel from rice bran fatty acids catalyzed by heterogeneous cesium-exchanged 12-tungstophosphoric acids20121165557NakkaL.MolinariJ. E.WachsI. E.Surface and bulk aspects of mixed oxide catalytic nanoparticles: oxidation and dehydration of CH3OH by polyoxometallates20091314215544155542-s2.0-7035031508710.1021/ja904957dYoriJ. C.GrauJ. M.BenítezV. M.SepúlvedaJ.Hydroisomerization-cracking of n-octane on heteropolyacid H3PW12O40 supported on ZrO2, SiO2 and carbon Effect of Pt incorporation on catalyst performance2005286171782-s2.0-1884438499510.1016/j.apcata.2005.03.003DengQ.ZhouW.LiX.PengZ.JiangS.YueM.CaiT.Microwave radiation solid-phase synthesis of phosphotungstate nanoparticle catalysts and photocatalytic degradation of formaldehyde20072621-21491552-s2.0-3384614480010.1016/j.molcata.2006.08.043FernandesS. A.CardosoA. L.Da SilvaM. J.A novel kinetic study of H3PW12O40: catalyzed oleic acid esterification with methanol via HNMR spectroscopy201296981032-s2.0-8485580095710.1016/j.fuproc.2011.12.025