# Methyl Esters Selectivity of Transesterification Reaction with Homogenous Alkaline Catalyst to Produce Biodiesel in Batch, Plug Flow, and Continuous Stirred Tank Reactors

Authors: N. F. Nasir; W. R. W. Daud; S. K. Kamarudin; Z. Yaakob
Journal: International Journal of Chemical Engineering
 (2014)

---

## Abstract

Selectivity concept is essential in establishing the best operating conditions for attaining maximum production of the desired product. For complex reaction such as biodiesel fuel synthesis, kinetic studies of transesterification reaction have revealed the mechanism of the reaction and rate constants. The objectives of this research are to develop the kinetic parameters for determination of methyl esters and glycerol selectivity, evaluate the significance of the reverse reaction in transesterification reaction, and examine the influence of reaction characteristics (reaction temperature, methanol to oil molar ratio, and the amount of catalyst) on selectivity. For this study, published reaction rate constants of transesterification reaction were used to develop mathematical expressions for selectivities. In order to examine the base case and reversible transesterification, two calculation schemes (Case  1 and Case  2) were established. An enhanced selectivity was found in the base case of transesterification reaction. The selectivity was greatly improved at optimum reaction temperature (60°C), molar ratio (9 : 1), catalyst concentration (1.5 wt.%), and low free fatty acid feedstock. Further research might explore the application of selectivity for specifying reactor configurations.


---

## Body

## 1. Introduction

Biodiesel, a major renewable fuel to replace fossil fuel, is clean, environmental friendly, and derived from readily available sources such as vegetable oils and animal fats. In the biodiesel production, transesterification reaction is used to convert triglycerides to esters. It can be defined as alcoholysis of oil and fat with an alcohol to form esters and glycerol [1]. Various types of catalysts such as acidic and basic catalysts, which are either homogeneous or heterogeneous, as well as enzymes, can catalyse the reaction.Transesterification of the oil often represented by triglycerides,T, with methanol, A, in the presence of a catalyst, usually alkaline, yields esters of fatty acids, E, and glycerol, G. Monoglyceride, M, and diglyceride, D, are the intermediates. The overall reaction is written typically as
(1)T+3A⟺catalyst3E+GHowever, (1) is misleading because the rates of the forward reactions are not equal. The transesterification reaction is known to occur in three stages, where D produced from methanolysis in the first stage reaction becomes the reactant for the second stage methanolysis and the M produced from the second stage reaction becomes the reactant for the third stage reaction. The third stage reaction produces G as a byproduct. All three reactions produce different types of methyl esters, collectively designated by E. The reactions are shown in (2) to (4), where k+1–k-3 are rate constants:
(2)T+A⟺k-1k+1E1+D(3)D+A⟺k-2k+2E2+M(4)M+A⟺k-3k+3E3+GA stoichiometric balance of the above three reversible reactions yield the correct form of the overall reaction, which is given below:(5)T+(1+(r+2-r-2)(r+1-r-1)+(r+3-r-3)(r+1-r-1))A→catalystE1+(r+2-r-2)(r+1-r-1)E2+(r+3-r-3)(r+1-r-1)E3kkkkkk+((r+3-r-3)(r+1-r-1))G,
where r+1 is the rate of the forward reaction in reaction 1, r-1 is the rate of the backward reaction in reaction 1, r+2 is the rate of the forward reaction in reaction 2, r-2 is the rate of the backward reaction in reaction 2, r+3 is the rate of the forward reaction in reaction 3, and r-3 is the rate of the backward reaction in reaction 3. Equation (5) reduces to (1) if and only if all the net rates of reactions are equal or (r+1-r-1)=(r+2-r-2)=(r+3-r-3).Kinetic studies of the chemical reactions in biodiesel production provide parameters that are used to predict the extent of the reactions at any time under any condition. Published kinetic data of biodiesel production are mostly dependent on the type of oil feed stock and the type of catalyst used as well as on the reaction conditions such as reaction temperature and catalyst concentration [2–4]. Previous evidence suggests that an important factor in the transesterification process is the degree of mixing between the alcohol, A, and the triglyceride, T, phases. As T and A phases are not miscible and they form two liquid layers upon their initial introduction into the reactor, mechanical mixing is typically applied to increase the contact between the reactants, resulting in an increase in mass transfer rate [5].Many works have studied the kinetics of the transesterification reaction mechanism and determined its reaction rate constants [2–11]. Leevijit et al. [3] who conducted batch transesterification reaction and evaluated the reaction rate constants k+1–k-3 found that the forward reactions are faster than the reverse reactions. The rate constants are presented as dataset 10 in Table 1. Similar kinetic study of palm oil methanolysis was conducted by Narvàez et al. [9] in order to investigate the effects of temperatures and concentration of NaOH catalyst. His results also suggest that the reaction kinetics for the forward reactions dominate over the reverse. From the result of their study on overall reaction of transesterification, the energy of activation for the forward reaction (16.4 kcal gmol−1) is greater than the energy of activation for the reverse reaction (13.6 kcal gmol−1).Table 1
Reaction rate constants data for transesterification reaction in biodiesel production.
Sr. number
1
2
3
4
5
6
7
8
9
10
k
+
1 (L/mol·s)
0.0500
0.0920
0.1650
0.0060
0.0063
0.0021
0.0031
0.1672
0.2251
0.634
k
-
1 (L/mol·s)
0.1100
0.1750
0.2710
0.0016
0.0007
0.0035
0.0000
0.9373
0.4711
0
k
+
2 (L/mol/s)
0.2150
0.5440
1.3050
0.0442
0.0092
0.0072
0.0123
0.3553
0.1317
7.1040
k
-
2 (L/mol/s)
1.2280
2.4350
4.6410
0.0127
0.0023
0.0102
0.1224
0.2674
0.3917
4.9120
k
+
3 (L/mol/s)
0.2420
0.3270
0.4340
0.0050
0.0070
0.0231
0.0092
1.9608
0.1190
7.8600
k
-
3 (L/mol/s)
0.0070
0.0110
0.0170
0.0003
0.0027
0.0035
0.0000
0.0077
0.0181
0.1210
Temperature (°C)
50
60
70
65
65
30
30
60
60
60
Molar ratioMR
6 : 1
6 : 1
6 : 1
9 : 1
3 : 1
9 : 1
9 : 1
6 : 1
6 : 1
6 : 1
Catalyst concentration (wt.%)
0.2
0.2
0.2
1.5
1.5
1.5
0.5
1.0
1.0
1.0
Oil feed stock
Soybean
Soybean
Soybean
Castor
Castor
Castor
Castor
Rapeseed
WasteSunflower
Palm
Reference
[5]
[5]
[5]
[28]
[28]
[28]
[28]
[27]
[27]
[3]A kinetic study on the transesterification ofJatropha curcas oil by Jain and Sharma [4] using a two-step acid-base catalyst process in a batch system found that the esterification reaction was slower than the transesterification reaction because the process was assisted by the esterification reaction. Apart from that, recent studies reported the kinetics of transesterification reaction in the absence of a catalyst using soybean oil, palm oil, and rapeseeds oil [6–8, 10]. The rate constants, conversion, and yield of E showed an increasing trend with the reaction temperature, but the E content in the reaction product decreased as the reaction temperature was increased.Reactor design uses information, knowledge, and experience from a variety of areas such as thermodynamics, chemical kinetics, fluid mechanics, heat transfer, mass transfer, and economics [12]. A considerable research has been devoted into the determination of the reaction rate constants of transesterification reactions of various types of edible and inedible vegetable oils. The published reaction rate constants at various conditions can be further used to determine the selectivity of E and G, providing opportunities for specifying reactor configurations using process system engineering tools such as the attainable region method. Selectivity can be described in a number of different ways, such as the production of the desirable component divided by the amount of limiting reactant converted or the production of desired component divided by the production of the undesired component [13]. The objectives of this present work are to develop the kinetic models of transesterification reaction and determine the selectivity of E and G using published reaction kinetic constants of the transesterification reactions and to investigate the effect of reaction conditions on them. In addition, this paper evaluates the significance of the reverse reactions in transesterification reaction on the E and G selectivities.


## 2. Methodology

### 2.1. Selectivity of Methyl Esters and Glycerol in Plug Flow Reactor

The amount of free fatty acids (FFAs) in oil feedstock could vary from 3 to 40% depending on types of oil used [14]. The saponification may occur in alkali-catalyzed transesterification reaction due to the reaction of acids with catalyst to produce soaps [15]. The transesterification reaction with the resulting soap formation tends to become viscous or forms a gel, and it is difficult to separate the glycerol from the mixture [16]. The saponification reaction is undesirable because it binds the catalyst into a form that does not contribute to accelerating the reaction [14]. For feedstock having a high free fatty acid content, pretreatment to the oil can be performed with methanol [17] or glycerol [18, 19] to decrease its acidity and thus reducing the potential for saponification to occur during the reaction. The presence of water in catalyzed transesterification reaction is also undesirable because water can consume the catalyst and reduce catalyst efficiency [20]. The saponification reaction can be taken into account while developing the kinetic model of alkali-catalyzed transesterification reaction [21]. However, the kinetics parameter for alkali-catalyzed transesterification in the present work is developed by assuming that triglycerides were pretreated earlier. Also, this work represents a preliminary design of reactors for a biodiesel plant. Therefore, some important issues such as the existence of competing reaction along with the transesterification such as saponification reaction and the presence of water in the oil, the transport process, heterogeneity of the reaction, and loss of polarity of an alkaline catalyst were not considered.The kinetic model was developed by first performing stoichiometric mole balance of each species in a plug flow reactor (PFR) using the reaction kinetics in (2) to (4). However, since a simpler mole balance of each species in a batch reactor produces the same solution as in Section 3, the batch reactor balance is used instead.Mole balance of T is as follows:(6)-dCTdt=k+1CTCA-k-1CDCE1.Mole balance of D is as follows:(7)-dCDdt=-(k+1CTCA-k-1CDCE1-k+2CDCA+k-2CMCE2).Mole balance of M is as follows:(8)-dCMdt=-(k+2CDCA-k-2CMCE2-k+3CMCA+k-3CGCE3).Mole balance of E1 is as follows:(9)-dCE1dt=-(k+1CTCA-k-1CDCE1).Mole balance of E2 is as follows:(10)-dCE2dt=-(k+2CDCA-k-2CMCE2).Mole balance of E3 is as follows:(11)-dCE3dt=-(k+3CMCA-k-3CGCE3).Mole balance of A (methanol) is as follows:(12)-dCAdt=k+1CTCA-k-1CDCE1+k+2CDCA-k-2CMCE2+k+3CMCA-k-3CGCE3.Mole balance of G is as follows:(13)-dCGdt=-(k+3CMCA-k-3CGCE3).Published reaction rate constants(k+1–kx-3) are listed in Table 1. In this work, methanol is assumed in excess while triglyceride is the limiting reactant.Closed form simultaneous solution of ordinary differential equations (6) to (13) by analysis is possible if only the forward reactions are considered significant but the backward reactions are neglected. On the other hand, if both forward and reverse reactions are equally significant, closed form solution for (6) to (13) cannot be found and they have to be solved numerically using a 4th order Runge-Kutta algorithm, which is available in MATLAB. The solution generates the molar concentrations of all species, Ci, as functions of the design variables, that is, the conversion of the limiting reactant, XT, and the molar ratio of the excess variable, MR. Conversion is defined as
(14)conversion=lossofreactantfeedofreactant.For reactantT, the conversion is
(15)XT=NT0-NTNT0=1-NT0/VNT0/V=1-CTCT0.Selectivity is defined as the ratio of the desired product to the amount of limiting reactant that has undergone chemical change [22]. That is,
(16)Selectivity=amountofdesiredproductamountoflimitingreactantthathasundergonechemicalchange.In this study, the total methyl ester is defined as(17)Totalmethylester=methylester1+methylester2+methylester3.The selectivity of the total methyl ester,SE, and of glycerol, SG, defined as the moles of total methyl ester or glycerol produced, respectively, per moles of the limiting reactant T that have reacted is then formulated as functions of the conversion of triglyceride, XT, as shown in
(18)SE=CE1+CE2+CE3CToXT,SG=CGCToXT.


### 2.2. Selectivity of Methyl Esters and Glycerol in Continuous Stirred Tank Reactor

The kinetic model for continuous stirred tank reactor (CSTR) is much simpler. The stoichiometric mole balance of each species in a CSTR, using the reaction kinetics in (2) to (4), yields a set of nonlinear equations in (19) to (26).Mole balance ofT is as follows:
(19)CTo-CT=VFT(k+1CTCA-k-1CDCE1)=τ(k+1CTCA-k-1CDCE1).Mole balance ofD is as follows:
(20)CDo-CD=-τ(k+1CTCA-k-1CDCE1-k+2CDCA+k-2CMCE2).Mole balance ofM is as follows:
(21)CMo-CM=-τ(k+2CDCA-k-2CMCE2-k+3CMCA+k-3CGCE3).Mole balance ofE1 is as follows:
(22)CE1o-CE1=-τ(k+1CTCA-k-1CDCE1).Mole balance ofE2 is as follows:
(23)CE2o-CE2=-τ(k+2CDCA-k-2CMCE2).Mole balance ofE3 is as follows:
(24)CE3o-CE3=-τ(k+3CMCA-k-3CGCE3).Mole balance ofA (methanol) is as follows:
(25)CAo-CA=τ(k+1CTCA-k-1CDCE1+k+2CDCA-k-2CMCE2kkkkkk+k+3CMCA-k-3CGCE3).Mole balance ofG is as follows:
(26)CGo-CG=-τ(k+3CMCA-k-3CGCE3).Closed form simultaneous solution of nonlinear equations (19) to (26) by analysis is not possible. They have to be solved numerically using a matrix inversion algorithm, which is available in MATLAB. The solution generates the molar concentrations of all species, Ci, as functions of the design variables, that is, the conversion of the limiting reactant, and the molar ratio of the excess variable, A.


## 2.1. Selectivity of Methyl Esters and Glycerol in Plug Flow Reactor

The amount of free fatty acids (FFAs) in oil feedstock could vary from 3 to 40% depending on types of oil used [14]. The saponification may occur in alkali-catalyzed transesterification reaction due to the reaction of acids with catalyst to produce soaps [15]. The transesterification reaction with the resulting soap formation tends to become viscous or forms a gel, and it is difficult to separate the glycerol from the mixture [16]. The saponification reaction is undesirable because it binds the catalyst into a form that does not contribute to accelerating the reaction [14]. For feedstock having a high free fatty acid content, pretreatment to the oil can be performed with methanol [17] or glycerol [18, 19] to decrease its acidity and thus reducing the potential for saponification to occur during the reaction. The presence of water in catalyzed transesterification reaction is also undesirable because water can consume the catalyst and reduce catalyst efficiency [20]. The saponification reaction can be taken into account while developing the kinetic model of alkali-catalyzed transesterification reaction [21]. However, the kinetics parameter for alkali-catalyzed transesterification in the present work is developed by assuming that triglycerides were pretreated earlier. Also, this work represents a preliminary design of reactors for a biodiesel plant. Therefore, some important issues such as the existence of competing reaction along with the transesterification such as saponification reaction and the presence of water in the oil, the transport process, heterogeneity of the reaction, and loss of polarity of an alkaline catalyst were not considered.The kinetic model was developed by first performing stoichiometric mole balance of each species in a plug flow reactor (PFR) using the reaction kinetics in (2) to (4). However, since a simpler mole balance of each species in a batch reactor produces the same solution as in Section 3, the batch reactor balance is used instead.Mole balance of T is as follows:(6)-dCTdt=k+1CTCA-k-1CDCE1.Mole balance of D is as follows:(7)-dCDdt=-(k+1CTCA-k-1CDCE1-k+2CDCA+k-2CMCE2).Mole balance of M is as follows:(8)-dCMdt=-(k+2CDCA-k-2CMCE2-k+3CMCA+k-3CGCE3).Mole balance of E1 is as follows:(9)-dCE1dt=-(k+1CTCA-k-1CDCE1).Mole balance of E2 is as follows:(10)-dCE2dt=-(k+2CDCA-k-2CMCE2).Mole balance of E3 is as follows:(11)-dCE3dt=-(k+3CMCA-k-3CGCE3).Mole balance of A (methanol) is as follows:(12)-dCAdt=k+1CTCA-k-1CDCE1+k+2CDCA-k-2CMCE2+k+3CMCA-k-3CGCE3.Mole balance of G is as follows:(13)-dCGdt=-(k+3CMCA-k-3CGCE3).Published reaction rate constants(k+1–kx-3) are listed in Table 1. In this work, methanol is assumed in excess while triglyceride is the limiting reactant.Closed form simultaneous solution of ordinary differential equations (6) to (13) by analysis is possible if only the forward reactions are considered significant but the backward reactions are neglected. On the other hand, if both forward and reverse reactions are equally significant, closed form solution for (6) to (13) cannot be found and they have to be solved numerically using a 4th order Runge-Kutta algorithm, which is available in MATLAB. The solution generates the molar concentrations of all species, Ci, as functions of the design variables, that is, the conversion of the limiting reactant, XT, and the molar ratio of the excess variable, MR. Conversion is defined as
(14)conversion=lossofreactantfeedofreactant.For reactantT, the conversion is
(15)XT=NT0-NTNT0=1-NT0/VNT0/V=1-CTCT0.Selectivity is defined as the ratio of the desired product to the amount of limiting reactant that has undergone chemical change [22]. That is,
(16)Selectivity=amountofdesiredproductamountoflimitingreactantthathasundergonechemicalchange.In this study, the total methyl ester is defined as(17)Totalmethylester=methylester1+methylester2+methylester3.The selectivity of the total methyl ester,SE, and of glycerol, SG, defined as the moles of total methyl ester or glycerol produced, respectively, per moles of the limiting reactant T that have reacted is then formulated as functions of the conversion of triglyceride, XT, as shown in
(18)SE=CE1+CE2+CE3CToXT,SG=CGCToXT.


