Evaluation of Thermo-oxidative Stability of Biodiesel - Energy & Fuels

Jun 27, 2017 - In this work, we propose Fourier transform infrared spectroscopy (FTIR) in association with principal components analysis (PCA) and hie...
5 downloads 6 Views 2MB Size
Subscriber access provided by Olson Library | Northern Michigan University

Article

Evaluation of Thermo-oxidative Stability of Biodiesel Elton L Savi, Leandro S. Herculano, Gustavo V. B. Lukasievicz, Alex S. Torquato, Mauro Luciano Baesso, Nelson Guilherme Castelli Astrath, and Luis C Malacarne Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b00696 • Publication Date (Web): 27 Jun 2017 Downloaded from http://pubs.acs.org on July 2, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Evaluation of Thermo-oxidative Stability of Biodiesel Elton L. Savi,† Leandro S. Herculano,‡ Gustavo V. B. Lukasievicz,‡ Alex S. Torquato,‡ Mauro L. Baesso,† Nelson G. C. Astrath,† and Luis C. Malacarne∗,† †Departamento de Física, Universidade Estadual de Maringá, 87020-900, Maringá, PR, Brazil ‡Departamento de Física, Universidade Tecnológica Federal do Paraná, 85884-000, Medianeira, PR, Brazil E-mail: [email protected] Abstract The biodiesel susceptibility to oxidation is one of the major problems concerning its use as an alternative to diesel fuel. Although well characterized, the oxidation process cannot be completely prevented since it can be affected by a large number of factors, as fuel composition, storage conditions, contaminants, temperature and presence of air and light. In this work, we propose the Fourier Transform Infrared Spectroscopy (FTIR) in association with Principal Components Analysis (PCA) and hierarchical cluster analysis (HCA) as a method for monitoring the extent of oxidation degree of biodiesel in the low rate phase, before the end of the induction period, in which the changes in the physical properties as viscosity, mass density, and chain composition of the sample remain almost undetectable. A detailed investigation of thermo-oxidation of biodiesel is reported for a mixture of 50/50 percent of soybean oil and animal fat biodiesels. The biodiesel degradation was accelerated maintaining the temperature of the sample at 110o C under

1 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

constant air flux for different times. Oil stability index, mass density, viscosity, gross caloric value, esters composition, UV-VIS and FTIR spectra were measured in order to analyze the degree of degradation of each sample. The multivariate analysis in the FTIR spectra clearly shows the discrimination of the samples in the first stage of the degradation process. These results could be useful as a method for monitoring the fuel quality during the storage, helping among fuel producers, suppliers, and users.

Introduction Biodiesel derived from vegetable oil or animal fats is an alternative source of energy to substitute petroleum fuels. The main advantage of biodiesel is its sustainable eco-friendly nature reducing the environmental impact in comparison with conventional diesel fuel, presenting a possible net negative carbon dioxide balance. 1 Biodiesel is non-toxic and presents high viscosity, providing better lubrication to engines than mineral diesel. Although biodiesel produces about 10% less energy than petroleum diesel, its actual performance in the engine is practically the same in terms of power output and torque. Biodiesel consists of esters produced by transesterification of vegetable oils or animal fat triglycerides. Oxidation process is one of the major issues influencing biodiesel oxidative stability, and can be affected by a large number of factors including fuel composition, storage conditions, contaminants, temperature and presence of air and light. Various aspects of biodiesel stability have been reported describing the significance of oxidation stability, its oxidation chemistry, the methods used for characterizing stability, the factors known to influence stability, and consequences of biodiesel oxidation in diesel engines. 2–4 The oxidation could be in principle observed by measuring a number of physical properties such as mass density, viscosity, gross caloric value, peroxide index, esters composition, etc. In terms of oxidative evaluation, the Rancimat method is one of the most popular methods used to determine the oxidation stability by accelerating the aging process of the sample by exposing it to heat and increased volumes of air. The method measures the time that 2 ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

