Investigacion Formativa Modelamiento

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA
FACULTAD DE INGENIERÍA DE PROCESOS
ESCUELA PROFESIONAL DE INGENIERÍA QUÍMICA
“MODELAMIENTO Y SIMULACIÓN DE UN REACTOR FLUJO PISTÓN EN

ESTADO TRANSITORIO CON INTERCAMBIO DE CALOR”
PRESENTADO POR:
 ALVA/LARREA, GILMAR FAVIAN
CARLO
 HUAMANI/HUACHACA, REYNALDO
 IBARRA/TINTAYA, HEIDY
 ROJAS/VARGAS, ANDREA VANESSA
 VILCHEZ/MOLINA, MARIA DEL PILAR
AREQUIPA – 2019
PLAN DE INVESTIGACIÓN FORMATIVA
1. RESPONSABLES DEL PROYECTO

2. PROBLEMÁTICA
Se propone el desarrollo de un programa para la simulación de reactores flujo pistón
en estado transitorio con intercambiador de calor. El modelo matemático se obtendrá a
partir de las leyes de conservación de materia, energía y momento lineal, pretende
relacionar las variables y parámetros relacionados con el desempeño del reactor; como
condiciones de alimentación, parámetros de diseño y parámetros ajustables del
proceso, lo que da lugar a modelos de parámetro distribuido, donde las variables son
influenciadas por el tiempo y la posición. La solución matemática de este tipo de
modelos es compleja y puede ser particular para cada caso, por ello se requiere utilizar
los métodos numéricos y la termodinámica de mezclas, esperando que tales métodos
sean los adecuados para poder desarrollar nuestro programa.
3. FUNDAMENTOS DEL PROYECTO DE INVESTIGACIÓN

3.1 DEFINICIÓN DEL PROBLEMA
La necesidad de un modelo dinámico es especialmente fuerte en procesos de
arranque, finalización o cuando se simulan procesos con perturbaciones, también es
importante en estudios de seguridad industrial y control de procesos.
Con el desarrollo de los métodos matemáticos y el auxilio de la computación, la
modelación matemática de los procesos adquiere gran importancia, ya que a través de
ella se puede obtener información del comportamiento del proceso sin necesidad de
realizar experimentos, que podrían ser costosos y sujetos a errores experimentales
difíciles de estimar.
La simulación surge de una necesidad de sintetizar procesos, logrando integrar
balances de materia, energía, dimensiones y costos, y así lograr una evaluación
económica preliminar, optimizando los mismos.
3.2 OBJETIVOS
3.2.1 OBJETIVO GENERAL
- Modelamiento y simulación de un reactor tipo Pistón en estado no
estacionario con intercambio de caló.
3.2.2 OBJETIVOS ESPECÍFICOS