## 2.2. Selectivity of Methyl Esters and Glycerol in Continuous Stirred Tank Reactor

The kinetic model for continuous stirred tank reactor (CSTR) is much simpler. The stoichiometric mole balance of each species in a CSTR, using the reaction kinetics in (2) to (4), yields a set of nonlinear equations in (19) to (26).Mole balance ofT is as follows:
(19)CTo-CT=VFT(k+1CTCA-k-1CDCE1)=τ(k+1CTCA-k-1CDCE1).Mole balance ofD is as follows:
(20)CDo-CD=-τ(k+1CTCA-k-1CDCE1-k+2CDCA+k-2CMCE2).Mole balance ofM is as follows:
(21)CMo-CM=-τ(k+2CDCA-k-2CMCE2-k+3CMCA+k-3CGCE3).Mole balance ofE1 is as follows:
(22)CE1o-CE1=-τ(k+1CTCA-k-1CDCE1).Mole balance ofE2 is as follows:
(23)CE2o-CE2=-τ(k+2CDCA-k-2CMCE2).Mole balance ofE3 is as follows:
(24)CE3o-CE3=-τ(k+3CMCA-k-3CGCE3).Mole balance ofA (methanol) is as follows:
(25)CAo-CA=τ(k+1CTCA-k-1CDCE1+k+2CDCA-k-2CMCE2kkkkkk+k+3CMCA-k-3CGCE3).Mole balance ofG is as follows:
(26)CGo-CG=-τ(k+3CMCA-k-3CGCE3).Closed form simultaneous solution of nonlinear equations (19) to (26) by analysis is not possible. They have to be solved numerically using a matrix inversion algorithm, which is available in MATLAB. The solution generates the molar concentrations of all species, Ci, as functions of the design variables, that is, the conversion of the limiting reactant, and the molar ratio of the excess variable, A.


## 3. Solution of the PFR and CSTR

### 3.1. Removal of Temporal Dependence

Equations (6) to (13) cannot be solved in their time dependent forms because of the single degree of freedom of the problem. In order to solve the equations simultaneously, the degree of freedom is reduced to zero by changing the independent variable from time to the conversion or concentration of the limiting reactant using the chain rule and the mole balance equation for T (6). Consider the following:
(27)dCidCT=(dCidt)(dtdCT).Hence, (7) to (13) can be rewritten with concentration of T as their independent variable as follows:
(28)dCDdCT=-(1-α2CDCA-α3CMCE2CTCA-α1CDCE1),dCMdCT=-(α2CDCA-α3CMCE2-α4CMCA+α5CGCE3CTCA-α1CDCE1),dCE1dCT=-1,dCE2dCT=-(α2CDCA-α3CMCE2CTCA-α1CDCE1),dCE3dCT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),dCAdCT=1+α2CDCA-α3CMCE2+α4CMCA-α5CGCE3CTCA-α1CDCE1,dCGdCT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),
where
(29)α1=k-1k+1,α2=k+2k+1,α3=k-2k+1,α4=k+3k+1,α5=k-3k+1.Similarly, (28) cannot be solved in their time dependent forms because of the single degree of freedom of the problem. In order to solve the equations simultaneously, the degree of freedom is reduced to zero by changing the independent variable from time to the conversion or concentration of the limiting reactant using algebraic manipulation and the mole balance equation for T (19). Consider the following:
(30)Cio-CiCTo-CT=-(Cio-Ciτ)(τCTo-CT).Hence, (28) can be rewritten with concentration of T as their independent variable as follows:
(31)CDo-CDCTo-CT=-(1-α2CDCA-α3CMCE2CTCA-α1CDCE1),CMo-CMCTo-CT=-(α2CDCA-α3CMCE2-α4CMCA+α5CGCE3CTCA-α1CDCE1),CE1o-CE1CTo-CT=-1,CE2o-CE2CTo-CT=-(α2CDCA-α3CMCE2CTCA-α1CDCE1),CE3o-CE3CTo-CT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),CAo-CACTo-CT=1+α2CDCA-α3CMCE2+α4CMCA-α5CGCE3CTCA-α1CDCE1,CGo-CGCTo-CT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1).Two main cases are considered: the base case where the reactions (2) to (4) are irreversible and the alternative case where the reactions are reversible.


### 3.2. Case  1: Base Case: Irreversible Reactions

In Case  1, the base case is considered by assuming (2) to (4) to be irreversible reactions. In other words, the forward reactions are so dominant that the reverse reaction can be neglected and the overall effect is that the reactions are irreversible. The results of the base case can then be compared with that of the alternative Case  2, where the reactions are all reversible, in order to verify the importance of the reverse reaction and its implication to the maximum attainable molar concentration of methyl ester. Since the reverse reactions are omitted, the value of α1, α3, and α5 is zero. Hence, (28) for PFR can be rewritten as follows:
(32)dCDdCT=-(1-α2CDCT),(33)dCMdCT=-(α2CD-α4CMCT),(34)dCE1dCT=-1,(35)dCE2dCT=-α2CDCT,(36)dCE3dCT=-α4CMCT,(37)dCAdCT=1+α2CD+α4CMCT,(38)dCGdCT=-α4CMCT.The solution of the set of ordinary differential equations (32) can be found analytically by using the integrating factor method [23]. The concentrations of all the species generated by the solution are then expressed in terms of the triglycerides conversion, XT, as shown in
(39)CT=CTo(1-XT),CD=CToα2-1[(CTCTo)-(CTCTo)α2]=CToα2-1[(1-XT)-(1-XT)α2],CM=CTo[α21-α2[(1-XT)1-α4-(1-XT)α2α2-α4]kkkkkk+α2(1-XT)α4(1-α4)(α2-α4)],CE1=CTo(1-(CTCTo))=CTo(1-(1-XT))=CToXT,CE2=CTo{1-α2α2-1[(1-XT)-(1-XT)α2α2]},CE3=α4CTo[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4],CA=CTo{((1-α2)α2+2α4-1(1-α4))MR+α2(1-α2)kkkkkk×[((1-α2)α2+2α4-1(1-α4))XTkkkkkkkkkk-(2α4-α2α2(α2-α4))[1-(1-XT)α2]kkkkkkkkkk-(1-α2)(1-α4)(α2-α4)kkkkkkkkk×[1-(1-XT)α4]((1-α2)α2+2α4-1(1-α4))]},CG=α4CTo[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4].For the CSTR, (31) can be manipulated algebraically and rewritten with conversion of T as their independent variable as follows:
(40)CD=CToXT(1-XT)(1-(1-α2)XT),(41)CM=α2CToXT2(1-XT)(1-(1-α4)XT)(1-(1-α2)XT),(42)CE1=CToXT,(43)CE2=α2CToXT2(1-(1-α2)XT),(44)CE3=α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT),(45)CA=CTo(MR-1+α2CToXT2(1-(1-α2)XT)kkkkkk+α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT)),(46)CG=α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT).


### 3.3. Case  2: Alternative Case: Reversible Reactions

In Case  2, an alternative case was considered by assuming that the reactions in (2) to (4) are reversible. Together with the base case, where the reactions are irreversible, it can be used to discover the significance of the reverse and forward reactions.Since the values of allαi’s in (28) for the PFR and (31) for the CSTR are not zero, no closed form solution to them can be found because they are too complex for analytical solution. The solution of both sets of equations can be achieved by numerical methods using a computer simulation software. Equations (28) for the PFR are solved by using the ordinary differential equation 45, (ODE 45) in MATLAB. Concentration of each species is obtained using Rungge-Kutta numerical integration algorithm. On the other hand, (31) for the CSTR are solved by using the simultaneous nonlinear equations solver in MATLAB. Concentration of each species is obtained using the multivariable Newton-Raphson algorithm.


### 3.4. Selectivity of Methyl Esters and Glycerol for Case  1 and Case  2

The net concentration of methyl esters is the sum of concentration of ester 1, ester 2, and ester 3. In Case  1, the net concentration of methyl esters in the PFR is given by(47)CE1+CE2+CE3=CTo{+α2[1-(1-XT)α4](1-α4)(α2-α4)α4]1+XT-α2α2-1[(1-XT)-(1-XT)α2α2]kkkkkkkkk+α4[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4]}.Similarly, the net concentration of methyl esters in the CSTR for Case  1 is given by(48)CE1+CE2+CE3=CE1o+CE2o+CE3o+[CToXT+α2CToXT2(1-(1-α2)XT)kkkkkk+α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT)].The selectivity of methyl esters in Case  1 for PFR is given by(49)SE={+α2[1-(1-XT)α4](1-α4)(α2-α4)α4XT]XT+1XT-α2α2-1[(1-XT)XT-(1-XT)α2α2XT]kk+α4[α21-α2[11-α4-[1-(1-XT)α2](α2-α4)α2XT]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4XT]}.The selectivity of glycerol in Case  1 for PFR is given by(50)SG=α4XT[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkik+α2[1-(1-XT)α4](1-α4)(α2-α4)α4].On the other hand, the selectivity of methyl esters in Case  1 for CSTR is given by(51)SE=[α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT)1+α2XT(1-(1-α2)XT)kkk+α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT)].The selectivity of glycerol in Case  1 for CSTR is given by(52)SG=α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT).For Case  2, there is no closed form expression for selectivity of methyl esters and glycerol. The selectivities of methyl esters and glycerol have to be calculated using (18) from the concentration of the species generated by the numerical solution.


