Peng-Robinson Equation

Aplikom Persamaan Peng-Robinson

Uraian Pada Comand Windows

Tc=405.5;%K

>> Pc=111.3;%atm

>> R=0.08026;%atm.L/mol.K

>> a=0.45724*(R^2*Tc^2/Pc)

a =

    4.3514

>> b=0.07780*(R*Tc/Pc)

b =

    0.0227

>> T=450;%K

>> Tc=405.5;%K

>> Tr=(T/Tc)

Tr =

    1.1097

>> w=0.1931;%w=faktor asentrik

>> k=0.37464+1.54226*w-0.26992*w^2

k =

    0.6624

>> x=sqrt(1.1094)

x =

    1.0533

>> y=0.6624-0.6977

y =

   -0.0353

>> A=(1-0.0353)

A =

    0.9647

>> P=56%atm;

P =

    56

>> T=450;%K

>> R=0.08026;%atm.L/mol.K

>> A=0.09647;

>> a=4.3514;

>> b=0.0227;

>> vdv=[P,(b*P-R*T),(a*A-3*b^2-2*R*T*b),-b*(b^2*P-R*T*b+a*A)];

>> v=roots(vdv)%L/mol;

v =

    0.6552

   -0.0389

    0.0059

5(a)

>> %sehingga volume molar yaitu = 0.6552, -0.0389,0.0059

>> %volume molar yang memungkinkan adalah 0.6552, sehingga faktor kompresibilitas yaitu:

>> v=0.6552;

>> R=0.08026;%atm.L/mol.K

>> T=450;%K

>> P=56;%atm

>> Z=(P*v/R*T)

Z =

  2.0572e+005

5b)  >> %Ulangi perhitungan untuk tekanan reduksi berikut : Pr = 1, 2, 4, 10, dan 20 !

>> Pr=[1 2 4 10 20]

Pr =

     1     2     4    10    20

>>  Pc=111.3;

>> P=Pr*Pc

P =

1.0e+003 *

0.1113    0.2226    0.4452    1.1130    2.2260

5c) %Bagaimana hubungan antara faktor kompresibilitas dan tekanan reduksi

(gambarkan dalam sebuah grafik) ?

P=[0.1113 0.2226 0.4452 1.1130 2.2260]

>> P1=0.1113;

>> P2=0.2226;

>> P3=0.4452;

>> P4=1.1130;

>> P5=2.2260;

>> z1=P1*v/(R*T)

z1 =

    0.0020

>> z2=P2*v/(R*T)

z2 =

    0.0040

>> z3=P3*v/(R*T)

z3 =

    0.0081

>> z4=P4*v/(R*T)

z4 =

    0.0202

>> z5=P5*v/(R*T)

z5 =

    0.0404

>> z=[0.0020 0.0040 0.0081 0.0202 0.0404];

>> Pr=[1 2 4 10 20];

>>plot(z,Pr)

Maka hasil grafiknya: