Liquid-Liquid Equilibria for Ternary Mixtures of (Water + Carboxylic Acid+ MIBK), Experimental, Simulation, and Optimization

Similar documents
PREDICTION OF VAPOR-LIQUID EQUILIBRIA OF BINARY MIXTURES USING QUANTUM CALCULATIONS AND ACTIVITY COEFFICIENT MODELS

Functions of Random Variables

ENGI 4421 Propagation of Error Page 8-01

Multiple Choice Test. Chapter Adequacy of Models for Regression

Research Journal of Chemical Sciences ISSN X Vol. 5(6), 64-72, June (2015)

Lecture 12 APPROXIMATION OF FIRST ORDER DERIVATIVES

We have already referred to a certain reaction, which takes place at high temperature after rich combustion.

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

Summary of the lecture in Biostatistics

A Helmholtz energy equation of state for calculating the thermodynamic properties of fluid mixtures

Module 1 : The equation of continuity. Lecture 5: Conservation of Mass for each species. & Fick s Law

Unsupervised Learning and Other Neural Networks

Evaluation of uncertainty in measurements

CHEMICAL EQUILIBRIA BETWEEN THE ION EXCHANGER AND GAS PHASE. Vladimir Soldatov and Eugeny Kosandrovich

Unimodality Tests for Global Optimization of Single Variable Functions Using Statistical Methods

General Method for Calculating Chemical Equilibrium Composition

1. A real number x is represented approximately by , and we are told that the relative error is 0.1 %. What is x? Note: There are two answers.

UNIVERSITY OF OSLO DEPARTMENT OF ECONOMICS

Quantitative analysis requires : sound knowledge of chemistry : possibility of interferences WHY do we need to use STATISTICS in Anal. Chem.?

Homework Assignment Number Eight Solutions

Manipulator Dynamics. Amirkabir University of Technology Computer Engineering & Information Technology Department

Derivation of 3-Point Block Method Formula for Solving First Order Stiff Ordinary Differential Equations

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

ESS Line Fitting

Econometric Methods. Review of Estimation

hp calculators HP 30S Statistics Averages and Standard Deviations Average and Standard Deviation Practice Finding Averages and Standard Deviations

Simple Linear Regression

Heat-Integrated Distillation Columns -Analysis and Modellingwith. Advanced Distillation as Supporting Subject

Analyzing Two-Dimensional Data. Analyzing Two-Dimensional Data

L5 Polynomial / Spline Curves

Chapter 5 Properties of a Random Sample

LINEAR REGRESSION ANALYSIS

Correlation and Regression Analysis

Solving Constrained Flow-Shop Scheduling. Problems with Three Machines

Maximum Likelihood Estimation

Descriptive Statistics

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

BAYESIAN INFERENCES FOR TWO PARAMETER WEIBULL DISTRIBUTION

Confidence Intervals for Double Exponential Distribution: A Simulation Approach

12.2 Estimating Model parameters Assumptions: ox and y are related according to the simple linear regression model

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

STA 108 Applied Linear Models: Regression Analysis Spring Solution for Homework #1

The number of observed cases The number of parameters. ith case of the dichotomous dependent variable. the ith case of the jth parameter

Objectives of Multiple Regression

Arithmetic Mean and Geometric Mean

Bayes Estimator for Exponential Distribution with Extension of Jeffery Prior Information

Lecture 3 Probability review (cont d)

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

Multivariate Transformation of Variables and Maximum Likelihood Estimation

SPECIAL CONSIDERATIONS FOR VOLUMETRIC Z-TEST FOR PROPORTIONS

X ε ) = 0, or equivalently, lim

Bayesian Inferences for Two Parameter Weibull Distribution Kipkoech W. Cheruiyot 1, Abel Ouko 2, Emily Kirimi 3

Chapter 8. Inferences about More Than Two Population Central Values

Vapor-Liquid Equilibria for Water+Hydrochloric Acid+Magnesium Chloride and Water+Hydrochloric Acid+Calcium Chloride Systems at Atmospheric Pressure

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

On generalized fuzzy mean code word lengths. Department of Mathematics, Jaypee University of Engineering and Technology, Guna, Madhya Pradesh, India

(Monte Carlo) Resampling Technique in Validity Testing and Reliability Testing

ENGI 3423 Simple Linear Regression Page 12-01

Multiple Regression. More than 2 variables! Grade on Final. Multiple Regression 11/21/2012. Exam 2 Grades. Exam 2 Re-grades

Carbonyl Groups. University of Chemical Technology, Beijing , PR China;

A New Method for Solving Fuzzy Linear. Programming by Solving Linear Programming

CS 2750 Machine Learning. Lecture 8. Linear regression. CS 2750 Machine Learning. Linear regression. is a linear combination of input components x


A Robust Total Least Mean Square Algorithm For Nonlinear Adaptive Filter

( ) = ( ) ( ) Chapter 13 Asymptotic Theory and Stochastic Regressors. Stochastic regressors model

Module 7. Lecture 7: Statistical parameter estimation

Bayes Interval Estimation for binomial proportion and difference of two binomial proportions with Simulation Study

Bias Correction in Estimation of the Population Correlation Coefficient

Chapter 8: Statistical Analysis of Simulated Data

Lecture 07: Poles and Zeros

Applying the condition for equilibrium to this equilibrium, we get (1) n i i =, r G and 5 i

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

Simulation Output Analysis

Lecture 7. Confidence Intervals and Hypothesis Tests in the Simple CLR Model

Journal of Polytechnic, 2017; 20 (1) : Journal of Polytechnic, 2017; 20 (1) :

STATISTICAL PROPERTIES OF LEAST SQUARES ESTIMATORS. x, where. = y - ˆ " 1

A Study of the Reproducibility of Measurements with HUR Leg Extension/Curl Research Line

Chapter 5 Elementary Statistics, Empirical Probability Distributions, and More on Simulation

Finite Difference Approximations for Fractional Reaction-Diffusion Equations and the Application In PM2.5

Introduction to local (nonparametric) density estimation. methods

Module 7: Probability and Statistics

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

{ }{ ( )} (, ) = ( ) ( ) ( ) Chapter 14 Exercises in Sampling Theory. Exercise 1 (Simple random sampling): Solution:

G S Power Flow Solution

Line Fitting and Regression

Chapter 4 Multiple Random Variables

A Method for Damping Estimation Based On Least Square Fit

BERNSTEIN COLLOCATION METHOD FOR SOLVING NONLINEAR DIFFERENTIAL EQUATIONS. Aysegul Akyuz Dascioglu and Nese Isler

Department of Agricultural Economics. PhD Qualifier Examination. August 2011

Study of Correlation using Bayes Approach under bivariate Distributions

Third handout: On the Gini Index

Chapter 9 Jordan Block Matrices

DKA method for single variable holomorphic functions

Numerical Simulations of the Complex Modied Korteweg-de Vries Equation. Thiab R. Taha. The University of Georgia. Abstract

Assignment 5/MATH 247/Winter Due: Friday, February 19 in class (!) (answers will be posted right after class)

Analysis of Lagrange Interpolation Formula

Linear Regression with One Regressor

ECONOMETRIC THEORY. MODULE VIII Lecture - 26 Heteroskedasticity

Comparison of Parameters of Lognormal Distribution Based On the Classical and Posterior Estimates

Transcription:

World Academy of Scece, Egeerg ad Techology teratoal Joural of Chemcal ad Molecular Egeerg Vol:7, No:6, 03 Lqud-Lqud Equlbra for Terary Mtures of (Water + Carboylc Acd+ MBK, Epermetal, Smulato, ad Optmzato D. Laad, A. Hassee, ad A. Merzougu Dgtal Ope Scece de, Chemcal ad Molecular Egeerg Vol:7, No:6, 03 waset.org/publcato/443 Abstract ths work, Epermetal te-le results ad solublty (bodal curves were obtaed for the terary systems (water + acetc acd + methyl sobutyl ketoe (MBK, (water + lactc acd+ methyl sobutyl ketoe at T = 94.5K ad atmospherc pressure. The cosstecy of the values of the epermetal te-les was determed through the Othmer-Tobas ad Hads correlatos. For the etracto effectveess of solvets, the dstrbuto ad selectvty curves were plotted. addto, these epermetal tele data were also correlated wth NRTL model. The teracto parameters for the NRTL model were retreved from the obtaed epermetal results by meas of a combato of the homotopy method ad the geetc algorthms. Keywords Lqud-lqud equlbra, homotopy methods, carboylc acd, NRTL. M. NTRODUCTON ANY attempts have bee made to descrbe the solvet etracto of carboylc acds from aqueous fermetato solutos. The (lqud lqud equlbrum (LLE measuremets ad phase behavour of terary systems cludg carboylc acds has bee the subject of much research recet years [] [5]. LLE data of the related systems are ot oly eeded for the desg of a effcet ad productve etracto system, but they are also dspesable calbrato ad verfcato of aalytcal models. the preset work, lqud lqud equlbrum data have bee obtaed for three dfferet systems, amely, (water + acetc acd + methyl sobutyl ketoe, (water + lactc acd+ methyl sobutyl ketoe at 94.5K ad at atmospherc pressure. The dstrbuto coeffcets ad separato factors were obtaed from epermetal results ad are also reported. The te les were determed ad were correlated by the methods of Othmer-Tobas ad Had o a mass-fracto bass. The epermetal results are compared wth values predcted by NRTL.. EPERMENTAL SECTON A. Chemcals Acetc acd, lactc acd, dchloromethae ad methyl sobutyl ketoe were purchased from Merck ad were of 99, 99, 98, ad 99% mass purty, respectvely. The chemcals were used wthout further purfcato. Deozed ad redstlled water was used throughout all epermets. B. Aalytcal Methods The Solublty curve was determed by the cloud pot method [6] usg a thermostated cell, equpped wth a magetc strrer ad sothermal flud jacket. The cell was kept a costat-temperature bath mataed at ±0. C. The cell was flled wth homogeeous water + carboylc acd mtures prepared by weghg, usg a Nahta YP40N balace wth a precso of 0 - g. The solvet was ttrated to the cell from a mcroburet wth a ucertaty of ±0.0cm 3. The ed pot was determed by observg the trasto from a homogeeous to a heterogeeous mture. Ths patter was coveet to provde the aqueous-rch sde of the curves. The data for orgac-rch sde of the curves were therefore obtaed by ttratg homogeeous solvet + carboylc acd bares wth water utl the turbdty had appeared. The mamum error the calculato of the compostos of the bmodal curve was estmated to be 0 -. Net, the refractve dees of these terary mtures are measured by usg a Nahta Modèle 690/ refractometer. Each measuremet was take o three occasos. For the te-le measuremet, a equlbrum cell was mmersed a thermostat cotrolled at the desred temperature (±0. o C. The pure compoets were added, ad the mture was strred for at least 3 h wth a magetc strrer. The twophase mture was allowed to settle for at least 3. Samples were take by syrge from the upper ad lower mtures. The refractve dees of both phases at equlbrum were measured to later determe ther compostos. D. Laad s wth the Departemet of dustral chemstry, Uversty Mohamed Khder Bskra, Algera (phoe: +3-06-6006663,e-mal: ayadh_jm@yahoo.com A. Hassee s wth the Departemet of dustral chemstry, Uversty Mohamed Khder Bskra, Algera (e-mal: hasee@yahoo.fr A. Merzougu s wth the Departemet of dustral chemstry, Uversty Mohamed Khder Bskra, Algera (,e-mal: merzougukarm@yahoo.com. MODELS THERMODYNAMC The NRTL model was used to correlate the epermetal data ths work, the NRTL model was epressed by: teratoal Scholarly ad Scetfc Research & ovato 7(6 03 384

World Academy of Scece, Egeerg ad Techology teratoal Joural of Chemcal ad Molecular Egeerg Vol:7, No:6, 03 Dgtal Ope Scece de, Chemcal ad Molecular Egeerg Vol:7, No:6, 03 waset.org/publcato/443 lγ = j= τ k= j j k j k j + j= jj τ j kj k k = ( α jτ j = k = τ = ep α j = α j The effectve teracto parameter τ g g R. T τ j j j = = j kj k j s defed by: The LLE epermetal data were used to determe the optmum NRTL bary teractos parameters. The thermodyamc model was ftted to epermetal data usg the Newto homotopy formulato [7]. V. HOMOTOPY CONTNUATON METHOD The phase equlbrum problem volves the separato of N C compoets of molar fractos z to two phases ad at a gve temperature T ad pressure P. The molar fractos of the compoets the two phases are deoted as, respectvely. these equatos; A T ( β z = 0 j ( ad β + (, Feed γ γ = 0 (3 = 0 (4 deote the mole fractos phases ad. β stads for the phase fracto of phase. The deotes the actvty coeffcets to be calculated by usg a approprate actvty coeffcet model such as the UNQUAC model whch s dscussed the prevous secto. ths work, a cove lear homotopy cotuato method s used; ths method s globally coverget method for the soluto of olear equatos;.e., ths method ca fd the soluto whle movg from a kow soluto or a startg pot to a ukow soluto. The homotopy fucto of a lear combato of two fuctos (, a easy fucto, ad F(, a dffcult fucto. γ H(,t = tf( + ( - t( (5 where t s a homotopy parameter that s gradually vared from 0 to as a path s tracked from ( (kow soluto to F( (ukow soluto ad s the vector of depedets varables, 0 s the tal of ad a soluto of (. For the Newto homotopy formulato, the fucto H(,t (4 wth (= F( - F( 0, ca takes a form as follows: H(, t = F( - ( - t F( 0 (6 Whe usg homotopy cotuato method, t s ot oly coveet but effectve as well, to talze the terary system wth pure methyl sobutyl ketoe ad water [7]. V. RESULTS AND DSCUSSON A. Epermetal Bodal ad Te-Le Measuremets The compostos defg the bodal curve of the terary systems, (water + acetc acd + methyl sobutyl ketoe, (water + lactc acd+ methyl sobutyl ketoe at 94.5K are lsted Table whch deotes the mass fracto of the th compoet. TABLE EPERMENTAL SOLUBLTY CURVE DATA Water ( + lactc Acd ( + MBK(3 Water ( + acetc Acd ( + MBK(3 ( ( (3 ( ( (3 0.038 0.04867 0.076 0.389 0.507 0.7453 0.983 0.406 0.47389 0.59 0.6508 0.7 0.853 0.8856 0.093 0.83 0.693 0.35 0.436 0.44773 0.4905 0.50359 0.44374 0.4608 0.3077 0.86 0.57 0.089 0.85997 0.76903 0.66547 0.5758 0.4383 0.37774 0.307 0.09305 0.0836 0.060 0.07775 0.04044 0.0353 0.0348 0.0544 0.09977 0.397 0.7393 0.7689 0.96 0.709 0.3048 0.47 0.38999 0.4093 0.45495 0.5357 0.5793 0.648 0.66787 0.744 0.866 0.90065 0.08568 0.4958 0.64 0.5353 0.856 0.347 0.3 0.3403 0.34444 0.34804 0.34365 0.34378 0.7 0.30399 0.6985 0.64 0.9446 0.855 0.07089 0.8599 0.75065 0.656 0.5753 0.53749 0.4689 0.4056 0.35559 0.3409 0.698 0.4704 0.07 0.696 0.67 0.08734 0.09849 0.0644 0.05484 0.0846 Table shows the epermetal te-le compostos of the equlbrum phases, for whch ad 3 refer to the mass fractos of the th compoet the aqueous ad solvet phases, respectvely. The bodal curves ad te-les are show Fgs. ad. TABLE EPERMENTAL TE-LNE RESULTS N MASS FRACTON FOR TERNARY SYSTEMS Water-rch phase (aqueous phase Solvet-rch phase (orgac phase (water 0,7666 0,65676 0,597 0,5785 0,4038 3 (solute (solvet (water (solute Water ( + Lactc Acd ( + MBK(3 0,9684 0,03656 0,0366 0,044 0,8344 0,0598 0,053 0,08687 0,354 0,07476 0,03588 0,0308 0,34376 0,07774 0,043 0,57 0,4595 0,0 0,0439 0,659 Water ( + Acetc Acd ( + MBK(3 0,034 0,043 0,0934 0,0376 0,095 0,03575 0,0639 0,097 3 (solvet 0,944 0,8796 0,8604 0,8053 0,790 0,9543 0,9489 0,873 0,8449 0,83085 6 0,04355 0,0796 0, 0,799 teratoal Scholarly ad Scetfc Research & ovato 7(6 03 385

World Academy of Scece, Egeerg ad Techology teratoal Joural of Chemcal ad Molecular Egeerg Vol:7, No:6, 03 0,7879 5 0,6884 0,689 0,585 0,884 0,39 0,706 0,05395 0,06566 0,0884 0,05 0,0955 0,3 0,3 0,5 0,497 0,8 0,048 0,4 48 0,7056 0,67 0,656 where D ad D are the dstrbuto coeffcets of water ad carboylc acd, respectvely. The dstrbuto coeffcets, D, for water ( = ad carboylc acd ( = were calculated as follows: Dgtal Ope Scece de, Chemcal ad Molecular Egeerg Vol:7, No:6, 03 waset.org/publcato/443 acetc acd,00,00,00 MBK water Fg. Terary dagram for LLE of {water ( + Acetc acd ( + MBK (3} at T = 94.5 K; epermetal solublty curve; epermetal te-le data; * epermetal te-le data at 308.5 K [8] Lacc Lactc acd acd,00,00,00 MBK water Fg. Terary dagram for LLE of {water ( + lactc acd ( + MBK (3} at T = 94.5 K; epermetal solublty curve; epermetal te-le data B. Dstrbuto Coeffcet ad Separato Factor The effectveess of etracto of carboylc acd by the solvet s gve by ts separato factor (S, whch s a dcato of the ablty of the solvet to separate carboylc acd from water. (7 where ad 3 D = (8 3 are the mass fractos of compoet the solvet-rch ad water-rch phases, respectvely. The dstrbuto coeffcets ad separato factors for each system are gve Table. Separato factors are greater tha for the systems reported here, whch meas that etracto of the studed carboylc acds from water wth MBK s possble. The separato factor s ot costat over the whole two-phase rego. Selectvty dagrams o a solvetfree bass are obtaed by plottg, 3 / 3 3 ( + ( versus / + for each carboylc acd Fg. 3 The selectvty dagram dcated that the performace of the MBK for etracto of acetc acd s hgher tha lactc acd. TABLE DSTRBUTON COEFFCENTS FOR WATER (D AND CARBOYLC ACD (D, AND SEPARATON FACTORS (S D D S Water ( + Lactc Acd ( + MBK(3 0,078 0,0505 0,06064 0,074 0,09899 0,007 0,0746 0,0958 0, 0,4943 0,788 0,408 0,30648 0,30905 0,443 0,3538 Water ( + acetc Acd ( + MBK(3,046,07 0,96497 0,94657 0,99595 0,976,043 6,0039 5,09654 5,93697 3,57434 49,90905 4,74 0,076 7,80945 6,6648 5,343 0,4398 0,8653 3,535 0,30584,0354 3,38545 3 /(3 +3 0,80 0,78 0,76 0,74 0,7 0,70 0,68 0,66 0,64 0,6 0,60 acetc lactque 0,0788 0,0975 0,895 0,6495 0,966 48 0,3 /( + Fg. 3 Solvet-free bass selectvty dagram of systems for ( Water-Lactc acd-mbk, ( Water-Acetc acd-mbk C. Othmer Tobas ad Had Correlatos ths study Othmer-Tobas [9], Had [0], correlatos were used to ascerta the relablty of the epermetal teratoal Scholarly ad Scetfc Research & ovato 7(6 03 386

World Academy of Scece, Egeerg ad Techology teratoal Joural of Chemcal ad Molecular Egeerg Vol:7, No:6, 03 Dgtal Ope Scece de, Chemcal ad Molecular Egeerg Vol:7, No:6, 03 waset.org/publcato/443 results for each system, where, mass fracto of water the aqueous phase; 3 ad, mass fracto of carboylc acd orgac ad aqueous phases, respectvely;, mass fracto of solvets orgac phase; a, b, a, ad b, the parameters of the Othme-Tobas correlato ad the Had correlato, respectvely, The Othmer-Tobas correlato s: l ( The Had correlato s l = a + bl = a' + b' l 3 (9 (0 The correlatos are gve Fgs. 4, 5, ad the costats of the correlatos are also gve Table V. The correlato factor (R beg appromately uty ad the learty of the plots dcate the degree of cosstecy of the measured LLE values ths study. TABLE V THE CORRELATON COEFFCENTS AND CORRELATON FACTORS FOR THE OTHMER-TOBAS AND HAND CORRELATONS System Othmer Tobas correlato Had correlato a b R a b R Water-Acetc 0,774,004 0,9935,3764 0,996 0,9943 acd-mbk Water-Lactc acd-mbk,3764,0898 0,9406 0,5993 0,675 0,967 L((- / 0,4 0,0-0,4-0,8 -, -,6 -,0 -,4 -,8-3, -4-3 - - 0 L((- / Fg. 4 Othmer Tobas plots of: ( Water-Acetc acd-mbk ad ( Water-Lactc acd-mbk terary systems L( / 0,0-0,5 -,0 -,5 -,0 -,5-3,0-3,5-4,0-4,0-3,6-3, -,8 -,4 -,0 -,6 -, -0,8-0,4 0,0 L( 3 / Fg. 5 Had plots of: ( Water-Acetc acd-mbk ad ( Water- Lactc acd-mbk terary systems V. OPTMZATON PROCEDURE The NRTL model was used to correlate the raw epermetal LLE values. the preset work, the value of the o-radomess parameter of the NRTL equato, α, was fed at 0.. The objectve fucto developed by Sorese [] was used to optmze the equlbrum models. The objectve fucto s the sum of the squares of the dfferece betwee the epermetal ad calculated data. The objectve fucto ca be defed as: mf m j = wk( k= j= = cal k ( j ep k ( j ( The observed results were used to determe the optmum NRTL ( g bary teracto eergy betwee a j par of j molecules or betwee each par of compouds (tables 5. The teracto parameters values obtaed wth geetc algorthm (A. The qualty of the correlato s measured by the root-mea square devato (RMSD. The RMSD value was calculated from the dfferece betwee the epermetal ad calculated mole fractos accordg to the followg equato: RMSD= 6 ( 3 k = j = = jk jk ( where s the umber of te-les, dcates the epermetal mole fracto, s the calculated mole fracto, ad the subscrpt dees compoets, j dees phases ad k =,,..., (te-les. The RMSD values the correlato by NRTL models for the systems studed at T = 394.5K are lsted Table V. teratoal Scholarly ad Scetfc Research & ovato 7(6 03 387

World Academy of Scece, Egeerg ad Techology teratoal Joural of Chemcal ad Molecular Egeerg Vol:7, No:6, 03 TABLE V NRTL ( α =0. BNARY NTERACTON PARAMETERS ( g j AND g j AND RMSD VALUES FOR LLE DATA OF THE TERNARY SYSTEMS AT T=94.5 K -j gj g j RMSD Water ( + Lactc Acd ( + MBK(3 - -4.04-8.796-3 -4.04-495.60-3 570.38 856.7449 Water ( + Acetc Acd ( + MBK(3-63.9883-38.045-3 80.936-38.6706-3 5.4 805.474 4.657.0-8.580.0 - Dgtal Ope Scece de, Chemcal ad Molecular Egeerg Vol:7, No:6, 03 waset.org/publcato/443 V. CONCLUSONS The epermetal te-le data of, (water + acetc acd + methyl sobutyl ketoe ad (water + lactc acd+ methyl sobutyl ketoe at T = 94.5K ad atmospherc pressure. The NRTL model was used to correlate the epermetal data. The RMSD values betwee observed ad calculated mass percets for the systems water-lactc-mbk, water-acetc acd- MBK were 4.657%, 8.580% respectvely t s apparet from the separato factors ad epermetal te-les that MBK s foud to be preferable solvet for separato of acetc ad lactc acd from aqueous solutos. The combg Homotopy cotuato ad eetc algorthms ca be appled to predct lqud-lqud equlbrum ad bary teracto parameters, respectvely for other lqud lqud systems as well as vapor lqud systems. REFERENCES [] M. Blg, Ş.Đ. Kırbaşlar,Ö. Ö zca, U. Dramur, J. Chem. Thermody. 37, 97-303, 005. [] S.Çehrel, Flud Phase Equlbra, 48, 4-8, 006. [3] L.Wag, Y.Cheg, L, J. Chem. Eg.Data, 5, 7-73, 007. [4] S. Şah, Ş. Đsmal Kırbaşlar, M.Blg, J. Chem. Thermody, 4, 97-0, 009. [5] H. haadzadeh, A. haadzadeh la, Kh. Bahrpama, R. Sarr, J. Chem. Thermody, 4, 67 73, 00. [6] J. J.Otero, J. F.Comesaňa, J. M.Correa, A.Correa, J. Chem. Eg. Data, 45, 898-90, 000. [7] J.W. Kovach, W.D.Seder, Comput. Chem. Eg.,, 593-605, 987. [8] M. ovdaraja, PL. Sabaratham ; Flud Phase Equlbra, 08, 69-9, 995. [9] T.F.Othmer, P.E.Tobas, d. Eg. Chem., 34, 693-696, 94. [0] D. B.Had, J. Phys. Chem., 34, 96-000, 930. [] J.M. Sorese, Ph.D. Thess, Techcal Uversty of Demark, Lygby, Demark, 980. teratoal Scholarly ad Scetfc Research & ovato 7(6 03 388