UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU

FACULTAD DE INGENIERIA QUÍMICA

REACTOR PFR EN ESTADO ESTACIONARIO PARA LLEVAR A CABO UNA

REACCIÓN ENDOTÉRMICA EN FASE GASEOSA

REACCIONES QUIMICAS II

PRESENTADO A:

ING. Salvador Ore Vidalon

POR:

NUÑEZ HERNANDEZ, Angel

HUANCAYO - PERU

2016
 Introducción
En un reactor PFR en estado estacionario para llevar a cabo una reacción
endotérmica en fase gaseosa, se presenta como uno de los pasos claves es la
desintegración catalítica en fase vapor de acetona a ceteno y metano que se
emplea en un tratamiento del diseño de una instalación para fabricar anhídrido
acético.
Planteamos una situación en la cual se desea evaluar cuál será el
comportamiento del PFR un reactor que consta de una banco de 1000 tubos
cedula 40 de una pulgada. Para unas condiciones iniciales de operación y
considerando tres casos: (1) El reactor se opera adiabáticamente, (2) El reactor
está rodeado por un intercambiador de calor y (3) El intercambiador de calor en
el caso 2 tiene ahora la temperatura de intercambio de calor variable. Esto con
el fin de obtener las datos apropiadas para la temperatura con la cual se pueda
obtener una conversión deseada analizando la influencia que tiene la
temperatura con respecto a la conversión, establecer métodos los cuales
permitan corregir las deficiencias, verificar el material del tubo y la cedula del
reactor para ver si son convenientes para tales temperaturas y presiones
sugiriendo las más adecuadas para su buen funcionamiento
 Marco teórico
1. . PFR en estado estacionario con cambiador de calor
En este trabajo consideraremos un PFR en el cual el calor se agrega o retira a
través de las paredes cilíndricas del mismo (figura 1), Para elaborar un modelo
de reactor, asumiremos que no hay gradientes radiales en el mismo y que el
flujo de calor a través de la pared por volumen unitario del reactor es como se
indica en la figura 1.

Figura 1: PFR con ganancia o pérdida de calor.

1.1. Deducción del balance de energía para un PFR
A continuación consideraremos sistemas de flujo que se operan en estado
estacionario.
El balance de energía en estado estacionario se obtiene igualando (𝑑𝐸̂𝑠𝑖𝑠𝑡𝑒𝑚𝑎 /
𝑑𝑡) a cero en la ecuación para dar:
𝑛̇ 𝑛

𝑄̇ − 𝑊𝑆 + ∑ 𝐹𝑖0 𝐻𝑖0 − ∑ 𝐹𝑖 𝐻𝑖 = 0 (1)

𝑖=1 𝑖=1

Efectuaremos un balance de energía para el volumen ∆𝑉 con 𝑊𝑆 = 0 La

ecuación (1) se transforma:

∆𝑄̇ + 𝐹𝑖 𝐻𝑖 |𝑉 − 𝐹𝑖 𝐻𝑖 |𝑉+∆𝑉 = 0 (2)

El flujo de calor al reactor ∆𝑄, se da en ténninos del coeficiente total de

transferencia de calor, U; el área de intercambio de calor, ∆𝐴 y la diferencia
entre la temperatura ambiente, Ta y la temperatura del reactor, T.
∆𝑄 = 𝑈∆𝐴(𝑇𝑎 − 𝑇) = 𝑈𝑎∆𝑉(𝑇𝑎 − 𝑇)

Donde a es el área de intercambio de calor por unidad de volumen del reactor.

Para el PFR
 𝐴 𝜋𝐷𝐿 4
