PRODUCTO INTEGRADOR DE APRENDIZAJE:

ANÁLISIS DE VARIABLES QUE AFECTAN LA

CONVERSIÓN DE LA REACCIÓN DE
SAPONIFICACIÓN EN UN CSTR.
Hernández Alan, Hérnandez Juan, Rodríguez Debany,
Santiago David, Torres César.
Resumen
El presente documento presenta el desarrollo del producto integrador de aprendizaje para el laboratorio de
reactores químicos. Muestra el análisis de distintas variables las cuales tienen un efecto en la variable de
respuesta, en este caso la conversión. Tomando como base los trabajos realizados a lo largo del semestre,
se utilizó la cinética de la reacción de saponificación entre el hidróxido de sodio y acetato de etilo
obtenida en la sesión A y B. Se analizaron varias situaciones como el tipo de propela, la variación de
flujos totales e individuales, la altura de la propela, variación en la agitación, etc. Todos estos análisis con
la finalidad de ver que tan significativo es el efecto de las variables en la conversión del NaOH. A partir
de la reacción de saponificación con un segundo orden y utilizando los datos cinéticos obtenidos por esta
misma reacción, se obtuvo que para el diseño del reactor se debe elegir un agitador de hélice marino de
tres palas, el flujo volumétrico total que nos dio mejores resultados fue el de 200 cc/min, de los cuales 50
cc/min son para el hidróxido de sodio y 150 para el acetato de etilo, de acuerdo con las relaciones
encontradas, se obtuvo que la altura apropiada de la propela es de 9.5 cm y de acuerdo con el análisis la
agitación por elegir sería de 250 rpm.

Introducción siglas en ingles CSTR (Continuous Stirred Tank

La reacción saponificación del acetato de etilo es Reactor).
de gran importancia en la industria debido al gran Debido a las propiedades del flujo, este reactor
número de aplicaciones que tienen sus productos. presenta las siguientes características:
Es utilizado, desde la fabricación de pinturas hasta o Dado que se presenta una agitación
sistemas de almacenamiento de calor. eficiente y el fluido dentro del recipiente
esta uniformemente mezclado (por tanto,
C H 3 COOC 2 H 5 + NaOH ⟶ C H 3 COONa +C 2 H 5 OH todos los elementos de fluido están
… (1) uniformemente distribuidos), todos los
elementos del fluido tienen la misma
Este tipo de reacciones se llevan a cabo probabilidad de abandonar el reactor en
generalmente en reactores de tipo tanque agitado cualquier tiempo.
junto con una camisa de refrigeración debido a o Existe una distribución de tiempos de
que su calor de reacción es de -62.1 kJ/mol [1], residencia. Dicha distribución puede
indicando que esta corresponde a una reacción representarse matemáticamente. El DTR
exotérmica. en conjunto con un modelo y datos
Las condiciones utilizadas comúnmente para cinéticos nos permite calcular
realizar esta reacción son entre 23ºC y 40ºC. Con conversiones en la salida del reactor. [3]
un pH entre 11.5 a 12.5.[2] o Las propiedades (concentración,
El reactor donde se lleva a cabo la reacción de temperatura) dentro del reactor son
saponificación es de mezcla completa, también uniformes. Se puede observar que la
conocido como: Reactor de mezcla completa de corriente de salida también presenta las
flujo continuo, reactor de tanque agitado y por sus

1
 mismas propiedades que el fluido dentro segundo componente es longitudinal y actúa en
del recipiente. dirección paralela al eje. El tercer componente es
o Debe existir un cambio en escalón desde el tangencial o rotacional, y actúa en dirección
valor de entrada al valor de salida de tangente a la trayectoria circular alrededor del eje.
[5]
cualquier propiedad del sistema.
o La velocidad de reacción es constante
dentro del reactor.
o Dado que la densidad en el sistema de Objetivo
flujo no es necesariamente constante. o Analizar las diferentes variables que tienen
Entonces, la densidad de las corrientes una influencia en la variable de respuesta,
puede cambiar entre la entrada y la salida. la conversión.
[4]

o Cumplir con la elección correcta de los

