Next Article in Journal
Reduction in Operating Costs and Environmental Impact Consisting in the Modernization of the Low-Power Cylindrical Wood Chipper Power Unit by Using Alternative Fuel
Next Article in Special Issue
A Review on Bio-Based Catalysts (Immobilized Enzymes) Used for Biodiesel Production
Previous Article in Journal
An Energy Management System-Based Control Strategy for DC Microgrids with Dual Energy Storage Systems
Previous Article in Special Issue
Development of Full-Cycle Utilization of Chlorella sorokiniana Microalgae Biomass for Environmental and Food Purposes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Universal Kinetic Model to Simulate Two-Step Biodiesel Production from Vegetable Oil

1
Department of Physical Chemistry and Applied Thermodynamics, EPS, Edificio Leonardo da Vinci, Campus de Rabanales, Universidad de Cordoba, Campus de Excelencia Internacional Agroalimentario ceiA3, 14071 Córdoba, Spain
2
Department of Rural Engineering, EPS, Edificio Leonardo da Vinci, Campus de Rabanales, Universidad de Cordoba, Campus de Excelencia Internacional Agroalimentario ceiA3, 14071 Córdoba, Spain
*
Author to whom correspondence should be addressed.
Energies 2020, 13(11), 2994; https://doi.org/10.3390/en13112994
Submission received: 9 May 2020 / Revised: 2 June 2020 / Accepted: 5 June 2020 / Published: 10 June 2020
(This article belongs to the Collection Feature Papers in Bio-Energy)

Abstract

:
To date, to simulate biodiesel production, kinetic models from different authors have been provided, each one usually applied to the use of a specific vegetable oil and experimental conditions. Models, which may include esterification, besides transesterification simulation, were validated with their own experimental conditions and raw material. Moreover, information about the intermediate reaction steps, besides catalyst concentration variation, is either rare or nonexistent. Here, in this work, a universal mathematical model comprising the chemical kinetics of a two-step (esterification and transesterification) vegetable oil-based biodiesel reaction is proposed. The proposed model is universal, as it may simulate any vegetable oil biodiesel reaction from the literature. For this purpose, a mathematical model using the software MATLAB has been designed. Using the mathematical model, the estimation of mass variation with time, of both reactants and products, as well as glyceride conversion and homogeneous catalyst concentration variation (instead of only alcohol/catalyst solution) are allowed. Moreover, analysis of the influence of some important variables affecting the reaction kinetics of biodiesel production (e.g., catalyst concentration), along with comparison and model validation with data from different authors may be carried out. In addition, Supplementary material with a collection of 290 rate constants, derived from 55 different experiments using different vegetable oils and conditions is provided.

1. Introduction

In the global effort to generate non-contaminant renewable fuels, biofuels comprising biodiesel, bioethanol, biobutanol and biogas, among others, play a significant role. Nowadays, many researchers work to improve processing efficiency and the integration of cellulose (to date, expensive) and organic waste recycling, as raw materials, in economically sound processes [1,2,3].
Biodiesel is a green fuel that can be used in diesel engines instead of fossil fuels. It is mainly produced through transesterification of vegetable oils/fats or animal fats, including further purification steps. Biodiesel is then used in diesel engines as a substitute for diesel fuel, straight or mixed with diesel fuel. Furthermore, it is well known that, during the process of combustion, biodiesel exhaust emissions, including carbon monoxide, hydrocarbons and particulates, are lower than those of diesel fuel, providing zero emissions of sulfur, lead and other heavy metals, besides the added advantage of its powerful lubricating properties [4].
The debate on food versus fuels focuses the discussion on the influence of biofuel prices on food prices and shortages, that are especially significant during a food crisis [5]. However, the crude oil price seems to be the main cause of food price increases [6]. Moreover, traditionally, land has been used to grow nonfood crops, e.g., tobacco or crops to produce alcoholic beverages, so its use to produce energy should be considered at the same level of acceptance. In any case, if land is not used to produce energy, food crop residues may still be used to provide biofuels, thus minimizing the polluting effect of these wastes [7,8].
In the search for new alternative low-cost raw materials, the reuse of edible frying oil may become an interesting raw material. Although this oil exhibits a high amount of free fatty acids (FFA), it may be used to produce biodiesel with minor pretreatment [9]. Other affordable commodities, which may be used as feedstock to produce biodiesel, are alternative vegetable oils, e.g., Jatropha curcas oil, also characterized by a high content of FFA. This represents a key problem for the common alkaline transesterification process because the alkaline catalyst may react with FFA, thus producing soap (saponification reaction), with consequent reduction of biodiesel yield and compounding separation of esters, glycerol and washing water. Soap formation also increases viscosity and leads to the formation of a gelatinous final product. In conclusion, the use of a raw material with an FFA content above 0.5% is not recommended for alkali-catalyzed transesterification reactions; therefore, a previous step consisting in FFA esterification with an acid catalyst (typically, H2SO4) is required. This previous step is followed by a second one, consisting in the transesterification of triacylglycerols by alcoholysis, usually with a basic catalyst.
During esterification, FFAs react with an alcohol to produce fatty acid ester and water. Subsequently, in the transesterification reaction, triacylglycerols react with an alcohol to produce other esters of the same fatty acids and glycerol. This reaction can be written as three consecutive reversible reactions with intermediate formation of diacylglycerols and monoacylglycerols. Finally, biodiesel is separated from glycerol by decantation or centrifugation [10]. Methanol is the most commonly used alcohol and, in this case, the reaction is known as methanolysis. There are physical and chemical properties that affect both esterification and transesterification reactions, i.e., the source of triacylglycerols, reaction temperature, type and amount of catalyst, type and amount of alcohol, stirring time, etc. [11]. Although this is a well-known process, the ratio of reactants affects the process in terms of conversion efficiency and this factor differs depending on the raw material chemical composition.
There are several kinetic studies about biodiesel production from vegetable oils under different reaction conditions. Esterification using H2SO4 as a catalyst has been carried out with methanol and either coconut or sunflower oil, and with ethanol and sunflower oil, at different temperatures and concentrations [12,13]. In the same way, methanolysis has been carried out using KOH or NaOH as catalyst, using rapeseed oil [14], sunflower oil [15,16], Pongamia pinnata oil [17], Brassica carinata oil [18], palm oil [19] or soybean oil [20,21] also at different temperatures and concentrations. However, only a few kinetic studies about enzyme- or lipase-assisted transesterification have been reported [22,23,24]. More recently, the use of new catalysts, e.g., barium cerate [25], metal oxide mixed Sr-Ce [26], Li/NaY (zeolite) [27], lithium-based chicken bone (Li-Cb) composite [28], calcium oxide from seashell [29], graphene [30] or solid acidic ionic liquid polymer [31] have been included in kinetic studies. Despite the potential of microalgal and microbial oil for biodiesel production [8], kinetic models have not yet been developed.
After a thorough analysis of previous studies about the kinetics of vegetable oil esterification and transesterification [12,13,14,15,16,17,18,19,20,21], it may be noticed that proposed mathematical models are generally similar and based on chemical stoichiometry, providing mass variation along reaction time. The input variable is the reaction rate constant, which depends on reaction conditions. Therefore, to work with these models, knowledge of these constants to further validate the model is mandatory. Current studies include kinetic studies on the synthesis of fuel additives from glycerol, biodiesel/hexanol blends or by means of neural networks [32,33,34].
The present work aims to provide a universal mathematical model, to simulate a two-step homogeneously catalyzed biodiesel reaction, that may be of application to any vegetable oil and reaction conditions. Moreover, intermediate conversion steps will be provided, besides any concentration variation of reactants and products. The present study is based on the characterization of reaction kinetics, while the objective is to show, through mathematical modeling, the influence of the most important variables affecting the reaction kinetics of biodiesel production, which are collected into constants of reaction rate. In this sense, models are an indisputable tool to both assist decision-making in the production process of biodiesel and predict reaction yield. However, to collect information about experiment conditions, previous laboratory tests must be carried out. To the best of our knowledge, a universal kinetic model that may be applied to any vegetable oil and reaction conditions is missing. Moreover, catalyst concentration variations are seldom mentioned, thus making difficult the comparison between models or model validation with data from literature.

