Optimization of Catalytic Glycerol Etherification with Ethanol in a

Mar 29, 2017 - ABSTRACT: In the past few decades, biodiesel has emerged as a promising renewable energy source. However, its production ... developmen...
0 downloads 15 Views 682KB Size
Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI

Article

Optimization of catalytic glycerol etherification with ethanol in a continuous reactor Caroline Ortega Terra Lemos, Leticia L. Rade, Marcos A.S. Barrozo, Lindoval Domiciano Fernandes, Lucio Cardozo-Filho, and Carla E. Hori Energy Fuels, Just Accepted Manuscript • DOI: 10.1021/acs.energyfuels.7b00194 • Publication Date (Web): 29 Mar 2017 Downloaded from http://pubs.acs.org on April 3, 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 24

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

Optimization of catalytic glycerol etherification with ethanol in a continuous reactor Caroline O. T. Lemosa, Leticia L. Radea, Marcos A. de S. Barrozoa, Lindoval D. Fernandesb, Lucio Cardozo-Filhoc,d, Carla E. Horia,*

a

Faculdade de Engenharia Química, Universidade Federal de Uberlândia, Av. João Naves de Ávila 2121, Campus Santa Monica - Bloco 1K, 38400-902, Uberlândia, MG, Brazil.

b

Departamento de Engenharia Química, Universidade Federal Rural do Rio de Janeiro. Rodovia BR 465, Km 07, s/n - Zona Rural, 23890-000, Seropédica, RJ, Brazil.

b

Departamento de Agronomia e Departamento de Engenharia Química, Universidade Estadual

de Maringa, Av. Colombo 5790 – Bloco D90 Jardim Universitário, 87020-900, Maringa, PR, Brazil. d

Centro Universitário da Fundação de Ensino Octávio Bastos, Av. Dr. Octávio Bastos, 2439 - Jd Nova São João, São João da Boa Vista, SP, Brazil.

* Corresponding author. Email address: [email protected]

KEYWORDS: Glycerol, catalytic etherification, continuous reactor.

ACS Paragon Plus Environment

1

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

Page 2 of 24

ABSTRACT: In the last decades, biodiesel has emerged as a promising renewable energy source. However, its production usually leads to the formation of large amounts of glycerol as byproduct. This paper evaluates the catalytic glycerol etherification with ethanol over a series of catalyst in continuous reactor. Amberlyst 15 presented the best performance in terms of glycerol conversion and yield of ethers. This result was attributed to its high pore size and acidity. Beta zeolite (Si/Al 12.5) and niobic acid, both with lower acid strength and pore size, were almost inactive to convert glycerol into ethers. A central composite design was developed to optimize the performance of Amberlyst 15. An increase in catalyst amount enhanced both glycerol conversion and yield of ethers. However, high temperatures and low molar ratio favored side reactions. The optimized value for glycerol conversion was 91% and for yield of ethers was 13%, at different reaction conditions.

ACS Paragon Plus Environment

2

Page 3 of 24

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

1. INTRODUCTION

In order to face the energy needs of this century, society demands the development of clean and renewable sources to replace fossil fuels.1,2 Biodiesel emerged as a promising renewable energy source 3,4. However, its main production route (transesterification of vegetable oils) leads to the formation of large amounts of glycerol as byproduct.5 Despite already being employed in food and beverage, pharmaceutical, cosmetics and paper industries, the production of glycerol is still much higher than its demand. Additionally, glycerol cannot be used directly as an additive in fuels due to its high polarity, hygroscopicity and decomposition in acrolein.6 Therefore, there is an urgent need to find new applications for the surplus glycerol. The catalytic glycerol processing may occur by hydrogenolysis, polymerization, etherification, oxidation, dehydration, acetylation and transesterification. The etherification reaction is the most promising because it generates oxygenated compounds which may be directly added to fuels.7,8 The most common glycerol etherification method is the reaction with alkenes catalyzed by acids. However, according to Ferreira et al.9 the glycerol etherification with primary and secondary alcohols is the most important route in terms of fuel industrial application.5,7,10,11 Klepacova et al.12 studied the etherification of glycerol with tert-butyl alcohol catalyzed by resins of Amberlyst type and by H-Y and H-Beta zeolites in a batch reactor. Amberlyst is very active reaching 100% of glycerol conversion with selectivity to di- and tri-ethers above 92%. Cannilla et al.13 also investigated the etherification of glycerol with butanol in a batch reactor with a tubular water perm selective membrane using Amberlyst 15 as solid acid catalyst. According to the authors, relevant glycerol conversion was reached at temperatures higher than 140 °C. However, at such temperatures, the removal of water through the membrane favored

