Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 601374, 6 pages

http://dx.doi.org/10.1155/2014/601374

Research Article

Numerical Analysis for the Synthesis of Biodiesel Using Spectral Relaxation Method

Z. G. Makukula,¹ S. S. Motsa,¹ and S. Shateyi²

¹School of Mathematics, Statistics and Computer Sciences, University of KwaZulu-Natal, Private Bag X01, Scottsville, Pietermaritzburg 3209, South Africa

²Department of Mathematics, University of Venda, Private Bag X5050, Thohoyandou 0950, South Africa

Received 28 March 2013; Accepted 6 November 2013; Published 23 January 2014

Academic Editor: Mufid Abudiab

Copyright © 2014 Z. G. Makukula et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Biodiesel is an alternative diesel fuel chemically defined as the mono-alkyl esters of long chain fatty acids derived from vegetable oils or animal fat. It is becoming more attractive as an alternative fuel due to the depleting fossil fuel resources. A mathematical model for the synthesis of biodiesel from vegetable oils and animal fats is presented in this study. Numerical solutions of the model are found using a spectral relaxation method. The method, originally developed for boundary value problems, is an iterative scheme based on the Chebyshev spectral collocation method developed by decoupling systems of equations using Gauss-Seidel type of techniques. The effects of the reaction rate constants and initial concentrations of the reactants on the amount of the final product are being investigated. The accuracy of the numerical results is validated by comparison with known analytical results and numerical results obtained using ode45, an efficient explicit 4th and 5th order Runge-Kutta method used to integrate both linear and nonlinear differential equations.

1. Introduction

Due to the continuous uncertainty and increasing environmental impact associated with the utilization of petroleum-based diesel fuel, the demands for alternative fuels have increased considerably in recent years [1]. Methyl and ethyl esters derived from vegetable oil or animal fat, known as biodiesel, have good potential as alternative diesel fuel [2]. Biodiesel is synthesized through a chemical process called transesterification. During this chemical process, triglycerides in vegetable oils or animal fats react with short chain alcohols such as methanol and ethanol in the presence of homogeneous basic catalysts such as sodium hydroxide and potassium hydroxide to form fatty acid methyl esters (biodiesel) and glycerin [3, 4].

Biodiesel has many advantages over petroleum-based diesel fuels. It is biodegradable and nontoxic and produce less particles, smoke, and carbon monoxide [2]. It is also renewable and usable in a variety of diesel engines with minimum or no modification necessary [3, 5]. It burns clean, and do not form engine deposits or generate harmful pollutants which results in a significant reduction of the types of pollutants that contribute to pollution and global warming and releases up to 85% fewer cancer-causing agents [6–10].

Experiments have been carried out to investigate different aspects in the production of biodiesel [3, 4, 11–15]. In this work, we carry out a numerical study of the transesterification process for the synthesis of biodiesel from vegetable oils. We develop a mathematical model from the experimental works of Santos et al. [2], Burnham et al. [6], and Gunvachai et al. [7], where the reaction equations and reaction rates are given. In this study, the chemical reaction equations and rates are reduced to a set of nonlinear first order differential equations that are solved using a spectral relaxation method. The method decouples nonlinear systems of equations using ideas similar to those of the Gauss-Seidel iterative scheme to give rise to a linear system which is solved sequentially using the Chebyshev spectral method. The method's applicability to initial value problems is also explored in this work. The effects of the reaction rate constants and initial concentrations of the reactants on the amount of the final product are being investigated. The accuracy of the numerical results is validated by comparison with known analytical results and numerical results obtained using ode45, an efficient explicit th and th order Runge-Kutta method used to integrate both linear and nonlinear differential equations.

2. Mathematical Formulation

Biodiesel () is produced by the transesterification of large, branched triglyceride () molecules into smaller, straight chain molecules with methanol () in the presence of an alkali or acid as a (nonreactive) catalyst. The process occurs in three consecutive reactions which in turn produce three by-products, di-glyceride (), mono-glyceride (), and glycerol () [2, 4, 6, 7, 15, 16]. The thermal transesterification reactions proposed in [12, 13] may be characterized as

