: DINÁMICA DE PROCESOS Y CONTROL

LABORATORIO

Alba Martha, Bermúdez Tania, Casallas Vanessa, García Laura, González Angie & Ñustes Nicolás.
marthaalgu@unisabana.edu.co, taniabero@unisabana.edu.co, vanessacaco@unisabana.edu.co,
lauragaroro@unisabana.edu.co, angiegoro@unisabana.edu.co & nicolasnuca@unisabana.edu.co.

Universidad de La Sabana
Jairo Enrique Rubiano Torres

30 de agosto del 2021

ABSTRACT Actualmente, es importante para la ingeniería

química llevar a cabo procedimientos que
In the following laboratory report it is desired to permitan estudiar el comportamiento de los
analyze the dynamic behavior of a continuous reactores químicos con el fin de evaluar las mejores
stirred tank reactor in the saponification reaction, condiciones de diseño.
taking into account the concentration and
temperature profiles, to carry out this procedure the Un modelo dinámico constituye una
dynamic model is carried out using the pertinent representación abstracta de un aspecto de la
equations and they are solved with the help of realidad. En su estructura intervienen, por una
mathematical tools, on the other hand, the parte, los elementos que caracterizan la realidad
comparison is carried out between the modelizada y, por otra parte, se puede evidenciar
temperatures of 40 ° C and 30 ° C of the reaction, las relaciones existentes entre ellos. (olmos, 2020)
which were the results obtained in the laboratory, all Este tema ocupa la modelización de relaciones
the analysis and development Mathematical is done entre variables dentro de un contexto de datos en
keeping in mind the concept of green engineering. el tiempo o series temporales. Estos datos
presentan un ordenamiento natural, ya que las
Palabras clave: Modelo dinámico, reacción de
observaciones están ordenadas de acuerdo con el
saponificación, temperatura, concentración, momento del tiempo en que han sido observadas.
reactor continuo. (Econo-Bolonia, 2015)
OBJETIVO Los procesos fisicoquímicos y fenómenos que se
Analizar el comportamiento dinámico de un reactor pueden evidenciar en el modelo dinámico son:
reacción química, difusión de masa, conducción de
continuo de tanque agitado a partir de la
calor, transferencia de calor por convección,
interpretación de los perfiles de concentración y
evaporación, mezcla turbulenta, transferencia de
temperatura. masa o energía y flujo de fluidos. Para llevar a cabo
1. Definir las ecuaciones que modelan el la solución de estos problemas se hace la
comportamiento dinámico de la reacción de evaluación del problema, mediante los siguientes
aspectos: evaluar los datos empíricos de los que se
saponificación en un reactor agitado de tanque
dispone y su exactitud; y evaluar los parámetros de
continuo.
los que se dispone y su exactitud. Para construir el
2. Aplicar softwares de programación Matlab y
modelo es necesario: desarrollar las ecuaciones
Python para resolver el sistema de ecuaciones del modelo, procedentes de principios de
del modelo dinámico. conservación (serán ecuaciones
3. Comprender el efecto de la temperatura en la algebraicas/diferenciales), procedentes de
reacción de saponificación mediante la ecuaciones constitutivas (serán en general
comparación de los resultados del modelo ecuaciones algebraicas), verificar que el modelo
dinámico a 30°C y 40°C. está correctamente especificado: Análisis de los
grados de libertad y por último verificar la
INTRODUCCIÓN consistencia del modelo: chequeo de unidades y
dimensiones. (Rodríguez, 2018)

