Computation of Isobaric Vapor-Liquid EquilibriumData for Binary and TernaryMixtures of Methanol , Water , and Ethanoic Acid from T , p , x , and HEm Measurements

Vapor-liquid equilibrium (VLE) data for the strongly associated ternary system methanol + water + ethanoic acid and the three constituent binary systems have been determined by the total pressure-temperature-liquid-phase composition-molar excess enthalpy of mixing of the liquid phase (p, T , x, HE m) for the binary systems using a novel pump ebulliometer at 101.325 kPa. The vapor-phase compositions of these binary systems had been calculated from Tpx and HE m based on the Q function of molar excess Gibbs energy through an indirect method. Moreover, the experimental T , x data are used to estimate nonrandom twoliquid (NRTL), Wilson, Margules, and van Laar model parameters, and these parameters in turn are used to calculate vapor-phase compositions. The activity coefficients of the solution were correlated with NRTL, Wilson, Margules, and van Laar models through fitting by least-squares method. The VLE data of the ternary system were well predicted from these binary interaction parameters of NRTL, Wilson, Margules, and van Laar model parameters without any additional adjustment to build the thermodynamic model of VLE for the ternary system and obtain the vapor-phase compositions and the calculated bubble points.


Introduction
New strategies for the correlation and accurate prediction of the vapor-liquid equilibrium (VLE) data play a vital role in the distillation and separation process in chemical industry.The common technique for obtaining VLE data is by direct measurement on the system.That is to say, when the VLE is established, and phases are sampled and analyzed.Consequently, the experimental technique must be rather highly accurate to ensure meaningful results in the operation of equilibrium stills.Actually, when the vapor-phase components are sampled and analyzed, the whole compositions of components in solution and vapor have been changed.Accordingly, the behavior of the systems has been changed with the amount of compositions.Moreover, it has been long realized that the analysis of vapor-liquid composition for the infinite dilute solution is very difficult.In addition, for VLE measurements of mixtures containing a highly volatile compound, the accurate measurement of the vapor-phase composition can be difficult.This fact, coupled with the necessity for much analytical work, tends to enhance interest in exploring new methods for the determination of equilibrium data that do not involve sampling and analysis of the vaporphase components.
Several methods have been explored for the calculation of component behavior from gross solution behavior.Van Ness and coworkers have suggested the classification of these methods into two categories, direct and indirect methods [1].The direct methods involve calculation of vapor compositions by integration of the coexistence equation, a firstorder differential equation derived from the Gibbs-Duhem equation relating phase compositions at equilibrium.Hala and coworkers have given a detailed discussion of the basic direct method [2], and Van Ness and coworkers have discussed techniques for handling nonconstant temperature or pressure conditions as well as nonideal vapor phase behavior [1,3].Moreover, in the total pressure method one can calculate y from T, p, x measurements using an indirect method discussed by Mixon et al. [4].
Of all the methods, the indirect methods involve first the measurement, by some appropriate means, the temperature and total pressure of the system, the liquid-phase compositions, and subsequent calculation of vapor-phase compositions there from.These methods usually involve ascertaining which of selected solution equations to the Gibbs-Duhem equation lead to the best fit to the experimental data, and of the determination of the parametric values producing the best fit.For example, Barker has developed a procedure based on the assumption that the excess free energy can be represented as a polynomial function of composition [5].
There is a basic difference in the degree of rigor associated with the direct method and the indirect method of Barker.In the former, one makes no assumptions about the solution behavior for the nature of molecular interactions.Solution behavior is determined directly from the experimental data.The Barker method necessitates the assumption of a particular model and the estimation of its parameters, this deficiency in the method of Barker has been recognized by Tao, who has presented another indirect method in which the necessity for the a priori assumption of a particular functional form for the excess free energy has been removed.Tao's procedure involves calculation of the activity coefficient essentially by integration of an equation resembling the coexistence equation.His procedure, though indirect, retains the rigor usually associated with the direct method [6].The method of Tao appears specific to binary systems and does not seem to be easily generalized.However, an indirect method such as that of Barker is readily and easily generalized to ternary and higher-ordered systems, but this method retains the disadvantage of lacking rigor as compared with the direct method.
This paper presents the vapor-liquid compositions calculated based on the measurements of VLE data for temperature-total pressure-liquid-phase composition-molar excess enthalpy energy of mixing of the liquid phase (T, p, x, H E m ) at 101.325 kPa according to the Q function of molar excess Gibbs energy through an indirect method.We know that the reaction of methanol carboxylated with carbon monoxide is the most common and important technology for synthesis of ethanoic acid in the chemical industrial process.Modeling the thermodynamic properties and correlating and predicting the phase equilibria of a mixture involving associating components forming hydrogen bonding such as carboxylic acid remain a challenging problem, since such systems show extremely nonideal behavior and the formation of monomer, dimer, and even trimer in vapor and liquid phase.In addition, for the VLE measurements of vaporphase components containing a highly volatile compound such as carboxylic acid, the accurate measurement of the vapor-phase composition can be difficult.Many attempts have been made to describe the vapor-liquid equilibria of carboxylic acid containing mixtures using the concept of multiscale association [7].Arlt reported the isothermal VLE data of a new apparatus for phase equilibria in reaction mixtures containing water with ethanoic acid and propanoic acid at (333 to 363) K [8].Xu and Chuang have developed a new correlation for the prediction of the vapor-liquid equilibrium of methyl acetate-methanol-water-ethanoic acid mixtures [9].Sawistowski and Pilavakis explored the vaporliquid equilibrium behavior of the quaternary system methyl acetate-methanol-water-ethanoic acid modeled by using the Margules equation in combination with Marek's method for the association of ethanoic acid [10].Moreover, Guan et al. investigated that the isobaric vapor-liquid equilibria for water + ethanoic acid + n-pentyl acetate, isopropyl acetate, N-methyl pyrrolidone, or N-methyl acetamide were correlated and predicted by both the nonrandom two-liquid (NRTL) and universal quasichemical activity coefficient (UNIQUAC) models used in combination with the Hayden-O'Connell (HOC) method [11][12][13][14].In our recent work, we have concluded that the VLE data for the associating systems containing the carboxylic acid system can be correlated and predicted [15,16].Although the VLE data of the mixture containing the associating systems were previously reported by the different research groups [17][18][19][20][21], respectively, the VLE data for the associating binary system containing the carboxylic acid have been still extensively studied because of the extensive association effects occurring in them and the difficulty of properly calculating the activity coefficients [22][23][24][25].Nominally, the system is binary but in practice, it is multicomponent.Because the carboxylic acid monomer undergoes partial dimerization and even higher polymerization.This association is attributed to the formation of hydrogen bonds and occurs in both the vapor and liquid phase.Therefore, the challenge for the VLE data of the associating systems has evoked more and more researchers to focus on new strategies for exploring them.The VLE data for methanol + water + ethanoic acid ternary system and the constituent binary systems are indispensable in distillation separation process to the product of methanol carboxylation through the correlation and prediction by the new method, while some of the isobaric VLE data for these binary and ternary systems are correlated and predicted earlier [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25].To provide the new correlation and prediction for some necessary basic thermodynamic data on the separation process of methanol carboxylation, therefore, it is very indispensable for these systems studied on the VLE data of the constituent binary and ternary systems using the new method.This paper reports a novel correlation and prediction for the VLE data for these systems containing the associating component that has been developed.The VLE data for methanol + water + ethanoic acid system and constituent binary systems were measured by the total pressure-temperature-liquid-phase composition-molar excess enthalpy of mixing of the liquid phase (p, T, x, H E m ) for the binary systems using the novel pump-ebulliometer at 101.325 kPa.Owing to the association of ethanoic acid molecules, Marek's method in combination with Hayden-O'Connell (HOC) model was used to deal with the associating properties of the liquid and vapor phase and the nonideality of vapor phase, respectively [26][27][28].However, the nonideality of liquid phase was corrected by the calculation of its activity coefficient, which was obtained based on NRTL, Wilson, Margules, and van Laar models as the function of T, x through the nonlinear fit of the least-squares method.NRTL, Wilson, Margules, and van Laar models were applied to correlate the VLE data for the three constituent binary systems, and the model parameters together with the deviations of temperature and vaporphase molar fractions calculated from T, p, x, H E m according to the Q function of molar excess Gibbs energy by the indirect method were obtained by the least-squares method.The VLE data of the ternary system were well predicted from these binary interaction parameters of NRTL, Wilson, Margules, and van Laar models without any additional adjustment to build the thermodynamic model of VLE for the ternary system and obtain the vapor-phase compositions and the calculated bubble points.The excess Gibbs free energies for these binary systems as the function of liquidphase composition and activity coefficient were calculated through the activity coefficient correlation to NRTL model parameters with the experimental data.

