UNIVERSIDAD AUTÓNOMA DE

NUEVO LEÓN

Facultad de Ciencias
Químicas
Ingeniería Química

Ingeniería de las reacciones químicas

Proyecto Integrador de Aprendizaje

”Diseño de reactor enzimático de lecho empacado
para transesterificación de 2-feniletilacetato”

Facilitador: Dr. Javier Rivera de la Rosa

Pérez Rivera José Uriel 1638636

Cd. Universitaria, San Nicolás de los Garza, Nuevo León a 6

de Diciembre del 2018
RESUMEN EJECUTIVO
Se diseñó un reactor tubular de lecho empacado (PBR) para la transesterificación de 2-
feniletilacetato con una enzima (Novozym 345), tomando como base de calculo la
producción de 120000 m3 en un año. Se calculó una conversión de 0.8 con a una
temperatura de operación de 60ºC y una caída de presión de P/P0=0.89, con 540 kg de
catalizador en un reactor de 1m3.

INTRODUCCIÓN
Los agentes saborizantes y fragancias son parte importante de la industria alimentaria y
farmacéutica, con un valor en el mercado alto, y que sigue creciendo. Pero existe una
problemática con estos productos, y es que para su obtención se tienen que llevar a cabo
reacciones a altas temperaturas, utilizando como catalizadores ácidos de Lewis, los cuales
son difíciles, y por lo tanto caros, de separar. Una alternativa para este tipo de reacciones
(esterificación/transesterificación) es el uso de enzimas empacadas como catalizador,
siendo estos efectivos al dar conversiones altas con disolventes orgánicos, y sin la necesidad
de altas temperaturas.
En este trabajo se llevará a cabo el diseño de un reactor de lecho empacado para la reacción
de transesterificación de 2-feniletilacetato, catalizada por una lipasa (Novozym 435)
utilizando hexano como solvente, a partir de 2-feniletanol y acetato de vinilo.

La reacción a estudiar es la siguiente:

El 2-feniletanol esta al 33% mol con hexano como solvente y se tomó como reactivo
limitante, ya que la estequiometria de la reacción es uno a uno, por lo que se dejó en exceso
el acetato de vinilo con una relación de 1:2 respecto al reactivo limitante. Para facilidad en
este proyecto se nombrará “A” al 2-feniletanol, “B” al acetato de vinilo, “C” al 2-
feniletilacetato, “D” al acetaldehído, y “H” al solvente.

Se utilizó una cinética de primer orden, ya que según artículos consultados [1] es la que mejor
ajusta a esta reacción. Se conocen los parámetros de energía de activación (E= 27.3 kJ/mol)
y el factor pre exponencial (k0= 857.52 s-1).

𝑘 𝑚𝑜𝑙
−𝑟 ′𝐴 = 𝜌 𝐶𝐴 [=] (1)
𝑏 𝑘𝑔𝑐𝑎𝑡 −𝑚𝑖𝑛

𝐸
𝑘 = 𝑘0 𝑒 −𝑅𝑇 (2)

ESQUEMA INICIAL
Se calculó el volumen del reactor como un PBR, tomando en cuenta el gráfico de
Levenspiel[2], suponiendo una conversión a la salida del 89%, con una temperatura a la
 entrada de 25ºC, una presión de 2 bar, una relación de flujos (B) de 2, y una producción
anual de 120,000 m3 de 2-feniletilacetato. Sin embargo, con el análisis de sensibilidad de
este proyecto se decidirá cuál son las mejores condiciones a la entrada del reactor, las
variables de salida de interés, y los parámetros requeridos para la transferencia de calor. El
reactor se dimensionó tomando en cuenta 100 tubos de diámetro de 2’’ y un largo de
4.15m. En la Figura 1 se aprecia el esquema inicial, el cual muestra las variables por calcular
en rojo y las variables propuestas en azul.

X, T, P, FA, FC
0.05 m

0.007 m

b=430 kg/m3
4.15 m
3
=0.349
V=1m
W

U, A, Ta, 𝒎ሶ

FA = 117.4 kmol/h
FB = 234.8 kmol/h FA= 1.174 kmol/h
T = 20ºC V0 = 0.73 m3/h
P = 2 bar
V0 = 73 m3/h Figura 2. Diagrama de un tubo del reactor.
Figura 1. Diagrama inicial de PBR.