Diseño de turbina. parámetros que afectan en la conversión y,
El encargado de diseñar un tanque agitado dispone por ende, en el diseño de un reactor CSTR.
de un gran e inusual número de opciones sobre el
tipo y localización del agitador, las proporciones
del tanque, el número y las proporciones de los Metodología
deflectores y otros factores. Cada una de estas Se analizaron diferentes incisos con la intención
decisiones afecta la velocidad de flujo del líquido, de obtener el efecto que causan estos en la
los patrones de velocidad y el consumo de conversión del NaOH en la reacción de
potencia. saponificación. A continuación, se describen
estos parámetros analizados, así como los datos
Impulsores de alta eficiencia. proporcionados para cada uno de ellos.
Se han desarrollado variantes de las turbinas de
aspas inclinadas para proporcionar un flujo axial 0.90 Conductividad vs tiempo
más uniforme y un mejor mezclado, además de 0.85

reducir la potencia requerida para determinado 0.80

flujo.
Conductividad (S/m)

0.75

Estos impulsores o agitadores son ampliamente 0.70

usados para mezclar líquidos de baja o moderada 0.65

viscosidad, pero no se recomiendan para líquidos 0.60

muy viscosos o para dispersar gases. 0.55

0.50
0.00 500.00 1000.00 1500.00 2000.00 2500.00 3000.00 3500.00 4000.00

Patrones de flujo. Tiempo (s)

La forma en que se mueve un líquido en un tanque

agitado depende de distintos factores como: el tipo a) Tipo de propelas
de impulsor, las características del líquido, Figura 1. Datos experimentales de conductividad vs tiempo de
los tipos de propelas.
especialmente su viscosidad; así como del tamaño
y las proporciones del tanque, deflectores y el
agitador. La velocidad del líquido en cualquier
punto del tanque tiene tres componentes, y el
patrón de flujo global en el mismo depende de las Propela 1
variaciones en estos tres componentes de
velocidad de un punto a otro. El primer Propela 2
Propela 3
componente de la velocidad es radial y actúa en
dirección perpendicular al eje del impulsor. El
2
 0.9 Conversión NaOH vs tiempo
0.8
b) Variación en flujos totales
0.7
0.10 Concentración de NaOH vs tiempo
0.6
0.09
0.08 300 cm3/min 200 cm3/min 0.5

ConversiónNaOH(XA)
0.07 Propela 1 Propela 2
0.4
0.06 150 cm3/min
C NAOH (M)

0.05 0.3
Propela 3
0.04
0.2
0.03
0.02 0.1
0.01
0.00
0.0
0 200 400 600 800 1000 1200 1400 1600 1800 0.00 500.00 1000.00 1500.00 2000.00
Tiempo (s) Tiempo (s)

d) Comparación de propelas
Figura 2. Datos experimentales de concentración y tiempo para
flujos volumétricos diferentes Figura 4. Datos experimentales de conversión y tiempo para
diferentes tipos de propelas.
c) Variación en los flujos individuales

Tabla 1. Cambios en los flujos volumétricos individuales para Propelas:

diferentes experimentos.
1. Agitador de turbina tipo hojas planas y
Experimento Flujo Flujo
volumétrico volumétrico disco.
individual de A individual de B
(cc/min) (cc/min) 2. Agitador de turbina de cuatro hojas
a 150 50 inclinadas.
b 100 100
c 50 150 3. Agitador de hélice marino de tres palas
A = Hidróxido de sodio, B = Acetato de etilo,
CA0-exp = 0.08 M
Conversión de NaOH vs tiempo
0.80
e) Altura de las propelas
0.70

0.60
Tabla 2. Alturas de propelas y conversiones experimentales.
0.50
Altura de propela
Conversión

0.40 Conversión (XA)

(cm)
0.30
b 13 0.8876
0.20

0.10
9 0.8000
0.00
5 0.7500
0 500 1000 1500 2000 2500 3000 3500 4000
Tiempo (s)

Figura 3. Conversiones experimentales para diferentes flujos

volumétricos individuales especificados en la tabla 1.

a
 Figura 6. Datos experimentales de concentración, conversión y
Prueba 2 a 150 rpm
Datos experimentales tiempo para un agitador de 250 rpm
0.45
0.12 1.0
0.9
0.10
0.4 0.8
g) Comparación de cambios en flujos
0.7
0.08
volumétricos individuales
Concentración NaOH (M)
Conductividad (S/m)

0.35 0.6

