Isobaric Vapor Liquid Equilibrium (VLE) Data of the Systems n-Butanol

Vapor liquid equilibria (VLE) of the binary systems n-butanol + butyric acid and n-butanol + acetic acid were determined at two pressures, 26.65 and 5...
0 downloads 0 Views 87KB Size
120

J. Chem. Eng. Data 2001, 46, 120-124

Isobaric Vapor Liquid Equilibrium (VLE) Data of the Systems n-Butanol + Butyric Acid and n-Butanol + Acetic Acid Laura A. L. Mun ˜ oz and M. Alvina Kra1 henbu 1 hl* Departamento Processos Quı´micossFaculdade de Engenharia Quı´micasUNICAMP, Caixa Postal 6066, CEP 13083-970, Campinas, SP, Brazil

Vapor liquid equilibria (VLE) of the binary systems n-butanol + butyric acid and n-butanol + acetic acid were determined at two pressures, 26.65 and 53.33 kPa. The equipment used was a flow ebulliometer, which is ideal for reactive systems. The quality of the measured P-T-x-y data was verified by applying the thermodynamic consistency test of Van Ness and Fredenslund. The binary interaction parameters for the determination of the liquid-phase activity coefficients, represented by the models Wilson, UNIQUAC, and NRTL, were adjusted by using the maximum likelihood method. The nonideality of the vapor phase was considered by using a chemical theory with the correlation of Hayden and O’Connell for the calculation of the second virial coefficient and the prediction of the chemical equilibrium dimerization constant.

Introduction This paper reports the experimental vapor liquid equilibrium (VLE) data for the binary systems n-butanol + butyric acid and n-butanol + acetic acid at two pressures, 26.65 and 53.33 kPa. The quality of the measured data was verified by the Van Ness method with the modification proposed by Fredenslund et al.1 The data of vapor liquid equilibria are correlated by fitting the binary interaction parameters of the models for the liquid-phase activity coefficients, by means of the maximum likelihood method, according to the procedure proposed by Stragevich.2 Experimental Section Materials. All materials were supplied by Merck. The purities of the chemicals were greater than 99%, as indicated in Table 1, and were used directly without further purification. All the liquids were subjected to gas chromatography (GC) analysis, and no extra peaks appeared. Apparatus and Procedure. The VLE measurements were carried out with a flow ebulliometer developed by Vilı´m et al.4 and Ha´la et al.,5 as shown in Figure 1, connected in series to the pressure controller. The temperatures in the separation camera were read with a precision thermometer calibrated by Ever Ready Thermometer Co. All the equilibrium compositions were determined by GC by means of equipment supplied from “Instrumentos Cientı´ficos C. G. Ltda.”, model CG-35, equipped with a thermal conductivity detector. The GC column was carbowax with a 1.5 m length stainless steel tube and 3.175 mm diameter. Even more than three analyses were made for both liquid and vapor phase. The accuracy of temperature measurements was estimated to be (0.05 K. The accuracy of pressure control was within (0.133 kPa. The accuracy of the equilibrium composition measurements was within (0.006 mole fraction. Results and Discussion VLE had been measured at 26.65 and 53.33 kPa for n-butanol + butyric acid and n-butanol + acetic acid. The * Corresponding author. E-mail: [email protected].

Table 1. Normal Boiling Temperature and Purity of the Chemicals purity/% (GC) TBlit/Ka TBexp/K a

n-butanol

butyric acid

acetic acid

min 99.5 390.9 390.8

>99 437.2 437.5

min 99.8 391.1 390.9

Reid et al.3

Figure 1. Flow ebulliometer still.

results are shown in Tables 2 and 3. The activity coefficient γ of pure liquid i in a nonideal mixture at temperature T and pressure P was calculated according to

γi )

yiP φˆ Vi

[

exp sat xiPsat i φi

10.1021/je000033u CCC: $20.00 © 2001 American Chemical Society Published on Web 12/21/2000

]

VLi (P - Psat i ) RT

(1)

Journal of Chemical and Engineering Data, Vol. 46, No. 1, 2001 121 Table 2. Experimental Data for the System n-Butanol + Butyric Acid P ) 26.65 kPa

P ) 53.33 kPa

T/K

x1

y1

T/K

x1

y1

397.63 396.05 394.85 393.55 391.55 389.75 387.25 385.75 383.05 375.45 372.25 369.75 367.95 364.05 361.45 357.55