TABLA ESTEQUIOMETRICA
Usando la simbología explicada anteriormente se mostrará la tabla estequiométrica, la
tabla auxiliar, y la tabla de propiedades termodinámicas.

Tabla 1. Estequiometria de la reacción.

Especie Entrada Rxn’ Salida Concentraciones (v = vo)
C8H10O (A) FA0 -FA0X FA = FA0 (1 – X) CA = CA0 (1 – X)
C4H6O2 (B) FB0 -FA0X FB = FA0 (θB – X) CB = CA0 (θB – X)
C10H12O2 (C) 0 FA0X FC = FA0 (θC + X) CC = CA0 (θC + X)
C2H4O (D) 0 FA0X FD = FA0 (θD + X) CD = CA0 (θD + X)
C6H14 (H) FH0 - FH = FA0 (θI) CH = CA0 (θI)
Tabla 2. Relaciones de flujo de entrada al reactor.
Símbol Fio (kmol/h) νi yio θi
o
A 1.175 -1 0.16 1
B 2.35 -1 0.33 2
C 0 1 0 0
D 0 1 0 0
H 3.5 0 0.5 3

Tabla 3. Datos termodinámicos.

Cpi ΔCp Σθi*Cpi H°fi (298.15 K) ΔH° rx’n
Símbolo
(J/mol.K) (J/mol.K) (J/mol.K) (kJ/mol) (kJ/mol)
A 252.64 -252.64 252.64 -137.7 137.7

B 169.5 -169.5 339 -350 350

C 123.36 123.36 0 -326 -326

D 98.15 98.15 0 -136.4 -136.4

H 196 0 588 -198.7 0

Total -200.63 1179.64 25.3

Como se ve en las tablas, se tiene una reacción endotérmica, por lo que es necesario realizar
el diseño de intercambio de calor para calentar.

TRANSFERENCIA DE CALOR
Se realizo el diseño de transferencia de calor con flujos paralelos[3] utilizando vapor saturado
a 6 bar. Las ecuaciones utilizadas son las siguientes.

𝑄ሶ1 = 𝐶𝑝𝑚𝑖𝑥 𝐹𝑇0 (𝑇 − 𝑇0 ) (3)

𝐶𝑝 𝑚(𝑇−𝑇0 )
𝑄ሶ2 = 𝑎𝑐𝑒𝑟𝑜 𝑡 (4)

𝑄ሶ𝑇 = 𝑄ሶ1 + 𝑄ሶ2 (5)

∆𝑇1 −∆𝑇2
∆𝑇𝑚𝑙 = ∆𝑇1 (6)
ln( )
∆𝑇2

∆𝑇1 = 𝑇𝑎 − 𝑇0 (7)

∆𝑇2 = 𝑇𝑎 − 𝑇 (8)

𝐴 = 𝜋𝐷𝑒𝑥𝑡 𝐿 (9)

𝑄ሶ
𝑈 = 𝐴∆𝑇𝑇 (10)
𝑚𝑙
 𝑄ሶ𝑇
𝑚ሶ𝑣 = (11)
𝜆

Se hizo la suposición inicial de que el reactor opera a 50ºC, y tiene una temperatura a la
entrada del reactor de 25ºC. Debido que es una reacción que tiene solventes orgánicos no
puede operar a temperaturas superiores a 69ºC porque habría cambio de fase. Con las
ecuaciones anteriores se calculo el coeficiente global de transferencia de masa y el vapor
necesario para llevar el reactor a 50ºC. Luego se calculo el vapor necesario para realizar esta
operación.

Vapor (kg/s) U (kW / K m2) A (m2) Ta (K) P (bar)

0.43 6.46 0.8 431.95 6

ECUACIONES DE DISEÑO

𝐸
𝑑𝑋 𝑘0
= (1 − 𝑋) 𝑒 −𝑅𝑇 (12)
𝑑𝑊 𝑣0 𝜌𝑏

𝐸
𝑈𝐴 𝑘
(𝑇𝑎 −𝑇)+(1−X) 0𝑒 −𝑅𝑇 [∆𝐻𝑟𝑥𝑛 𝑜 +∆𝐶𝑝(𝑇−𝑇𝑅 )]
𝑑𝑇 𝐹𝐴0 𝑣0
=
𝑑𝑊 𝜌𝑏 (𝐶𝑝0 +∆𝐶𝑝𝑋)
(13)