## 3.1. Removal of Temporal Dependence

Equations (6) to (13) cannot be solved in their time dependent forms because of the single degree of freedom of the problem. In order to solve the equations simultaneously, the degree of freedom is reduced to zero by changing the independent variable from time to the conversion or concentration of the limiting reactant using the chain rule and the mole balance equation for T (6). Consider the following:
(27)dCidCT=(dCidt)(dtdCT).Hence, (7) to (13) can be rewritten with concentration of T as their independent variable as follows:
(28)dCDdCT=-(1-α2CDCA-α3CMCE2CTCA-α1CDCE1),dCMdCT=-(α2CDCA-α3CMCE2-α4CMCA+α5CGCE3CTCA-α1CDCE1),dCE1dCT=-1,dCE2dCT=-(α2CDCA-α3CMCE2CTCA-α1CDCE1),dCE3dCT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),dCAdCT=1+α2CDCA-α3CMCE2+α4CMCA-α5CGCE3CTCA-α1CDCE1,dCGdCT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),
where
(29)α1=k-1k+1,α2=k+2k+1,α3=k-2k+1,α4=k+3k+1,α5=k-3k+1.Similarly, (28) cannot be solved in their time dependent forms because of the single degree of freedom of the problem. In order to solve the equations simultaneously, the degree of freedom is reduced to zero by changing the independent variable from time to the conversion or concentration of the limiting reactant using algebraic manipulation and the mole balance equation for T (19). Consider the following:
(30)Cio-CiCTo-CT=-(Cio-Ciτ)(τCTo-CT).Hence, (28) can be rewritten with concentration of T as their independent variable as follows:
(31)CDo-CDCTo-CT=-(1-α2CDCA-α3CMCE2CTCA-α1CDCE1),CMo-CMCTo-CT=-(α2CDCA-α3CMCE2-α4CMCA+α5CGCE3CTCA-α1CDCE1),CE1o-CE1CTo-CT=-1,CE2o-CE2CTo-CT=-(α2CDCA-α3CMCE2CTCA-α1CDCE1),CE3o-CE3CTo-CT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1),CAo-CACTo-CT=1+α2CDCA-α3CMCE2+α4CMCA-α5CGCE3CTCA-α1CDCE1,CGo-CGCTo-CT=-(α4CMCA-α5CGCE3CTCA-α1CDCE1).Two main cases are considered: the base case where the reactions (2) to (4) are irreversible and the alternative case where the reactions are reversible.


## 3.2. Case  1: Base Case: Irreversible Reactions

In Case  1, the base case is considered by assuming (2) to (4) to be irreversible reactions. In other words, the forward reactions are so dominant that the reverse reaction can be neglected and the overall effect is that the reactions are irreversible. The results of the base case can then be compared with that of the alternative Case  2, where the reactions are all reversible, in order to verify the importance of the reverse reaction and its implication to the maximum attainable molar concentration of methyl ester. Since the reverse reactions are omitted, the value of α1, α3, and α5 is zero. Hence, (28) for PFR can be rewritten as follows:
(32)dCDdCT=-(1-α2CDCT),(33)dCMdCT=-(α2CD-α4CMCT),(34)dCE1dCT=-1,(35)dCE2dCT=-α2CDCT,(36)dCE3dCT=-α4CMCT,(37)dCAdCT=1+α2CD+α4CMCT,(38)dCGdCT=-α4CMCT.The solution of the set of ordinary differential equations (32) can be found analytically by using the integrating factor method [23]. The concentrations of all the species generated by the solution are then expressed in terms of the triglycerides conversion, XT, as shown in
(39)CT=CTo(1-XT),CD=CToα2-1[(CTCTo)-(CTCTo)α2]=CToα2-1[(1-XT)-(1-XT)α2],CM=CTo[α21-α2[(1-XT)1-α4-(1-XT)α2α2-α4]kkkkkk+α2(1-XT)α4(1-α4)(α2-α4)],CE1=CTo(1-(CTCTo))=CTo(1-(1-XT))=CToXT,CE2=CTo{1-α2α2-1[(1-XT)-(1-XT)α2α2]},CE3=α4CTo[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4],CA=CTo{((1-α2)α2+2α4-1(1-α4))MR+α2(1-α2)kkkkkk×[((1-α2)α2+2α4-1(1-α4))XTkkkkkkkkkk-(2α4-α2α2(α2-α4))[1-(1-XT)α2]kkkkkkkkkk-(1-α2)(1-α4)(α2-α4)kkkkkkkkk×[1-(1-XT)α4]((1-α2)α2+2α4-1(1-α4))]},CG=α4CTo[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4].For the CSTR, (31) can be manipulated algebraically and rewritten with conversion of T as their independent variable as follows:
(40)CD=CToXT(1-XT)(1-(1-α2)XT),(41)CM=α2CToXT2(1-XT)(1-(1-α4)XT)(1-(1-α2)XT),(42)CE1=CToXT,(43)CE2=α2CToXT2(1-(1-α2)XT),(44)CE3=α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT),(45)CA=CTo(MR-1+α2CToXT2(1-(1-α2)XT)kkkkkk+α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT)),(46)CG=α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT).


## 3.3. Case  2: Alternative Case: Reversible Reactions

In Case  2, an alternative case was considered by assuming that the reactions in (2) to (4) are reversible. Together with the base case, where the reactions are irreversible, it can be used to discover the significance of the reverse and forward reactions.Since the values of allαi’s in (28) for the PFR and (31) for the CSTR are not zero, no closed form solution to them can be found because they are too complex for analytical solution. The solution of both sets of equations can be achieved by numerical methods using a computer simulation software. Equations (28) for the PFR are solved by using the ordinary differential equation 45, (ODE 45) in MATLAB. Concentration of each species is obtained using Rungge-Kutta numerical integration algorithm. On the other hand, (31) for the CSTR are solved by using the simultaneous nonlinear equations solver in MATLAB. Concentration of each species is obtained using the multivariable Newton-Raphson algorithm.


## 3.4. Selectivity of Methyl Esters and Glycerol for Case  1 and Case  2

The net concentration of methyl esters is the sum of concentration of ester 1, ester 2, and ester 3. In Case  1, the net concentration of methyl esters in the PFR is given by(47)CE1+CE2+CE3=CTo{+α2[1-(1-XT)α4](1-α4)(α2-α4)α4]1+XT-α2α2-1[(1-XT)-(1-XT)α2α2]kkkkkkkkk+α4[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkkkkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4]}.Similarly, the net concentration of methyl esters in the CSTR for Case  1 is given by(48)CE1+CE2+CE3=CE1o+CE2o+CE3o+[CToXT+α2CToXT2(1-(1-α2)XT)kkkkkk+α2α4CToXT3(1-(1-α4)XT)(1-(1-α2)XT)].The selectivity of methyl esters in Case  1 for PFR is given by(49)SE={+α2[1-(1-XT)α4](1-α4)(α2-α4)α4XT]XT+1XT-α2α2-1[(1-XT)XT-(1-XT)α2α2XT]kk+α4[α21-α2[11-α4-[1-(1-XT)α2](α2-α4)α2XT]kkkkkkkk+α2[1-(1-XT)α4](1-α4)(α2-α4)α4XT]}.The selectivity of glycerol in Case  1 for PFR is given by(50)SG=α4XT[α21-α2[XT1-α4-[1-(1-XT)α2](α2-α4)α2]kkkkik+α2[1-(1-XT)α4](1-α4)(α2-α4)α4].On the other hand, the selectivity of methyl esters in Case  1 for CSTR is given by(51)SE=[α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT)1+α2XT(1-(1-α2)XT)kkk+α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT)].The selectivity of glycerol in Case  1 for CSTR is given by(52)SG=α2α4XT2(1-(1-α4)XT)(1-(1-α2)XT).For Case  2, there is no closed form expression for selectivity of methyl esters and glycerol. The selectivities of methyl esters and glycerol have to be calculated using (18) from the concentration of the species generated by the numerical solution.


## 4. Results and Discussion

### 4.1. Comparison of Selectivities for Cases  1 and  2

