Fluid Phase Equilibria 221 (2004) 1 6 Vapor liquid equilibria for the binary system 2,2 dimethylbutane + 1,1 dimethylpropyl methyl ether (TAME) at 298.15, 318.15, and 338.15 K Armando del Río a, Baudilio Coto b, Juan A.R. Renuncio a, Concepción Pando a, a Departamento de Química Física I, Universidad Complutense, E-28040 Madrid, Spain b ESCET, Universidad Rey Juan Carlos, E-28933 Móstoles, Madrid, Spain Received 16 January 2004; accepted 2 April 2004 Available online 7 June 2004 Abstract Vapor liquid equilibrium (VLE) data are reported for the binary mixtures formed by 2,2-dimethylbutane and the branched ether 1,1- dimethylpropyl methyl ether (tert-amyl methyl ether or TAME). A Gibbs-van Ness type apparatus was used to obtain total vapor pressure measurements as a function of composition at 298.15, 318.15, and 338.15 K. Deviations from Raoult s law are positive and very small. These VLE data are analyzed together with excess enthalpy (Hm E ) data previously reported at 298.15 K using the UNIQUAC model. The modified UNIFAC (Dortmund) model is used to predict VLE data. 2004 Elsevier B.V. All rights reserved. Keywords: Binary liquid mixtures; Vapor liquid equilibria; Gibbs energy; Branched ether 1. Introduction In recent years, the interest in oxygenated gasoline has resulted in a large number of studies on the thermophysical properties of branched ethers and their binary and ternary mixtures. The purpose of this paper is to report phase equilibria for the binary system formed by 2,2-dimethylbutane and 1,1-dimethylpropyl methyl ether (tert-amyl methyl ether or TAME). 2. Experimental section 2,2-Dimethylbutane was purchased from Sigma, with a purity >99 mol%. TAME (Fluka 97 mol% purity) was refluxed over 0.3 nm molecular sieves for several hours and then fractionated. The middle distillate used in the present work (approximately 50% of the initial amount) had a purity Abbreviations: TAME, 1,1-dimethyl-propyl methyl ether (tert-amyl methyl ether); VLE, vapor liquid equilibrium Corresponding author. Fax: +34-91-3944135. E-mail address: pando@quim.ucm.es (C. Pando). better than 99.6%, as measured by a gas chromatographic analysis. Chemicals were handled under a dry nitrogen atmosphere and were degassed by reflux distillation for several hours following a procedure described elsewhere [1]. VLE data were measured using a Gibbs-van Ness type static apparatus [2] which was described previously by Coto et al. [1]. Binary liquid solutions of known composition were prepared in a test cell by volumetric injection of degassed liquids using calibrated pistons. The equilibrium cell volume is ca. 100 ml. Values for the radii of the left and right pistons are 6.000 ± 0.001 and 6.006 ± 0.001 mm, respectively. The liquid inside the piston-injector is kept at a pressure higher than the cell pressure to avoid a vapor phase. The volume injected of each component can be determined from the position of the piston which is registered by a length gauge (ACU-RITE A R/1, 2 m accuracy). Therefore, the accuracy of the injected volume is ±0.0002 ml. The amount injected (in mol) is calculated using the compressibility of the liquid and its density at atmospheric pressure. The accuracy of the mole fraction is estimated to be ca. 0.0001 in the dilute regions and ca. 0.0003 in the middle of the concentration range. Cell and piston-injectors were immersed in a water bath in which temperature was controlled within ±0.002 K. 0378-3812/$ see front matter 2004 Elsevier B.V. All rights reserved. doi:10.1016/j.fluid.2004.04.001
2 A. del Río et al. / Fluid Phase Equilibria 221 (2004) 1 6 Table 1 Properties of pure components used in this study T(K) V m (cm 3 mol 1 ) κ T 10 9 (Pa 1 ) B (cm 3 mol 1 ) p (kpa) p (kpa) (literature values range) 2,2-Dimethylbutane 298.15 136.83 [3] 2.01 [4] 1889.0 42.87 42.51 42.55 [5,6] 318.15 141.05 [3] 2.01 [4] 1566.0 86.83 86.79 86.85 [5,6] 338.15 145.70 [3] 2.01 [4] 1127.8 160.93 160.92 160.93 [5,6] TAME 298.15 133.44 [7] 1.26 [7] 2285.0 10.02 10.01 10.10 [8 10] 318.15 135.86 [7] 1.80 [7] 1875.0 24.02 23.99 24.13 [9,10] 338.15 138.56 [7] 2.34 [7] 1714.0 50.83 50.82 50.88 [11,12] The temperature was measured with a quartz thermometer, Testo 781, with an accuracy of 0.01 K. The total vapor pressure was measured when phase equilibrium was reached and maintained for at least 15 min using a differential MKS Baratron pressure gauge with a resolution of 0.08% of the reading. Pressure accuracy is estimated to be 0.01 kpa. Table 1 lists values for the molar volumes, isothermal compressibilities and second virial coefficients of pure components used in the data reduction and values for the vapor pressures of 2,2-dimethylbutane and TAME obtained in this study and those found in the literature at the three temperatures studied. The source of molar volumes and isothermal compressibilities is given in Table 1. Second virial coefficients were calculated by means of the Hayden and O Connell method [13]. The vapor pressures measured in this work are in good agreement with literature values and this is considered an indication of the chemicals purity. azeotrope), and the slopes of the pressure versus x 1 diagram are extremely high in the dilute regions. VLE data obtained were shown to be consistent with previous literature data. More details about this comparison ant the data-reduction procedure are found in Ref. [1]. Table 2 lists values for the 2,2-dimethylbutane liquid composition, x 1, and the total pressure, p, at the three temperatures studied, together with the calculated total pressure, p calc, the hydrocarbon vapor composition, y 1,calc, the molar excess Gibbs energy, G E m, for 2,2-dimethylbutane(1) + TAME(2), the logarithm of the activity coefficients γ 1 and γ 2, the coefficients P i for the G E m representation by Eq. (1), and their uncertainties, and the standard deviations between experimental and calculated values of x 1, σ x, and p, σ p. Fig. 1 shows plots of the total pressure against liquid and vapor compositions for the 2,2-dimethylbutane + TAME system at the three temperatures studied. Both 3. Results and discussion VLE measurements for 2,2-dimethylbutane + TAME were carried out at 298.15, 318.15, and 338.15 K. Results were analyzed using a modified Barker method and the maximum likelihood principle [14]. The temperature, T, and the amount of substance for components 1 and 2 were considered to be the independent variables in the data reduction. The material balance was included to take into account the amounts of substance present in the vapor phase. The molar excess Gibbs energy, G E m, of the liquid phase was described by the Redlich Kister equation 3 G E m = RTx 1(1 x 1 ) P i (2x 1 1) i (1) i=0 where P i are adjustable parameters, and x 1 is the mole fraction of 2,2-dimethylbutane. The vapor phase is described using the virial equation. In order to test the quality of the experimental data obtained using the experimental technique and the data-reduction procedure described above, VLE data for the ethanol + cyclohexane system were previously measured at 298.15 K by Coto et al. [1]. This system shows very high positive deviations from ideality (including an Fig. 1. VLE data for the 2,2-dimethylbutane(1) + TAME(2) system: ( ) 298.15 K; ( ) 318.15 K; ( ) 338.15 K; ( ) calculated values using Eq. (1).
A. del Río et al. / Fluid Phase Equilibria 221 (2004) 1 6 3 Table 2 Vapor liquid equilibrium data and coefficients and standard deviations for G E m representation by Eq. (1) for 2,2-dimethylbutane(1) + TAME(2) at 298.15, 318.15, and 338.15 K x 1 p (kpa) p calc (kpa) y 1,calc G E m (J mol 1 ) ln γ 1 ln γ 2 T = 298.15 K a 0.0000 10.02 10.02 0.0000 0 0.0000 0.0518 12.07 12.06 0.2103 22 0.1673 0.0002 0.0986 13.87 13.87 0.3458 41 0.1591 0.0008 0.1509 15.84 15.84 0.4586 60 0.1475 0.0025 0.2014 17.69 17.70 0.5422 78 0.1352 0.0051 0.2484 19.38 19.38 0.6045 92 0.1237 0.0085 0.2964 21.06 21.06 0.6573 105 0.1123 0.0127 0.3466 22.79 22.78 0.7037 116 0.1012 0.0180 0.3989 24.48 24.54 0.7449 126 0.0907 0.0243 0.4460 26.16 26.10 0.7772 133 0.0820 0.0306 0.4954 27.77 27.72 0.8072 138 0.0736 0.0381 0.5491 29.47 29.45 0.8361 142 0.0649 0.0476 0.6250 31.77 31.84 0.8714 142 0.0529 0.0647 0.6752 33.31 33.37 0.8920 139 0.0448 0.0798 0.7333 35.04 35.09 0.9135 132 0.0351 0.1031 0.7944 36.83 36.84 0.9341 118 0.0246 0.1374 0.8412 38.21 38.14 0.9490 103 0.0166 0.1732 0.8928 39.67 39.55 0.9650 79 0.0087 0.2251 0.9407 40.98 40.84 0.9800 50 0.0031 0.2885 1.0000 42.47 42.47 1.0000 0 0.0000 T = 318.15 K b 0.0000 24.02 24.02 0.0000 0 0.0000 0.0472 27.49 27.52 0.1660 19 0.1466 0.0004 0.0906 30.64 30.65 0.2828 35 0.1338 0.0013 0.1527 35.04 35.02 0.4119 55 0.1176 0.0036 0.1943 37.91 37.88 0.4809 67 0.1080 0.0056 0.2441 41.28 41.25 0.5503 80 0.0978 0.0084 0.2987 44.89 44.88 0.6140 92 0.0878 0.0122 0.3461 47.96 47.99 0.6614 101 0.0799 0.0159 0.3954 51.10 51.19 0.7045 108 0.0724 0.0203 0.4501 54.55 54.70 0.7463 115 0.0645 0.0261 0.4950 57.65 57.55 0.7768 119 0.0582 0.0317 0.5439 60.81 60.61 0.8068 121 0.0515 0.0390 0.5889 63.40 63.36 0.8317 122 0.0453 0.0471 0.6342 66.10 66.10 0.8547 121 0.0390 0.0570 0.6880 69.24 69.29 0.8797 117 0.0315 0.0718 0.7426 72.39 72.45 0.9030 109 0.0239 0.0910 0.7915 75.19 75.23 0.9226 98 0.0173 0.1127 0.8417 77.98 78.04 0.9416 84 0.0110 0.1408 0.8934 80.91 80.90 0.9606 63 0.0056 0.1768 0.9468 83.96 83.86 0.9800 35 0.0015 0.2236 1.0000 86.83 86.83 1.0000 0 0.0000 T = 338.15 K c 0.0000 50.83 50.83 0.0000 0 0.0000 0.0378 55.68 55.68 0.1185 15 0.1284 0.0003 0.0750 60.26 60.27 0.2141 27 0.1131 0.0013 0.1119 64.85 64.73 0.2943 38 0.1003 0.0026 0.1518 69.34 69.43 0.3687 48 0.0887 0.0043 0.2024 75.31 75.30 0.4492 59 0.0770 0.0068 0.2461 80.10 80.32 0.5093 67 0.0690 0.0091 0.2938 85.86 85.77 0.5669 74 0.0621 0.0117 0.3381 91.03 90.81 0.6144 80 0.0570 0.0140 0.4043 98.16 98.30 0.6768 87 0.0508 0.0177 0.4699 105.61 105.71 0.7302 93 0.0456 0.0217 0.5233 111.65 111.69 0.7685 96 0.0415 0.0258 0.5810 118.86 118.13 0.8056 97 0.0367 0.0317 0.6292 123.12 123.32 0.8336 97 0.0324 0.0383 0.6759 127.92 128.32 0.8586 96 0.0278 0.0470 0.7326 133.94 134.27 0.8865 91 0.0218 0.0613 0.7872 139.80 139.88 0.9113 83 0.0158 0.0804
4 A. del Río et al. / Fluid Phase Equilibria 221 (2004) 1 6 Table 2 (Continued ) x 1 p (kpa) p calc (kpa) y 1,calc G E m (J mol 1 ) ln γ 1 ln γ 2 0.8371 144.98 144.91 0.9327 72 0.0105 0.1036 0.8949 151.12 150.65 0.9565 54 0.0050 0.1394 0.9455 156.40 155.63 0.9772 32 0.0015 0.1801 1.0000 160.93 160.93 1.0000 0 0.0000 a P 0 = 0.223 ± 0.002; P 1 = 0.068 ± 0.005; P 2 = 0.059 ± 0.007; P 3 = 0.04 ± 0.01; σ x = 0.0001; σ p = 56 Pa. b P 0 = 0.180 ± 0.002; P 1 = 0.050 ± 0.004; P 2 = 0.042 ± 0.006; P 3 = 0.01 ± 0.01; σ x = 0.0001; σ p = 71 Pa. c P 0 = 0.134 ± 0.004; P 1 = 0.039 ± 0.009; P 2 = 0.057 ± 0.001; P 3 = 0.01 ± 0.03; σ x = 0.0002; σ p = 33 Pa. experimental values and those calculated using Eq. (1) are shown. Deviations from Raoult s law are positive and small. Since the mixtures behave almost ideally, the correction for real behavior of the vapor phase is not necessary. Fig. 2 shows plots of G E m values against x 1 for the three temperatures studied. Values for the excess Gibbs energy of the 2,2-dimethylbutane + TAME mixtures are positive and decrease as temperature increases for a given mole fraction. Maximum G E m values range from 142 J mol 1 at 298.15 K to 97 J mol 1 at 338.15 K. These maxima appear at approximately the same mole fraction (x 1 0.6) for the three temperatures studied. A similar behavior was reported for hexane + TAME mixtures [15]. Figs. 3 5 show plots of ln γ 1 and ln γ 2 against x 1 at 298.15, 318.15, and 338.15 K, respectively. As could be expected for an almost ideal system, values for the activity coefficients are close to one and show a smooth variation with composition. The temperature effect on the activity coefficients is small. The limiting activity coefficients were calculated using the method of Maher and Smith [16]. Values of 1.17, 1.16 and 1.15 were obtained for γ 1 at 298.15, 318.15, and 338.15 K, respectively. Values of 1.91, 1.89 and 1.87 were obtained for γ 2 at 298.15, 318.15, and 338.15 K, respectively. Fig. 6 shows plots of G E m /RT against temperature in the 298.15 338.15 K range for x 1 = 0.25, 0.50 and 0.75, respectively. G E m /RT decreases linearly with temperature; the temperature effect on G E m is in agreement with the moderately endothermic values of the molar excess enthalpy (Hm E) reported at 298.15 K for these mixtures: Knezevic-Stevanovic et al. [17] reported a maximum value of 224.1 J mol 1 occurring at x 1 = 0.500. Values for Hm E calculated from the slopes of the G E m /RT against T plots at 298.15 K are more endothermic than the experimental excess enthalpies: a value of 411 J mol 1 is obtained for x 1 = 0.500. More accurate values for Hm E at this temperature could be obtained if similar plots could be drawn for a temperature range beginning at a temperature lower than 298.15 K. The UNIQUAC model [18] was used to correlate the VLE data reported in this paper at 298.15, 318.15, and 338.15 K and the Hm E data previously obtained at 298.15 K. Uncertainties in temperature, pressure and composition were considered to be smaller than those of vapor pressures and excess enthalpies. The interaction parameters of the UNIQUAC model, A ji, were considered to vary linearly with temperature according to A ji = A ji,1 + A ji,2 (T 298.15) (2) A least-square procedure was used to minimize deviations between experimental and calculated vapor pressure and enthalpies. The resulting values for the interaction parameters Fig. 2. G E m against x 1 values for the 2,2-dimethylbutane(1) + TAME(2) system: ( ) 298.15 K; ( ) 318.15 K; ( ) 338.15 K; ( ) calculated values using Eq. (1). Fig. 3. Activity coefficients against x 1 values for the 2,2-dimethylbutane(1) + TAME(2) system at 298.15 K: ( ) lnγ 1 ;( ) lnγ 2.
A. del Río et al. / Fluid Phase Equilibria 221 (2004) 1 6 5 Table 3 Vapor liquid equilibrium for 2,2-dimethylbutane(1) + TAME(2): correlation using the UNIQUAC model and prediction using the modified UNI- FAC (Dortmund) model a T (K) UNIQUAC Modified UNIFAC (Dortmund) σ p (kpa) σ rel σ p (kpa) σ rel 298.15 0.21 0.5 0.24 0.6 318.15 0.13 0.2 0.27 0.3 338.15 0.39 0.2 0.61 0.4 a σ p is the standard deviations between experimental and calculated vapor pressures and σ rel is the percent ratio of σ p and the maximum value of the vapor pressure. Fig. 4. Activity coefficients against x 1 values for the 2,2-dimethylbutane(1) + TAME(2) system at 318.15 K: ( ) lnγ 1 ;( ) lnγ 2. Fig. 5. Activity coefficients against x 1 values for the 2,2-dimethylbutane(1) + TAME(2) system at 338.15 K: ( ) lnγ 1 ;( ) lnγ 2. are: A 12,1 = 140.3K, A 12,2 = 0.6580, A 21,1 = 101.7K, A 21,2 = 0.5410. A value of 1.6 J mol 1 was obtained for the standard deviation between experimental and calculated molar excess enthalpies, σ H. This leads to a value of 0.7 for the percent ratio of σ H and the maximum molar excess enthalpy value. Values for the standard deviation between experimental and calculated vapor pressures, σ p, and the percent ratio of σ p and the maximum value of vapor pressure, σ rel, are listed in Table 3. The correlation is accurate: values for σ p are lower than 0.4 kpa and values of σ rel are lower than 0.5. The modified UNIFAC (Dortmund) [19,20] model was used to predict VLE data for the 2,2-dimethylbutane + TAME mixtures at 298.15, 318.15, and 338.15 K. Values for σ p and σ rel are listed in Table 3. These values are slightly higher than those obtained using the UNIQUAC model. We may conclude that the UNIQUAC model provides a simultaneous description of VLE and Hm E data for the 2,2-dimethylbutane + TAME system and that vapor pressures are accurately predicted for this system using the modified UNIFAC (Dortmund) model. List of symbols A ji UNIQUAC model interaction parameters B second virial coefficient G Gibbs energy H enthalpy p vapor pressure P i Redlich Kister equation parameter V volume x mole fraction in the liquid phase y mole fraction in the vapor phase Greek symbols γ activity coefficient κ isothermal compressibility σ standard deviation Superscripts E excess property Fig. 6. G E m /RT against temperature values for the 2,2-dimethylbutane(1) + TAME(2) system: ( ) x 1 = 0.25; ( ) x 1 = 0.50; ( ) x 1 = 0.75. Subscripts calc calculated i, j components of binary systems
6 A. del Río et al. / Fluid Phase Equilibria 221 (2004) 1 6 H excess enthalpy m molar property p vapor pressure rel relative to maximum value x liquid composition of component 1 Acknowledgements This work was funded by the Spanish Ministry of Science and Technology, PPQ-2002-04143-C02-02. A. del Río acknowledges the Comunidad Autonoma de Madrid for its support through a predoctoral grant. References [1] B. Coto, C. Pando, R.G. Rubio, J.A.R. Renuncio, J. Chem. Soc., Faraday Trans. 91 (1995) 273 278. [2] R.E. Gibbs, H.C. van Ness, Ind. Eng. Chem. Fundam. 11 (1972) 410 413. [3] TRC Thermodynamic Tables, Thermodynamic Research Center, Texas A&M University System, College Station, TX. [4] D.-Y. Peng, G.C. Benson, B.C.-Y. Lu, J. Solut. Chem. 28 (1999) 505 519. [5] J. Gmehling, U. Onken, W. Arlt, Vapor liquid equilibrium data collection, DECHEMA, Chem. Data Ser. I, Part 6b DECHEMA, Frankfurt Main, Germany, 1980. [6] R.C. Reid, J.M. Prausnitz, B.E. Poling, The Properties of Gases and Liquids, McGraw-Hill, Singapore, 1988. [7] U.P. Govender, T.M. Letcher, J. Chem. Eng. Data 41 (1996) 147 150. [8] B. Coto, F. Mössner, C. Pando, R.G. Rubio, J.A.R. Renuncio, J. Chem. Soc., Faraday Trans. 92 (1996) 4435 4440. [9] I. Cervenkova, T. Boublik, J. Chem. Eng. Data 29 (1984) 425 427. [10] M. Palczewska-Tulinska, D. Wyrzykowska-Stankiewicz, J. Cholinski, K. Zieborak, Fluid Phase Equilibria. 54 (1990) 57 68. [11] F. Mössner, B. Coto, C. Pando, R.G. Rubio, J.A.R. Renuncio, J. Chem. Eng. Data 41 (1996) 537 542. [12] M.A. Krähenbühl, J. Gmehling, J. Chem. Eng. Data 39 (1994) 759 762. [13] J.G. Hayden, J.P. O Connell, Ind. Eng. Chem. Process Des. Dev. 14 (1975) 209 216. [14] R.G. Rubio, J.A.R. Renuncio, M. Diaz Peña, Fluid Phase Equilibria. 12 (1983) 217 234. [15] A. del Río, B. Coto, C. Pando, J.A.R. Renuncio, Ind. Eng. Chem. Res. 41 (2002) 1364 1369. [16] P.J. Maher, B.D. Smith, Ind. Eng. Chem. Fundam. 18 (1979) 354 357. [17] A. Knezevic-Stevanovic, Z.-L. Jin, G.C. Benson, B.C.-Y. Lu, J. Chem. Eng. Data 40 (1995) 340 342. [18] D.S. Abrams, J.M. Prausnitz, AIChE J. 21 (1975) 116 128. [19] U. Weidlich, J. Gmehling, Ind. Eng. Chem. Res. 26 (1987) 1372 1381. [20] J. Gmehling, L. Jiding, M. Schiller, Ind. Eng. Chem. Res. 32 (1993) 178 193.