ACS Paragon Plus Environment

3

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

Page 4 of 24

glycerol dehydration. Ozbay et al. 14, 15 investigated the etherification of glycerol with tert-Butyl alcohol in batch and continuous flow reactors. Both studies showed the advantages of operating in a continuous mode, as this type of operation provides better contact between the catalyst and the reaction mixture and constant water removal. Pariente et al.7 analyzed the performance of different types of heterogeneous acid catalysts, including sulfonic resins, zeolites and graphited silicas in the glycerol etherification with ethanol. The best results were achieved using sulfonic-acid polystyrene resins of the Amberlyst family and zeolites with Si/Al ratios around 25. Yuan et al.5 studied the glycerol etherification with ethanol over a series of catalyst such as H-ZSM5, H-β and tungstophosphoric acid (HPW). HPW showed the highest conversion of glycerol (97%) at 160 °C, ethanol/glycerol of 6:1. Until now, the glycerol etherification reaction with ethanol was only carried out in batch type reactors. This feature decreases the conversion of glycerol into ethers due to the high water content in the reaction mixture at the end of the reaction. The presence of water may facilitate the reversible reaction of the ethers hydrolysis.5 Moreover, in batch reactors the energy costs are high and catalyst must be removed from the product, among other problems. Considering the relevance of the glycerol etherification route, the aim of this work is to study glycerol etherification with ethanol in a continuous reactor over three catalysts: commercial strong acid ion-exchange resins (Amberlyst 15), beta zeolite and niobic acid (HY-340). A lightoff curve was used to select the catalyst that showed the best performance of glycerol conversion/yield to ethers in a temperature range of 150-300°C. For this catalyst, a design of experiments (DOE) was performed to evaluate the effect of temperature, catalyst amount and ethanol:glycerol molar ratio as well as their interactions in the glycerol conversion and yield to ethers. Finally, a canonical analysis was applied in order to optimize the response variables.

ACS Paragon Plus Environment

4

Page 5 of 24

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

2. EXPERIMENTAL SECTION

2.1.

Materials

Glycerol and ethanol (purity 99.5%) supplied by Synth were used as reactants. Analytical standards for GC analysis were supplied by Sigma-Aldrich. The macroreticular sulfonic resins Amberlyst A15 provided by Dow Chemical Company was received in dry form and dried at 110 °C for 1 hour. The niobic acid (HY-340) was supplied by Companhia Brasileira de Metalurgia e Mineração (CBMM) in powder form and calcined in air flow at 300 °C for 2 hours to remove adsorbed impurities and water. The beta zeolite sample Si/Al of 12.5 was prepared according to methodology described in the U.S. patent.16 In this synthesis, a gel was prepared with the following molar composition: 25 SiO2: 1.0 Al2O3: 1.5 Na2O: 4.5 TEAOH: 240 H2O. To this end, 18.32 g of sodium aluminate were added in a mixture containing 148.78 g of 40% TEAOH and 95.66 g of deionized water. After stirring for 10 minutes, 337.23 g of colloidal silica (Ludox HS40, Sigma-Aldrich) were added and resulting suspension was stirred for 30 minutes and then heated at 150 °C for 10 days. The solid formed was recovered by filtration, dried at 100 °C for one night and calcined at 540 °C for 8 hours. For calcination it was used a heating rate of 0.2 °C/min, with two levels of 30 minutes each at 150 °C and 350 °C. Finally, the ion exchange was carried out using NH4Cl followed by calcination at 550 °C for 8 hours.

2.2.

Catalyst Characterization

Powder X-ray diffraction (XRD) of the samples was recorded on a Rigaku Miniflex II powder X-ray diffratometer using nickel filtered CuKα radiation (λ = 1.540 Å). The measurements for