- Formular Modelos Matemáticos para la resolución de la problemática.
- Generar los perfiles de temperatura, caída de presión y conversión
respecto al tiempo.
- Comparar los resultados con diferentes simulaciones utilizando softwares
comerciales.
3.3 HIPÓTESIS
A partir del modelamiento de un reactor flujo pistón en estado estacionario
reconociendo nuestro sistema haciendo balance hallamos la función de transferencia
global, al reconocer las variables dependientes e independientes haremos la simulación
usando el software MATLAB en el cual nosotros podemos generar diferentes perfiles.
En este caso el objeto de estudio es la hidrolisis de la sacarosa. Generamos entonces
perfiles de temperatura, caída de presión y conversión con respecto al tiempo.
3.4 JUSTIFICACIÓN
3.4.1 ALCANCES
La simulación de procesos se ha convertido en una herramienta imprescindible para
los ingenieros, ya que por medio de esta es posible integrar balances de materia,
energía, dimensiones y costos, y así lograr una evaluación económica preliminar,
optimizando los mismos. Por ello se hace necesario desarrollen modelos más rigurosos
que representen de manera adecuada la fenomenología del proceso.
Los programas de simulación se han implementado de forma generalizada en la
industria química, llegando a convertirse en herramienta indispensable para el diseño,
modificación y evaluación de equipos y plantas, dada la gran versatilidad que ofrecen
en cuanto a la manipulación de variables involucradas en los procesos, y el ahorro
económico que proporcionan al permitir realizar el análisis de diversas alternativas
para la mejora de un proceso, sin necesidad de materializar los cambios propuestos
antes de obtener conclusiones confiables acerca de la conveniencia de implementar o
no, la alternativa. Durante años, la simulación ha ganado terreno dentro del diseño y
análisis de procesos, al emplear algunos de los simuladores comerciales que se
encuentran en el mercado, para verificar la operabilidad, controlabilidad y el
desempeño ante los cambios de carga de un proceso potencial.
La mayoría de las industrias de procesos dentro de sus plantas poseen reactores
químicos. Los reactores son el corazón de la planta, ya que en estos ocurren los
principales procesos de transformación de la materia prima y por lo tanto, un
funcionamiento óptimo y eficiente de éste garantiza en buena parte el éxito de un
proceso.
Una de las formas más de procesamiento de materias primas en la industria es llevada
a cabo en reactores tubulares, sin embargo estos reactores presentan ciertas desventajas
provenientes de su comportamiento dinámico ya que por su forma geométrica y su tipo
de flujo pueden presentar puntos calientes que generan daños en los equipos y
dificultades en el mantenimiento y control de la calidad final del producto deseado, por
tanto se hace necesario el conocimiento vasto de la dinámica del proceso para poder
tomar medidas útiles en caso de inconvenientes en planta en este contexto se pretende
modelar un reactor de flujo tapón en estado no estacionario.
3.4.2 RESTRICCIONES
4. MARCO TEORICO
Un reactor de flujo de pistón es un reactor ideal ya que en el citado reactor se
considera que todos los elementos del fluido se mueven con velocidad uniforme a lo
largo de líneas paralelas. Este tipo de reactores opera en continuo, por lo que las
variables del proceso en un punto determinado del reactor tienen siempre el mismo
valor, cambiando longitudinalmente de un punto a otro del reactor.
Sin embargo, en la práctica muchos sistemas no cumplen las condiciones de
idealidad. Entre las causas principales de la desviación del comportamiento ideal de
un reactor destacan las siguientes:
- Difusión axial y radial en reactores tubulares.
- Retro-mezcla en reactores tubulares.
- Presencia de zonas estáticas y de muy baja velocidad dentro del sistema
de reacción. Este punto es especialmente importante en reactores de
tanque agitado.
- Canalización.
- Desviación del flujo del fluido respecto del comportamiento ideal.
La desviación que experimenta un reactor con respecto al comportamiento ideal se
puede estudiar mediante la determinación de la curva de distribución de tiempos de
residencia. Los métodos experimentales más habituales para llevar a cabo este
estudio están basados en la adición de una sustancia inerte, o trazador, al sistema de
reacción y en la medida a la salida del mismo de la concentración de esa sustancia
inerte con el tiempo.
Dos de las técnicas más comunes en el estudio del comportamiento de un reactor son
las siguientes:
- Test de impulso. Esta técnica experimental consiste en añadir una
determinada cantidad de sustancia inerte, o trazador, en un tiempo muy
corto, a la entrada del reactor, midiendo a continuación la variación de
la concentración de sustancia inerte con el tiempo a la salida del reactor.
- Test de escalón. En esta técnica se añade una cantidad constante de
inerte a la entrada del sistema de reacción, midiendo la variación de la
concentración de inerte con el tiempo a la salida del sistema. En este
caso particular, se denomina curva F, o curva escalón normalizada, a la
curva representativa de la concentración del trazador a la salida del
reactor, dividida por la concentración a la entrada, frente al tiempo,
Consideraremos un reactor tubular en el que se añade o elimina calor a través de las