2. Materials and Methods

2.1. Analysis of Process Variables

In this work, a kinetic model based on reaction stoichiometry is proposed. One mole of triacylglycerol (TAG) and three moles of alcohol (A) react to produce three moles of fatty acid esters (FAE) and 1 mole of glycerol (Gly). The general reaction is named transesterification.
To reduce the order of the model equation, while describing the intermediate steps, the previous equation is divided into three consecutive reversible reactions. At each step, one mole of FAE, diacylglycerols (DAG), monoacylglycerols (MAG) and Gly is released.
In the case of a raw material with a large presence of FFA, a previous esterification reaction is needed. In this situation, one mole of FFA reacts with one mole of A to produce one mole of FAE, plus water (W).
While the catalysts used for esterification are acids (H+), the preferred transesterification catalysts are alkalis (OH-). In fact, during biodiesel production, esterification is carried out under acidic medium (H2SO4) while transesterification usually takes place in a basic medium (NaOH or KOH); both processes being completely independent. Moreover, both reactions cannot take place simultaneously. This ends up with two processes that cannot meet either physically or mathematically.
The reaction rate varies proportionally with the involved species concentrations; the proportionality constant is known as the speed constant, ki. This constant depends on the type of oil, alcohol type, temperature, presence and type of catalyst. It is important to note that the rate constant, k, only depends on the catalyst, whatever the concentration may be.
Although stirring speed is usually not considered among variables, it should be included among parameters affecting processing. Variations of the amount of oil, alcohol, catalyst and temperature affect the reaction rate, although the rate constant remains unchanging.
Previous kinetic studies proposed that oil and alcohol type, alcohol-to-oil molar ratio, catalyst type and concentration, and reactor stirring speed were parameters to take into consideration [12,13,14,15,16,17,18,19,20]. Then, performing the reaction at different temperatures provided information to further calculate activation energy and, thus, rate constants.
A review of the literature shows that, when mathematical models are proposed and validated with experimental data, the catalyst concentration is seldom or never included in the reaction rate equation. To gain knowledge about the implications of this decision, and considering first order (with respect to OH) transesterification reaction, it may be found that by doubling the catalyst concentration, reaction rate would also be doubled. That is, catalyst concentration may increase or decrease reaction rate and, therefore, reaction time, which is really important in the production of biodiesel. Furthermore, as the rate constant has usually been considered independent of reagent concentration, then, alcohol-to-oil molar ratio should also be independent. However, a dependence between rate constant and alcohol-to-oil molar ratio may be found, as the last one depends on catalyst type and concentration [11].

2.2. Mathematical Analysis