ACS Paragon Plus Environment

5

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

Page 6 of 24

beta zeolite were obtained in the 2θ range of 2 to 50° (steps of 0.05°) with an account time of 3 seconds. XRD data for niobic acid were recorded in a 2θ range from 20 to 80° at a scan rate of 0.05°/step and scan time of 2 s/step. The total acidity of the materials was measured by temperature programmed desorption of NH3 (TPD-NH3), in a Quantachrome Chembet-3000 equipment. Firstly, the samples were treated at 300 °C under N2 flow rate of 20 mL/min by 1 h. The adsorption of NH3 was carried out at 100 °C for 30 minutes. Physisorbed NH3 was purged during 2 h under N2 flow rate of 20 mL/min. The heating rate was then adjusted to 10 °C/min, from 100 °C to 700 °C, to perform the desorption of chemisorbed NH3. Textural analysis of the catalysts was carried out using N2 gas adsorption–desorption isotherms at -196.15 °C using Quantachrome Nova 1200 equipment. Surface area measurements were made according to the BET method and the pore size distribution was calculated using the tmethod (micropores) and BJH method (mesopores) with data determined from the desorption branch of the isotherms. The validity of the method was based on the BET constant C associated to the interaction of adsorbate/adsorbent, and c>0 was employed as isotherms adjustment region.

2.3.

Reaction Procedure

Etherification reactions were carried out in a fixed bed using a tubular reactor made of stainless steel tubing (316 L 3/8 in OD, inner diameter 6.8 mm, 30 cm). The feeding of substrate in the system was controlled by a high pressure liquid pump (LabAlliance, Series III) which was adjusted to a flow rate of 0.40 mL/min. The residence time was calculated as the time that the substrate remains in contact with the catalyst, which ranged from 12 to 80 seconds for 0.12 and 0.80 g of catalyst respectively. After filling the system with the substrate, the process of heating

ACS Paragon Plus Environment

6

Page 7 of 24

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

the furnaces started. The temperature was monitored by thermocouples connected to the reactors, with a precision of 2 °C. Finally, when the temperatures were stabilized at the desired conditions and after waiting twice the time necessary to fill all the reactor system, samples were collected. The etherification of glycerol with ethanol yields the following desired products: monoethylglycerols (MEG), di-ethylglycerols (DEG) and tri-ethylglycerols (TEG). Reaction products and remaining glycerol were diluted with ethanol and analyzed by gas chromatography (Shimadzu-2014) using a CP-WAX 52 CB column (30 m x 0.25 mm x 0.25 µm) and a flame ionization detector (FID) according to the methodology of Melero et al.11 The retention times were 6.9, 5.9, 6.1, 9.6 minutes for MEG, DEG, TEG and glycerol, respectively. The injection and detector temperatures were 250 °C. Analyses were carried out in split mode of 1:100 with temperature program from 50 to 210 °C (with a heating step of 70 °C/min) and at 210 °C for 10 minutes isothermally. The 1,4-butanediol analytical standard was used after a calibration with commercial MEG and glycerol. In the case of DEG and TEG, non-commercially available, response factors were extrapolated from that of MEG assuming similar behavior. Glycerol conversion (X), selectivity of MEG/DEG/TEG (S) and yield of MEG/DEG/TEG (Y) were calculated using the following equations:

=

       ⁄      ⁄

× 100

(1)

=

      ⁄   ⁄

× 100

(2)

=

         ⁄ ×      ⁄

100

(3)

ACS Paragon Plus Environment

7

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

Page 8 of 24

All the solid materials were tested in a temperature range of 150 to 300 °C, molar ratio ethanol/glycerol of 10/1 and 0.4 g of catalyst in order to evaluate the best catalyst to convert glycerol into ethers. The catalyst which presented the best performance was then chosen to have its performance optimized using a design of experiments methodology.

2.4.

Statistical analysis using design of experiments

The design of experiments (DOE) was used for quantifying the effects of the most significant variables, as well as their interactions on the glycerol conversion and yield of ethers. A Central Composite Design (CCD) with response surface method was used to determine the optimal reactions conditions. Three independent variables were selected:

temperature (XT),

ethanol/glycerol molar ratio (XM.R.) and catalyst amount (XC.A.). Each variable was investigated at five levels: -α, -1, 0, +1 and +α (coded levels), as shown in Table 1, where α is the axial point of 1.414 (orthogonal design). The selection of the levels was based on results obtained in preliminary tests. According to the experimental results for glycerol conversion and yield of ethers, a regression or statistical model was fitted. The significant terms in the model were found by analysis of variance (ANOVA) for each response. The statistically non-significant terms (p > 0.10) were removed from the initial models by stepwise selection and the experimental data were refitted to produce the final model. The central point experiment was repeated four times in order to determine the variability of the results and assess the experimental error. In all experiments performed, the feeding of substrate in the system was adjusted to a flow rate of 0.40 mL/min.

2.5.

Optimization using canonical analysis

ACS Paragon Plus Environment

8

Page 9 of 24

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

Canonical correlation analysis was used to study the interrelationships among sets of multiple dependent and independent variables for glycerol conversion and yield of ethers. An empirical equation was fitted to the experimental data which correlates the predict response with independent variables and first order parameters vectors and a matrix of the quadratic parameters. The maximum value for the response will be a set of conditions such that the derivatives are simultaneously zero. This point is called as stationary point of the fitted surface and can be a maximum, minimum or a saddle point, depending of the arrows sign (!i). Thus, the response function can be expressed in terms of new variables whose axes correspond to the principal axes of the contour system. The form of the function in terms of these variables is called canonical form and is given by Equation 4. The software Maple 17 was used to perform this analysis.

" = "# + !% &% ' + !' &' ' + … + !) &) '

(4)

3. RESULTS AND DISCUSSION

3.1.

Characterization of catalysts

XRD patterns of beta zeolite exhibited only the typical features, showing the characteristic diffraction peaks of beta zeolite, at 2θ = 7.5 – 8° and 22.4°, thus indicating high crystallinity without impurities.17 XRD patterns of niobic acid calcined at 300 °C showed an amorphous structure and this result is in agreement with the data reported by Lebarbier et al.18 Table 2 summarizes the most relevant physicochemical properties for the Amberlyst 15, beta zeolite and niobic acid. The acidity of the three samples was evaluated by NH3 adsorption and

ACS Paragon Plus Environment

9

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

Page 10 of 24

the results are presented in Table 2. The acidity of Amberlyst 15 was significantly higher than the other samples, followed by beta zeolite and niobic acid. Catalysts showed different surface areas, ranging from 53 to 434 m2/g and porosities between 0.0042 and 0.40 cm3/g. Amberlyst 15 presented a mesopore volume and mesopore size much higher than the other catalysts studied. Beta zeolite and niobic acid showed to be structured both by micro and mesopores. The presence of mesopores should favor the access to the active sites.

3.2.

Light off results

The light off curves for each type of catalyst are presented in Figure 1. The catalysts studied were active in different temperature ranges. The highest glycerol conversion value (85 %) was reached using Amberlyst 15 as solid acid at 200 °C. However, it should be noted that even with a low conversion of 53 % at 175 °C the ether yield was higher than the yield obtained at higher temperature. This behavior was repeated for the other catalysts studied, since the increase in the reaction temperature favors the formation of undesirable compounds such as acrolein and diethyl ether. Moreover, at high temperatures, water forms at a higher rate, favoring the inverse reaction and competing for the acid sites. In the presence of water, thermodynamic constraints and problems related to the stability of the catalyst may inhibit the reaction12,

19, 20

. Beta zeolite

reached a maximum conversion of 13 % at 300 °C. For this catalyst, glycerol conversion increases linearly with increasing temperature, however the ether yield remained very low over the range evaluated. One possible reason for this small conversion is the lower acidity presented by this sample compared to Amberlyst 15. Another explanation for this low efficiency observed for Beta zeolite is its small pore size (micropores around 6 Å). The average molecule size of glycerol é 4.3 Å and its smallest projection is around 4 Å. The addition of ethoxy groups to

ACS Paragon Plus Environment

10

Page 11 of 24

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

