SOLUCIONARIO DEL SEGUNDO EXAMEN DE

REACCIONES QUÍMICAS II

1. 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 y graficar .

SOLUCIÓN:
 352 ° C+ 273.15 K =625.15 K
 kmol kmol benceno
2 ∗( 0.011 )=0.022
S s
kmol kmol aire
2 ∗( 0.989 ) =1.978
S s

kmol aire kmol O 2

1.978 ∗( 0.21 )=0.41538
s s

kmol aire kmol N 2

1.978 ∗( 0.79 )=1.5626
s s
Productos ingresantes
C 6 H 6 =0.022; PM =72

O 2=0.4153 ; PM =32

N 2=1.5626 ; PM =28

Calculo del PM Promedio

PM Promedio =PM Ben∗y Ben+ PM O 2∗ y O 2+ PM N 2∗y N 2 …(1)

Calculo dela composición molar de cada compuesto ingresante
fi
y i= …(2)
Ft
0.022
y Ben = =0.011
2
0.41538
y O 2= =0.20769
2
1.56262
y N 2= =0.78131
2
Reemplazando en (3):
kg
PM Promedio =72∗0.011+32∗0.20769+ 28∗0.78131=29.38
kmol
1. BALANCE DE ENERGÍA

dT −∆ H R∗F A
= …(3) o

dx FT ∗C p

Del enunciado tenemos los siguientes datos:

∆ H R=−1285 409 kJ / kmol

Cp=0.992kJ /kg·K
De cálculo:
kg
PM Promedio =29.38
kmol
Reemplazando en (3):

kJ kmol
−−1285409 ∗0.022
dT kmol s
= =485.14 K …(4)
dx kmol kJ kg
2 ∗0.992 ∗29.38
s kg∗K kmol

Ordenando e integrando:

dT =485.14∗dx
T xA

∫ dT =¿ 485.14 ∫ dx ¿
T0 0

T =485.14 x+T o

2. BALANCE DE MATERIA

d ZA −F A
= o
… (6 )
d x A ρC ∆ T ∗AC ∗r p

De enunciado:

F A =2 kmol /s
o

ρC ∆ T =1300 kg /m3

Cálculo del área transversal ( AC ):

Ac=π∗r 2
 D 2
Ac=π∗( )
2

π∗D 2
Ac =
4
π
Ac= ∗(2.54∗10−2 )2
4

2.54
=0.0254 m
100

π∗0.02542
Ac =
4

Ac=5.06707∗10−4 m 2

Reemplazando en (6):
kmol
∗3600 s
s
0.022
dz 1h
=
dx kg kmol
1300
m3 [
∗5.067∗10−4 m2∗r P
kg cat∗h ]
dz 120.235
= …(7)
dx rp

Cálculo de r p :

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

Calculo de k

13636
lnk =28−
T
−13636
28 T
k =e ∗e
 −13636
k =1.446∗1012∗e T

Calculo de las fracciones en función de la conversión:

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

COMPONENTES COEFICIENTES INICIAL CAMBIO FINAL

C6 H6 1  0.022  −0.022∗x  0.02−0.022∗x
 N 2 No inter  FN 2  No inter  FN 2
 
−b
∗0.022∗x
 O 2 4  0.41538 a  0.41538−4∗(0.022)∗x
 C 4 H 2 O3 1   F C 0=0  0.022∗x  0.022∗x
 H 2O 2   F D 0 =0  2∗0.022∗x  2∗0.022∗x
 CO 2 1   F E 0 =0  0.022∗x  0.022∗x
 CO 1   F F 0=0  0.022∗x  0.022∗x
∑ ¿ 0.022+ F N 2+ 0.41538=F T 0

F A 0 −F A 0∗x 0.022−0.022∗x
y A= =
FT 0 2
F B 0−4∗F A 0∗x 0.41538−4∗(0.022)∗x
yB = =
FT 0 2

Reemplazando en (7):
dz 120.235
=
dx k∗p A∗p B

