Effect of Neutralization on Yield of Corn Oil Methyl Ester Compared to Other Vegetable Oils
This research started by producing corn oil methyl ester using homogeneous base catalyst (NaOH pellets) in lab-scale. However, the recovered yield of the corn methyl ester was very low (< 50%) because it was difficult to separate the methyl ester from the glycerol phase (Fig. 1). There was no distinct separation phase for the corn methyl ester compared to that of vegetable (soybean) oil methyl ester as shown in Fig. 1. The corn methyl ester had significant soap formation and was mixed in the aqueous phase. It was important to investigate the best way to address the poor separation of corn methyl ester.
Poor separation of corn oil methyl ester can be explained by looking at some properties of the different oils tabulated in Table 2. Corn oil has the highest saponification value and ester value among three oils (Table 2). A low saponification value usually indicates easy soap formation at a higher temperature. Canola oil had the lowest saponification value (Table 2) causing more soap formation during transesterification reaction. The soap increases the methyl ester solubility in glycerol and results in emulsification of the ester and glycerol, which causes difficulties in the separation of the esters, thus reduces the recovered yield.
Table 2
The iodine, acid and saponification value of corn, canola and soybean oils including distribution of fatty acids
Oils | Saponification value (mg KOH g− 1) | Ester value (mg KOH g− 1) | Distribution of Fatty acid, % | |
Palmitic (16:0) | Stearic (18:0) | Oleic (18:1) | Linoleic (18:2) | Linolenic (18:3) | Others |
Corn | 193.3 | 182.4 | 13.2 | 1.6 | 28.2 | 54.7 | N/A | 2.2 |
Canola | 186.9 | 179.5 | 4.5 | 0.2 | 63.1 | 1.9 | 19.5 | 10.8 |
Soybean | 192.9 | 179.6 | 11.5 | 3.7 | 21.0 | 55.8 | 7.3 | 0.8 |
Saponification values include the neutral fatty acids and free fatty acid content present in the oil. On the other hand, the ester value represents the amount of neutral fat in the oil, which is directly associated with acid and saponification values. A high ester value indicates a high amount of ester with a low molecular weight fatty acid indicates the availability of short-chain fatty acid. Table 2 shows that the ester value is higher in corn oil than in other oils, and it may hinder the separation of glycerol and ester. Saponification is a side reaction happening along with transesterification. Saponification occurs when the free hydroxide of a catalyst breaks the ester bonds between fatty acids and glycerol in a triglyceride, resulting in more free fatty acids and glycerol. The ester bond is more prone to break down with short chain fatty acid. The sodium used for a catalyst is then bound with the fatty acid and unusable, thus complicating the separation and recovery of esters.
The soybean oil had a similar saponification value to the corn oil and a similar ester value to the canola oil (Table 2). Nevertheless, it is already evident from Fig. 1 that the soybean oil had less soap formation during reaction, when compared to the corn oil resulting in better methyl ester yield. Soybean oil had the lowest amount of oleic acid among others, and free fatty acid is often expressed as percent oleic acid. The presence of low oleic acid could be the reason for less soap formation during saponification towards transesterification of soybean oil.
To improve the recovered yield, the effect of stopping or neutralizing the reaction of the different corn, canola, and soybean methyl esters was studied. These are commonly used oils in the US and the expectation is that the yield of methyl ester production will be similar since the refining steps are also similar. The results illustrated in Fig. 2 revealed that neutralization has a positive effect on the yield of methyl ester. The recovered yield increased with the use of acid to stop the reaction. The methyl ester yield was 52–55% for corn and canola oil samples that were not neutralized. This yield significantly increased to 75% when the reaction was neutralized using acid. However, it was also visible that soybean methyl ester yield in reaction with no neutralization was around 72%, which was not significantly different with corn and canola methyl ester yield neutralized with acid. Neutralizing the reaction of soybean methyl ester significantly increased the yield to 88%. and was comparable to the maximum expected theoretical yield of 90%. The small variation might be an indication of some esters were washed during washing.
This finding is interesting because the recovered yield is dependent on how well the organic phase separated from the aqueous phase and also the formation of soap that usually dissolves in the aqueous phase. The higher recovered yield seen with soybean oil is likely due to less soap formation. Moreover, the added HCl is more strongly attracted to the metal ion on the sodium soap than the fatty acid chain. So, the metal ion combines with the Cl from the HCl to produce NaCl, and the hydrogen freed from the HCl converts the fatty acid chain to free fatty acid. In this way, the soap formation was reduced, and an increasing recovered yield was observed when acid was added after the transesterification reaction.
Most previous studies on the transesterification of these three mentioned oils were mainly focused on the property of methyl ester from different oils and their combustion performance. The findings of this present study is quite similar to Karademir and Karademir [25] who measured the efficiency of biodiesel production from soybean, corn and canola oil. Their result also showed a better ester conversion rate in soybean compared to corn or canola. But their yield measurement was based on total ester rate and the percent of linoleic acid conversion through gas chromatography instrument rather than the recovered yield. Moreover, the previous studies did not use any neutralization step during their experiment. The neutralization step of stopping transesterification reaction and the recovered yield on the basis of weight differentiates the present study from previous studies. Overall, the differences observed with the recovered yield of methyl ester justified why optimizing the process steps for each type of vegetable oil used is important. Hence, it was essential to optimize the process for corn oil which was the focus of this research.
Modelling the Factors Influence the Yield of Corn Oil Methyl Ester
With the finding of the importance of neutralization towards increasing the yield during transesterification, it was important to identify the amount of acid needed and its interaction with the total reaction time. The experimental runs corresponding to the factorial design along with the values of the yield for each run are presented in Table 1, where the predicted value is based on the statistical model. The model to predict the yield of methyl ester was first developed by considering both the linear and quadratic terms in Eq. 2. For developing a quadratic model, the present experimental design was expanded to a central composite design by the addition of three new runs (at time 1 h and stopping the reaction with 2.6 mL HCl). A statistical analysis was carried out on these experimental values, and the main effects and interaction effects of the factors were determined. Initially, both time (A) and square of time (A2) in the model were not statistically significant at p-value < 0.05. Therefore, a reduced model was developed using forward selection method as seen in Eq. (3) and the ANOVA results of the developed model are shown in Table 3.
$$\text{Yield (%) = 73.45-16.77A+0.31B-0.003}{\text{B}}^{\text{2}}\text{+0.22 AB}$$
3
Table 3
ANOVA results for the reduced quadratic model
Source of variation | Degree of freedom | Adjusted Sum of squares | Adjusted Mean Square | F-Value | P-Value |
Model | 4 | 3421.9 | 855.5 | 43.2 | < 0.0001 |
Linear | 2 | 2432.9 | 1216.4 | 61.4 | < 0.0001 |
Time (A) | 1 | 253.7 (6.0%) | 253.7 | 12.8 | 0.001 |
Acid amount (B) | 1 | 2179.2 (51.7%) | 2179.2 | 110.0 | < 0.0001 |
Acid Amount× Acid Amount (B2) | 1 | 539.2 (12.8%) | 539.2 | 27.2 | < 0.0001 |
Interactions (A×B) | 1 | 449.8 (10.7%) | 449.8 | 22.7 | < 0.0001 |
Error | 40 | 792.2 (18.8%) | 19.8 | | |
Lack-of-Fit | 10 | 339.6 (8.1%) | 34.0 | 2.3 | 0.042 |
Pure Error | 30 | 452.6 (10.7%) | 15.1 | | |
Total | 44 | 4214.1 | | | |
According to the ANOVA result (Table 3), the reaction time (A), the acid amounts used to stop the reaction (B), square term of acid amount (B2) and interaction term (A×B) significantly affected the yield at 95% confident level. The p-values (< 0.0001) indicate that all of the factors were important in modeling the transesterification reaction. The ANOVA table shows that almost 82% of the source of variation (SOV) in the process was due to variation in the levels of the selected factors. Acid amount (B) contributed more than 50% to the variation in the model. Both the squared term of acid amount (B2) and interaction term (A×B) had similar contribution to predict the yield of methyl ester (12.8% and 10.7%, respectively). Time (A) had contributed only 6% in the variation of the developed model. The source of variation for each factor is represented in the Pareto chart in Fig. 3. to present the magnitude and the importance of the effects. On the Pareto chart, bars that cross the reference line are statistically significant.
However, the lack of fit was slightly significant at p-value = 0.042 in Table 3, which was very close to being insignificant at p-value > 0.05. The significant error term was probably due to the elimination of square term of time (A2) from the initial model and higher relative deviation (greater than 5%) occurred in experimental run 1 to 6 (Table 1). The developed model indicates that for the reactions with no acid addition, the separation between methyl ester and glycerol becomes more difficult, as the reaction time increases.
The significance of each coefficient in Eq. (3) was evaluated by the p-value shown in Table 4. The smaller the magnitude of the P-value, the more significant is the corresponding coefficient. From Table 4, it can be seen that all terms in the model were found to be statistically significant. The acid amount (B) term had the most significant effect, followed by the square of acid amount (B2) and interaction (AB) term. The reaction time (A) had the least impact and a similar trend is also seen at 5% significance level with the Pareto chart at Fig. 3.
Table 4
Regression coefficients of the predicted model
Term | Coefficient | P-Value |
Constant | 74.80 | 0.000 |
Time (A) | -2.91 | 0.001 |
Acid Amount (B) | 9.84 | 0.000 |
Acid Amount × Acid Amount (B2) | -8.27 | 0.000 |
Time × Acid Amount (AB) | 5.48 | 0.000 |
Reaction time (A), amount of acid (B), and time-acid amount interaction effects (AB) were fitted by multiple regression analysis to develop a linear model. This linear model was developed to improve on the quadratic models. The adequacy of all the statistical models (linear, quadratic, quadratic with forward selection) was compared by model summary statistics (Table 5). The linear model maximizes the value of R2 and R2adj and minimizes the standard deviation compared to the other quadratic models. The high R2 value refers to the acceptable goodness of fit of the linear model. The R2 value 0.893 of the linear model revealed that 89.3% of the variation in the response was due to the difference observed in the factors. The linear model to predict methyl ester yield (response function) for the significant main effects and interactions is in Eq. (4).
$$\text{Y}\text{i}\text{e}\text{l}\text{d} \left(\text{%}\right) = 70.66-2.91\text{A}+9.84\text{B}+5.48 \text{A}\text{B}$$
4
Table 5
Model summary statistics of different model
Model type | Standard Deviation | R2 | R2adj | R2pred |
Linear | 3.885 | 0.893 | 0.843 | 0.758 |
Quadratic | 4.484 | 0.814 | 0.793 | 0.751 |
Quadratic with forward selection | 4.450 | 0.812 | 0.793 | 0.763 |
Main Effect and Interaction Effect on Corn Oil Methyl Ester Yield
Figure 4 and 5 illustrate the main effects and interaction effects of the treatment combinations on the yield of corn oil methyl ester. The model results indicate that both the main effects’ such as reaction time and acid amounts were statistically significant. Moreover, the interaction between the reaction time and amount of acid added was also significant. In the main effects plot (Fig. 4), reaction time at 0.5 h and 1 h were associated with the highest mean yield. The methyl ester yield at 1.5 h was significantly different than the yield at 0.5 h and 1 h. The decreasing trend of the yields with increasing time was also evident in Table 4. The negative coefficient of factor time (A) in Table 4 indicates a reduction in methyl ester yield with the increasing reaction time.
Highest yield was achieved by stopping the transesterification reaction with 5.2 mL of 6 M HCl. However, the highest yield was not statistically different for 2.6 mL and 3.9 mL of acid addition to the reaction. Lowest yield of methyl ester was found when no acid was added. Even with a little acid added to stop the reaction, the yield increases significantly (from 55.9% with no acid to 69.4% with 1.3 mL of acid). The positive coefficient of factor acid amount (B) in Table 4 indicates an increasing trend of methyl ester yield with the increasing acid amount. The negative coefficient of the square of acid amount (B2) indicates a concave curvature relationship with yield (Table 4; Fig. 4).
In the interaction plot (Fig. 5), the highest yield (79.1%) of methyl ester was achieved with the treatment condition of 1 h reaction time and 5.2 mL acid amounts used to stop the reaction. The significant differences of means were not illustrated here in Fig. 5, but it was already presented in Table 1 and Fig. 4. From Fig. 5, it was evident that the methyl ester yield of corn oil depends both on reaction time and the amount of acid added to stop the reaction. Figure 5 exhibited a trend of increasing yield with the increasing time and acid amount. The only exception happened when no acid was added to stop the reaction. When no acid was used, the yield of methyl ester reduced significantly from 69.6% at 0.5 h to 52.3% and 45.9% at 1.0 h and 1.5 h, respectively. From practical observation, the separation between glycerol and methyl ester became very difficult when no acid was added. Moreover, the interaction effect only occurs between 0.5-1 h. As the reaction time extended beyond 1 h, the interaction was not seen, and methyl ester yield was also slightly reduced.
Quality of Methyl Ester
Characteristics of the produced methyl ester from the experimental runs were tabulated and compared with standard methyl ester values and other previous studies as shown in Table 6. The selected characteristics were water content, viscosity, total glycerin content, acid value, cloud point, and pour point. For sucrose esters of fatty acid (bioresin) production using methyl esters, water content of the methyl ester is considered to be the most important parameter. This is because high water content slows down the catalysts by not only participating in the formation of emulsions but also causing hydrolysis or hydrolytic oxidation during the esterification reaction [26]. The water content of the produced corn oil methyl ester was well below the limit of ASTM D6751 (Table 6). Some of the previous studies have reported a very high water content of the methyl ester and might be an indication of using a different drying agent (diatomaceous earth) than that used in present study [27]. Both diatomaceous earth and MgSO4 have been reported to be good drying agents for oil but maybe the vacuum drying may have helped lower the moisture content. The amount of drying agent can also impact the water content in the final product. The kinematic viscosity and total glycerin content were also found to be within the limits of ASTM standards and comparable to the previous studies as well (Table 6). The low total glycerin content indicated that the conversion was good and very little impurity was present. The acid value found in this study was very similar to [28]. The low acid value may be due to the difference of acid value in the corn oil used.
Table 6
Characterization of corn oil methyl ester
Characteristics | ASTM D6751 (limit) | Corn oil methyl ester (present study)* | El Boulifi et al. [8] | Moser and Vaughn [28] | Mata et al. [27] |
Water content, ppm | < 500 | 242.05 ± 42.01 | 323 | 400 | 892 |
Viscosity, mm2 s− 1 | 1.9–6.0 | 4.55 ± 0.18 | 4.33 | 4.14 | 4.55 |
Total glycerin, % | < 0.24 | 0.05 ± 0.01 | 0.18 | nd | nd |
Acid value, mg KOH g− 1 | < 0.5 | 0.49 ± 0.08 | 0.04 | 0.49 | 0.25 |
Cloud point, °C | Report | -1.96 ± 0.51 | -3.6 | -5 | nd |
Pour point, °C | Report | -4.41 ± 2.06 | -6 | -6 | nd |
*n = 30 (15 sample with two replicates) displayed as average ± standard deviation, nd: not determined |
Cloud point and pour point was a main indicator of cold flow properties of methyl esters. However, ASTM D6751 requires that cloud point and pour point be reported. The observed cloud point and pour point values for the present study were − 1.96°C and − 4.41°C, respectively. Other studies have reported lower cloud point and pour point for their methyl esters, indicating good cold flow properties (Table 6). In our study, this characteristic has less attention as it was mainly related with fuel properties. Overall, it was evident that despite having a low yield in some of the experimental runs, the quality was still within acceptable limits.