glycerol molecule will increase the molecule size (see molecules below). For MEG, the smallest projection would be around 5 Å and the diffusion would still be viable. For DEG and TEG, the diffusion limitations would be very severe, as previously reported by Yuan et al.5 and Klepacova et al.12. On the other hand, Amberlyst 15 is characterized by an average pore diameter much larger than beta zeolite. The niobic acid was quite inactive in the glycerol conversion and this fact may be related to its low acidity when compared to the other tested catalysts. Suitable textural properties of an open porous accessible structure and the high acidity characteristic of Amberlyst 15 lead to a best catalytic activity and production of desired products. Thereby, this catalyst was chosen to be used in the next steps of the study.

3.3.

Development of regression model and response surface

The experimental results obtained at different operating conditions used in the glycerol etherification with ethanol are presented in Table 3. As it can be seen, the glycerol conversion at different reaction conditions showed a large sensitivity to the variables studied (ranging from 0 to 91%). The highest yield of ethers (9%) was obtained by entry 7 at 167 °C, 14:1 ethanol:glycerol molar ratio and 0.68 g of catalyst. The temperature influences the distribution of the desired products. A higher temperature accelerates glycerol conversion and enhances the selectivity toward products DEG and TEG, whereas the selectivity to MEG decreased, as it can be seeing comparing entries 1 and 2, 3 and 4, 5 and 6 and 7 and 8. The influence of ethanol/glycerol molar ratio on the products selectivity cannot be easily evaluated. However the yield of ethers enhances with the increase of the molar ratio, according to entries 1 and 3, 5 and 7 and, 6 and 8. The glycerol conversion and yield of ethers are good

ACS Paragon Plus Environment

11

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

Page 12 of 24

results considering that the residence time (from 12 until 80 seconds) is much lower than the literature data presented until now.5, 7, 11 The arithmetical averages and the standard deviations were calculated for the central point experiments: glycerol conversion of 61% ± 1% and yield of ethers of 4% ± 1%. The results obtained were then analyzed using analysis of variance (ANOVA) and is the results are shown in Table 4. The ANOVA summary showed that the models obtained for conversion and yield of ethers had both a good correlation with the response variable. For glycerol conversion the quadratic term of ethanol/glycerol molar ratio and the three interaction between variables terms were not significant at a confidence interval of 90%. On the other hand, for yield of ethers only the interaction temperature/molar ratio and the quadratic term of catalyst amount were not significant (p-value higher than 0.1). From the matrix generated by the experimental data (Table 3) and after eliminating the insignificant parameters for CCD, a multiple regression was developed for each response variable. The empirical equations obtained, expressed in coded factors, is given by Equations 5 and 6.

 = 59 + 15.98. + 12.7112 + 26.2145 − 7.50. ' − 745' (r2=0.92)

(5)

 = 3.67 − 0.42. + 1.1012 + 2.4145 − 0.38.45 + 0.881245 − 0.71. ' + 0.5412 ' (r2=0.98)

(6)

Figures 2 and 3 show the response surfaces for glycerol conversion and ethers yield predicted by the statistical models presented on Equations 5 and 6. Statistical analysis of the studied

ACS Paragon Plus Environment

12

Page 13 of 24

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

experimental range identifies the catalyst amount (CA) as the most important factor in the glycerol conversion response (Eq. 5). A higher amount of catalyst leads to a consequently longer residence time and, in this way, the glycerol conversion will be beneficiated. The second most important factor is the temperature (T) followed by the ethanol/glycerol molar ratio (MR). That means that an increase in any of these variables produces an increase in the conversion of glycerol. However the enhancement of this response with temperature is more significant at high values of the molar ratio and the corresponding improvement with the molar ratio is also more substantial at high temperatures. Nevertheless, under high temperatures, a lot of non-desired compounds are produced by side reactions. Equation 5 shows two quadratic terms which indicates a non-linear behavior of the system: at higher catalyst amount and temperature, these parameters will start to have a negative effect on the glycerol conversion due to its degradation process and problems related to mass transfer. In terms of ethers yield, it is observed from Equation 6 that this response is affected by the linear terms of all studied parameters. Molar ratio and catalyst amount have positive effects on ethers yield, but temperature has a negative effect. The glycerol etherification is enhanced when the amount of ethanol present in the reaction medium is increased. In addition, the quadratic effects of temperature and molar ratio have a significant influence on the ethers yield, which indicates that the increase or decrease in these variables do not produce a constant variation on response due to curvature effects.

