INGENIERIA DE LAS REACCIONES QUIMICAS II

REACTOR CATALÍTICO DE LECHO FIJO

EVALUACION

Para la producción de anhídrido maléico se alimenta a un FIXED BED CATALITIC

REACTOR, 2 kmol/s de una mezcla que está formado de 1.1 % de benceno y 98.9 % de
aire. El proceso es adiabático y la reacción es:
C 6 H 6 + 4 O2 → C 4 H 2 O3 +2 H 2 O+CO2 +CO
A+ 4 B→ C+ 2 D+ E+ F
Los reactantes se precalientan inicialmente hasta la temperatura de 352°C.
Se asume que la reacción es de seudo primer orden, considerando un exceso de oxígeno,
donde
kmol de A
r p =k p A p B ; [
kg cat . · h ]
kmol
k=
[ kg cat .· h · atm2 ]
13636
lnk =28− ; T =[ K ]
T

DATOS ADICIONALES:
Diámetro del reactor, 2.54 cm
Densidad global del catalizador, 1300 kg/m3
Calor de reacción, - 1 285 409 kJ/kmol
Presión total, 1 atm
Capacidad calorífica, 0.992 kJ/kg·K
Formular la solución numérica para determinar la temperatura y profundidad del lecho
catalítico en función de la conversión

A. BALANCE DE ENERGÍA
Realizamos un balance de energía para un FBCR
 dT −F Ao ∆ H R
= se tienen dos componentes:
d xA F t Cp
Aire y benceno

kmol kmol benceno 3600 s kmol benceno

2 ∗( 0.011 )=0.022 x =79.2
S s h h

kmol kmol aire 3600 s kmol aire

2 ∗( 0.989 ) =1.978 x =7120.8
S s h h

kmol aire kmol N 2

 1.978 ∗( 0.79 )=1.5626 ( f nitrogeno )
s s

kmol aire kmol O 2

 1.978 ∗( 0.21 )=0.4154 ( f oxigeno )
s s

Relacionamos
La corriente total con la corriente de cada compuesto
fi
y i=
Ft
Ft = 2 kmol/s

Component Fracción molar Pesos moleculares

e
Benceno 0.011 72
Nitrogeno 0.7813 28
Oxígeno 0.2077 32

PESO MOLECULAR PROMEDIO

kg
PM Promedio =(72(0.011)+ 32(0.2077)+ 28(0.7813))=29.32
kmol

kmol C 6 H 6
 −F Ao =¿ 0.022
s

kmol total
 −F t=2
s
 CAPACIDAD CALORÍFICA
kJ
Cp=0.992
kgK

dT −F Ao ∆ H R
=
d xA F t Cp

kJ kmol
−(−1285 409) (0.022)
dT kmol s
=
d xA kmol kJ kg
2 (0.992 )(29.32 )
s KgK kmol

B. BALANCE DE MATERIA

dW F Ao
=
d x A −r p

dz F Ao
=
d x A ρAc(−r ¿¿ p) ¿
π
Ac= D 2
4
π
Ac= (2.54 cm)2
4
Ac=5.0671cm2

Calculo de k
Del problema tenemos

13636
lnk =28−
T
 −13636
elnk=e28∗e T

−13636
k =1.446∗1012∗e T
Cálculo de r p :

kmol de A
r p =k p A p B ; [ kg cat . · h ]
r p =k P A PB
Calculo de la presión parcial
pi
y i=
P
Donde la presión total es 1 atm
y i= pi

ESTEQUIMETRIA

C 6 H 6 + 4 O 2 → C 4 H 2 O 3 +2 H 2 O+CO2 +CO
A+ 4 B→ C+ 2 D+ E+ F

COMPONENTES INICIA CONVERSIÓN FINAL Fracción

