TALLER DE INGENIERIA DE REACCIONES QUIMICAS

1. Una reacción química irreversible en fase liquida A→B es llevada a cabo en un CSTR. Se han tomado
diferentes datos de tiempo de residencia para estudiar la ley de velocidad de reacción. A entra puro al
reactor a una concentración de 7.5 moles/litro en cada uno de los experimentos
Run 1 2 3 4
τ (min) 1 15 30 100
CA (mol/lit) 3.2 0.72 0.46 0.21
a. Escribir el balance de moles para un CSTR donde una reacción química de orden n toma lugar
(en términos del tiempo de residencia)
b. Obtener el orden de reacción y la constante de velocidad de
reacción. RESPUESTA
a) Para este tipo de ejercicios, balance general de moles significa encontrar una ecuación general con la
cual puedas encontrar el orden de la reacción. Se que ahora mismo no suena muy claro, pero en el
inciso b podrán entender mejor que carajo
PARTIMOS DESDE LA ECUACIÓN DE DISEÑO DE UN CSTR
𝐹𝑎0𝑥
𝑉𝐶𝑆𝑇𝑅 =
−𝑟𝑎
𝑉𝑐𝑠𝑡𝑟 = 𝐶𝑎0𝑣0𝑥
𝑘𝐶𝑎𝑛
𝑁𝑜𝑠 𝑝𝑖𝑑𝑒𝑛 𝑢𝑠𝑎𝑟 𝑡𝑖𝑒𝑚𝑝𝑜 𝑑𝑒 𝑟𝑒𝑠𝑖𝑑𝑒𝑛𝑐𝑖𝑎 (τ)
V
τ=
v0
𝑉𝑐𝑠𝑡𝑟 𝐶𝑎0𝑥
=
𝑣0 𝑘𝐶𝑎𝑛
τ = 𝐶𝑎0𝑛𝑥
𝑘𝐶 𝐶 𝑥
𝑎0
𝑘𝐶 = 𝑎
𝑛

𝑎
τ
A este punto ya se podría dejar así, todo lo que no conocemos (la k y el orden de reacción) están a la
izquierda, pero a la derecha tenemos la conversión x, la cual no la especifican en el ejercicio, por lo que
hay que cambiarla: Haciendo balance de masa (lo que entra - lo que reacciona = lo que sale)
𝐸𝑛𝑡𝑟𝑎𝑛 𝐶𝑎0, 𝑠𝑒 𝑐𝑜𝑛𝑣𝑖𝑒𝑟𝑡𝑒𝑛 𝐶𝑎0𝑥 𝑦 𝑠𝑎𝑙𝑒𝑛 𝐶𝑎
𝐶𝑎 = 𝐶𝑎0 − 𝐶𝑎0𝑥
𝐶𝑎 − 𝐶𝑎0 = −𝐶𝑎0𝑥
𝐶𝑎0 − 𝐶𝑎 = 𝐶𝑎0𝑥
𝑃𝑜𝑟 𝑙𝑜 𝑞𝑢𝑒 𝑎𝑙 𝑟𝑒𝑒𝑚𝑝𝑙𝑎𝑧𝑎𝑟:
𝐶𝑎0 − 𝐶𝑎
𝑘𝐶𝑎𝑛 =
τ
Esta es una ecuación con la que podemos encontrar la k y el orden de reacción n, pues todo lo que está
a la derecha lo conocemos: Ca0 lo da el ejercicio, al igual que Ca y τ
 b) Ahora si nos piden encontrar esa K y ese orden de reacción n
𝑃𝑎𝑟𝑎 𝑟𝑒𝑎𝑐𝑡𝑜𝑟 𝐶𝑆𝑇𝑅:
𝐶𝐴0 − 𝐶𝐴
= 𝑘𝐶𝐴𝑛
𝑟
Tenemos que linealizar:
𝐶𝐴0 − 𝐶𝐴
ln( ) = ln 𝑘 + 𝑛ln(𝐶𝐴)
𝑟
Esta es una función lineal de tipo Y=mx+b, donde b = ln(k), x=ln(Ca), y n = m, por lo que si graficamos en
polymath o en Excel y pedimos que nos de la función, tenemos los datos:

De Excel o polymath podemos sacar que b= -0,2886

Y como sabemos que b=ln(k)
Entonces:
−0,2886 = ln 𝑘
𝑒−0,2886 = 𝑘
0,7493 = k
Ahora para el orden de reacción:
n=m Y como vemos en la gráfica,
m=1,497
Por lo que podríamos decir que el orden de reacción n es 1,5

Y listo ese es todo el primer punto. No𝑦es necesario hacerlo a pc, también pueden sacar la pendiente a mano
2−𝑦1
después de tabular los datos con 𝑚 = y después armar su función y=mx+b y encontrar b despejando:
𝑥2−𝑥1
−2,618667 − 1,458615
𝑚= = 1,497
−1,560648 − 1,163151
𝑦 = 𝑚𝑥 + 𝑏
𝑝𝑎𝑟𝑎 𝑦 = 1,458615, x = 1,163151
1,458615 = 1,497 ∗ 1,163151 + b
𝑏 = −0,2886
−0,2886 = ln 𝑘
𝑒−0,2886 = 𝑘
0,7493 = k
Son dos formas diferentes de hacer lo mismo y dan iguales resultados
El acetato de etilo es un disolvente que se usa mucho y que se puede formar por la esterificación en fase de
vapor de ácido acético y etanol. La reacción se estudió usando una resina microporosa como catalizador en
un reactor de lecho empacado. La reacción es de primer orden respecto al etanol y de pseudo-orden cero
respecto al ácido acético, con una velocidad de alimentación equimolar de ácido acético y etanol, la velocidad
de reacción especifica es de 1.2 dm3/g cat. min. la velocidad de alimentación molar total es de 10 mol/min, la
presión inicial es de 10 atm. La temperatura es de 118 °C y el parámetro de caída de presión, α, es igual a 0.01
g-1
a. Calcule el peso máximo de catalizador que se podría usar manteniendo una presión en la salida de más de
1 atm.
b. Determine el peso de catalizador necesario para alcanzar una conversión de 90%. determine el cociente
del peso de catalizador necesario para lograr el ultimo 5% (85 a 90%) de conversión, entre el peso
necesario para lograr el primer 5% (0 a 5%) de conversión en el reactor
RESPUESTA
Esta pregunta es de reactores PBR, recuerden que las fórmulas para estos reactores SOLAMENTE EN FASE
GASEOSA son dos: Una para la caída de presión y otra para el diseño del PBR (reitero que solo para PBR en
fase gaseosa se usan esas dos, en fase liquida solo se usa solo la ecuación de diseño porque no hay caída de
presión)

Ecuación de diseño de un reactor PBR:

𝑑𝑥
𝐹𝑎0 = −𝑟´𝑎
𝑑𝑊
𝐷𝑜𝑛𝑑𝑒 𝑤 𝑒𝑠 𝑒𝑙 𝑝𝑒𝑠𝑜 𝑑𝑒𝑙 𝑐𝑎𝑡𝑎𝑙𝑖𝑧𝑎𝑑𝑜𝑟
Ecuación de caída de presión:
𝑑𝑦 α 𝑇
=− (1 + ex)
𝑑𝑊 2y 𝑇0
𝑃
𝐷𝑜𝑛𝑑𝑒 𝑦 = 0
𝑃
Ahora si procedemos a hacer el ejercicio
Lo primero es tener clara la reacción que se está dando:
Ácido acético + Etanol -> Acetato de etilo + agua (según busqué en Google)

Podríamos decir que: A+B produce C+D

Ahora para el inciso a:

a) Nos están pidiendo calcular el peso del catalizador dada cierta presión. Si analizamos las dos
ecuaciones para los PBR, vemos que la fórmula que sirve es la de caída de presión, puesto a que la
ecuación de diseño utiliza una conversión X la cual no tenemos, pero sí tenemos la caída de presión P
que nos da el ejercicio, por lo que usaremos esa.
 𝑑𝑦 α 𝑇
=− (1 + ex)
𝑑𝑊 2y 𝑇0
𝑆𝑢𝑝𝑜𝑛𝑒𝑚𝑜𝑠 𝑒𝑠𝑡𝑎𝑑𝑜 𝑖𝑠𝑜𝑡𝑒𝑟𝑚𝑖𝑐𝑜
𝑑𝑦 α
= − (1 + ex)
𝑑𝑊
2y
𝑃𝑎𝑟𝑎 𝑒:
𝑒 = 𝑦𝑎0(𝛿)
𝑐𝑜𝑚𝑜 𝑙𝑎 𝑟𝑒𝑎𝑐𝑐𝑖ó𝑛 𝑒𝑠 A + B → C + D
𝛿 =1+1−1−1=0
𝑒=0
𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠:
𝑑𝑦 α
=−
𝑑𝑊 2y
−2ydy = αdW
𝑦 𝑤
∫ −2ydy = α ∫ 𝑑𝑊
1 0
1 − y2 = αW
1 − y2
=W
α 𝑃
Sabemos que 𝑦 =
𝑃0
El ejercicio dice que P0 = 10atm
Y P no debe ser inferior a 1atm, por lo que:
1atm 2
1−( )
10atm
0,01 g−1 = W
99g = W
El peso máximo que puede tener el catalizador teniendo una presión de salida de 1atm son 99 gramos

b) Ahora nos dan la conversión y nos piden encontrar el peso del catalizador, por lo que si tenemos
que usar la ecuación de diseño:
𝑑𝑥
𝐹𝑎0 = −𝑟´𝑎
𝑑𝑊
El ejercicio nos dice que la reacción es de primer orden respecto al etanol y de orden casi 0 respecto al
Ácido acetico

−𝑟´𝑎 = 𝑘𝐶𝑎𝐶𝑏0
−𝑟´𝑎 = 𝑘𝐶𝑎

Estamos en fase gaseosa:

1−x 𝑃
𝐶𝑎 = 𝐶𝑎0
1 + ex 𝑃0
𝑑𝑥 1−x 𝑃
𝐹𝑎0 = 𝑘 (𝐶𝑎0 )
𝑑𝑊 1 + ex 𝑃0
𝑃
𝑆𝑎𝑏𝑒𝑚𝑜𝑠 𝑞𝑢𝑒 = 𝑦
𝑃0
 𝑑𝑥 1−x
𝐹𝑎0 = 𝑘 (𝐶𝑎0 ∗ 𝑦)
𝑑𝑊 1 + ex
𝐷𝑒𝑙 𝑝𝑢𝑛𝑡𝑜 𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟, 𝑠𝑎𝑏𝑒𝑚𝑜𝑠 𝑞𝑢𝑒:
1 − y2 = αW
1
𝑦 = (1 − αW)2
𝑑𝑥 1−x 1

𝐹𝑎0 = 𝑘𝐶𝑎0 ( ) (1 − αW)2