Modeling Section
There is an added complexity when working with carboxylic acids because they associate in the vapor and liquid phases.This association can be represented by assuming that the organic acid exists as monomer and dimer according to the Marek's method [26,27].This fact, coupled with the necessity for much analytical work, tends to enhance interest in exploring the new strategies for correlation and prediction of the VLE data for the systems containing the associating carboxyl acid.
According to the Marek's chemical theory, there are monomer and dimer carboxylic acid molecules in both liquid and vapor phases.The equilibrium constant of vapor-phase dimerization, C V E , is calculated by the expression and the equilibrium constant of liquid-phase dimerization, C L E , is calculated by the expression ( In equations mentioned above, y * D and x * D are the mole fractions of dimers of ethanoic acid molecules in both vapor and liquid phases, respectively, and C V E can be defined as the function of temperature by the expression obtained from the literature [27] C where C V E was presented in kPa −1 and T in K.When the dimers of ethanoic acid molecules mainly exist, the binary system for the methanol or water (1) + ethanoic acid (2) is nominally binary; however, it is actually ternary for the methanol or water (1) + monomer ethanoic acid (M) + dimer ethanoic acid (D) system.In this system, the mole fractions of vapor-liquid equilibrium phases are y * 1 , y * M , y * D , x * 1 , x * M , and x * D , respectively.Therefore, the VLE relations of the nonassociating and associating components are individually calculated by the expressions In ( 4), the actual mole fractions can be denoted by the apparent mole fractions of easily determined components (y 1 , y 2 , x 1 , and x 2 multiplied by a modified coefficient.Likewise, the measured apparent vapor pressure saturated of ethanoic acid (p s 2 ) multiplied by a modified coefficient can also denote the actual vapor pressures saturated of monomer and dimer ethanoic acid (p s M , p s D ).Consequently, the apparent compositions and vapor pressures saturated substitute for the actual ones, and the VLE relation of the nonassociating component, methanol or water (1), is expressed by where and for ethanoic acid (2), which is the associating component, its relation is expressed by where In ( 5) to (11), α 1 and α 2 can be regarded as modified coefficients for the deviation from ideality in vapor phase on account of association.The fugacity of coefficients Φ 1 and Φ 2 is not negligible, and their values were obtained through the HOC model [28].From another point of view, θ s 2 can be viewed as a modified coefficient for the vapor pressure of the associating component.Moreover, β 1 and β 2 denote modified coefficients of the deviation from ideality in the liquid phase by reason of the existence of association, and γ 1 and γ 2 express the deviation from ideality in liquid  phase because of other factors.Herein, in liquid phase, the modified coefficients of deviation from ideality, β 1 and β 2 , can be incorporated to the activity coefficients, γ 1 and γ 2 , respectively.So ( 5) and ( 8) are obtained by the expressions