𝑎= = 2 =
𝑉 𝜋𝐷 𝐿 𝐷
4
Donde D es el diámetro del reactor. Sustituyendo ∆𝑄̇, en la ecuación (2),
dividiendo la ecuación (2) entre ∆𝑉 y en el límite, cuando ∆𝑉 → 0 , obtenemos
𝑑 ∑(𝐹𝑖 𝐻𝑖 )
𝑈𝑎(𝑇𝑎 − 𝑇) − =0
𝑑𝑉
Expandiendo
𝑑𝐹𝑖 𝑑𝐻𝑖
𝑈𝑎(𝑇𝑎 − 𝑇) − ∑ 𝐻𝑖 − ∑ 𝐹𝑖 𝐻 =0 (3)
𝑑𝑉 𝑑𝑉 𝑖

Efectuando un balance molar para la especie 𝑖, tenemos

𝑑𝐹𝑖
= 𝑟𝑖 = 𝑣𝑖 (−𝑟𝐴 ) (4)
𝑑𝑉

Diferenciando la ecuación de entalpía con respecto a V

𝑑𝐻𝑖 𝑑𝑇
= 𝐶𝑃 𝑖 (5)
𝑑𝑉 𝑑𝑉

Sustituyendo las ecuaciones (4) y (5) en la ecuación (3), obtenemos

𝑑𝑇
𝑈𝑎(𝑇𝑎 − 𝑇) − ∑ 𝑣𝑖 𝐻𝑖 (−𝑟𝐴 ) − ∑ 𝐹𝑖 𝐶𝑃 𝑖 =0
𝑑𝑉
Reordenando, nos da la siguiente expresión,

𝑑𝑇 𝑟𝐴 ∆𝐻𝑅𝑥 − 𝑈𝑎(𝑇𝑎 − 𝑇)
= (6)
𝑑𝑉 ∑ 𝐹𝑖 𝐶𝑃 𝑖

Esta forma del balance de energía también se aplicará a reacciones múltiples

𝑟𝐴 ∆𝐻𝑅𝑥 : Calor generado
𝑈𝑎(𝑇𝑎 − 𝑇): Calor retirado
A continuación expresaremos 𝑟𝐴 en función de las concentraciones para
sistemas líquidos o flujos ·molares para sistemas gaseosos.
Utilizaremos esta forma del balance de energía para reactores de membrana y
también la ampliaremos para reacciones múltiples.
Podríamos escribir la ecuación (6) en términos de la conversión, recordando
que 𝐹𝑖 = 𝐹𝐴0 (Θ𝑖 + 𝑣¡ 𝑋) y sustituyendo la expresión en el denominador de la
ecuación (6).
 𝑑𝑇 𝑈𝑎(𝑇𝑎 − 𝑇) + 𝑟𝐴 ∆𝐻𝑅𝑥
= (7)
𝑑𝑉 𝐹𝐴0 (∑ Θ𝑖 C𝑃𝑖 + ΔC𝑃 𝑋)

Balance de energía para PFR

Para un reactor empacado 𝑑𝑊 = 𝜌𝑏 𝑑𝑉, donde 𝜌𝑏 es la densidad volumétrica
del catalizador,

𝑈𝑎
𝑑𝑇 (𝑇𝑎 − 𝑇) + 𝑟′𝐴 ∆𝐻𝑅𝑥
𝜌 (8)
= 𝑏
𝑑𝑊 ∑ 𝐹𝑖 𝐶𝑃 𝑖

Balance de energía para reactor empacado

Tras haber realizado la derivación de tales ecuaciones será más fácil aplicarlas
con precisión a problemas de ingeniería de reacciones químicas con efectos de
calor. La ecuación diferencial que describe el cambio de temperatura con el
volumen a lo largo del reactor (es decir, la distancia) es

𝑑𝑇
= 𝑔(𝑋, 𝑇) (A)
𝑑𝑉

Y debe acoplarse con el balance molar,

𝑑𝑋 −𝑟𝐴
= = 𝑓(𝑋, 𝑇) (B)
𝑑𝑉 𝐹𝐴0

Para resolverlas simultáneamente. Si la temperatura del enfriador varía a lo

largo del reactor, es necesario agregar el balance del enfriador, que es