𝑑𝑃 4𝑣0 𝑇𝑃0 (1−∅)(1+𝜀𝑋) 150(1−∅)𝜇 4𝑣0 𝜌0

=− [ + 1.75 ] (14)
𝑑𝑊 𝜌𝑏 𝑔𝑐 𝐷𝑃 𝜙3 𝑇0 𝑃𝜋𝐷2 𝐷𝑃 𝜋𝐷2

𝜀 = 𝑦𝐴 0 𝛿 (15)

𝑏 𝑐 𝑑
𝛿 = ∑ Υ𝑖 = 1 − 𝑎 + 𝑎 + 𝑎 (16)

ANÁLISIS DE SENSIBILIDAD
Primero se resolverá el sistema de ecuaciones simultaneo con las condiciones de entrada al
reactor propuestas en el diagrama inicial. Posteriormente se realizará un análisis de
sensibilidad y se tomará una decisión de qué condiciones son las mas convenientes para
este reactor.

Sustituyendo en las ecuaciones las condiciones propuestas.

𝑘𝐽
27300
− 𝑚𝑜𝑙
𝑑𝑋 857.52𝑠 −1 (8.314
𝑘𝐽
)𝑇
= (1 − 𝑋) 3 𝑒 𝑚𝑜𝑙−𝐾
𝑑𝑊 𝑚 𝑘𝑔
(0.0002 ) (430 3 )
𝑠 𝑚
 𝑘𝐽
𝑘𝑊 27300
𝑚𝑜𝑙
(6.46 )(0.8 𝑚2 ) 857.52𝑠 −1
− 𝑘𝐽
𝐾 − 𝑚2 (431.95 𝐾 − 𝑇) + (1 − 𝑋) 𝑒 (8.314𝑚𝑜𝑙−𝐾)𝑇 [25300
𝑘𝐽
+ (−200.63
𝑘𝐽
)(𝑇 − 298.15 𝐾)]
𝐹𝐴0 𝑚3 𝑘𝑔 𝑚𝑜𝑙 𝑚𝑜𝑙 − 𝐾
𝑑𝑇 (0.0002 ) (430 3 )
𝑠 𝑚
=
𝑑𝑊 𝑘𝑔 𝑘𝐽 𝑘𝐽
(430 3 ) (1179.64 + (−200.63 )𝑋)
𝑚 𝑚𝑜𝑙 − 𝐾 𝑚𝑜𝑙 − 𝐾

𝑏 𝑐 𝑑
𝛿 = ∑ Υ𝑖 = 1 − + + =1-1+1+1=2
𝑎 𝑎 𝑎

𝜀 = 𝑦𝐴 0 𝛿=(0.16)(2)=0.32

𝑚3 𝑚3 𝑘𝑔
𝑑𝑃 4 (0.0002 ) 𝑇(200000𝑃𝑎)(1 − 0.349))(1 + 0.32𝑋) 150(1 − 0.349)𝜇 4 (0.0002 ) (929 3 )
𝑠 𝑠 𝑚
=− [ + 1.75 ]
𝑑𝑊 𝑘𝑔 0.007 𝑚 𝜋(0.05 𝑚)2
(430 3 )(1)(0.007𝑚)(0.349)3 (293.15 𝑘)𝑃𝜋(0.05 𝑚)2
𝑚

Para resolver este sistema de ecuaciones, se lleva a cabo un método númerio (Euler),
dando una variable de paso de 0.5.

A continucación se hará el análisis de las condiciones para el reactor.

Variables de operación con parámetros propuestos

Gráfico 1. Conversión, Temperatura y Caída de presión en función del catalizador.

Como se muestra en el gráfico, podríamos obtener una conversión de 0.8 con 5.8kg de
catalizador a una temperatura de operación de 58ºC. Sin embargo, hace falta hacer un
análisis de sensibilidad para explorar otras opciones de operación.
Gráfico 2. Efecto de la temperatura de entrada al reactor en las variables de operación.

Como se muestra en el Gráfico 2, existe una fuerte correlación de la temperatura de entrada