3.4.

Optimization

Aiming to have a larger set of data close to the region of the highest values of the responses (glycerol conversion and yield of ethers), additional experiments were performed in this region.

ACS Paragon Plus Environment

13

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

Page 14 of 24

Table 5 shows the results of these new experiments. The experimental values and coded reaction parameters of Table 5 were coupled to the values of central composite design and a new multiple regression was fitted for each response variable. Thus, the canonical analysis for the new fitted models was applied and the stationary points have been found. Canonical forms of these new models are presented in Equations 7 and 8. It can be seen that the stationary points are maximum points for both response variables (negative values for the eigenvalues).

"   9 = −28.38&% − 22.42&' − 11&: + 93.44

(7)

"9 = −6.78&% − 1.42 − 0.88&: + 12.77

(8)

The maximum predicted glycerol conversion was 93% at 193 °C, glycerol/ethanol molar ratio of 10.98/1 and 0.61 g of catalyst. The experimental value of glycerol conversion at this reaction condition was 91%. For yield of ethers the maximum predicted value was 13% at 170 °C, molar ratio of 17.34/1 and 0.57 g of Amberlyst 15. The experimental result obtained at this operating conditions was also 13%. Values of ethers yield are not usually presented in previous studies. However, Melero et al.11 presented the ethers yield as optimized response variable and the maximum value obtained was 40% using ethanol and modified-SBA-15 with arenesulfonic as catalyst, after 4 hours of batch reaction. Therefore, maximum values of the response variables optimized in this work are good results since the residence time is much smaller than the residence time presented on the current literature data. Furthermore, the use of continuous mode reactors proposed in this study favors the intensification and scalability of the process.

ACS Paragon Plus Environment

14

Page 15 of 24

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. CONCLUSIONS

Amberlyst 15 showed better performance in terms of glycerol conversion and yield of ethers than beta zeolite and niobic acid. According to the central composite design, the conversion of glycerol was influenced positively by catalyst amount followed by temperature and molar ratio. For yield of ethers, the catalyst amount was the most important positive factor followed by molar ratio, however, temperature had a negative effect. The canonical analysis led to a maximum value of glycerol conversion of 91% and a maximum yield of ethers of 13%, even with a residence time around 60 seconds in both cases.

AUTHOR INFORMATION *Carla Eponina Hori Faculdade de Engenharia Química, Universidade Federal de Uberlândia, Av. João Naves de Ávila 2121, Campus Santa Monica - Bloco 1K, 38400-902, Uberlândia, MG, Brazil. Email address: [email protected]

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. All the authors contributed equally.

ACS Paragon Plus Environment

15

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

Page 16 of 24

ACKNOWLEDGMENT

The authors wish to acknowledge the financial support of CAPES, CNPq, Vale S.A. and FAPEMIG and CBMM for providing the niobic acid sample (HY-340).

ABBREVIATIONS

HPW, tungstophosphoric acid; DOE, design of experiments; GC, gas chromatography; CBMM, Companhia Brasileira de Metalurgia e Mineração; XRD, powder X-ray diffraction; TPD-NH3, temperature programmed desorption of NH3; MEG, mono-ethylglycerols; DEG, diethylglycerols; TEG, tri-ethylglycerols; FID, flame ionization detector; CCD, central composite design; ANOVA, analysis of variance; CA, catalyst amount; T, temperature; MR, molar ratio.

REFERENCES

(1) Galan, M. I., Bonet, J., Sire, R., Reneaume, J. M., Ples, A. E. Bioresour. Technol. 2009, 100, 3775–3778. (2) Sivaiah, M. V., Robles-Manuel, S., Valange, S., Barrault, J. Catal. Today 2012, 198, 305– 313. (3) Izquierdo, J. F., Montiel, M., Pale, I., Outo, P. R., Gala, M., Jutglar, L., Villarrubia, M., Izquierdo, M., Hermo, M. P., Ariza, X. Renew. Sust. Energ. Rev. 2012, 16, 6717–6724. (4) Di Serio, M., Casale, L., Tesser, R., Santacesaria, E. Energy Fuels 2010, 24, 4668-4672.

