Problem D3 - Three Phase Equilibrium - Bubble Point   |POLVER05_0
f(x11) = x11*(beta+(1-beta)*k11/k12)-z1    #Rearrangement of Equation (9) for component i = 1
f(x12) = x12*k12-x11*k11    #Rearrangement of Equation (10) for component i = 1
f(x21) = x21*(beta+(1-beta)*k21/k22)-z2    #Rearrangement of Equation (9) for component i = 2
f(x22) = x22*k22-x21*k21    #Rearrangement of Equation (10) for component i = 2
f(t) = x11-y1+x21-y2    #Rearrangement of Equation (11)
f(beta) = (x11-x12)+(x21-x22)    #Rearrangement of Equation (12)
p1 = 10^(7.62231-1417.9/(191.15+t))    #Antoine equation for component i = 1
p2 = 10^(8.10765-1750.29/(235+t))    #Antoine equation for component i = 2
A = 1.7    #Numerator onstant in Equation (15)
B = 0.7    #Numerator constant in Equation (16)
gamma11 = 10^(A*x21*x21/((A*x11/B+x21)^2))    #Equation (15) for component i =1 in liquid phase j = 1
gamma21 = 10^(B*x11*x11/((x11+B*x21/A)^2))    #Equation (16) for component i =2 in liquid phase j = 1
gamma12 = 10^(A*x22*x22/((A*x12/B+x22)^2))    #Equation (15) for component i =1 in liquid phase j = 2
gamma22 = 10^(B*x12*x12/((x12+B*x22/A)^2))    #Equation (16) for component i =1 in liquid phase j = 2
k11 = gamma11*p1/760    #Equation (13) for component i = 1 in liquid phase j = 1
k21 = gamma21*p2/760    #Equation (13) for component i = 2 in liquid phase j = 1
k12 = gamma12*p1/760    #Equation (13) for component i = 1 in liquid phase j = 2
k22 = gamma22*p2/760    #Equation (13) for component i = 2 in liquid phase j = 2
z1 = 0.2    #Mole fraction of component i = 1 in feed
z2 = 0.8    #Mole fraction of component i = 2 in feed
y1 = k11*x11    #Equation (10) for mole fraction of i = 1 in vapor phase
y2 = k21*x21    #Equation (10) for mole fraction of i = 2 in vapor phase
x11(0)=0
x12(0)=1
x21(0)=1
x22(0)=0
t(0)=100
beta(0)=0.8