0.0000 0.0160 0.0472 0.0736 0.1069 0.1526 0.1991 0.2314 0.2927 0.4609 0.5476 0.6315 0.7061 0.8108 0.9072 1.0000

0.0000 0.0584 0.1609 0.2108 0.3017 0.3898 0.4933 0.5408 0.6411 0.8401 0.8877 0.9409 0.9606 0.9818 0.9959 1.0000

416.62 413.65 411.35 407.95 397.85 395.35 392.55 389.75 385.25 383.95 383.05 373.70

0.0000 0.0483 0.0908 0.1605 0.3751 0.4488 0.5059 0.5820 0.7094 0.7440 0.7701 1.0000

0.0000 0.1448 0.2546 0.4038 0.7209 0.7954 0.8592 0.9203 0.9704 0.9816 0.9908 1.0000

Table 3. Experimental Data for the System n-Butanol + Acetic Acid P ) 26.65 kPa

Figure 2. Thermodynamic consistency test for n-butanol + butyric acid at 26.65 kPa.

P ) 53.33 kPa

T/K

x1

y1

T/K

x1

y1

352.61 353.45 354.55 355.75 356.95 357.65 358.15 358.55 358.75 359.15 359.95 360.25 360.35 360.35 360.15 359.75 359.35 358.65 357.55

0.0000 0.0825 0.1460 0.2135 0.2832 0.3363 0.3675 0.4001 0.4206 0.4555 0.5643 0.6055 0.6741 0.7304 0.7639 0.8451 0.8827 0.9416 1.0000

0.0000 0.0435 0.0845 0.1266 0.1872 0.2271 0.2681 0.3074 0.3370 0.3761 0.5017 0.5779 0.6647 0.7424 0.7909 0.8849 0.9371 0.9767 1.0000

371.15 371.85 372.95 373.45 374.05 374.75 375.25 376.25 377.05 377.45 377.55 377.55 377.55 377.45 376.95 376.35 375.75 374.65 373.75

0.0000 0.0528 0.1118 0.1538 0.1901 0.2368 0.2835 0.3667 0.4806 0.5693 0.6122 0.6295 0.7053 0.7265 0.8218 0.8809 0.9202 0.9753 1.0000

0.0000 0.0358 0.0738 0.0991 0.1287 0.1675 0.1965 0.2941 0.3822 0.5202 0.6086 0.6283 0.7062 0.7547 0.8977 0.9507 0.9767 0.9992 1.0000

The fugacity coefficients of the vapor phase were calculated by means of the chemical theory with the correlation of Hayden and O’Connell6 for the calculation of the second virial coefficient and the prediction of the chemical equilibrium dimerization constant. The vapor pressures of the pure components were expressed by the Antoine equation in the form

log Psat i (mmHg) ) Ai -

Bi t (°C) + Ci

(2)

The experimental vapor pressure measured in this work are given in Table 4, and the Antoine constants of the chemicals were obtained by fitting these experimental data. In Table 10 are given all the parameters of the pure properties used in this investigation. A test was done to have an idea of ester formation. The method of the area percentage in the chromatographic analysis was used. In such a method it is not required to specify data of the pure components involved in the system; the area of the chromatogram is just taken into consideration. In Tables 5 and 6 the results are shown and it is verified that there was a small amount of ester formed in both systems. In general it was less than 0.05% in the liquid and vapor phases for the both pressures. The amount of ester formed is not significant, and it can be ignored in the calculations.

Figure 3. Thermodynamic consistency test for n-butanol + butyric acid at 53.33 kPa. Table 4. Experimental Vapor Pressures for the Pure Components acetic acid

n-butanol

butyric acid

P/kPa

T/K

P/kPa

T/K

P/kPa

T/K

20.00 25.72 26.67 28.04 30.00 37.08 50.00 53.84 60.00 70.00 74.65 79.99

345.15 351.45 352.35 353.65 355.35 361.05 369.45 371.55 374.75 379.35 381.35 383.55

26.66 27.73 30.00 34.16 39.97 40.47 50.00 53.45 60.00 66.63 80.00 86.81 95.00

357.50 358.20 360.10 363.10 366.80 367.00 372.20 373.90 376.80 379.50 384.40 386.60 389.10

20.82 23.31 26.13 28.89 32.06 37.49 42.06 47.62 50.00 54.26 58.13 60.00 63.85 66.66 70.00 75.00 79.99