passes until oxidation takes place at a high rate - the induction time or oxidation stability index (OSI). However, the oxidation process characterization involves the study of full modification in the samples properties and composition. As this process is non-linear, the simple measurement of a specific property do not necessarily describe the degree of degradation. The early stages of oxidation typically follow first order reaction kinetics, in which quality parameters such as mass density, acid values, and viscosity remain almost unchanged. At the final phase, where the oxidation takes high rate, the changes in the physical parameters are significant and easily detected. Maintaining and monitoring the biodiesel fuel quality during storage presents a concern among fuel producers, suppliers, and users. In this direction, the availability of quick tests for the determination of the degradation degree represents an important step towards the improvement of biodiesel handling and trading. Raman, IR, and UV-Vis spectroscopies have already become very widespread techniques in the field of characterization and analysis of biofuels. For instance, 5,6 FTIR was used in conjunction with principal components analysis (PCA) and hierarchical cluster analysis (HCA) to classify biodiesels produced from different vegetable oils. The results showed that it is possible to develop a reliable and rapid methodology to identify used vegetable oils with raw materials in biodiesel production by applying the unsupervised techniques of multivariate analysis to the data. HCA and PCA were also used in the analysis of gas chromatography as a tool to classify/analyze biodiesels for chemical information. 7,8 Multivariate analysis was also used for predicting water content in biodiesel, 9 to monitor and quantify transesterification reaction, 10 in H-NMR spectra of fatty esters and biodiesels to infer the corresponding cetane number, 11,12 to associate and discriminate samples for potential use in forensic and environmental applications, 13 and for biodiesel/diesel blend classification. 14,15 Recently, the change of the dielectric property was found to be consistent with the Rancimat conductivity, infrared peak area, ultraviolet absorbance, and the residual mass during the degradation process. 16 In this work, FTIR spectroscopy in association with multivariate analysis is used to

3 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

differentiate the degradation degree of a commercial biodiesel in the low rate phase, before the end of the induction period. In this phase, the changes in the physical properties such as viscosity, mass density, and chain composition of the sample remain almost undetectable. The accelerated degradation is induced at the temperature of 110o C under constant air flux for different time periods. A set of physical properties is determined in order to analyze the degree of degradation of each sample. Measurements of OSI, mass density, viscosity, caloric power, esters composition, UV-VIS/FTIR spectra are performed and the results used to analyze the oxidation stability of biodiesels. The PCA shows that the first three components are sufficient for a complete discrimination of the degradation dynamics, especially in the low rate phase of the degradation process. So far as we know, similar work has not yet been reported.

Experimental Section Samples of commercial biodiesel were provided by BSBios Inc. with composition of 50/50 percent of soybean oil and waste animal fat from meat-processing industries. This biodiesel presents mostly fatty acid group with 16- and 18-carbon chain lengths (17.4% of palmitic acid, 31.3% of oleic acid, 30.6% of linoleic acid). The effects of the oxidation process were investigated maintaining ten samples under accelerated degradation conditions for different times. The thermo-degradation was carried out at 110o C under constant air flux of 10L/h for 1, 2, 3, 5, 7, 9, 12, 15 and 18 hours using the Rancimat method (892 Professional Rancimat, Metrohm). The high temperature accelerated the biodiesel degradation due to the enhancement in the rate of oxidation. These samples were used as stock for analysis. The oil stability index (OSI) or induction time was measured by the Rancimat method. The method is based on the plot of conductivity versus time and the inflection point is used to define the OSI time, according to the technical norm EN14112. The mass density and the viscosity were recorded by a densimeter (Aton-paar DMA 5000)

4 ACS Paragon Plus Environment

Page 4 of 18

Page 5 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

and a rheometer (Brookfield DV-III) coupled to a heat controller, respectively, according to the technical norms ASTM D4052 and D445. The gross calorific value or combustion power of the samples was obtained using a calorimeter (IKA Works C-2000) following the ASTM D240 norm. The tests were performed with approximately 0.5g of sample, which was analyzed in dynamic mode at 40o C and 15 bars of oxygen pressure with an accuracy of 0.1o C. Fatty acid methyl esters of the biodiesel were separated by gas chromatography (GC) (Perkin Elmer Clarus 680) equipped with a flame ionization detector (FID) and a fused-silica capillary column (100m, 0.25mm i.d. and 0.25µm of polysiloxane, Select FAME CP-7420, Agilent Technologies). The operation parameters were as follows: detector temperature, 275o C; injection port temperature, 240o C. Initially, the column temperature was held at 80o C for 1min. Subsequently, the temperature was increased to 160o C at a rate of 20o C/min, and then to 198o C at a rate of 1o C/min. After this period, it was once again raised to 250o C at a rate of 5o C/min for 5min, totalizing 58min of chromatographic run. The carrier gas flow (helium) was 1.1mL/min. The sample split ratio was 1 : 100. For identification, the retention times of the fatty acids were compared to those of standard methyl esters (Sigma, St. Louis, MO) and commercial oils of known composition, such as chia and soybean oils. Retention times and peak area percentages were automatically computed by a software (TotalChrom, Perkin Elmer). The main fatty acids compositions of all samples are shown in Table 1. The UV-Vis spectra were recorded by a spectrophotometer (T90+UV/Vis, PG Instruments Limited) with a precision of 0.3 nm from 300 nm to 700 nm with a step of 0.5 nm at 25o C. The samples were held in a 0.2 mm thick quartz cuvette. A FTIR spectrometer (PerkinElmer model Spectrum 400) was used to detect the infrared absorption of biodiesel using a Universal Attenuated Total Reflectance (UATR-FTIR) sensor, from 4, 000 cm−1 to 600 cm−1 , with a resolution of 4 cm−1 and 32 scans. To remove noise the spectra were then treated using first derivative procedure with a first-order polynomial. Mean-centered data were used as pre-processing tools for multivariate analysis. To obtain

