Esterification Optimization of Crude African Palm Olein Using Response Surface Methodology and Heterogeneous Acid Catalysis

In this work, the effect of zeolite montmorillonite KSF in the esterification of free fatty acids (FFAs) of crude African palm olein (Eleaias guinnesis Jacq) was studied. To optimize the esterification of FFAs of the crude African palm olein (CAPO), the response surface methodology (RSM) that was based on a central composite rotatable design (CCRD) was used. The effects of three parameters were investigated: (a) catalyst loading (2.6–9.4 wt %), (b) reaction temperature (133.2–166.2 ◦C), and (c) reaction time (0.32–3.68 h). The Analysis of variance (ANOVA) indicated that linear terms of catalyst loading (X1), reaction temperature (X2), the quadratic term of catalyst loading (X2 1), temperature reaction (X2 2), reaction time (X 2 3), the interaction catalyst loading with reaction time (X ∗ 1 X3), and the interaction reaction temperature with reaction time (X∗ 2 X3) have a significant effect (p < 0.05 with a 95% confidence level) on Fatty Methyl Ester (FAME) yield. The result indicated that the optimum reaction conditions to esterification of FFAs were: catalyst loading 9.4 wt %, reaction temperature 155.5 ◦C, and 3.3 h for reaction time, respectively. Under these conditions, the numerical estimation of FAME yield was 91.81 wt %. This result was experimentally validated obtaining a difference of 1.7% FAME yield, with respect to simulated values.