This study has shown the development of kinetic parameters for determining methyl esters and glycerol selectivity in a biodiesel transesterification reaction; it has also established several mathematical equations for the same. Furthermore, a particularly important parameter for designing a reactor, that is, product selectivity, was determined for each case and condition. For comparing the base case and an alternative case, datasets 2 and 3 were evaluated using each calculation schemes. The difference in the selectivity between Case  1 and Case  2, which were investigated under the same conditions, is very significant and is attributed to the reversibility factor of the reaction, as illustrated in Figures1 and 2. The results of the base case show that when all the reactions are irreversible, it implies that the maximum attainable molar concentration and selectivity of methyl ester have been reached. This verifies the significant role of the forward reactions in (2)–(4).Figure 1
Selectivity plots of Cases  1 and  2 at reaction temperature 60°C, molar ratio 6 : 1, and catalyst concentration 0.2 wt.%.Figure 2
Selectivity plots of Cases  1 and  2 at reaction temperature 70°C, molar ratio 6 : 1, and catalyst concentration 0.2 wt.%.For a simple reversible reaction, we can see the effects of thermodynamic parameter like temperature on the reversibility of the reaction whereas the reversibility of a reaction is controlled by reaction equilibrium constant. The reaction equilibrium constant,K, is defined as the forward rate constant divided by the reverse rate constant and affected by temperature [12]. This transesterification reaction is a complex reaction which consists of three reversible and consecutive reactions. Therefore, for simplification and preliminary design purpose, selectivity is effectively used to determine the effect of temperature and role of catalyst on the reversibility. Figures 1 and 2 show significant changes of the selectivity plots for conversion of triglycerides from 0.7 to 0.9 and the plots finally reach the equilibrium at the end of the reaction. In order to access the role of temperature and amount of catalyst on the reversibility of transesterification reaction, Table 2 was constructed by observing the methyl esters selectivity for conversion of triglycerides from 0.7 to 0.9. It presents the difference in selectivity for Case  1 and Case  2. Datasets (1–3) are used to determine the effects of different temperatures on the reversibility of the transesterification reaction which represent a temperature of 50, 60, and 70, respectively. Conversion of triglyceride at 0.9 is selected for analysis. From this data, we can see that the forward reaction alone is able to provide better selectivity than reverse reaction at any temperature. Datasets (6-7) are used to analyze the role catalyst on the reversibility of transesterification reaction whereas datasets 6 and 7 represent the amount of catalyst of 1.5% and 0.5%, respectively. A significant role of the forward reaction once again was determined at a higher catalyst concentration. The results obtained also indicate that a higher catalyst concentration controls the forward reactions and gives preferably product selectivity.Table 2
Selectivity at 0.7–0.9 triglycerides conversion for datasets (1–3) and (6-7).
Conversion,XT
Dataset
Methyl esters selectivity
Case  1
Case  2
0.7
1
2.610
2.526
2
2.643
2.511
3
2.630
2.470
6
2.610
2.613
7
2.493
2.313
0.8
1
2.766
2.675
2
2.784
2.666
3
2.769
2.629
6
2.764
2.766
7
2.682
2.480
0.9
1
2.895
2.832
2
2.903
2.816
3
2.892
2.788
6
2.893
2.879
7
2.852
2.655According to Narvàez et al. [9] and Noureddini and Zhu [5], the forward reactions in a transesterification reaction are the most important reactions. However, Diasakou et al. [7] claimed that the reverse reactions during a transesterification reaction do not influence the forward reactions. In this study, a lower concentration of methyl ester was found for Case  2, which affects the selectivity and methyl ester conversion, which is likely due to the reverse nature of the reaction [4]. This result suggests that for establishing the optimal chemical process for biodiesel fuel production, it can be assumed that when irreversible reactions are employed, better selectivity can be obtained because of the absence of reverse reactions.


### 4.2. Effects of Reaction Characteristics on Selectivity

The reactor design stage for a biodiesel plant utilizes methyl ester selectivity for the material balances and stream costs. Attaining the appropriate operating conditions through product selectivity can eliminate material and energy recycling, thus reducing costs. The effects of the reaction temperature, molar ratio of methanol to oil, and the catalyst concentration on the selectivity of methyl esters and glycerol are analyzed for both Case  1 and Case  2 in the following subsections.


#### 4.2.1. Effect of Reaction Temperature on Selectivity

An investigation on the effect of temperature on methyl ester selectivity was conducted using datasets (1–3), where the molar ratio of alcohol to oil and the catalyst concentration are fixed. The Arrhenius equation was used for this investigation:(53)k=k°exp⁡⁡(-ERT).According to (53), the reaction rate, k, is an integral part of the selectivity expressions. It depends on the energy of activation (E) obtained from the Arrhenius equation [22]. For reactions at different temperatures, fixed molar ratios, and catalyst concentrations, the activation energy is constant. In this study, the kinetic data used to analyze the temperature effect were adapted from Noureddini and Zhu [5]. Because the data published in their research utilize the average reaction rate constants at 50°C and constant activation energy, (53) can be manipulated to obtain the other reaction rate constants (at 60°C and 70°C) listed in Table 1. The results of the activation energy investigated in their study are shown in Table 3.Table 3
Activation energy.
Reaction
Activation energy (cal/mol)
T
→
D
13145
D
→
T
9932
D
→
M
19860
M
→
D
14639
M
→
G
6421
G
→
M
9588Figures3 and 4 show the effects of the reaction temperature on the selectivity of methyl esters and glycerol for Cases  1 and  2, respectively. The reaction temperature had a small but substantial effect on the selectivity. The figures show that the selectivity of methyl esters is high at a high temperature in the early conversion of T. However, in both figures, the selectivities of E for both Cases  1 and 2 at 70°C were found to decrease toward the completion of the reaction. As the boiling point of methanol is 64.7°C, an additional 5°C has the effect of lowering the product concentration and affecting the selectivity of E at the equilibrium. This temperature effect changes the slope of selectivity from a concave to a convex shape, which is also related with the reaction rate constant at different temperatures.Figure 3
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  1.Figure 4
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  2.According to Darnoko and Cheryan [2], the transesterification reaction rates were increased at 65°C as compared with those at 50°C. The results of their research also showed an enhanced methyl ester yield at 65°C compared with that at 50°C, which was because of the higher viscosity of oil at 50°C. In addition, 55°C is a stable temperature that can maintain the oil liquidity. This result indicates that an optimum temperature promotes the formation of methyl esters and, thus, contributes to a higher selectivity of biodiesel products. It can also be suggested that the reaction rate constant, k, which is greatly dependent on the temperature and selectivity of E, is greatly improved at the optimum temperature.The findings confirm that methyl ester selectivity effects can be realistically obtained at the optimum temperature, which is 60°C in this study. As the reactor optimization requires the relationship between the conversion and selectivity in both cases (Case  1 and Case  2), the selectivity within the studied temperature indicates that the biodiesel reactor may work more efficiently in converting T to E at 60°C.


#### 4.2.2. Effect of Molar Ratio of Methanol to Oil on Selectivity

An investigation into the effects of the methanol-to-oil molar ratio involves datasets 4 and 5. Selectivity plots of methyl esters and glycerol are shown in Figures5 and 6. For Cases  1 and  2, the selectivity of methyl esters and glycerol was relatively affected by the molar ratio of methanol to oil. It can be seen from the figures that the higher molar ratio of methanol to oil encourages the formation of transesterification reaction products, thus contributing to a higher selectivity.Figure 5
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  1.Figure 6
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  2.The most striking result to emerge from the figures is the change in the shape of the selectivity plot from concave to convex for a molar ratio of 9 : 1. An increasing selectivity of methyl ester can be found in an early conversion of triglyceride, but toward equilibrium, the selectivity is found to slightly decrease. A possible explanation for this might be that the methyl ester concentration has both quadratic as well as cubic effects that can be explained by (42)–(44) and Figure 7.Figure 7
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 4.According to (42), the selectivity of methyl ester 1 is constant at 1. Equation (43) is a quadratic equation, which gives a parabolic shape to the plot of methyl ester 2. Equation (44) clearly illustrates the selectivity of methyl ester 3 when the equation is a cubic one. In more detail, the power of two is more dominant at an XT of less than 0.3, which means that the formation of E2 is more active under this condition. However, we can see an active formation of E3 at an XT of 0.5–1 when the power of three is dominant. At the same time, E1 is formed linearly with a triglyceride conversion. The stepwise transesterification reaction in (2)–(4) shows that each E is formed consecutively. This means that to achieve the maximum yield, the maximum selectivity at the maximum T conversion is necessary, and each E must be formed during the process. Therefore, it can be assumed that the selectivity of methyl ester is more greatly affected by the quadratic and cubic functions at a higher molar ratio of methanol to oil than at a lower one.In addition, Freedman et al. [11] reported that the best methyl ester conversion is achieved at a molar ratio of 6 : 1 with alkaline catalyst and that higher molar ratios do not increase the productivity but do increase the cost of alcohol recovery. At a high molar excess of methanol, a pseudo-first-order reaction appeared more likely than a second-order reaction. However, high molar ratio of methanol is used for noncatalytic and supercritical methanol reaction [6, 11, 20, 24]. Narvàez et al. [9] stated that when an excess amount of methanol is used, it can positively affect the reversible characteristic of the reaction and the displacement of the reaction toward the product. As a consequence, excess methanol does not necessarily contribute to an increase in the concentration of methyl ester and glycerol.