Conversión
0.06 0.5
concentración conversión
0.3 0.4 Figura 8. Datos experimentales de la prueba 1, dónde el flujo
0.04
0.3 volumétrico es de 150 cc/min en cada reactivo.
0.25 0.2
0.02
0.1
0.2
0.00 0.0
00 220 440 6
60 8
80 10
100 12
120 14014 16016 18018 200 20
Tiempo
Tiempo(s)(s) Figura [Link] experimentales de la prueba 2, dónde el flujo
volumétrico es de 150 cc/min para el acetato de etilo y de 75
cc/min para el hidróxido de sodio.
f) Variación en agitación
Datos experimentales a 250 rpm Figura 7. Datos experimentales de concentración, conversión y
0.12 1.0
De igual forma seunutilizaron
tiempo para datos
agitador de 350 rpm cinéticos

0.9
obtenidos en las sesiones A y B. El experimento
0.10
0.8 elegido fue el SA2 experimental. Los datos se
0.08
0.7 muestran a continuación (CA0 representa la
Concentración NaOH (M)

0.6
concentración inicial en el hidróxido de sodio,
Conversión

0.06 concentración conversión 0.5

mientras que CB0 es la concentración inicial del
acetato de etilo):
0.4
0.04
0.3

0.2
0.02
0.1

0.00 0.0
0 10 20 30 40 50 60 70 80
Tiempo (s)
Tabla 3. Parámetros cinéticos utilizados.

Figura 5. Datos experimentales de concentración, conversión y Parámetro Valor

tiempo para un agitador de 150 rpm
Volumen del reactor 2.3 L
Velocidad de reacción (k) L
0.327
Datos experimentales a 350 rpm
mol s
0.12 1.0

0.9
Temperatura 30 °C
0.10
0.8
CA0 0.1 M
0.7 CB0 0.1 M
Concentración NaOH (M)

0.08
0.6
CA 0 1
Conversión

0.06 0.5
(θ)
concentración conversión 0.4 CB0
Orden de reacción 2
0.04
0.3

0.2
0.02
0.1

0.00
0 50 100 150 200 250
0.0 Con estos datos proporcionados se procedió a
Tiempo (s) realizar cálculos para la distribución en tiempos
de residencia en todos los incisos, a excepción
del inciso “a” y “e” donde solo fueron calculados
Prueba 1
las dimensiones ideales del reactor, el cálculo de
0.75
la DTR es necesario para poder encontrar las
conversiones con los modelos de segregación y
de mezclado máximo, para de esta forma
0.7

encontrar si se trata de un macro o microfluido.

Conductividad (S/m)

0.65

0.6

0.55 4
0.5
0 5 10 15 20 25
Tiempo (s)
Cálculo de DTR Con los datos cinéticos obtenidos de la sesión
Primero se obtuvo la función de la distribución SA2 de datos experimentales se procede a
de tiempos de residencia E(t) a través de la analizar dos modelos, el totalmente segregado y
siguiente ecuación. el de mezclado máximo con el fin de obtener una
C (t ) predicción de la conversión.
E ( t )= ∞
… (2) Modelo totalmente segregado
∫ C ( t ) dt Para obtener la conversión se desarrolla la
0
ecuación 9 con lo que se obtiene la ecuación 10,
De acuerdo con la ecuación 2, C(t) es la donde Xexp(t) es la conversión experimental la
concentración, por lo que se procede a integrar cual se calcula con la ecuación 11 para una
esta concentración. En los casos donde se reacción de segundo orden y con las
proporcionaron conductividades se encontró una concentraciones de ambos reactivos iguales.
relación entre esta y la concentración, como se
muestra en las ecuaciones 3 y 4. Dónde k i dX
representa la conductividad. =X exp ( t)E (t) … (9)
dt
C i=m k i … (3)
Haciendo la sumatoria de todos los datos de
mk ( t ) conversión se obtiene la ecuación 10 integrando
E ( t )= ∞
… (4) con el método del trapecio.
m∫ k ( t ) dt
0 ∞
X =∫ X exp ( t ) E ( t ) dt … (10)
Es necesario un método para integrar 0
numéricamente la conductividad o la
concentración respecto al tiempo (k(t)) por lo 1
X ( t )=1− … (11)
que se seleccionó el método del trapecio, el cual C A 0 kt +1
se describe a continuación. Modelo de mezclado máximo
b b
Para encontrar la conversión a la salida del
I =∫ f ( x ) dx ≅∫ f I ( x ) dx reactor se aplica la siguiente fórmula (13) para
… (5)
a a una reacción de segundo orden.
∞
f (a)+ f (b)
I =( b−a ) … (6) ∫ E ( t ) dt =1−F (t) … (12)
2 t