Experimental Section
3.1.Materials.Methanol (99.8 mass %) and ethanoic acid (99.8 mass %) were purchased from Sigma.The purities of the chemicals are provided by the manufacture's specifications.Ultra sound was used to dispel the solvent air in the materials, which were dried on molecular sieve (pore diameter 30 nm from Fluka) to remove all possible traces of moisture before use, but no other treatments were applied.The densities and refractive indices at 298.15 K and normal boiling points at 101.325 kPa of the pure component were compared with the literature values of Riddick et al. [29].
The results show that the measured values are approximately in agreement with those of the literature, as presented in Table 1.The measurement method of the composition dependence of densities and refractive indices has been previously reported [30].

Apparatus and Procedure.
A new type of magnetic pump ebulliometer described in detail by Qiu et al. [31] was used to measure the boiling points with different liquid phase compositions at the 101.325kPa.The experimental main apparatus for pump-ebulliometer is schematically shown in Figure 1.The recirculation still is entirely constructed from borosilicate glass.The main parts are a magnetic recirculating bump (F), a magnetic stirrer (G), a thermo-well (filled with a conducting oil) provided to enable good measurement of equilibrium temperature (H), a part-heating tube (resistance heating wire inserted) (K), and a tube window to observe in VLE status (I).The apparatus is an all-glass dynamic recirculation still with total volume of about 1.00 × 10 −4 m 3 .During the run, to avoid the over heating, the still was submerged in a constant temperature bath at about 3 • C below the equilibrium boiling point, which was obtained by the Nichrome wire in the part-heating tube (K) to partially heat the known mass of the material.The equilibrium pressure was measured using a Fischer digital manometer with a precision of ±0.01 kPa.The pressure measurement for the manometer had an uncertainty of ±0.07 kPa, as provided by the manufacturer's specifications.The total uncertainty of calibration and pressure measurement is estimated to be ±0.15 kPa because of the uncertainty of the calibration and measurement errors.Since the barometric pressure changed slightly, the experimental temperatures of the systems were automatically calibrated to that at 101.325 kPa with selfadjusted pressure system.The temperature was measured using a Heraeus QuaT100 quartz thermometer with a thermosensor, with an accuracy of ±0.01 K.The calibration of the thermometers was carried out at the accredited calibration laboratory (Quality and Technique Bureau, Anhui).The total uncertainty of calibration and temperature measurement is evaluated to be ±0.085K because of the uncertainty of the calibration, the probe's position, and the pressure fluctuations.The equilibrium temperature, T, was measured by means of the thermosensor inserted into the thermowell (filled with conducting oil) (H).In each experiment, a known mass of the material was introduced from the injector into the still from feed inlet (B) and heated to equilibrium boiling point of the system at a fixed pressure of 101.325 kPa by an automatic pressure regulation system.The liquid mixtures of required composition were prepared gravimetrically, with the use of a Sartorius electronic analytic balance (model ER-182A) with an accuracy of ±0.0001 g.The values of mole fraction were reproducible to ±0.0001 and have uncertainty of 0.01%.The ebulliometer was charged with the mixture of desired composition, and the boiler was then heated by nichrome wire wound around the boiler.After the liquid mixture started boiling, the bubbles along with the drops of liquid slowly spurted on the thermowell one by one through the tube window observed in VLE status (I).When the VLE state was attained by adjusting the pressure to 101.325 kPa, it remained constant for 20 min to ensure the stationary state, and then the boiling temperature was measured.
A flow microcalorimeter (model 2107, LKB produkter, Bromma, Sweden) was applied to measure the molar excess enthalpies, H E m , of the mixtures.The electrical-calibration apparatus and its operating procedure have been described elsewhere by Francesconi and Comelli [34].Two automatic burets (ABU, Radiometer, Copenhagen, Denmark) were used to pump liquids through the mixing cell of the calorimeter.The performance of the calorimeter was checked by measuring the H E m of the test mixture hexane + cyclohexane system reported by Gmehling [35].Agreement with literature is better than 0.5% over the central range of mole fraction of hexane.
The uncertainties in mole fraction and H E m are estimated to be 0.0005 and 0.5%, respectively.The liquid-phase mole fraction of component i, x i , could be calculated from the known mass of the material added to the still.The vaporphase mole fraction of component i, y i , was calculated from the experimental T, p, x, and H E m data based on the Q function (the function of molar excess Gibbs energy) by the indirect method.