𝑑𝑇𝑎 𝑈𝑎 (𝑇𝑎 − 𝑇)
= = 𝑓(𝑋, 𝑇) (C)
𝑑𝑉 𝑚̇ 𝑐 𝐶𝑃𝑐

Es posible emplear diversos métodos numéricos para resolver las ecuaciones

diferenciales acopladas, (A), (B) Y (C).
1.2. Balance para el fluido enfriador en la transferencia de calor
El fluido en la transferencia de calor será un enfriador para reacciones
exotérmicas y un medio de calentamiento para reacciones endotérmicas. Si el
flujo del fluido en la transferencia de calor es suficientemente alta con respecto
al calor liberado (o absorbido) por la mezcla de reacción, entonces la
temperatura del fluido en la transferencia de calor será constante a lo largo del
reactor.
En el material siguiente desarrollaremos ecuaciones básicas para un enfriador
utilizado para retirar calor de reacciones exotérmicas; sin embargo, tales
ecuaciones serán aplicables a reacciones endotérmicas en las cuales se
emplea un medio de calentamiento para suministrar calor.

Por convención, 𝑄̇ es el calor agregado al sistema. A continuación realizaremos

el balance para el enfriador en la sección anular entre R1 y R2 Y entre 𝑉 y 𝑉 +
Δ𝑉. El flujo de masa del enfriador es 𝑚̇ 𝐶 Consideraremos el caso en el que el
radio externo del canal enfriador R2 está aislado.
Caso A: flujo a Co-corriente
El reactivo y el enfriador fluyen en la misma dirección.

El balance de energía para el enfriador en el volumen entre 𝑉 y 𝑉 + Δ𝑉 es

𝐶𝑎𝑙𝑜𝑟 𝑎𝑔𝑟𝑒𝑔𝑎𝑑𝑜 𝑝𝑜𝑟
𝐸𝑛𝑒𝑟𝑔𝑖𝑎 𝑞𝑢𝑒 𝐸𝑛𝑒𝑟𝑔𝑖𝑎 𝑞𝑢𝑒
[ ]−[ ] + [𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖𝑜𝑛 𝑎 𝑡𝑟𝑎𝑣𝑒𝑠 ] = 0
𝑒𝑛𝑡𝑟𝑒 𝑒𝑛 𝑉 𝑠𝑎𝑙𝑒 𝑒𝑛 𝑉 + Δ𝑉
𝑑𝑒 𝑙𝑎 𝑝𝑎𝑟𝑒𝑑 𝑖𝑛𝑡𝑒𝑟𝑛𝑎
𝑚̇𝐶 𝐻𝐶 |𝑉 + 𝑚̇ 𝐶 𝐻𝐶 |𝑉+Δ𝑉 + 𝑈𝑎(𝑇 − 𝑇𝑎 )Δ𝑉 = 0
Donde 𝑇𝑎 es la temperatura del enfriador y T es la temperatura de la mezcla de
reacción en el tubo interno.
Dividiendo entre Δ𝑉 Y tomando el límite cuando Δ𝑉 → 0
𝑑𝐻𝐶
−𝑚̇ 𝐶 + 𝑈𝑎(𝑇 − 𝑇𝑎 ) = 0 (9)
𝑑𝑉

El cambio de entalpía del enfriador puede escribirse como sigue

𝑑𝐻𝐶 𝑑𝑇𝑎
= 𝐶𝑃𝑐 (10)
𝑑𝑉 𝑑𝑉
La variación de temperatura del enfriador 𝑇𝑎 a lo largo del reactor es

𝑑𝑇𝑎 𝑈𝑎(𝑇 − 𝑇𝑎 )
= (11)
𝑑𝑉 𝑚̇ 𝐶 𝐶𝑃𝑐

A continuación se muestran los perfiles de temperatura del fluido en la

transferencia de calor típicos para reacciones exotérmicas y endotérmicas

Caso B: Flujo a contracorriente