paredes cilíndricas del reactor. Al modelar el reactor supondremos que no hay
gradiente radial en el reactor y que el transporte de calor a través de la pared por
unidad de volumen del reactor es como se muestra en la figura.
Figura 1. Reactor tubular con ganancia o pérdida de calor.
Recordando la ecuación e ignorando cualquier trabajo efectuado sobre el fluido de

reacción, obtenemos:
Diferenciamos respecto al volumen V y reunimos términos para llegar a:
Recordando que –ra=FA0(dX/dV) , y sustituyendo la ecuación por (dQ/dV),

podemos reacomodar la ecuación para obtener:
La ecuación diferencial que describe el cambio de temperatura con el volumen de

reactor (es decir, distancia reactor abajo),
Se debe acoplar con el balance de moles:

ECUACIONES DE TRANSFERENCIA DE MASA
Se toma que en el reactor se está realizando una reacción en fase gas:
Se define la velocidad de reacción:
1.
La ecuación para un reactor tubular de flujo pistón, la ecuación de diseño es la

siguiente:
2.
La ecuación de diseño expresada en volumen de reactor:
Donde es la cantidad de A que va desapareciendo por reacción de unidad de
volumen y de tiempo, que se va a relacionar con por medio de la densidad global

del lecho:
La conversión máxima va a ser fijada por condicionamiento del operador.

ECUACIONCES DE TRANSFERENCIA DE CALOR
Partiendo de la hipótesis de flujo pistón, el balance de calor para un elemento
incremental de reactor se realiza de el mismo modo a los reactores homogéneos:
3.
Dando valores:
4.
Donde se puede relacionar con el volumen y el peso del catalizador:

Sustituyendo el valor de de la ecuación de diseño (2):
5.
Llevando este valor a la ecuación de balance de calor (4) y despejando queda:
6.
Ecuación que da la variación de temperatura en el interior del reactor en función del

avance de la conversión y en la cual tenemos que:
CALCULO DEL COEFICIENTE GLOBAL DE TRANSFERENCIA DE

CALOR
Hasta este punto se ha realizado del diseño del reactor partiendo de un supuesto valor
de U. Una fórmula para el cálculo del verdadero valor de U en tubos metálicos de poco
espesor es:
El cálculo del factor de ensuciamiento ( ) es problemático por lo cual se recurre a
diferentes tablas en las cuales dan diferentes valores de , según el fluido circulante.
CAUDAL DE FLUIDO EXTERNO

Una vez que ya han sido calculadas y confirmadas todas las variables de diseño, solo
queda calcular la cantidad de fluido externo a utilizar. Se hace en función al calor
transferido, que ya ha sido calculado anteriormente en el balance de calor. Si el fluido
solo realiza un cambio de calor sensible, seria:
Donde es la cantidad de fluido externo de capacidad calorífica que experimenta
una variación de temperatura , que anteriormente se había considerado
despreciable.
Y si existe cambio de fase, entonces:
4.1 MODELOS MATEMÁTICOS

Dado que las ecuaciones que se plantearon anteriormente son EDOs, donde la
variable independiente es el volumen, se resolverán por el método numérico de
RUNGE KUTTA de cuarto orden.
MÉTODO DE RUNGE-KUTTA DE CUARTO ORDEN PARA SISTEMAS

DE DOS ECUACIONES DIFERENCIALES ORDINARIAS.
Para aproximar la solución al (PVI)

X´ = F(V,X,T)
T´ = G(V,X,T)
X(Vo) = Xo , X(Vf) = ?
T(Vo) = To , T(Vf) = ?
Proporcionar las funciones F(V,X,T) y G(V,X,T).

