Ingeniería de las reacciones químicas.

Lizeth Dayana León López 317045

PRUEBA 3
Ejercicio individual
El propilenglicol (C3H8O2) se obtiene, a temperatura ambiente, por hidrólisis del óxido de propileno
(OP, C3H6O), en presencia de ácido sulfúrico. Un ingeniero químico está encargado de supervisar el
funcionamiento de un CSTR adiabático en el que se produce propilenglicol a escala industrial por el
método descrito. Desgraciadamente el reactor está comenzando a deteriorarse, puesto que el ácido
sulfúrico ataca al acero con el que está construido, y hay que reemplazarlo por uno nuevo. Buscando
en el almacén de la empresa, el ingeniero químico encuentra un CSTR de 1.14 m3 de volumen útil
(con un revestimiento interior de vidrio que impedirá que el ácido sulfúrico ataque al acero con el que
está construido el reactor) y que se podría utilizar para producir el propilenglicol. En la instalación
actual se están alimentando 19.54 kmol/h de OP. El alimento consta de dos corrientes:
Corriente A: una mezcla equivolumétrica de OP (1.32 m3/h) y de metanol (1.32 m3/h); y
Corriente B: una corriente de agua que contiene un 0.1% de ácido sulfúrico, con un caudal
volumétrico 6.6 m3 /h.
El caudal molar de alimentación del metanol es 32.63 kmol/h y el del agua 364.5 kmol/h. La
temperatura de las dos corrientes es 14.5°C, antes de mezclarlas, pero hay un incremento instantáneo
de 9.5°C en la corriente que resulta debido al elevado calor de mezcla. En consecuencia, la
temperatura de entrada del alimento al reactor es 24°C. La reacción es de primer orden respecto a la
concentración de OP y de orden aparente cero respecto al agua que se utiliza en exceso. La variación
de la constante cinética con la temperatura es de la forma:

 −75362  −1
k = 16.96 1012  exp   h con la energía de activación expresada en J/mol
 RT 
En la operación hay una importante restricción: el OP es una sustancia de bajo punto de ebullición (a
1 atm. abs. hierve a 34.3 °C). Por lo tanto, con la mezcla que se está utilizando, el reactor no debe
sobrepasar los 57 °C para evitar las pérdidas de OP a través del sistema de ventilación del reactor.
Datos
FA,0 = 19540 mol / h CpA = 146.5 J v0, A = 1.32 m3 / h
mol.K
FB ,0 = 364500 mol / h CpB = 75.35 J v0, B = 6.6 m3 / h
mol.K
FC ,0 = 0 mol / h CpC = 192.6 J v0,C = 0 m3 / h
mol.K
FM ,0 = 32630 mol / h CpM = 81.63 J v0, M = 1.32 m3 / h
mol.K
Tmax = 330 K H rxn = −85300 J / mol VR = 1.14 m3
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

Reacción.

C3 H 6O + H 2O ⎯⎯⎯
H 2 SO4
→ C3 H8O2 El metanol (M) actúa como un inerte.
La presencia del catalizador (H2SO4) desencadena la
A + B ⎯⎯⎯
H 2 SO4
→C reacción, sin embargo, este se alimenta en muy baja
concentración respecto al agua, reactivo altamente en
exceso. Por tanto, se considera que la cantidad de
catalizador tiende a cero.

Ley de velocidad de reacción rA = −kC A

FA
CA =
vT ,0 FA,0  (1 − X A )
rA = −k
FA,0  (1 − X A ) vT ,0
CA =
vT ,0

 −75362  −1
Donde k = 16.96 1012  exp   h
 RT 
Tabla estequiométrica

Sustancia Inicio Cambio por rxn Final

A FA,0 −FA,0 X A FA,0 ( A − X A )
B FB,0 −FA,0 X A FA,0 (B − X A )
C - + FA,0 X A FA,0 (C + X A )
M FM ,0 - FA,0 (M − X A )
Total FT ,0 −FA,0 X A FT ,0 − FA,0 X A