391.50 394.30 397.50 399.80 402.60 407.10 409.90 413.10 415.20 417.50 419.30 420.50 422.00 423.30 425.20 427.30 428.90

All the binary systems were tested for their thermodynamic consistency by means of the Van Ness method modified by Fredenslund et al.1 The mean deviation between experimental and calculated y values should be lower than 1 × 10-2. In this method, the heat of mixing is neglected.

122

Journal of Chemical and Engineering Data, Vol. 46, No. 1, 2001

Table 5. Chromatographic Analyses in Percent Area for the System n-Butanol + Butyric Acid P ) 26.65 kPa liquid phase

P ) 53.33 kPa vapor phase

n-butanol

butyric acid

n-butyl butyrate

0.00 1.55 4.59 7.16 10.42 14.89 19.29 22.63 28.68 45.53 54.05 62.48 70.02 80.62 90.46 100.00

100.00 98.40 95.38 92.82 89.55 85.07 80.69 77.35 71.31 54.46 45.93 37.50 29.98 19.35 9.52 0.00

0.00 0.05 0.03 0.02 0.03 0.03 0.02 0.02 0.01 0.01 0.02 0.02 0.00 0.03 0.01 0.00

liquid phase

n-butanol

butyric acid

n-butyl butyrate

0.00 5.68 15.68 20.61 29.57 38.31 48.62 53.37 60.71 83.69 76.01 93.91 95.92 99.98 99.99 100.00

100.00 94.19 84.30 79.38 70.42 61.68 51.37 46.62 39.28 16.29 23.93 6.06 4.05 0.00 0.00 0.00

0.00 0.13 0.02 0.01 0.01 0.01 0.00 0.01 0.02 0.02 0.07 0.03 0.03 0.02 0.01 0.00

vapor phase

n-butanol

butyric acid

n-butyl butyrate

n-butanol

butyric acid

n-butyl butyrate

0.00 6.16 11.45 15.64 36.98 44.17 49.88 59.54 70.35 73.84 76.51 100.00

100.00 93.80 88.52 84.31 63.00 55.81 50.12 40.46 29.65 26.14 26.14 0.00

0.00 0.04 0.03 0.05 0.02 0.02 0.00 0.00 0.00 0.01 0.00 0.00

0.00 18.00 30.68 39.69 71.63 79.06 85.58 91.82 96.96 98.10 99.06 100.00

100.00 82.00 69.31 60.28 28.36 20.92 14.42 8.18 3.04 1.89 0.94 0.00

0.00 0.00 0.02 0.03 0.01 0.02 0.00 0.00 0.00 0.01 0.00 0.00

Table 6. Chromatographic Analyses in Percent Area for the System n-Butanol + Acetic Acid P ) 26.65 kPa liquid phase

P ) 53.33 kPa vapor phase

liquid phase

vapor phase

n-butanol

acetic acid

n-butyl acetate

n-butanol

acetic acid

n-butyl acetate

n-butanol

acetic acid

n-butyl acetate

n-butanol

acetic acid

n-butyl acetate

0.00 10.12 17.63 25.36 33.07 38.85 42.61 45.49 47.61 51.14 61.84 65.76 72.13 77.21 80.21 87.21 90.39 95.27 100.00

100.00 89.88 82.37 74.63 66.90 61.14 57.38 54.50 52.39 48.85 38.15 34.23 27.86 22.77 19.79 12.78 9.60 4.72 0.00

0.00 0.00 0.00 0.01 0.03 0.01 0.01 0.01 0.00 0.01 0.01 0.01 0.01 0.02 0.00 0.02 0.01 0.01 0.00

0.00 5.83 10.35 15.35 22.27 26.89 31.43 35.71 38.88 42.99 55.74 63.14 71.23 78.21 82.53 90.55 94.88 98.09 100.00

100.00 94.62 89.64 84.64 77.61 73.10 68.56 64.28 61.11 56.98 44.24 36.84 28.72 21.76 17.43 9.41 5.09 1.87 0.00

0.00 0.00 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.02 0.02 0.02 0.05 0.04 0.04 0.04 0.03 0.04 0.00

0.00 6.52 13.98 18.54 22.70 27.97 33.01 42.01 53.65 62.32 66.39 68.01 74.97 76.87 82.36 87.36 91.03 97.96 100.00

100.00 93.47 86.02 81.46 77.30 72.03 66.99 57.99 46.34 37.68 33.60 31.99 25.03 23.12 17.64 12.64 8.91 1.98 0.00