𝑑𝑊 1 + ex
𝐸𝑠𝑡𝑎 𝑝𝑎𝑟𝑡𝑒 𝑒𝑠 𝑢𝑛 𝑝𝑜𝑐𝑜 𝑗𝑜𝑑𝑖𝑑𝑎
𝑑𝑥 𝑘𝐶𝑎0 1
= (1 − αW)2𝑑𝑊
1−x 𝐹𝑎0
(1 + ex
)𝑥 𝑤
(1 + 𝑘𝐶𝑎0 1
ex) 𝑑𝑥 = ∫ (1 − αW)2 𝑑𝑊
∫ 𝐶𝑎0𝑣0 0
1−x
0
𝐼𝑛𝑡𝑒𝑔𝑟𝑎𝑙𝑒𝑠 𝑦𝑎 𝑟𝑒𝑠𝑢𝑒𝑙𝑡𝑎𝑠 𝑒𝑛 𝑒𝑙 𝑓𝑜𝑔𝑙𝑒𝑟
1 𝑘 2 3
(1 + e) ln − 𝑒𝑥 = [1 − (1 − αW)2]
1−x 𝑣0 3α
𝑅𝑒𝑐𝑢𝑒𝑟𝑑𝑒𝑛 𝑞𝑢𝑒 𝑛𝑜𝑠 𝑝𝑖𝑑𝑒𝑛 𝑤 𝑐𝑢𝑎𝑛𝑑𝑜 𝑥 = 0,9
𝑌𝑎 𝑠𝑎𝑏𝑒𝑚𝑜𝑠 𝑞𝑢𝑒 𝑒 = 0 𝑝𝑜𝑟 𝑒𝑙 𝑝𝑢𝑛𝑡𝑜 𝑎
dm3
1 (1,2 g cat ∗ min) 2 3
(1 + 0) ln −0=
1 − 0,9 𝑣0 [1 − (1 − αW)2]
3 ∗ 0,01g−1
𝐹𝑎𝑙𝑡𝑎 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 𝑐𝑎𝑢𝑑𝑎𝑙, 𝑝𝑒𝑟𝑜 𝑝𝑜𝑑𝑒𝑚𝑜𝑠 𝑠𝑎𝑐𝑎𝑟𝑙𝑜 𝑐𝑜𝑛 𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑃0𝑣0 = 𝑛𝑅𝑇0
𝑣0 𝑛𝑅𝑇0 𝑚𝑜𝑙
= → 𝑛=5 (𝑙𝑎 𝑚𝑖𝑡𝑎𝑑 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 𝑡𝑜𝑡𝑎𝑙), 𝑇 = 391,15K, P = 10 atm
𝑃0𝑚𝑜𝑙 𝑚𝑛
𝑖 atmL
( 𝑚𝑖𝑛
5 ) (0,08205746 𝑚𝑜𝑙 𝐾 ) (391,15K)
𝑣𝑜 = = 16,0484 L
10 atm 𝑚𝑖𝑛

Ahora si volvemos a la ecuación para despejar w, ay mi madre el bicho

(1,2 dm3 ) ∗ 1L
1 g cat ∗ min 1dm3 2 3
(1 + 0) ln −0=
1 − 0,9 L [1 − (1 − αW)2]
16,0484 3 ∗ 0,01g−1
3
2,302585 = 4,985 [1 − (1 − αW)2]
3
0,46191 = 1 − (1 − αW)2
3
1 − 0,46191 = (1 − αW)2
2
0,538093 = 1 − αW
1 − 0,66156
=W
0,01g−1
33,8439g = W
Piden encontrar el cociente entre el peso para una conversión de 0% a 5% y el peso a una conversión entre
85% y 90%, asi que lo único que hay que hacer es cambiar los limites de integración en la integral de la
conversión:

Para x de 0% a
5%: 0,05
(1 + ex)
∫ 𝑘 𝑤 1
𝑑𝑥 = ∫ (1 − αW)2 𝑑𝑊
1−x 𝑣0 0
0
dm3
1 (1,2 g cat ∗ min) 2 3
ln −1
)W) 2]
1 − 0,05 L [1 − (1 − (0,01g
= 16,0484 𝑚𝑖𝑛 3(0,01g−1)
3
0,0103 = 1 − (1 − (0,01g−1)W)2
0,9931 = 1 − (0,01g−1)W
0,6871g = W para x de 0% a 5%
Para x de 85% a 90%:
Ya Conocemos el peso del catalizador para x=90%, asi que encontramos el peso cuando x=85%
0,85
(1 + 𝑊
ex) 𝑘 1
∫ 𝑑𝑥 = ∫ (1 − αW)2 𝑑𝑊
1−x 𝑣0 0
0
dm3
1 (1,2 g cat ∗ min) 2 3
ln −1
)W) 2]
1 − 0,85 L [1 − (1 − (0,01g
= 16,0484 𝑚𝑖𝑛 3(0,01g−1)
3
0,38057 = 1 − (1 − (0,01g −1
)W)2

0,72665 = 1 − (0,01g−1)W
27,334g = W Para x de 0% a 85%

Para x de 85% a 90% solo hay que restar los dos pesos:
𝑊𝑥=0,85→ 0,9 = 𝑊𝑥=0,9 − 𝑊𝑥=0,85 = 33,8439g − 27,334g = 6,509244g
𝐸𝑙 𝑐𝑜𝑐𝑖𝑒𝑛𝑡𝑒 𝑠𝑒𝑟𝑖𝑎:
𝑊𝑥=0→ 0,05 0,6871g
𝑊𝑥=0,85→ 0,9 = = 0,10555
6,509244g

Joa me duele la cabeza

La reacción de primer orden, A→B, es llevada a cabo en un PBR. La constante de velocidad de reacción y la ley
de velocidad de reacción está dada por:
𝑚𝑜𝑙
−𝑟𝐴´ = 𝑘𝑝𝐴
𝑘𝑔 𝑐𝑎𝑡 ℎ
 𝑚𝑜𝑙
𝑘 = 0.80
𝑎𝑡𝑚 𝑘𝑔 𝑐𝑎𝑡 ℎ
Al reactor es alimentado A con 40 % inertes a 350 ⁰C y 2 atmosfera. El flujo del alimento es 40 moles/h. El
parámetro de la caída de presión es α = 0.005 kg-1.
a. Calcule la conversión alcanzada con un peso de catalizador de 150 kg.
b. Calcule la caída de presión.

Respuesta

B) Es mas fácil hacer el b primero, se hace igual al punto anterior

𝑑𝑦 α 𝑇
=− (1 + ex)
𝑑𝑊 2y 𝑇0
𝑒 = (1 − 1)(ya0) = 0
𝑑𝑦 α
=−
𝑑𝑊 2y
1
𝑦 = (1 − αW)2
1
𝑦 = (1 − 0,005kg−1 ∗ 150kg)2
𝒚 = 𝟎, 𝟓
𝑷 = 𝟎, 𝟓(𝟐 𝐚𝐭𝐦) = 𝟏𝐚𝐭𝐦

A) Ahora si el punto difícil

𝐹𝑎𝑜 = 𝐹0 ∗ 0,6 (𝐸𝑛 𝑙𝑎 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑒𝑠𝑡á 𝐴 𝑦 𝑢𝑛 40% 𝑑𝑒 𝑖𝑛𝑒𝑟𝑡𝑒𝑠)

𝐹𝑎0 = 40 𝑚𝑜𝑙 mol
∗ 0,6 = 24 h
ℎ
𝑑𝑥
𝐹𝑎0 = −𝑟´𝑎
𝑑𝑊
𝑑𝑥
𝐹𝑎0 = 𝑘𝑝𝐴
𝑃𝑜𝑟𝑑𝑊𝑔𝑎𝑠𝑒𝑠 𝑖𝑑𝑒𝑎𝑙𝑒𝑠:
𝑃𝑎 ∗ 𝑣 = 𝑛𝑎 ∗ 𝑅 ∗ 𝑇
𝑛𝑎𝑅𝑇
𝑃𝑎 =
𝑃𝑎 = 𝐶𝑎𝑣𝑅𝑇
𝑑𝑥
𝐹𝑎0 = 𝑘𝐶𝑎𝑅𝑇
𝑑𝑊 1−x 𝑃
𝐹𝑎𝑠𝑒 𝑔𝑎𝑠𝑒𝑜𝑠𝑎: 𝐶𝑎 = 𝐶𝑎0
𝑑𝑥 1 − x1𝑃+ ex 𝑃0
𝐹𝑎0 = 𝑘𝑅𝑇𝐶𝑎0
𝑑𝑊 1 + ex 𝑃0
𝑅𝑇𝐶𝐴0 = 𝑝𝐴0
𝑃 1
= 𝑦 = (1 − αW)2
𝑃0
 𝑑𝑥 1−x 1

𝐹𝑎0 = 𝑘𝑝𝐴0 (1 − αW)2

𝑑𝑥𝑑𝑊 𝑘𝑝𝐴0 1 + ex 1
= (1 − αW)2𝑑𝑊
1−x 𝐹𝑎0
1 𝑥+ ex
1 + ex 𝑘𝑝𝐴0 𝑤 1
∫ 𝑑𝑥 = ∫ (1 − αW) 𝑑𝑊
2
0 1−x 𝐹𝑎0 0

1 𝑘𝑝𝐴0 2 3
(1 + e) − 𝑒𝑥 = [1 − (1 − αW)2]
ln 1−x 𝐹𝑎0 3α

e=0
w = 150 kg
𝑝𝑎0 = 𝑃 ∗ 𝑦𝑎0 = 2atm ∗ 0,6 = 1,2 atm
𝑚𝑜𝑙
0,80 ∗ 1,2 atm 2
1 𝑎𝑡𝑚 𝑘𝑔 𝑐𝑎𝑡 ℎ −1 3
ln mol ∗ [1 − (1 − 0.005 kg ∗ 150 kg)2]
1=− x 3(0,005 kg−1)
24
1 h
ln =
1−x
4,6666
1
4,6666
1−x= e
1
= 106,3426754 despejando x
1 − x
x = 0,9905

Este fue casi lo mismo que el anterior asi que no fui tan detallista. Espero que se entienda

4. Una reacción química, 2A →B, es llevada a cabo en un reactor por lotes para estudios cinéticos. De los
siguientes valores de concentración, determine el orden de reacción y la constante de velocidad de
reacción.
t (min) 0 5 13 20 30
CA (mol/lit) 5 0.93 0.41 0.27 0.18

Respuesta

Este se parece al primero. Nos piden encontrar el orden de la reacción y el valor de k a apartir de unos datos,
entonces la forma de proceder es la misma:

Partimos del balance de moles para un reactor Batch

𝑑𝑁𝑎
= −𝑟𝑎𝑉
𝑑𝑡
Asumimos Volumen constante

𝑑𝑁𝑎
= −𝑟𝑎
𝑉 ∗ 𝑑𝑡
𝑁𝑎
= 𝐶𝑎
𝑉
 𝑑𝐶𝑎
= 𝑘𝐶𝑛
𝑑𝑡 𝑎

Linealizamos
𝑑𝐶𝑎
ln = ln 𝑘 + 𝑛 ln 𝐶𝑎
𝑑𝑡
𝑦 = 𝑚𝑥 + 𝑏
𝑏 = ln 𝑘
𝑚=𝑛
𝑥 = ln 𝐶𝑎

Hay que obtener los valores de dCa/dt, y esto lo haremos con un método llamado diferencias finitas
Sabemos que un diferencial es casi lo mismo que un delta, asi que decimos que:

𝑑𝐶𝑎 𝐶𝑎2 − 𝐶𝑎1

=
𝑑𝑡 𝑡2 − 𝑡1

(Nota: Para sacar ln(dCa/dt), tomen el valor como positvo, pues no existe el ln de un numero negativo)

De aquí podemos ver que:

𝑚 = 𝑛 = 2,7893. El orden de la reacción es ese, pero podemos decir que es 3

𝑏 = −0,1093 = ln 𝑘
𝑒−0,1093 = 𝑘 = 0,89646

Tambien se puede hacer a mano como el primer punto, pero los valores van a dar un poco diferentes porque
la linea de tendencia no está tan marcada

−4,7105307 + 0,20579491
𝑚= = 2,74
−1,71479843 + 0,07257069
−4,7105307 = (−1,71479843)2,74 + b
𝑏 = −0,01198
𝑘 = 0,988