Introduction
Fossil fuels such as oil, coal, and gas represent the main energy source worldwide.However, these turn out to be limited resources and contribute to environmental pollution and global warming.This scenario demands policies that stimulate the use of renewable and sustainable energy sources in order to mitigate these effects [1,2].In this context, biodiesel has been glimpsed as a viable replace respect to engine conventional fuels.Among the several biodiesel alternatives, the use of Fatty Methyl Ester (FAME) in combustion engines can be found [3].Biodiesel has advantages over conventional diesel fuel, such as it can be obtained from oils and animal fats; also, it is biodegradable, non-toxic, and environmentally sustainable, avoiding incorporate new CO 2 to the carbon cycle [4,5].
Currently, the biggest obstacle to the substitution of conventional diesel by the biodiesel is its high costs, where the feedstock represents 75% of its manufacturing cost [6].Therefore, it is necessary to use lower-cost feedstocks in its manufacture.Among the alternatives for producing biodiesel at a reasonable cost, are the cooking vegetable oil, agricultural waste, animal fats, and unrefined vegetable oils [7,8].For example, the cost of biodiesel made from used cooking oil and animal fats waste is estimated between $0.09-$0.20 USD per pound [9].The cost of biodiesel from crude African palm oil has been estimated between $0.49 and $0.53 USD per liter [10].Unfortunately, those feedstocks have a high content of free fatty acids (FFAs) which present soap formation during the process of conventional direct base-catalyzed transesterification, making their use not viable [11].As an alternative, the use of a multistep process allows employing feedstocks with high FFAs concentrations by first carrying out the acid-catalyzed pre-esterification of the FFAs prior to the base-catalyzed triglycerides transesterification [12].
The usual esterification procedure implies the use of p-toluene sulphonic acid, sulfuric, hydrochloric, and phosphoric acids as a catalyst.However, the use of these catalysts has several disadvantages: (a) during the esterification reaction water is generated, favoring the reversibility of the reaction, which lowers the yield of methyl esters; (b) these acids are highly corrosive, which reduces the lifetime of the equipment used in the process; (c) toxic wastes that are harmful to the environment are generated; and, (d) costs of purification of biodiesel are high increasing the biodiesel cost [13].Another method that can be used for esterification of FFAs is the supercritical process.Its main advantages consist of reduce mass-transfer limitations and improve the phase solubility.Furthermore, the reaction rate is augmented, completing it in shorter periods and simpler separations and purification of FAME [14].Nevertheless, to present satisfactory FAME yield the supercritical method needs high molarity of alcohol to oil and high pressure and temperature conditions for the reaction.This implies high processing costs and causes the degradation of the fatty acids esters formed [15].
A very attractive alternative to carry out the esterification of FFAs is the use of heterogeneous catalyst systems, which include cationic-exchanged resins.
Examples of these resins are Amberlyst-15, metal oxides, heteropolyacids supported on silica and carbon, aluminophosphates, and silicoalumino-phosphate molecular sieves and clays [16].The clays have presented a substantial interest in several organic syntheses as catalysts [17].Among these, the montmorillonite possesses both Lewis and Brönsted acidities sites.Due to these properties, there is particular interest in using montmorillonite as catalysts in the biodiesel manufacture, as they have several advantages: (a) availability and low cost; (b) environmental compatible; (c) easy recovered by filtration and reusable; (d) high selectivity and thermal stability; and, (e) allow the development of continuous processes to increase the production volume of biodiesel, which reduces manufacturing costs.In particular, the montmorillonite KSF (commercial name for clays activated with acid treatment) presents high concentrations of Lewis and Brönsted acid sites, thus oil esterification can be employed with high FFAs concentrations [18].
The esterification reaction of FFAs is represented by the following reversible chemical equation This reaction is affected by temperature changes, reaction time, and the change of catalyst and alcohol concentrations.These variations can enhance or reduce the FFAs to FAME conversion.For this reason, it is important to evaluate the effect of the process variables to determine the optimal reaction conditions of FFAs to FAME conversion and to maximize the yield and reduce costs in the manufacture of biodiesel.It is also necessary to develop the experimental design Response Surface Methodology (RSM) and a quadratic model to calculate the best fit.The RSM is a group of mathematical and statistical tools of important use in the modeling, analysis and optimization problems.The response of interest is influenced by various factors in RMS [11].In this methodology, a suitable experiment can significantly reduce the number of tests, and provide the information needful for process optimization.The RSM uses polynomial models to characterize the response of the surface [19].
The RSM has been widely reported as a suitable tool for the modeling and optimization of biodiesel yield, for example: Jahirul et al. [11] used RSM based in a Box-Behnken design to determine a method to optimize the biodiesel production process from beauty-leaf oil.The authors considered catalyst concentration, methanol to oil molar ratio and reaction temperature as the main factors in the conversion process.FFAs and FAME concentration prediction and reaction conditions optimization were achieved through linear and full quadratic regression models.Ali et al. [20] investigated on fuel properties and engine performance of blended palm biodiesel-diesel using diethyl ether as an additive.The independent variables considered were the engine revolution per minute and the diethyl ether blend.A design of experiment based on RSM was employed to optimize the fuel properties viscosity, density, acid value, heating value, pour point, and cloud point.Corral Bobadilla et al. [21] used RSM to improve the production of biodiesel by finding the best combination of input variables used in the transesterification reactions.Various biodiesel optimization scenarios were proposed based on the objective to improve the biodiesel production and higher heating value, and decrease the viscosity, density, and turbidity.Omar et al. [22] used the RSM to study the relationship of different variables in the waste palm cooking oil heterogeneous transesterification to biodiesel by Sr/ZrO 2 catalyst.
According to Ramos et al. [23], soybean, canola, sunflower, and African palm (Elaeis guinnensis Jacq) oils are the most widely used for the production of biodiesel in the world.The African palm is an oil-producing species, its expansion around the world is due to its easy adaptation to tropical environments.The high oil content that has its fruit allows for yields of 5.5 ton ha −1 , generating the decrease of production costs.This makes the crude African palm oil (CAPO) a viable option for its use as raw material in the manufacture of low-cost biodiesel [24].The main producing CAPO countries in the world are Malaysia, Indonesia, Thailand, Nigeria, and Colombia, among others.In Mexico, the African palm cultivation takes place in the states of Chiapas, Campeche, Tabasco, and Veracruz with 54,434.00ha planted in 2010 [24].There is a significant possibility to project 2.5 million ha of African palm plantation in these areas.Based on these conditions, the African palm may be considered as a potential alternative to increasing the manufacture of biodiesel in Mexico [25].
The objective of this study is to evaluate the efficiency of the montmorillonite KSF for the esterification of the FFAs present in the CAPO.Furthermore, the determination the optimal reaction conditions that yield the optimization of FFAs to FAME conversion at a laboratory level, through the application of a central composite rotatable design (CCRD) and RSM techniques.
The esterification reaction was carried out by using a batch-type tube reactor (200 mL and external diameter = 38.1 mm) made of glass (Vineland).The esterification reaction used a molar ratio methanol/oil 12:1 and a palm oil quantity of 20 g for each experiment.Initially, oil, methanol, and montmorillonite KSF were charged into the tube reactor.The tube reactor was immersed in an oil bath.Previously, the bath oil was preheated to the reaction temperature of each experiment.A mercury thermometer was used to measure the temperature of the oil bath.The stirring and heating of the mixture of reactants and heating oil bath were performed using a magnetic hot plate stirrer with a digital temperature controller (Thermoscientific, Waltham, MA USA).After a fixed reaction period, the tube reactor will be transferred into another oil bath (30 • C), in order to quench the reaction immediately (5 min) due to the short extern diameter of the reactor.The mixture reaction products containing clay, methanol, water, FAME, and unreacted triglycerides.The mixture was filtered through Whatman filter paper #40 to recover the clay.To remove methanol and water from the mixture obtained by filtration, this was subjected to vacuum evaporation using a Yamato rotary evaporator Re-801, under the following operating conditions: Rotation speed of 220 rpm, a vacuum pressure of 480 mm Hg and temperature of 80 • C in water bath.Finally, the mixture that was obtained from the evaporation contains FAME and unreacted triglyceride was transferred into test tubes to determine acid.The experimental conversion (yield) of FFAs to FAME Y exp can be calculated by Equation (1): where V 0 and V 1 are the acid values of feed and products, respectively.The acid value was determined by the method developed by Aricetti et al. [26].This novel method presents the main advantages of being less expensive and nontoxic since it uses ethanol and water as solvents.Besides, the evaluation of this method exhibited a similar accuracy and it is comparable to the official method AOCS Cd 3d-63 [27], since there is no statistically significant difference between both.