5 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

a simplified and comprehensive two-dimensional view of the relationships that exist among various thermo-degraded samples, PCA 17,18 and HCA were carried out in this work. Both PCA and HCA were performed with Mathematica 7.0 software by using the Multivariate Statistics Package and Hierarchical Clustering Package. For PCA and HCA the Agglomerate method to find clustering hierarchically, with the Chessboard distance, was used.

Results and Discussion The oxidation induction time for the undegraded sample was measured following the standard technical norm EN14112 based on the analysis by the Rancimat instrument at 110o C. This value was calculated as the maximum change in the rate of oxidation defined by the maximum of the second derivative of the conductivity versus time curve, as showed in the inset of Fig. 1. At about 12h, the oxidation starts at a high rate - a period whose termination is defined by the detection of secondary decomposition products. The OSI time measured using the Rancimat method as a function of temperature is plotted in Fig. 1 (a) and in logarithm scale in (b). The results follow the model ln(OSI) = B0 + B1 T −1 , where T is the absolute temperature in K and B0 and B1 are constants. 19 The parameter B1 is related to the reaction activation energy, Ea (kJ/mol), by Ea = Rg B1 , in which Rg is the gas constant. The regression analysis shows a linear correlation (R2 = 0.99), yielding Ea = 77.8kJ/mol, which is in agreement with the literature value for fatty acids and ester obtained by isothermal differential scanning calorimetry. 20 The result provides a correlation for estimating OSI as a function of temperature, although an extrapolation of OSI estimation to reasonable storage temperatures is more problematic. 19 Figure 2 presents the values of mass density, kinematic viscosity, and gross calorific value for all the samples. All the properties present almost no variation for the samples under thermo-degradation time of up to 9 hours. After 9 hours, the changes in mass density and viscosity become apparent. At 12 hours, the measured values are outside the limit

6 ACS Paragon Plus Environment

Page 6 of 18

Page 7 of 18

Temperature (ºC) 70

80

90

100

110

120

130

(a) second order derivative

150

( s/cm)

OSI (h)

150

100

(Conductivity)

T=110ºC

75

0 0

4

8

12

Time (h)

50

0

ln (OSI)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

4

(b)

2 0 2.5

2.6

2.7 -3

1/T (10

2.8

2.9

-1

K )

Figure 1: Rancimat test for cell conductivity with time displaying oxidations stability time for the fresh biodiesel. (a) The OSI time measured using the Rancimat method for the undegraded sample as a function of temperature. (b) OSI in logarithm scale. Inset shows Rancimat test for cell conductivity in time at 110o C, following the EN14112 norm. values allowed for commercial biodiesel, according to the technical norm EN14214:2012. A decreasing of gross calorific value after 9 hours of thermo-degradation is also observed. As biodiesel is oxidized, there is a loss of hydrogen molecules or addition of oxygen molecules to the fuel that contributes to the fast decreasing in the combustion power or the amount of energy released per mass unit. The same slow change in the early stage of oxidation followed by a fast change in the second stage of oxidation is observed as for the density and viscosity results. The slow change in the physical properties during the initial period of oxidation is a wellknown result. 21 During the initial stages of the oxidation of biodiesel, the methylene allyl and bis allyl groups are more active and interact more with the double bonds to produce a carbon free radical. The resulting radicals interact with oxygen causing the formation of peroxides for the propagation step. Peroxides propagate through the continuous abstraction of hydrogen from a carbon to form another carbon radical and a hydroperoxide (ROOH). The newly formed carbon free radicals are again combined with oxygen and the propagation 7 ACS Paragon Plus Environment

