Dataof (Aceticacid+water) at (412.6,443.2,483.2) K
Dataof (Aceticacid+water) at (412.6,443.2,483.2) K
2078–2082
DOI:10.1021/acs.jced.5b01104
Abstract
Experimental – data were measured for the (acetic acid + water) system at the
temperatures of (412.6, 443.2 and 483.2) K between the pressures of 196 kPa and
1902 kPa over the entire range of concentrations. Experiments were carried out in a
static-analytical equipment with quantitative analysis by GC. The experimental data
were correlated with the Perturbed-Chain Statistical Associating Fluid Theory
considering two association sites for acetic acid and four association sites for water.
Keywords
Published by American Chemical Society. This is the Author Accepted Manuscript issued with:
Creative Commons Attribution Non-Commercial No Derivatives License (CC:BY:NC:ND 4.0).
The final published version (version of record) is available online at DOI:10.1021/acs.jced.5b01104
Please refer to any applicable publisher terms of use.
* Corresponding author. Tel.: (+44) 0121 414 5351. E-mail address: [email protected].
1
1 Introduction
Acetic acid and its derivatives are used in several industries making it one of the major
worldwide commodity chemicals. The global demand for acetic acid in 2013 was
around 10.5 Mt and it is projected to reach 15.5 Mt by 2020 with a market revenue of
12,190 million USD [1,2]. More than 65% of the acetic acid produced in the world ends
in the form of polymers, mainly as derivatives of vinyl acetate and cellulose acetate
[3]. Other usages include the production of pharmaceuticals (e.g. acetyl salicylic acid)
[4], as a reaction solvent of purified terephtalic acid (PTA) for further production of
polyethylene terephthalate (PET) [3], and in the food industry, as an acidulant,
preservative and flavouring agent [5]. Its most common household usage is in the form
of vinegar and as a chemical descaler.
Experimental vapour–liquid equilibria (VLE) data for the system (acetic acid + water)
are the most readily available of the carboxylic acids; more than 60 different articles
were found in the open literature. The oldest work seems to date back to 1921 [7] while
the most recent one to 2012 [8]. Most of the available data are isobaric, the majority
are made at sub- and atmospheric conditions. Wichterle, et al. [9,10] and Gmehling
and Onken [11] have compiled most of these work. None of the available literature, to
our knowledge, have reported azeotropic behaviour of the mixture.
Most of the isothermal data available are at temperatures below 373 K, a temperature
below the normal boiling point of acetic acid (391.05 K), mainly because of the
corrosive nature of the compounds; the common apparatus for these kind of
measurements is made of glass. Special alloys are required to handle the increasing
corrosive conditions encountered at higher temperatures. An exemption to the low
temperature data is that of Freeman and Wilson [12,13] who presented – – data at
2
= (372.77, 412.57, 462.06 and 502.86) K and pressures ranging from 56 kPa to
2778 kPa.
The Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) [14,15] EoS was
selected in this work to model the experimental data since it has been developed to
handle association interactions explicitly. In the PC-SAFT framework, the association
scheme for the molecules involved has to be defined, for which the classification of
Huang and Radosz [16] is commonly used. Based on previous studies [15,17-19] it
was decided in this work to model acetic acid with the 2B association scheme, while
water with the 4C association scheme. Pure component parameters were refitted in
each case. Additionally, the EoS results were compared with those obtained from the
UNIQUAC-HOC [20,21] activity coefficient model, where appropriate.
2 Experimental
2.1 Materials
Table 1 summarizes the chemical compounds used in the experiments as well as their
corresponding purities. Gas-Chromatography (GC) analysis of the purchased acetic
acid revealed two main peaks, one corresponding to the carboxylic acid while the other
3
to water in accordance with its hydrophilic characteristics. A third small peak
accounting for 0.028% of the mass sample was an unidentified impurity which was
considered to be part of the water content. Chemicals were used without further
purification and were only subjected to a degassing process as described in Section
2.3.
4
Combined standard uncertainties in pressure, ( ), after considering calibration,
repeatability and pressure drop during sampling were estimated as ( ) = 1 kPa.
The GC was equipped with a TCD detector and a packed column 3 ft x 1/8” SS
Porapak N 80/100 mesh. GC operating conditions were: inlet temperature: 503 K;
oven temperature: 423 K, column flow: 40 mL/min (1 min), 60 mL/min (4 min) at 10
mL/min; detector temperature: 473 K; carrier gas: Helium CP grade (BOC, certified
purity ≥ 0.99999 mole fraction). 0.2 μL volume samples were injected with an
autosampler. The GC calibration procedure of Raal and Mühlbauer [24] for liquid
sample injections was adopted in this work. A linear relationship of area ratio vs. mass
fraction ratio was obtained from the low to the high acetic acid concentration regions.
Uncertainties in composition are presented for each experimental point in Section 4.
2.3 Procedure
The experimental procedure has been described previously [22,25] and is briefly
summarized as follows: the equilibrium cell was cleaned thoroughly before each
experimental run and subjected to a pressure test with compressed nitrogen. Acetic
acid and water were degassed for 1 h in an ultrasonic bath. A mixture of the
compounds was loaded into the equilibrium cell and the system was vacuumed down
between = 1.5 kPa and = 2.0 kPa at room temperature. The desired equilibrium
temperature was reached by increasing or decreasing the air bath temperature. The
system was allowed to reach equilibrium state under constant stirring, which was
assumed when temperature and pressure did not vary within ±0.05 K and ±0.5 kPa,
respectively, for at least 5 min. A minimum of five samples (20 μL each) of the liquid
phase were withdrawn and collected in 250 μL vial inserts. Pressure was increased by
adding water into the cell by means of the liquid-pump in order to stablish a new
experimental point.
5
3 Modelling
= + + (1)
=1+ + + (2)
The five pure component parameters needed to characterize a pure component are:
the number of segments per chain ( ), the segment diameter ( ), the depth of the pair
potential ( ), the association energy ( ) and the association volume ( ). For
mixtures, conventional combining rules for and are used:
( + )
= (3)
2
= (1 − ) (4)
where is the binary interaction parameter. Wolbach and Sandler [15] combining
rules were used for the cross-association interactions:
+
= (5)
2
6
= (6)
( + )/2
Pure component parameters were obtained by fitting vapour pressure ( ) and liquid
density ( ) data simultaneously with the procedure described in a previous
publication [22], with the following objective function:
, − , , − ,
= (7)
, ,
The optimum was obtained by regressing the experimental bubble point pressures
( ) with Eq. (8) as the objective function:
−
= (8)
4.1 Experimental
The experimental – data obtained at = (412.6, 443.2 and 483.2) K as well as the
uncertainties in composition are shown in Table 2. Relative errors in the measured
water vapour pressures against reference data from NIST [23] were: 0.07, 0.12, and
0.40% at = (412.6, 443.2 and 483.2) K, respectively.
7
their studies of (acetic acid + water) in a Type SS-316 still that a presence of 2%
content gave no interference to the phase behaviour when testing a 85 wt.% acetic
acid solution at = 503 K and = 2.17 MPa.
Table 2. Experimental equilibrium data for the {acetic acid (1) + water (2)} system at
pressure , liquid mole fraction and temperatures = (413.2, 443.2 and 483.2) K.a
= 412.6 K = 443.2 K = 483.2 K
/kPa ( ) /kPa ( ) /kPa ( )
196 0.976 0.005 426 0.983 0.005 1012 0.977 0.005
210 0.951 0.005 437 0.975 0.005 1052 0.955 0.005
230 0.889 0.005 461 0.955 0.005 1150 0.895 0.005
266 0.721 0.005 504 0.903 0.005 1272 0.822 0.005
290 0.579 0.005 551 0.823 0.008 1401 0.689 0.005
311 0.433 0.005 605 0.697 0.005 1515 0.580 0.005
326 0.347 0.005 654 0.568 0.001 1652 0.444 0.005
333 0.282 0.005 704 0.409 0.003 1690 0.407 0.005
341 0.186 0.005 731 0.315 0.001 1750 0.335 0.005
348 0.102 0.010 760 0.195 0.005 1796 0.275 0.005
349 0.085 0.005 768 0.178 0.005 1821 0.234 0.005
353 0.014 0.005 774 0.148 0.005 1855 0.178 0.005
356 0.000 780 0.133 0.005 1870 0.153 0.005
786 0.105 0.005 1881 0.100 0.005
789 0.055 0.005 1890 0.078 0.005
792 0.000 1896 0.048 0.005
1900 0.034 0.005
1902 0.000
a Combined standard uncertainties, , are ( ) = 0.1 K, ( ) = 1 kPa. ( ) are displayed
in each temperature column.
4.2 Modelling
The pure component parameters for both compounds are presented in Table 3, with
acetic acid and water modelled as a 2B and a 4C association molecule, respectively.
The effect of different parameters for water on the modelling has been discussed
before [27].
8
d Δ = ∑ , for either or .
The adjusted were -0.0099, -0.0012 and -0.0016 for the (412.6, 443.2 and 483.2)
K isotherms, respectively. These values are one order of magnitude lower than the
ones reported in previous publications [18,19].
Figure 1 shows a comparison of the new data against those of the literature values of
Freeman and Wilson [12] at = 412.6 K. Both data sources for the liquid compositions
agree relatively well. It is necessary, however, to point out some aspects of the
literature values. They exhibit a maximum pressure of = 359 kPa in the dilute acetic
acid region after which the pressure decreases sharply, reaching a water vapour
pressure of 350.2 kPa. There is no obvious reason for this maximum to happen. It is
likely that repeatability was not easy to achieve because of thermal and composition
gradients not eliminated due to the lack of stirring. It is reasonable to assume that
deviations in pressure must be larger than the reported ones since the stated vapour
pressure of water at 412.6 K gives a relative error of 1.6% (NIST data as reference) in
contrast to that of 0.07% of this work. One can conclude that the experimental values
of the liquid compositions obtained in the present work are more reliable than those
reported by the previous authors, although the accuracy reported in the literature is
actually lower ( ( ) = 0.001 vs. an average of ( ) = 0.005 in this work).
Figure 1 also presents the modelling results with PC-SAFT, and as a comparison,
modelling with the UNIQUAC [20] activity coefficient coupled with the Hayden-
O’Connell term [21] (UNIQUAC-HOC), which considers the non-idealities of the
vapour phase of organic acids, is presented as well. The following default pure
compound and binary interaction parameters were used in the UNIQUAC-HOC model:
association parameter ( ) = 4.5, 1.7 and 2.5 for the acetic acid, water and cross-
interactions, respectively; = 0.7446, = 0.0042, =-615.264 and =
196.899. As can be seen, in general the liquid phase concentrations are better
correlated by PC-SAFT, even in the dilute acetic acid region (Figure 2) with the
exception of the concentrations of < 0.014.
9
The isotherms at 443.2 K and 483.2 K are presented in Figures 3 and 4, respectively.
It can be observed that PC-SAFT can correlate satisfactorily the bubble pressures at
both temperatures.
As an additional plot, the isotherm reported by Freeman and Wilson [12] at = 502.85
K is presented in Figure 5. PC-SAFT is able to model accurately the compositions of
both phases with = -0.012.
10
360
340
320
300
280
p/kPa
260
240
PC-SAFT, kij = -0.0099
220 UNIQUAC-HOC
200
180
0.0 0.2 0.4 0.6 0.8 1.0
x1, y1
Figure 1. Vapour–liquid equilibrium diagram for the {acetic acid (1) + water (2)} system
at = 412.6 K. Experimental data: (x) Freeman and Wilson [12]; (●) this work.
11
360
350
340
p/kPa
330
PC-SAFT, k ij = -0.0099
UNIQUAC-HOC
320
310
0.00 0.05 0.10 0.15 0.20 0.25 0.30
x1, y1
Figure 2. Vapour–liquid equilibrium diagram for the {acetic acid (1) + water (2)}
system in the =(0 to 0.3) region at = 412.6 K. Experimental data: (x)
Freeman and Wilson [12]; (●) this work.
12
800
700
p/kPa
600
400
0.0 0.2 0.4 0.6 0.8 1.0
x1, y1
Figure 3. Vapour–liquid equilibrium diagram for the {acetic acid (1) + water (2)}
system at = 443.2 K. Symbols: experimental data from this work.
13
2000
1800
1600
p/kPa
1400
1000
x1, y1
Figure 4. Vapour–liquid equilibrium diagram for the {acetic acid (1) + water (2)}
system at = 483.2 K. Symbols: experimental data from this work.
14
2800
2600
2400
2200
p/kPa
2000
1800
PC-SAFT, kij = -0.012
1600
1400
x1, y1
Figure 5. Vapour–liquid equilibrium diagram for the {acetic acid (1) + water (2)} system
at = 502.85 K. Symbols: experimental data from Freeman and Wilson [12].
15
5 Conclusions
Experimental pressures and liquid phase compositions were measured for the (acetic
acid + water) system at the temperatures of (412.6, 443.2 and 483.2) K. PC-SAFT,
with water modelled as a 2B and acetic acid as a 4C association molecule, was able
to correlate satisfactorily the experimental data including high-pressure data. New
pure component parameters for acetic acid and water are presented for the
corresponding association schemes.
Acknowledgments
L.A. Román-Ramírez gratefully acknowledges the National Council for Science and
Technology (CONACyT) of Mexico for their financial support through a PhD fellowship.
The authors would like to thank the technical assistance of Mr Robert W. Sharpe in
setting up the experimental equipment.
16
References
[1] Grand View Research Acetic Acid Market Analysis and Segment Forecasts to 2020.
http://www.grandviewresearch.com/industry-analysis/acetic-acid-market (7 August 2014).
[3] Cheung, H.; Tanke, R. S.; Torrence, G. P., Acetic Acid. In Ullmann's Encyclopedia of
Industrial Chemistry [online]; Wiley-VCH: 2011; pp. 209-237.
http://dx.doi.org/10.1002/14356007.a01_045.pub2 (30 July 2014).
[4] Burdick, D. L.; Leffler, W. L., Petrochemicals in Nontechnical Language. In [online]; 3rd ed.;
PennWell: 2001; pp. 256-261.
http://app.knovel.com/hotlink/toc/id:kpPNTLE002/petrochemicals-in-
nontechnical/petrochemicals-in-nontechnical (30 July 2014).
[6] Chiusoli, G. P.; Maitlis, P. M., Metal-Catalysis in Industrial Organic Processes. In [online];
Royal Society of Chemistry: 2006; pp. 119-120.
http://app.knovel.com/hotlink/toc/id:kpMCIOP001/metal-catalysis-in-industrial/metal-
catalysis-in-industrial (July 2014).
[7] Pascal, P.; Dupuy, E.; Garnier, M., Study of Binary and Ternary Mixtures Found in the
Manufacture of Synthetic Acetic Acid. Bull. Soc. Chim. Fr. 1921, 29, 9-21.
[8] Xin, H.; Wang, X.; Li, J. Q., Isobaric Vapor–Liquid Equilibria for the System Containing
Acetic Acid at 101.3 kPa. Advanced Materials Research 2012, 560-561, 79-85.
[9] Wichterle, I.; Linek, J.; Hála, E., Vapor–Liquid Equilibrium Data Bibliography. Elsevier:
Amsterdam, 1973.
[10] Wichterle, I.; Linek, J.; Hála, E., Vapor–Liquid Equilibrium Data Bibliography: Supplement.
Elsevier: Amsterdam, 1976.
[11] Gmehling, J.; Onken, U., Vapor–Liquid Equilibrium Data Collection: Aqueous–Organic
Systems. Chemistry Data Series. DECHEMA: Frankfurt, Germany, 1977; Vol. 1. Part 1.
17
[13] Freeman, J. R.; Wilson, G. M., High Temperature PVT Properties of Acetic Acid/Water
Mixtures. In AIChE Symposium Series 244, Benson, M. S.; Zudkevitch, D., Eds. American
Institute of Chemical Engineers: New York, 1985; Vol. 81, pp 1-13.
[14] Gross, J.; Sadowski, G., Perturbed-Chain SAFT: An Equation of State Based on a
Perturbation Theory for Chain Molecules. Ind. Eng. Chem. Res. 2001, 40, 1244-1260.
[15] Gross, J.; Sadowski, G., Application of the Perturbed-Chain SAFT Equation of State to
Associating Systems. Ind. Eng. Chem. Res. 2002, 41, 5510-5515.
http://dx.doi.org/10.1021/ie010954d.
[16] Huang, S. H.; Radosz, M., Equation of State for Small, Large, Polydisperse, and
Associating Molecules. Ind. Eng. Chem. Res. 1990, 29, 2284-2294.
[18] Chen, Y.; Afef, A.; Fabrice, M.; Roland, S.; Jeday, M. R., Thermodynamic Modeling of
Mixtures Containing Carboxylic Acids Using the PC-SAFT Equation of State. Ind. Eng. Chem.
Res. 2012, 51, 13846-13852. http://dx.doi.org/10.1021/ie301930q.
[19] Janecek, J.; Paricaud, P., Influence of Cyclic Dimer Formation on the Phase Behavior of
Carboxylic Acids. II. Cross-Associating Systems. J. Phys. Chem. B 2013, 117, 9430-9438.
http://dx.doi.org/10.1021/jp4012125.
[20] Abrams, D. S.; Prausnitz, J. M., Statistical Thermodynamics of Liquid Mixtures: A New
Expression for Excess Gibbs Energy of Partly or Completely Miscible Systems. AlChE J. 1975,
21, 116-128. http://dx.doi.org/10.1002/aic.690210115.
[21] Hayden, J. G.; O'Connell, J. P., A Generalized Method for Predicting Second Virial
Coefficients. Industrial and Engineering Chemistry Process Design and Development 1975,
14, 209-216. http://dx.doi.org/10.1021/i260055a003.
[22] Román-Ramírez, L. A. Vapour – Liquid Equilibria of Acetic Acid + Water and Propanoic
Acid + Water: Experimental Measurement and Thermodynamic Modelling. PhD, The
University of Birmingham, Birmingham, United Kingdom, 2015.
[23] NIST. Chemistry Webbook. NIST Standard Reference Database 69. [online];
http://webbook.nist.gov/chemistry/ (October 2012).
[24] Raal, J. D.; Mühlbauer, A. L., Phase Equilibria: Measurement and Computation. Taylor &
Francis: Washington D.C., USA, 1998.
[25] Román-Ramírez, L. A.; García-Sánchez, F.; Santos, R. C. D.; Leeke, G. A., Vapour–
Liquid Equilibrium of Propanoic Acid + Water at 423.2, 453.2 and 483.2 K from 1.87 to 19.38
bar. Experimental and Modelling with PR, CPA, PC-SAFT and PCP-SAFT. Fluid Phase
Equilib. 2015, 388, 151-159. http://dx.doi.org/10.1016/j.fluid.2015.01.004.
18
[26] Othmer, D. F.; Silvis, S. J.; Spiel, A., Composition of Vapors from Boiling Binary Solutions.
Pressure Equilibrium Still for Studying Water–Acetic Acid System. Ind. Eng. Chem. 1952, 44,
1864-1872. http://dx.doi.org/10.1021/ie50512a041.
[27] Liang, X.; Tsivintzelis, I.; Kontogeorgis, G. M., Modeling Water Containing Systems with
the Simplified PC-SAFT and CPA Equations of State. Ind. Eng. Chem. Res. 2014, 53, 14493-
14507. http://dx.doi.org/10.1021/ie501993y.
[28] DIPPR, 801 Database. Data Compilation of Pure Compound Properties. AIChE: 2012.
19
2000
T = 483.2 K
1500
p/kPa
1000
T = 443.2 K
500
T = 412.6 K
x, y acetic acid
for Table of Contents use only.
– data of (acetic acid + water) at = (412.6, 443.2, 483.2) K.
Luis A. Román-Ramírez
Gary A. Leeke
20