En este caso, la mezcla de reacción y el enfriador fluyen en dirección opuesta.
En la entrada del reactor, V = O, los reactivos ingresan a la temperatura 𝑇0 y el
enfriador sale a la temperatura 𝑇𝑎2 . Al final del reactor, los reactivos y los
productos salen a la temperatura T, mientras el enfriador entra a 𝑇𝑎0 .

De nuevo escribiremos un balance de energía sobre un volumen diferencial del

reactor para llegar al volumen del reactor

𝑑𝑇𝑎 𝑈𝑎(𝑇𝑎 − 𝑇)
= (12)
𝑑𝑉 𝑚̇ 𝐶 𝐶𝑃𝑐

En la entrada V = O  X = O Y 𝑇𝑎 = 𝑇𝑎2 .
En la salida V = Vf  𝑇𝑎 = 𝑇𝑎0
Observamos que la única diferencia entre las ecuaciones (11) y (12) es el signo
de menos (𝑇𝑎 − 𝑇) = −(𝑇𝑎 − 𝑇) .

Planteamiento del problema

Producción de anhídrido acético
Jeffreys, en un tratamiento del diseño de una instalación para fabricar anhídrido
acético, indica que uno de los pasos claves es la desintegración catalítica en
fase vapor de acetona a ceteno y metano:

CH3COCH3 CH2CO + CH4

Dice además que tal reacción es de primer orden con respecto a acetona y que
la velocidad de reacción específica puede expresarse como sigue

34.222
𝐿𝑛𝑘 = 34.34 −
𝑇

Donde k se da en s - 1 y Ten kelvins. En este diseño se desea alimentar 7850

kg de acetona por hora a un PFR. El reactor consta de un banco de 1000 tubos
cédula 40 de una pulgada.

Consideraremos tres casos:

A. CASO 1 El reactor se opera adiabáticamente.
B. CASO 2 El reactor está rodeado por un intercambiador de calor, en el
cual el coeficiente de transferencia de calor es de 100 J/m2. s. K, Y la
temperatura de intercambio de calor, Ta, es constante, en 1150 K.
C. CASO 3 El intercambiador de calor en el caso 2 tiene ahora la
temperatura de intercambio de calor variable.
La temperatura y la presión de entrada son iguales para todos los casos en
1035 K Y 162 kPa (1.6 atm), respectivamente. Grafique la conversión y la
temperatura a lo largo del reactor.

Solución

Sea A = CH3COCH3, B = CH2CO y C = CH4. Reescribiendo la reacción

simbólicamente se obtiene
AB+C
𝑑𝑋 𝑟
1. Balance molar: = −𝐹𝐴
𝑑𝑉 𝐴0

2. Ley de velocidad: −𝑟𝐴 = 𝑘𝐶𝐴

3. Estequiometria (reacción en fase gaseosa sin caída de presión):

𝐶𝐴0 (1 − 𝑋)𝑇0
−𝑟𝐴 =
(1 + 𝜀𝑋)𝑇

𝜀 = 𝑦𝐴0 𝛿 = 1(1 + 1 − 1) = 1

4. Combinando se obtiene:
𝑘𝐶𝐴0 (1 − 𝑋) 𝑇0
−𝑟𝐴 =
(1 + 𝑋) 𝑇

𝑑𝑋 −𝑟𝐴 𝑘 1 − 𝑋 𝑇0
= = ( )
𝑑𝑉 𝐹𝐴0 𝑣0 1 + 𝑋 𝑇
Para resolver esta ecuación diferencial, primero es necesario usar el balance
de energía para determinar T en función de X.

5. Balance de energía:

CASO 1. OPERACIÓN ADIABÁTICA

En el caso en que no se realiza trabajo sobre el sistema, Ws = 0, y la operación
es adiabática, 𝑄̇ = O (es decir, U = O), la ecuación (8-36) se transforma en
0 (𝑇 )
𝑑𝑇 (−𝑟𝐴 )(−[∆𝐻𝑅𝑋 𝑅 + ∆𝐶𝑃 (𝑇 − 𝑇𝑅 )])
=
𝑑𝑉 𝐹𝐴0 (∑ Θ𝑖 𝐶𝑃𝑗 + 𝑋∆𝐶𝑃 )

