Implementación de la Isomerización del Butano

%--------------------ISOMERIZACIÓN DEL BUTANO-----------------%
% N-BUTANO I-BUTANO %


%BALANCE GLOBAL

%DATOS
X=input ('Introduce un valor de X = ');
T=433+43.4*X %°K
XNB1=0.995;
F4=20; %kmol/h
XNB6=input ('XNB6 = ');
nkr=48;%kmol/h
Ep=0.05;

%BALANCE

XEM1=1-XNB1;
XEM6=1-XNB6;

F1=((F4*(1-XNB6))/(XNB1-XNB6)); %kmol/h

F6=((1-XNB1)/(1-XNB6))*F1; %kmol/h%EVALUACIÓN DE LA Keq
%Se reportan dos valores conocidos para dicha reacción en fase líquida, el deltaH
%de reacción, además de tener una Keq conocida para reacción dada a 60°C
%que en esteproyecto le llamaremos kcs, y utiliazndo la ecuación de Vant
%Hoff para tener la relación de keq que usamos.

rAE=31.1; %velocidad de reacción especificakcs=3.03; %valor conocido de la Keq a 60°C
deltaH=-6900; %kJ/kmol
Ts=330;%°K
R=8.31446; %kJ/kmol°K

K=rAE*exp((deltaH/R)*((1/Ts)-(1/T)))

Kc=kcs*exp((deltaH/R)*((T-333)/(333*T)))keq=Kc/K;

Xeq=Kc/(1+Kc)

%BALANCE DE RECIRCULACIÓN

FNB2=((((F4^2*XEM6)/(nkr-XEM6))-F4-(F4/keq))/((F4/nkr)+((XEM6*F4)/(nkr-XEM6))-1));

FNB3=FNB2-F4;

FIB3=F4;FEM2=(XEM6*(FNB2-F4))/1-XEM6;

FEM3=FEM2;

F3=FIB3+FNB3+FEM3;

F2=FNB2+FEM2+F1;

XEM2=FEM2/F2;

XIB3=((1-XEM2)/(1+(F3/nkr)+(1/Xeq)));
XNB3=FEM3/F3;
XEM3=1-XIB3-XNB3;

%DIVISOR

FEM6=Ep*FEM2;FEM5=FEM2;

FNB6=Ep*FNB3;

FNB5=FNB3;

F5=FNB5+FEM5;

F7=F5-F6;

FEM7=FEM5-FEM6;
FNB7=FNB5-FNB6;

%ANÁLISIS ECONOMICO DE LOS EQUIPOS

%REACTOR

%DENSIDAD DE...