Cálculos preliminares

vT ,0 = v0, A + v0, B + v0,C + v0, M

• Flujo volumétrico total de alimentación
vT ,0 = 9.24 m3 / h
• Capacidad calorífica de alimentación
N
Cp i i A A B B C C M M
ˆ =  Cp =  Cp +  Cp +  Cp +  Cp
i =1

ˆ = FA,0 Cp + FB,0 Cp + FC ,0 Cp + FM ,0 Cp
Cp A B C M
FA,0 FA,0 FA,0 FA,0
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

ˆ = FA,0 Cp + FB,0 Cp + FC ,0 Cp + FM ,0 Cp
Cp A B C M
FA,0 FA,0 FA,0 FA,0

ˆ = 1688.4 J
Cp
mol  K
ˆ = 146.5 J + 364500 mol / h  75.35 J  + 0 + 32630 mol / h  81.63 J 
Cp    
mol  K 19540 mol / h  mol  K  19540 mol / h  mol  K 

ˆ = 1688.4 J
Cp
mol  K

VR 1.14m3
• Tiempo de residencia = =
vT ,0 9.24m3 / h
 = 0.1234h
Deducciones matemáticas

• Balance molar E + G – S = Ac
En un reactor CSTR E+G–S=0
Fi ,0 + rV
i R − Fi = 0
Para A FA,0 + rAVR − FA = 0
De la tabla estequiométrica FA = FA,0  (1 − X A )
Reemplazando rAVR + FA,0 X A = 0
Reemplazando la ley de F  (1 − X A )
velocidad de reacción −k A,0 VR + FA,0 X A = 0
vT ,0
Reorganizando V
−k (1 − X A ) R + X A = 0
vT ,0
VR −k (1 − X A ) + X A = 0
Con  =
vT ,0
Con un poco de álgebra k
X A, BM =
1 + k (1)

• Balance de energía E–S=A

dE
Eent − Esal =
dt
En estado estable Eent − Esal =0
Flujos de energía que afectan el sistema Q − W + Fent Eent − Fsal Esal = 0
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

La velocidad de energía total que se adiciona N N

y se remueve del sistema está descrita por la Q − W +  Fi Ei |ent − Fi Ei | sal = 0

i =1 i =1
sumatoria del aporte de cada sustancia que
participa.

La velocidad de flujo de trabajo representa N N

la energía necesaria para el ingreso y la W = − FPV

i i | ent +  Fi PVi | sal +Weje
i =1 i =1
salida de la masa del sistema y el trabajo de
eje.

El enunciado no hace referencia a un N N

agitador o turbina, por lo tanto, Weje = 0

W = − FPV
i i | ent +  Fi PVi | sal
i =1 i =1

Reemplazando
N N N N
Q +  Fi PVi |ent − Fi PVi | sal + Fi Ei |ent − Fi Ei | sal = 0
i =1 i =1 i =1 i =1
Factorizando
N N
Q +  Fi ( Ei + PVi ) |ent − Fi ( Ei + PVi ) | sal = 0
i =1 i =1
La energía interna representa el mayor aporte de energía en el sistema
N N
Q +  Fi (U i + PVi ) | ent − Fi (U i + PVi ) | sal = 0
i =1 i =1

Como H = U + PV
N N
Q +  Fi ( H i ) | ent − Fi ( H i ) | sal = 0
i =1 i =1

Entrada N

 F (H ) |
i =1
i i ent
= FA,0 H A,0 + FB,0 H B,0 + FC ,0 H C ,0 + FM ,0 H M ,0

Salida N

 F (H ) |
i =1
i i sal
= FA H A + FB H B + FC H C + FM H M

Reemplazando los flujos de cada sustancia a partir de la tabla estequiométrica y reagrupando

N N

 F (H ) |