Como sólo entra A,

∑ Θ𝑖 𝐶𝑃𝑗 = 𝐶𝑃𝐴 )

Entonces la ecuación se transforma en

0 (𝑇 )
𝑑𝑇 (−𝑟𝐴 )(−[∆𝐻𝑅𝑋 𝑅 + ∆𝐶𝑃 (𝑇 − 𝑇𝑅 )])
=
𝑑𝑉 𝐹𝐴0 (𝐶𝑃𝐴 + 𝑋∆𝐶𝑃 )

6. Cálculo de parámetros para el balance molar por cada tubo:

7.850 𝑘𝑔/ℎ 1
𝐹𝐴0 = × 1000 𝑇𝑢𝑏𝑜𝑠 =0.135 kmol/h=0.0376 mol/s
58 𝑔/𝑚𝑜𝑙

𝑃𝐴0 162 𝑘𝑃𝑎 𝑘𝑚𝑜𝑙

𝐶𝐴0 = = = 0.0188 = 18.8 𝑚𝑜𝑙/𝑚3
𝑅𝑇 𝑘𝑃𝑎. 𝑚3 𝑚3
8.31 (1035𝐾)
𝑘𝑚𝑜𝑙. 𝐾
𝐹 5𝑚3
𝑣0 = 𝐶𝐴0 = 2.037𝑑𝑚3 /𝑠 𝑉 = 1000 𝑡𝑢𝑏𝑜𝑠 = 5𝑑𝑚3
𝐴0
7. Cálculo de parámetros para el balance de energía:
0 (𝑇 ):A
a. [∆𝐻𝑅𝑋 𝑅 298 K, los calores estándar de formación son:

0 (𝑇 )𝑎𝑐𝑒𝑡𝑜𝑛𝑎
∆𝐻𝑅𝑋 𝑅 = −216.67 𝑘𝐽/𝑚𝑜𝑙

0 (𝑇 )𝑐𝑒𝑡𝑒𝑛𝑜
𝑘𝐽
∆𝐻𝑅𝑋 𝑅 = −61.09
𝑚𝑜𝑙
0 (𝑇 )𝑚𝑒𝑡𝑎𝑛𝑜
∆𝐻𝑅𝑋 𝑅 = −74.81 𝑘𝐽/𝑚𝑜𝑙
0 (𝑇 ) 0 0 0
∆𝐻𝑅𝑋 𝑅 = 𝐻𝐵 (𝑇𝑅 ) + 𝐻𝐶 (𝑇𝑅 ) − 𝐻𝐴 (𝑇𝑅 )

= (−61.09) + (−74.81) − (−216.67)𝑘𝐽/𝑚𝑜𝑙

= 80.77𝑘𝐽/𝑚𝑜𝑙
b. ΔCp : Las capacidades caloríficas medias son:

CH3COCH3: CPA = 163 J/mol K

CH2CO: CPB = 83 J/mol K
 CH4: CPC = 71 J/mol.K

Δ𝐶𝑝 = 𝐶𝑝𝐵 + 𝐶𝑝𝐶 − 𝐶𝑝𝐴 = (83 + 71 − 163)𝐽/𝑚𝑜𝑙 · 𝐾

Δ𝐶𝑝 = −9 𝐽/𝑚𝑜]. 𝐾

Gráfica del programa Polymath. Para operación adiabática, no importa si

alimentamos o no todo a un tubo con V = 5 m3 o distribuimos el flujo a
los 1000 tubos cada uno con V =5 dm3. Los perfiles de temperatura y
conversión son idénticos, porque no hay intercambio de calor.

Tabla 1: Programa Polymath para la operación adiabática