ACS Paragon Plus Environment

16

Page 17 of 24

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

(5) Yuan, Z., Xia, S., Chen, P., Hou, Z., Zheng, X. Energy Fuels 2011, 25, 3186-3191. (6) Srinivas, M., Raveendra, G., Parameswaram, G., Sai Prasad, P. S., Lingaiah, N. J. Mol. Catal. A. Chem. 2016, 413, 7–14. (7) Pariente, S., Tanchoux, N., Fajula, F. Green Chem. 2009, 11, 1256-1261. (8) Gholami, Z., Abdullah, A. Z., Lee, K. T. Appl. Catal., A: Gen. 2014, 479, 76–86. (9) Ferreira, M. O., Cardozo-Filho, L., Silva, C., Sousa,E. M. B. D. Lat. Am. Appl. Res. 2014, 44, 47-56. (10) Frusteri, F., Arena, F., Bonura, G., Cannila, C., Spadaro, L., Di Blasi, O. Appl. Catal., A: Gen. 2009, 367, 77–83. (11) Melero, J. A., Vicente, G., Paniagua, M., Morales, G., Munoz, P. Bioresour. Technol. 2012, 103, 142-151. (12) Klepacova, K., Mravec, D., Bajus, M. Appl. Catal., A: Gen. 2005, 294, 141–147. (13) Cannilla, C., Bonura, G., Frusteri, L., Frusteri, F. Chem. Eng. J. 2015, 282, 187–193. (14) Ozbay, N., Oktar, N., Dogu, G., Dogu, T. Ind. Eng. Chem. Res. 2012, 51, 8788–8795. (15) Ozbay, N., Oktar, N., Dogu, G., Dogu, T. Top. Catal. 2013, 56, 1790-1803. (16) Wadlinger, R.L., Kerr, G.T., Rosinski, E.J. U.S. patent 3,308,069, March 1967. (17) Gonzalez, M. D., Salagre, P., Linares, M., Garcia, R., Serrano, D., Cesteros, Y. Appl. Catal., A: Gen. 2014, 473, 75–82. (18) Lebarbier, V., Houalla, M., Omfroy, T. Catal. Today 2012, 192, 123-129.

ACS Paragon Plus Environment

17

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

Page 18 of 24

(19) Ayoub, M., Khayoon, M.S., Abdullah, A.Z. Bioresour. Technol. 2012, 112, 308–312. (20) Cannilla, C., Bonura, G., Frusteri, L., Frusteri, F. Environ. Sci. Technol. 2014, 48, 60196026.

ACS Paragon Plus Environment

18

Page 19 of 24

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

TABLES AND FIGURES

Table 1. Central composite design with the three independent variables evaluated. Levels -α

-1

0

+1



(XT) Temperature (°C)

160

167

185

203

210

(XM.R.) Molar ratio

4:1

6:1

10:1

14:1

16:1

(XC.A.) Catalyst amount (g)

0

0.12

0.40

0.68

0.80

Table 2. Physicochemical properties for the solid materials. Catalyst

Acidity (mmol/g)

Surface (m2/g)

area Pore (cm3/g)

volume Pore size (Å)

Meso

Micro

Meso

Micro

Amberlyst 15

4.70

53*

0.40*

-

300*

-

Beta zeolite

0.29

434

0.05

0.15

15

6

Niobic acid

0.10

111

0.079

0.0042

30

9

*Provided by the supplier.

Table 3. Experimental conditions and results for glycerol etherification with ethanol over Amberlyst 15 obtained using central composite design.

ACS Paragon Plus Environment

19

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

Entry

T (°C)

MR

CA (g)

X (%)

Page 20 of 24

Y (%)

S (%) MEG

DEG

TEG

1

167

6

0.12

5

1

95

3

2

2

203

6

0.12

25

1

94

6

0

3

167

14

0.12

12

2

97

2

1

4

203

14

0.12

45

1

88

9

3

5

167

6

0.68

26

5

95

3

2

6

203

6

0.68

83

3

92