L
C6 H6 1.1 -1.1 x 1.1 -1.1 x 1.1−1.1 x
100
O2 20.77 -4(1.1)x 20.77 – 4.4x 20.77 – 4.4 x
100
N2 78.13 0 78.13 78.13
100
C 4 H 2 O3 0 1.1x 1.1 x 1.1 x
100
H2O 0 2(1.1x) 2.2x 2.2 x
100
CO 2 0 1.1x 1.1 x 1.1 x
100
CO 0 1.1 1.1 x 1.1 x
100
TOTAL 100 1

dz F Ao
=
d x A ρAc(−r ¿¿ p) ¿
 −13636
k =1.446∗1012∗e T
kg
ρ=1300 3
m
Ac=5.0671 x 10−4 m 2

ya tenemos los datos calculados

r p =k p A p B

98.9−4.4 x 1.1∗( 1−x )

r p =k ( (
100
∗ )(
100
) )
Vamos a reemplazar en la ecuación general

kmol benceno 3600 s

0.022 x
dZ s h
=
d xA (28−13636
1300 kg T ) 98.9−4.4 x 1.1∗( 1−x ) kmol de A
m (
3
∗5.0671(10−4
) m 2
∗e ) ∗
100( ∗ )(
100 ) kg∗h

ECUACIONES FINALES
TÉRMINO INDEPENDIENTE

(xA)

TÉRMINOS DEPENDIENTES

T(K)

z (m)

kJ kmol
−(−1285 409)
(0.022)
dT kmol s
 d x = kmol kJ kg
A
2 (0.992 )(29.32 )
s KgK kmol
 
kmol benceno 3600 s
0.022 x
dZ s h
=
d xA (28−13636
1300 kg T ) 98.9−4.4 x 1.1∗( 1−x ) kmol de A
( m 3
−4 2
)
∗5.0671(10 ) m ∗e ∗ ( 100 )(
∗
100 ) kg∗h

clc, clear all;

disp('======================================================')
disp(' REACTOR CATALITICO DE CAMA FIJA (FBCR) ')
disp('======================================================')
disp('**********************************************************')
disp('**********************************************************')
disp('* CONVERSIÓN : xA')
disp('* LONGITUD : Z (m)')
disp('* TEMPERATURA: T(K)')
% Desarrollamos el sistema EDO
%datos del problema
%P=1;
%d=1300
%D=0.0254m
%Ac=(pi)*((D/2)^2)
%FA0=0.022;
%H=-1285409
%Ft=2
%Cp=0.992;
disp('=============================')
disp('C6H6+4O2? C4H2O3 +2H2O+CO2+CO')
disp('=============================')
P=1;d=1300;D=0.0254;Ac=(pi)*((D/2)^2); FA0=0.022;H=-
1285409;Ft=2;Cp=0.992;
dT=@(xA,T,z) ((0.022*3600)/(1300*(5.067*10^(-4))*(((exp(28-
13636/z))*((1.1-1.1*xA)/100)*((20.77-4.4*xA)/100)))));
dz=@(xA,T,z) (-H*FA0)/(Ft*Cp*29.35);
xA=[0,0];
T=[0;0];
z=[625;0];
h=0.01;
disp('*************************')
disp(' Cuadro de resultados ')
disp('*************************')
disp('==================================')
disp(' T(°K) z(m) xA ')
disp('==================================')
for i=1:85
k1=dT(xA(i),T(i),z(i));
c1=dz(xA(i),T(i),z(i));
k2=dT(xA(i)+h/2,T(i)+h*k1/2,z(i)+h*c1/2);
c2=dz(xA(i)+h/2,T(i)+h*k1/2,z(i)+h*c1/2);
k3=dT(xA(i)+h/2,T(i)+h*k2/2,z(i)+h*c2/2);
c3=dz(xA(i)+h/2,T(i)+h*k2/2,z(i)+h*c2/2);
k4=dT(xA(i)+h,T(i)+h*k3,z(i)+h*c3);
c4=dz(xA(i)+h,T(i)+h*k3,z(i)+h*c3);
T(i+1)=double(T(i)+(1/6)*h*(k1+2*k2+2*k3+k4));
z(i+1)=double(z(i)+(1/6)*h*(c1+2*c2+2*c3+c4));
xA(i+1)=xA(i)+h;
fprintf(' %8.4f %8.4f %8.2f\n',z(i),T(i),xA(i));
end
plotyy(T,xA,T,z),grid on
[hAx,hLine1,hLine2] = plotyy(T,xA,T,z);
ylabel(hAx(1),'CONVERSION xA'),ylabel(hAx(2),'TEMPERATURA(K)')
xlabel('z (m)')
grid on