1
De antemano, también se debe tener en cuenta En este reporte, se considera la reacción de
que existen subprocedimientos importantes para saponificación de orden global de 2 y de primer
la elaboración del modelo dinámico tales como: orden respecto al acetato de etilo y al hidróxido de
establecer los volúmenes donde se van a realizar sodio.
los balances de materia, energía y momento,
determinar qué variables van a caracterizar el −𝑟 = 𝑘 𝐶 𝐴 𝐶𝐵
sistema: entradas, salidas y estados internos, Ecuación 2 Velocidad de reacción
plantear los balances de conservación: materia,
energía y momento, plantear las ecuaciones de Donde a representa el hidróxido de sodio y b será
transporte de calor, masa y momento entre los el acetato de etilo.
diferentes volúmenes, Establecer las leyes
termodinámicas, relación entre diferentes METODOLOGÍA
volúmenes donde se aplican las ecuaciones, definir La reacción de saponificación de acetato de etilo
las restricciones del proceso y definir y aplicar las
con hidróxido de sodio ocurrió en un reactor
suposiciones del modelo. (Rodríguez, 2018)
continuo de tanque agitado, donde las
Por otra parte, un modelo matemático es un tipo de concentraciones iniciales de ambos reactivos eran
modelado que se basa en la lógica matemática, de 0.1 M. Para la toma de datos experimentales de
cuyos elementos son esencialmente variables y la reacción, se hizo uso de la unidad de reacción
funcionales, y las relaciones entre ellos vienen CEM MkII del fabricante Armfield.
expresadas a través de relaciones matemáticas,
tales como: operadores lógicos, ecuaciones, El volumen total del reactor es de 2 L y el diámetro
inecuaciones, entre otros; las cuales corresponden es de 0.153 m. El volumen y presión de operación
a las relaciones que se pueden encontrar en el del reactor fue de 1 L y 0.74 atm respectivamente.
mundo real que se modelizan. (olmos, 2020) El caudal para ambos reactivos fue de 50 ml/min,
impulsados por dos bombas peristálticas, que
La saponificación consiste en la formación de sales
envían las sustancias hacia el tanque del reactor.
de sodio y potasio de ácidos carboxílicos a partir
de la hidrólisis alcalina de esteres (Vollhardt & Adicionalmente, el reactor cuenta con un sensor de
Schore, 2007). En general, la hidrólisis de un éster temperatura ubicado en la parte superior de este,
en una disolución básica lleva por nombre que permite medir la temperatura de la solución al
saponificación. El presente reporte se basa en la
interior del reactor a medida que transcurre la
reacción del acetato de etilo (CH3COOC2H5 o
reacción. Junto a la termocupla, está un
AcEt) con una base fuerte tal como el hidróxido de
sodio (NaOH) para producir acetato de sodio conductímetro, encargado de medir y enviar la
(CH3COONa) y etanol. información al software quien almacena los datos
del proceso. El seguimiento a una reacción de
La saponificación de acetato de etilo con hidróxido saponificación del acetato de etilo se logra por
de sodio es una de las reacciones más estudiadas y medio del registro de la conductividad de NaOH.
se encuentra reportada en su mayoría como una Cabe recordar que:
reacción no catalítica de segundo orden global.
“La conductimetría es un método que se utiliza para
𝐶𝐻3 𝐶𝑂𝑂𝐶2 𝐻5 + 𝑁𝑎𝑂𝐻 → 𝐶𝐻3 𝐶𝑂𝑂𝑁𝑎 + 𝐶2 𝐻5 𝑂𝐻 medir la conductividad de una disolución,
𝐴+𝐵 → 𝐶+𝐷 determinada por su carga iónica, o salina, de gran
movilidad entre dos puntos de diferente potencial.
Ecuación 1 Reacción de saponificación del acetato de etilo
La conductividad eléctrica es un fenómeno de
transporte en el cual la carga eléctrica (Electrones
La anterior reacción de saponificación del acetato
o iones) se mueve a través de un sistema.” (Laidler,
de etilo se caracteriza por ser poco exotérmica y los
Keith y John Meiser, 2002).
valores en la literatura que indican que la constante
de equilibrio para esta reacción es muy grande. En el caso particular donde se trabaja con
Adicionalmente, dado que las concentraciones de soluciones electrolíticas la conductividad se
los reactantes usadas para este reporte son bajas conoce como conductividad electrolítica cuya
(0.1 M para ambos reactivos) y las velocidades de unidad en el SI es siemens partido metro. Sin
las reacciones inversas son despreciables frente a embargo, por razones prácticas se suele utilizar la
las de reacción directa, se puede asumir que la unidad siemens partido centímetro (S/cm). La
reacción es irreversible (Muhktar, 2015).
2
conductividad electrolítica se puede ver afectada Dado que a un reactor CSTR lo que entra, menos lo
por dos factores principalmente: La concentración que sale más lo que se genera no genera
de la solución y el carácter (débil o fuerte) del acumulación al ser un reactor continuo el balance
electrolito en sí. Por consiguiente, la velocidad de de masa general se expresa como:
una reacción en donde participan iones puede ser
evaluada por mediciones de conductividad de la 𝐹𝑖𝑛 − 𝐹𝑜𝑢𝑡 + 𝑟𝐴 = 0
mezcla reaccionante en función del tiempo. Una
solución de hidróxido de sodio y acetato de etilo Escribiendo el flujo como el producto entre la
experimenta con el tiempo una disminución concentración la velocidad de alimentación v o y la
notable de la conductividad pues los iones OH-, de expresión de la velocidad de reacción:
gran movilidad iónica, son sustituidos
𝐶𝐴𝑜 𝑣𝑜 − 𝐶𝐴 𝑣𝑜 + (−𝑟𝐴 𝐶𝐴 𝐶𝐵 )𝑉 = 0
progresivamente por los iones acetato, de baja
movilidad iónica. Por lo tanto, pueden usarse las Dividiendo por la velocidad de alimentación vo
medidas de conductividad de NaOH para seguir el toda la expresión para obtener el balance en
curso de esta reacción, para finalmente obtener términos del tiempo de residencia 𝜏. Teniendo en
una expresión matemática que relacione la cuenta que no se alimenta acetato de sodio ni
concentración de NaOH en función de la
etanol al reactor, las ecuaciones para el balance de
conductividad medida.
moles finalmente se simplifican en:
Por otro lado, el equipo cuenta con un software que
𝐶𝐴𝑜 − 𝐶𝐴 + (−𝑟𝐴 𝐶𝐴 𝐶𝐵 )𝜏 = 0
contiene un diagrama mímico del equipo, que
permite controlar los caudales para cada reactivo, 𝐶𝐵𝑜 − 𝐶𝐵 + (−𝑟𝐴 𝐶𝐴 𝐶𝐵 )𝜏 = 0
temperaturas y velocidad de agitación según los
−𝐶𝑐 + (𝑟𝐴 𝐶𝐴 𝐶𝐵 )𝜏 = 0
requerimientos del usuario. Para el presente
reporte, la velocidad de operación de definió en 70 −𝐶𝐷 + (𝑟𝐴 𝐶𝐴 𝐶𝐵 )𝜏 = 0
rpm.
Ecuación 3 Ecuaciones para el balance de masa por
La reacción de saponificación en el reactor se componente de un reactor CSTR en función de t