0.00 0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.01 0.00 0.00 0.00 0.00 0.00 0.06 0.05 0.00

0.00 4.44 9.07 12.10 15.60 20.12 23.44 34.27 43.62 42.41 66.04 67.90 75.05 79.36 89.14 94.14 96.78 98.84 100.00

100.00 95.56 90.93 87.90 84.36 79.87 76.55 65.73 56.35 57.58 33.94 32.09 24.94 20.61 10.86 5.86 3.16 1.15 0.00

0.00 0.00 0.00 0.00 0.04 0.01 0.01 0.00 0.02 0.01 0.02 0.01 0.01 0.02 0.00 0.00 0.06 0.01 0.00

Table 7. Mean Deviations of Y from the Consistency Test system

|∆y| P26.65

|∆y| P53.33

n-butanol + butyric acid n-butanol + acetic acid

0.0101 0.0081

0.0112 0.0099

Table 7 presents the obtained mean deviations for each set. The graphic representation of the consistency results together with the experimental data for each pressure is in Figures 2-5. It can be seen that there is a good distribution between positive and negative differences of ∆y, indicating no systematic errors. The VLE experimental data of all binary systems were used to obtain the interaction parameters of three activity coefficient models: Wilson, NRTL, and UNIQUAC.3 The fitted parameters were obtained through an algorithm based on the principle of the maximum likelihood method developed by Stragevitch.2 The fitted parameters together with the mean values of the absolute deviations of the equilibrium temperature (∆T) and the mole fraction of the vapor phase (∆y) are shown in Tables 8 and 9. For the NRTL model, the parameter Rij was considered to be a constant equal to 0.2. It was seen that the value of this parameter causes no significant influence on the results. The calculated results for the

Table 8. Correlation Parameters and Absolute Mean Deviations in Boiling Temperatures and Vapor-Phase Mole Fractions for the System n-Butanol + Butyric Acid model

A12a

A21a

|∆y|

|∆T|

Wilson UNIQUAC NRTL

P ) 26.65 kPa -426.76 506.20 307.95 -264.96 -689.15 854.56

0.0137 0.0125 0.0098

0.62 0.71 1.37

Wilson UNIQUAC NRTL

P ) 53.33 kPa -503.27 700.97 372.04 -306.93 990.78 -802.15

0.0207 0.0192 0.0170

0.59 0.74 1.38

Wilson UNIQUAC NRTL

Global (26.65 and 53.33 kPa) -440.27 504.74 0.0173 295.88 -265.43 0.0161 -764.97 988.39 0.0138

0.60 0.71 1.41

a The binary adjustable parameters for various models are as follows: Wilson, Aij ) (λij - λii); NRTL, Aij ) (gij - gjj); UNIQUAC, Aij ) (uij - ujj).

Wilson, UNIQUAC, and NRTL models are illustrated in Figures 6-9. It can be seen from Tables 8 and 9 that all the adjustments are equivalent. The resulting mean deviations of the

Journal of Chemical and Engineering Data, Vol. 46, No. 1, 2001 123 Table 9. Correlation Parameters and Absolute Mean Deviations in Boiling Temperatures and Vapor-Phase Mole Fractions for the System n-Butanol + Acetic Acid model

A12a

A21a

|∆y|

|∆T|

Wilson UNIQUAC NRTL

P ) 26.65kPa 566.02 -191.96 -234.83 397.44 -588.54 923.62

0.0044 0.0043 0.0043

0.46 0.45 0.34

Wilson UNIQUAC NRTL

P ) 53.33 kPa 733.69 -222.73 -273.77 498.38 -705.70 1208.1

0.0090 0.0088 0.0071

0.65 0.64 0.50

Wilson UNIQUAC NRTL

Global (26.65 and 53.33 kPa) 656.82 -208.80 0.0071 -254.31 446.40 0.0071 -641.65 1044.5 0.0066

0.51 0.54 0.39

a The binary adjustable parameters for various models are as follows: Wilson, Aij ) (λij - λii); NRTL, Aij ) (gij - gjj); UNIQUAC, Aij ) (uij - ujj).

Figure 5. Thermodynamic consistency test for n-butanol + acetic acid at 53.33 kPa.

Table 10. Pure Component Properties Used in This Work parameter

n-butanol

butyric acid

acetic acid

Aia Bia Cia Tcb/K Pcb/bar Vcb/cm3.mol-1 $b Vb/cm3‚mol-1 RDc/Å ETAc Rd q (q ) q′)d

