Nombre: Guardado Rodrı́guez Valeria 71-x

Problemas de Reactores.

6.21-Tiscareño. Un reactor de tanque agitado de 10m3 con calentamiento

se alimenta con una solución acuosa a 28◦ C que contiene 2.5 M de A. Las dos
reacciones relevantes son:

2A → B → P roductos

y las velocidades de reacción, por mol de reacción, están dadas por:

r1 = k1 CA (1)
r2 = k2 CB (2)
donde k1 = 0,21min−1 y k2 = 0,05min−1 a 72◦ . temperatura a la que debe
salir la corriente del reactor. Los calores de reacción para las reacciones (1) y
Kcal Kcal
(2) son 25 moldeA y 30 moldeB , respectivamente. Para el calentamiento se dispone
de vapor saturado a 112◦ con un calor latente de vaporización de 530 cal g . Si se
especifica un rendimiento de B de 0.5, que equivale a CB1 = 0,625M , determine:
a)CA1 y fA1
m3
b)El flujo volumétrico de la alimentación en min ;y
Kg
c)El flujo de vapor alimentado en min
Solución: a)Primero obtenemos el balance de masa, como son dos ecuacio-
nes eso implica que debemos obtener dos balances de masa. Escogemos como
elementos independientes a A y B:
Con respecto a A:
CA0 − CA
τ= (3)
−rA
Con respecto a B:
CB0 − CB
τ= (4)
−rB
Ahora, obtenemos rA y rB de las ecuaciones de r1 y r2 . Primero debemos
observar en cuales reacciones participan cada uno. Y de eso partimos para ob-
tener la velocidad de reacción por mol de compuesto. Donde se multiplica la
velocidad de reacción por el coeficiente estequiometrico del reactivo.

rA = (υA )(r1 ) = (−2)(k1 CA ) = −2k1 CA (5)

rB = (υB )(r1 ) + (υB )(r2 ) = (1)(k1 CA ) + (−1)](k2 CB ) = k1 CA − k2 CB (6)

Después, sustituimos en las ecuaciones (3) y (4):

CA0 − CA CA0 − CA
τ= = (7)
−(−2k1 CA ) 2k1 CA

CB0 − CB C B0 − C B
τ= = (8)
−(k1 CA − k2 CB ) −k1 CA + k2 CB

1
 Ahora podemos observar que en ambas ecuaciones tenemos dos incognitas y
las podemos resolver. Ası́ que ahora despejaremos las ecuaciones para igualarlas
a cero. Con respecto a A:

τ (2k1 CA ) = CA0 − CA

2τ k1 CA − CA0 + CA = 0 (9)
Con respecto a B:

τ (−k1 CA + k2 CB ) = CB0 − CB

−τ k1 CA + τ k2 CB − CB0 + CB = 0 (10)
Conocemos el valor de k1 , k2 , CA0 , CB0 y CB , entonces las podemos resolver
para obtener el valor de CA y τ . Resolvemos con Newton-Raphson para dos
ecuaciones.
Derivamos primero las ecuaciones con respecto a τ y a CA :

δ(2τ k1 CA − CA0 + CA )
= 2k1 CA (11)
δτ
δ(2τ k1 CA − CA0 + CA )
= 2τ k1 + 1 (12)
δCA
Haciendo lo mismo con la ecuación (10):

δ(−τ k1 CA + τ k2 CB − CB0 + CB )
= −k1 CA + k2 CB (13)
δτ
δ(−τ k1 CA + τ k2 CB − CB0 + CB )
= −τ k1 (14)
δCA
Abriendo el programa de Newton-Raphson en MATLAB, que es el que se
muestra a continuación:

2
Corriendo el programa y dandole los valores iniciales correspondientes:

Donde x es τ y y CA que nos da un valor de:

τ = 3,333 y CA = 1,041.
Ahora vamos a sacar el valor de la fracción conversión:
CA0 − CA 2,5 − 1,041
fa = = = 0,5836
CA0 2,5

3
 b)
Al resolver los balances de masa de manera simultanea pudimos obtener el
valor del tiempo espacial, y a partir de este podemos obtener el flujo volumetrico:
V
τ=
v̇
V
v̇ =
τ
Sustituyendo el volumen por 10m3 y τ por 3.333 min obtenemos:

10 m3
v̇ = = 3,00
3,333 min

c) A partir del balance de energı́a, podemos obtener el calor y de ahı́ obtener

el flujo de vapor en la entrada.
X
0 = Q − v̇ρC¯p (T − Tref ) − ∆Hr VR rr (15)

Despejamos el calor (Q) y lo calculamos, suponiendo la ρ = 1000 Kg ¯

m3 y Cp =
cal
1000 Kg◦ C como las del agua:

Q = v̇ρC¯p (T − Tref ) + (∆Hr1 r1 VR + ∆Hr2 r2 VR )