#### 4.2.3. Effect of Catalyst Concentration on Selectivity

Apart from the temperature and molar ratio, the parameters that affect the selectivity are the same as those influencing the reaction rates, namely, the catalyst. The effects of the catalyst concentration are studied through datasets 6 and 7, and the selectivity plots can be seen in Figures8 and 9. When analyzing the figures, it can be seen that the selectivity of both products is higher at a higher catalyst concentration.Figure 8
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  1.Figure 9
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  2.Toward the equilibrium of the reaction, the selectivity was observed to increase at a higher catalyst concentration, thus contributing to an increased amount of methyl ester formation. Vicente et al. [25] stated that the largest methyl ester conversions are obtained when a large catalyst concentration is employed (1.3%) at a mild temperature (20–50°C). This agrees with the resulting selectivity curves plotted, which indicates that the catalyst concentration has a positive effect on the selectivity of methyl esters and glycerol. However, a greater application of an alkali catalyst can lead to the production of large amounts of soap, and extra cost is required to remove the catalyst from the reaction system. Srivastava and Prasad [26] have stated that a catalyst concentration of 0.5–1% can yield up to a 99% methyl ester conversion.In Figures8 and 9, there is a clear increasing trend of the selectivity of methyl ester for a 1.5 wt.% catalyst concentration during the early conversion of the reactant, which decreases at the end of the reaction. The same transition trend has been found when analyzing the influence of the molar ratio on methyl ester selectivity. A further detailed analysis showed that the quadratic and cubic effects are more significant at a higher catalyst concentration than at a lower catalyst concentration. As shown in Figure 10, the selectivity of E2 has a quadratic plot, whereas the selectivity of E3 has a parabolic shape. Referring to (42)–(44), we can see that E1 is linearly formed with a T conversion. E2 is produced more at a T conversion of lower than 0.35, where we can see that the power of two in (43) is more dominant at this stage. At a T conversion of over 0.5, the effect of the power of three in (44) is dominant, causing a higher formation of E3. Because the transesterification reaction is a reversible and consecutive equation, the formation of three Es is important to achieve a better selectivity. With an optimum amount of catalyst, the reactor becomes more efficient in converting the reactant into the products.Figure 10
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 6.


### 4.3. Comparison of Methyl Ester Selectivities from Different Oil Feedstocks

To compare the selectivities of methyl esters and glycerol from different feedstocks, datasets (8–10) were used. Figure11 presents the selectivity of methyl esters derived from rapeseed oil, waste sunflower oil, and palm oil for Case  1 and Case  2. As can be seen from the graph, the methyl esters and glycerol were more likely to form if palm oil is used as the feedstock compared to the other two oil feedstocks. The present finding also suggests that the selectivity of methyl ester was much lower when using low cost feedstock such as waste sunflower oil. Waste oil contains large amounts of free fatty acids (FFAs) (about 20%) which affect the yield of fuel from that process. The amount of FFAs in palm oil is 0.64% [3], but for rapeseed oil, the author used fresh oil [27] which is estimated less than 1%.Figure 11
Selectivity plots of methyl ester from various feedstock for Cases  1 and  2.As discussed before, these free fatty acids will react with alkali catalysts to produce soaps that will inhibit the transesterification reaction. Soaps will allow emulsification that causes the phase separation of methyl ester and glycerol to be less sharp. It produces water that can hydrolyze the triglycerides and contribute to the formation of more soap [14]. The effect of reversibility of transesterification can be seen significantly when analysing the selectivity of methyl ester derived from palm oil. In Figure 11, methyl ester derived from palm oil shows a higher selectivity plot for base case. This finding suggests that forward reaction dominates the reaction and contributes to a better selectivity when low FFAs feedstock is used in biodiesel production. As there are many low cost feedstocks available for biodiesel production, a pretreatment to these feedstocks can be carried out in order to ensure the value of FFAs is within the allowable range.


## 4.1. Comparison of Selectivities for Cases  1 and  2

This study has shown the development of kinetic parameters for determining methyl esters and glycerol selectivity in a biodiesel transesterification reaction; it has also established several mathematical equations for the same. Furthermore, a particularly important parameter for designing a reactor, that is, product selectivity, was determined for each case and condition. For comparing the base case and an alternative case, datasets 2 and 3 were evaluated using each calculation schemes. The difference in the selectivity between Case  1 and Case  2, which were investigated under the same conditions, is very significant and is attributed to the reversibility factor of the reaction, as illustrated in Figures1 and 2. The results of the base case show that when all the reactions are irreversible, it implies that the maximum attainable molar concentration and selectivity of methyl ester have been reached. This verifies the significant role of the forward reactions in (2)–(4).Figure 1
Selectivity plots of Cases  1 and  2 at reaction temperature 60°C, molar ratio 6 : 1, and catalyst concentration 0.2 wt.%.Figure 2
Selectivity plots of Cases  1 and  2 at reaction temperature 70°C, molar ratio 6 : 1, and catalyst concentration 0.2 wt.%.For a simple reversible reaction, we can see the effects of thermodynamic parameter like temperature on the reversibility of the reaction whereas the reversibility of a reaction is controlled by reaction equilibrium constant. The reaction equilibrium constant,K, is defined as the forward rate constant divided by the reverse rate constant and affected by temperature [12]. This transesterification reaction is a complex reaction which consists of three reversible and consecutive reactions. Therefore, for simplification and preliminary design purpose, selectivity is effectively used to determine the effect of temperature and role of catalyst on the reversibility. Figures 1 and 2 show significant changes of the selectivity plots for conversion of triglycerides from 0.7 to 0.9 and the plots finally reach the equilibrium at the end of the reaction. In order to access the role of temperature and amount of catalyst on the reversibility of transesterification reaction, Table 2 was constructed by observing the methyl esters selectivity for conversion of triglycerides from 0.7 to 0.9. It presents the difference in selectivity for Case  1 and Case  2. Datasets (1–3) are used to determine the effects of different temperatures on the reversibility of the transesterification reaction which represent a temperature of 50, 60, and 70, respectively. Conversion of triglyceride at 0.9 is selected for analysis. From this data, we can see that the forward reaction alone is able to provide better selectivity than reverse reaction at any temperature. Datasets (6-7) are used to analyze the role catalyst on the reversibility of transesterification reaction whereas datasets 6 and 7 represent the amount of catalyst of 1.5% and 0.5%, respectively. A significant role of the forward reaction once again was determined at a higher catalyst concentration. The results obtained also indicate that a higher catalyst concentration controls the forward reactions and gives preferably product selectivity.Table 2
Selectivity at 0.7–0.9 triglycerides conversion for datasets (1–3) and (6-7).
Conversion,XT
Dataset
Methyl esters selectivity
Case  1
Case  2
0.7
1
2.610
2.526
2
2.643
2.511
3
2.630
2.470
6
2.610
2.613
7
2.493
2.313
0.8
1
2.766
2.675
2
2.784
2.666
3
2.769
2.629
6
2.764
2.766
7
2.682
2.480
0.9
1
2.895
2.832
2
2.903
2.816
3
2.892
2.788
6
2.893
2.879
7
2.852
2.655According to Narvàez et al. [9] and Noureddini and Zhu [5], the forward reactions in a transesterification reaction are the most important reactions. However, Diasakou et al. [7] claimed that the reverse reactions during a transesterification reaction do not influence the forward reactions. In this study, a lower concentration of methyl ester was found for Case  2, which affects the selectivity and methyl ester conversion, which is likely due to the reverse nature of the reaction [4]. This result suggests that for establishing the optimal chemical process for biodiesel fuel production, it can be assumed that when irreversible reactions are employed, better selectivity can be obtained because of the absence of reverse reactions.


## 4.2. Effects of Reaction Characteristics on Selectivity

The reactor design stage for a biodiesel plant utilizes methyl ester selectivity for the material balances and stream costs. Attaining the appropriate operating conditions through product selectivity can eliminate material and energy recycling, thus reducing costs. The effects of the reaction temperature, molar ratio of methanol to oil, and the catalyst concentration on the selectivity of methyl esters and glycerol are analyzed for both Case  1 and Case  2 in the following subsections.


### 4.2.1. Effect of Reaction Temperature on Selectivity