Characterization of MT-KSF Catalyst
X-ray diffraction analysis (XRD) is a widely used method to characterize catalyst in different materials.XRD analysis was performed in the laboratory using an APD 2000 PRO X-Ray diffractometer with CuKα radiation (λ = 1.5405Å).The patterns for of Montmorillonite-KSF (MT-KSF) sample was recorded in the 2θ range of 5 • to 75 • with a scanning speed of 0.025 deg s −1 and time steps of 10 s.The obtained crystalline phases of the sample were identified by comparison with Joint Committee of the Power Diffraction Standards (JCPDS) cards.
Scanning electron microscopy analysis (SEM) is used to study the morphological features and surface characteristics of MT-KSF sample.Chemical composition was determined by energy-dispersive X-ray spectroscopy at different points on the surface of the samples.SEM-EDS analysis was performed using HITACHI S-3400N Microscope instrument.The working conditions of SEM-EDS were 10 kV for the energy dispersive analysis, with a current intensity of 30 pA, and a work distance of 10.6 mm.The infrared spectroscopy analysis of the MT-KSF sample was conducted using an Agilent Cary 660 infrared spectrometer with ATR imaging with diamond crystal design (Gladi Pike).The equipment was operated within a range from 4000 to 400 cm −1 , with a resolution of 4 cm −1 in a KBr wafer.
X-ray diffraction: XRD pattern of the MT-KSF sample is illustrated in Figure 1.It can be observed the natural clay presence of the main reflections at 2θ = 8.90 in a KBr wafer.
X-ray diffraction: XRD pattern of the MT-KSF sample is illustrated in Figure 1.It can be observed the natural clay presence of the main reflections at 2θ = 8.90°, 17.71°, 19.85°, 27.95°, 37.60°, and 61.88°, which correspond to the crystalline structure of the MT-KSF reported in the JCPDS standards (00-025-1349 card).Minor amounts of quartz at 2θ = 20.91°,26.56°, 39.43°, 45.58°, 50.05°, 59.89°, 68.24° (JCPDS 01-085-0504) and Calcite at 2θ = 28.04°were also found.Scanning Electron Microscopy and Energy Dispersive X-ray Spectroscopy (SEM/EDS): SEM image of MT-KSF sample is given in Figure 2. It can be observed that the surface morphology demonstrates a layered surface with some flakes, which is the typical structure of montmorillonite.Moreover, the surface of the clay does not have homogeneous dispersion and the most crystals display characteristic monoclinic symmetry, where many of them are fluffy-shaped.This confirms the results obtained by XRD analysis.Scanning Electron Microscopy and Energy Dispersive X-ray Spectroscopy (SEM/EDS): SEM image of MT-KSF sample is given in Figure 2. It can be observed that the surface morphology demonstrates a layered surface with some flakes, which is the typical structure of montmorillonite.Moreover, the surface of the clay does not have homogeneous dispersion and the most crystals display characteristic monoclinic symmetry, where many of them are fluffy-shaped.This confirms the results obtained by XRD analysis.The Chemical composition of the MT-KSF sample was determined by the EDS.The results of the punctual analysis at different zones Figure 2a,b are showed in Table 1.The total weight percentage indicates that the chemical elements constituting the montmorillonite KSF are: oxygen, magnesium, aluminum, silicon, sulfur, potassium, calcium, and iron.On the other hand, it was found that oxygen, silicon, aluminum, and sulfur are the most abundant elements in the clay, with 72.49%, 12.12%, 5.77%, The Chemical composition of the MT-KSF sample was determined by the EDS.The results of the punctual analysis at different zones Figure 2a,b are showed in Table 1.The total weight percentage indicates that the chemical elements constituting the montmorillonite KSF are: oxygen, magnesium, aluminum, silicon, sulfur, potassium, calcium, and iron.On the other hand, it was found that oxygen, silicon, aluminum, and sulfur are the most abundant elements in the clay, with 72.49%, 12.12%, 5.77%, and 5.32% weight, respectively.In addition, the average elemental composition of clay (Table 2) tested that the ratio Si/Al is about 2.62 characteristic to MT-KSF.The EDS analysis indicates that the MT-KSF presents high percentages of oxygen, silicon, aluminum, and sulfur in the clay, contributing to the formation of Al-O-Si, Al-O-H, and S-OH, Si-O-H chemical bonds that constitute the skeleton of the clays [28][29][30].The presence of these chemical bonds was confirmed by infrared spectroscopy analysis.These chemical bonds forming Lewis and Brönsted acid sites, the molar ratio of Si/Al relate the amount of acid sites present in the clay.Low values of this ratio (Si/Al) indicate that the MT-KSF has high catalytic activity [28-32].