dz 120.235
=
dx 12
−13636
T 0.022−0.022∗x 0.41538−4∗(0.022)∗x
1.446∗10 ∗e ∗( )∗( )
2 2
Reemplazando el vapor de r p
dZ 2 kmol /s∗3600 s/h
=
dx (28−13636
T ) 1.1∗( 1−x )
( 1300 kg /m3∗1.0178∗10−3 m2 )∗e ∗ ( 98.9−4.4
100
x
)∗(
100 ) kmol de A
kg kg∗h

SOLUCIÓN NUMÉRICA POR MATLAB

clc, clear,clf
disp('------------METODO DE RUNGE-KUTTA----------------')
P=1;
d=1300;
D=2.54/100;
Ac=(pi)*((D/2)^2);
F0=0.022;
AHR=-1285409;
Ft=2;
Cp=0.992;
dy=@(x,y,z) 0.022/(1300*(5.067*10^(-4))*(((exp(28-13636/z))*((1.1-
1.1*x)/100)*((20.77-4.4*x)/100))*(1/3600)));
dz=@(x,y,z) (-AHR*F0)/(Ft*Cp*29.35);
x=[0,0];
y=[0;0];
z=[625;0];
h=0.01;
disp('_______________________________')
disp(' x Z(ft) T(K)\n')
disp('_______________________________')
for i=1:80
k1=dy(x(i),y(i),z(i));
c1=dz(x(i),y(i),z(i));
k2=dy(x(i)+h/2,y(i)+h*k1/2,z(i)+h*c1/2);
c2=dz(x(i)+h/2,y(i)+h*k1/2,z(i)+h*c1/2);
k3=dy(x(i)+h/2,y(i)+h*k2/2,z(i)+h*c2/2);
c3=dz(x(i)+h/2,y(i)+h*k2/2,z(i)+h*c2/2);
k4=dy(x(i)+h,y(i)+h*k3,z(i)+h*c3);
c4=dz(x(i)+h,y(i)+h*k3,z(i)+h*c3);
y(i+1)=y(i)+h*(k1+2*k2+2*k3+k4)/6;
z(i+1)=z(i)+h*(c1+2*c2+2*c3+c4)/6;
x(i+1)=x(i)+h;
fprintf('%8.2f %8.4f %8.2f\n',x(i),y(i),z(i));
end
plotyy(y,x,y,z),title('Conversión, T°K y Z(ft)'),
xlabel('Z(ft)'),ylabel('conversion'),grid minor
text(2,0.15,'T°(k)');
text(6,0.075,'Conversión');

RESULTADO, CORRIENDO EL PROGRAMA

------------METODO DE RUNGE-KUTTA----------------

_______________________________

x Z(ft) T(K)\n

_______________________________

0.00 0.0000 625.00

0.01 1.0062 629.86

0.02 1.8681 634.71

0.03 2.6084 639.57

0.04 3.2460 644.43

0.05 3.7964 649.28

0.06 4.2728 654.14

0.07 4.6861 658.99

0.08 5.0455 663.85

0.09 5.3589 668.71

0.10 5.6327 673.56

0.11 5.8724 678.42

0.12 6.0828 683.28

0.13 6.2679 688.13

0.14 6.4311 692.99

0.15 6.5752 697.85

0.16 6.7027 702.70

0.17 6.8158 707.56

0.18 6.9163 712.42

0.19 7.0058 717.27

0.20 7.0857 722.13

0.21 7.1570 726.98

0.22 7.2209 731.84

0.23 7.2783 736.70

0.24 7.3298 741.55

0.25 7.3761 746.41

0.26 7.4180 751.27

0.27 7.4558 756.12

0.28 7.4899 760.98

0.29 7.5209 765.84

0.30 7.5491 770.69

0.31 7.5747 775.55

0.32 7.5980 780.40

0.33 7.6193 785.26

0.34 7.6387 790.12

0.35 7.6565 794.97

0.36 7.6728 799.83

0.37 7.6877 804.69

0.38 7.7015 809.54

0.39 7.7141 814.40

0.40 7.7258 819.26

0.41 7.7366 824.11

0.42 7.7466 828.97

0.43 7.7558 833.83

0.44 7.7644 838.68

0.45 7.7723 843.54

0.46 7.7797 848.39