al reactor con la conversión y la temperatura de operación, lo cual no sucede con la caída
de presión, porque como se ve no existe una diferencia entre las 3 líneas. Basándonos en
este gráfico es conveniente operar con una temperatura de entrada al reactor de 30ºC.

Gráfico 3. Efecto del diámetro de del catalizador en las variables de operación.

Para el Gráfico 3, se utilizó la temperatura de entrada seleccionada anteriormente (30ºC).
Este gráfico nos muestra algo que las ecuaciones ya nos mostraban, y es que la única
variable que depende del diámetro de la partícula es la caída de presión, sin embargo, es
importante estudiar cual es su efecto porque la caída de presión se traduce en potencia de
la bomba. En base al gráfico se escoge un diámetro de partícula de 0.008m ya que es con la
que menos cae la presión.

Gráfica 4. Efecto de la relación de flujo de A y B a la entrada del reactor sobre las variables de operació

Para el Gráfico 4 se utilizan los parámetros ya seleccionados anteriormente (T 0=30ºC y

Dp=0.008m). Como se observa, si nuestra relación de flujo es 1:1 nuestra caída de presión
es menor y nuestra conversión aumenta. La temperatura de operación permanece
prácticamente intocable. Tomando de referencia este análisis se decide optar por la
relación ya mencionada.

Una vez seleccionadas nuestros parámetros de temperatura a la entrada, diámetro de

partícula del catalizador y relación de fujos, se realiza nuevamente el gráfico para
seleccionar cuanto catalizador se debe utilizar en el lecho.

RESULTADOS Y DIAGRAMA FINAL

 Variables de operacion a condiciones
seleccionadas
1.2
1
X, T/T0, P/P0

0.8
0.6 Conversión

0.4 Temperatura

0.2 Caída de Presión

0
0 1 2 3 4 5 6 7 8 9 10
W( kg)

Gráfico 5. Variables de operación a las condiciones seleccionadas en el análisis de sensibilidad.

El gráfico 5 muestra las variables de operación a las condiciones seleccionadas en el análisis

anterior. Lo que falta ahora es seleccionar el catalizador necesario. Tomando en cuenta que
el hexano, nuestro componente mas volátil en este sistema, ebulle a 69ºC, es prudente
hacer nuestra selección de acuerdo a la temperatura de operación, así que se opta por
seleccionar nuestras variables de operación cuando la temperatura del reactor es 60ºC. Por
lo tanto nuestra conversión será de 0.8, la presión a la salida de 1.78 bar, y el catalizador
necesario para realizar esta operación es 5.4 kg en un solo tubo. El diagrama final con todas
las condiciones seleccionadas es el siguiente.
FA = 23.48 kmol/h
FB = 23.48 kmol/h
FC = 93.92 kmol/h
X=0.8 0.05 m
T = 60ºC
P = 1.78 bar
V = 62 m3/h

0.008 m

4.15 m
V = 1 m3
W = 540 kg-cat
U = 6.48 kW/m2-ºC
A = 78 m2
Ta = 158.8 ºC
𝒎ሶ = 290 kg/s

FA= 117.4 kmol/h

FB= 117.4 kmol/h
T=30ºC FA= 1.174 kmol/h
P=2 bar V0 = 0.62 m3/h
V0 = 62 m3/h
CONCLUSIÓN

Se tenían muchas restricciones para este proceso, como el de la temperatura, ya que el

objetivo era que se realizara en fase liquida para que la enzima lograra llevar la reacción de
forma eficiente. A pesar de eso se logró encontrar condiciones mejor que las propuestas
para la misma conversión, tomando en cuenta la caída de presión y cambio de la
temperatura.
REFERENCIAS
1. Sanchez, A. and Gomez, E. (2018). Kinetic study of the enzymatic synthesis of 2-
phenyl-ethyl acetate in discontinuous tank reactor. 1st ed. [ebook] Madrid, España.
Available at: http://pubs.acs.org [Accessed 1 Feb. 2018].

2. Scott Fogler, H. (2008). Elementos de ingeniería de las reacciones químicas.

3. Ciudad de México, México: PEARSON. Propiedades térmicas de algunas
aleaciones. 2018, obtenido de: [http://materias.fi.uba.ar/6731/Tablas/Tabla5.pdf].