Problema 29

Código de programación Matlab

% PROBLEMA 29

clc
clear
format compact; format short g;

% Se tiene un compresor de desplazamiento positivo que maneja una mezcla de
% gases que se dan en la tabla, el flujo de gas es de 1800 ft^3/min a
% condiciones estándar. La presión inicial es de 85 lb/in^2 y 85°F y la
% presión final es de 2400 ln/in^2. Calcule lo siguiente:
% a) Número de etapas
% b) Presiones y temperaturas inter-etapa, considere un enfriamiento completo
% y una DP=6 PSI.
% c) Potencia real con una eficiencia del 85%

T=input(‘Valor de la temperatura [°K]’)
%T=302.594
T2=input(‘Valor de temperatura [°R]’)
%T2=545
Pd=input(‘Valor de la presión de descarga ‘);
%Pd=2400
Ps=input(‘Valor de la presión de succión ‘);
%Ps=85
dp=input(‘Valor de la diferencia de presiones’)
%dp=6
Q=input(‘Valor del gasto [ft/seg]’)
%Q=30
ef=input(‘Valor de la eficiencia ‘)
%ef=0.8
z=0.89
Pman=Ps-11.3
patm=14.7
Lo=1.22 %tabla

z=0.89
R=1545
disp(‘Proporciona los valores correspondientes según la sustancia’)

n=input(‘Cuantas sustancias hay en el sistema ‘)
%n=3
g=zeros(3,6);
for j=1:n

g(j,1)=input(‘Moles de la sustancia ‘);
g(j,2)=g(j,1)/100
g(j,3)=input(‘Peso Molecular de la sustancia ‘);
g(j,4)=g(j,2)*g(j,3)
g(j,5)=input(‘Presión Critica de la sustancia ‘);
g(j,6)=g(j,2)*g(j,5)
fprintf(‘\n Moles Yi PM Pc Pc*yi
\n’)

end
PM=g(:,4);
pm=sum(PM)
Pc=g(:,6);
pc=sum(Pc)
n=input(‘Cuántas sustancias hay en el sistema ‘)
h=zeros(3,3);
for j=1:n

h(j,1)=input(‘Valor de la temperatura crítica de la sustancia’);
h(j,2)=g(j,2)*h(j,1);
h(j,3)=input(‘Valor de la constante A’);
h(j,4)=input(‘Valor de la constante B’);
h(j,5)=h(j,3)+h(j,4)*T2;
h(j,6)=g(j,2)*h(j,5);
fprintf(‘\n Tc Tc fraccional A B Mcp
Mcp fraccional \n’)
h

end
tc=h(:,2);
tcp=sum(tc)
mc=h(:,6);
mcp=sum(mc)

n=input(‘Cuantas sustancias hay en el sistema ‘)
%n=3;
m=zeros(3,3);
for j=1:n
m(j,1)=input(‘Valor de VIS B] ‘);
m(j,2)=input(‘Valor de VISTO ‘);
m(j,3)=10^[(m(j,1)*(1/T-1/m(j,2)))];
m(j,4)=m(j,3)*.000671968;
m(j,5)=g(j,2)*m(j,4);
fprintf(‘\n VIS B VISTO VISCOSIDAD[CP]
VISCOSIDAD[lb/ft.s] \n’)
m
end
vis=m(:,5);
visp=sum(vis)

% Razón de compresión total
Rc=Pd/Ps
% Calculo de la razón de compresión por etapas
n=2
Rc1=(Rc)^(1/2)
n=3
Rc2=(Rc)^(1/3)

% Caso ideal sin enfriamiento
Pd1=Rc2*Ps
Pd2=Rc2*Pd1
Pd3=Rc2*Pd2

%Calculo de Ps y Pd reales
Pd1r=Pd1+(dp/2)
Ps2r=Pd1-(dp/2)
Pd2r=Pd2+(dp/2)
Ps3r=Pd2-(dp/2)
Pd3r=Pd

%Calculo de las razones de compresión reales

Rc1r=Pd1r/Ps
Rc2r=Pd2r/Ps2r
Rc3r=Pd3r/Ps3r
Rcprom=(Rc1r+Rc2r+Rc3r)/3

%Calculo de temperaturas de succión y descarga
k=mcp/(mcp-1.987)
o=(k-1)/k
Td1r=T2*Rc1r^o
Td2r=T2*Rc2r^o
Td3r=T2*Rc3r^o

%Para el calculo de la potencia

Prm=Ps/pc
Trm=T2/tcp
dc=(Ps*144*pm)/(R*T2*z)
bb=(Ps*144*k)/((k-1)*dc)
hh=[(Rc1r)^o+ (Rc2r)^o+(Rc3r)^o-3]

wf=bb*hh
Qinicial=Q*(14.7/(11.3+Pman))*((460+Ps)/520)
w=Qinicial*dc

%POTENCIA TEÓRICA
disp(‘El valor de la potencia esta dado en [HP]’)
PT=(wf*w)/550
disp(‘El valor de la potencia real está dado en [HP]’)
Preal=(PT/ef)*Lo