De acuerdo con la ecuación 7 se obtiene la E (λi )

fracción del material que sale del reactor y que
ha residido en él, dentro de t1 y t2. Y de acuerdo
X i−1= X i−∆ λ
[
1−F( λi)
2
X i−k C A 0 (1−X i ) …
]
con este concepto se da la ecuación 8. (13)
t2 Donde λ es el tiempo.
∫ E ( t ) dt … (7) Conversión experimental
t1 De acuerdo con los datos proporcionados, ya
fueran concentración (14) o conductividad (15)
∞
se emplearon las siguientes ecuaciones para
∫ E ( t ) dt=1 … (8) obtener las conversiones experimentales.
0
C
X exp=1− A … (14)
Las ecuaciones 7 y 8 de igual forma son CA 0
obtenidas a través del método del trapecio. K o−K t
X exp= … (15)
K o−K ∞
5
Como se mencionó anteriormente, las ecuaciones Resultados y discusiones
anteriores fueron empleadas en el calculo de los a) Tipo de propelas
incisos, a excepción del inciso e, para este inciso Como se muestra en la figura 1, se llevó a cabo
se utilizaron relaciones encontradas en la un experimento en continuo, cambiando los tipos
referencia [5] para determinar la altura ideal de de propelas, de acuerdo con esta misma figura,
la propela, a continuación, se explican las los tipos de propelas que se utilizaron para este
variables en la figura 10 y 11. experimento son los de las figuras 12, 13 y 14.

Figura 12. Propela 1. Agitador de turbina tipo hojas planas y

disco.

Figura 10. Dimensiones de un CSTR.

Figura 11. Dimensiones de la propela de un CSTR.

Figura 13. Propela 2. Agitador de turbina de cuatro hojas
Las relaciones encontradas se presentan a inclinadas.
continuación.
Da 1
= … (16)
Dt 3

H
=1 … (17)
Dt

J 1
= … (18)
Dt 12
Figura 14. Propela 3. Agitador de hélice marino de tres palas.
E 1
= … (19)
Dt 3 A través de la ecuación 19, ya que tenemos datos
de conductividad, obtenemos la figura 15, donde
W 1 podemos notar el cambio de conversión respecto
= … (20)
Da 5 al tiempo, de igual manera se observa el
intervalo de cada una de estas propelas.
L 1
= … (21)
Dt 4

6
 a) Tipos de propelas (conversión) b)Conversiones a diferentes flujos
1 1.2 150 cc/min 200 cc/min 300 cc/min
0.9
0.8
1
Conversión NAOH

0.7
0.6

Conversión NaOH
0.8
0.5
0.4 Propela 2
0.6
0.3
0.2
Propela 1 0.4
0.1
0
-500.00 500.00 1500.00 2500.00 3500.00 0.2

Tiempo (s)
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Figura 15. Conversión a partir de los datos proporcionados en Tiempo (s)
el inciso a). Figura 16. Conversiones para diferentes flujos volumétricos de
acuerdo con el inciso b).
De acuerdo con la figura 15, podemos notar
como en la propela 1 se llega hasta Se analizaron los datos del flujo de 150 cc/min,
aproximadamente 0.58 en la conversión, para la donde se obtuvieron los siguientes resultados
propela 2 se nota un cambio en la conversión y tanto para las predicciones de las conversiones
como esta va en aumento, por último, en la para los modelos de segregación y mezclado
propela 3 la conversión va en aumento, por lo máximo y se compararon con el experimental.
que esta propela es la que ayuda mas al
mezclado en el CSTR. 150 cc/min
Los cálculos de DTR no se llevaron a cabo para 1.2 experimental segregado
este inciso debido a que los datos para las
1 mezclado maximo
propelas no estaban separados adecuadamente,
por lo que esto interferiría en la predicción de las 0.8
conversión

conversiones con los modelos de segregación y 0.6