Response Surface Methodology (RSM)
RSM is a technique based on the design of experiments (DOE) that is used to establish the relationship between independent variables and one or more dependent variables [33].It is useful for modeling, analysis, and optimization industrial process and products [34,35].RMS consist of statistical methods that use a polynomial model that relies on a low degree polynomial function Equation (2): where Y mod represents the relationship of the response variables, f is a function of cross-products of the polynomial's terms, X 1 , X 2 , X 3 . . . .X k are the independent variables, and ε is the noise or error.The quadratic model is a polynomial function for predicting the optimal points is expressed as follows: where Y mod represents the value of the predicted response, β 0 , β i , β ii , and β ij are the regression coefficient (β 0 is the intercept coefficient (offset), β i is the linear effect term, β ii is the pure quadratic Energies 2018, 11, 157 terms, and β ij is the interaction effect term), X i , X j , and X 2 i are the independent variables, and k is the total number of independent variables, ε is a random error.The second order model equation, RSM and contour plots were used to predict the optimum value and to analyze the interaction between response (Y) and the variables of the process.The second order model equation is recommended since the reaction rate of esterification of FFAs is governed by the second order reversible reaction [36,37].
The value of p (p-value) is a statistical parameter, which indicates the significant effect of an independent variable on the dependent variable [5].If p-value < 0.05, then the variable is significant in the process.The p-value (or Prob > F), computed by an analysis of variance (ANOVA), is used to determine the probability of generating a result equal or greater than the observed.If it exceeds the model's F and the model has no terms that exceed the significance level (e.g., p = 0.05) it will suffice with a confidential interval of 1-p [21].
In order to generalize the capacity prediction of the model the Mean Absolute Percentage Error (MAPE, which is the computed average of errors (%) by estimating the predictions of a variable), Root Mean Square Error (RMSE, used to compute the differences between estimated and observed values), and the adjusted determination coefficient (R adj 2 , a statistical measure that indicates the proportion of variation explained by the estimated regression line and that determine the predictive capacity of the model rigorously between 0 and 1) are calculated [38]: in this case, Y exp represents the experimental response, Y mod is the model output, − Y exp is the average experimental value, m represents the experiments produced, and k the number of independent variables.