Calculation of Vapor-Phase Mole Fraction y for the Binary Systems from T px and H E
m .The vapor-phase mole fraction y i of the components was calculated from and the fugacity of coefficients Φ i of the components was obtained by the expression where ϕ V i is the fugacity coefficient of component i in the vapor mixture, ϕ s i is the fugacity coefficient of component i at saturation, V L i is the molar volume of component i in the liquid phase, R is the universal gas constant, T is the experimental temperature, p is the total pressure (101.325kPa), and p s i is the vapor pressure of pure component i.The Antoine equations were applied to calculate the values of these vapor pressures.The Antoine constants A i , B i , C i , and the values of T c , p c , V c , Z c and ω were obtained from Shi et al. [36], as shown in Table 2.The Poynting correction factor was also included in the calculation of fugacity at the reference state.The liquid molar volumes were evaluated from the modified Rackett equation [37].
According to the Gibbs-Duhem equation, any extensive molar thermodynamic property of a given phase, such as the Gibbs and Helmholtz energies, the enthalpy, and molar volume, must satisfy the following differential relation: x i dM i = 0. ( 16) In ( 16) M is a general molar property, x i the molar fraction of component i in the phase under consideration, x j / = i the pertinent set of compositions, C the number of components, and M i the partial contribution of component i to M.
When considering VLE, the molar excess Gibbs energy, G E , can be evaluated from measurable (T, p, x) data using activity coefficient relations.Replacing M by G E /RT in the following equation yields the well-known relation where H E m and V E m are the molar excess enthalpy and volume of mixing of the liquid phase.According to the thermodynamic principles, the activity coefficients γ i of the components were calculated from the expression as follows: where Q is the function of molar excess Gibbs energy and combining ( 17) and ( 18) yields Application to a binary system gives Simultaneous solution of ( 18) and (20) yields Equation ( 18) substituted into (21) reduced it to For a binary system comprised of species 1 and 2 at VLE state, at constant pressure, (dp/dx 1 ) = 0, substitution into (22) reduced it to The right-hand side, (H E m /RT 2 )(dT/dx i ), of ( 23) cannot be neglected.Proper use of ( 23) requires the availability Table 2: Antoine coefficients of the compounds and published parameters [36] used to calculate fugacity coefficients: critical temperature T c , critical pressure p c , critical volume V c , critical compression Z c , and acentric factor ω of pure compounds. of heats of mixing as a function of composition and temperature.The activity coefficients γ i of the components as functions of the excess Gibbs energy are as follows: From ( 14), this equation is rearranged to obtain Because y i = 1, for the binary system, the equation may be summed to give In (26), solved for y i by difference method.Suppose that [0, 1] is subdivided into n subintervals [x k , x k+1 ] of equal step size h = 1/n by using x k = kh for k = 0, 1, . . ., n.In k difference point, we obtain Meanwhile, ( 23) may be shown as follows: And ( 27) is linearized to obtain: The number of n linear equation from ( 28) is solved for ΔQ| k by chasing method where t is relaxation factor, finally, y ik (k = 1 ∼ n) is obtained by difference method.The block diagram for calculation procedure was detailedly shown in Figure 2.