Reporte de EDO (RKF 45)
Ecuaciones diferenciales introducidas por el usuario
[1] d(X)/d(V) = -ra/Fao
[2] a(T)/d(V) = -ra*(-deltaH)/(Fao*(Cpa+X*deICp))
Ecuaciones explícitas introducidas por el usuario
[1] Fao = .0376
[2] Cpa = 163
[3] delCp =-9
[4] Cao=18.8
[5] To = 1035
[6] deltaH = 80770+deICp*(T-298)
[7] ra = -Cao*3.5S*exp(34222*(1/To-1/T))*(1-X) *(To/T)/(1 +X)

Figura E-l: Perfiles de conversión y temperatura (a) y velocidad de

reacción (b).

Cada vez que la temperatura desciende, k disminuye; por lo tanto, la velocidad,

-rA, desciende a un valor no significativo.

CASO 2. INTERCAMBIO DE CALOR A TEMPERATURA CONSTANTE

A continuación veremos lo que ocurre a medida que agregamos calor a la
mezcla de reacción. Véase figura 2.
 Figura 2: PFR con intercambio de calor.

1. Balance molar:
𝑑𝑋 𝑟
(Paso 1) = −𝐹𝐴
𝑑𝑉 𝐴0

Usando el algoritmo: (paso 2) la ley de velocidad y el (paso 3)

estequiometria, para el caso adiabático discutido previamente en el
caso 1 (paso 4), procederemos a combinar para obtener la velocidad de
reacción como sigue:

(1−𝑋) 𝑇
(Paso 4) −𝑟𝐴 = 𝑘𝐶𝐴0 (1+𝑋) 𝑇0
5. Balance de energía:
0 (𝑇)])
𝑑𝑇 𝑈𝑎(𝑇𝑎 − 𝑇) + (−𝑟𝐴 )([−∆𝐻𝑅𝑋
=
𝑑𝑉 𝐹𝐴0 (∑ Θ𝑖 𝐶𝑃𝑗 + 𝑋∆𝐶𝑃 )

Para el sistema de reacción de acetona,

0 (𝑇)])
𝑑𝑇 𝑈𝑎(𝑇𝑎 − 𝑇) + (−𝑟𝐴 )([−∆𝐻𝑅𝑋
=
𝑑𝑉 𝐹𝐴0 (𝐶𝑃𝐴 + 𝑋∆𝐶𝑃 )
6. Evaluación de parámetros:
a) Balance molar. Por cada tubo, 𝜐𝑜 = 0.002 m3/s. La concentración
de acetona es de 18.8 mol/m3 , de modo que el flujo molar de entrada
es

𝑚𝑜𝑙 −3
𝑚3 𝑚𝑜𝑙
𝐹𝐴0 = 𝐶𝐴0 𝜐0 = (18.8 3 )(2 × 10 ) = 0.0376
𝑚 𝑠 𝑠

El valor de la k a 1035 K es de 3.58 𝑠 −1 ; en consecuencia,

tenemos

1 1
𝑘(𝑇) = 3.58 exp[34.222( − )]𝑠 −1
1035 𝑇

b) Balance de energía. Por el caso adiabático, caso l, tenemos Δ𝐶𝑝 y

𝐶𝑃𝐴 . El área de transferencia de calor por unidad de volumen de
tubería es

𝜋𝐷𝐿 4 4
𝛼= 2 = = = 150𝑚 −1
𝜋𝐷 𝐷 0.0266
( 4 )𝐿
𝑈 = 110𝐽/𝑚2 . 𝑠 · 𝐾
 Combinando el coeficiente total de transferencia de calor con el área,
se obtiene
𝑈𝑎 = 16 ,500 𝐽/𝑚3 . 𝑠 . 𝐾

Ahora utilizaremos las ecuaciones (E8-5.1) a (E8-5.6), así como las

