Isobaric Vapor-Liquid Equilibrium data for Binary Mixtures of n-butylamine and Triethylamine with Cumene at 97.3 kpa

Similar documents
Isobaric Vapor-Liquid Equilibria of Mesitylene + 1- Heptanol and Mesitylene +1-Octanol at 97.3 kpa

A simple expression for radial distribution functions of pure fluids and mixtures

Mean Activity Coefficients of Peroxodisulfates in Saturated Solutions of the Conversion System 2NH 4. H 2 O at 20 C and 30 C

A Cubic Equation of State for Reservoir Fluids

An Improved Model for Calculating Heats of Dilution and Equilibrium Constants for High Temperature Aqueous Electrolyte Solutions

Characterizing Pure and Undefined Petroleum Components

Modeling real gas equations of state in high density combustion

A multi-fluid nonrandom lattice fluid model: General derivation and application to pure fluids

NMR spin-lattice relaxation time and activation energy in some molecular systems

General Equilibrium. What happens to cause a reaction to come to equilibrium?

( ) ( ) Volumetric Properties of Pure Fluids, part 4. The generic cubic equation of state:

Chapter 8 Thermodynamic Relations

Isobaric Vapour-Liquid Equilibrium of Binary Mixture of 1, 2-Di-chloroethane with 1-Heptanol at Kpa

Size Exclusion Chromatography and Nuclear Magnetic Resonance Techniques for Predicting Asphalt Yields and Viscosities from Crude Oil Analyses

Chapter 3. Volumetric Properties of Pure Fluids

What s New in ChemSep TM 6.8

A Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil

KINETICS OF IRON OXIDE DIRECT REDUCTION BY COAL E.R. ABRIL 1

PREDICTION OF THE DENSITIES AND PRESSURES OF ETHANE ON THE COEXISTENCE

New Mexico Institute of Mining and Technology

Determination of the reaction order

Flash point of organic binary mixtures containing alcohols: experiment and prediction

Chemistry (Physical chemistry) Lecture 10.

Process engineers are often faced with the task of

Analysis of discretization in the direct simulation Monte Carlo

Generalization of a Double-Point Method to Determine the Intrinsic Viscosity in a Polymer-Solvent Mixture

PHASE EQUILIBRIA OF CARBON DIOXIDE AND METHANE GAS-HYDRATES PREDICTED WITH THE MODIFIED ANALYTICAL S-L-V EQUATION OF STATE

JF Physical Chemistry JF CH 1101: Introduction to Physical Chemistry.

Thermodynamic Properties of Supercritical Fluids: Example of n-hexane

A Stochastic Analysis of Liquid Mixing in Bubble Column

Replacement of gas phase with liquid, for hexamine production

Pure Component Phase Diagram. Definitions. Definitions (cont.) Class 17 Non-Ideal Gases

Review of classical thermodynamics

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Temperature Control of Batch Suspension Polyvinyl Chloride Reactors

CEE 670 TRANSPORT PROCESSES IN ENVIRONMENTAL AND WATER RESOURCES ENGINEERING. Kinetics Lecture #1

Isobaric Vapor Liquid Equilibria of Systems containing N-Alkanes and Alkoxyethanols

Introduction (1) where ij denotes the interaction energy due to attractive force between i and j molecules and given by; (2)

Vapor liquid equilibria for the binary system 2,2 dimethylbutane + 1,1 dimethylpropyl methyl ether (TAME) at , , and 338.

QCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines

Modeling of Corrosion Kinetics of Mild Steel in Hydrochloric Acid in the Presence and Absence of a Drug Inhibitor

Boiling points for five binary systems of sulfolane with aromatic hydrocarbons at kpa

19.3 SPECTROSCOPY OF ALDEHYDES AND KETONES

Robust Flight Control Design for a Turn Coordination System with Parameter Uncertainties

Pneumatic Performance Study of a High Pressure Ejection Device Based on Real Specific Energy and Specific Enthalpy

Kinetics of Esterification of Ethylene Glycol with Acetic Acid Using Cation Exchange Resin Catalyst

Lab.6. Extraction. Upper phase. water. Bottom phase. chloroform

PRACTICAL DATA CORRELATION OF FLASHPOINTS OF BINARY MIXTURES BY A RECIPROCAL FUNCTION: THE CONCEPT AND NUMERICAL EXAMPLES

2. Mass transfer takes place in the two contacting phases as in extraction and absorption.

Chapter 14. The Concept of Equilibrium and the Equilibrium Constant. We have for the most part depicted reactions as going one way.