Calculation of Vapor-Phase Mole Fraction y for the Binary Systems from Model.
There are many methods concerning the correlation and prediction of VLE data.The model-free approach data treatment of vapor-liquid equilibrium is also one of the best strategies for the correlation and prediction of VLE data.Wisniak's group has developed that the novel model-free computation techniques and limiting conditions have been applied to VLE data for azeotropic systems [42].Moreover, Segura and coworkers reported that a modelfree approach dealt with VLE data in application of ternary systems [38,43].For the three binary systems of this study, the activity coefficients were correlated with the NRTL [39], Wilson [40], Margules [41], and van Laar [44] equations, respectively.The interaction parameters optimized were achieved by the objective function (OF) minimized using the least-squares fitting as follows: where N is the number of experimental points x i,cal and x i,exp are the liquid-phase mole fraction of component i calculated and experimental values from the ( 12) or ( 13) and from measured data, respectively.Because carboxylic acids are always present in an associated form, like a dimer or trimer, in both the vapor and liquid phases even at low pressures, a significant deviation in fugacity coefficient may exist using the ideal gas assumption.To illustrate the deviation from ideal behavior, Marek's chemical theory and HOC model were applied to deal with the associating component and modify the deviation from idealities of vapor phase [26][27][28], respectively.The Poynting correction factor was also included in the calculation of fugacity at the reference state.The liquid molar volumes were evaluated from the modified Rackett equation [37].However, under isobaric conditions, the most volatile component cannot exist in the liquid state, only as superheated vapor.Hence, there is no way to calculate or measure this property for the molar volumes of the pure liquids.Therefore, the correct procedure for isobaric measurements is to calculate the overall values of A i j and A i j , as adjustable parameters,    [32,33]; •, liquidphase experimental temperature; •, liquid-phase mole fraction x 1 from the literature [32,33]; -, NRTL correlation temperature.
and not the values of the interaction excess energy.Hence, for Wilson model, A i j and A i j were reported as adjustable parameters.
The activity coefficients computed on the basis of NRTL model were used to evaluated dimensionless excess Gibbs function at 101.325 kPa for three binary systems over the overall range of composition.Liquid-phase mole fraction x 1 , experimental boiling point temperature T exp , calculated bubble point temperature T cal , vapor-phase mole fraction y 1 from the model, activity coefficients γ 1 and γ 2 using NRTL equation correlation, fugacity coefficients ϕ V 1 and ϕ V 2 , and the molar excess enthalpies of mixing of the liquid phase H E m are included in Table 3.The T-x 1 -y 1 diagrams for the methanol (1) + water (2), methanol (1) + ethanoic acid (2), and water (1) + ethanoic acid (2) three binary systems at 101.325 kPa are shown in Figures 3, 4, and 5, respectively.The plot of excess Gibbs energy function G E /RT versus liquid-phase mole fraction x 1 for the three binary systems is shown in Figure 6.All the mixtures exhibit deviations from ideality with a range that may be attributed to interactions leading to the formation of various associated aggregates.Observed nonideal behavior is indicative of the magnitude of experimental activity coefficients γ i , as well as of the variation of excess Gibbs energy function, G E /RT, with composition, as depicted in Figure 6.The obtained absolute maximum values of G E /RT for the methanol (1) + water (2), methanol (1) + ethanoic acid (2), and water (1) + ethanoic acid (2) three binary systems are 0.1703, 0.0180, and 0.0892, respectively.The values of excess Gibbs energy Table 3: VLE data for the methanol (1) + water (2), methanol (1) + ethanoic acid (2), and water (1) + ethanoic acid (2).Three binary systems at 101.325 kPa: liquid-phase mole fraction x 1 , experimental boiling point temperature T exp , calculated bubble point temperature T cal , vapor-phase mole fraction y 1 , activity coefficients γ 1 and γ 2 using NRTL equation correlation, fugacity coefficients ϕ V 1 and ϕ V 2 , and molar excess enthalpies of mixing of the liquid phase H E m .function G E /RT are positive for methanol (1) + water (2) and water (1) + ethanoic acid (2) binary systems.However, for methanol (1) + ethanoic acid (2) system, the values of those are negative in the overall range of mole fraction.G E /RT values follow the order methanol (1) + water (2) > water (1) + ethanoic acid (2) > methanol (1) + ethanoic acid (2).The absolute maximum value of G E /RT is approximate at an equimolar fraction in three binary systems.Figures 3-5 show that the comparison of the predicted vapor-phase and experimental liquid-phase compositions with those of the literature [10,11,32,33].Comparing with the values of vaporphase and liquid-phase component from the literatures, the values of those from the paper are in good agreement with the literature, as shown in Figures 3-5.The results have demonstrated that the methods for T px and H E m are appropriate for representing the experimental data of the three binary systems.The optimum model interaction parameter of liquid activity coefficient and the absolute average deviations for the different models, and from T px and H E m are listed in Table 4. Herein, we obtained the results by the four different types of correlations for the prediction of activity coefficients in these systems, which reveal that the deviations of NRTL, Wilson, Margules, and van Laar equations are reasonably small in Table 4.For comparison, the mean deviations obtained by Gmehling and Onken [32,33,45,46] are also shown in the Table 4.It can seem that two sets of deviation values are comparable.Since the superiority of one method over the others is not always obvious, practice must rely on experience and analogy.The comprehensive comparisons of four of the methods (NRTL, Wilson, van Laar, and Margules) were made in Table 4. From the data analysis, the temperature deviations between the experimental and calculated values of four different types of model are very similar in the three binary systems, and the vapor-phase mole fraction deviations between calculated values from T px and H E m , and from the model are very similar.Therefore, the activity coefficient models are appropriate for representing the experimental data of the three binary systems.In Table 4, the absolute average deviations dT of the difference between boiling point temperature from experiment and bubbling point temperature from calculation by NRLT model parameters for the three binary systems are 0.28 • C, 0.51 • C, and 0.33 • C, respectively.And the absolute average deviations dy of the difference between vapor-phase mole fraction from T px and H E m calculation and from NRTL model calculation are 0.0072, 0.0054, and 0.0096, respectively.

