function DISENO_DEL_REACTOR_MATLAB_CLORURO_DE_ALILO
clc, clear all
disp('======================================================
=======================================================')
disp(' DISEÑO DE
REACTOR ISOTERMICO ')
disp(' REACTORCATALITICO DE LECHO
FIJO O EMPACADO MULTITUBULAR ')
disp(' PRODUCCIÓN DE
CLORURO DE ALILO ')
disp('======================================================
=======================================================')
%EL REACTOR OPERA ISOTERMICAMENTE (T CONSTANTE=510°C)
%TEMPERATURA DE OPERACION DEL REACTOR=510°C
%PRESION DE OPERACION DEL REACTOR=3.04 bar
%TEMPERATURA DE ALIMENTACION DEL FLUIDO EN FASE GASEOSA AL
REACTOR=510°C
%PRESION DE ALIMENTACION DEL FLUIDO EN FASE GASEOSA AL
REACTOR=3.04 bar
%NUMERO DE TUBOS:900
%DIAMETRO INTERNO DE LOS TUBOS:0.0375
syms w x P
%...........................................................
...............
%.........................VARIABLES DEL
DISEÑO.............................
%...........................................................
...............
NT=900;%REACTOR MULTITUBULAR N°TUBOS
Di=0.0375;%DIAMETRO INTERNO DE LOS TUBOS (m)
Ac=(3.1415*(Di^(2)))/4;%AREA TRANVERSAL DE LA TUBERIA.(m^2)
FC=33.5406/3600;%FLUJO MOLAR DEL CLORO (Kmol/s)
FP=(33.5406*2)/3600; %FLUJO MOLAR DEL PROPILENO (Kmol/s)
FT=FC+FP;%FLUJO MOLAR TOTAL (Kmol/s)
PM=(FC/FT)*70.9+(FP/FT)*42;%PESO MOLECULAR PROMEDIO DE LA
MEZCLA %Kg/Kmol
yoC=FC/FT;%FRACCION MOLAR DEL CLORO EN LA ALIMENTACION
yoP=FP/FT;%FRACCION MOLAR DEL PROPILENO EN LA ALIMENTACION
po=304; %PRESION INICIAL (KPa)
T=273+510;%TEMPERATURA DE OPERACION DEL REACTOR(°K)
R=8.3145;%CONSTANTE UNIVERSAL DE LOS GASES(KJ/Kmol K)
Dp=0.003;%DIAMETRO DE LA PARTICULA DEL CATALIZADOR(m)
pc=2000;%DENSIDAD DE LA PARTICULA DEL CATALIZADOR(Kg/m^3)
pg=2.413;%DENSIDAD DEL FLUIDO EN FASE GASEOSA A
510°C(Kg/m^3)
E=0.5;%Void fraction
u=0.000024;%VISCOSIDAD DEL FLUIDO EN FASE GASEOSA (Kg/m s)
G=(FT*PM)/(Ac*NT);%VELOCIDAD MASICA SUPERFICIAL (Kg/sm^2)
e=0;
%...........................................................
...............
% ECUACION DE DISEÑO
%...........................................................
...............
%..CINETICA DE LAS TRES
REACCIONES.........................................
%--PARA LA PRIMERA REACCION
% C3H6+Cl2-->CH2=CHCH2Cl+HCl
k1=0.322*exp(-85000/(T*R));%Kmol/[Link]^2
%--PARA LA SEGUNDA REACCION
%C3H6+Cl2-->CH2=CClCH3+HCl
k2=(1.83*10^(-5))*exp(-43800/(T*R));%Kmol/[Link]^2
%--PARA LA TERCERA REACCION
%C3H6+2Cl2-->C3H4Cl2+2HCl
k3=(1.27*10^(-3))*exp(-82100/(T*R));%Kmol/[Link]^3
dxdw=((k1+k2)*(((yoC*po)^2)*((P/po)^2)*((1-x)*(1-x)))+
(k3*(((yoC*po)^3)*((P/po)^3)*((1-x)^2)*(6.3-(1/2)*x))))/
(FC);
%...........................................................
........
%COMO EL REACTOR OPERA A T CONSTANTE (ISOTERMICO) LA
ECUACION DE ERGUN NOS
%AYUDA A CALCULAR LA PRESION A LA SALIDA DEL REACTOR
%..........ECUACION DE ERGUN (CAIDA DE
PRESION).....................
%...........................................................
........
b0=(((G*(1-E))/(pg*1000*Dp*(E^(3))))*(((150*(1-E)*u)/Dp)
+1.75*G));
a=(2*b0)/(Ac*(1-E)*pc*po);
dPdw=-(a/2)*(po/(P/po));
%.........................METODO NUMERICO RUNGE
KUTTA......................
%...........................................................
...............
xi=0;
Pi=304;
w0=0;
wf=20;
Po(1)=Pi;
xi=0;
n=200;
h=(wf-w0)/n;
i=0;
disp(' RESULTADOS DEL
DISEÑO DEL REACTOR
')
disp(' N° CONVERSION (Kg CAT./TUBO) P salida(Kpa) Caida
de P(Kpa) ALTURA L.(m) Kg [Link] V Reac.(m^3)
ALTURA Re.(m)')
disp('------------------------------------------------------
------------------------------------------------------------
-------')
while i<n;
w=w0;
x=xi;
P=Pi;
c1=eval(dxdw);
k1=eval(dPdw);
P=Pi+h*k1/2;
x=xi+h*c1/2;
w=w0+h/2;
c2=eval(dxdw);
k2=eval(dPdw);
P=Pi+h*k2/2;
x=xi+h*c2/2;
w=w0+h/2;
c3=eval(dxdw);
k3=eval(dPdw);
P=Pi+h*k3;
x=xi+h*c3;
w=w0+h;
c4=eval(dxdw);
k4=eval(dPdw);
xi=xi+h/6*(c1+2*c2+2*c3+c4);
Pi=Pi+h/6*(k1+2*k2+2*k3+k4);
w0=w0+h;
i=i+1;
wi(i+1)=w0;
yi(i+1)=xi;
Po(1+i)=Pi;
%1)CAIDA DE PRESION en KPa
DP=po-Pi;
%2)ALTURA DEL LECHO en m
H=w0/(Ac*(1-E)*pc);
% VOLUMEN DEL REACTOR en m^3
Vt=Ac*H*NT;
%3)CANTIDAD DE CATALIZADOR REQUERIDO Kg
Ca=w0*NT;
%4)ALTURA DEL REACTOR en m
HT=H*1.2;
fprintf(' %2d %10.4f %14.6f %16.6f %12.6f %13.2f %17f %12.4f
%9.2f\n',i,xi,w0,Pi,DP,H,Ca,Vt,HT)
end
disp('====================================================')
disp(' DATOS ADICIONALES ')
disp('====================================================')
disp('NUMERO DE TUBOS: 900')
disp('DIAMETRO DEL TUBO:0.0375 m')
disp('DENSIDAD DE LA PARTICULA DEL CATALIZADOR:2000 Kg/m^3')
disp('DIAMETRO DE LA PARTICULA DEL CATALIZADOR:0.003 m')
%DIAMETRO DEL REACTOR en m
At=Ac*NT;
D=sqrt((At*4)/3.1415);
disp('DIAMETRO DEL REACTOR (m):')
disp(D)
figure (1)
plot(wi,Po,'-r'), grid on, title('CAIDA DE PRESIÓN EN
FUNCIÓN A LA MASA DEL
CATALIZADOR'),ylabel('Presión(kPa)'),xlabel('w (kg/tubo)')
legend('Caida de presión kPa')
figure (2)
plot(wi,yi,'-b'), grid on, title('CONVERSIÓN EN FUNCIÓN A LA
MASA DEL CATALIZADOR'),ylabel('Conversión (x)'),xlabel('w
(kg/tubo)')
set(gca,'ytick',0:0.2:1)
disp('======================================================
=========================================================')
disp(' EFECTOS
ENERGETICOS
')
disp('======================================================
=========================================================')
% PASO N°1..............................................
%ENTALPIA DE FORMACION H°rx,298°K
%COMPONENTES(H°f,298°K (J/mol)......
P=20.41; %PROPILENO
Cl=0.000; %CLORO
CA=-5.6; %CLORURO DE ALILO
HC=-92.310; %ACIDO CLORHIDRICO
CP=-24.700;%2-CLORO PROPENO
DP=-73.470;%DICLORO PROPENO
% PASO N°2..............................................
%CALCULANDO LA ENTALPIA DE REACCION ESTANDAR H°rx,298°K
%H°rx,298°K(productos)-%H°rx,298°K(reactivos)
%--PARA LA PRIMERA REACCION
% C3H6+Cl2-->CH2=CHCH2Cl+HCl
Hf1=(CA+HC)-(P-0);
%--PARA LA SEGUNDA REACCION
%C3H6+Cl2-->CH2=CClCH3+HCl
Hf2=(CP+HC)-(P-0);
%--PARA LA TERCERA REACCION
%C3H6+2Cl2-->C3H4Cl2+2HCl
Hf3=(DP+2*HC)-(P-0);
%VALORES DE Cp (J/molK) en funcion de la T
%FORMULA Cp/R=a0+a1*T+a2*T^2+a3*T^3+a4*T^4 (J/molK)
%COMPONENTE HCl
a0C=3.827;
a1C=-2.936*10^(-3);
a2C=0.879*10^(-5);
a3C=-1.031*10^(-8);
a4C=0.439*10^(-11);
%COMPONENTE PROPILENO
a0P=3.834;
a1P=3.893*10^(-3);
a2P=4.688*10^(-5);
a3P=-6.013*10^(-8);
a4P=2.283*10^(-11);
%COMPONENTE CLORO
a0L=3.0560;
a1L=5.3708*10^(-3);
a2L=-0.8098 *10^(-5);
a3L=0.5693*10^(-8);
a4L=-0.15256 *10^(-11);
%COMPONENTE 2-CLORO PROPENO
a0R=3.5571;
a1R=5.3211*10^(-3);
a2R=4.8087 *10^(-5);
a3R=2.2693*10^(-8);
a4R=-0.11256 *10^(-11);
%COMPONENTE DICLORO PROPENO
a0D=3.5780;
a1D=10.9208*10^(-3);
a2D=-0.9098*10^(-5);
a3D=8.5693*10^(-8);
a4D=13.45256*10^(-11);
%FORMULA Cp=a0+a1*T+a2*T^2+a3*T^3 (cal/molK)
%COMPONENTE CLORURO DE ALILO
a0A=0.604;
a1A=4.617*10^(-2);
a2A=-7.942*10^(-5);
a3A=17.42*10^(-10);
%CALCULANDO LA ENTALPIA DE REACCION A UNA TEMPERATURA DE
510°C
%H°rx,808°K=H°rx,298°K+[Link]
%--Cp del acido clorhidrico (HCl)
syms T
f=(a0C+a1C*T+a2C*(T^(2))+a3C*(T^(3))+a4C*(T^(4)))*8.314;
a=298;%°K
b=783;%°K
CpCl=double(int(f,T,a,b));
%--Cp del propileno (C3H6)
syms T
f1=(a0P+a1P*T+a2P*(T^(2))+a3P*(T^(3))+a4P*(T^(4)))*8.314;
a=298;%°K
b=783;%°K
CpP=double(int(f1,T,a,b));
%--Cp del cloro (Cl2)
syms T
f2=(a0L+a1L*T+a2L*(T^(2))+a3L*(T^(3))+a4L*(T^(4)))*8.314;
a=298;%°K
b=783;%°K
CpL=double(int(f2,T,a,b));
%--Cp del cloruro de alilo (CH2=CHCH2Cl)
syms T
f3=(a0A+a1A*T+a2A*(T^(2))+a3A*(T^(3)))/0.238846;
a=298;%°K
b=783;%°K
CpA=double(int(f3,T,a,b));
%--Cp del 2-cloro propeno (CH2=CHCH2Cl)
syms T
f4=(a0R+a1R*T+a2R*(T^(2))+a3R*(T^(3))+a4R*(T^(4)))*8.314;
a=298;%°K
b=783;%°K
CpR=double(int(f4,T,a,b));
%--Cp del diclopropileno (C3H4Cl2)
syms T
f5=(a0D+a1D*T+a2D*(T^(2))+a3D*(T^(3))+a4D*(T^(4)))*8.314;
a=298;%°K
b=783;%°K
CpD=double(int(f5,T,a,b));
disp('********************************ENTALPIA DE REACCION A
510°C EN KJ/mol********************************')
disp('(1) REACCION PRINCIPAL C3H6+Cl2->CH2=CHCH2Cl+HCl')
HR1=Hf1+(((CpCl+CpA)-(CpP+CpL))/1000);
disp(HR1)
disp('(2) REACCION SECUNDARIA C3H6+Cl2->CH2=CClCH3+HCl')
HR2=Hf2+(((CpCl+CpA)-(CpR+CpL))/1000);
disp(HR2)
disp('(3) REACCION SECUNDARIA C3H6+2Cl2->C3H4Cl2+2HCl')
HR3=Hf3+(((2*CpCl+CpA)-(CpD+2*CpL))/1000);
disp(HR3)
disp('**********************CALOR GENERADO POR LA REACCION
EN EL REACTOR A 510°C EN KJ/h********************')
QR=(HR1+HR2+HR3)*(FT)*3600*1000*(1-0.958);
disp(QR)
if QR<0
disp('LAS REACCIONES QUE SE PRODUCEN EN EL REACTOR SON
EXOTERMICAS')
else
if QR>0
disp('LAS REACCIONES QUE SE PRODUCEN EN EL REACTOR
SON ENDOTERMICAS')
end
end
