Estimating Binary Interaction Parameter by Regression

Wilson equation is commonly used to predict non ideality in binary mixture vapor liquid equilibrium. This article shows how to estimate binary interaction parameters used in wilson equation from experimental data by regression in excel spreadsheet.

Determine binary interaction parameters used in wilson equation for a mixture of acetone and chloroform from T x-y experimental data available at 101.33 kPa.

VLE data acetone chloroform mixture

Obtain pure component properties of acetone and chloroform from literature mainly vapor pressure data and liquid molar volume.

acetone chloroform example data

Featured Resources :

Hydrocarbon Processing the premier magazine providing job-help information to technical and management personnel in petroleum refining, gas processing, petrochemical/chemical and engineer/constructor companies throughout the world – since 1922.

Based on modified Raoult’s Law following relation is obtained.

 yi.P = xii.Pisat 
 yi   = xii.Ki

From above equation experimental γ1 is obtained for acetone.

 γ1(Exp) = y1/ ( x1.K1 )
 K1 = [ e(A - B/ (T + C))] / P

An initial value of binary interaction parameter is assumed.

 A12 = 200 cal/gmol
 A21 = 200 cal/gmol

Liquid phase γ1 is obtained from above interaction parameter and using wilson equation.

 Y12 = V2/V1.e-A12/RT
 Y21 = V1/V2.e-A21/RT
 lnγ1 = -ln(x1 + (1-x1)*Y12) + (1-x1)[ Y12/(x1 + Y12.(1-x1)) - Y21/(1 - x1 + Y21.x1) ]

Square of difference of γ1_experimental and γ1_calculated is obtained for all data points. An objective function is defined as summation of all these differences.

binary vle example data

Click on Solver in Data Ribbon (Excel 2010) to open dialog box for Solver parameters and input data as shown below.

binary vle solver box parameter

Minimize the objective function by changing values of A12, A21. Uncheck Make Unconstrained Variables Non-Negative, as these variables can take negative values. Click solve to start regression and new values of A12 and A21 are calculated.

 A12 = 157.9 cal/gmol
 A21 = -570.3 cal/gmol

Above values are then used to calculate y1 values and results are plotted to check the deviation.

plot of estimated parameters

Determine binary interaction parameters used in wilson equation for a mixture of benzene and acetonitrile from P x-y experimental data available at 318.15 °K.

VLE data benzene acetonitrile mixture

Use above steps and change formula for Y12, Y21 as temperature is fixed. After doing regression binary interaction parameters are obtained and result is plotted as following.

XY plot benzene acetonitrile mixture

Spreadsheet for Binary Interaction Parameter