Correlation and Prediction of VLE of Ternary System.
The binary interaction parameters of the NRTL, Wilson, Margules, and van Laar model given in Table 4 were used to correlate and predict the VLE data of the ternary system.VLE data for methanol (1) + water (2) + ethanoic acid (3) at 101.325 kPa included liquid-phase mole fraction x 1 , x 2 , and x 3 , experimental boiling point temperature T exp , calculated bubble point temperature T cal , calculated vaporphase mole fraction y 1 , y 2 , and y 3 , activity coefficients γ 1 , γ 2 , and γ 3 , the average deviation in the bubble temperatures of the ternary system using NRTL equation correlation listed in Table 5.The absolute average and maximum deviation

Conclusions
VLE data for the ternary system methanol + water + ethanoic acid and three constituent binary systems at 101.325 kPa: methanol + water, methanol + ethanoic acid, and water + ethanoic acid were determined at different liquid-phase compositions using a novel pump ebulliometer.The equilibrium composition of the vapor phase was calculated from T, p, x, and H E m based on the Q function of excess Gibbs free energy by the indirect method.The experimental data were correlated using the NRTL, Wilson, Margules, and van Laar equations.It was shown that the deviations of NRTL, Wilson, Margules, and van Laar equations are reasonably small.The VLE data of ternary system were predicted by NRTL, Wilson, Margules, and van Laar equation.The calculated bubble points accorded well with the experimental data.The results show that the calculated bubble point is fitted by the models, which satisfy the need for the design and operation of separation process in chemical industry.Moreover, the method will provide theoretical guidance for the research of VLE data of strongly associating system of vapor and liquid phase in nonideal behavior.
, B i , C i , T c , p c , V c z, and p, T, x measured T∼ x, H m E ∼ x by the cubic spline fit Calc.

Figure 2 :
Figure 2: Block diagram for the calculation y by T, p, x, and H E m .

Table 1 :
Physical properties of the pure compounds: densities ρ, refractive indexes n D at 298.15 K and normal boiling points T b .

Table 4 :
Correlation parameters for activity coefficients, average deviation for studied systems and vapor-phase composition mean absolute deviation from the literatures.Margules and van Laar interaction parameters (dimensionless).dT = |T exp − T cal |/N ; N: number of data points; T cal : calculated bubble point from model, K; T exp : experimental boiling point temperature, K. dy = |y cal − y mod |/N ; N: number of data points; y cal : calculated vapor-phase mole fraction from Tpx; y mol : calculated vapor-phase mole fraction from model.c The values in parentheses from the literatures [38-41].