Energy & Fuels

process continues. This is a continuous chain process that concludes in the termination step; the reaction terminates when two free radicals react with each other to yield stable products. In the induction period, the concentration of hydroperoxide remains very low, given the oxidation stability of the fatty oil or biodiesel under stress conditions. After this period the ROOH level increases rapidly and properties of fatty oils and biodiesel fuels are

0.94

density

0.92

(a)

0.90

0.88

10

(b) viscosity

8

calorific value

(kJ/g)

2

Kinematic (mm /s)

3

Mass (g/cm )

altered directly.

Gross

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 18

6

4

40

39

(c)

38

37 0

2

4

6

8

10

12

14

16

18

Degradation time (h)

Figure 2: (a) Mass density, (b) kinematic viscosity, and (c) gross calorific value or combustion power as a function of the thermo-degradation time. A great change in the values of density, viscosity and gross calorific values is observed only after the end of the induction period. Table 1 shows the GC-FID analysis. It reveals a fatty acid composition of 17.4% of palmitic acid (C16), 11.7% of stearic acid (C18), 30.6% of oleic acid (C18:1n9c), 31.3% of linoleic acid (C18:2n6c), 3.5% of α-linolenic (C18:3n3), residual quantities of other esters for the undegraded biodiesel, and the corresponding modification for the thermo-degraded samples. A similar behavior is observed; a small change in composition in the first stage of oxidation followed by a large decreasing of linoleic and α-linolenic acids and a corresponding increasing in palmitic, stearic and oleic composition in the second stage. The results are 8 ACS Paragon Plus Environment

Page 9 of 18

depicted in Fig. 3. Table 1: Main esters composition of all thermo-degraded samples obtained by CG-FID analysis. Analitic 0h 17.4 11.7 30.6 31.3 3.5

Palmitic (C16) Stearic (C18) Oleic (C18:1n9c) Linoleic (C18:2n6c) α-Linolenic (C18:3n3)

1h 17.5 11.7 30.6 31.3 3.5

2h 17.6 11.7 30.6 31.2 3.5

Normalized Area (%) 3h 5h 7h 9h 17.5 17.5 17.6 17.5 11.8 11.8 11.8 11.9 30.6 30.6 30.7 30.8 31.2 31.2 31.0 30.9 3.5 3.4 3.4 3.4

12 h 18.8 12.8 32.4 27.7 2.6

15 h 23.0 15.7 35.7 17.5 0.9

18 h 27.7 18.8 36.4 8.2 0.2

100

Concentration change (%)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

50

C16: Palmitic acid C18: Stearic acid C18:1n9c: Oleic acid

0

C18:2n6c: Linoleic acid C18:3n2:

-linolenic acid

-50

-100 0

2

4

6

8

10

12

14

16

18

Degradation Time (h)

Figure 3: Principal changes in the hydrocarbon chain composition for the thermo-degraded samples obtained by gas chromatography with flame ionization detection - CG-FID. After the end of the induction period a large decreasing of linoleic and linolenic acids and a corresponding increasing in palmitic, stearic and oleic composition is observed. Although slight changes are noticeable in the physical properties in the first stage of the oxidation process, no clear information regarding the biodiesel degree of oxidation is ascribed. In this direction, a quick test to monitor the quality of biodiesel fuel during storage, especially in the first stage of oxidation, represents an important factor for the improvement of biodiesel 9 ACS Paragon Plus Environment

Energy & Fuels

handling, trading and use. The UV-Vis absorption results are shown in Fig. 4. The thermo-degraded time influence on the evolution of the absorption bands is prominent. The strong absorption in the region of 400−500 nm is predominantly a result of carotenoids present in the biodiesel samples. 22,23 The decreasing of absorption in the first stage of degradation, as shown in the inset for 450 nm and 480 nm, indicates a bleaching by carotenoid loss in the samples. Despite the reduction of the optical absorption in the first stage of oxidation, the use of UV-Vis spectroscopy for identification of oxidative stability depends on the feedstock of the corresponding biodiesel. For this particular composition, the absorption measurements in the range of 400 − 500 nm seems to represent an easy, fast and non-expensive method to monitor the fuel quality during the storage period. However, the degradation of carotenoids in the sample could not necessarily give a direct correlation with the biodiesel oxidation, especially since the biodiesel degradation can be affected by a large number of factors, not only by thermal degradation. 3.0