i =1
i i ent
− Fi ( H i ) | sal = FA,0 ( H A,0 − H A ) +  B ( H B,0 − H B ) + C ( H C ,0 − H C ) + M ( H M ,0 − H M )  ...
i =1

... − ( H C − H B − H A ) FA,0 X A
N N

 F (H ) |
i =1
i i ent
− Fi ( H i ) | sal = FA,0 ( H A,0 − H A ) +  B ( H B,0 − H B ) + C ( H C ,0 − H C ) + M ( H M ,0 − H M )  ...
i =1

... − H Rxn FA,0 X A

Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

N N N

 Fi (Hi ) |ent − Fi (Hi ) | sal = FA,0 i ( Hi,0 − Hi ) − H Rxn FA,0 X A

i =1 i =1 i =1

Retomando N N
Q +  Fi ( H i ) | ent − Fi ( H i ) | sal = 0
i =1 i =1

Reemplazando N
Q + FA,0 i ( H i ,0 − H i ) − H Rxn FA,0 X A = 0
i =1

La entalpía H i (T ) = H i (T0 ) + Cpi (T − Ti ,0 )

H i (T0 ) − H i (T ) = Cpi (Ti ,0 − T )

Reemplazando N
Q + FA,0 i Cp(Ti ,0 − T ) − H Rxn FA,0 X A = 0
i =1
Todas las sustancias ingresan a la misma
temperatura (T0). ˆ (T − T )
Q + FA,0Cp
Finalmente XA = 0

H Rxn FA,0
La reacción es exotérmica, por lo tanto, ˆ (T − T )
−Qsal + FA,0Cp (2)
el calor es un flujo de salida. X A, BE = 0

H Rxn FA,0

a) ¿Puede utilizarse el CSTR del almacén para sustituir el que se está deteriorando,
funcionando de forma adiabática? Si se puede emplear, ¿cuál será la conversión de OP
en propilenglicol en estado estacionario?
Retomando (1) k
X A, BM =
1 + k
Reemplazando datos k (0.1234h) (3)
X A, BM =
1 + k (0.1234h)
Con constante de
 −75362 J 
k = 16.96 1012  exp  mol  h −1
velocidad de reacción
 8.314 J  (T ) 
 mol  K 
Retomando (2) −Qsal + FA,0Cp ˆ (T − T )
X A, BE = 0

H Rxn FA,0
En una operación ˆ (T − T )
Cp
adiabática Q = 0 X A, BE = 0

H Rxn
Reemplazando datos con 1688.4 J / mol  K (297 K − T ) (4)
una T0 = 297 K X A, BE =
−85300 J / mol
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

Se grafica la conversión del balance molar y la conversión del balance de energía en un intervalo de
293 a 353 K, obteniendo la Figura 1, donde la intersección representa las condiciones de operación
del reactor en estado estacionario adiabático; se observa que la temperatura es mayor a la temperatura
máxima admitida.

Figura 1. Operación adiabática.

Se calcula la temperatura de operación.

El reactor opera a la temperatura en la que
la conversión en el balance molar es igual X A, BM = X A, BE
a la conversión en el balance de energía
(estado estacionario).

Igualando a cero y reemplazando (3) y (4)

k (0.1234h−1 ) 1688.4 J / mol  K (297 K − T )
−1
− =0
1 + k (0.1234h ) −85300 J / mol

Con constante de velocidad de reacción

 −75362 J 
k = 16.96 10  exp
12  mol  h −1
 8.314 J  (T ) 
 mol  K 
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

Se resuelve utilizando la herramienta fsolve de Matlab®, encontrando una temperatura de operación