7.32672 1302.39 175.114 563.1 44.2 275.0 0.593 91.97 3.225 2.20 3.4543 3.0520

7.76035 1780.62 201.590 628.0 52.7 290.0 0.683 92.43 3.134 4.50 3.5512 3.1520

9.00908 2749.36 330.948 592.7 57.9 171.0 0.447 57.54 2.595 4.50 2.2024 2.0720

a Adjusted with the collected experimental data: t/°C and P/mmHg. b Reid et al.3 c Fredenslund et al.1 d Gmehling et al.7

Figure 6. System n-butanol + butyric acid at 26.65 kPa.

Figure 4. Thermodynamic consistency test for n-butanol + acetic acid at 26.65 kPa. Figure 7. System n-butanol + butyric acid at 53.33 kPa.

variables can be attributed to experimental errors and ester formation that was neglected. However, when Figures 6-9 are analyzed, it appears that a significant displacement appears for the system n-butanol + acetic acid between the experimental data and the calculated curve. This is due to the smaller temperature interval (∼8 K) for this system compared to that for the other system (∼40 K).

Conclusion According to the results of the consistency test, it can be said that the reported data are good enough to represent the systems, since they are right at the cutoff at which Fredenslund et al.1 consider data to be consistent. Analyses of experimental data of all binary systems using the Wilson, NRTL, and UNIQUAC equations show that all

124

Journal of Chemical and Engineering Data, Vol. 46, No. 1, 2001

Figure 8. System n-butanol + acetic acid at 26.65 kPa.

C ) Antoine coefficient ETA ) association parameter gij-gjj ) NRTL binary parameters of components i and j P ) pressure Pc ) critical pressure Psat ) saturated vapor pressure q ) UNIQUAC area parameter r ) UNIQUAC volumetric parameter R ) ideal gas constant RD ) radius of gyration Vc ) critical volume VL ) liquid molar volume T ) temperature Tb ) boiling temperature Tc ) critical temperature uij-ujj ) UNIQUAC binary parameters of components i and j x ) mole fraction of liquid phase y ) mole fraction of vapor phase Rij ) NRTL parameter nonrandom γi ) activity coefficient of component i λij-λii ) Wilson binary parameters of components i and j φsat ) saturated vapor-phase fugacity coefficient $ ) acentric factor Literature Cited

Figure 9. System n-butanol + acetic acid at 53.33 kPa.

three models were generally satisfactory. It can be assumed that no reaction between the compounds has taken place in the systems, since the amount of ester formed can be neglected.

(1) Fredenslund, A.; Gmehling, J.; Rasmussen, P. Vapor-liquid equilibria using UNIFACsA group contribution method; Elsevier Scientific Publishing Company: New York, 1977. (2) Stragevitch, L. Equilı´brio lı´quido-lı´quido de misturas de na˜o eletro´litos. Ph.D. Thesis, Faculdade de Engenharia Quı´micas UNICAMPsCampinas, Brazil, 1997. (3) Reid, R. C.; Prausnitz, J. M.; Poling, B. E. The properties of gases and liquids; McGraw-Hill: New York, 1987. (4) Vilı´m, O.; Ha´la, E.; Fried, V.; et al. Collect. Czech. Chem. Commum. 1954, 19, 1330-1334. (5) Ha´la, E.; Pick, J.; Fried, V.; et al. Vapour-Liquid Equilibrium; Oxford: Pergamon: 1967. (6) Hayden, J. G.; O’Connell, J. P. Ind. Eng. Chem. Process Des. Dev. 1975, 14 (3), 209-216. (7) Gmehling, J.; Onken, U. Vapor-liquid equilibrium data collection chemistry data series; Dechema: Frankfurt, 1977; pp 1-7.

Acknowledgment The authors thank Mr. D. R. Queiroz for collaboration on the measurements. Nomenclature

A ) Antoine coefficient B ) Antoine coefficient

Received for review January 28, 2000. Accepted October 7, 2000. The authors gratefully acknowledge the financial support of CNPqsConselho Nacional de Desenvolvimento Cientı´fico e Tecnolo´gico (Brasil), FAPESPsFundac¸ a˜o de Amparo a` Pesquisa do Estado de Sa˜o Paulo (Brasil), and FAEPsFundo de Apoio ao Ensino e a` Pesquisa (UNICAMP, Brasil).

JE000033U