mezclado máximo.
La propela que se elegiría sería la 3, ya que 0.4
después del estancamiento presentado por quitar 0.2
la propela 2, esta mejora la conversión debido a
0
que hay un mejor mezclado. 0 200 400 600 800 1000 1200 1400 1600 1800 2000
Tiempo (s)
b) Variación en los flujos totales
De acuerdo con los datos mostrados en la figura Figura [Link]ón de conversiones para un flujo de 150
2, se calcularon las conversiones experimentales cc/min.
a partir de la ecuación 18, con lo que se obtuvo
la figura 17. Como se observa en el grafico anterior, el
modelo de segregación es el más adecuado a los
datos analizados antes de los mil segundos,
después de este ningún modelo lo predice
adecuadamente. La conversión máxima
alcanzada es de 0.59.
Pasando al análisis de los datos de 200 cc/min
7
 c) Variación en los flujos individuales
200 cc/min
Conversión contra tiempo c)
experimental
segregado 0.80
mezclado máximo 0.60

Conversión
1.2 0.40
1
0.20
Conversión

0.8
0.6 0.00
0 500 1000 1500 2000 2500 3000 3500 4000
0.4
0.2 Tiempo (s)
0
0 0 0 0 0 0 0 0 0 0 0 Figura 20. Conversiones contra el tiempo para diferentes flujos
20 40 60 80 100 120 140 160 180 200
Tiempo (s) individuales.

Figura [Link]ón de conversiones para un flujo de 200

cc/min. De acuerdo con el gráfico obtenido de la
conversión (figura 20) podemos observar que la
mayor conversión sucede en el marcado como c
Con los datos de este flujo, se alcanza una que es donde existe un flujo para A de 50 y para
conversión máxima de 0.68, podemos observar b de 150, donde A es hidróxido de sodio y B es
que el modelo de segregado se ajusta hasta un acetato de etilo, esto se confirma con la figura
tiempo de 1400 segundos. Sin embargo, a partir 21, a pesar de que ningún modelo se ajusta
de este punto ningún modelo se ajusta. correctamente.
Ahora se calcula para los datos cuando existe un
flujo de 300 cc/min
Comparación conversiones
300 cc/min
experimental segregado
experimental mezclado máximo
segregado 1.00
mezclado máximo
0.80
1.2
Conversión

0.60
1
0.40
0.8
Conversión

0.20
0.6
0.00
0.4 0 500 1000 1500 2000 2500 3000 3500 4000 4500

0.2 Tiempo (s)

0
0 200 400 600 800 100012001400160018002000 Figura 21. Comparación de conversiones con modelos cinéticos.

Tiempo (s) d) Comparación de propelas

De acuerdo con los tres tipos de propelas
Figura [Link]ón de conversiones para un flujo de 300
cc/min. mencionas y mostrados en la figura 12, 13 y 14,
y con los datos proporcionados se calculan las
En este experimento, el modelo que mas se conversiones para cada una de ellas como se
ajusta es el mezclado máximo, además de que muestra en la figura 4. A partir de aquí
los datos experimentales son los que dan una analizamos cada una comparándolas con los
mejor conversión cercana al 1. modelos de segregación y mezclado máximo.
8
 Para obtener la altura total de la propela, es
Conversión NaOH vs tiempo necesario restar la altura (H) a la altura de las
aspas (E), dando como resultado 9.54 cm.
Propela 1
f) Variación en agitación
Propela 2
Propela 3
Se realizó una grafica comparando el tiempo y la
Mezclado máximo conversión para cada una de estas revoluciones
1.0 Segregado por minuto, obteniendo la figura 23, donde se
0.9 observa que en donde hay 250 rpm llega más
rápido a una conversión alta, por lo que esta
0.8
agitación es la que se elige.
Conversión NaOH (XA)

0.7
0.6 Variación de agitación
0.5
150 RPM 250 RPM 350 RPM
0.4
1.00
0.3
0.90
0.2
0.80
0.1
0.70
0.0
0.00 500.00 1000.00 1500.00 2000.00 Conversión 0.60

Tiempo (s) 0.50