= 450 nm

2.0

Absorbance

= 480 nm 2.5

2.0

1.5

1.0

0.5

Absorbance

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 18

0.0

1.5

0

3

6

9

12

15

18

Degradation Time (h) 1.0

0.5

0 h

7 h

1 h

9 h

2 h

12 h

3 h

15 h

5 h

18 h

0.0 400

450

500

550

600

650

700

Wavelength (nm)

Figure 4: UV-Vis spectra of thermo-degraded biodiesel for different times. The arrow represents the increasing degradation time (0, 1, 2, 3, 5, 7, 9, 12, 15 and 18 h). The inset shows the decreasing of absorption in the first stage of degradation for the 450 nm and 480 nm, indicating a bleaching by carotenoid loss.

10 ACS Paragon Plus Environment

Page 11 of 18

FTIR spectroscopy gives the characteristics of functional groups that constitute the components in the samples. For qualitative analysis, the mid-IR region (4000 − 400 cm−1 ) gives a fingerprint of functional groups present in the biodiesel. However, instead of a detailed analysis of each specific peak, a global analysis considering the whole spectrum can be obtained with the chemometric method. For instance, the mid-IR region spectroscopy is found to be efficient in the classification of oils and can be used as a rapid and simple discriminant analysis for oils and fats. 24 Here, we show that FTIR spectroscopy associated with chemometric analysis is also a powerful tool for monitoring or investigating biodiesel oxidation. It allows for the assessment of the relative oxidative state of a biodiesel sample, including the first stage, in which the main physical properties present slow variation. FTIR spectra are presented in Fig. 5. The curves are quite similar to one another, making it difficult to identify differences among the samples by simple inspection. This characteristic situation can commonly be addressed using principal components analysis. Additionally, hierarchical cluster analysis can also be used to monitor the biodiesel degree of oxidation. 100 Transmittance (%)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

0 h

90

1 h 2 h 3 h

80

5 h 7 h 9 h

70

12 h 15 h 18 h

60 500

1000

1500

2000

2500

3000

3500

4000

-1

Wavenumber (cm )

Figure 5: FTIR spectra of thermo-degraded samples for different times at 110o C and constant air flux.

Figure 6 shows the results from PCA analysis. The PCA is a statistical technique that uses an orthogonal transformation to convert a large amount of observations in a set of non 11 ACS Paragon Plus Environment

Energy & Fuels

correlated values, known as principal components. The PCA scores plot, PC1 vs. PC2 (PC1 contains 77% of the variance between the data sets and PC2 12% of variance), gives a clear characteristic discrimination from samples in the first stage process of degradation (negative values of PC1) from samples in the second stage of degradation (samples with large variation in physical properties). The two stages appear to be well separated from one another, suggesting that the degradation after the first stage has significant changes in chemical bonds compared to that observed in the second stage of degradation. In addition, the plot of PC1 versus PC3 (5% of variance) gives a increasing discrimination for the sample in the first stage of degradation, in which the change in the physical properties as viscosity, density, and chain composition of the sample remain almost undetectable. Note that the three first components contain 94% of the total variance, showing that with this orthogonal transformation of the PCA method almost all information of FTIR spectra are captured in the first components. A simple spectroscopic analysis taking a specific peak, as performed in Ref., 16 was not able to discriminate the degradation evolution in the low rate phase.

PC2 (12%)

2,0 18h

1,0

12h 3h

0,0

5h

9h

1h 2h 0h 7h

-1,0 Undegraded biodiesel

Degraded biodiesel

15h

-2,0

9h 0,5 7h

5h

0,0

3h 2h

-0,5

1h 0h

-1,0 -2

-1

Increasing degradation

1,0

PC3 (5%)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 18

12h 15h

18h

0

1

2

3

4

PC1 (77%)

Figure 6: PCA scores plot (P C2 × P C1 and P C3 × P C1) obtained from FTIR spectra of thermo-degraded biodiesel for different times at 110o C. Negative values of P C1 and increasing values of P C3 discriminate the samples in the low rate phase of the degradation process. Figure 7 shows the dendrogram plot of the statistical distance for the correlation matrix 12 ACS Paragon Plus Environment