where , , and are reaction rate constants. The alternative mathematical model proposed by Kusdiana and Saka [14] ignores the two intermediate reactions and assumes the reaction

In this study we consider the three-stage model defined by (1)–(3). Following [4, 6], we assume that there is only one biodiesel product produced in all three reactions so that . In [6], an adaptive Runge-Kutta integrator was used to simulate the set of ODEs that represent the reaction network defined by (1)–(3). In this study, we derive a system of differential equations that model the change in the concentrations of the reactants and the product species with time from (1)–(3) by using the notation,

where the square bracket denotes the concentration of the given chemical at time . These equations are

where the prime denotes the time derivative, subject to the initial conditions

In (6), the positive terms in the differential equations are those that contribute to an increase in the th variable , that is to say those terms that contribute positively to the concentration of the different products. The negative terms are those contributing to the decay in the concentration of the different species. Effectively, we start with definite known concentrations of the reacting species with no by-products. As , the concentrations of the reacting species are expected to approach zero and those of the by-products to approach some limiting values,   .

3. Method of Solution

In this section, (6) are solved using the spectral relaxation method (SRM). Principally, the method has been developed for similarity boundary layer problems with at least one the essential profiles such as velocity, temperature, or concentration decaying exponentially. The governing systems of equations are decoupled in a manner parallel to the Gauss-Seidel relaxation method normally used for solving linear algebraic system of equations. Worth noting also is that the method bears some similarities with the Jacobi and Gauss-Seidel waveform relaxation methods earlier developed by Lelarasmee et al. [17] for solutions of initial value problems (IVPs) in very large-scale electrical networks. However, these methods may result into nonlinear differential equations which require further linearization using appropriate linearization schemes while the SRM gives rise into a sequence of linear differential equations which can be directly solved using standard methods. Spectral methods [18, 19] are then used in this work to solve the resulting scheme hence, the name spectral relaxation method. Spectral methods have been chosen in this work because of their well-documented advantages including high accuracy. In developing the sequence of linear equations using the SRM, the following guidelines are used.

(1)The equations have to be arranged in such a way that the equation with the least number of unknowns is placed first on the list.(2)Variables say can be used to assign the ordered equations in [1], such that , is an unknown function in the th equation. It is chosen to be the unknown function associated with the highest order derivative of the th equation.(3)To generate the iterative scheme, in each equation, only terms linear in are to be evaluated at the current iteration level, . All other terms, linear and nonlinear, in the other functions are to be evaluated from the previous iteration, , with an exception from the second equation, where updated solutions for obtained from the previous equation will be used. Nonlinear terms in are also evaluated from the previous iteration. Terms made up of a product of and its derivatives, the derivative terms are to be evaluated from .

In view of the guidelines above, we develop the iterative scheme for (6) as follows:

To solve the iteration schemes (8) we use the Chebyshev spectral collocation method. For brevity, we omit the details of the spectral methods and refer interested readers to [18, 19]. Before applying the spectral method, it is convenient to transform the domain on which the governing equation is defined to the interval on which the spectral method can be implemented. We use the transformation to map the interval to , where is a finite time. The basic idea behind the spectral collocation method is the introduction of a differentiation matrix which is used to approximate the derivatives of the unknown variables at the collocation points as the matrix vector product

where is the number of grid points known as collocation points, , and is the vector function at the collocation points . Applying the Chebyshev spectral collocation method in (8) gives

where

From the above equations, is a diagonal matrix and ,  , are the values of the functions when evaluated at the collocation points. To start the SRM iterative scheme (20), we use the initial conditions (7) as initial solutions. The scheme is repeated until the solution converges to the numerical solution obtained using the ode45.

4. Results and Discussion