Using previous equations and in agreement with the above conditions and hypothesis, a series of expressions based on mass balance and their derivatives with respect to time are proposed. These expressions allow the inclusion of the change in concentration of both reagents and final products with time. These expressions can be represented by the following differential Equations (1):
dFFA dt = k 1 [ FFA ] [ A ] [ H + ] + k 2 [ FAE ] [ W ] [ H + ] dA dt = k 1 [ FFA ] [ A ] [ H + ] + k 2 [ FAE ] [ W ] [ H + ] k 3 [ TAG ] [ A ] [ OH ] + k 4 [ FAE ] [ DAG ] [ OH ] k 5 [ DAG ] [ A ] [ OH ] + k 6 [ FAE ] [ MAG ] [ OH ] k 7 [ MAG ] [ A ] [ OH ] + k 8 [ FAE ] [ Gly ] [ OH ] dFAE dt = k 1 [ FFA ] [ A ] [ H + ] k 2 [ FAE ] [ W ] [ H + ] + k 3 [ TAG ] [ A ] [ OH ] k 4 [ FAE ] [ DAG ] [ OH ] + k 5 [ DAG ] [ A ] [ OH ] k 6 [ FAE ] [ MAG ] [ OH ] + k 7 [ MAG ] [ A ] [ OH ] k 8 [ FAE ] [ Gly ] [ OH ] dW dt = k 1 [ FFA ] [ A ] [ H + ] k 2 [ FAE ] [ W ] [ H + ] dTAG dt = k 3 [ TAG ] [ A ] + k 4 [ FAE ] [ DAG ] [ OH ] dDAG dt = k 3 [ TAG ] [ A ] [ OH ] k 4 [ FAE ] [ DAG ] [ OH ] k 5 [ DAG ] [ A ] [ OH ] + k 6 [ FAE ] [ MAG ] [ OH ] dMAG dt = k 5 [ DAG ] [ A ] [ OH ] k 6 [ FAE ] [ MAG ] [ OH ] k 7 [ MAG ] [ A ] [ OH ] + k 8 [ FAE ] [ Gly ] [ OH ] dGly dt = k 7 [ MAG ] [ A ] [ OH ] k 8 [ FAE ] [ Gly ] [ OH ] dH + dt = k 1 [ FFA ] [ A ] [ H + ] + k 2 [ FAE ] [ W ] [ H + ] dOH dt = k 3 [ TAG ] [ A ] [ OH ] + k 4 [ FAE ] [ DAG ] [ OH ] k 5 [ DAG ] [ A ] [ OH ] + k 6 [ FAE ] [ MAG ] [ OH ] k 7 [ MAG ] [ A ] [ OH ] + k 8 [ FAE ] [ Gly ] [ OH ]
where, k1, k3, k5 and k7 are reaction effective rate constants and k2, k4, k6 and k8 are reverse reaction constants. The chemical elements located between brackets show the concentration of the different species.

2.3. Mathematical Model

To solve the set of differential Equation (1) and mathematical calculations, the software MATLAB 5.3 Academic Use (The MathWorks Inc., Natick, MA, USA) was used and toolbox “Ordinary Differential Equations” (ODE) solver was selected. For this purpose, allocation of states representing Equation (1) was performed. This tool allows ordinary differential equations to be converted into a system of equations of state, where each state corresponds to the concentration of each reagent or product involved in the kinetics. Thus, the model is transformed into a third-degree equation system (Expression (2)), depicting the same number of equations and unknown variables.
k1 x1 x2 x9 + k2 x3 x4 x9 = 0;
k1 x1 x2 x9 + k2 x3 x4 x9k3 x2 x5 x10 + k4 x3 x6 x10k5 x2 x6 x10 + k6 x3 x7 x10k7 x2 x7 x10 + k8 x3 x8 x10 = 0;
k1 x1 x2 x9k2 x3 x4 x9 + k3 x2 x5 x10k4 x3 x6 x10 + k5 x2 x6 x10k6 x3 x7 x10 + k7 x2 x7 x10k8 x3 x8 x10 = 0;
k1 x1 x2 x9k2 x3 x4 x9 = 0;
k3 x2 x5 x10 + k4 x3 x6 x10 = 0;
k3 x2 x5 x10k4 x3 x6 x10k5 x2 x6 x10 + k6 x3 x7 x10 = 0;
k5 x2 x6 x10k6 x3 x7 x10k7 x2 x7 x10 + k8 x3 x8 x10 = 0;
k7 x2 x7 x10k8 x3 x8 x10 = 0;
k1 x1 x2 x9 + k2 x3 x4 x9 = 0;
k3 x2 x5 x10 + k4 x3 x6 x10k5 x2 x6 x10 + k6 x3 x7 x10k7 x2 x7 x10
+k8 x3 x8 x10 = 0;
where allocation of states corresponds to Equation (3):
[FFA] = x1 [DAG] = x6
[A] = x2 [MAG] = x7
[FAE] = x3 [Gly] = x8
[W] = x4 [H+] = x9
[TAG] = x5  [OH] = x10
To simulate the previous model, a rigid system solver (ODE15 and ODE45, included in ODE solver, from MATLAB 5.3 Academic Use) based on a numerical differentiation formula of variable order and multiple width or step was used. ODE15 function is recommended to solve rigid differential equations, while ODE45 function is time consuming due to the rigidity of ordinary differential equations. This rigidity is due to the presence of some terms in equations that can lead to a rapid and unstable variation of solution, unless an extremely small step size is selected. In general, function ODE45 is more accurate, although it consumes more computing time to perform numerical integration. To develop the mathematical model, experimental data from literature were used.

3. Results and Discussion

The proposed model may be positively validated with any previous work from literature. Results from proposed model simulation, with data from literature, are shown below.

