No Title |POLVER05_1 |2 P=100 R=1.987 # ideal gas constant x1=0.85 x2=0.15 # Antoine equations for vapor pressure f(P2sat) = ln(P2sat)-16.5+3850/(T-47) P2sat(0) = 50 # initial guess for saturation pressure of 2 f(P1sat) = ln(P1sat)-17+3600/(T-54) P1sat(0) = 90 # initial guess for saturation pressuer of 1 #Modified Raoult's law f(y2) = y2*P-x2*gamma2*P2sat y2(0) = 0.5 # initial guess for y2 f(y1) = y1*P-x1*gamma1*P1sat y1(0) = 0.6 #initial guess for y1 f(T) = 1-y1-y2 T(0) = 350 # initial guess for temperature # Wilson parameters a12=440 a21=1250 V1=77 V2=18 f(G21) = G21-(V1*EXP(-a21/R/T))/V2 G21(0) = 1 # initial guess for G21 f(G12) = G12-(V2*EXP(-a12/R/T))/V1 G12(0) = 1 # initial guess for G12 f(W) = W- G12/(x1+x2*G12) +G21/(x1*G21+x2) W(0) = 0.5 # initial guess for W #Activity coefficient calculations f(gamma2) = ln(gamma2)+ln(x1*G21+x2)+x1*W gamma2(0) = 1 f(gamma1) = ln(gamma1) + ln(x1+x2*G12) -x2*W gamma1(0) = 1