In this section we present SRM results of the governing equations (6) subject to the initial conditions (7). In experiments, the contents of , , , and in samples of reaction product are analyzed periodically using, for example, gas chromatography and the mass fractions of each component calculated based on the concentration of some internal standard [4, 20]. However in this work we use concentration-time curves to determine the concentration of each reacting component and products at any given time. We also present results showing the effect of the reaction rate constants on the main product, biodiesel, and the effect of the concentration of the main reactant, methanol, on the concentration of biodiesel produced. The accuracy of the SRM in this work was validated using ode45, a fourth and fifth order Runge-Kutta scheme whose tolerance level was set to be seven decimal digits (). A further comparison with similar results in the literature was made. Figures 1, 2, 3, and 4 show the conversion of the reaction and yield of biodiesel and the other by-products with reaction time. In all figures, the solid line represents the ode45 solution while the filled circles, squares, and diamonds represent the SRM solutions. Convergence of the SRM was reached after the th iteration when in all simulations. Figure 1 shows the variation of with reaction rate constant at different concentrations of . The concentration of was taken to be and while . The concentration of produced shows to be directly proportional to as it increases with increase in . The concentration of produced doubles when is also doubled. The SRM solutions agree with the ode45 results.

(a)

(b)

Figure 1: Concentration-time curves showing the effect of on the concentration of biodiesel produced for different concentrations of . The solid lines represent ode45 solutions while the figures represent SRM solutions for , , .

(a)

(b)

Figure 2: The effect of on the concentration-time curves for biodiesel produced at different concentrations of . The solid lines represent ode45 solutions while the figures represent SRM solutions for , , .

(a)

(b)

Figure 3: Concentration-time curves for at different concentrations of and that of and when , , . The solid lines represent ode45 solutions while the figures represent SRM solutions.

Figure 4: The variation of , , and with time when , , , . The solid lines denote ode45 solutions while the figures represent SRM solutions.

In Figure 2, the effect of is shown to decrease . In the figure, , , and . The concentration of biodiesel produced in this case is shown to be equal to that of the methanol used. A good agreement between the two solutions is observed.

The variation of with the concentration of the main reactant, , is shown in Figure 3. Increasing increases that of produced. In the same figure, it is shown that when ,, and , about of gets used up in the reaction, while methanol gets exhausted after about 120 days. The results here are consistent with the analytical results in our earlier study [21]. The two sets of solutions continue to show a very good agreement.

Concentration-time curves for the three by-products are shown in Figure 4. Di-glyceride is produced in stronger concentrations followed by and lastly .

5. Conclusion

A spectral relaxation method has been used successfully in this work to solve a system of differential equations governing the thermal transesterification reactions in the production of biodiesel from vegetable oils. The method proved useful in solving initial value problems of the type in this study. The results were consistent with numerical results obtained using ode45 and those in the works of Sibanda et al. [21]. In the study we investigated the effect of the reaction rate constants and that of the concentration of methanol as the main reactant on the final concentration of biodiesel produced. For optimal results, the study showed that the first reaction rate constant should be kept large while the value of should be small. The results showed that that is increasing increases biodiesel production while increasing decreases biodiesel production. The amount of biodiesel produced increases with increase in the initial concentration of methanol used. Methanol is used up in the reaction after sometime. Experimental research is costly compared to theoretical studies. It is however essential to use methods that give accurate results for trustworthy findings. The spectral relaxation method promises to be a viable tool for theoretical studies in the biosciences.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

N. U. Soriano Jr., R. Venditti, and D. S. Argyropoulos, “Biodiesel synthesis via homogeneous Lewis acid-catalyzed transesterification,” Fuel, vol. 88, no. 3, pp. 560–565, 2009.

F. F. P. Santos, J. Q. Malveira, M. G. A. Cruz, and F. A. N. Fernandes, “Production of biodiesel by ultrasound assisted esterification of Oreochromis niloticus oil,” Fuel, vol. 89, no. 2, pp. 275–279, 2010.

K. Gunvachai, M. G. Hassan, G. Shama, and K. Hellgardt, “A new solubility model to describe biodiesel formation kinetics,” Process Safety and Environmental Protection B, vol. 85, no. 5, pp. 383–389, 2007.