3.1. Rate Constants

As mentioned above, as a first step prior to experimentation, rate constant determination was required. In this work, a set of experimental data from different authors corresponding to rate constants, depending on experimental conditions (the most important variables affecting kinetics of biodiesel production) were studied. A small collection of rate constants for different vegetable oils with their respective variables is provided as Supplementary material. These constants are nominated as k1, k2, k3, k4, k5, k6, k7 and k8, following the criteria established in the proposed model and constitute a small database for virtual experimentation with this kind of mathematical model.

3.2. Concentration of Reactants and Products

To validate the model, three experiments provided by different authors from the literature, with different conditions (experiment no. 1 [15], experiment no. 2 [16] and experiment no. 3 [12]), have been analyzed and applied to the proposed model.

3.3. Model Validation with Experiment No. 1

Kinetics of sunflower oil methanolysis at 60 °C, 400 rpm stirring speed and a catalyst concentration (NaOH) of 0.5% (w/w) are shown in Table 1 [15]. For model validation purposes, simulation results from this experiment, together with experimental data from Table 1, are shown in Figure 1. As may be seen, the model suited experimental data, providing information that corresponded to the simulated reality (see Section 3.7 for statistical analysis). Figure 1 exhibits the concentration vs. time of each compound, showing that reaction is completed after the first half-hour.
Table 2 shows rate constants used in the simulation, according to data from Table 1 (complete data are included as Supplementary Material).

3.4. Validation of the Model with Experiment No. 2

In this section, the aims were both to validate the model with other experimental data and to check the effect of the catalyst concentration variation in the proposed model. In this sense, Table 3 shows kinetics from methanolysis with sunflower oil at 25 °C, 600 rpm stirring speed and a catalyst concentration (KOH) of 1.5% (w/w), according to Vicente et al. [16].
For model validation purposes, simulation results are shown in Figure 2, together with experimental data shown in Table 3. It may be inferred that the model appears to fit to the experimental data from experiment no. 2, provided by [16], thus validating the proposed model (see Section 3.7 for statistical analysis).
Table 4 shows the effect of catalyst concentration (KOH) in the production of sunflower oil biodiesel, using methanol at 25 °C and 600 rpm stirring speed, according to results from [16]. Table 5 shows the rate constants used for simulation. Simulation results and experimental data (from Table 4) are shown in Figure 3. It may be inferred that the model appears to fit to experimental data from experiment no. 2, provided by [16], thus validating the proposed model (see Section 3.7 for statistical analysis). Figure 3 shows the effect of catalyst concentration variation included in the proposed model.

3.5. Validation of the Model with Experiment No. 3

Table 6 shows results from M. L. Pisarello et al. [12] after sunflower oil esterification with ethanol at 70 °C and a concentration of catalyst (H2SO4) of 0.2% (v/v). On the other hand, Table 7 shows rate constants used for simulation, according to data conditions provided in Table 6.
For model validation purposes, simulation results from this experiment are shown in Figure 4, together with experimental data from Table 6. It may be seen that the theoretical data calculated from the model (representing concentration variation of W and FFA with time) are in good agreement with experimental data from [12]. However, an exception is provided by fatty acid ethyl esters (FAEE), for which the comparison between last simulated and experimental values seem to differ: this may be explained by an experimental calculation error (see Section 3.7 for statistical analysis).

3.6. Simulation to Analyze the Effect of Catalyst Concentration

Considering the proposed model, change in concentration of transesterification reactants and products, after 1, 5 and 10 min of reaction time, have been simulated. Results are shown in Figure 5, while simulation conditions are listed in Supplementary material, considering an initial concentration of methanol and triacylglycerols of [A] = 12 mol/L and [TAG] = 2 mol/L, respectively. As expected, simulation shows that, for the same time instants, the higher the amount of catalyst, the higher the reaction rate. It may also be inferred that, for a catalyst concentration around 1% (w/w), the conversion of [TAG] to products ([FAME], [DAG], [MAG] and [Gly]) is faster, which leads to conclude that, for this production, concentration of catalyst is optimal.

3.7. Statistical Analysis of the Model

A summary of the statistics of the experiments is shown in Table 8 and Table 9. As may be seen from Table 9, model and experimental values are in good agreement. The model goodness of fit has been calculated by means of the coefficient of determination (R-squared) that shows the proportion of variation of results that can be explained by the model. In other words, a value of 1.0 indicates a highly reliable model, while a value close to 0.0 means a failing model. In this sense, the correlation between experimental data and simulated values (Figure 1, Figure 2, Figure 3 and Figure 4) results in 0.92184506, 0.99999914, 0.99232725 and 0.98553624, respectively. As may be seen, R2 value is, in all cases, above 92%, thus showing the goodness and reliability of the model.

4. Conclusions

Mathematical models are versatile tools for a low-cost study, reducing raw material consumption and experimentation instrument needs, among others. In this sense, models for biodiesel production are increasingly widespread. Models found in literature consider the use of a solution of catalyst and alcohol; therefore, experiments with separated reagent concentration variations are not modeled. In this work, as a novelty, the proposed model considers the presence and variations of homogenous catalysts and alcohol, separately, allowing further simulations.
Comparing studies from literature with the proposed model, it may be seen that the reaction depends on a number of variables, including the type of vegetable oil, reaction stirring speed, type and concentration of catalyst, type of alcohol, alcohol-to-oil molar ratio, reaction time and temperature. These parameters are crucial and differentiating elements in the transformation of oils into biodiesel.
To increase the usefulness of kinetic models from literature, a universal model of biodiesel production that allows implementing control strategies is provided. In this sense, the model could provide information about best reaction time to achieve maximum biodiesel production, optimum catalyst concentration for different reaction times, etc., allowing automated production, thus reducing production costs.
The biggest drawback presented by previous models found in literature is the need, before starting simulation, of the rate constant determination, which means that experimental work is needed. For this reason, a collection of rate constants, derived from different experiments using different vegetable oils and conditions is provided as Supplementary material. This collection allows simulation of a potential production plant, thus helping to find out production deviations.