corrió dos veces a 30°C y 40°C para analizar el

Adicionalmente, se consideró un proceso
efecto del cambio de la temperatura en la reacción.
isotérmico y no adiabático. Por ello, el balance de
Para ello, el equipo Armfield cuenta con un tanque
energía del equipo debe tenerse en cuenta para la
de agua y una resistencia que se calienta para
solución del modelo dinámico.
suministrar energía en forma de calor al reactor. La
temperatura de deseada en el reactor se define en 𝐶𝑎𝑜 − 𝐶𝑎 = 𝑘(𝐶𝑎𝑜 − 𝑉)
el software y se selecciona la opción “automático”
𝑑𝑇
para la operación del controlador. La temperatura 𝑚𝑡 𝐶𝑝𝑖 = ∑ 𝐹𝑖𝑜 (𝐻𝑖𝑜 − 𝐻𝑖 ) + 𝑉 ∑(∆𝐻𝑟𝑥 )[−𝑟𝑖 (𝑇, 𝑋𝑖 ]
𝑑𝜃
que se define actúa como un setpoint para el 𝑖 𝑖
controlador y el equipo calentara o enfriara la + 𝐴𝑡 𝑈(𝑇𝑟 − 𝑇)
resistencia según considera pertinente para
Puesto que este tipo de reactor opera en estado
mantener la temperatura optima en el reactor. 𝑑𝑇
estable entonces = 0 y por ende el balance de
𝑑𝜃
Por otro lado, mediante el planteamiento de las energía a solucionar implica:
ecuaciones o expresiones matemáticas que
0= 𝑄𝑂 𝜌𝑓 𝐶𝑝− (𝑇𝑜 − 𝑇) + 𝑉 ∑(∆𝐻𝑟𝑥 )[−𝑟𝑖 (𝑇, 𝑋𝑖 ]
modelan el sistema dinámico del proceso es
𝑖
posible analizar el comportamiento de la + 𝐴𝑡 𝑈(𝑇𝑟 − 𝑇)
temperatura y las composiciones de los reactivos y
productos en un periodo de tiempo determinado. Ecuación 4 Balance de energía para un reactor CSTR
El modelo matemático resulta de los balances
Dado que el reactor no es adiabático se debe
molares por componente y el balance de energía
calcular el coeficiente de transferencia de calor (U).
del reactor. Para ello, teniendo en cuenta las
para ello es necesario emplear correlaciones como
ecuaciones de diseño de un rector continuo de
el número de Reynolds, Nusselt y Prandtl, esto se
tanque agitado se obtienen las siguientes
debe a que la ecuación del Nusselt nos permite
expresiones:

3
obtener el cálculo de h para lograr cumplir con la 1 1 𝐿 1
= + +
ecuación del cálculo de U. 𝑈 ℎ1 𝜆 ℎ2

La ecuación del número de Nusselt se presenta a Ecuación 8 Cálculo de U

continuación, en esta ecuación es posible
RESULTADOS Y ANÁLISIS DE RESULTADOS
encontrar una relación con el número de Reynolds
y Prandtl, una vez obtenido el número de Nusselt el Con base en los resultados experimentales
cual en este caso se emplea la ecuación reportados por el software del equipo CSTR de
relacionada con un cilindro por que el equipo Armfield y el modelo dinámico que parte de las
presente en el laboratorio cuenta con esta forma. ecuaciones de balance molar y el balance de
Por lo tanto, es necesario establecer las ecuaciones energía es posible obtener los perfile de
para las demás correlaciones presentes en la temperatura y composición molar para cada una de
ecuación. El número de Nusselt es un número las especias que intervienen la reacción.
adimensional, el cual se encuentra estrechamente
relacionado con el número de Péclet y ambos se
usan para describir la relación de la energía térmica
que se convence al fluido con respecto a la energía
térmica conducida dentro del fluido
(Connor,2019).
4/5
ℎ𝐷 0.62 𝑅𝑒1/2 𝑃𝑟 1/3 𝑅𝑒 5/8
𝑁𝑢𝑐𝑖𝑙 = = 0.3 + [1 + ( )]
𝑘 [1 + (0.4/Pr )2/3 ]1/4 282000

Ecuación 5 Número de Nusselt

De acuerdo con la ecuación 1 se encuentran 2

variables “desconocidas” las cuales hacen
referencia al número de Reynolds y Prandtl; el
número de Reynolds es uno de los números Ilustración 1 Perfil de temperatura para la reacción a 30°C con
adimensionales más importantes de la mecánica Python
de fluidos y, a su vez, de los más básicos. Este
número indica la relación entre las fuerzas de Teniendo en cuenta que la rección es de carácter
inercia y las fuerzas viscosas de un fluido (Osma, exotérmico y la entalpía de la reacción es negativa,
era de esperarse que la temperatura aumente con
2020), por otro lado, el número de Prandtl se
el tiempo. De acuerdo con la ilustración 1 es
define como la relación de difusividad impulso a la
posible verificar el comportamiento creciente de la
difusividad térmica (Connor,2019). A continuación,
temperatura. No obstante, considerando que la
se presentan las ecuaciones para el cálculo de rección es poco exotérmica, el cambio de T en el
ambos números. tiempo es bajo. Con ello, la temperatura aumenta
𝜌𝑉𝐷 solo 4°C para un tiempo de operación de 175 min.
𝑅𝑒 = La tasa de calentamiento es mayor al inicio de la
𝜇
reacción dado que la cantidad de reactivos al inicio
Ecuación 6 Número de Reynolds es mayor y, por ende, la velocidad de reacción será
más rápida liberando más energía en forma de
𝜇𝐶𝑝 calor. Por el contrario, a medida que disminuye la
𝑃𝑟 =
𝑘 cantidad de NaOH y AcEt, la formación de
Ecuación 7 Número de Prandtl productos será más lenta y con ello, la energía que
se desprende será menor.
Una vez obtenidas estas correlaciones, es posible
aplicar la ecuación para el cálculo del coeficiente
total de transferencia de calor, como se mencionó
anteriormente, mediante el cálculo del número de
Nusselt se realiza el despeje para obtener el valor
para h, y con esto emplear la siguiente ecuación y
de esta forma obtener el valor de U.
4
 Diametro (m) 0,15300
Altura (m) 0,10000
Atrans (m2) 0,04807
Q (ml/min) 50,00000
Q (m3/seg) 8,33E-07
A trans. (m2) 0,01839
ρ (kg/m3) 991,00000
µ (kg/m*s) 0,00100
Vel (m/s) 0,00005
Re 6,86969
Cp (j/kg*K) 4186
K (w/m*K) 0,607
Pr 6,89621
Nu 3,26907
h (W/m2*K) 12,96946
1/U 0,24945
U (W/m2*K) 4,00888
Ilustración 2 Perfil de composición molar para la reacción a
30°C con Python Ilustración 3 Cálculo del coeficiente global de transferencia U

De acuerdo con la Ilustración 2 se observa el En la Ilustración 3, es posible encontrar el cálculo

consumo de los reactivos y la formación de los de cada una de las correlaciones para la obtención
productos con el paso del tiempo. El del coeficiente total de transferencia de calor, este
comportamiento dinámico para el consumo del valor se encuentra influenciado por el espesor y la
NaOH en la reacción será el mismo que el conductividad térmica de los materiales a través de
comportamiento del AcEt. Esto se debe a que los cuales se realiza la transferencia de calor.
ambos reactivos influyen en la misma magnitud en
Cuanto más grande el coeficiente, más fácil se
la expresión de la velocidad de la reacción,
transfiere el calor desde su fuente hacia el producto
recordando que se considera de primer orden
respecto a cada uno de los reactivos y con una que está siendo calentado. En este caso este
relación estequiométrica 1:1. coeficiente se tiene en cuenta para ambas
variaciones de temperatura (30°C y 40°C), esto se
Por otro lado, partiendo del análisis de la solución debe a que cada una de las correlaciones no
del modelo dinámico para el reactor continuo de presenta una relación con la temperatura del
tanque agitado y la velocidad de reacción, la proceso, estas se enfocan más en propiedades
composición molar de los reactivos disminuirá con propias del fluido o del material con el cual se
el tiempo hasta consumirse casi en su totalidad. Si realiza la transferencia de calor.
embargo, conociendo la naturaleza de la reacción,
donde las reacciones que intervienen son Los códigos de Python empleados en la solución
equilibrios químicos, la composición molar de los de los modelos matemáticos del balance tanto de
reactivos y productos después de un largo periodo masa como energía son:
de tiempo se mantendrán constantes.
Figura 1 Códigos de Python para el balance de materia por
componente y de energía para el reactor CSTR

import numpy as np
from scipy.optimize import fsolve
import matplotlib.pyplot as plt

# Constantes

Cao = 0.1
Cbo = 0.1
Cco = 0
Cdo = 0
To=295

5
# ARREGLOS
y = Ca.append(respuesta[0])
Ca = [] # Hidróxido de sodio m = Cb.append(respuesta[1])
Cb = [] # Acetato de etilo z = Cc.append(respuesta[2])
Cc = [] # Acetato de sodio s = Cd.append(respuesta[3])
Cd = [] # Etanol l = T.append(respuesta[4])
T=[]
tao = np.linspace(0, 180, 50) print(Ca)
print(Cb)
def f(w, t): print(Cc)
Ca = w[0] print(Cd)
Cb = w[1] print(T)
Cc = w[2]
Cd = w[3] plt.plot(tao, Ca, 'r--', label=r'NaOH')
T = w[4] plt.plot(tao, Cb, 'b-', label=r'Acetato de etilo')
plt.plot(tao, Cc, 'y-', label=r'Acetato de sodio')
k0 = 2.0831 # L/mol*min plt.plot(tao, Cd, 'm--', label=r'Etanol')
V = 1 # Litros plt.ylabel('Concentración (mol/L)')
Fi = 0.16 # L/min plt.plot(tao, T, 'm--', label=r'Temperatura')
Ti = 295 # K plt.ylabel('Temperatura (K)')
Ea = 11.3 # kcal/mol*K plt.xlabel('Tiempo (min)')
R = 1987 # kcal/mol*K plt.legend(loc='best')
Cai = 0.1 plt.title("Perfil de temperatura del reactor CSTR")
rho = 970 # constante plt.title("Perfil de concentraciones del reactor
CpA = 0.013384 # kcal/mol*k CSTR")
CpB = 0.04038 # j/mol*K plt.show()
deltaH = -210.47 #kcal/mol
At=0.048
U=4.13 En cuanto a los resultados experimentales de los
datos obtenidos por medio del software Armfield,
k = k0 * np.exp((Ea / R) * ((1 / Ti) - (1 / T))) estos pueden ser convertido a datos relacionados
# Balance a las concentraciones gracias a las curvas de
calibración del conductímetro empleado,
f1 = Cao - Ca - (k * Ca * Cb * t) adicionalmente de acuerdo con las practicas del
f2 = Cbo - Cb - (k * Ca * Cb * t) manual se para calcular las concentraciones
f3 = -Cc + (k * Ca * Cb * t) asociada a las medidas de conductividad debe
f4 = -Cd + (k * Ca * Cb * t) llevarse a cabo una corrección de temperatura que
f5= (Fi*rho*((CpA+CpB)/2)*(Ti-T))+(V*deltaH*(- afecta las medidas cuyas relaciones se muestran
k*Ca*Cb*t))+(At*U*(Ti-T)) como:

return [f1, f2, f3, f4, f5]

# x0 = [0, 0] Ilustración 4 Relaciones para la corrección de la conductividad

# Valor inicial de acuerdo con la temperatura
x0 = [0.1, 0.1, 0, 0, 295] Posteriormente dado que se toma al hidróxido de
sodio como el reactivo limitante la concentración
for i in tao: # tao de este y conversión a lo largo del tiempo y a
medida que se varia la temperatura se expresa
# print(i) como:
respuesta = fsolve(f, x0, args=(i,))
x0 = respuesta (𝐾 − 𝐾∞ )
𝐶𝑁𝑎𝑂𝐻 = 𝐶𝑁𝑎𝑂𝐻𝑜
print(respuesta) (𝐾𝑂 − 𝐾∞ )
6
 𝐶𝑁𝑎𝑂𝐻𝑜 − 𝐶𝑁𝑎𝑂𝐻 del reactor, el controlador manda una señal para
𝑋𝑁𝑎𝑂𝐻 =
𝐶𝑁𝑎𝑂𝐻𝑜 cambiar la temperatura de la resistencia y con ellos
busca corregir los valores medidos. La resistencia
Para el otro reactivo, el acetato de etilo: cambia su temperatura varias veces durante el
proceso hasta que se logra mantener la
𝐶𝐴𝑐𝐸𝑡 = 𝐶𝐴𝑐𝐸𝑡𝑜 − (𝐶𝑁𝑎𝑂𝐻𝑜 − 𝐶𝑁𝑎𝑂𝐻 )
temperatura desea. Entonces, al inicio del proceso,
𝐶𝐴𝑐𝐸𝑡𝑜 − 𝐶𝐴𝑐𝐸𝑡 los cambios en la temperatura serán más drásticos
𝑋𝐴𝑐𝐸𝑡 = y a medida que pasa el tiempo, el equipo logra
𝐶𝐴𝑐𝐸𝑡𝑜
controlar la variable hasta conseguir un estado
Ecuación 9 Set de ecuaciones para la concentraciones y estacionario, donde la temperatura del reactor será
conversiones de los reactivos en el reactor CSTR experimental
constante.

Ilustración 5 Perfil de composición molar de los resultados

experimentales a 30°C Ilustración 7 Perfil de composición molar de los resultados
experimentales a 40°C

En la ilustración 5 se debe tener en cuenta que la

cantidad de datos tomados para el presente
estudio fueron significativos cuantitativamente
hablando, las mínimas diferencias en la
concentración que se obtuvo transformando la
lectura original del sistema en conductividad a
concentración molar puede explicar la naturaleza
vibratoria de la gráfica. En este orden de ideas, se
presume detallar que la línea de color gris
determina la formación de etanoato de sodio y
etanol a una tasa equitativa. En este orden de ideas,
Ilustración 6 Perfil de temperatura de los resultados
experimentales a 30°C
se puede comprobar el mismo comportamiento
entre el hidróxido de sodio y el etanoato de etilo
Se consideró un proceso isotérmico, donde la que se ven representados en esta ocasión por la
temperatura deseada definida en el software era línea naranja.
de 30°C. Según la ilustración 4 es posible observar
un comportamiento oscilatorio en los datos de la
temperatura. Esto se debe al tiempo que tarda el
sistema en dar una respuesta o ajustarse a las
condiciones definidas por el operador. Una vez es
definido el setpoint, el controlador calienta o
apaga la resistencia con el fin de calentar el agua
del tanque e impulsarla por una bomba hacia el
serpentín al interior del reactor. Según los datos
registrados podría afirmarse que se trata de un
sistema de control feedback, quien actúa a partir
de un proceso que ocurrió. Es decir, que según la
temperatura que mide la termocupla en el interior

7
 Ilustración 8 Perfil de temperatura de los resultados
experimentales a 40°C

Por otra parte, en la ilustración 6 se puede observar

el sistema térmico del reactor CSTR después de
haber fijado el set-point en 40°C, igualmente se
puede establecer que funciona en modalidad
feedback dado que las señales son generadas
principalmente por los cambios en la temperatura,
ocasionados por el paso del fluido de servicio que
circula en el interior del reactor. A su vez, se puede
explicar que la línea pasa de un estado de turbidez
a un estado constante a medida que se aproxima a
la temperatura objetivo.

Ilustración 10 Perfil de temperatura para la reacción a 40°C con

Python

CONCLUSIONES.

La reacción de saponificación del acetato de etilo

con hidróxido de sodio para producir acetato de
sodio se caracteriza por ser de carácter exotérmico
es decir, libera energía para poder favorecer la
producción del acetato de sodio como producto
de interés, por ende el aplicar perturbaciones
sobre el sistema asociada al suministro de energía
por medio de un sistema de calentamiento permite
evidenciar que se alcanzan menores
Ilustración 9 Perfil de composición molar de los resultados concentraciones de productos cuanto mayor es la
experimentales a 40°C
temperatura de operación.
Según el desarrollo experimental y el modelo Los modelos matemáticos del balance de energía
matemático planteado en el software de y materia desarrollados y solucionados por medio
programación Python, se puede establecer que la de una herramienta numérica como Python
formación de etanoato de sodio y etanol se permiten comparar los resultados ideales frente a
comporta de manera equimolar a partir de la los reales obtenidos de los datos registrados por
reacción de saponificación que se lleva a cabo de medio del software del reactor los cuales reflejan la
enfrentar etanoato de etilo e hidróxido de sodio, en realidad del comportamiento del sistema teniendo
este orden de ideas, se pretende enunciar que la en cuenta las perturbaciones que pueden
tasa de consumo de reactivos es equitativa con la presentarse como la velocidad y lectura de
tasa de producción de los compuestos bombeo de los reactivos, la velocidad de agitación,
mencionados anteriormente. la medida del conductímetro y la eficiencia del
En este orden de ideas, se puede analizar de los calentamiento empleado.
datos obtenidos a las diferentes temperaturas que
el comportamiento de la reacción es muy similar,
de este modo también es correcto afirmar que el REFERENCIAS.
aumento de temperatura genera producción de
etanoato de sodio, aunque en menor proporción. Laidler, Keith y John Meiser. Fisicoquímica. México:
Editorial CECSA, 2002. 987pp.

Duran, M. (2017). Cinética en la hidrólisis del

acetato de etilo midiendo conductividad. Informe
de Informática Aplicada a la Química. Recuperado

8
de:
http://www.mikeciencia.com/practicas/practica-
02.pdf

Rodríguez, M. (2018). Modelado y simulación en

Ingeniería Química. Retrieved 19 September 2021,
from
http://www.diquima.upm.es/old_diquima/docencia/
modysim2000/docs/tema2_fe

Olmos. (2020). INTRODUCCIÓN A LOS MODELOS

DINÁMICOS. Retrieved 19 September 2021, from
https://www.uv.es/olmos/Ecuaciones%20diferenciale
s.pdf

Econo-Bolonia. (2015). Modelos Din´amicos.

Retrieved 19 September 2021, from
https://ocw.ehu.eus/file.php/23/MDINAMICOS.pdf

Connor, N. (2019). ¿Qué es el número de Nusselt?

Recuperado de https://www.thermal-
engineering.org/es/que-es-el-numero-de-nusselt-
definicion/

Osma, M. (2020). ¿Qué es el número de Reynolds?

Recuperado de
https://www.aerodinamicaf1.com/2020/10/que-es-
el-numero-de-reynolds/

Connor, N. (2019). ¿Qué es el número de Prandtl?

Recuperado de https://www.thermal-
engineering.org/es/que-es-el-numero-de-prandtl-
definicion/

TLV. (s.f). Coeficiente total de transferencia de calor.

Recuperado de
https://www.tlv.com/global/LA/steam-theory/overall-
heat-transfer-coefficient.html

Çengel, Y. A., & Ghajar, A. J. (2011). Transferencia de

calor y masa: fundamentos y aplicaciones, cuarta
edición. McGraw-Hill Interamericana.