de 339.5 K y, reemplazando este dato, se obtiene una conversión en estado estacionario de 0.8420,
sin embargo, debido a que la temperatura de operación máxima admitida es de 330 K no se puede
utilizar el CSTR en operación adiabática porque habría pérdidas del producto deseado.
b) En el mismo almacén se encuentra también disponible un serpentín que tiene un área
útil de transferencia de calor de 5 m2. Si el agua de refrigeración que se utiliza circula
a través del serpentín con un caudal suficientemente grande para admitir que su
temperatura permanece constante e igual a 29.5 °C, y para este caudal y las condiciones
del reactor se sabe que el coeficiente global de transmisión de calor vale 1500
kJ/(h·m2·K), ¿podría instalarse este serpentín en el reactor del almacén para llevar a
cabo la reacción, con las restricciones antes citadas, en el caso de que se incremente a
27 °C la temperatura de entrada del alimento? ¿Qué conversión y temperatura de
estado estacionario resultarían?
Datos
T0 = 297 K Ta = 302.5 K A = 5m 2 U = 1500 103 J
h  m2  K

Retomando (1) k
X A, BM =
1 + k
Reemplazando datos k (0.1234h) (3)
X A, BM =
1 + k (0.1234h)
Con constante de (4)
 −75362 J 
velocidad de reacción
k = 16.96 10  exp
12  mol  h −1
 8.314 J  (T ) 
 mol  K 
Retomando (2) ˆ (T − T )
−Qsal + FA,0Cp
X A, BE = 0

H Rxn FA,0
Intercambio de calor con
el serpentín Qsal = U  A  (T − Ta )
Reemplazando ˆ (T − T )
−U  A  (T − Ta ) + FA,0Cp
X A, BE = 0

H Rxn FA,0
Reemplazando datos
J mol J
1500 103  5m 2  (T − 302.5) K + 19540 11688.4  (297 − T ) K (5)
X A, BE = hm K
2
h mol  K
J mol
−85300 19540
mol h

Se grafica la conversión del balance molar y la conversión del balance de energía en un intervalo de
293 a 353 K, obteniendo la Figura 2, donde la intersección representa las condiciones de operación
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

del reactor en estado estacionario con intercambio de calor. Se observa que es menor a la temperatura
máxima admitida.

Figura 2. Operación con intercambio de calor.

Se calcula la temperatura de operación.

El reactor opera a la temperatura en la que X A, BM = X A, BE
la conversión en el balance molar es igual
a la conversión en el balance de energía.

Igualando a cero X A, BM − X A,BE = 0

Ecuaciones (3) y (5)

Se resuelve utilizando la herramienta fsolve de Matlab®, encontrando una temperatura de operación

(T) de 329.1 K, teniendo en cuenta que la temperatura de operación máxima admitida es de 330 K se
puede utilizar el CSTR con intercambio de calor (serpentín). Sin embargo, debe tenerse en cuenta que
se utilizó un calor de reacción constante en el intervalo de temperatura, debido a que el enunciado no
presenta la temperatura a la cual se reporta este dato, de modo que con un dato más específico y
corrigiendo por temperatura, es posible que la temperatura de operación supere los 330 K, haciendo
inviable la utilización del reactor. Adicionalmente, la más mínima perturbación del sistema (p.e.
Cambios momentáneos que conduzcan al incremento de la temperatura de alimento y del tiempo de
Ingeniería de las reacciones químicas.
Lizeth Dayana León López 317045

residencia) podría ocasionar una temperatura de operación mayor a 330K y, con esto, ebullición de
la mezcla y pérdida de materia a través del sistema de ventilación del reactor.
Se calcula la conversión
Reemplazando T = 329.1 K
 −75362 J 
k = 16.96 1012  exp  mol  h −1
en (4)
 8.314 J  (T ) 
 mol  K 
−1
k = 18.4824h
Reemplazando k en (5) 18.4824h −1 (0.1234h −1 )
XA =
1 + 18.4824h −1 (0.1234h −1 )

Conversión X A = 0.6951

Tabla de resultados
Operación Temperatura Conversión Límite máximo de
temperatura
Adiabática 339.5 K 0.8420
Con intercambio de 329.1 K 0.6951 330 K
calor (serpentín)