Supplementary Materials

The following are available online at https://www.mdpi.com/1996-1073/13/11/2994/s1, Table S1: Reaction inputs (oil, reaction type and conditions) to determine effective rate constants, Table S2: Effective rate constants resulting from esterification (SFO: sunflower oil; CO: coconut oil; PoO: pongamia oil; ERO: Ethiopian rape oil; PO: palm oil; SO: soybean oil; RO: rapeseed oil; Me: methanol; Et: ethanol; n.d.: no data).

Author Contributions

Conceptualization, M.P.D. and F.J.L.-G.; software, A.A.; validation, A.A.; formal analysis, M.P.D., F.J.L.-G. and A.A.; writing—original draft preparation, A.A.; writing—review and editing, M.P.D.; supervision, M.P.D.; project administration, M.P.D.; funding acquisition, M.P.D. All authors have read and agreed to the published version of the manuscript.

Funding

Authors wish to acknowledge the support of the Spanish Ministry of Economy and Competitiveness (grant no. ENE2013-47769R).

Acknowledgments

Authors are thankful to David Leiva-Candia for great statistical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, H.S.; Li, S.Z. Optimization of continuous solid-state distillation process for cost-effective bioethanol production. Energies 2020, 13, 854. [Google Scholar] [CrossRef] [Green Version]
  2. Mahmud, N.; Rosentrater, K.A. Techno-economic analysis (TEA) of different pretreatment and product separation technologies for cellulosic butanol production from oil palm frond. Energies 2020, 13, 181. [Google Scholar] [CrossRef] [Green Version]
  3. Al-Addous, M.; Alnaief, M.; Class, C.; Nsair, A.; Kuchta, K.; Alkasrawi, M. Technical possibilities of biogas production from olive and date waste in Jordan. BioResources 2017, 12, 9383–9395. [Google Scholar]
  4. Zahan, K.A.; Kano, M. Biodiesel production from palm oil, its by-products, and mill effluent: A review. Energies 2018, 11, 2132. [Google Scholar] [CrossRef] [Green Version]
  5. Filip, O.; Janda, K.; Kristoufek, L.; Zilberman, D. Food versus fuel: An updated and expanded evidence. Energy Econ. 2019, 82, 152–166. [Google Scholar] [CrossRef] [Green Version]
  6. Shrestha, D.S.; Staab, B.D.; Duffleld, J.A. Biofuel impact on food prices index and land use change. Biomass Bioenergy 2019, 124, 43–53. [Google Scholar] [CrossRef]
  7. Gu, Y.Q.; Cai, F.F.; Zhu, Z.; Dai, Z.Q.; Chen, C.; Liu, G.Q. Improving the methane production from zucchini stem by response surface methodology and different pretreatments. Ind. Crop. Prod. 2020, 150, 112402. [Google Scholar] [CrossRef]
  8. Leiva-Candia, D.E.; Pinzi, S.; Redel-Macias, M.D.; Koutinas, A.; Webb, C.; Dorado, M.P. The potential for agro-industrial waste utilization using oleaginous yeast for the production of biodiesel. Fuel 2014, 123, 33–42. [Google Scholar] [CrossRef]
  9. Dorado, M.P.; Cruz, F.; Palomar, J.M.; Lopez, F.J. An approach to the economics of two vegetable oil-based biofuels in Spain. Renew. Energy 2006, 31, 1231–1237. [Google Scholar] [CrossRef]
  10. Van Gerpen, J. Biodiesel processing and production. Fuel Process. Technol. 2005, 86, 1097–1107. [Google Scholar] [CrossRef]
  11. Dorado, M.P.; Ballesteros, E.; Mittelbach, M.; Lopez, F.J. Kinetic parameters affecting the alkali-catalyzed transesterification process of used olive oil. Energy Fuels 2004, 18, 1457–1462. [Google Scholar] [CrossRef]
  12. Pisarello, M.L.; Dalla Costa, B.; Mendow, G.; Querini, C.A. Esterification with ethanol to produce biodiesel from high acidity raw materials. Kinetic studies and analysis of secondary reactions. Fuel Process. Technol. 2010, 91, 1005–1014. [Google Scholar] [CrossRef]
  13. Berrios, M.; Siles, J.; Martin, M.A.; Martin, A. A kinetic study of the esterification of free fatty acids (FFA) in sunflower oil. Fuel 2007, 86, 2383–2388. [Google Scholar] [CrossRef]
  14. Kadi, M.A.; Akkouche, N.; Awad, S.; Loubar, K.; Tazerout, M. Kinetic study of transesterification using particle swarm optimization method. Heliyon 2019, 5, e02146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bambase, M.E.; Nakamura, N.; Tanaka, J.; Matsumura, M. Kinetics of hydroxide-catalyzed methanolysis of crude sunflower oil for the production of fuel-grade methyl esters. J. Chem. Technol. Biotechnol. 2007, 82, 273–280. [Google Scholar] [CrossRef]
  16. Vicente, G.; Martinez, M.; Aracil, J.; Esteban, A. Kinetics of sunflower oil methanolysis. Ind. Eng. Chem. Res. 2005, 44, 5447–5454. [Google Scholar] [CrossRef]
  17. Karmee, S.K.; Chandna, D.; Ravi, R.; Chadha, A. Kinetics of base-catalyzed transesterification of triglycerides from Pongamia oil. J. Am. Oil Chem. Soc. 2006, 83, 873–877. [Google Scholar] [CrossRef]
  18. Vicente, G.; Martinez, M.; Aracil, J. Kinetics of Brassica carinata oil methanolysis. Energy Fuels 2006, 20, 1722–1726. [Google Scholar] [CrossRef]
  19. Narvaez, P.C.; Rincon, S.M.; Sanchez, F.J. Kinetics of palm oil methanolysis. J. Am. Oil Chem. Soc. 2007, 84, 971–977. [Google Scholar] [CrossRef]
  20. Noureddini, H.; Zhu, D. Kinetics of transesterification of soybean oil. J. Am. Oil Chem. Soc. 1997, 74, 1457–1463. [Google Scholar] [CrossRef]
  21. Bashiri, H.; Pourbeiram, N. Biodiesel production through transesterification of soybean oil: A kinetic Monte Carlo study. J. Mol. Liq. 2016, 223, 10–15. [Google Scholar] [CrossRef]
  22. Fedosov, S.N.; Brask, J.; Pedersen, A.K.; Nordblad, M.; Woodley, J.M.; Xu, X.B. Kinetic model of biodiesel production using immobilized lipase Candida antarctica lipase B. J. Mol. Catal. B 2013, 85–86, 156–168. [Google Scholar] [CrossRef]
  23. Firdaus, M.Y.; Guo, Z.; Fedosov, S.N. Development of kinetic model for biodiesel production using liquid lipase as a biocatalyst, esterification step. Biochem. Eng. J. 2016, 105, 52–61. [Google Scholar] [CrossRef]
  24. Firdaus, M.Y.; Brask, J.; Nielsen, P.M.; Guo, Z.; Fedosov, S. Kinetic model of biodiesel production catalyzed by free liquid lipase from Thermomyces lanuginosus. J. Mol. Catal. B 2016, 133, 55–64. [Google Scholar] [CrossRef]
  25. Sahani, S.; Roy, T.; Sharma, Y.C. Clean and efficient production of biodiesel using barium cerate as a heterogeneous catalyst for the biodiesel production; kinetics and thermodynamic study. J. Clean. Prod. 2019, 237, 117699. [Google Scholar] [CrossRef]
  26. Banerjee, S.; Sahani, S.; Sharma, Y.C. Process dynamic investigations and emission analyses of biodiesel produced using Sr-Ce mixed metal oxide heterogeneous catalyst. J. Environ. Manag. 2019, 248, 109218. [Google Scholar] [CrossRef] [PubMed]
  27. Li, Z.; Ding, S.X.; Chen, C.; Qu, S.K.; Du, L.X.; Lu, J.; Ding, J.C. Recyclable Li/NaY zeolite as a heterogeneous alkaline catalyst for biodiesel production: Process optimization and kinetics study. Energy Convers. Manag. 2019, 192, 335–345. [Google Scholar] [CrossRef]
  28. AlSharifi, M.; Znad, H. Development of a lithium based chicken bone (Li-Cb) composite as an efficient catalyst for biodiesel production. Renew. Energy 2019, 136, 856–864. [Google Scholar] [CrossRef]
  29. Sun, Y.; Zhang, J.P.; Sun, Z.; Zhang, L. Biodiesel production using calcium-based catalyst from venus shell: Modeling of startup production in an industrial reactor. Environ. Prog. Sustain. Energy 2019, 38, e13053. [Google Scholar] [CrossRef]
  30. Borah, M.J.; Sarmah, H.J.; Bhuyan, N.; Mohanta, D.; Deka, D. Application of Box-Behnken design in optimization of biodiesel yield using WO3/graphene quantum dot (GQD) system and its kinetics analysis. Biomass Conv. Bioref. 2020. [Google Scholar] [CrossRef]
  31. Panchal, B.; Chang, T.; Qin, S.J.; Sun, Y.Z.; Wang, J.X.; Bian, K. Optimization and kinetics of tung nut oil transesterification with methanol using novel solid acidic ionic liquid polymer as catalyst for methyl ester synthesis. Renew. Energy 2020, 151, 796–804. [Google Scholar] [CrossRef]
  32. Kulkarni, R.M.; Britto, P.J.; Narula, A.; Saqline, S.; Anand, D.; Bhagyalakshmi, C.; Herle, R.N. Kinetic studies on the synthesis of fuel additives from glycerol using CeO2-ZrO2 metal oxide catalyst. Biofuel Res. J. 2020, 7, 1100–1108. [Google Scholar] [CrossRef] [Green Version]
  33. Thomas, J.J.; Manojkumar, C.V.; Sabu, V.R.; Nagarajan, G. Development and validation of a reduced chemical kinetic model for used vegetable oil biodiesel/1-hexanol blend for engine application. Fuel 2020, 273, 117780. [Google Scholar] [CrossRef]
  34. Galvan, D.; Cremasco, H.; Mantovani, A.C.G.; Bona, E.; Killner, M.; Borsato, D. Kinetic study of the transesterification reaction by artificial neural networks and parametric particle swarm optimization. Fuel 2020, 267, 117221. [Google Scholar] [CrossRef]