An investigation on the effect of temperature on methyl ester selectivity was conducted using datasets (1–3), where the molar ratio of alcohol to oil and the catalyst concentration are fixed. The Arrhenius equation was used for this investigation:(53)k=k°exp⁡⁡(-ERT).According to (53), the reaction rate, k, is an integral part of the selectivity expressions. It depends on the energy of activation (E) obtained from the Arrhenius equation [22]. For reactions at different temperatures, fixed molar ratios, and catalyst concentrations, the activation energy is constant. In this study, the kinetic data used to analyze the temperature effect were adapted from Noureddini and Zhu [5]. Because the data published in their research utilize the average reaction rate constants at 50°C and constant activation energy, (53) can be manipulated to obtain the other reaction rate constants (at 60°C and 70°C) listed in Table 1. The results of the activation energy investigated in their study are shown in Table 3.Table 3
Activation energy.
Reaction
Activation energy (cal/mol)
T
→
D
13145
D
→
T
9932
D
→
M
19860
M
→
D
14639
M
→
G
6421
G
→
M
9588Figures3 and 4 show the effects of the reaction temperature on the selectivity of methyl esters and glycerol for Cases  1 and  2, respectively. The reaction temperature had a small but substantial effect on the selectivity. The figures show that the selectivity of methyl esters is high at a high temperature in the early conversion of T. However, in both figures, the selectivities of E for both Cases  1 and 2 at 70°C were found to decrease toward the completion of the reaction. As the boiling point of methanol is 64.7°C, an additional 5°C has the effect of lowering the product concentration and affecting the selectivity of E at the equilibrium. This temperature effect changes the slope of selectivity from a concave to a convex shape, which is also related with the reaction rate constant at different temperatures.Figure 3
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  1.Figure 4
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  2.According to Darnoko and Cheryan [2], the transesterification reaction rates were increased at 65°C as compared with those at 50°C. The results of their research also showed an enhanced methyl ester yield at 65°C compared with that at 50°C, which was because of the higher viscosity of oil at 50°C. In addition, 55°C is a stable temperature that can maintain the oil liquidity. This result indicates that an optimum temperature promotes the formation of methyl esters and, thus, contributes to a higher selectivity of biodiesel products. It can also be suggested that the reaction rate constant, k, which is greatly dependent on the temperature and selectivity of E, is greatly improved at the optimum temperature.The findings confirm that methyl ester selectivity effects can be realistically obtained at the optimum temperature, which is 60°C in this study. As the reactor optimization requires the relationship between the conversion and selectivity in both cases (Case  1 and Case  2), the selectivity within the studied temperature indicates that the biodiesel reactor may work more efficiently in converting T to E at 60°C.


### 4.2.2. Effect of Molar Ratio of Methanol to Oil on Selectivity

An investigation into the effects of the methanol-to-oil molar ratio involves datasets 4 and 5. Selectivity plots of methyl esters and glycerol are shown in Figures5 and 6. For Cases  1 and  2, the selectivity of methyl esters and glycerol was relatively affected by the molar ratio of methanol to oil. It can be seen from the figures that the higher molar ratio of methanol to oil encourages the formation of transesterification reaction products, thus contributing to a higher selectivity.Figure 5
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  1.Figure 6
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  2.The most striking result to emerge from the figures is the change in the shape of the selectivity plot from concave to convex for a molar ratio of 9 : 1. An increasing selectivity of methyl ester can be found in an early conversion of triglyceride, but toward equilibrium, the selectivity is found to slightly decrease. A possible explanation for this might be that the methyl ester concentration has both quadratic as well as cubic effects that can be explained by (42)–(44) and Figure 7.Figure 7
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 4.According to (42), the selectivity of methyl ester 1 is constant at 1. Equation (43) is a quadratic equation, which gives a parabolic shape to the plot of methyl ester 2. Equation (44) clearly illustrates the selectivity of methyl ester 3 when the equation is a cubic one. In more detail, the power of two is more dominant at an XT of less than 0.3, which means that the formation of E2 is more active under this condition. However, we can see an active formation of E3 at an XT of 0.5–1 when the power of three is dominant. At the same time, E1 is formed linearly with a triglyceride conversion. The stepwise transesterification reaction in (2)–(4) shows that each E is formed consecutively. This means that to achieve the maximum yield, the maximum selectivity at the maximum T conversion is necessary, and each E must be formed during the process. Therefore, it can be assumed that the selectivity of methyl ester is more greatly affected by the quadratic and cubic functions at a higher molar ratio of methanol to oil than at a lower one.In addition, Freedman et al. [11] reported that the best methyl ester conversion is achieved at a molar ratio of 6 : 1 with alkaline catalyst and that higher molar ratios do not increase the productivity but do increase the cost of alcohol recovery. At a high molar excess of methanol, a pseudo-first-order reaction appeared more likely than a second-order reaction. However, high molar ratio of methanol is used for noncatalytic and supercritical methanol reaction [6, 11, 20, 24]. Narvàez et al. [9] stated that when an excess amount of methanol is used, it can positively affect the reversible characteristic of the reaction and the displacement of the reaction toward the product. As a consequence, excess methanol does not necessarily contribute to an increase in the concentration of methyl ester and glycerol.


### 4.2.3. Effect of Catalyst Concentration on Selectivity

Apart from the temperature and molar ratio, the parameters that affect the selectivity are the same as those influencing the reaction rates, namely, the catalyst. The effects of the catalyst concentration are studied through datasets 6 and 7, and the selectivity plots can be seen in Figures8 and 9. When analyzing the figures, it can be seen that the selectivity of both products is higher at a higher catalyst concentration.Figure 8
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  1.Figure 9
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  2.Toward the equilibrium of the reaction, the selectivity was observed to increase at a higher catalyst concentration, thus contributing to an increased amount of methyl ester formation. Vicente et al. [25] stated that the largest methyl ester conversions are obtained when a large catalyst concentration is employed (1.3%) at a mild temperature (20–50°C). This agrees with the resulting selectivity curves plotted, which indicates that the catalyst concentration has a positive effect on the selectivity of methyl esters and glycerol. However, a greater application of an alkali catalyst can lead to the production of large amounts of soap, and extra cost is required to remove the catalyst from the reaction system. Srivastava and Prasad [26] have stated that a catalyst concentration of 0.5–1% can yield up to a 99% methyl ester conversion.In Figures8 and 9, there is a clear increasing trend of the selectivity of methyl ester for a 1.5 wt.% catalyst concentration during the early conversion of the reactant, which decreases at the end of the reaction. The same transition trend has been found when analyzing the influence of the molar ratio on methyl ester selectivity. A further detailed analysis showed that the quadratic and cubic effects are more significant at a higher catalyst concentration than at a lower catalyst concentration. As shown in Figure 10, the selectivity of E2 has a quadratic plot, whereas the selectivity of E3 has a parabolic shape. Referring to (42)–(44), we can see that E1 is linearly formed with a T conversion. E2 is produced more at a T conversion of lower than 0.35, where we can see that the power of two in (43) is more dominant at this stage. At a T conversion of over 0.5, the effect of the power of three in (44) is dominant, causing a higher formation of E3. Because the transesterification reaction is a reversible and consecutive equation, the formation of three Es is important to achieve a better selectivity. With an optimum amount of catalyst, the reactor becomes more efficient in converting the reactant into the products.Figure 10
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 6.


## 4.2.1. Effect of Reaction Temperature on Selectivity

An investigation on the effect of temperature on methyl ester selectivity was conducted using datasets (1–3), where the molar ratio of alcohol to oil and the catalyst concentration are fixed. The Arrhenius equation was used for this investigation:(53)k=k°exp⁡⁡(-ERT).According to (53), the reaction rate, k, is an integral part of the selectivity expressions. It depends on the energy of activation (E) obtained from the Arrhenius equation [22]. For reactions at different temperatures, fixed molar ratios, and catalyst concentrations, the activation energy is constant. In this study, the kinetic data used to analyze the temperature effect were adapted from Noureddini and Zhu [5]. Because the data published in their research utilize the average reaction rate constants at 50°C and constant activation energy, (53) can be manipulated to obtain the other reaction rate constants (at 60°C and 70°C) listed in Table 1. The results of the activation energy investigated in their study are shown in Table 3.Table 3
Activation energy.
Reaction
Activation energy (cal/mol)
T
→
D
13145
D
→
T
9932
D
→
M
19860
M
→
D
14639
M
→
G
6421
G
→
M
9588Figures3 and 4 show the effects of the reaction temperature on the selectivity of methyl esters and glycerol for Cases  1 and  2, respectively. The reaction temperature had a small but substantial effect on the selectivity. The figures show that the selectivity of methyl esters is high at a high temperature in the early conversion of T. However, in both figures, the selectivities of E for both Cases  1 and 2 at 70°C were found to decrease toward the completion of the reaction. As the boiling point of methanol is 64.7°C, an additional 5°C has the effect of lowering the product concentration and affecting the selectivity of E at the equilibrium. This temperature effect changes the slope of selectivity from a concave to a convex shape, which is also related with the reaction rate constant at different temperatures.Figure 3
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  1.Figure 4
Selectivity of methyl ester at various temperatures and fixed molar ratio and catalyst concentration for Case  2.According to Darnoko and Cheryan [2], the transesterification reaction rates were increased at 65°C as compared with those at 50°C. The results of their research also showed an enhanced methyl ester yield at 65°C compared with that at 50°C, which was because of the higher viscosity of oil at 50°C. In addition, 55°C is a stable temperature that can maintain the oil liquidity. This result indicates that an optimum temperature promotes the formation of methyl esters and, thus, contributes to a higher selectivity of biodiesel products. It can also be suggested that the reaction rate constant, k, which is greatly dependent on the temperature and selectivity of E, is greatly improved at the optimum temperature.The findings confirm that methyl ester selectivity effects can be realistically obtained at the optimum temperature, which is 60°C in this study. As the reactor optimization requires the relationship between the conversion and selectivity in both cases (Case  1 and Case  2), the selectivity within the studied temperature indicates that the biodiesel reactor may work more efficiently in converting T to E at 60°C.