0.47 7.7866 853.25

0.48 7.7930 858.11

0.49 7.7990 862.96

0.50 7.8047 867.82

0.51 7.8099 872.68

0.52 7.8148 877.53

0.53 7.8195 882.39

0.54 7.8238 887.25

0.55 7.8279 892.10

0.56 7.8318 896.96

0.57 7.8354 901.81

0.58 7.8389 906.67

0.59 7.8422 911.53

0.60 7.8453 916.38

0.61 7.8482 921.24

0.62 7.8510 926.10

0.63 7.8537 930.95

0.64 7.8562 935.81

0.65 7.8587 940.67

0.66 7.8610 945.52

0.67 7.8632 950.38

0.68 7.8654 955.24

0.69 7.8674 960.09

0.70 7.8694 964.95

0.71 7.8713 969.80

0.72 7.8732 974.66

0.73 7.8750 979.52

0.74 7.8767 984.37

0.75 7.8784 989.23

0.76 7.8801 994.09

0.77 7.8817 998.94

0.78 7.8833 1003.80

0.79 7.8849 1008.66

GRÁFICA RELACIÓN DE LA CONVERSIÓN DE T° Y Z:

2.Para la producción del ácido sulfúrico es importante calcular los efectos externos
en función de la variación de las concentraciones en la superficie del catalizador y
la corriente global.
Tomando en consideración la reacción catalítica:
SO2 + 0.5 O2  SO3
Donde el catalizador platino soportado sobre alúmina tiene la forma de un cilindro
de 1/8 x 1/8 de pulgada.
DATOS:
Temperatura = 458 °C
Presión = 790 mm Hg
Velocidad de masa superficial = 245 lb/hr-ft 2
Velocidad de alimentación de SO2 = 0.0475 mol-g/min
Velocidad de alimentación del aire = 0.681 mol-g/min
Velocidad de reacción = 0.094 mol-g SO3 producido/h-g de catalizador
Área externa superficial/masa de catalizador = 5.12 ft 2/lb
Viscosidad de la mezcla reaccionante = 0.032 cp
Densidad del pellet = 112.8 lb/ft3
Porosidad del lecho = 0.40
Difusividad del SO3 en la mezcla reaccionante = 0.484 cm2/seg
Determinar la variación de la concentración.

SOLUCIÓN:
SO2 + 0.5 O2  SO3

A=1.8 pulg∗1.8 pulg=0.15 ft

2/ 3
u
Cb−Cs=
rp( )ρD
(1)
G
Jd∗at ( )
g

5.12 ft 2
at = masa de catalizador
lb
 G=245 lb/hr−ft 2

0.094 lbmol
rp=0.094 molgr /gcath=
lbcat∗h
lb
ρ=112.8
ft 3
Hallamos Jd
Empezamos hallando los Re
dpG
ℜ=
u

u=0.032 cp ( 0.0011 cpkg/ms )=3.2∗10 −5

kg/ms ( 0.451lbkg )( 3.28
1m 3600 h
ft )( 1 h )
u=0.0774 lb/ft∗¿

Hallando dp

3
A= π D2 A=π∗dp 2 D=0.15 ft
2

dp=0.183 ft
Entonces:

0.183 ft∗240 lb/hft 2

ℜ=
lb
0.774
ft∗h
ℜ=567.44

Ya que Re> 10 Aplicamos

0.458 −0.407
Jd= ( ℜ) =0.0867
0.40
Hallando Sc

u=0.0774 lb /ft∗h
 lb
ρ=112.8
ft 3

cm2 0.03282 ft 2 3600 s

D=0.484
s ( 1 cm2 )( h )
ft 2
D=1.874
s

0.0774 lb/ft∗h
Sc= 2
=3.66∗10− 4
lb ft
112 3 ∗1.874
ft h

Sc=3.66∗10−4
lb 3/ 2
0.094 ∗( 3.66∗10− 4 )
lbcat∗h
Cb−Cs=
ft 2
∗240 lb/hft 2
lb
0.0867∗5.12
lb
112.8 3
ft
lbmol
Cb−Cs=5.1759∗10−4
ft 3