Notesale: Turn your study into money

Document Preview

Extracts from the notes are below, to see the PDF you'll receive please use the links above

---

clc;
clear all;
format long
Z= [0
...
30458 0
...
62 1223
...
54 544
...
88];%Critical pressures are in psia
w= [0
...
2033 0
...
732;
P= 2681;% pressure is in psia, guess the bubble point
Kij=[0 0
...
03; 0
...
01; 0
...
0096 0];
Nc =numel(Z);%number of components
Tr= (T)
...
/Pc;
counter=0;
%-----calculate the parameters of PR EOS (b,m,alpha, a_c,a_T)----for ii=1:Nc
b(ii)=0
...
37464+1
...
26992*(w(ii))^2;
alpha(ii)=(1+m(ii)*(1-sqrt(Tr(ii))))^2;
a_c(ii)=(0
...

difference=1;
while abs(difference)>0
...
/Pc;
for ii=1:Nc
K(ii)= (1/Pr(ii))*(exp(5
...
This is the
Wilson equation than 1
...

end
while error>0
...
5;%initial guess for n_g
delta_n_g=1;% to declair variation
while delta_n_g>10^(-9)
f_old= 0;

for jj=1:Nc
f_old = f_old+(Z(jj)*(K(jj)-1))/(1+n_g_guess*(K(jj)-1));
end
n_g = n_g_guess-0
...
0000001*f_old/(f_old-f_new);
delta_n_g=abs(n_g_new-n_g_guess);
n_g_guess=n_g_new;
end
result=n_g_guess;
for ll=1:numel(K)
x(ll)=(Z(ll))/(1+n_g_new*(K(ll)-1));
end
y=K
...
5);
g=@(Z_vap)Z_vap^3-(1-B_vap)*Z_vap^2+(A_vap-2*B_vap-3*B_vap^2)*Z_vap(A_vap*B_vap-B_vap^2-B_vap^3);
Z_vap=fzero(g,0

---