S. Joelianingsih, H. Maeda, S. Hagiwara et al., “Biodiesel fuels from palm oil via the non-catalytic transesterification in a bubble column reactor at atmospheric pressure: a kinetic study,” Renewable Energy, vol. 33, no. 7, pp. 1629–1636, 2008.

M. Canakci and J. H. Van Gerpen, “Comparison of engine performance and emissions for petroleum diesel fuel, yellow grease biodiesel, and soybean oil biodiesel,” Transactions of the American Society of Agricultural Engineers, vol. 46, no. 4, pp. 937–944, 2003.

S. C. Burnham, D. P. Searson, M. J. Willis, and A. R. Wright, “Inference of chemical reaction networks,” Chemical Engineering Science, vol. 63, no. 4, pp. 862–873, 2008.

K. Gunvachai, M. G. Hassan, G. Shama, and K. Hellgardt, “A new solubility model to describe biodiesel formation kinetics,” Process Safety and Environmental Protection B, vol. 85, no. 5, pp. 383–389, 2007.

K. Murugesan, T. Ganapathy, and R. P. Gakkhar, “Performance optimization of Jatropha biodiesel engine model using Taguchi approach,” Applied Energy, vol. 86, no. 11, pp. 2476–2486, 2009.

P. Pogorevc, B. Kegl, and L. Skerget, “Diesel and biodiesel fuel spray simulations,” Energy and Fuels, vol. 22, no. 2, pp. 1266–1274, 2008.

J. van Gerpen, “Biodiesel processing and production,” Fuel Processing Technology, vol. 86, no. 10, pp. 1097–1107, 2005.

D. G. B. Boocock, S. K. Konar, V. Mao, C. Lee, and S. Buligan, “Fast formation of high-purity methyl esters from vegetable oils,” JAOCS, Journal of the American Oil Chemists' Society, vol. 75, no. 9, pp. 1167–1172, 1998.

D. Darnoko and M. Cheryan, “Kinetics of palm oil transesterification in a batch reactor,” JAOCS, Journal of the American Oil Chemists' Society, vol. 77, no. 12, pp. 1263–1267, 2000.

M. Diasakou, A. Louloudi, and N. Papayannakos, “Kinetics of the non-catalytic transesterification of soybean oil,” Fuel, vol. 77, no. 12, pp. 1297–1302, 1998.

D. Kusdiana and S. Saka, “Methyl esterification of free fatty acids of rapeseed oil astreated in supercritical methanol,” Journal of Chemical Engineering of Japan, vol. 34, no. 3, pp. 383–387, 2001.

Y. Liu, L. Wang, and Y. Yan, “Biodiesel synthesis combining pre-esterification with alkali catalyzed process from rapeseed oil deodorizer distillate,” Fuel Processing Technology, vol. 90, no. 7-8, pp. 857–862, 2009.

M. Slinn and K. Kendall, “Developing the reaction kinetics for a biodiesel reactor,” Bioresource Technology, vol. 100, no. 7, pp. 2324–2327, 2009.

E. Lelarasmee, A. E. Ruehli, and A. L. Sangiovanni-Vincentelli, “The waveform relaxation method for time-domain analysis of large scale integrated circuits,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 1, no. 3, pp. 131–145, 1983.

C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer, Berlin, Germany, 1988.

L. N. Trefethen, Spectral Methods in MATLAB, SIAM, Philadelphia, Pa, USA, 2000.

B. Freedman, E. H. Pryde, and W. F. Kwolek, “Thin layer chromatography/flame ionization analysis of transesterified vegetable oil,” Journal of the American Oil Chemists' Society, vol. 61, no. 7, pp. 1215–1220, 1984.

P. Sibanda, S. Motsa, and Z. Makukula, “Analysis of the synthesis of biodiesel using the homotopy analysis method,” Journal of Advanced Research in Applied Mathematics, vol. 2, no. 3, pp. 23–37, 2010.