*Sat, 19 Dec 2015*

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.

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

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

`y`

_{i}.P = x_{i}.γ_{i}.P_{i}^{sat}`y`

_{i}= x_{i}.γ_{i}.K_{i}

From above equation experimental γ_{1} is obtained for acetone.

`γ`

_{1}(Exp) = y_{1}/ ( x_{1}.K_{1})`K`

_{1}= [ 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(x_{1}+ (1-x_{1})*Y12) + (1-x_{1})[ Y12/(x_{1}+ Y12.(1-x_{1})) - Y21/(1 - x_{1}+ Y21.x_{1}) ]

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.

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

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.

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.

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.

Spreadsheet for Binary Interaction Parameter

Heat Loss from Insulated Pipe PT Flash Calculation using PR EOS