Scribd
Upload
20 views

%trabajamos Con Los Componentes Metanol y 1-Hepteno %componente Metanol

The document calculates vapor-liquid equilibrium for a mixture of methanol and 1-heptene at 353.15 K. It defines constants for the Antoine equation for each component and solves for the vapor pressures (pvpa and pvpb). It then calculates vapor phase compositions (y1 and y2) and total pressure (P) across a range of liquid phase compositions (x1). The results are displayed in a table and plotted on graphs.

Uploaded by

Branco R. R. Hermoza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
20 views

%trabajamos Con Los Componentes Metanol y 1-Hepteno %componente Metanol

The document calculates vapor-liquid equilibrium for a mixture of methanol and 1-heptene at 353.15 K. It defines constants for the Antoine equation for each component and solves for the vapor pressures (pvpa and pvpb). It then calculates vapor phase compositions (y1 and y2) and total pressure (P) across a range of liquid phase compositions (x1). The results are displayed in a table and plotted on graphs.

Uploaded by

Branco R. R. Hermoza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 3

clear,clc

T=353.15
%Trabajamos con los componentes metanol y 1-hepteno
%Componente metanol
aa=-8.54796;
ab=0.76982;
ac=-3.10850;
ad=1.54481;
tca=512.6;
pca=8090;
ra=1-T/tca;
syms pvpx
pvpa=solve(log(pvpx/pca)-(1-ra)^1*(aa*ra+ab*ra^1.5+ac*ra^3+ad*ra^6),pvpx)
%Componente 1-hepteno
ba=-8.26875;
bb=3.02688;
bc=-26.18709;
bd=4.33049;
tcb=537.3;
pcb=2830;
rb=1-T/tcb;
syms pvpy
pvpb=solve(log(pvpy/pcb)-(1-rb)^1*(ba*rb+bb*rb^1.5+bc*rb^3+bd*rb^6),pvpy)
x1=0:0.05:1+0.00001;
x2=1-x1;
P=x1.*pvpa+x2.*pvpb;
y1=x1.*pvpa./P;
y2=x2.*pvpb./P;
y1exp=[0 0.5696 0.644 0.6734 0.6892 0.6991 0.7058 0.7108 0.7143 0.710
0.7147 0.718 0.7209 0.7236 0.7265 0.7296 0.733 0.7572 0.78 0.8329 1];
res(:,1)=x1'; %columna 1 de la matriz
res(:,2)=y1'; %columna 2 de la martiz
res(:,3)=y2';
res(:,4)=P';
res(:,5)=y1exp';
disp('Tabla')
disp('
x1
y1
y2
P
y1(exp)')
disp(res)
figure(1)
plot(x1,P,y1,P,y1exp,P)
figure(2)
plot(x1,y1,x2,y2,x1,y1exp)

T=

353.1500

pvpa =

8090/exp(8555459402732829/2251799813685248)

pvpb =

2830/exp(5608164014944055/1125899906842624)

Tabla
x1
0

y1
0

y2

y1(exp)

1.0000 19.4332

0.0500

0.3291

0.6709 27.5160

0.5696

0.1000

0.5087

0.4913 35.5989

0.6440

0.1500

0.6219

0.3781 43.6818

0.6734

0.2000

0.6997

0.3003 51.7646

0.6892

0.2500

0.7565

0.2435 59.8475

0.6991

0.3000

0.7997

0.2003 67.9304

0.7058

0.3500

0.8338

0.1662 76.0132

0.7108

0.4000

0.8614

0.1386 84.0961

0.7143

0.4500

0.8840

0.1160 92.1790

0.7100

0.5000

0.9031

0.0969 100.2618

0.7147

0.5500

0.9193

0.0807 108.3447

0.7180

0.6000

0.9332

0.0668 116.4276

0.7209

0.6500

0.9454

0.0546 124.5104

0.7236

0.7000

0.9560

0.0440 132.5933

0.7265

0.7500

0.9655

0.0345 140.6762

0.7296

0.8000

0.9739

0.0261 148.7590

0.7330

0.8500

0.9814

0.0186 156.8419

0.7572

0.9000

0.9882

0.0118 164.9248

0.7800

0.9500

0.9944

0.0056 173.0076

0.8329

1.0000

1.0000

0 181.0905

1.0000

You might also like