ESCUELA DE INGENIERÍA,
ARQUITECTURA Y DISEÑO
PROGRAMA DE INGENIERÍA QUÍMICA
CURSO MODELAMIENTO Y SIMULACIÓN
Práctica de laboratorio:
Simulación de Modelos Fenomenológicos con dos ecuaciones diferenciales ordinarias acopladas
Utilizando Matlab y Python.
1. OBJETIVOS
Objetivo General: Desarrollar una simulación computacional para analizar la dinámica de un
sistema químico compuesto por dos reacciones acopladas.
Objetivos Específicos
• Identificar y formular modelos fenomenológicos que describan reacciones químicas en
función del tiempo.
• Aplicar métodos numéricos para resolver sistemas de ecuaciones diferenciales ordinarias
de primer orden.
• Interpretar el comportamiento dinámico de sistemas químicos modelados mediante
simulaciones computacionales.
• Representar gráficamente las soluciones obtenidas y analizarlas desde el punto de vista
químico e ingenieril.
2. Marco Teórico:
Velocidad de Reacción: La velocidad de una reacción química se define como el cambio de la
concentración de un reactivo o producto por unidad de tiempo. Esta puede expresarse en forma
general como:
dC
v =
dt
Donde C representa la concentración [mol/L] y t el tiempo [s]. La velocidad puede depender de
varios factores: naturaleza de los reactivos, temperatura, presión, presencia de catalizadores y
concentración.
Cinética de Reacción: La cinética química estudia cómo y a qué velocidad se producen las
reacciones químicas. Según el orden de la reacción, la dependencia de la velocidad con respecto
a la concentración de reactivos se clasifica así:
➢ Cinética de orden cero La velocidad no depende de la concentración del reactivo:
v=k
Donde k es la constante de velocidad con unidades [mol/L·s]. Ocurre, por ejemplo, cuando el
reactivo se encuentra en exceso o cuando hay saturación catalítica.
Universidad Tecnológica de Bolívar – Programa Ingeniería Química.
1
➢ Cinética de primer orden: La velocidad es proporcional a la concentración del reactivo:
v=kC
La unidad de k en este caso es [1/s]. Es típico en reacciones de descomposición o en procesos de
transferencia limitados por una sola fase.
➢ Cinética de segundo orden: La velocidad depende del cuadrado de la concentración de un
reactivo o del producto de dos concentraciones diferentes:
v = k C2
v = k C1 C2
La unidad de k será [L/mol·s]. Se aplica a reacciones bimoleculares.
Equilibrio Químico y Reacciones Reversibles: Una reacción química reversible puede proceder
en ambas direcciones: hacia los productos y hacia los reactivos. Se alcanza el equilibrio cuando
las velocidades de las reacciones directa e inversa se igualan:
A⇌B
La constante de equilibrio K se define como:
[B]eq
K=
[A]eq
El valor de K permite estimar la posición del equilibrio. Cuando K≫1, predomina la formación
de productos; si K≪1, se favorecen los reactivos.
Sistemas Dinámicos: Un sistema dinámico es aquel cuyo comportamiento cambia con el tiempo.
En Ingeniería Química, los sistemas dinámicos modelan fenómenos como reacciones, transporte,
mezclado o acumulación. Se clasifican según el número y orden de las ecuaciones diferenciales
que los describen.
➢ Sistemas de primer orden: Se describen mediante ecuaciones diferenciales de primer orden:
dy
= f(y, t)
dt
Son comunes en procesos con una sola etapa dominante, como un tanque perfectamente mezclado
o una reacción simple.
➢ Sistemas de segundo orden: Involucran derivadas de segundo orden, típicos en sistemas
donde hay aceleración o fuerzas:
d2 y dy
2
+ a + by = f(t)
dt dt
Ejemplos: dinámica de fluidos, control de procesos entre otros.
Universidad Tecnológica de Bolívar – Programa Ingeniería Química.
2
➢ Sistemas de orden superior: Combinan múltiples ecuaciones o derivadas de orden mayor a
uno. Su complejidad requiere métodos computacionales para obtener soluciones analíticas o
numéricas. Aparecen en columnas de destilación, reactores en serie o procesos acoplados
térmicos y químicos.
3. Materiales y Equipos
• Computador con sistema operativo Windows, macOS o Linux.
• Software MATLAB.
• Acceso a documentación oficial de MATLAB y MathWorks.
• Cuaderno de apuntes, bolígrafos, calculadora científica (opcional).
4. Metodología
Descripción del problema a simular:
En esta práctica de laboratorio se propone la simulación de un modelo matemático que describe
la cinética de dos reacciones químicas acopladas en un reactor. La primera reacción es A → B, la
cual es de primer orden con respecto a la especie A y está caracterizada por una constante de
velocidad k = 0.1 s⁻¹. La segunda reacción es B → C, también de primer orden respecto a B, con
una constante de velocidad k' = 0.05 s⁻¹. Las concentraciones se expresan en moles por litro
(mol/L), y las condiciones iniciales del sistema son CA(0) = 1 mol/L y CB(0) = 0 mol/L, donde
CA y CB representan las concentraciones de las especies A y B, respectivamente, al tiempo t = 0.
El objetivo de esta simulación es analizar cómo varían las concentraciones de A y B a lo largo del
tiempo mediante la solución de un sistema de dos ecuaciones diferenciales ordinarias.
A→B→C
Las ecuaciones diferenciales que modelan el comportamiento del sistema son:
dCA
= −k · CA
dt
dCB
= k · CA − k′ · CB
dt
Se resolverá este sistema de ecuaciones usando el método numérico ode45 de MATLAB, el cual
emplea un algoritmo de Runge-Kutta adaptativo de orden 4/5.
Procedimiento general de simulación:
1. Crear un archivo .m en MATLAB para definir el sistema de ecuaciones.
2. Definir los parámetros del modelo (constantes cinéticas, condiciones iniciales).
3. Implementar el solucionador ode45.
4. Generar una gráfica de CA(t) y CB(t) vs. tiempo.
5. Documentar cuidadosamente el código con comentarios por línea.
6. Completar las tablas con los resultados obtenidos.
7. Responder el cuestionario de análisis y entregar el informe.
8. En el archivo de Google Collab se presenta la simulación utilizando Python.
Universidad Tecnológica de Bolívar – Programa Ingeniería Química.
3
Figura 1: Primera Parte Código Matlab
Universidad Tecnológica de Bolívar – Programa Ingeniería Química.
4
5. Análisis:
Tabla 1. Parámetros del sistema
Parámetro Valor Unidades Descripción
K 0.1 s−1 Velocidad de A → B
'k′ 0.05 s−1 Velocidad de B → C
CA(0) 1.0 mol/L Concentración inicial de A
CB(0) 0.0 mol/L Concentración inicial de B
Tabla 2. Resultados de la simulación
Tiempo (s) CA (mol/L) CB (mol/L)
... ... ...
Gráficas requeridas:
• Gráfico de Concentración vs. Tiempo: mostrar CA(t) y CB(t) en el mismo plano.
• Ejes etiquetados: tiempo [s], concentración [mol/L].
• Título: “Evolución temporal de A y B”.
• Leyenda para cada curva.
6. Cuestionario
1. ¿Qué representan físicamente las ecuaciones diferenciales planteadas?
2. ¿Qué significado tienen los parámetros k y k′ en este sistema?
3. ¿Por qué se utiliza ode45 en esta práctica y qué ventajas ofrece?
4. ¿Qué comportamiento observas en la concentración de la especie B a lo largo del tiempo?
5. ¿Cómo se comporta el sistema si duplicas el valor de k? ¿Y si duplicas k′?
6. ¿Qué diferencias habría si se tratara de una reacción reversible?
7. ¿Qué relación tiene este modelo con procesos reales en la industria química?
8. ¿Cuál es la importancia de simular estos sistemas en lugar de resolverlos solo analíticamente?
9. Simular otro sistema dinámico que se represente por medio de mínimo dos ecuaciones
diferenciales ordinarias. Puede utilizar los ejemplos planteados en clase, en el texto de
Corripio o que encuentre en el cual se realicen balances de materia, energía o cantidad de
movimiento. Simularlo en ambos programas Matlab y Python.
7. Bibliografía
• Fogler, H. S. (2016). Elements of Chemical Reaction Engineering. 5th ed. Prentice Hall.
• Chapra, S. C. & Canale, R. P. (2015). Numerical Methods for Engineers. 7th ed. McGraw-
Hill.
• Smith, J. M., Van Ness, H. C., & Abbott, M. M. (2005). Introduction to Chemical
Engineering Thermodynamics. McGraw-Hill.
• MathWorks. (2023). MATLAB Documentation. Disponible en:
[Link]
• Seborg, D. E., Edgar, T. F., Mellichamp, D. A., & Doyle, F. J. (2010). Process Dynamics and
Control. 3rd ed. Wiley.
• Smith, C. A., & Corripio, A. B. (2006). Principles and Practice of Automatic Process Control
(3rd ed.). Wiley. ISBN: 978-0471431909.
Universidad Tecnológica de Bolívar – Programa Ingeniería Química.
5