Figure 1. Simulation and validation of proposed model with experimental data from [15]; A: alcohol concentration, FAME: fatty acid methyl ester concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, Gly: glycerol concentration.
Figure 1. Simulation and validation of proposed model with experimental data from [15]; A: alcohol concentration, FAME: fatty acid methyl ester concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, Gly: glycerol concentration.
Energies 13 02994 g001
Figure 2. Simulation and validation with experimental data from [16]; A: alcohol concentration, FAME: fatty acid methyl ester concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, Gly: glycerol concentration.
Figure 2. Simulation and validation with experimental data from [16]; A: alcohol concentration, FAME: fatty acid methyl ester concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, Gly: glycerol concentration.
Energies 13 02994 g002
Figure 3. Simulation and model validation with experimental data from [16]; FAME: fatty acid methyl ester concentration.
Figure 3. Simulation and model validation with experimental data from [16]; FAME: fatty acid methyl ester concentration.
Energies 13 02994 g003
Figure 4. Simulation and validation of the proposed model with experimental data from [12]; W: water; FFA: free fatty acid concentration; FAEE: fatty acid ethyl ester concentration.
Figure 4. Simulation and validation of the proposed model with experimental data from [12]; W: water; FFA: free fatty acid concentration; FAEE: fatty acid ethyl ester concentration.
Energies 13 02994 g004
Figure 5. Simulation of catalyst concentration for production of fatty acid methyl ester (FAME), diacylglycerol concentration (DAG), monoacylglycerol concentration (MAG) and glycerol concentration (Gly) at 60 s, 300 s and 600 s of reaction time; TAG: triacylglycerol concentration, A: alcohol concentration.
Figure 5. Simulation of catalyst concentration for production of fatty acid methyl ester (FAME), diacylglycerol concentration (DAG), monoacylglycerol concentration (MAG) and glycerol concentration (Gly) at 60 s, 300 s and 600 s of reaction time; TAG: triacylglycerol concentration, A: alcohol concentration.
Energies 13 02994 g005
Table 1. Kinetics of sunflower oil methanolysis at 60 °C, 400 rpm stirring speed and a catalyst concentration (NaOH) of 0.5% (w/w) [15]; t: time, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAME: fatty acid methyl ester concentration, Gly: glycerol concentration, A: alcohol concentration.
Table 1. Kinetics of sunflower oil methanolysis at 60 °C, 400 rpm stirring speed and a catalyst concentration (NaOH) of 0.5% (w/w) [15]; t: time, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAME: fatty acid methyl ester concentration, Gly: glycerol concentration, A: alcohol concentration.
t (min)A (mol L−1)TAG (mol L−1)DAG (mol L−1)MAG (mol L−1)FAME (mol L−1)Gly (mol L−1)
0610000
13.190.720.030.020.080.17
22.950.550.050.030.440.23
32.680.280.0890.0851.680.53
42.600.210.080.081.910.55
62.550.180.080.082.220.62
102.450.170.070.062.280.68
202.440.160.060.062.350.69
302.440.150.040.042.380.72
452.430.140.030.032.410.73
602.420.130.010.012.40.73
752.420.110.010.012.410.73
902.420.11002.40.73
1052.420.1100.012.40.72
1202.420.11002.40.73
Table 2. Reaction conditions and rate constants used from sunflower oil biodiesel simulation, according to data conditions from [15].
Table 2. Reaction conditions and rate constants used from sunflower oil biodiesel simulation, according to data conditions from [15].
ParametersValues
Temperature, T (°C)60
CatalystNaOH
Catalyst concentration (%, w/w)0.5
Methanol-to-oil (molar ratio)6:1
Stirring speed (rpm)400
k3 (L mol−1·s−1)0.00149167
k4 (L mol−1·s−1)0.00015667
k5 (L mol−1·s−1)0.00580000
k6 (L mol−1·s−1)0.00214167
k7 (L mol−1·s−1)0.00814000
k8 (L mol−1·s−1)0.00064167
Table 3. Kinetics of sunflower oil methanolysis at 25 °C, 600 rpm stirring speed and a catalyst concentration (KOH) of 1.5% (w/w) [16]; t: time, A: methanol concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAME: fatty acid methyl ester concentration, Gly: glycerol concentration.
Table 3. Kinetics of sunflower oil methanolysis at 25 °C, 600 rpm stirring speed and a catalyst concentration (KOH) of 1.5% (w/w) [16]; t: time, A: methanol concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAME: fatty acid methyl ester concentration, Gly: glycerol concentration.
t (min)A (mol L−1)TAG (mol L−1)DAG (mol L−1)MAG (mol L−1)FAME (mol L−1)Gly (mol L−1)
051.620000
141.150.050.030.890.09
23.60.750.110.061.30.2
33.40.440.210.11.550.31
43.30.410.090.041.70.39
53.20.390.080.041.780.49
63.10.360.080.041.850.55
830.330.070.031.950.62
102.90.310.050.022.050.69
122.830.270.040.022.120.71
152.790.230.040.022.180.74
202.710.20.030.012.250.76
252.650.180.020.012.350.77
302.60.180.020.012.390.79
452.580.180.0102.430.81
602.570.170.0102.480.81
902.550.17002.490.81
1202.550.17002.490.81
Table 4. Effect of catalyst concentration (KOH) in the production of sunflower oil fatty oil methyl ester (mol L−1) and methanol (25 °C y 600 rpm) [16].
Table 4. Effect of catalyst concentration (KOH) in the production of sunflower oil fatty oil methyl ester (mol L−1) and methanol (25 °C y 600 rpm) [16].
Time (min)Catalyst Concentration % (w/w)
0.511.5
0000
10.6131.1221.234
21.3951.3881.605
31.6511.7491.8552
41.7711.9011.955
61.8831.9742.181
7.51.9632.0722.255
102.0012.112.329
122.0172.1542.357
152.0512.222.381
202.0622.2362.422
302.0712.2412.447
452.0822.2562.469
602.0912.2452.472
902.0972.2512.488
1202.1262.2582.489
Table 5. Reaction conditions and rate constants used in sunflower oil biodiesel simulation, according to conditions provided by [16].
Table 5. Reaction conditions and rate constants used in sunflower oil biodiesel simulation, according to conditions provided by [16].
ParametersValues
Temperature, T (°C)252525
CatalystKOHKOHKOH
Catalyst concentration (%, w/w)0.51.01.5
Methanol-to-oil (molar ratio)6:16:16:1
Stirring speed (rpm)600600600
k3 (L mol−1 s−1)0.014583330.020964360.02591036
k4 (L mol−1 s−1)0.052083330.102725370.12955182
k5 (L mol−1 s−1)0.031250000.175052410.26260504
k6 (L mol−1 s−1)0.029166670.228511530.30462185
k7 (L mol−1 s−1)0.045833330.028301890.02240896
k8 (L mol−1 s−1)0.003333330.011530400.00053922
Table 6. Free fatty acids (FFA), fatty acid ethyl ester (FAEE) and water (W) concentrations in the reaction of sunflower oil with ethanol at 70 °C and a concentration of catalyst (H2SO4) of 0.2% (v/v) [12]; t: time.
Table 6. Free fatty acids (FFA), fatty acid ethyl ester (FAEE) and water (W) concentrations in the reaction of sunflower oil with ethanol at 70 °C and a concentration of catalyst (H2SO4) of 0.2% (v/v) [12]; t: time.
t (min)FFA (mol·L−1)FAEE (mol·L−1)W (mol·L−1)
00.3800
100.220.140.34
600.110.290.50
1200.060.380.51
Table 7. Reaction conditions and rate constants used in the simulation according to data provided by [12]; n.d.: no data.
Table 7. Reaction conditions and rate constants used in the simulation according to data provided by [12]; n.d.: no data.
ParametersValues
Temperature, T (°C)70
CatalystH2SO4
Catalyst concentration (%, w/w)0.2
Ethanol-to-oil (molar ratio)15:1
Stirring speed (rpm)n.d.
k1 (L mol−1·s−1)0.00016667
k2 (L mol−1·s−1)0.00041667
Table 8. Statistical analysis of simulation of experiments no. 1, 2 and 3; A: alcohol concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAEE: fatty acid ethyl ester concentration, Gly: glycerol concentration; FFA: free fatty acid concentration; W: water content.
Table 8. Statistical analysis of simulation of experiments no. 1, 2 and 3; A: alcohol concentration, TAG: triacylglycerol concentration, DAG: diacylglycerol concentration, MAG: monoacylglycerol concentration, FAEE: fatty acid ethyl ester concentration, Gly: glycerol concentration; FFA: free fatty acid concentration; W: water content.
Simulation of Experiment No. 1Simulation of Experiment No. 2Simulation of Experiment No. 3
Standard DeviationVarianceStandard DeviationVarianceStandard DeviationVariance
A (mol·L−1) 0.03613720.0013059
TAG (mol·L−1)0.02737740.00079980.02262730.0005120
DAG (mol·L−1)0.00251420.00000740.00273230.0000075
MAG (mol·L−1)0.00312500.00001210.00309660.0000096
FAEE (mol·L−1)0.04390170.00195470.01770230.00031340.03705760.0013733
Gly (mol·L−1)0.00265130.00001000.00494640.0000245
FFA (mol·L−1) 0.36928910.1363744
W (mol·L−1) 0.00285000.0000081
Table 9. Statistical analysis of experiment no. 2 (including the effect of catalyst concentration, from 0.5 to 1.5%, w/w); FAEE: fatty acid ethyl ester concentration.
Table 9. Statistical analysis of experiment no. 2 (including the effect of catalyst concentration, from 0.5 to 1.5%, w/w); FAEE: fatty acid ethyl ester concentration.
FAEE (mol·L−1)
Standard DeviationVariance
Catalyst concentration, H2SO4(%, w/w)0.50.00154110.0000024
10.00390860.0000153
1.50.00037560.0000001

Share and Cite

MDPI and ACS Style

Alcantara, A.; Lopez-Gimenez, F.J.; Dorado, M.P. Universal Kinetic Model to Simulate Two-Step Biodiesel Production from Vegetable Oil. Energies 2020, 13, 2994. https://doi.org/10.3390/en13112994

AMA Style

Alcantara A, Lopez-Gimenez FJ, Dorado MP. Universal Kinetic Model to Simulate Two-Step Biodiesel Production from Vegetable Oil. Energies. 2020; 13(11):2994. https://doi.org/10.3390/en13112994

Chicago/Turabian Style

Alcantara, A., F. J. Lopez-Gimenez, and M. P. Dorado. 2020. "Universal Kinetic Model to Simulate Two-Step Biodiesel Production from Vegetable Oil" Energies 13, no. 11: 2994. https://doi.org/10.3390/en13112994

APA Style

Alcantara, A., Lopez-Gimenez, F. J., & Dorado, M. P. (2020). Universal Kinetic Model to Simulate Two-Step Biodiesel Production from Vegetable Oil. Energies, 13(11), 2994. https://doi.org/10.3390/en13112994

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop