Ejercicio Nº 4:

Se va ha hacer una oxidación catalítica en fase gaseosa, entre una sustancia A y B (oxígeno) en un reactor tubular
empacado con catalizador, operará isotérmicamente, la temperatura de operación es 260ºC. Las sustancias A y B
se alimentan al reactor en proporciones estequiométricas. La sustancia A se alimenta a razón de 0.30 lbmol/s y
presión de 10 atm. Se propone usar 10 bancos de tubería de diámetro de 11/2 pulgada y cédula 40, rellena con
catalizador, se tienen 100 tubos por banco. En consecuencia, el flujo molar en cada tubo será de 3 x10-4 lbmol/s.
Las propiedades del fluido reaccionante puede considerarse idénticas a las del aire a éstas temperatura y presión.
Considerar que el oxígeno que se alimenta proviene del aire, trayendo nitrógeno como inerte. La densidad de las
partículas de catalizador de 1/4 de pulgada es de 120lb/ft3 y la fracción de vacío en el lecho es de 0.45. La ley de
velocidad es:

La reacción es : A + 1/2 B --------> C

1 2
3 3
-rA = k PA  PB

lbmol
k = 0.0141
atm  lb hr

a) Grafique conversión y presión en función de peso de catalizador.

b) Averiguar el peso de catalizador necesario para lograr una conversión del 60%. ¿Cuál será la caída de presión en
el lecho al alcanzar esta conversión?
c) Calcule la cantidad de catalizador necesaria para alcanzar la misma conversión, pero sin tener en cuenta la caída
de presión.
En este caso ¿es importante tenerla en cuenta? ¿Por qué?

Para el cálculo de la caída de presión en un lecho empacado se usará la ecuación de Ergun (ver
pág.177-H.Scott Fogler)

dy -α T
=  ( 1 + ε  X) 
dW 2 y To

dy -α
Para operación isotérmica tenemos: =  ( 1 + ε  X)
dW 2 y

2  βo
α=
Donde : Ac ρc ( 1 - ϕ)  Po

bo :constante que depende de las propiedades del lecho empacado y las condiciones de entrada.
f : porosidad = volumen de vacío / volumen global del lecho = fracción de vacío.
r: densidad del catalizador sólido
Ac: área de la sección transversal de

-5 1
α := 3.656 10 
gm
DATOS
J
lbmol := 453.5924mol kJ := 1000J Rg := 8.3144
mol K
P0 := 10atm T0 := 260 °C

0.3 lbmol 0.15 lbmol 79

FA0 := FB0 := FI0 := FB0 FT0 := FB0 + FA0 + FI0
1000 s 1000 s 21

FB0 FI0 FA0 P0 FT0 3

-3m
yB0 := yI0 := yA0 := CT0 := FV0 := = 2.013  10
FT0 FT0 FT0 Rg T0 CT0 s

lb 1
ρc := 120 ϕ := 0.45 Dp :=  in
3 4
ft

A + 1/2 B --------> C

Especie Simbolo Inicial Cambio Final

A FA0 -FA0 X FA = FA0 ( 1 - X)

-b
FB = FA0  ΘB -  X
b
B FB0 = FA0 ΘB  FA0 X
a  a 

FC = FA0  ΘC +  X
c c
C FC0 = FA0 ΘC  FA0 X
a  a 

FI = FA0 ΘI
Inerte I FI0 = FA0 ΘI
Coeficientes estequimetricos Valores de Θ

a := 1 c := 1 yA0 yB0
ΘA := =1 ΘC := 0 ΘB := = 0.5
yA0 yA0
1
b := yI0
2 ΘI := = 1.881
yA0

c b a
δ := - - = -0.5 ε := δ yA0 = -0.148
a a a

CONCENTRACIONES

FA FA0 ( 1 - X) CA0 ( 1 - X) CA0 ( 1 - X) P

CA = = = = 
FV  P0  T  P0 ( 1 + ε  X) P0
FV0       ( 1 + ε  X)  ( 1 + ε  X)
 P  T0   P

C A0  ΘB -  X CA0  ΘB -  X

b b
FB
CB = =
 a  =  a  P
FV P0 ( 1 + ε  X) P0
 ( 1 + ε  X)
P

CA0 ( 1 - X) P PA0
PA = CA Rg T =   Rg T CA0 =
( 1 + ε  X) P0 Rg T0