0.40
Figura 22. Conversión a la salida del reactor.
0.30
Como se observa en el gráfico el modelo de 0.20
segregación predice la conversión a la salida del 0.10
reactor para las propelas. Sin embargo, antes de
0.00
la salida podemos observar como hay una gran 0 50 100 150 200 250
diferencia entre ambos modelos y los datos Tiempo (s)
experimentales, lo que se puede deber a la
solubilidad de los reactivos. Para solucionar esto Figura 23. Conversiones a diferentes agitaciones.
se podría cambiar la concentración de los
g) Cambio en flujos volumétricos
reactivos, de modo que la concentración de
individuales.
acetato de etilo sea mayor a la de hidróxido de
Se obtuvo una grafica de conversiones
sodio.
experimentales como se muestra en la figura 24
e) Altura de la propela
De acuerdo con las relaciones mostradas en la
metodología, se obtienen los siguientes datos:
Tabla 4. Dimensiones para el reactor CSTR

Volumen total 2.3 L

Diametro total (Dt) 0.1431 m
Altura (H) 0.1431 m
Diametro aspas (Da) 0.0477 m
Grosor de deflectores (J) 0.0119 m
Altura de aspas (E) 0.0477 m
Grosor de aspa (W) 0.0095 m
Longitud de aspa (L) 0.0358 m
9
 Cambio en las
Conversiones experimentales concentraciones.
e) Altura de la 9.54 cm
experimental 2 experimntal 1
propela
1.2
f) Variación de 250 rpm
1 la agitación
0.8 g) Flujo Se alcanza una menor
Conversión

individual conversión en el
0.6 experimento 2, por lo
0.4 que es recomendable
que se tenga mayor
0.2
concentración de
0
0 5 10 15 20 25
hidróxido de sodio.
Tiempo (s)

Figura 24. Conversiones experimentales de datos brindados.

Conclusión
Como se puede observar a través de este
Las conversiones llegan a un punto y se documento, es posible decir que se logró el
detienen, esta es la conversión máxima que se correcto análisis del comportamiento de las
logra, lo cual puede ser debido a las diferentes variables puestas a prueba mediante la
concentraciones, el experimento 1 comienza bien experimentación para así comparar su
en los primeros segundos, pero luego tarda en comportamiento respecto a la conversión
subir esa conversión. También se observa que el obtenida en cada experimento.
experimento 2 alcanza una menor conversión, y De acuerdo con el análisis realizado se concluye
esto podría ser por las concentraciones iniciales que para el diseño de un reactor donde se tenga
de los reactivos, ya que el acetato de etilo podría mejor aprovechamiento de las conversiones, es
ser el reactivo limitante. recomendable que se utilice un agitador de hélice
marino de tres palas, debe estar a una altura de
Resumen de resultados 9.54 cm y contar con una agitación de 250 rpm,
Tabla 5. Resultados elegidos para el diseño de CSTR. el reactor debe contar con un flujo de 200
cc/min, y que en este 50 cc/min sean de
Parámetro analizado Elección hidróxido de sodio y 150 cc/min sean de acetato
a) Tipo de Propela 3- Agitador de etilo, sin embargo, que exista una mayor
propela de hélice marino de 3 concentración de hidróxido de sodio, debido a
palas. que el acetato de etilo es el reactivo limitante de
b) Flujo total Flujo de 200 cc/min, la reacción.
propone una más alta Cabe resaltar que el análisis realizado de los
conversión. diferentes factores, los cuales influyen en la
c) Flujo Flujo de 50cc/min conversión de los reactivos, fue realizado de
individual para hidróxido de manera individual, por lo que sería
sodio y de 150cc/min recomendable también realizar pruebas en
para el acetato de conjunto respecto a los parámetros y
etilo dimensiones seleccionados para el reactor para
d) Comparación Ningún modelo determinar si realmente estas propiedades en
de las propelas predice conjunto favorecen la conversión en el reactor.
adecuadamente.

10
 Referencias

[1] Perry, J. H. (1949). Chemical Engineering

Handbook.
[2] Çitak, A., & Kivrak, A.
(2019). Determination of the Expression Rate of
Ethyl Acetate Hydrolysis Reaction Depending on
the Temperature [PDF]. Recuperado el 10
octubre 2020, de
[Link]
file/651595.
[3] Fogler, H. (2008). Elementos de ingeniería de
las reacciones químicas (4th ed., pp. 867-933).
Pearson Educación.
[4] Dr. Rogelio Cuevas García. (2008).
Reactores homogéneos. 12 de diciembre de
2020, de Ingeniería de reactores Sitio web:
[Link]
ncemateria_10657.pdf
[5] McCabe, W. L. (2007). Operaciones
Unitarias en Ingeniería Química. Ciudad de
México: McGraw-Hill Interamericana.

11