EXPT. 3 DETERMINATION OF CHROMIUM AND MANGANESE IN A MIXTURE

THALES Project No. 65/1181

Density and surface tension of liquid ternary Ni-Cu-Fe alloys

EXTRACTIVE DISTILLATION OF ETHYLBENZENE AND STYRENE USING SULFOLANE AS SOLVENT: LOW PRESSURE ISOBARIC VLE DATA

PREDICTION OF CONCRETE COMPRESSIVE STRENGTH

Chemical Engineering Thermodynamics II ( ) 02 - The Molar Gibbs Free Energy & Fugacity of a Pure Component

METU Chem. Eng. Dept. Ch.E. 320 Chem. Eng. Lab I EXPERIMENT 33 MOLECULAR WEIGHT DETERMINATION OF POLYMERS BY CAPILLARY FLOW METHOD

- - I Thermodynamic I. Evaluation of Surface Tension of Molten Salt Mixtures based on Thermodynamic Databases. } BulkPha~

Definitions. Pure Component Phase Diagram. Definitions (cont.) Class 16 Non-Ideal Gases

13.4 REDOX STOICHIOMETRY

Physics of Relaxation. Outline

Calibration of Piping Assessment Models in the Netherlands

Is classical energy equation adequate for convective heat transfer in nanofluids? Citation Advances In Mechanical Engineering, 2010, v.

Extraction of Phenol from Industrial Water Using Different Solvents

STATIC SIMULATION PROGRAM OF COPPER SOLVENT EXTRACTION CONFIGURATIONS USING MICROSOFT EXCEL SOLVER

Final Exam: know your section, bring your ID!

Surfactants in cloud droplet activation: mixed organic-inorganic particles

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012

Numerical simulation of a one-dimensional shock tube problem at supercritical fluid conditions

Agung Ari Wibowo, S.T, M.Sc THERMODYNAMICS MODEL

BINARY RANKINE CYCLE OPTIMIZATION Golub, M., Koscak-Kolin, S., Kurevija, T.

ULTRASOUND PIEZOCERAMIC TRANSDUCER FOR UNDERWATER ACOUSTICS IN SONOCHEMISTRY

MATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II

Correspondence should be addressed to Harmandeep Singh Gill; harman

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION

IV Transport Phenomena. Lecture 21: Solids and Concentrated Solutions

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry

Possibility of Refuse Derived Fuel Fire Inception by Spontaneous Ignition

Copyright 2018 Society of Photo-Optical Instrumentation Engineers (SPIE). One print or electronic copy may be made for personal use only.

Homework Set 4. gas B open end

DETERMINATION OF MATERIAL PARAMETERS OF A TEXTILE REINFORCED COMPOSITE USING AN INVERSE METHOD

A Chemical Engineering Approach to Cellulose Substitution Kinetics

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED

Neuro-Fuzzy Control of Chemical Reactor with Disturbances

Supplementary Information. Infrared Transparent Visible Opaque Fabrics (ITVOF) for Personal Cooling

11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules

DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS

Neuro-Fuzzy Modeling of Heat Recovery Steam Generator

Ranking of Bangladeshi Coals Based on Fuzzy Set Theory

Maximum Sauter Mean Diameter and Terminal Velocity of Drops in a Liquid-liquid Spray Extraction Column

Performance of an irreversible Diesel cycle under variable stroke length and compression ratio

ABSTRACT. and since the phase transition of asymmetric one-component fluid or mixture

Phase equilibria properties of binary and ternary systems containing isopropyl ether + isobutanol + benzene at K.

23.1 Tuning controllers, in the large view Quoting from Section 16.7:

Inter-fibre contacts in random fibrous materials: experimental verification of theoretical dependence on porosity and fibre width

Damage Evaluation of Core Concrete by AE

DETERMINATION OF CARRIER TEMPERATURE FROM JUNCTION I(V) MEASUREMENTS

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2011

On scaling laws at the phase transition of systems with divergent order parameter and/or internal length : the example of DNA denaturation

Transcription:

International Journal of Chemial and Biologial Engineering 3:3 1 Isobari Vapor-Liquid Equilibrium data for Binary Mixtures of n-butylamine and Triethylamine with Cumene at 97.3 kpa Baljinder K. Gill, V. K. Rattan, and Seema Kapoor 1 Abstrat Isobari vapor-liquid equilibrium measurements are reported for the binary mixtures of n-butylamine and Triethylamine with Cumene at 97.3 kpa. The measurements have been performed using a vapor reirulating type (modified Othmer's) equilibrium still. The binary mixture of n-butylamine + Cumene shows positive deviation from ideality. Triethylamine + Cumene mixture shows negligible deviation from ideality. None of the systems form an azeotrope. The ativity oeffiients have been alulated taking into onsideration the vapor phase nonideality. The data satisfy the thermodynami onsisteny test of Herington. The ativity oeffiients have been satisfatorily orrelated by means of the Margules, NRTL, and Blak equations. The ativity oeffiient values obtained by the UNIFAC model are also reported. Keywords Binary mixture, Cumene, n-butylamine, Triethylamine, Vapor-liquid equilibrium. I. INTRODUCTION EPARATION of liquid mixtures by distillation is one of Sthe most important proesses in hemial industries. For the design of distillation olumns, knowledge of vaporliquid equilibrium data is of utmost importane. Due to omplex vapor-liquid equilibrium problems arising from numerous new industrial proesses, there is a need for the aurate vapor-liquid equilibrium determinations experimentally. In addition to this, experimental data are required to update and improve the data bank used to fit the model parameters of various theoretial models. Very limited work has been reported on vapor-liquid equilibrium study of binary mixtures ontaining umene as one of the omponents. Isobari vapor-liquid equilibrium data for binary mixtures of N-Methylaetamide and N,N-Dimethylaetamide with Cumene has been studied and reported in [1]. In the present work, experimental vapor-liquid equilibrium data for binary mixtures of n-butylamine and Triethylamine with umene are reported. The measurements were performed under isobari onditions at a pressure of 97.3 kpa using a modified version of the reirulating type equilibrium still that has been Baljinder K. Gill is with the Department of Chemial Engineering, Beant College of Engineering and Tehnology, Gurdaspur-14351, India (orresponding author; phone: 91 981599884; fax: +91-1874-1463; e- mail: bkg-7@hotmail.om) Dr. V. K. Rattan is with the Univ. Institute of Chemial Engineering & Tehnology, Panjab University, Chandigarh-1614, India (e-mail: vkrattan@pu.a.in). Dr.Seema Kapoor is with the Univ. Institute of Chemial Engineering & Tehnology, Panjab University, Chandigarh-1614, India (e-mail: seemakap_4@sify.om). desribed earlier [], [3]. The binary system n-butylamine + Cumene has a wide boiling range i.e. 75 K and the other binary system Triethylamine + Cumene has a boiling range of 6.6 K. None of the systems form an azeotrope. The ompounds studied have a wide range of appliations and are of great industrial importane. Cumene is used to manufature other hemials suh as phenol, aetone, aetophenone, and methyl styrene. It is used as a thinner in paints, laquers, and enamels. Also, it is a omponent of highotane motor fuels. Natural soures of isopropylbenzene inlude rude petroleum and oal tar. n-butylamine is used as an intermediate in the synthesis of dyes, drugs, rubber additives, emulsifiers, tanning agents and insetiides. It is also used as a vulanizing aelerator for rubber and as a uring agent for polymers. Triethylamine is ommonly employed in organi synthesis as a base, most often in the preparation of esters and amides with ayl hlorides. It is also used in the synthesis of pestiides, pharmaeutials, paints and oatings. Other appliations of Triethylamine inlude uring, hardening and orrosion inhibition for polymers and its use as a propellant. II. EXPERIMENTAL Chemials: n-butylamine and Triethylamine were obtained from C.D.H (P) Ltd., India and Cumene was obtained from Merk-Shuhardt, Germany. All hemials were AR grade materials and had purities (by hromatographi analysis, as given by the manufaturer in area perent) of 98. %, 99. % and 99. % respetively. The hemials were purified using standard proedures [4] and stored over moleular sieves. The purity of the hemials was heked by measuring the refrative indies for the pure ompounds and omparing them with the values reported in the literature. The results are listed in Table I. Apparatus and Proedure: The vapor-liquid equilibrium data were obtained by using a modified version of the equilibrium still. The equilibrated mixtures were analyzed using a Baush and Lomb Abbe-3L refratometer. The apparatus, modifiations, and analytial tehniques have already been desribed earlier [5]. All the measurements were made at a onstant temperature with the help of a irulating-type ryostat (type MK7, MLW, Germany) maintained at a temperature within ±. K. The estimated unertainties in the measurements of mole fration were ±., in refrative index were ±., and in temperature were ±. K. 14

International Journal of Chemial and Biologial Engineering 3:3 1 III. RESULTS AND DISCUSSION The liquid-phase ativity oeffiients ( ) were alulated from the experimental data using the equations [6] below, whih take into aount the vapor phase nonideality: 1 ( P y1 / P1 )exp[{( B11 V1 )( P P1 )/ RT} ( P1y )/ RT] (1) ( Py / P x)exp[{( B V )( PP )/ RT} ( P1y1 )/ RT] () 1 B 1 B 11 B (3) where x 1, x and y 1, y are the equilibrium mole frations of omponents 1 and in the liquid and vapor phases, respetively; T and P are the boiling point and the total pressure; V 1 and V are the molar liquid volumes; B 11 and B are the seond virial oeffiients of the pure omponents; and B 1 is the ross seond virial oeffiient. Table II gives the physial onstants of the pure omponents. The pure omponent vapor pressures ( P ) were alulated aording to the Antoine equation: Log ( P /.133 ) A B /( C T 73.15) (4) The Antoine s onstants A, B, and C are reported along with physial onstants of pure omponents in Table II. The experimental vapor-liquid equilibrium data ( T, x 1, and y 1 ) at 97.3 kpa along with the alulated ativity oeffiients for n-butylamine + Cumene are presented in Table III and for Triethylamine + Cumene are presented in Table IV. The Yen and Woods [7] method was used for the estimation of liquid molar volumes. The Pitzer and Curl equation modified by Tsonopoulos [8] was used in the evaluation of seond virial oeffiient as well as ross virial oeffiients in this work. The data for the systems were assessed for thermodynami onsisteny by applying the Herington area test [9]. It shows that the experimental data are thermodynamially onsistent. The ativity oeffiients were orrelated with Margules, NRTL [1], and Blak equations. The adjustable parameter α 1 for the NRTL orrelation equation was set equal to.4 for the n-butylamine + Cumene system and was set equal to.41 for the Triethylamine + Cumene system. The estimation of parameters for the three orrelation equations is based on minimization of ln( 1 / ) as an objetive funtion using the nonlinear least square method of Nagahama, Suzuki, and Hirata as used by Rattan et al [11]. The orrelation parameters A 1, A, A and deviation in vapor phase omposition for 3 the n-butylamine + Cumene system and the Triethylamine + Cumene system are listed in Table V and Table VI respetively. The Blak equation gave the best fit with.375 as the average absolute deviation in the vapor phase omposition of n-butylamine for the n-butylamine + Cumene system. The Margules equation gave the best fit with.359 as the average absolute deviation in the vapor phase omposition of triethylamine for the Triethylamine + Cumene system. TABLE I REFRACTIVE INDEX, nd AT 98.15 K n D Compound Exptl. Lit. n-butylamine 1.39839 1.3987 [4] Triethylaminee 1.39819 1.398 [4] Cumene 1.4889 1.4889 [4] TABLE II PHYSICAL CONSTANTS OF THE PURE COMPOUNDS Constant n-butylamine Triethylamine Cumene Moleular wt 73.14 [1] 11.193 [1] 1. [1] Boiling point at 11.3 kpa (K) 349.5 [13] 36.5 [13] 45.6 [13] 1.3987 [4] Refrative index, nd at 98.15 K 1.398 [4] 1.4889 [4] 531.9 [13] 535. [13] 631.13 [14] P (kpa) 4198.9 [13] 39.3 [13] 38.1 [13] V 1 6 (m 3 mol -1 ) 77. [15] 389. [13] 48. [14] Aentri fator,.39 [13].319 [1].35 [1] Dipole moment, (Debyes) 1.37 [4].87 [4].39 [4] Constants of Antoine s equation, eq.4 A 6.94519[16] 7.18658 [16] 6.9316 [16] B 11157.81 [16] 1341.3 [16] 1457.318 [16] C 7.8[ 16]. [16] 7.37 [16] TABLE III VAPOR LIQUID EQUILIBRIUM DATA OF THE BUTYLAMINE (1) + CUMENE () SYSTEM y1 ln 1 ln 349..988.9984.48.3445 349.73.9698.9961.356.318 351.8.9188.989.836.897 35.19.885.984.113.5198 354.35.815.974.1199.576 356.51.7478.968.93.17938 36.5.6493.946.351.1499 361.75.6137.9338.34.13186 363.7.5661.917.4591.11156 368.95.466.887.5548.819 373.56.39.851.6965.5818 378.74.3144.86.9516.3879 384.97.46.747.11485.77 389.38.1997.6868.187.91 394.5.1583.65.15379.173 397.58.1317.5688.168.1415 46.1.793.4196.18965.119 413.85.46.548.1998.79 418.3.16.1487.397.461 4.87.4.99.99.17 143

International Journal of Chemial and Biologial Engineering 3:3 1 TABLE IV VAPOR LIQUID EQUILIBRIUM DATA OF THE TRIETHYLAMINE (1) + CUMENE () SYSTEM y1 ln 1 ln 36.45.9686.9949 -.5.443 363.69.913.9866 -.66.448 364.56.89.989 -.31.5186 365.88.855.976 -.18.3999 367.67.7979.966 -.11.375 368.75.7684.9534 -.46.375 37.96.79.937 -.93.681 373.5.6575.95 -.48.375 375.77.5944.8973 -.447.1987 379.33.5197.8636 -.54.1637 38.5.4676.8355 -.499.181 384.86.4195.85 -.676.585 388.16.3647.764 -.346.45 391.95.31.715 -.143 -.637 395.45.68.663.314 -.933 41.17.194.5631.1646 -.69 46.66.135.4531.758 -.46 411.93.888.3388.48 -.965 416.9.54.9.639 -.34 418.84.336.1546.8411 -.35 41.96.13.646.14 -.8 Mole fration of n -Butylamine in vapor phase 1..9.8.7.6.5.4.3..1. 44 4 4 Mole fration of n -Butylamine in liquid phase Fig. 1. VLE of the n-butylamine + Cumene system at 97.3 kpa. TABLE V CORRELATION PARAMETERS FOR ACTIVITY COEFFICIENT AND DEVIATION IN VAPOR-PHASE COMPOSITION FOR THE BUTYLAMINE (1) + CUMENE () SYSTEM Correlations A 1 A A 3 Deviation (Δy) Margules.96.3195.881.394 NRTL.5677 -.334 -.411 Blak.146.365.14.375 TABLE VI CORRELATION PARAMETERS FOR ACTIVITY COEFFICIENT AND DEVIATION IN VAPOR-PHASE COMPOSITION FOR THE TRIETHYLAMINE (1) + CUMENE () SYSTEM Correlations A 1 A A 3 Deviation (Δy) Margules.11.533.874.359 NRTL -.4686.6619 -.366 Blak.948.385.1.511 Temperature(K) 38 36 34 3 3 Composition Fig.. Temperature vs. Composition urves for the n-butylamine + Cumene system at 97.3 kpa. Fig. 1 shows the experimental vapor liquid equilibrium data for the n-butylamine + Cumene binary mixture. In Fig., the Temperature vs. Composition urves are drawn for the n- Butylamine + Cumene system at 97.3 kpa. Fig. 3 shows the plot of ln value of ativity oeffiients as obtained by UNIFAC method [17] vs. omposition for the n- Butylamine + Cumene system. The graph learly indiates positive deviation from ideal behavior for the binary system studied. The mixture does not form an azeotrope. Fig. 4 shows the plot of y vs. y for the n-butylamine + Cumene system using NRTL equation. 144

International Journal of Chemial and Biologial Engineering 3:3 1.5 1..4.3..1. M ole fration of n-butylamine in liquid phase Mole fration of triethylamine in vapor phase.9.8.7.6.5.4.3..1. Mole fration of triethylamine in liquid phase Fig. 3. Plot of ln ln vs. omposition for the n-butylamine + Cumene system at 97.3 kpa., UNIFAC. Fig. 5. VLE of the Triethylamine + Cumene system at 97.3 kpa. 43. Mole fration of n-butylamine in vapor phase (NRTL) 1..9.8.7.6.5.4.3..1. Mole fration of n-butylamine in vapor phase (exptl.) Temperature (K) 4. 41. 4. 39. 38. 37. 36. 35. Composition Fig. 4. Plot of y vs. y for the n-butylamine + Cumene system using NRTL equation Fig. 6. Temperature vs. Composition urves for the Triethylamine + Cumene system at 97.3 kpa. 145

International Journal of Chemial and Biologial Engineering 3:3 1 ln1, ln (Ativity oeffiients).1.5. -.5 IV. CONCLUSION Vapor-Liquid equilibrium data at P=97.3 kpa has been obtained for the two binary systems; n-butylamine + Cumene and Triethylamine + Cumene. Both mixtures are nonazeotropi in nature. Whereas the n-butylamine + Cumene binary mixture shows positive deviation from ideality, the Triethylamine + Cumene binary mixture shows negligible deviation from ideality. This work will be of great use in improving the databank for estimation of model parameters for mixtures formed by umene with amines, and thus will enhane the preditability of the group ontribution model. Mole fration of Triethylamine in vapor phase (NRTL) -.1 1..9.8.7.6.5.4.3..1 Mole fration of Triethylamine in liquid phase Fig. 7. Plot of ln ln vs. omposition for the Triethylamine + Cumene system at 97.3 kpa., UNIFAC.. Mole fration of Triethylamine in vapor phase (exptl.) Fig. 8. Plot of y vs. y for the Triethylamine + Cumene system using NRTL equation Fig. 5 shows the experimental vapor liquid equilibrium data for the Triethylamine + Cumene binary mixture. In Fig. 6, the Temperature vs. Composition urves are drawn for the Triethylamine + Cumene system at 97.3 kpa. Fig. 7 shows the plot of ln value of ativity oeffiients as obtained by UNIFAC method [17] vs. omposition for the Triethylamine + Cumene system. The graph shows nearly ideal behavior for the binary system studied. The mixture does not form an azeotrope. Fig. 8 shows the plot of y vs. y for the Triethylamine + Cumene system using NRTL equation. REFERENCES [1] B. K. Gill, V. K. Rattan, and S. Kapoor, Vapor liquid equilibrium data for N-Methylaetamide and N,N-Dimethylaetamide with Cumene at 97.3 kpa, Journal of Chemial and Engineering Data, 9, vol. 54, pp. 1175-1178. [] B. N. Raju, R. Ranganathan, and M. N. Rao, Vapor-liquid equilibrium still for partially misible systems, Indian Chemial Engineer, 1965, vol. 7, pp. T33 T37. [3] B. Kumar, and K. S. N. Raju, Vapor-liquid equilibrium data for the systems -methoxyethanol-ethylbenzene, -methoxyethanol-pxylene, and -ethoxyethanol-p-xylene, Journal of Chemial and Engineering Data, 1977, vol., pp.134 137. [4] J. A. Riddik, W. B. Bunger, and T. K. Sakano, Organi Solvents: Physial Properties and Methods of Purifiation. 4th ed. Wiley- Intersiene: New York, 1986. [5] B. K. Sood, O. P. Bagga, and K. S. N. Raju, Vapor-liquid equilibrium data for systems ethylbenzene-anisole and p-xylene anisole, Journal of Chemial and Engineering Data, 197, vol. 17, pp. 435 438. [6] H. C. Van Ness, and M. M. Abbott, Classial Thermodynamis of Non-eletrolyte Solutions. MGraw-Hill: New York, 198. [7] C. L. Yen, and S. S. Woods, A Generalized equation for omputer alulation of liquid densities, AIChE Journal, 1966, vol. 1, pp. 95 99. [8] C. Tsonopoulos, An empirial orrelation of seond virial oeffiients, AIChE Journal, 1974, vol., pp. 63 7. [9] E. F. G. Herington, Tests for the onsisteny of experimental isobari vapor-liquid equilibrium data, Journal of Institute of Petroleum, 1951, vol. 37, pp. 457 47. [1] H. Renon, and J. M. Prausnitz, Loal ompositions in thermodynami exess funtions for liquid mixtures, AIChE Journal, 1968, vol. 14, pp. 135 144. [11] V. K. Rattan, S. Kapoor, and S. Singh, Isobari vapor-liquid equilibria of 1-butanol- p-xylene system, International Journal of Thermophysis, 6, vol. 7, pp. 85 91. [1] R. M. Stephenson, and S. Malanowski, Handbook of the thermodynamis of organi ompounds. Elsevier Publiations, 1987. [13] R. C. Reid, J. M. Prausnitz, and B. E. Poling, The Properties of Gases & Liquids. 4th ed. MGraw-Hill: New York, 1987. [14] J. A. Riddik, W. B. Bunger, and T. K. Sakano, Organi Solvents: Physial Properties and Methods of Purifiation. 3rd ed. Wiley- Intersiene: New York, 197. [15] D. R. Lide, CRC Handbook of Physis and Chemistry. 85th ed. CRC Press, 4-5. [16] T. Boublik, V. Fried, and E. Hala, The Vapor Pressures of Pure Substanes. Elsevier: New York, 1975. [17] J. Gmehling, J. Lohmann, and R. Wittig, Vapor-liquid equilibria by UNIFAC Group Contribution. 6. Revision and extension, Industrial Engineering Chemistry Researh, 3, vol. 4, pp.183 188. 146