Experimental Design
In this work, a central composite rotatable design (CCRD) [33] was used to study the effect of three factors on the FAME yield: the catalyst loading (X 1 ), reaction temperature (X 2 ), and reaction time (X 3 ).The study was carried out in a total of 20 experiments with triple repetition, previously randomized.The experimental design considers eight factorial points (2 3 ), six axial points, and six central points, according to the established by Montgomery [32].As is enlisted in Table 3, the experimental results of the esterification process reported a FAME yield range (Y exp ) from 34.98 to 88.91 wt %, depending on the conditions of experiments.The DOE consisted of eight factorials ( 23 ), six central, and six axial points, respectively.The distance of star points from the center point is provided by α = (2 n ) 1/4 , where the number of independent factors is given by n [32], of three factors α = 1.68.The range and coded levels of the esterification process variables studied are listed in Table 4.The yield of FAME was the response of the experimental design.

Modeling Results
To model the FAME yield, The Equation ( 3) was fitted using the data presented in Table 3.The experimental data were analyzed using the statistical package STATISTICA V.7 to construct the second-order polynomial model.The Levenberg-Marquardt algorithm was used to compute the parameters of the second-order model for nonlinear multiple regression [39].This analysis includes all of the independent variables and their binary interactions, indicating the real influence of each variable in Equation ( 7): We have considered 95% of confidence to evaluate the capacity of predict model.The remaining terms (the variables with a p-value > 0.05 had no significant effect on the response variable) despite not respect to the experimental data.The result of the determination coefficient adjusted (Radj 2 = 0.9723) determines that the sample variation of 97.23% for FAME yield is related to the independent variables, where only 2.77% of the total variations are not explained by the mathematic model.In general, the statistical results validate the accuracy and prediction capacity of the second-order model, guaranteeing the reliability of the model.Figure 4 presents the standardized residual values with zero mean and unit variance, of the regression model (Figure 3).As it can be appreciated, all of the residual values are included in the ±3 interval, determining the absence of outlier points.Therefore, the model predicts with certainty the experimental data, agreement with Montgomery [29].Figure 4 presents the standardized residual values with zero mean and unit variance, of the regression model (Figure 3).As it can be appreciated, all of the residual values are included in the ±3 interval, determining the absence of outlier points.Therefore, the model predicts with certainty the experimental data, agreement with Montgomery [29].

Analysis of Variance (ANOVA)
The ANOVA was used to evaluate the statistical significance model p-value by regression and Mean Square of Residual Error values.Table 6 contains the results of this analysis.

Analysis of Variance (ANOVA)
The ANOVA was used to evaluate the statistical significance model p-value by regression and Mean Square of Residual Error values.Table 6 contains the results of this analysis.According to Table 6, the ANOVA results also determine that the most significant linear terms on FAME yield (with p-values < 0.05) are the reaction temperature (X 2 ), the quadratic term of reaction time (X 2 3 ), and the quadratic reaction temperature (X 2 2 ).Likewise, catalyst concentration (X 1 ), the quadratic term of reaction time (X 2 1 ), the interaction catalyst loading*reaction time (X * 1 X 3 ), and the interaction reaction temperature*reaction time (X * 2 X 3 ) were also significant on FAME yield.Finally, the interactions catalyst loading*reaction temperature (X * 1 X 2 ) and the linear term of reaction time (X 3 ) had no significant effect on the FAME yield (p-value > 0.05).As it can be seen, although X 3 does not present significance on its own, its importance when interacting with the two other variables determines that it is a required element for the model result.Therefore, according to the ANOVA, the three input variables have a relevance importance on the model response.