## 4.2.2. Effect of Molar Ratio of Methanol to Oil on Selectivity

An investigation into the effects of the methanol-to-oil molar ratio involves datasets 4 and 5. Selectivity plots of methyl esters and glycerol are shown in Figures5 and 6. For Cases  1 and  2, the selectivity of methyl esters and glycerol was relatively affected by the molar ratio of methanol to oil. It can be seen from the figures that the higher molar ratio of methanol to oil encourages the formation of transesterification reaction products, thus contributing to a higher selectivity.Figure 5
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  1.Figure 6
Selectivity of methyl ester at various molar ratios, constant temperature of 65°C, and catalyst concentration of 1.5 wt.% for Case  2.The most striking result to emerge from the figures is the change in the shape of the selectivity plot from concave to convex for a molar ratio of 9 : 1. An increasing selectivity of methyl ester can be found in an early conversion of triglyceride, but toward equilibrium, the selectivity is found to slightly decrease. A possible explanation for this might be that the methyl ester concentration has both quadratic as well as cubic effects that can be explained by (42)–(44) and Figure 7.Figure 7
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 4.According to (42), the selectivity of methyl ester 1 is constant at 1. Equation (43) is a quadratic equation, which gives a parabolic shape to the plot of methyl ester 2. Equation (44) clearly illustrates the selectivity of methyl ester 3 when the equation is a cubic one. In more detail, the power of two is more dominant at an XT of less than 0.3, which means that the formation of E2 is more active under this condition. However, we can see an active formation of E3 at an XT of 0.5–1 when the power of three is dominant. At the same time, E1 is formed linearly with a triglyceride conversion. The stepwise transesterification reaction in (2)–(4) shows that each E is formed consecutively. This means that to achieve the maximum yield, the maximum selectivity at the maximum T conversion is necessary, and each E must be formed during the process. Therefore, it can be assumed that the selectivity of methyl ester is more greatly affected by the quadratic and cubic functions at a higher molar ratio of methanol to oil than at a lower one.In addition, Freedman et al. [11] reported that the best methyl ester conversion is achieved at a molar ratio of 6 : 1 with alkaline catalyst and that higher molar ratios do not increase the productivity but do increase the cost of alcohol recovery. At a high molar excess of methanol, a pseudo-first-order reaction appeared more likely than a second-order reaction. However, high molar ratio of methanol is used for noncatalytic and supercritical methanol reaction [6, 11, 20, 24]. Narvàez et al. [9] stated that when an excess amount of methanol is used, it can positively affect the reversible characteristic of the reaction and the displacement of the reaction toward the product. As a consequence, excess methanol does not necessarily contribute to an increase in the concentration of methyl ester and glycerol.


## 4.2.3. Effect of Catalyst Concentration on Selectivity

Apart from the temperature and molar ratio, the parameters that affect the selectivity are the same as those influencing the reaction rates, namely, the catalyst. The effects of the catalyst concentration are studied through datasets 6 and 7, and the selectivity plots can be seen in Figures8 and 9. When analyzing the figures, it can be seen that the selectivity of both products is higher at a higher catalyst concentration.Figure 8
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  1.Figure 9
Selectivity of methyl ester at various catalyst concentrations, constant temperature of 30°C, and molar ratio 9 : 1 for Case  2.Toward the equilibrium of the reaction, the selectivity was observed to increase at a higher catalyst concentration, thus contributing to an increased amount of methyl ester formation. Vicente et al. [25] stated that the largest methyl ester conversions are obtained when a large catalyst concentration is employed (1.3%) at a mild temperature (20–50°C). This agrees with the resulting selectivity curves plotted, which indicates that the catalyst concentration has a positive effect on the selectivity of methyl esters and glycerol. However, a greater application of an alkali catalyst can lead to the production of large amounts of soap, and extra cost is required to remove the catalyst from the reaction system. Srivastava and Prasad [26] have stated that a catalyst concentration of 0.5–1% can yield up to a 99% methyl ester conversion.In Figures8 and 9, there is a clear increasing trend of the selectivity of methyl ester for a 1.5 wt.% catalyst concentration during the early conversion of the reactant, which decreases at the end of the reaction. The same transition trend has been found when analyzing the influence of the molar ratio on methyl ester selectivity. A further detailed analysis showed that the quadratic and cubic effects are more significant at a higher catalyst concentration than at a lower catalyst concentration. As shown in Figure 10, the selectivity of E2 has a quadratic plot, whereas the selectivity of E3 has a parabolic shape. Referring to (42)–(44), we can see that E1 is linearly formed with a T conversion. E2 is produced more at a T conversion of lower than 0.35, where we can see that the power of two in (43) is more dominant at this stage. At a T conversion of over 0.5, the effect of the power of three in (44) is dominant, causing a higher formation of E3. Because the transesterification reaction is a reversible and consecutive equation, the formation of three Es is important to achieve a better selectivity. With an optimum amount of catalyst, the reactor becomes more efficient in converting the reactant into the products.Figure 10
Selectivity plots: (a) methyl ester 1, (b) methyl ester 2, and (c) methyl ester 3 for dataset 6.


## 4.3. Comparison of Methyl Ester Selectivities from Different Oil Feedstocks

To compare the selectivities of methyl esters and glycerol from different feedstocks, datasets (8–10) were used. Figure11 presents the selectivity of methyl esters derived from rapeseed oil, waste sunflower oil, and palm oil for Case  1 and Case  2. As can be seen from the graph, the methyl esters and glycerol were more likely to form if palm oil is used as the feedstock compared to the other two oil feedstocks. The present finding also suggests that the selectivity of methyl ester was much lower when using low cost feedstock such as waste sunflower oil. Waste oil contains large amounts of free fatty acids (FFAs) (about 20%) which affect the yield of fuel from that process. The amount of FFAs in palm oil is 0.64% [3], but for rapeseed oil, the author used fresh oil [27] which is estimated less than 1%.Figure 11
Selectivity plots of methyl ester from various feedstock for Cases  1 and  2.As discussed before, these free fatty acids will react with alkali catalysts to produce soaps that will inhibit the transesterification reaction. Soaps will allow emulsification that causes the phase separation of methyl ester and glycerol to be less sharp. It produces water that can hydrolyze the triglycerides and contribute to the formation of more soap [14]. The effect of reversibility of transesterification can be seen significantly when analysing the selectivity of methyl ester derived from palm oil. In Figure 11, methyl ester derived from palm oil shows a higher selectivity plot for base case. This finding suggests that forward reaction dominates the reaction and contributes to a better selectivity when low FFAs feedstock is used in biodiesel production. As there are many low cost feedstocks available for biodiesel production, a pretreatment to these feedstocks can be carried out in order to ensure the value of FFAs is within the allowable range.


## 5. Conclusions

This paper has presented a comprehensive kinetic modeling of a transesterification reaction for biodiesel with homogeneous alkaline catalyst. The purpose of the study was to develop the kinetic models of transesterification reaction, determine the selectivity of methyl ester and glycerol, and evaluate the significance of the reverse reactions in the transesterification reaction. This study has shown that a process model solution was found for two types of common biodiesel reactors. The concentration of each species during transesterification was found, leading to a selectivity determination during the biodiesel production. In addition, the reversible and irreversible effects of a transesterification reaction were discovered. The forward reaction is more dominant than the reverse reaction, and the maximum selectivity can be found if a base case of transesterification reaction is assumed to have occurred in the biodiesel reactor system. For the various reaction characteristics under investigation, the selectivities are influenced by the temperature, molar ratio of methanol to oil, catalyst concentration, and type of oil used as the feedstock. Furthermore, an analysis of the selectivity plots indicates an enhanced effect of the reaction characteristics on the selectivity. This research has shown that the selectivity of methyl ester is maximized at 60°C, a 9 : 1 molar ratio of methanol to oil, and a catalyst concentration of 1.5%. This investigation into the effects of the molar ratio of methanol to oil and the catalyst concentration on the selectivity can be made more interesting if several more curves are plotted at a molar ratio of methanol to oil between 3 : 1 to 9 : 1 and a catalyst concentration between 0.5 and 1.5%. However, a limitation of this research remains in the availability of the reaction rate constants under various reaction conditions. In addition, a comparison should be made through mutual reaction conditions. The findings from this study will enhance the concept of selectivity and its applications in the evaluation of complex transesterification reactions when the goal is to attain the maximum production of methyl esters and glycerol. It is recommended that further research is undertaken in exploring the selectivity data for specifying the reactor configurations and may be used in the synthesis of reactor network for biodiesel production.



---

*Source: 931264-2014-08-25.xml*