6

2

7

167

14

0.68

87

9

87

8

5

8

203

14

0.68

91

7

83

10

7

9

160

10

0.40

12

2

92

7

1

10

210

10

0.40

67

2

92

7

1

11

185

4

0.40

36

3

94

4

2

12

185

16

0.40

76

6

90

6

4

13

185

10

0.00

0

0

0

0

0

14

185

10

0.80

81

7

86

9

5

15

185

10

0.40

60

4

93

5

2

16

185

10

0.40

61

3

91

6

3

17

185

10

0.40

62

4

91

6

3

18

185

10

0.40

59

4

91

6

3

Note: T, temperature; MR, ethanol/glycerol molar ratio; CA, catalyst amount; X, conversion of glycerol; Y, yield to ethers; S, selectivity of MEG, DEG and TEG.

Table 4. Analysis of variance (ANOVA) for response surface quadratic model for, A) glycerol conversion; and B) Yield of ethers. Source

Sum squares

of Degrees freedom

of Mean square

F-value

P-value

ACS Paragon Plus Environment

20

Page 21 of 24

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

A Model

14524.55

9

1613.84

15.87

0.000000a

T

3064.95

1

3064.95

30.15

0.000580a

MR

1939.74

1

1939.74

19.08

0.002387a

CA

8245.13

1

8245.13

81.10

0.000018a

T2

449.85

1

449.85

4.42

0.068558a

MR2

4.52

1

4.52

0.04

0.838186b

CA2

391.86

1

391.86

3.85

0.085225a

T.MR

200.00

1

200.00

1.97

0.198331b

T.CA

8.00

1

8.00

0.08

0.786197b

MR.CA

220.50

1

220.50

2.17

0.179043b

Error

813.31

8

101.66

Model

100.02

9

11.11

37.03

0.000000a

T

2.08

1

2.08

6.93

0.029896a

MR

14.61

1

14.61

48.70

0.000115a

CA

69.60

1

69.60

232.00

0.000000a

T2

4.01

1

4.01

13.37

0.006411a

MR2

2.35

1

2.35

7.83

0.023255a

CA2

0.01

1

0.01

0.03

0.834784b

T.MR

0.12

1

0.12

0.40

0.536571b

T.CA

1.12

1.

1.12

3.73

0.088757a

MR.CA

6.12

1

6.12

20.40

0.001951a

Error

2.40

8

0.30

B

a

Significant at 90% confidence interval.

b

Not significant at 90% confidence interval.

ACS Paragon Plus Environment

21

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

Page 22 of 24

Table 5. Experimental glycerol conversion and yield of ethers of the reactions added to central composite design. Entry

T (°C)

MR

CA (g)

Experimental value S (%) (%) MEG

DEG

TEG

Glycerol conversion 19

213

7.94

0.83

77

92

7

1

20

211

9.82

0.87

89

82

12

6

21

209

11.69

0.90

91

89

7

4

22

205

15.44

0.98

88

86

10

4

23

184

19.44

0.64

65

88

8

4

24

190

15.71

0.61

82

87

8

5

25

167

7.08

1.58

5

96

2

2

26

163

13.91

1.87

9

95

3

2

27

161

17.32

2.01

12

90

6

4

28

159

20.73

2.15

8

94

3

3

Yield ethers

of

Note: T, temperature; MR, ethanol/glycerol molar ratio; CA, catalyst amount.

(a)

(b)

ACS Paragon Plus Environment

22

Page 23 of 24

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

Figure 1. Glycerol conversion (a) and yield of ethers (b) at different temperatures. Reactions conditions: Molar ratio ethanol/glycerol=10/1; Catalyst amount=0.40 g.

ACS Paragon Plus Environment

23

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)

Page 24 of 24

(b)

Figure 2. (a) Response surface for glycerol conversion against ethanol/glycerol molar ratio and temperature; (b) Response surface for glycerol conversion against catalyst amount and ethanol/glycerol molar ratio.

(a)

(b)

Figure 3. (a) Response surface for yield of ethers against catalyst amount and temperature; (b) Response surface for yield of ethers against catalyst amount and ethanol/glycerol molar ratio.

ACS Paragon Plus Environment

24