======================================================
REACTOR CATALITICO DE CAMA FIJA (FBCR)
======================================================
**********************************************************
**********************************************************
* CONVERSIÓN : xA
* LONGITUD : Z (m)
* TEMPERATURA: T(K)
=============================
C6H6+4O2? C4H2O3 +2H2O+CO2+CO
=============================
*************************
Cuadro de resultados
*************************
==================================
T(°K) z(m) xA
==================================
625.0000 0.0000 0.00
629.8564 1.0062 0.01
634.7128 1.8681 0.02
639.5692 2.6084 0.03
644.4256 3.2460 0.04
649.2820 3.7964 0.05
654.1384 4.2728 0.06
658.9948 4.6861 0.07
663.8512 5.0455 0.08
668.7076 5.3589 0.09
673.5640 5.6327 0.10
678.4204 5.8724 0.11
683.2768 6.0828 0.12
688.1332 6.2679 0.13
692.9896 6.4311 0.14
697.8460 6.5752 0.15
702.7024 6.7027 0.16
707.5588 6.8158 0.17
712.4152 6.9163 0.18
717.2716 7.0058 0.19
722.1280 7.0857 0.20
726.9844 7.1570 0.21
731.8407 7.2209 0.22
736.6971 7.2783 0.23
741.5535 7.3298 0.24
746.4099 7.3761 0.25
751.2663 7.4180 0.26
756.1227 7.4558 0.27
760.9791 7.4899 0.28
765.8355 7.5209 0.29
770.6919 7.5491 0.30
775.5483 7.5747 0.31
780.4047 7.5980 0.32
785.2611 7.6193 0.33
790.1175 7.6387 0.34
794.9739 7.6565 0.35
799.8303 7.6728 0.36
804.6867 7.6877 0.37
809.5431 7.7015 0.38
814.3995 7.7141 0.39
819.2559 7.7258 0.40
824.1123 7.7366 0.41
828.9687 7.7466 0.42
833.8251 7.7558 0.43
838.6815 7.7644 0.44
843.5379 7.7723 0.45
848.3943 7.7797 0.46
853.2507 7.7866 0.47
858.1071 7.7930 0.48
862.9635 7.7990 0.49
867.8199 7.8047 0.50
872.6763 7.8099 0.51
877.5327 7.8148 0.52
882.3891 7.8195 0.53
887.2455 7.8238 0.54
892.1019 7.8279 0.55
896.9583 7.8318 0.56
901.8147 7.8354 0.57
906.6711 7.8389 0.58
911.5275 7.8422 0.59
916.3839 7.8453 0.60
921.2403 7.8482 0.61
926.0967 7.8510 0.62
930.9531 7.8537 0.63
935.8095 7.8562 0.64
940.6658 7.8587 0.65
945.5222 7.8610 0.66
950.3786 7.8632 0.67
955.2350 7.8654 0.68
960.0914 7.8674 0.69
964.9478 7.8694 0.70
969.8042 7.8713 0.71
974.6606 7.8732 0.72
979.5170 7.8750 0.73
984.3734 7.8767 0.74
989.2298 7.8784 0.75
994.0862 7.8801 0.76
998.9426 7.8817 0.77
1003.7990 7.8833 0.78
1008.6554 7.8849 0.79
1013.5118 7.8864 0.80
1018.3682 7.8879 0.81
1023.2246 7.8894 0.82
1028.0810 7.8909 0.83
1032.9374 7.8924 0.84