m3 kg cal cal
Q = (3 )(1000 3 )(1000 ◦ )(72−28)◦ C+((∗25000 )(0,2186M )(10000Lt)
min m kg C moldeA
cal
+(30000 )(0,0312M )(10000lt))
moldeB
Obteniendo un valor de, Q = 67990 kcal
min . Ahora tenemos que el calor es igual
a:
Q = Lv m
Despejamos la masa:
m = Q/Lv
Y la calculamos:
67990 kcal
min g kg
m= = 128283,01 = 128,28
0,530 kcal
g
min min

4
 6.23-Tiscareño. Se desea analizar el efecto de utilizar uno u otro de los
inertes I1 e I2 en la producción de C en fase gaseosa. El proceso dispone de un
reactor tubular de 22 m3 y el sistema de reacción es:
−15,300 cal
mol
0,5 cal 0,5
2A + 12 B → C r1 = 1,9x107 segltmol0,5 exp 1,987 molK T CA CB
−18,700 cal
mol
0,5 cal 0,5
C + 12 B → D r2 = 8,5x107 segltmol0,5 exp 1,987 molK T CC CB

Los calores de reacción a 200◦ C son -52.4 y -50.9 kcal/mol de reacción, res-
pectivamente; éstos no deben considerarse constantes. La alimentación consiste
en 550 lt/seg a 1.5 atm y 200◦ C. Su composición molar es 1 % de A, 1 % de B y
98 % de inerte. El reactor es adiabático y se pueden despreciar las caidas de pre-
cal
sión. Se cuenta con las siguientes capacidades calorificas en mol◦C , y tampoco
deben considerarse constantes:

CP A = 10,4 + 0,020T CP D = 8,2 + 0,010T

CP B = 7,2 + 0,002T CP I1 = 7,6 + 0,014T
CP C = 12,8 + 0,024T CP I2 = 18,6 + 0,014T

donde la temperatura debe estar en ◦ C.

Utilice 200 pasos de integración con el método básico de Runge-Kutta de
cuarto orden y determine fA , SA→C , RA→C y T en la salida si:
a)Se emplea el inerte I1 ; y
b)Se emplea el inerte I2 .
c)Explique sus resultados.
a) Comenzamos con los balances de masa, A y B como nuestros elementos
independientes:
dFA
= rA (16)
v̇dt
dFB
= rB (17)
v̇dt
Tenemos los balances del reactor:

FA = FA0 − 1

FB = FB0 − 0,51 − 0,52

FC = FC0 + 1 − 2
FD = FD0 + 2
Como tenemos A y B como elementos independientes, debemos dejar las
ecuaciones de C y D en función de estas. De la ecuación de A despejamos 1 y
de la ecuación de B, despejamos 2 .

1 = FA0 − FA
FB0 − FB − 0,51 FB0 − FB − 0,5(FA0 − FA )
2 = =
0,5 0,5
Sustituimos en las ecuaciondes de C y D:
FB0 − FB − 0,5(FA0 − FA )
FC = FC0 + FA0 − FA − ( )
0,5

5
 También necesitamos las ecuaciones de rA y rB :
0,5
rA = −k1 CA CB
0 0,5
rB = −0,5k1 CA CB ,5 − 0,5k2 CC CB
Sustituimos en las ecuaciones (16) y (17):

dFA 0,5
= v̇(−k1 CA CB ) (18)
dt
dFB 0 0,5
= v̇(−0,5k1 CA CB ,5 − 0,5k2 CC CB ) (19)
dt
También tenemos que dejar las concentraciones en función del flujo molar de
cada uno, donde tenemos que:
FA
CA =
v̇
FB
CB =
v̇
FC
CC =
v̇
Las sustituimos en las ecuaciones (18) y (19).

dFA FA FB0,5
= v̇(−0,5k1 ) (20)
dt v̇ v̇ 0,5

dFB FA FB0,5 FC FB0,5

= v̇(−0,5k1 0,5
− 0,5k2 ) (21)
dt v̇ v̇ v̇ v̇ 0,5
Sustituyendo el valor de FC por el balance en el reactor ya dependiendo de
los elementos independientes:
F −FB −0,5(FA0 −FA )
dFB FA FB0,5 FC0 + FA0 − FA − ( B0 0,5 ) FB0,5
= v̇(−0,5k1 −0,5k2 )
dt v̇ v̇ 0,5 v̇ v̇ 0,5
(22)
Sustituyendo los valores de k1 y k2 :
0,5
dFA −15,300 FA F
= v̇(−0,5(1,9x107 exp 1,987T ) B
) (23)
dt v̇ v̇ 0,5
−15,300 0,5
dFB FB
dt = v̇(−0,5(1,9x107 exp 1,987T ) Fv̇A v̇ 0,5

FB0 −FB −0,5(FA0 −FA )

7
−18,700 FC0 + FA0 − FA − ( 0,5 ) FB0,5
−0,5(8,5x10 exp 1,987T ) ) (24)
v̇ v̇ 0,5