Page 13 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

for all FTIR spectra. One can observe the increasing statistical distance (vertical amplitude in the plot) from the sample before the fast degradation process (0h to 9h) to the full degraded times (12h-18h). The results are in total agreement with the PCA scores plot, and is an additional tool for monitoring the dynamic evolution of thermo-degradation process. Note that this method also takes a global information from the FTIR spectrum without the need for analysis of specific peaks. The cluster analysis technique for data classification was used. Data elements are partitioned into groups called clusters that represent proximate collections of data elements based on a distance function. The partition data into exactly n lists of similar elements is found by clustering hierarchically. Increasing n leads to an increasing classification in the degradation. For instance, the gray area in Fig. 7 represents high light level n = 4. Additional sub-options methods are available to allow for more control over the clustering. For a given set of data and distance function, the choice of the best number of clusters may be unclear. A Significance Test can be used to determine statistically significant clusters to help choose an appropriate number. A possible value of Significance Test is Silhouette. The Silhouette test subdivides the data into successively more clusters looking for the first minimum of the silhouette statistic.

Figure 7: Dendrogram plot of a Hierarchical Cluster Analysis obtained from FTIR spectra of thermo-degraded biodiesel showing the statistical distance between samples. The increasing in the hierarchical distance is observed in function of thermo-degradation time.

13 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Conclusions To summarize, the chemometric techniques in association with FTIR spectroscopy proved to be a useful tool to monitor the evolution of the degradation process in biodiesel by taking an overall description of the evolution of the degradation. Measurements of OSI time, mass density, viscosity, gross calorific value, FTIR and UV-Vis spectra were employed to study thermo-degraded biodiesel. Upon comparing the samples at different thermo-degradation times, it becomes evident that it is not possible to determine the degradation evolution in the low rate phase only by monitoring the physical changes on the samples. The application of principal component analysis allows us to locate the thermo-degradation classification of the biodiesel samples, giving a clear discrimination on the stage of degradation. The spectral fingerprints generated by FTIR in addition to the chemometric techniques of PCA and HCA are sufficient information to differentiate the biodiesel stages of degradation and may help maintaining and monitoring fuel quality of biodiesel during the storage. The present paper suggests that it is possible to develop a methodology to identify the degree of degradation for biodiesel sample helping among fuel producers, suppliers, and users.

Acknowledgement The authors acknowledge the support from the Brazilian agencies CAPES, CNPq, and Fundação Araucária, and thank BSBios Inc. for providing the samples.

References (1) Demirbas, A. Political, Economic and Environmental Impacts of Biofuels: A Review. Appl. Energy 2009, 86, S108-S117. (2) Siddharth, J.; Sharma, M. P. Review of Different Test Methods for the Evaluation of Stability of Biodiesel. Renew. Sustain. Energy Rev. 2010, 14, 1937-1947. 14 ACS Paragon Plus Environment

Page 14 of 18

Page 15 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(3) Pullen, J.; Saeed, K. An Overview of Biodiesel Oxidation Stability. Renew. Sustain. Energy Rev. 2012, 16, 5924-5950. (4) Yaakob, Z.; Narayanan, B. N.; Padikkaparambil, S.; Unni K. S.; Akbar M. P. A Review on the Oxidation Stability of Biodiesel. Renew. Sustain. Energy Rev. 2014, 35, 136-153. (5) Mueller, D.; Ferrão, M. F.; Marder, L.; da Costa, A. B.; Schneider, R. C. S. Fourier Transform Infrared Spectroscopy (FTIR) and Multivariate Analysis for Identification of Different Vegetable Oils Used in Biodiesel Production. Sensors 2013, 13, 4258-4271. (6) Veras, G.; Gomes, A. A.; da Silva, A. C.; de Brito, A. L. B.; de Almeida, P. B. A; de Medeiros, E. P. Classification of Biodiesel Using NIR Spectrometry and Multivariate Techniques. Talanta 2010, 83, 565-568. (7) Schale, S. P.; Le, M. T.; Pierce, K. M. Predicting Feedstock and Percent Composition for Blends of Biodiesel with Conventional Diesel using Chemometrics and Gas Chromatography and Mass Spectrometry. Talanta 2012, 94, 320-327. (8) Pasupuleti, D.; Eiceman, G. A.; Pierce, K. M. Classification of Biodiesel and Fuel Blends using Gas Chromatography-Differential Mobility Spectrometry with Cluster Analysis and Isolation of C18:3 me by Dual Ion Filtering. Talanta 2016, 155, 278288. (9) Felizardo. P.; Baptista, P.; Menezes, J. C.; Correia M. J. N. Multivariate Near Infrared Spectroscopy Models for Predicting Methanol and Water Content in Biodiesel. Anal. Chim. Acta 2007, 595, 107-113, (10) Fontalvo-Gomez, M.; Colucci, J. A.; Velez, N.; Romanach, R. In-Line Near-Infrared (NIR) and Raman Spectroscopy Coupled with Principal Component Analysis (PCA) for In Situ Evaluation of the Transesterification Reaction. Appl. Spectrosc. 2013, 67, 1142-1149.