ecuaciones (E8-5.10 Y (E8-5.11), junto con la ecuación E8-5.12, en la
tabla E8-5.3 del programa Polymath (tabla E8-5.4), para determinar los
perfiles de conversión y temperatura que se muestran en la figura E8-
5.3.

𝑉𝑓 = 0.001𝑚3 = 𝑑𝑚3

Tabla 2: Programa Polymath para un PFR con intercambio de calor

Resultados del Polyrnath

Ejemplo 8-S Producción de anhídrido acético con intercambio de calor
(Ta constante) 08-16-2004. Rcv5.1.232
Reporte ODE (RKF45)
Ecuaciones diferenciales introducidas por el usuario
[1] d{X)/d{V) = -ralFao
[2] d{T)/d{V) = (Ua'{Ta-T)+ra'deitaH)/{Fao*{Cpa+X*deICp))
Ecuaciones explicitas introducidas por el usuario
[1] Fao=.0376
[2] Cpa= 163
[3] delCp = -9
[4] Cao = lS.8
[5] To=1035
[6] deltaH = S0770+deICp'{T-29S)
[7] ra = -Cao*3.5S*exp{34222'{lIT0-11T))*{l-X)*{To/T)/{l+X)
[8] Ua = 16500
[9]Ta=1150

Figura E-2 Perfiles de temperatura y conversión en un PFR. Los perfiles

de temperatura y conversión (a) y velocidad de reacción (b) en un PFR
con intercambio de calor a temperatura constante, Ta.
CASO 3. INTERCAMBIO DE CALOR A TEMPERATURA VARIABLE
También emplearemos aire como corriente de calentamiento a co-
corriente entrando a temperatura de 1250 K y un flujo molar de (0.11
mol/s). La capacidad calorífica del aire es de 34.5 J/mol ' K

Solución
Para flujo a co-corriente
𝑑𝑇𝑎 𝑈𝑎(𝑇 − 𝑇𝑎 )
=
𝑑𝑉 𝑚̇𝐶 𝐶𝑃𝑐

𝑚𝑜𝑙 𝑗
𝑚̇ 𝐶 𝐶𝑃𝑐 = (0.11 ) (34.5 ) = 3.83 𝑗/𝑠. 𝐾
𝑠 𝑚𝑜𝑙. 𝐾

El código Polymath se modifica reemplazando Ta = 1150K en las tablas

ES 5.3 y ES-5.4 por la ecuación (ES-5 .7) , así como agregando los
valores numéricos para 𝑚̇ 𝐶 y Cp, '

Resultados del Polymath

Ejemplo 8-5 Producción de anhídrido acético con intercambio de calor
(Ta variable)
Ecuaciones diferenciales introducidas por el usuario
[1] d{X)/d{V) = -ra/Fao
[2] d{T)/d{V) = (Ua*{Ta-T)+ra*deltaH)/{Fao*{Cpa+X*deICp))
[3] d(Ta)/d(V)=Ua*(T-Ta)/mc/Cpc
Ecuaciones explicitas introducidas por el usuario
[1] Fao=.0376
[2] Cpa= 163
[3] delCp = -9
[4] Cao = l8.8
[5] To=1035
[6] deltaH = 80770+deICp'{T-298)
[7] ra = -Cao*3.5S*exp{34222'{lIT0-11T))*{l-X)*{To/T)/{l+X)
[8] Ua = 16500
[9] mc=.111
[10] Cpc=34.5

Los perfiles de temperatura y conversión se muestran en la figura E-3

.Figura E-3 (a) Perfiles de temperatura y (b) conversión en un PFR a
temperatura variable con intercambio de calor, Ta .

Programando en otro lenguaje de programación MATHCAB:

PRODUCTION OF ACETIC ANHYDRIDE
REACTOR OPERATED ADIABATICALLY

CAo  18.8 To  1035 Tr  298 Fao  38.3 Hrxo  80770

6 2
Cpa( T)  26.63  0.183 T  45.86 10 T

6 2
Cpb ( T)  20.04  0.095 T  30.95 10 T

6 2
Cpc( T)  13.39  0.077 T  18.71 10 T

  13.39  20.04  26.63

  0.077  0.095  0.183

   18.71 10 6   30.95 10 6  45.86 10 6

a  26.63 a  0.183 
a   45.86 10
6 

Cp ( T)  Cpc( T)  Cpb ( T)  Cpa( T)

k( T)  8.2 10  exp 
14 34222 
 
 T 

k( T)  CAo  ( 1  X )  To
rA ( X T) 
(1  X)T

rA ( X T)
f1( X T) 
Fao

rA ( X T)   
Hrxo    ( T  Tr) 
 
    T3  Tr3
2
 T  Tr
2

f2( X T) 
 2 3 
Fao   a  a  T  a  T  X  Cp ( T) 
2

 f1 Y0 Y1 
D( VY)   
 f2 Y0 Y1 
 

 0 
ic    V0  0 Vf  5 N  1000
 1035 

S  rkfixed ( ic V0 Vf N D)

 0  1  2
v  S x  S t  S
 PRODUCTION OF ACETIC ANHYDRIDE
REACTOR OPERATED WITH HEAT EXCHANGE

U  110 avol  150 Ta  1150

rA ( X T)
f3( X T) 
Fao
U avol ( Ta  T)  rA ( X T)  Hrxo    ( T  Tr) 
  2
 T  Tr    T3  Tr3
2

f4( X T) 
 2 3 
 2
Fao  a  a  T  a  T  X  Cp ( T) 
 f3 Y0 Y1 
D( VY)   
 f4 Y Y  
 0 1 

 0 
ic    N  100 V0  0 V1  1
 1035 

S  rkfixed ( ic V0 V1 N D)

 0  1  2
v1  S X1  S T1  S
 Resultados y discusión de resultados
Se observó en el caso (1) que para dicha reacción endotérmica adiabática, la
reacción prácticamente finaliza (se extingue) tras 3.5 m3, debido a la gran
caída de temperatura, y que se logra muy poca conversión más allá de este
punto. Una manera de aumentar la conversión sería agregar un diluyente como
el nitrógeno, que pudiera suministrar el calor sensible para esta reacción
endotérmica. Sin embargo, si se agrega demasiado diluyente, la concentración
y la velocidad serán muy bajas.
Por otra parte, si se agrega muy poco diluyente, la temperatura descenderá y
prácticamente extinguirá la reacción.
Un banco de 1000 tubos cédula 40 de una pulgada con 1.79 m de largo
corresponde a 1.0 m3 y da una conversión del 20%. El ceteno es inestable y
tiende a explotar, lo cual es un buen motivo para mantener baja la conversión.
Sin embargo, el material de la tubería y la cédula del tubo deben verificarse
para ver si son convenientes para tales temperaturas y presiones.
Se observa en el caso (2) que la temperatura del reactor pasa por un mínimo a
lo largo del reactor. En la parte delantera del reactor, la reacción se lleva a
cabo muy rápidamente, tomando energía del calor sensible del gas y
provocando que la temperatura de este gas descienda, porque el
intercambiador de calor no puede suministrar energía a velocidad igualo mayor.
Esta caída de temperatura, acoplada con el consumo de reactivos, hace más
lenta la velocidad de reacción a lo largo del reactor. Por tal velocidad de
reacción más lenta, el intercambiador de calor suministra energía a velocidad
mayor de la que la reacción retira energía de los gases; como resultado, la
temperatura aumenta.
 Bibliografía
1. Smith J. (1990) Ingeniería Cinética Química, Editorial Cecsa; México
2. Fogler H. S. (2001) Elementos de Ingeniería de las Reacciones Químicas.
Editorial Pearson Educación, cuarta edición.
3. Levenspiel O. (2004) Chemical reaction engineering, editorial John Wiley &
Sons. 3ra edición New York.
4. Castaneda R. (1999) Chemical reactor in porous media, editorial Wiley 4ta
edicion