PASO 1: hacer H = (Vf - Vo )/N
PASO 2: Hacer I = 1
PASO 3: Mientras I <=N repetir los pasos 4 a 15
PASO 4: Hacer K11 = F(Vi,Xi,Ti)
PASO 5: Hacer K12 = G(Vi,Xi,Ti)
PASO 6: Hacer K21 = F(Vi + h/2 , Xi + h/2*K11, Ti + h/2*K12)
PASO 7: Hacer K22 = G(Vi + h/2 , Xi + h/2*K11, Ti + h/2*K12)
PASO 8: Hacer K31 = F(Vi + h/2 , Xi + h/2*K21, Ti + h/2*K22)
PASO 9: Hacer K32 = G (Vi + h/2, Xi + h/2*K21, Ti + h/2*K22)
PASO 10: Hacer K41 = F( Vi + h , Xi +h*K31, Ti + h*K32)
PASO 11: Hacer K42 = G(Vi + h , Xi +h*K31, Ti + h*K32)
PASO 12: Hacer Xi+1 = Xi + h/6*(K11 + 2*K21 + 2K31 + K41)
PASO 13: Hacer Ti+1 = Ti + h/6*(K12 + 2*K22 + 2*K32 + K42)
PASO 14: Hacer Vo = Vo + h
PASO 15: Hacer I = I + 1
PASO 16: Imprimir Xi, Ti y terminar.
5. BIBLIOGRAFÍA
- H. Scott Fogler (2001). Elementos de ingeniería de las reacciones
químicas. 3ra Ed. Pearson Educación: México. pp 458-459
- Henao., Carlos Andrés (2006). Simulación y Evaluación de los

procesos químicos, Editorial Universidad Bolivariana. Medellín.
- Daniel F. Rodríguez, Brenda M. Restrepo, Jorge A. Velásquez (2009).

Simulador dinámico de reactores flujo tapón, Editorial Universidad
Pontificia Bolivariana. Medellín.

Investigacion Formativa Modelamiento

Cargado por

Copyright:

Formatos disponibles

Investigacion Formativa Modelamiento

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Investigacion Formativa Modelamiento

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE SAN AGUSTÍN DE AREQUIPA

FACULTAD DE INGENIERÍA DE PROCESOS

ESCUELA PROFESIONAL DE INGENIERÍA QUÍMICA

“MODELAMIENTO Y SIMULACIÓN DE UN REACTOR FLUJO PISTÓN EN

1. RESPONSABLES DEL PROYECTO

3. FUNDAMENTOS DEL PROYECTO DE INVESTIGACIÓN

3.2.2 OBJETIVOS ESPECÍFICOS

Consideraremos un reactor tubular en el que se añade o elimina calor a través de las

Figura 1. Reactor tubular con ganancia o pérdida de calor.

Recordando la ecuación e ignorando cualquier trabajo efectuado sobre el fluido de

Diferenciamos respecto al volumen V y reunimos términos para llegar a:

Recordando que –ra=FA0(dX/dV) , y sustituyendo la ecuación por (dQ/dV),

La ecuación diferencial que describe el cambio de temperatura con el volumen de

Se debe acoplar con el balance de moles:

Se define la velocidad de reacción:

La ecuación para un reactor tubular de flujo pistón, la ecuación de diseño es la

La ecuación de diseño expresada en volumen de reactor:

Donde es la cantidad de A que va desapareciendo por reacción de unidad de

volumen y de tiempo, que se va a relacionar con por medio de la densidad global

La conversión máxima va a ser fijada por condicionamiento del operador.

Donde se puede relacionar con el volumen y el peso del catalizador:

Llevando este valor a la ecuación de balance de calor (4) y despejando queda:

Ecuación que da la variación de temperatura en el interior del reactor en función del

CALCULO DEL COEFICIENTE GLOBAL DE TRANSFERENCIA DE

El cálculo del factor de ensuciamiento ( ) es problemático por lo cual se recurre a

CAUDAL DE FLUIDO EXTERNO

una variación de temperatura , que anteriormente se había considerado

4.1 MODELOS MATEMÁTICOS

MÉTODO DE RUNGE-KUTTA DE CUARTO ORDEN PARA SISTEMAS

Para aproximar la solución al (PVI)

Proporcionar las funciones F(V,X,T) y G(V,X,T).

- Henao., Carlos Andrés (2006). Simulación y Evaluación de los

- Daniel F. Rodríguez, Brenda M. Restrepo, Jorge A. Velásquez (2009).

También podría gustarte