Optimization
Optimization seeks to find the desired catalyst concentration, reaction temperature, and reaction time to obtain a required FAME yield.The second-order model (Equation (7)) illustrated the relationship between the response and experimental levels of each factor.This equation can be implemented to deduce the optimum conditions by response surface and the contour plots.To compute the optimum FAME yield, the empirical model was plotted in a three-dimensional surface that represent the FAME yield response as a function of two variables within the experimental range considered (Figure 5). Figure 5a illustrates the interaction of catalyst loading (2.6-9.4 wt %) and reaction temperature (133.2-166.2• C) on FAME yield, respectively.It can be observed that catalyst concentrations that are lower than 6 wt % and reaction temperatures lower than 150 • C, yield to FAME concentrations of 72 wt % or smaller.On the contrary, when catalyst concentration increases higher than 8 wt % and reaction temperatures are higher than 150 • C, FAME concentration increases significantly.This increase may be caused by an increment in the catalyst concentration, which also increases the number of activation sites.Since, the results of the EDX analysis show that montmorillonite KSF contains 5.75% and 12.12% by weight of aluminum and silicon, respectively.The presence of aluminum and silicon will generate Brönsted acid sites in the crystal structure of the clay.On the other hand, the presence of Fe, Mg, K, Ca, and S atoms in the montmorillonite structure will form bonds with oxygen to generate Lewis acid sites.These acid sites decrease the activation energy, which accelerates the rate of esterification reaction favoring the synthesis of FAME.According to Arrhenius' law, when the reaction temperature increases, the velocity of reactants reaction increases, and improve the FAME yield [42].For catalysis, concentrations higher than 8 wt % and reaction temperatures higher than 160 • C, FAME concentration of 88.91% was reached.These results are consistent to those that were reported in previous studies by Boffito et al. for the esterification of FFAs of waste cooking oil (WCO) using ionic resins as catalytic agents.For reaction times of 8 h, the reaction temperature of 63 • C, and a relation methanol WCO (16:100 wt %), FAME conversions of 90 wt % were achieved [43].
Temperature is one the most important variable affecting the conversion of FFAs.The FAME yield, along with the reaction time and different reaction temperatures, is illustrated in Figure 6b.In to optimize the FAME yield, the effects of different reaction temperatures (133.2-166.2• C) and reaction times (0.32-3.68 h) on the FAME yield was investigated.As is illustrated in Figure 6b, with increasing reaction temperature, the FAME yield increased quickly since the temperature increase caused higher molecule motion speed and mass transfer rate.For elemental reactions at reaction temperature increments of 10 • C, the reaction rate is doubled [40], which enhances FAME concentration.For reaction times of 3 h or longer, and temperatures of 160 • C, FAME yield of 90 wt % were achieved.Nandiwale et al. studied the esterification of levulininc acid to obtain ethyl levulitate, using an acid silicate (H-ZSM-5) compound as a catalytic agent.For a reaction time of 7 h and a reaction temperature of 130 • C, FAME concentration of 95 wt % was achieved [44].
According to Arrhenius' law, when the reaction temperature increases, the velocity of reactants reaction increases, and improve the FAME yield [42].For catalysis, concentrations higher than 8 wt % and reaction temperatures higher than 160 °C, FAME concentration of 88.91% was reached.These results are consistent to those that were reported in previous studies by Boffito et al. for the esterification of FFAs of waste cooking oil (WCO) using ionic resins as catalytic agents.For reaction times of 8 h, the reaction temperature of 63 °C, and a relation methanol WCO (16:100 wt %), FAME conversions of 90 wt % were achieved [43].
Temperature is one the most important variable affecting the conversion of FFAs.The FAME yield, along with the reaction time and different reaction temperatures, is illustrated in Figure 6b.In to optimize the FAME yield, the effects of different reaction temperatures (133.2-166.2°C) and reaction times (0.32-3.68 h) on the FAME yield was investigated.As is illustrated in Figure 6b, with increasing reaction temperature, the FAME yield increased quickly since the temperature increase caused higher molecule motion speed and mass transfer rate.For elemental reactions at reaction temperature increments of 10 °C, the reaction rate is doubled [40], which enhances FAME concentration.For reaction times of 3 h or longer, and temperatures of 160 °C, FAME yield of 90 wt % were achieved.Nandiwale et al. studied the esterification of levulininc acid to obtain ethyl levulitate, using an acid silicate (H-ZSM-5) compound as a catalytic agent.For a reaction time of 7 h and a reaction temperature of 130 °C, FAME concentration of 95 wt % was achieved [44].The response surface and the two dimensions contour plots are generally used to find the optimal values of the variables process.This technique allows for inferring the optimum FAME yield, and its respective operational parameters, from the FAME curve gradient.The contour plots are illustrated in Figure 6.Each contour curve represents an infinite number of two test variable combination (reaction temperature (X 2 ) and reaction time (X 3 )) with the other one maintained at their respective zero level (catalyst loading (X 1 )).The contour plots in Figure 6 indicate that there is a significant mutual interaction between reaction temperature and reaction time.Figure 6 predicts that the optimal values of the test are reaction temperature 155.5 • C, reaction time 3.3 h, and catalyst concentration 9.4 wt %.The mathematical model predicts that under optimal conditions obtained 91.81 wt % FAME yield.Further increase in the reaction temperature or reaction time would not have positive effects on the FAME yield.
In order to validate the optimal experimental condition that was generated by the mathematical model, an experiment was conducted in the laboratory considering the following conditions: reaction temperature 155 • C, reaction time 3 h 20 min (equivalent to 3.3 h) and catalyst concentration 9.5 wt %.The experimental evaluation generated a 90.1% FAME yield, with a difference of 1.71% between the modeled output and experimental result.This shows, the efficiency of the computational methodology to estimate the reaction optimum conditions.
In order to validate the optimal experimental condition that was generated by the mathematical model, an experiment was conducted in the laboratory considering the following conditions: reaction temperature 155 °C, reaction time 3 h 20 min (equivalent to 3.3 h) and catalyst concentration 9.5 wt %.The experimental evaluation generated a 90.1% FAME yield, with a difference of 1.71% between the modeled output and experimental result.This shows, the efficiency of the computational methodology to estimate the reaction optimum conditions.

