Research Article: Process Optimization For Biodiesel Production From Corn Oil and Its Oxidative Stability
Research Article: Process Optimization For Biodiesel Production From Corn Oil and Its Oxidative Stability
Research Article: Process Optimization For Biodiesel Production From Corn Oil and Its Oxidative Stability
Research Article
Process Optimization for Biodiesel Production from Corn Oil and
Its Oxidative Stability
Copyright © 2010 N. El Boulifi 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.
Response surface methodology (RSM) based on central composite design (CCD) was used to optimize biodiesel production
process from corn oil. The process variables, temperature and catalyst concentration were found to have significant influence
on biodiesel yield. The optimum combination derived via RSM for high corn oil methyl ester yield (99.48%) was found to be
1.18% wt catalyst concentration at a reaction temperature of 55.6◦ C. To determine how long biodiesel can safely be stored, it is
desirable to have a measurement for the stability of the biodiesel against such oxidation. Storage time and oxygen availability have
been considered as possible factors influencing oxidative instability. Biodiesel from corn oil was stored for a period of 30 months,
and the physico-chemical parameters of samples were measured at regular interval of time. Results show that the acid value (AV),
peroxide value (PV), and viscosity (ν) increased while the iodine value (IV) decreased. These parameters changed very significantly
when the sample was stored under normal oxygen atmosphere. However, the ν, AV, and IV of the biodiesel sample which was stored
under argon atmosphere were within the limit by the European specifications (EN 14214).
as a by-product by using new technology that will make Factorial design of experiments gives more information
ethanol production more efficient. This corn oil can be per experiment than unplanned approaches; it allows to
converted into a biofuel, such as biodiesel [13]. see interactions among experimental variables within the
Biodiesel is being commercialized as a substitute or range studied, leading to better knowledge of the process and
blending stock of fossil-diesel. However, biodiesel is less therefore reducing research time and costs [22].
resistant to oxidation than typical fossil fuel and therefore
doping of biodiesel in fossil-diesel will affect the stability
2. Experimental Section
of fuel significantly [14]. Thus, as the demand and the
production of biodiesel grows fast, the development of 2.1. Equipment. Reactions were carried out in a batch stirred
methods to ensure the quality of the biodiesel industry reactor of 500 cm3 volume, equipped with a reflux condenser
and standardization becomes an urgent topic for market and a mechanical stirrer. The impeller speed was set at
introduction of biodiesel. 600 rpm to avoid external mass transfer limitation [6].
As corn biodiesel chemically is a mixture of long-chain
fatty acid methyl esters (FAMEs), it is more susceptible
2.2. Materials. Corn oil was supplied by Koipe Spain. The
to autoxidation and thus has a higher level of chemical
quality control of the corn oil used in this study is presented
reactivity than fossil-diesel. This oxidation instability is
in Table 1. The characteristics were determined according
dependant on the number and location of methylene-
to AOCS official method. Methanol of 99.8% purity was
interrupted double bounds in the FAMEs. Thus carbons
supplied by Panreac (Spain). The catalyst used was potassium
that are simultaneously allylic to two olefinic groups will
hydroxide purchased from Merck (Barcelona, Spain).
be extremely susceptible to the initiation of peroxidation.
An early study [15] measured the relative rate of oxidation
for the methyl esters of oleic (18 : 1), linoleic (18 : 2) and 2.3. Production of Corn Oil FAMEs. The corn oil FAMEs were
linolenic (18 : 3) acids to be 1 : 12 : 25. produced via a transesterification process in which a strong
A number of reports have appeared on the storage and alkali catalyst was used. This latter is frequently used in the
oxidative stability of biodiesel synthesized from different transesterification reaction, primarily due to its significant
vegetable and frying oils (including sunflower, soybean advantages in terms of the smaller quantity of catalyst used
and rapeseed oil) [16–18]. However, there is no available and the shorter reaction time required [6].
information about storage stability of biodiesel from corn oil. Corn oil was used as the raw oil to be transesterified
These stability studies were carried out under different stor- with methanol in a reacting tank. The temperature values
age condition, such as elevated temperature and exposure are below the boiling point of methanol (63◦ C), to prevent
to light, air, water, and other contaminants. In other studies the methanol in the reactant mixture from evaporating.
biodiesel was subjected to accelerated methods of oxidation, The molar ratio of the methanol and corn oil was set at
including EN 14112 (Rancimat method) and pressurized 6 : 1. The potassium hydroxide was stirred with methanol
differential scanning calorimetry [19, 20]. The oxidation for 10 minutes using an electric-magnetic stirrer to form
stability does not only depend on the raw material on the potassium methoxide, which was then poured into the
production process, but also on how the oil pressed and reacting tank and mixed with the corn oil. The total
refined [21]. reaction time was 60 minutes. Almost total conversion to
In a previous work [16], long-term storage tests on corn oil FAMEs was achieved quickly after a few minutes
biodiesel from different sources have been conducted under from the start of the reaction, depending on the reaction
different storage condition. Some chemicals properties have conditions.
been severely affected by oxidation of the FAMEs, when
the biodiesel was exposed to air and daylight. In order to 2.3.1. Fatty Acid Methyl Esters (FAMEs) Purification. At the
understand whether normal oxygen atmosphere availability end of a run the reaction mixture was allowed to cool down.
can adversely affect these properties, the study of the The upper phase consisted of FAME, and the lower phase was
change in chemical and physical proprieties of biodiesel glycerol.
in absence of oxygen (under argon atmosphere) becomes Once the glycerol and FAME phases have been separated,
necessary. the last one was purified by gentle washing with distilled
In this paper, the storage stability of biodiesel made from water to remove residual catalyst, glycerol, and soaps. The
corn oil was investigated over a storage time of 30 months pH of washing water was initially very high 10.22 due to
under argon atmosphere conditions, whose properties were the dissolved KOH. After 3 successive rinses with water, the
compared to a corn biodiesel sample which stored under washing water became clear and its pH was 7.9. The washing
normal oxygen atmosphere. The focus of the study was on process was continued (twice more) until a pH of about 7
the influence of storage time on the biodiesel properties, was achieved. Finally, the methyl ester phase was distilled to
such as PV, AV, IV and ν. This work also discusses the remove the residual water.
results of experiments carried out to evaluate the different The final water content of the corn oil FAME was less
variables affecting the alkaline methanolysis of corn oil. The than 0.01%. Water in the sample can promote microbial
optimum value for the variables affecting the process will be growth, lead to tank corrosion, participate in the formation
determined by application of factorial design and response of emulsions, as well as cause hydrolysis or hydrolytic
surface methodology. oxidation [23].
International Journal of Chemical Engineering 3
2.4. Analytical Method. Reaction products were monitored Table 2: Experimental matrix and experimental results.
by capillary column gas chromatography, using a Hewlett-
Run XC XT C (%) T (◦ C) Y (%)
Packard 5890 series II equipped with a flame ionization
detector (FID). The injection system was split-splitless. The 1 −1 −1 34.4 0.86 94.40
carrier gas was helium at a flow rate of 1 mL/minutes. The 2 +1 −1 55.6 0.86 98.25
analytical procedures and operating conditions have been 3 −1 +1 34.4 1.18 96.50
described in detail in a previous work [24]. The internal 4 +1 +1 55.6 1.18 99.48
standard technique was used to quantify the amount of the 5 0 0 45 1.02 98.25
chemical species. 6 0 0 45 1.02 97.85
7 0 0 45 1.02 98.30
2.5. Preparation the FAME Samples. Two biodiesel samples, 8 0 0 45 1.02 98.40
three litres each of were stored at ambient temperature 9 −α 0 30 1.02 97.50
for 30 months at two different storage conditions: the 10 +α 0 60 1.02 98.70
samples were stored in closed glass exposed to daylight,
11 0 −α 45 0.80 96.9
one under normal oxygen conditions and the other under
12 0 +α 45 1.24 96.4
atmosphere argon. During storage, samples were taken out
periodically and different quality parameters (PV, AV, IV, and Note: T: temperature; C: catalyst concentration; X: coded value; Y: biodiesel
yield.
ν) were monitored. The analyses were done according to the
following procedures: PV (AOCS Cd 8-53), AV (AOCS-Ca
5a-40), ν (ISO 3104), and IV (AOCS cd 1-25). response, yield of ester, may be thus expressed as a function
of the significant factors.
3. Results and Discussion
3.1.1. Linear Stage. A linear stage was considered in the first
3.1. Statistical Analysis. The synthesis of biodiesel was stud- step. Table 2 shows the experiments corresponding to the
ied using factorial design of experiments. The experimental 22 factorial design (experiments 1–4) and four experiments
design applied to this study was a full two-level factorial in the centre point to evaluate the experimental error
design 22 (two factors each, at two-levels) and extended to (experiments 5–8). A statistical analysis was carried out
response surface methodology (RSM). on these experimental values, and the main effects and
The response selected, Y, was the yield of methyl ester. interaction effects of the variables were calculated. The
The factors chosen were reaction temperature, XT , and initial analysis of the main effects and its interaction for the chosen
catalyst concentration, XC . Initial alcohol/oil molar ratio was response and the test of statistical significance are given in
fixed at 6 : 1. Table 3.
Selection of the levels was carried out based on the Temperature (XT ), catalyst concentration (XC ), and
results obtained in preliminary studies [25]. Temperature catalyst concentration-temperature interaction effects (XTC )
levels were selected according to reactants properties, thus were fitted by multiple regression analysis to a linear model.
the lower value was set at 34.4◦ C and the higher at 55.6◦ C. The response function for the significant main effects and
The levels of catalyst concentration were chosen on the basis interactions is
of preliminary experiments where the amount of catalyst was
progressively increased, and the ester yield monitored versus Y = 97.16 + 1.7XC + 0.83XT − 0.21XTC , r = 0.98.
time. The levels chosen were 0.86 and 1.18 wt.%, based on (1)
the total mass reaction.
The statistical analysis was thereafter applied. The exper- The statistical analysis of experimental results revealed
imental matrix for the factorial design is shown in Table 2. that the most significant factor is the catalyst concentration,
All the runs were performed at random. Four experiments while it also shows a significant value for curvature for the
were carried out at the centre point level, coded as “0”, for chosen responses. These data indicate the nonlinearity of the
experimental error estimation. model and thus justifies planning a more complex design to
The use of analysis and factorial design of experiments fit the data to a second-order model.
allowed us to express the amount of ester produced as a
polynomial model (if the levels of the factors are equally 3.1.2. Nonlinear Stage. To better predict the effect of vari-
spaced, then orthogonal polynomials may be used). The ables, a quadratic model was investigated. Here, the 22
4 International Journal of Chemical Engineering
Response: Y (%)
Number of experiments: 4
Degree of freedom: 3
Results of statistical analysis
Y (average) 97.16
Main effects and interactions
XT = 1.66, XC = 3.41, XTC = −0.43
Statistical Significance of t-test
Confidence level: 95%, Standard deviation, S = 0.24
Student’s t value, t = 3.182, Confidence interval: ±0.38
Significant effects and Interactions: XT (+), XC (+), XTC (−)
Statistical significance of curvature
YC = 98.2, Y = 97.16 Curvature: YC − Y − = 1.04
S = 0.11, t = 3.182 Confidence curvature interval: ±0.54
Curvature: significant
Response equation:
Y = 97.16 + 1.7XC + 0.83XT − 0.21XTC r = 0.98
2
2
2
94
Y = β0 + βi Xi + βii Xi2 + βi j Xi X j , (2)
i=1 i=1 i=
/ j
92
where Y is the response (methyl ester yield), Xi and X j 30
35
are the uncoded independent variables, and β0 , βi , βii , and 40
T 45 1.2
βi j are intercept, linear, quadratic and interaction constant (◦ 1.1
C) 50
55 1
coefficients, respectively. 0.9 )
The coefficients of (2) were determined by multiple 60
0.8 C (%
regression analysis. This analysis includes all the independent
Figure 1: Experimental yield versus temperature and catalyst
variables and their interactions, regardless of their signifi-
concentration.
cance levels. The best-fitting response surfaces found can be
written as follows:
Y = 98.2 + 1.06XC + 0.32XT − 0.21XT XC + 0.82XT2 significant in the range studied. This effect has a positive
(3)
− 0.1XC2 , r = 0.97. influence on the response. As the temperature increases, the
solubility of methanol in the oil increases and so does the
The statistical model was obtained from coded levels. speed of reaction. As a matter of fact, at low temperatures,
Equation (3) can be represented as dimensional surfaces and methanol is not soluble at all in the oil; when the stirring
contour plots, as shown in Figure 1. These show the ester is started an emulsion appears. The reaction takes place at
yield predicted for the experimental range of temperature the interface of the droplets of alcohol in the oil and then as
and initial catalyst concentration. The influence of these soon as the first FAMEs are formed, the alcohol solubilizes
variables on the ester yield will now be discussed. The progressively because the esters are mutual solvents for the
influence of the main factors and interactions will be derived alcohol and the oil.
from (1) and (3).
3.1.4. Influence of Catalyst Concentration. From the statistical
3.1.3. Influence of Temperature. For both linear and non- analysis it can be concluded that, within the experimental
linear models, the temperature influence is statistically range, initial catalyst concentration is the most important
International Journal of Chemical Engineering 5
80
Residual
50 0.2
20 −0.3
5
−0.8
1
0.1 −1.3
−2.2 −1.2 −0.2 0.8 1.8 2.8 3.8 94 95 96 97 98 99 100
Predicted
Figure 2: Normal probability plot for the methyl esters yield.
Figure 3: Residual plot of methyl esters yield for the second-order
model.
Table 4: Quality control of corn oil biodiesel fuels used in this study before and after storage time of 30 months compared to EN 14214
0.8 7
4
0.5 3
0.4 2
1
0.3
0
0.2 0 10 20 30
Month
0.1
Figure 8: Variation in the kinematic viscosity of corn oil biodiesel
stored under normal oxygen atmosphere and under argon
0
atmosphere.
0 10 20 30
Month
Figure 6: Variation in the acid value of corn oil biodiesel stored storage time (Figure 8). The sample which was stored under
under normal oxygen atmosphere and under argon atmosphere.
normal oxygen atmosphere appeared to have the highest rate
of reduction in IV, implying that the largest conversion of
7 polyunsaturated fatty acids to hydroperoxides took place. In
contrast to this result, the sample which was stored under
6 argon atmosphere was observed to have low rate of decrease
Viscosity at 40◦ C (mm2 /s)
to this study, and from an economical point of view, the [10] F. Ferella, G. Mazziotti, I. De Michelis, V. Stanisci, and F.
maximum yield of ester (98.75%) can be obtained, working Vegliò, “Optimization of the transesterification reaction in
with an initial catalyst concentration (1.1%) and operation biodiesel production,” Fuel, vol. 89, pp. 36–42, 2010.
temperature (55◦ C). These models are useful to determine [11] S. V. Ghadge and H. Raheman, “Process optimization for
the optimum operating conditions for the industrial process biodiesel production from mahua (Madhuca indica) oil using
using a minimal number of experiments with the consequent response surface methodology,” Bioresource Technology, vol.
economical benefit. 97, no. 3, pp. 379–384, 2006.
Biodiesel synthesised, which consists of long-chain [12] A. K. Tiwaria, A. Kumara, and H. Raheman, “Biodiesel
production from jatropha oil (Jatropha curcas) with high free
FAME, generally suffers from lower oxidation stability.
fatty acids: an optimized process,” Biomass and Bioenergy, vol.
Results from this study suggest that for a remarkably stable 31, no. 8, pp. 569–575, 2007.
biodiesel and in order to avoid oxidation, special precautions
[13] GS AgriFuels to Convert Corn Oil into Biodiesel at Ethanol
must be taken during long storage, such as storage under Facilities, GS AgriFuels Corporation, New York, NY, USA,
argon atmosphere. Nevertheless this action delays oxidation 2006.
but it doesnot prevent it. The results of this work have [14] R. O. Dunn and G. Knothe, “Oxidative stability of biodiesel in
allowed develop a methodology to overcome the obvious blends with jet fuel by analysis of oil stability index,” Journal of
problems of long storage stability of biodiesel. Long-term the American Oil Chemist’s Society, vol. 80, no. 10, pp. 1047–
storage study gives us a better understanding of the effect of 1048, 2003.
the normal oxygen atmosphere on the stability of biodiesel. [15] F. D. Gunstone and T. P. Hilditch, “The union of gaseous
The use of corn oil as an alternative raw material and oxygen with methyl oleate, linoleate, and linolenate,” Journal
renewable feedstock to produce biodiesel which fulfills the of the Chemical Society, pp. 836–841, 1945.
specification of EU standards for biodiesel (EN 14214) is [16] A. Bouaid, M. Martı́nez, and J. Aracil, “Long storage stability
of great interest to build an integrated and self-sustained of biodiesel from vegetable and used frying oils,” Fuel, vol. 86,
biorefinery. no. 16, pp. 2596–2602, 2007.
[17] A. Bouaid, M. Martı́nez, and J. Aracil, “Production of biodiesel
from bioethanol and Brassica carinata oil: oxidation stability
Acknowledgment study,” Bioresource Technology, vol. 100, no. 7, pp. 2234–2239,
2009.
Financial support from the (CICYT), Spanish Project CTQ [18] P. Bondioli, A. Gasparoli, L. D. Bella, and T. Silvia, “Evaluation
2006-10467/PPQ is gratefully acknowledged. of biodiesel storage stability using reference methods,” Euro-
pean Journal of Lipid Science and Technology, vol. 104, no. 12,
pp. 777–784, 2002.
References [19] B. R. Moser, “Comparative oxidative stability of fatty acid alkyl
[1] F. Staat and E. Vallet, “Vegetable oil methyl ester as a diesel esters by accelerated methods,” Journal of the American Oil
substitute,” Chemistry & Industry, vol. 21, pp. 856–863, 1994. Chemist’s Society, vol. 86, no. 7, pp. 699–706, 2009.
[2] M. P. Dorado, E. Ballesteros, J. M. Arnal, J. Gomez, and F. J. [20] R. O. Dunn, “Effect of oxidation under accelerated conditions
L. Gimenez, “Testing waste olive oil methyl ester as a fuel in a on fuel properties of methyl soyate (biodiesel),” Journal of the
diesel engine,” Energy & Fuels, vol. 17, no. 6, pp. 1560–1565, American Oil Chemist’s Society, vol. 79, no. 9, pp. 915–920,
2003. 2002.
[21] M. H. Chahine and R. F. Macneill, “Effect of stabilization of
[3] N. Usta, “An experimental study on performance and exhaust
crude whale oil with tertiary-butylhidroquinone and other
emissions of a diesel engine fuelled with tobacco seed oil
antioxidants upon keeping quality of resultant deodorized
methyl ester,” Energy Conversion and Management, vol. 46, no.
oil. A feasibility study,” Journal of the American Oil Chemist’s
15-16, pp. 2373–2386, 2005.
Society, vol. 51, no. 3, pp. 37–41, 1974.
[4] S. P. Singh and D. Singh, “Biodiesel production through the
[22] G. Box and J. Hunter, “Response surface methods,” in Statistics
use of different sources and characterization of oils and their
for Experiments Part IV: Building Models and Using Them,
esters as the substitute of diesel: a review,” Renewable and
chapter 5, John Wiley & Sons, New York, NY, USA, 1978.
Sustainable Energy Reviews, vol. 14, no. 1, pp. 200–216, 2010.
[23] G. Knothe, “Analyzing biodiesel: standards and other meth-
[5] J. M. Marchetti, V. U. Miguel, and A. F. Errazu, “Possible ods,” Journal of the American Oil Chemist’s Society, vol. 83, no.
methods for biodiesel production,” Renewable and Sustainable 10, pp. 823–833, 2006.
Energy Reviews, vol. 11, pp. 1300–1311, 2007.
[24] T. Garcia, A. Coteron, M. Martı́nez, and J. Aracil, “Optimiza-
[6] G. Vicente, M. Martı́nez, and J. Aracil, “Integrated biodiesel tion of the enzymatic synthesis of isopropyl palmitate using a
production: a comparison of different homogenous catalysts central composite design,” Transactions of Chemical Engineers,
systems,” Bioresource Technology, vol. 92, pp. 297–305, 2004. vol. 73, pp. 140–144, 1995.
[7] Y. C. Sharma and B. Singh, “Development of biodiesel from [25] G. Vicente, A. Coteron, M. Martı́nez, and J. Aracil, “Appli-
karanja, a tree found in rural India,” Fuel, vol. 87, no. 8-9, pp. cation of the factorial design of experiments and response
1740–1742, 2008. surface methodology to optimize biodiesel production,”
[8] G. Knothe, “Current perspectives on biodiesel,” Information, Industrial Crops and Products, vol. 8, no. 1, pp. 29–35, 1998.
vol. 13, no. 12, pp. 900–903, 2002. [26] A. Monyem, M. Canakci, and J. H. Van Gerpen, “Investigation
[9] H. Shi and Z. Bao, “Direct preparation of biodiesel from rape- of biodiesel thermal stability under simulated in-use condi-
seed oil leached by two-phase solvent extraction,” Bioresource tions,” Applied Engineering in Agriculture, vol. 16, pp. 373–378,
Technology, vol. 99, no. 18, pp. 9025–9028, 2008. 2000.
International Journal of Chemical Engineering 9
Rotating
Machinery
International Journal of
The Scientific
Engineering Distributed
Journal of
Journal of
Journal of
Control Science
and Engineering
Advances in
Civil Engineering
Hindawi Publishing Corporation Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
Journal of
Journal of Electrical and Computer
Robotics
Hindawi Publishing Corporation
Engineering
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
VLSI Design
Advances in
OptoElectronics
International Journal of
International Journal of
Modelling &
Simulation
Aerospace
Hindawi Publishing Corporation Volume 2014
Navigation and
Observation
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014
in Engineering
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014
Engineering
Hindawi Publishing Corporation
http://www.hindawi.com Volume 2010
Hindawi Publishing Corporation
http://www.hindawi.com
http://www.hindawi.com Volume 2014
International Journal of
International Journal of Antennas and Active and Passive Advances in
Chemical Engineering Propagation Electronic Components Shock and Vibration Acoustics and Vibration
Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation
http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014