( 1 - X) P PA0
PA =   Rg T Si ademas se considera isotermico
( 1 + ε  X) P0 Rg T0

( 1 - X)  yA0
PA = P
( 1 + ε  X)

 Θ - b  X  y
 B 
 a  A0
PB = P
( 1 + ε  X)

Ley de Velocidad

1 2
3 3 dX dX -rA
-rA = k PA  PB FA0 = -rA =
gmol A /g catalizador hr dW dW FA0
 lbmol 1
k := 0.0141 α = 0.037
atm  lb hr kg
2
dP -α P0
=   ( 1 + ε  X)
dW 2 P

2
1 3
  Θ - b  X  y 
 ( 1 - X)  yA0    B a  A0 
3
-rA = k   P    P
 ( 1 + ε X)   ( 1 + ε  X) 
2
1
 yA0  3
 P  ( 1 - X)   ΘB -  X
3 b
-rA = k 
 ( 1 + ε  X)   a 

2
1
yA0 3
 P ( 1 - X)   ΘB -  X
dX k 3 b
= 
dW FA0 ( 1 + ε  X)  a 

2
yA0 3
 P ( 1 - X)   
dX k 1
= 
dW FA0 ( 1 + ε  X) 2

2
( 1 - X)  yA0 3
 P  
dX k 1
= 
dW FA0 ( 1 + ε  X) 2
2
dP -α P0
=   ( 1 + ε  X)
dW 2 P

0 X
M :=   M= 
 10  P
  k 2
 

mol
atm s kg (1 - M0) yA0  M   1  3 
   

F A0 (1 + ε M0) 1  2  
D( W , M ) :=  
mol

 s 
 -α 2 
  P0  
1  
  atm  
 (1 + ε M )
kg
  
2 M 0
 1 

Vector := rkfixed( M , 0 , 27.5 , 55 , D)

0 1 2 0
Wvec := Vector [=] kg
0 0 0 10
1 0.5 0.02639 9.9083594 1
Xvec := Vector
2 1 0.0519465 9.8162209
3 1.5 0.0766935 9.723559 2
Pvec := Vector [=]atm
4 2 0.1006542 9.6303473
5 2.5 0.1238514 9.5365587
6 3 0.1463074 9.4421653
Vector = 7 3.5 0.1680438 9.3471382
8 4 0.1890817 9.2514474
9 4.5 0.2094417 9.1550621
10 5 0.2291439 9.0579501
11 5.5 0.2482077 8.960078
12 6 0.2666522 8.861411
13 6.5 0.2844957 8.7619128
14 7 0.3017561 8.6615456
15 7.5 0.3184508 ...

10

8 0.8

6 0.6
Pvec = 5.24
Pvec Xvec 42

4 0.4

2 0.2

0 0
0 10 20 30
Wvec
b) Averiguar el peso de catalizador necesario para lograr una conversión del 60%. ¿Cuál
será la caída de presión en el lecho al alcanzar esta conversión?

De acuerdo a los resultados, para una conversion del 60%

 
W60 := Wvec -  Xvec - 0.6 
Wvec - Wvec ( ) kg
( )
42 41
 42  42 Xvec - Xvec 
  42 41 

W60 = 20.559 kg Masa de catalizador

 Xvec - 0.6  (
  Pvec - Pvec ) 
P60 := Pvec - ) Xvec - Xvec  atm
(
42 41
 42 42
  42 41 

P60 = 5.379 atm

ΔP := P0 - P60 = 4.621 atm Caida de Presión

c) Calcule la cantidad de catalizador necesaria para alcanzar la misma conversión, pero sin
tener en cuenta la caída de presión.

2
( 1 - X)  yA0 3
 P0  
dX k 1
=  
dW FA0 ( 1 + ε  X) 2

dX
dW =
 2
 
3
 k ( 1 - X)  yA0   
1
 F  ( 1 + ε  X)  P0   2  
 A0 

0.6
 1
Wc :=  dX
  2
  
3
  k ( 1 - X)  yA0   
1
  F  ( 1 + ε  X)  P0   2  

  A0 
0
Wc = 16.213 kg Sin considerar la caida de presión

En este caso ¿es importante tenerla en cuenta? ¿Por qué?

ΔW := W60 - Wc = 4.345 kg
En este caso si es importante consideara la caida de presión
De lo contrario se cometeria un error del 21.1% y la conversion
ΔW seria del 53.4%
= 0.211
W60