Conclusions
This study is focused in optimize the Fatty Methyl Ester yield of the effect of montmorillonite KSF in the esterification of FFAs of CAPO.The operation parameters catalyst loading, reaction time, and reaction temperature were optimized using a response surface methodology that was based on a central composite rotatable design to achieve the maximum Fatty Methyl Ester yield.According to the statistical results, the second-order model generated is a feasible tool for estimation of Fatty Methyl Ester yield with minimums error levels (MAPE = 2.13, RMSE = 1.43, and R 2 = 97.23%).Additionally, ANOVA results demonstrate that the three operation parameters present relative importance for the model output.
According to the experimental results, the montmorillonite KSF could be used as an effective heterogeneous catalyst for the conversion of FFAs of CAPO to biodiesel.The optimum reaction

Conclusions
This study is focused in optimize the Fatty Methyl Ester yield of the effect of montmorillonite KSF in the esterification of FFAs of CAPO.The operation parameters catalyst loading, reaction time, and reaction temperature were optimized using a response surface methodology that was based on a central composite rotatable design to achieve the maximum Fatty Methyl Ester yield.According to the statistical results, the second-order model generated is a feasible tool for estimation of Fatty Methyl Ester yield with minimums error levels (MAPE = 2.13, RMSE = 1.43, and R 2 = 97.23%).Additionally, ANOVA results demonstrate that the three operation parameters present relative importance for the model output.
According to the experimental results, the montmorillonite KSF could be used as an effective heterogeneous catalyst for the conversion of FFAs of CAPO to biodiesel.The optimum reaction conditions were a catalyst loading 9.4 wt %, reaction temperature 155.5 • C and reaction time, 3.3 h.Under these conditions, was possible to obtain 91.81 wt % Fatty Methyl Ester yield.After the esterification steps, the FFAs to less than 1.1 wt %, this allows for the use of alkaline catalysis.This would allow for developing a two-step process for the synthesis of biodiesel, with high efficiency and low energy cost.

Figure 3 .
Figure 3. Scatter diagram of experimental FAME yield and predicted FAME yield.

Figure 4 .
Figure 4. Diagram of standardized residuals of predicted FAME yield.

Figure 4 .
Figure 4. Diagram of standardized residuals of predicted FAME yield.

Figure 5 .
Figure 5. FAME yield response surface of the free fatty acids (FFAs) esterification as a function: (a) The catalyst loading and temperature reaction; and (b) reaction temperature and reaction time.

Figure 6 .
Figure 6.Graph of contour illustrating the optimum reaction temperature and reaction time for the esterification FFAs of CAPO.

Figure 6 .
Figure 6.Graph of contour illustrating the optimum reaction temperature and reaction time for the esterification FFAs of CAPO.

Table 2 .
Average chemical composition of MT-KSF sample by EDS analysis.

Table 4 .
Experimental range and levels coded of independent variables.

Table 6 .
ANOVA for the model regression.

Table 6 .
ANOVA for the model regression.