UNIQUAC Activity Coefficient & Excess Gibbs Energy Binary Mixture Benzene and 2,2,4-Trimethylpentane

#UNIQUAC #thermodynamics #VLE #matlab #activity
Activity Coefficient Binary Mixture Using UNIQUAC (MATLAB).
Focus on Activity Coefficient and Excess Gibbs Energy Calculation.
Benzene(1) and 2,2,4-Trimethylpentane(2)
Verhoef, A., De Ridder, E., Degrève, J., Van Der Bruggen, B., 2011. Determination of activities in membrane processes: The UNIQUAC model expressed in mole and mass fractions. AIChE J. 57, 1889-1896. doi.org/10.100...
Weissman, S., Wood, S.E., 1960. Vapor-Liquid Equilibrium of Benzene-2,2,4-Trimethylpentane Mixtures. The Journal of Chemical Physics 32, 1153-1160. doi.org/10.106...
clc;clear;close all;
%data
xd=[0.0819;0.2192;0.3584;0.3831;0.5256;0.5924;0.8478;0.9872];
yd=[0.1869;0.4065;0.5509;0.5748;0.6786;0.7293;0.8741;0.9863];
P=[201;236.86;266.04;270.73;293.36;302.07;324.66;327.39];
P1sat=327.05; %mmHg at 55degC
P2sat=178.08; %mmHg at 55degC
%Calculate gamma from data
g1=yd.*P./(xd.*P1sat);
g2=(1-yd).*P./((1-xd).*P2sat);
lng1=log(g1);
lng2=log(g2);
%Calculate GE from data GE=RT(x1*ln(g1)+x2*ln(g2))
R=1.987; %cal/mol/K
T=55+273.15;
GEd=R.*T.*(xd.*lng1+(1-xd).*lng2);

[@MRojas-nh7vk]

Nice job, maybe do you have an example of ternary system? for example extraction

[@rahmatsunarya]

Thank you for your comment. If you have an idea about it, please let me know at [e-mail]