15 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(11) Ghesti, G. F.; de Macedo, J. L.; Braga, V. S.; de Souza, A. T. C. P.; Parente, V. C. I.; Figuerêdo, E. S.; Resck, I. S.; Dias, J. A.; Dias, S. C. L. Application of Raman Spectroscopy to Monitor and Quantify Ethyl Esters in Soybean Oil Transesterification. J. Am. Oil Chem. Soc. 2006, 83, 597-601. (12) Nadai, D. V. ; Simões, J. B.; Gatts, C. E. N.; Miranda, P. C. M .L. Inference of the Biodiesel Cetane Number by Multivariate Techniques.Fuel 2013, 105, 325-330. (13) Huppa, A. M.; Marshall, L. J.; Campbell, D. I.; Smith, R. W.; McGuffin, V. L. Chemometric Analysis of Diesel Fuel for Forensic and Environmental Applications. Anal. Chim. Acta 2008, 606, 159-171. (14) Silva, G. W. B. ; Gomes, A. A.; da Silva, P. ; Costa, G. B.; Fernandes, D. D. S.; Fontes, M. M.; Veras, G. Biodiesel/Diesel Blends Classification with Respect to Base Oil Using NIR Spectrometry and Chemometrics Tools. J. Am. Oil. Chem. Soc. 2012, 89, 1165-1171. (15) Pimentel, M. F.; Ribeiro, G. M .G. S.; da Cruz, R. S.; Stragevitch, L.; Pacheco Filho, J. G. A.; Teixeira, L. S. G. Determination of Biodiesel Content when Blended with Mineral Diesel Fuel using Infrared Spectroscopy and Multivariate Calibration. Microchem. J. 2006, 82, 201-206. (16) Zhou, J.; Xiong, Y.; Gong, Y.; Liu, X. Analysis of the oxidative degradation of biodiesel blends using FTIR, UV-Vis, TGA and TD-DES methods. Fuel 2017, 202, 23-28. (17) Jolliffe, I. T. Principal components analysis. Springer-Verlag: New York, 2002, 2nd ed. (18) Beebe, K. R.; Pell, R. J.; Seasholtz, M. B. Chemometrics: a practical guide. Willey: NewYork, 1996. (19) Dunn, R. O. Effect of Temperature on the Oil Stablility Index (OSI) of Biodiesel, Energy & Fuels 2008, 22, 657-662. 16 ACS Paragon Plus Environment

Page 16 of 18

Page 17 of 18

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(20) Dunn, R. O. Thermal Oxidation of Biodiesel by Pressurized Differential Scanning Calorimetry: Effects of Heating Ramp Rate, Energy & Fuels 2012, 26, 6015-6024. (21) Jain, S.; Sharma, M.P. Stability of biodiesel and its blends: A review. Renewable and Sustainable Energy Reviews 2010, 14, 667-678. (22) Lichtenthaler, H. K.; Buschmann, C. Chlorophylls and Carotenoids: Measurement and Characterization by UV-VIS Spectroscopy. Current Protocols in Food Analytical Chemistry 2001 F:F4:F4.3. (23) Stefanovich, A. F.; Karel, M. Kinetics of Beta-Carotene Degradation at Temperatures Typical of Air Drying of Foods. J. Food Process. and Pres. 1982, 6, 227-242. (24) Yang, H.;. Irudayaraj, J.; Paradkar, M. M. Discriminant analysis of edible oils and fats by FTIR, FT-NIR and FT-Raman spectroscopy. Food Chem. 2005, 93, 25-32.

17 ACS Paragon Plus Environment

Energy & Fuels

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical TOC Entry

18 ACS Paragon Plus Environment

Page 18 of 18