CC BY
79
CHEMISTRY AND ECOLOGY
ENTHALPY OF PHASE TRANSITION AND PREDICTION OF PHASE EQUILIBRIA IN SYSTEMS OF GLYCOLS AND GLYCOL ETHERS
Z. N. Esinaa’*, A. M. Miroshnikovb, and M. R. Korchuganovaa
aKemerovo State University, st. Krasnaya 6, Kemerovo, 650043 Russia,
*phone/fax: +7 (3842) 54-27-70, e-mail: [email protected] bKemerovo Institute of Food Science and Technology, bul'v. Stroitelei 47, Kemerovo, 650056 Russia
(Received March 7, 2014; accepted in revised form March 17, 2014)
Abstract: The PCEAS model was used to study the liquid-solid and liquid-vapor phase transitions at constant pressure in systems containing glycols and glycol ethers. This method is based on minimizing the excess Gibbs energy over the solvation parameter, which takes into account the processes of association of molecules in various phases. To compute the diagrams, the data on enthalpy and phase transition temperatures of pure components are required, while the information about the interactions in the binary system is not necessary. We present analytical expressions for the enthalpy of vaporization and enthalpy of melting of glycols and glycol ethers obtained with the theory of similarity using molecular weight, critical temperature, temperature of the triple point, and the number of carbon atoms in the molecule as the parameters. In the absence of information about the critical temperature, the enthalpy of vaporization may be calculated using the boiling point value. It is shown that the prediction of the enthalpy of melting and enthalpy of vaporization allows us to calculate of the phase diagram, as well as the azeotropic and eutectic parameters in water-glycol ether and glycol ether-alkane systems.
Keywords: modeling, glycol, glycol ether, PCEAS, enthalpy of vaporization, melting enthalpy, thermodynamic similarity, liquid-solid equilibrium, liquid-vapor equilibrium, eutectics, azeotrope
UDC 544.015; 517.958; 544.275.5; 544.275
DOI 10.12737/4142
INTRODUCTION
Glycols and glycol ethers are widely used as additives in food and cosmetics industries. Therefore, systems containing glycols and glycol ethers in mixtures with water, alkanes, and salts are actively studied [1-3]. These systems are the subjects of the current work because experimental and theoretical foundation for their application is lacking. There should be thermodynamic models that would be able to predict properties of pure components and phase equilibrium in systems containing glycols and glycol ethers with high accuracy.
For some of the members of the homology series, enthalpies of vaporization and melting are unknown or have not been measured with sufficient accuracy. In work [4] it has been noted that calculation of the parameters of eutectics and azeotrope mixtures used in technological processes of crystallization and rectification is complicated by the absence of accurate data on enthalpies of phase transitions of individual components.
Boiling temperature Tb and melting temperature Tm, pressure P, density p, enthalpy of vaporization AHvap, and enthalpy of melting AHm are determinative for computation of thermodynamic properties of pure substances and solutions. Prediction of the compound properties implies that the properties are determined basing on the compound structural formula. The theory
connecting the structure of a molecule with its macroscopic parameters is not complete yet; therefore, it is necessary to summarize empirical data on the properties of various compounds. Method of thermodynamics similarity, which is a part of general similarity theory, forms theoretical basis for such a summary [5].
Enthalpy of vapor formation at normal temperature of boilingis called the enthalpy of vaporization. Solid-liquid and liquid-vapor phase transitions proceeding at normal temperature of melting are characterized by enthalpy of melting and enthalpy of vaporization. Methods of group contribution [6, 7] and quantum chemistry [8] are used for predictions.
Enthalpy of phase transition of a pure component may be found from pressure-temperature dependence curve of the two-phase equilibrium, which is expressed as the Clausius-Clapeyron equation:
dP AH _ (i)
,/'/ T II ''
Provided that the evaporation heat does not depend on the temperature, vapor is considered an ideal gas, the volume of liquid is small compared to the gas volume, and the following equation may be written for the liquid-vapor coexistence curve:
In P - -A//vap IRT + const
If the data on vapor pressure are absent, enthalpy of vaporization may be calculated according to empirical equations, most of which are based on the data on critical parameters of the compounds [9, 10]. Critical properties of many gases and liquids are known or may be computed [11, 12]. In [13], an empirical formula is reported:
A# vap = A
Ta-T
Ta-Tu
where AHtr is the enthalpy of vaporization in the triple point; Tcr , critical temperature; Ttr, temperature of the triple point; zcr , specifies the number of similarity.
In work [13], authors arguethat a wide range of liquids, with the exception of quantum liquids, has a universal behavior on the evaporation curve. The group of substances with similar chemical structure is characterized by similarity parameters that vary upon transition to a different group. Similarity relations are associated with the concept of the critical state. According to the theory of similarity, a system of units is chosen in function of the nature of a subject under study and then the transition is made to non-dimensional, the so-called corrected, values. To perform such a change, that is, to transfer from one system to another, a linear transformation of values is to be done:
R N
where a(N)/R is similarity, a B(N) = (3(Ar)(—1)A -a more accurate
the determinative parameter of non-dimensional coefficient; y( Y) is the function introduced for modeling of the enthalpy of vaporization for the members of a homology series with even and odd numbers of carbon atoms; and P(N) and Y(N) are expressed in J/mol.
Table 1 reports the a(N) coefficient of the equation
(3) as a function of the number of carbon atoms in a molecule for glycols and glycol ethers, P(N) = y(N) = 0.
Table 1. Parameters of the equation for calculation of the enthalpy of vaporization of the members of the homology series of glycols and glycol ethers
Homology series a(7V), J/(mol K)
Glycols #(# + 1.1) 1.7 (A,r — 1)
Glycol ethers #(# + 1.1) 1.15(#-1)
Figure 1 presents the experimental data and the results of prediction of the enthalpy of vaporization for a homology series of glycol ethyl ethers.
where c is a constant value called the similarity
parameter and x, and x' are the respective values.
Laws revealed for one system are correct for a group of similar systems, which allows obtaining information on the properties of unknown systems.
The method of thermodynamic similarity may be briefly characterized as a method of isolating a number of parameters that adequately define the thermal physical properties of a substance. In the cases when the data on critical parameters are absent, it is reasonable to use models allowing to present the enthalpy as a function of the boiling temperature Tb and other reliably defined parameters.
Kistyakovsky proposed a rule to define the entropy of vapor formation:
AS^=^ = 8J5 + R]nTb,
(2)
which allows to calculate the entropy of vapor formation.
Equation (2) for various types of organic compounds rewritten by Vetere in the form of a function of Tb and molecular mass M, as well as other correlations, are reported in [7, 10].
MATERIALS AND METHODS
In the work, molar enthalpy of vaporization AHev at normal boiling point of a pure component for homology series of alkanes, alcohols, acids, glycols, and glycol ethers is expressed as a function of the number of carbon atoms N, relative molecular mass M, normal boiling temperature Tb, and the universal gas constant R:
Fig. 1. Enthalpy of evaporation of the homology series of glycol ethyl ethers: 1 - AHev, according to formula (3), and 2 - AHev - the experimental data [14].
The following relation is proposed for these members of different homology series:
AH”-B„(N) a„(N)MnT,
rev " n
Ml:-Bm{N) am {N)MmTr where £.(AO = (3 ,.(#)(-!)" + y,(A0, i = n.
(4)
b m
To relate the enthalpy of evaporation for different members of the same homology series with the number of carbon atoms N1 and N2, the following ratio may be used:
AHZ -B(N2) a(N2)M(N2)T*
(5)
where a(N2)/a(Ni) is the ratio of the determinative parameters of similarity.
b
Equation (1) may also be used to determine AHb, but the data on the dependence of melting temperature on the pressure are available for a limited range of pure components.
Enthalpy of melting of a pure component may be presented for a molecule with the number of carbon atoms of N as a function of molar mass M, melting temperature Tm, and the universal gas constant R:
A/7 =
a(N) MRTm
R
N
- A(N),
(6)
where A(N) = b(NX~lf + c(N) ■ Here a(N)/R is the determinative parameter of similarity, a non-dimensional coefficient; b(N) and c(N) are expressed in J/mol.
Table 2 provides the coefficients of the equation (3), a(N), b(N), and c(N), as functions of the number of carbon atoms in the molecule. This presentation is not the only one possible and depends on the accuracy of experimental data used to approximate the enthalpy.
AH*-A(N2) a(N2)M (N2)Tr
m
2 n7
AH™ -AIN,) aiN^AKN^
(8)
Relation (9) allows to calculate of the enthalpy of melting using the known enthalpy of vaporization value:
Mfm-b(N)(-l)N-c(N) a{N)Tm
A^ev-p^X-l)" -y(N) a(N)Tb
(9)
Table 3 presents the relation between the parameters of similarity a(N)/a(N) for such a comparison.
Table 3. Relation between the parameters of similarity for the members of the homology series of glycols and glycol ethers
Homology series a(N)/a(N)
Glycols (1 < N < 10) (N2 +l)(Ar-1) 15.59iV(iV + 1.1)
Glycol ethers (1 < N < 10) (W2 +1)(W-1) 23.04JV(JV + 1.1)
Table 2. Parameters of the equation used to calculate the enthalpy of melting of the members of homology series of glycols and glycol ethers
Homology series a(N), J/(mol K) b(N), J/mol c(N), J/mol
Glycols and glycol ethers (1</V<10) #2+l 26.5 0 -320N+10755
Figure 2 presents the experimental data and predicted values of the enthalpy of melting for the homology series of glycols.
A.Hm, J/mol 26000 -
24000
22000
20000
18000
16000
114000
12000
10000
0 2 4 6 8 10 N
Fig. 2. Enthalpy of melting of glycols: 1 - AHm according to the formula (3) and 2 - AHm, the experimental data [14].
The similarity of the enthalpy of melting taken up to a constant A(N) of the relevant members of homology series of glycols and glycol ethers is considered:
Mi:-An(N) a„(N)M„T,r (7)
AC-4W am(N)MnlT™ ’
where .4,.(TV) +c,.(7V), series number.
The following ration is true for similarity of the enthalpy of melting between the different members of the same homology series with carbon atom numbers N1 and N2:
The advantages of enthalpy of melting and enthalpy of evaporization calculations using the formula (4)-(5) and (7)-(8) is the decrease of the effect of systematic measurement errors in the course of molecular mass and melting temperature determination of the homology series members since they form a ratio in the formula
(4)-(5) and (7)-(8). Relation (9) does not include molecular mass, which excludes the influence of molecular mass measurement error. Since ratio Tm /Tb is in the right side of the relation, the contribution of systematic error of temperature measurement to the error of enthalpy calculation decreases.
RESULTS AND DISCUSSION
Empirical models for calculation of enthalpy of evaporation and enthalpy of melting based on the data of phase transition temperatures, molecular mass, and the number of carbon atoms in the molecule were proposed. Upon serial application of the similarity theory, critical temperature and triple point temperature are introduced into the system as reducing parameters, allowing for a decrease of the error of determination of the phase transition enthalpy
of a pure component to 1-2%.
The method of thermodynamic potentials was used as the mathematical method to study the phase equilibrium in real solutions [15, 16]. Phase chart eutectic and azeotropic systems (PCEAS) [17] is the model of state equation based on the minimization of excess Gibb’s energy over the solvation parameter X, characterizing the ratio of the number of A molecules to the number of B molecules in a compound.
The result of the work is the computation of liquid-solid and liquid-vapor phase equilibrium in binary systems containing glycol ethers using the PCEAS software. Enthalpy of melting and vaporization for glycol ethers necessary for calculation were found using the models proposed in the work. Phase equilibria taking into account association and solvation of molecules were studied under normal atmospheric pressure.
Figures 3 and 4 present the liquidus and solidus in the systems of 2-methoxyethanol-water and 2-methoxy-ethanol-cyclohexane.
T,K___________
260 1 —
250 240 230 220 210 200 190 180 170
0.1 0.2 0.3 0,4 0.5 0.6 0.7 0.8 0.9 1.0
Xj. mass fraction
Fig. 3. Liquid-solid phase diagram of the
2-methoxyethanol-water system.
Table 4 provides for the results of calculations of eutectics and azeotrope parameters in systems based on glycol ethers.
T, K 270 v 260 \ i 250 \
240 \ I
230 \ j
220 210 200 190 180 170
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
xj, mass fraction
Fig. 4. Liquid-solid phase diagram of the
2-methoxyethanol-cyclohexane system.
Table 4. Computed and experimental eutectics and azeotrope parameters in systems containing glycol ethers
Solution Xeut, mass., comp. 8 f m f ° . f qaz, mass., comp. t °C ‘•az? comp. qaz, mass., exp. [18] t °C ‘•az? exp. [18]
2 -methoxyethanol-water 0.959 -102.87 0.21 99.77 0.212 99.75
2 -methoxyethanol-cyclohexane 0.452 -111.66 0.09 77.56 0.080 79.80
2 -methoxyethanol-heptane 0.404 -102.40 0.18 91.88 0.230 92.50
2 -methoxyethanol-octane 0.809 -88.73 0.329 105.59 0.480 110.00
methyl carbitol-water 0.948 -91.09 non- azeotr. non- azeotr. non- azeotr. non- azeotr.
methyl carbitol -undecane 0.959 -85.65 0.340 173.59 0.400 (v.) 178.70
2 -ethoxyethanol-water 0.909 -80.57 0.297 99.74 0.300 99.30
2 -ethoxyethanol-heptane 0.291 -97.74 0.147 93.69 0.140 94.50
2 -ethoxyethanol-octane 0.801 -78.54 0.294 110.11 0.380 116.00
2 -ethoxyethanol-nonane 0.606 -79.91 0.462 120.85 0.50 128.00
Figure 5 presents the dependence of boiling temperature on mass fraction of the first component in the 2-methoxyethanol-water system.
r.K
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Xj. mass fraction
Fig. 5. Liquid-vapor phase diagram of the 2-methoxyethanol-water system.
Results of the calculations reported in Table 4 agree with the experimentally obtained parameters of eutectics and azeotropes [18].
CONCLUSIONS
Predicting thermodynamic parameters of glycols and glycol ethers based on the similarity theory allows to calculate phase diagrams of binary systems for which no experimental data on enthalpy of melting or enthalpy of evaporization of pure components is available. The method can be recommended when selecting components and composition in order to obtain the desired properties of the solution (temperature and enthalpy of melting and vaporizatio n), which will greatly reduce the time and expenses spent on experiments.
REFERENCES
1. Grenner, A., Kontogeorgis, G. M., Nicolas von Solms, and Michelsen, M. L., Application of PC-SAFT to glycol containing systems - PC-SAFT towards a predictive approach, Fluid Phase Equilibria, 2007, vol. 261, no. 1-2, pp. 248-257; http://dx.doi.org/10.10167j.fluid.2007.04.025.
2. Chen, L.-F., Soriano, A. N., and Li, M.-H., Vapour pressures and densities of the mixed-solvent desiccants
(glycols + water + salts), J. Chem. Thermodynamics, 2009, vol. 41, no. 6, pp. 724-730;
http://dx.doi.org/10.1016/jjct.2008.12.003.
3. Garrido, N. M., Folas, G. K., and Kontogeorgis, G. M., Modelling of phase equilibria of glycol ethers mixtures
using an association model, Fluid Phase Equilibria, 2008, vol. 273, no. 1-2, pp. 11-20;
http://dx.doi.org/10.1016/j.fluid.2008.08.006.
4. Kenisarin, M., and Mahkamov, K., Solar energy storage using phase change material, Renewable and Sustainable Energy Reviews, 2007, vol. 11, no. 9, pp. 1913 - 1965; http://dx.doi.org/10.1016/j.rser.2006.05.005.
5. Filippov, L. P., Prognozirovanie teplofizicheskikh svoistv zhidkostei i gazov (Predicting Thermal Physical Properties of Liquids and Gases), Moscow: Energoatomizdat, 1988.
6. Basarova, P., and Svoboda, V., Prediction of the enthalpy of vaporization by the group contribution method, Fluid Phase Equilibria, 1995, vol. 105, no. 1, pp. 27-47; http://dx.doi.org/10.1016/0378-3812(94)02599-V.
7. Marrero, J., and Gani, R., Group-contribution based estimation of pure component properties. Fluid Phase Equilibria, 2001, vol. 183-184, pp. 183-208; http://dx.doi.org/10.1016%2FS0378-3812%2801%2900431-9.
11. Gharagheizi, F., Determination of normal boiling vaporization enthalpy using a new molecular-based model, Fluid Phase Equilibria, 2012, vol. 317, pp. 43-51; http://dx.doi.org/10.1016%2Fj.fluid.2011.12.024.
12. Reid, R., Prausnitz, J. M., and Sherwood, T. K., Svoistva zhidkostei i gazov (The Properties of Gases and Liquids), Leningrad: Khimiya, 1982.
13. Kikic, I., and Vetere, A., Evaluation of several literature equations to predict the vaporization enthalpies at the
normal boiling point, Fluid Phase Equilibria, 2011, vol. 309, pp. 151-154;
http://dx.doi.org/10.1016%2Fj.fluid.2011.06.026.
14. Nesterov, I. A., Nesterova, T. N., Nazmutdinov, A. G., et al, Prediction of critical temperatures of the liquid-vapor equilibrium for organic compounds, Zh. Fiz. Khim. (J. Physical Chemistry), 2006, vol. 80, no 11, pp. 2032-2040.
15. Nesterova, T. N., and Nesterov, I. A. Kriticheskie temperatury i davleniya organicheskikh soedinenii. Analiz sostoyaniya bazy dannykh i razvitie metodov prognozirovaniya (Critical temperatyre and pressure values of organic compounds. Analysis of database conditions and development of methods of predicition), Samara: Izdatel'stvo Samarskogo nauchnogo tsentra RAN, 2009.
16. Meyra, A. G., Kuz, V. A., and Zarragoicoechea, G. J., Universal behavior of the enthalpy of vaporization: an
empirical equation, Fluid Phase Equilibria, 2004, vol. 218, pp. 205-207;
http://dx.doi.org/10.1016%2Fj.fluid.2003.12.011.
17. Dyment, O. N., Kazansky, K. S., and Miroshnikov, A. M., Glycols and other derivatives of propylene oxide and ethylene, Moskow: Khimiya, 1976.
18. Esina, Z. N., Korchuganova, M. R., and Murashkin, V. V., Mathematical modeling of the liquid-solid phase transition, Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel’naya tekhnika i informatika (Tomsk State University Bulletin. Management, computer engineering, and informatics), 2011, vol. 3, no. 16, pp. 13-23.
19. Esina, Z. N., Murashkin, V. V., and Korchuganova, M. R., Physical and mathematical model of the solid-gas equilibrium at constant pressure, Vestnik Kemerovskogo gosudarstvennogo universiteta (Kemerovo State University Bulletin), 2013,vol. 2-1, no. 54, pp. 72-80.
20. Esina, Z. N., Murashkin, V. V., and Korchuganova, M. R., «Phase Chart Eutectic and Azeotropic System (PCEAS)»: State registration certificate of an application software № 2012618394 (RU) - № 2012616324; filed 25.07.2012; registered in the catalog of the application software, Moscow, September 17, 2012.
21. Ogorodnikov, S. K., Lesteva, T. M., and Kogan, V. B., Azeotropnye smesi. Spravochnik (Azeotropic Mixtures. Reference Book), Kogan, V. B., Ed., Leningrad: Khimiya, 1971.
