UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA QUÍMICA

ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA QUÍMICA

PROYECTO DE TESIS

I. GENERALIDADES

1. Título:

MODELAMIENTO Y SIMULACIÓN DE UN REBOILED ABSORBER

DE MULTICOMPONENTES EN ESTADO ESTACIONARIO POR EL
MÉTODO DE CORRECCIONES SIMULTÁNEAS EN LA UNIDAD DE
FRACCIONAMIENTO DE CRAQUEO CATALÍTICO

2. Tipo de investigación

Investigación Básica Desarrollo Tecnológico

Investigación Aplicada Innovación

3. NOMBRE DEL TESISTA(S)

NOMBRE DEL TESISTA 1 Robert Eduardo Graus Peláez

Código o número matrícula: 1051601211

Teléfono: 968 685 806
Email: [email protected]
Dirección: Campiña de Moche Lt:13. Curva del Sun.
NOMBRE DEL TESISTA 2 Jorge Pablo Escobar Esquivel

Código o número matrícula: 1051602111

Teléfono: 954119657
Email: [email protected]
Dirección: Av. Bolivia Mz 05 Lt 06
 4. NOMBRE DEL ASESOR
4.1. ASESOR PRINCIPAL

Apellidos y Nombres: Moreno Eustaquio Walter

Categoría y Dedicación Docente: Profesor Auxiliar-Tiempo Completo

Código: 3309
Título Profesional: Ingeniero Químico
Grado Académico: Maestro en Ciencias Económicas
Departamento Académico: Ingeniería Química
Teléfono: 044 280431 - 969819999
Email: [email protected]
Dirección: Calle Isabel de Bobadilla Mz C2 Lote 2, Monserrate
V Etapa

5. LUGAR DE EJECUCIÓN DEL PROYECTO

- Laboratorio de Computación e Informática (LABCI)

- Facultad de Ingeniería Química
- Universidad Nacional de Trujillo

6. PRESUPUESTO

Presupuesto requerido en nuevos soles S/. 5795,00

7. DURACIÓN DEL PROYECTO

INICIO: 03 de Enero del 2018

Duración del proyecto en meses referidos al año
de ejecución vigente.
TERMINO: 25 de Junio del 2018

8. CRONOGRAMA DE EJECUCIÓN DEL PROYECTO

Dedicación semanal
Etapas Inicio Término
(Horas)
a. Recolección de datos 03/01/2018 15/03/2018 16
b. Experimentación y análisis de resultados 16/03/2018 15/05/2018 16
c. Redacción del informe 06/05/2018 25/06/2018 16
 9. RECURSOS DISPONIBLES

a. Personal:
- Autores del Proyecto de Investigación.
- Asesor del Proyecto de Investigación.

b. Materiales y Equipos:

Precio Precio
Descripción Cantidad
Unitario, S/. Total, S/.
Materiales y Equipos: 7400,00
Software Matlab Home R2015a version 8.5.0 01 200,00 400,00
Portátil ASUS X555LB i7-5500 8GB 1TB 15.6'' 02 3500,00 7000,00
Servicios: 1 560,00
Servicio de Internet UNT 01 360,00 360,00
Biblioteca de Ingeniería Química (Libros y Revistas) 06 200,00 1 200,00
*Monto Total (S/.): 8960,00
* El Monto de los recursos disponibles está considerado para dos personas durante la elaboración
del proyecto.
c. Locales:
- Laboratorio de Computación e Informática (LABCI)
- Facultad de Ingeniería Química
- Universidad Nacional de Trujillo

10. PRESUPUESTO

*Monto
Partida Descripción
(S/.)
2.3.2 1.1 2 Viáticos y Asignaciones 3700,00
- Movilidad 1000,00
- Alimento 1800,00
- Hospedaje 900,00
2.3.1 99.1 3 Bienes de Consumo 1135,00
- Material de escritorio: 30,00
Papel Bond A-4 Atlas, el Millar 20,00
Lápices, Lapiceros, Resaltador, otros 10,00
- Material de Impresión: 1 000,00
Libros de Consulta 700,00
Revistas o Papers 300,00
- Fotocopias y Transparencias 105,00
2.6.7 1.6 2 Gasto por la Compra de Bienes 450,00
- Impresora Epson L355 Con Sistema Continuo 450,00
2.3.2 2.1 Servicios de Energía eléctrica 360,00
- Consumo de Energía 360,00
2.3.1 99.1 99 Otros bienes 150,00
Monto Total (S/.): 5795,00
* El Monto de los recursos no disponibles está considerado para dos personas durante la
elaboración del proyecto.
 11. FINANCIAMIENTO

a. Con Recursos Universitarios: Los recursos disponibles serán financiados por la

Universidad Nacional de Trujillo, Facultad de Ingeniería Química, LABCI.
b. Con Recursos Externos: No se considera el financiamiento externo.
c. Autofinanciamiento: Los recursos no disponibles serán financiados por los
elaboradores del presente proyecto con un monto de S/. 5 795,00.

II. PLAN DE INVESTIGACIÓN

1. ANTECEDENTES

La recuperación de hidrocarburos ligeros a partir de mezclas gaseosas multicomponente

por absorción con un aceite absorbente es una operación importante en refinerías de
petróleo y plantas de gasolina natural. La fuerza motriz para la transferencia de un
componente ligero desde la fase gaseosa a la fase líquida es proporcional a la diferencia
entre la presión parcial del componente en la fase gaseosa y la presión de vapor de
equilibrio del componente en solución con el aceite pobre. La absorción según la
volatilidad y las proporciones de líquido y vapor hace que los componentes indeseables
estén presentes con los componentes deseados. Para eliminar los componentes indeseables
del aceite enriquecido; las secciones de Stripping y reboilers son frecuentemente
agregados a una unidad básica de absorbedor.
Para sistemas comerciales, la concentración de componentes recuperables en el vapor de
alimentación es apreciable. La absorción de estos componentes provoca un cambio
significativo en las tasas de líquido y vapor en la columna. El calor de absorción es
desarrollado y se recupera como calor sensible por el aceite enriquecido. Por lo tanto,
ocurren variaciones de temperatura significativas en la columna. Los enfriadores
intermedios son frecuentemente agregados para mantener bajas temperaturas en la
columna, eliminando el calor de absorción. Las temperaturas más bajas otorgan una mayor
selectividad a la absorción. Los aceites absorbedores ligeros también se utilizan debido a
su mayor capacidad de absorción por unidad de volumen.

La operación en estado estacionario de una columna reboiled absorption (absorción con

ebullición) de multicomponentes industrial con multi-alimentación, previamente
investigada por Petryschuk (1964) ha sido simulada por dos modelos matemáticos
alternativos: un modelo de etapa de equilibrio y un modelo de eficiencia de bandeja basado
en el concepto de la eficiencia de Murphree. [1]

Ravicz ha desarrollado relaciones de transferencia de masa que combina las eficiencias de

Murphree para etapas no ideales, para predecir los productos de un simple absorbente.
Este método ha sido programado para una computadora digital a gran escala (IBM 704).
[2].

Kremser desarrolló relaciones para predecir la eficiencia de la absorción, el cual requiere

la suposición del concepto de etapa teórica. Estas ecuaciones están en términos del factor
de Absorción o Stripping en cada etapa teórica. [ 3]
Edmister ha desarrollado un factor efectivo que es riguroso para una columna de dos
placas, pero una aproximación para columnas con más placas. Él supone que el factor de
absorción o de separación efectiva es una función de los factores de absorción del terminal.
[4]

El desarrollo de estos modelos requirió el uso de un conjunto de procedimientos de

estimación de propiedades físicas. Cada modelo reproduce de manera cuantitativa el
comportamiento de estado estable de la columna dentro de la exactitud de los datos
disponibles de la planta y, cuando está sujeto a variaciones de parámetros, el modelo de la
etapa de equilibrio verifica el comportamiento cualitativo asociado con la columna. Los
modelos desarrollados por Petryschuk han demostrado ser preferibles y mejora sobre los
modelamientos de columna. [1]

Hasta aquí se desarrolló los métodos de cálculo para sistemas de Absorber-Stripping, en

los cuales el balance de energía no ha sido considerado. Esto nos lleva a tener un límite en
las variables a especificar para el modelamiento de este sistema.

Se establece un modelo teórico para una etapa de equilibrio en el cual se hace un balance
de Materia, Equilibrio, Energía y se considera la Sumatoria de fracciones molares, a lo
cual Wang y Henke los denomina como ecuaciones MESH. [5]
Existe una gran variedad de métodos iterativos para la resolución de ecuaciones
algebráicas no lineales (como son las ecuaciones MESH) que consiste en la solución de
las ecuaciones por tanteo y/o linealizacion por las técnicas de Newton Raphson. Se ha
propuesto varios intentos en la resolución de las ecuaciones MESH como: Cálculo de etapa
a etapa y ecuación a ecuación de Lewis – Matheson (1932) y Thiele – Geddes (1933),
Método Theta de Holland; llegando a concluir que su solución en forma manual es tediosa.
Resulta sencillo utilizar un ordenador digital. [6, 7, 8]

En un estudio clásico, Friday y Smith analizaron de forma sistemática una serie de

métodos de tanteo para la resolución de las ecuaciones MESH, considerando
cuidadosamente la elección de la variable de salida para cada ecuación. Demostraron que
no hay una técnica única capaz de resolver todos los tipos de problemas. Para el caso de
separadores donde la alimentación está constituida por componentes de volatilidades
parecidas (puntos de ebullición próximos), recomendaron un método de Amundson-
Pontinen modificado, denominado método del punto de burbuja (BP). En el caso de
alimentaciones con componentes de muy diferente solubilidad o volatilidad (puntos de
ebullición alejados), el método BP mostró encontrarse sometido a fallos, sugiriéndose el
método denominado de suma de flujos (SR). Para casos intermedios el método de tanteo
puede no converger; en este caso Friday y Smith sugieren la utilización del método de
Newton-Raphson, o bien una combinación de éste con un método de tanteo. El método de
Newton-Raphson permite una considerable flexibilidad en la elección de las variables
especificadas y, por lo general, es capaz de resolver todos los problemas. [9]

Este método está basado en la resolución de todas las ecuaciones MESH, o una
combinación de las mismas, por técnicas de correcciones simultáneas (SC).
El desarrollo de los procedimientos SC implica varias decisiones importantes:
1) ¿Qué variables se deben emplear?
2) ¿Qué ecuaciones se deben utilizar?
3) ¿Cómo se deben ordenar las variables?
4) ¿Cómo se deben ordenar las ecuaciones?
5) ¿Cómo se debe proporcionar la flexibilidad en las especificaciones?
6) ¿Cuáles derivadas de las propiedades físicas se deben mantener?
7) ¿Cómo deben linealizarse las ecuaciones?
8) Si se utilizan técnicas de linealización de Newton o cuasi-Newton, ¿Cómo debe
actualizarse el Jacobiano?
9) Las correcciones de las incógnitas que se calculan en cada iteración ¿deben
modificarse para amortiguarlas o acelerar la solución o mantenerlas dentro de ciertos
límites?
10) ¿Qué criterio de convergencia debe especificarse? [10]

Naphtali and Sandholm desarrollan un nuevo enfoque para los cálculos de separación que
tiene muchas ventajas prácticas. Las ecuaciones de conservación de masa y energía y de
equilibrio se agrupan por etapas y luego se linealizan. El conjunto resultante de ecuaciones
tiene una estructura bloque-tridiagonal que permite la solución mediante una técnica
simple. Así, se diseña una nueva forma de aplicar la técnica de Newton-Raphson a los
cálculos de separación. El método ha sido probado en una serie de problemas elegidos
para exhibir características que hacen que falle otra técnica de solución. La mayoría de los
problemas se resuelven con diez iteraciones o menos. [11] Fredenslund, Gmehling y
Rasmussen dan un programa de cálculo para este método. [12]

La aplicación de esta metodología nos permite que las especificaciones para cualquier tipo
de problema sean bastante flexibles. Presión, composiciones, flujos y etapa de localización
son necesarias para todas las alimentaciones.
La condición térmica de cada alimentación puede darse en función de la entalpía,
temperatura o fracción molar vaporizada. Si una alimentación consta de dos fases, éstas
pueden enviarse a la misma etapa, o bien el vapor puede dirigirse a la etapa superior
respecto a la diseñada para alimentación. La presión y la eficacia de las etapas puede
diseñarse especificando los valores de cabeza y cola. El resto de los valores se obtiene por
interpolación lineal. Por defecto, se supone un comportamiento adiabático para las etapas
intermedias, a no ser que se especifiquen los valores de calor y temperatura en cada etapa.

En 1974, Boston y Sullivan presentaron un algoritmo diseñado para reducir el tiempo

dedicado a calcular las propiedades termodinámicas al hacer cálculos de columnas. Se
emplean dos conjuntos de modelos de propiedades termodinámicas: (1) un conjunto
empírico simple aproximado que se usa frecuentemente para converger cálculos de bucle
interno, y (2) un conjunto riguroso utilizado con menos frecuencia en el bucle externo .
Las ecuaciones MESH siempre se resuelven en el bucle interno con el conjunto
aproximado. Los parámetros en las ecuaciones empíricas para el conjunto de bucle interno
se actualizan con poca frecuencia en el bucle externo utilizando las ecuaciones rigurosas.
La característica distintiva de Boston-Sullivan son los bucles internos y externos, de ahí
el nombre de método de INSIDE – OUT. [13]

En 1979, Lars J. Christiansen, Michael L. Michelsen and Aage Fredenslund aplicaron el

método el desarrollo del proceso y diseño de columnas de destilación asociadas con la
 separación de líquidos de gas natural (LGN). Esto logró una convergencia rápida y segura
incluso cerca de la región crítica.

En el 2001, se publicó en Computers & Chemical Engineering la introducción del

concepto de mezcla pseudo binaria el algoritmo de Ishii – Otto para cálculos en multiples
etapas. Este algoritmo resuelve las ecuaciones del modelo de la etapa de forma simultánea
y estructura las ecuaciones linealizadas de manera que las manipulaciones de la matriz
tratan casi por completo con matrices tridiagonales. Esto reduce en gran medida el tiempo
de computación y los requisitos de almacenamiento de la computadora. Las características
de convergencia para sistemas altamente no ideales se mejoran notablemente al emplear
el concepto de PBM para tener en cuenta la dependencia de la composición de los valores
K de equilibrio. Se conservan las ventajas y la simplicidad de la formulación original del
algoritmo de Ishii y Otto. La solidez del método se ilustra con una serie de difíciles
problemas de separación que se sabe que son altamente sensibles a las suposiciones
iniciales para los perfiles de temperatura. [15]

2. FUNDAMENTACIÓN CIENTIFÍCA, TÉCNICA O HUMANÍSTICA

La unidad de craqueo catalítico de fluidos (FCCU) es una de las instalaciones de refinación

de petróleo más importantes en el proceso. La unidad de recuperación de vapor (VRU) es
la parte de separación de FCCU, que procesa el gas rico y el destilado de la parte superior
del fraccionador para producir gasolina y gas licuado de petróleo con el gas seco como
producto secundario. Una planta de concentración de gas de las FCCU o una unidad de
recuperación de vapor generalmente consta de: Reboiled Absorber - Stripper,
Debutanizador, Despropanizador, Gasoline splitter. Solo el reboiled absorber presenta
nuevos problemas únicos para las operaciones de FCCU. [16]

El Reboiled Absorber se refiere a la operación de un absorbedor con ebullición para

separar al menos un hidrocarburo de bajo punto de ebullición de otro hidrocarburo que
también tiene bajo punto de ebullición. En otro aspecto, se refiere a la separación de
metano de una mezcla que incluye hidrocarburos C2. En muchas separaciones, se
encuentra que el uso de absorbedores con ebullición es el tipo de operación más eficiente
para la separación de uno o más hidrocarburos de bajo punto de ebullición. Dichos
absorbedores consisten generalmente en una sección absorbente y una sección de
stripping, es decir, la entrada de alimentación está aproximadamente en el medio de la
torre con la sección encima del nivel de entrada de alimentación que es la sección
absorbente y la porción inferior es la sección de stripping. Se agrega un absorbente líquido
a la parte superior de la sección de absorción mientras que la parte inferior de la sección
de stripping se opera como un hervidor. Este método de operación da como resultado una
cantidad casi constante de líquido que fluye hacia abajo en ambas secciones de la torre.
Para la operación más eficiente, encontramos que la cantidad de líquido que fluye hacia
abajo en la sección de absorción debe ser tan grande como sea posible y la cantidad de
líquido que fluye hacia abajo en la sección de stripping debe ser tan pequeña como sea
posible. [17]

Para fijar los parámetros de operación de este reboiled absorber, su modelamiento y

simulación involucra el desarrollo de las ecuaciones MESH en cada etapa de la columna.
Desafortunadamente, estas ecuaciones son relaciones algebraicas no lineales que
 interaccionan entre sí fuertemente, lo que hace que los métodos de resolución se vuelvan
complejos y tediosos.

Para la solución de estas ecuaciones, se aplica el método de Newton-Raphson, que está

basado en las resoluciones de todas las ecuaciones MESH por técnicas de corrección
simultánea (SC). Este método nos permite una considerable flexibilidad en la elección de
variables especificadas y, por lo general, es capaz de resolver todos los problemas.

El método de SC de Naphtali y Sandholm se desarrolla con detalle, dado que utiliza

muchas de las técnicas matemáticas de métodos de tanteo. Se llega a desarrollar las
siguientes ecuaciones MEH:

Balance de materia:

𝑀𝑖𝑗 = 𝑙𝑖𝑗 (1 + 𝑠𝑗 ) + 𝑣𝑖𝑗 (1 + 𝑆𝑗 ) − 𝑙𝑖𝑗−1 − 𝑣𝑖𝑗+1 − 𝑓𝑖𝑗 = 0

Equilibrio de Fases:

∑𝑐𝑘=1 𝑣𝑘𝑗
𝐸𝑖𝑗 = 𝑘𝑖𝑗 𝑙𝑖𝑗 − 𝑣𝑖𝑗 = 0
∑𝑐𝑘=1 𝑙𝑘𝑗

Balance de Energía:
𝑐 𝑐 𝑐 𝑐 𝑐

𝐻𝑗 = ℎ𝐿𝑗 (1 + 𝑠𝑗 ) ∑ 𝑙𝑖𝑗 + ℎ𝑉𝑗 (1 + 𝑆𝑗 ) ∑ 𝑣𝑖𝑗 − ℎ𝐿𝑗−1 ∑ 𝑙𝑖𝑗−1 − ℎ𝑉𝑗+1 ∑ 𝑣𝑖𝑗+1 − ℎ𝐹𝑗 ∑ 𝑓𝑖𝑗 − 𝑄𝑗 = 0
𝑖=1 𝑖=1 𝑖=1 𝑖=1 𝑖=1

Si N y todas las 𝑓𝑖𝑗 , 𝑇𝐹𝑗 , 𝑃𝐹𝑗 , 𝑃𝑗 , 𝑠𝑗 , 𝑆𝑗 y 𝑄𝑗 , las ecuaciones M, E y H son no lineales con
variables de salida 𝑣𝑖𝑗 , 𝑙𝑖𝑗 , y 𝑇𝑗 , donde i=1 hasta c y j=1 hasta N. Las ecuaciones se
resuelven de forma simultánea por el método iterativo de Newton–Raphson, donde se
generan sucesivos conjuntos de variables de salida hasta que las funciones MEH alcancen
una cierta tolerancia en torno de cero. En las iteraciones, a los valores de las funciones
distintos de cero se les denomina discrepancia o error. Se agrupan las funciones y variables
de salida por etapas de cabeza a cola. Como se verá posteriormente, esto se hace así para
obtener una estructura de bloque tridiagonal en la matriz jacobiana de las derivadas
parciales con el fin de poder aplicar el algoritmo de Thomas. Sea:

donde 𝑋𝑗 es el vector de variables de salida para la etapa j dispuestos en el siguiente orden:

y 𝐹𝑗 es el vector de funciones para la etapa j dispuesto en el orden:

La iteración de Newton-Raphson se desarrolla resolviendo (que en forma matricial

resulta):
 donde ∆𝑋 son las correcciones de las variables de salida. Estas correcciones se utilizan
para calcular la siguiente aproximación del conjunto de variables de salida mediante:

La cantidad de (𝜕𝐹̅̅̅̅ ⁄𝜕𝑋) es el siguiente jacobiano o matriz (N X N) de bloques de

derivadas parciales de todas las funciones con respecto a todas las variables de salida:

Cada bloque 𝐴𝑗 , 𝐵𝑗 𝑜 𝐶𝑗 representa una submatriz de derivadas parciales de las funciones

de la etapa j.

La suma de los cuadrados de discrepancia, basada en las suposiciones iniciales para todas
las variables de salida, se calcula y compara mediante un criterio de convergencia:

Para conseguir que los valores de todas las discrepancias sean el mismo orden de
magnitud, es necesario dividir las funciones del balance de energía 𝐻𝑗 por un factor de
escala del orden del calor latente de vaporización.

El criterio de convergencia se calcula a partir de:

3. JUSTIFICACIÓN

El modelamiento de una columna de destilación de multicomponentes, cabe mencionar

los casos de separaciones en las que intervienen mezclas liquidas con una fuerte no
idealidad (por ejemplo, en destilación extractiva) o en aquellos casos donde el separador
presenta un comportamiento de absorbedor o agotamiento, y de fraccionador (por ejemplo,
un absorbedor con ebullición ); es bastante complejo debido a la gran cantidad de
parámetros y variables involucradas y también, por lo general, se requiere resolver el
conjunto de ecuaciones no lineales (ecuaciones MESH) y ecuaciones diferenciales que
intervienen en el modelamiento.

Por ello, el disponer de técnicas prácticas para modelar operaciones de separación de

contacto vapor – liquido más o menos realistas y el desarrollar sistemas de control eficaces
y fiables, es muy importante a fin de conseguir un funcionamiento eficaz y seguro de los
sistemas de separaciones multicomponentes Industriales.
 Las disponibilidades de grandes computadoras digitales hacen posibles las soluciones
rigurosas de los modelos de etapas de equilibrio, para separación de multicomponentes y
etapas múltiples, hasta un grado de exactitud limitado exclusivamente por la precisión de
los datos de entalpía y equilibrio de fase. El tiempo y costo para obtener dichas soluciones
es muy bajo, en comparación con el costo de las soluciones manuales.

4. PROBLEMA

¿Cómo realizar el modelamiento y simulación de un Reboiled Absorber en estado

estacionario con todos los parámetros que demande aplicando el método de correcciones
simultáneas de Naphtali y Sandholm?

5. HIPÓTESIS

El modelamiento y simulación de un Reboiled Absorber en estado estacionario se

realizará aplicando las ecuaciones de MESH (balances de materia, relaciones de
equilibrio, sumatoria de fracciones molares y balances de entalpía) y el método de
Correcciones Simultáneas de Naphtali y Sandholm mediante programación con Matlab
R2015a 8.5.0, lo que permitirá evaluar el comportamiento del proceso propuesto.

6. OBJETIVOS

a. Objetivo General:

Modelar y simular un Reboiled Absorber de multicomponente en estado estacionario

aplicando el método de Correcciones Simultáneas de Naphtali y Sandholm.

b. Objetivo Específico:

1. Formular las ecuaciones desarrolladas MEH (balances de materia, relaciones de

equilibrio y balances de entalpía) y ordenar dichas funciones y variables de salida
para obtener una estructura de bloque tridiagonal en la Matriz Jacobiana de las
derivadas parciales.
2. Elaborar el programa computacional usando el lenguaje de programación Matlab
R2015a 8.5.0.
3. Aplicar el programa computacional a los casos de estudio propuesto en el texto
Fouad M. Khoury (2004) a los Ejemplos: 9-6 Pag 226. [19]
4. Contrastar los resultados obtenidos con el programa computacional
‘‘ReboiledAbsorberSC-NS’’ con los reportados mediante Aspen Hysys 8.1,
ChemCad 6.1.3. y Fouad M. Khoury del texto antes mencionado.

7. METODOLOGÍA DE TRABAJO

Para el desarrollo del programa computacional se aplicará:

a. Las ecuaciones MEH desarrolladas por Naphtali y Sandholm:

- Balances de materiales de los componentes para cada etapa

 - Relaciones de Equilibrio de los componentes para cada etapa.
- Balances de Energía en cada plato

b. Procedimiento de cálculo y resolución del modelo matemático para separaciones

multietapas de multicomponentes de mezclas no ideales:

1) Se especifican los flujos y condiciones de alimentación, todos los flujos de calor

que se llegasen a dar en cada etapa excepto los de cabeza y fondo, las corrientes
laterales, presión de la columna, eficiencia de cada plato y numero de etapas.
2) La iteración de Newton de inicia al proporcionar suposiciones razonables para
todas las incógnitas las cuales se pueden generar a partir de las suposiciones de
Temperaturas (T1, TN) flujos de líquido (L1, LN) y vapor (V1, VN) de la cabeza y
de fondo. Los valores restantes de Tj, Lj y Vj se obtienen por interpolación lineal.
3) Una estimación para las fracciones molares de cada etapa se obtiene por
cualquiera de estas dos técnicas: la primera se obtiene realizando cálculos de
flash para una combinación de alimentaciones a una presión media y una
relación V/L que se aproxime a la relación entre los flujos de destilado y de colas
obteniendo así dichas composiciones en fracciones molares para las fases de
vapor y liquido en equilibrio, la otra técnica consiste en la primera iteración del
método BP. Así por medio de cualquiera de estas dos técnicas obtendremos las
suposiciones iniciales lj y vj.
4) Se calcula la suma de los cuadrados de las funciones de discrepancia, basadas en
las suposiciones iniciales para todas las variables de salida y se compara
mediante un criterio de convergencia.
5) Se procede a la segunda iteración en las cuales los flujos de vij, lij y los valores
de temperatura Tj son cambiados las correcciones de Newton – Raphson.
Alternativamente puede ultilizarse un factor de avance escalar no negativo t en
las correcciones.
6) Al haber convergido el problema, se calculan los valores de Vj, Lj, Q1 y QN.

Este procedimiento de cálculo, se supone que la caída de presión de un plato a

otro es despreciable en relación con la presión total.
 MODELAMIENTO Y SIMULACIÓN DE UN
REBOILED ABSORBER DE
MULTICOMPONENTES EN ESTADO
ESTACIONARIO POR EL MÉTODO DE
CORRECCIONES SIMULTÁNEAS EN LA UNIDAD
DE FRACCIONAMIENTO DE CRAQUEO
CATALÍTICO

 Balance de materia
Desarrollar las ecuaciones
 Relaciones de equilibrio
MEH
 Balances de energía

Especificar Fj, zij, condiciones

de alimentación (TFj, PFj, HFj),
Pj, ηj, N, todas las Qj y QN una
variable para la etapa de
cabeza y de fondo.
Ajustar k=1 (para iniciar la primera
iteración)
Iniciar variables de tanteo Tj,
Vj, Lj

Calcular la suma de los

cuadrados de la función de
discrepancia

¿Es 𝜏 < 𝜀? Calcular Vj y Lj

Ajustar k=k+1 No

Calcular correcciones de Calcular Q1 de H1 y QN de HN

Newton Raphson. sino se especifican.

Salida
Calcular t óptimo para
minimizar 𝜏. Calcular nuevos Resolución simultánea de las ecuaciones
valores de vij, lij, Tj
no lineales por el método de Newton -
Raphson mediante programación con
Matlab.

Figura1. Esquema general del modelamiento y simulación de una Separación de

multicomponentes.
8. POSIBLE CONTRIBUCIÓN E IMPACTO

 El proyecto de esta investigación contribuirá en la industria de hidrocarburos, destilería

y otros afines, permitiendo reducir el tiempo empleado en el análisis de los procesos
de destilación, se logrará simplificar los procesos productivos, haciendo posible una
experimentación controlada y así evitar el gasto innecesario en la compra de un
software con un costo muy elevado. Además, demostrar que modelar y simular un
programa computacional es una alternativa viable en lo que respecta al control de
gastos de Software.

 Crear un banco de programas de simulación para nuestra Facultad de Ingeniería

Química.

9. REFERENCIAS BIBLIOGRAFICAS

[1] D. R. Prowse & A. I. Johnson (1972). MultiComponent Reboiled Absorber Modeling

and Parametric Studies. The Canaddim Journal of Chemical Engineering, Vol. 50,
April, 1972.
[2] A. E. Ravicz (1959). Non-Ideal Stage Multicomponent Absorber Calculations by
Automatic Digital Computern University of Michigan, 1959.
[3] A. Krernser (1930). Nat. Pet. News, 22, No. 21, p. 48. 1930.
[4] W. C. Edmister (1947). Pet. Eng., 18, No. 13, pp. 130-144. 1947.
[5] J. C. Wang & G. E. Henke (1966). Hydrocarbon Processing, 45 (8), 155-163. 1966.
[6] W. K. Lewis & G. L. Matheson (1932) Ind. Eng. Chem., 24, 496-498. 1932.
[7] E. W. Thiele & R. L. Geddes (1933). Ind. Eng. Chem., 25, 290. 1933.
[8] C. D. Holland (1963). Multicomponent Distillation. Prentice-Hall, Inc., Englewood
Cliffs, N. J. 1963.
[9] J. R. Friday & B. D. Smith (1964). AIChe J., 10, 698-707. 1964..
[10] Robert H. Perry & Don W. Green (1994). Manual del Ingeniero Químico. Sexta
Edición MC Graw Hill.1994
[11] Leonard M. Naphtali and Donald P. Sandholm (1971). Muhicomponent Separation
CaIcuIations by Linearization. AlChE Journal, January, 1971.
[12] Fredenslund A., Gmehling J. & Rasmussen P. (1977). Vapor-Liquid Equilibria using
UNIFAC, A Group Contribution Method. Elsevier, Amsterdam. 1977
[13] Boston J. F. & S. L. Sullivan (1974). Can. J. Chem. Engr., 52, 52-63. 1974.
[14] Lars J. Christiansen, Michael L. Michelsen and Aage Fredenslund (1979). Naphtali-
Sandholm Distillation Calculations for NGL Mixtures near the critical region.
Computers & Chemical Engineering, Vol. 3, pp. 535-542, 1979.
[15] Yoshikazu Ishi & Fred D. Otto (2001). An efficient simultaneous correction
procedure for multicomponent, multistage separation calculations for non-ideal
systems. Computers and Chemical Engineering 25, 1285-1298. 2001.
[16] Norman P. Lieberman (2009). Troubleshooting process operations. Fourth Edition.
PennWell Corporation. 2009.
[17] Robert L. McIntire and Russell O. Shelton (1955). Reboiled Absorber Operation.
Phillips Petroleum Company. 1955.
[18] J. D. Seader, Ernest J. Henley, D. Keith Roper 1998. Separation Process Principles:
chemical and biochemical operations. Third Edtion. John Wiley & Sons, Inc.
[19] Fouad M. Khoury 2004. Multistage Separation Processes. Third Edtion. CRC Press,
2000 N.W. Corporate Blvd., Boca Raton, Florida 33431.
[20] Giorgio Soave, 1972. Equilibrium constants from a modified Redkh-Kwong
equation of state. Chemical Engineering Science, 1972.
[21] Otto Redlich and J. N. S. Kwong (1948). On the Thermodynamics of Solutions. An
equation of state. Fugacities of gaseous solutions Shell Development Company,
Emeryville, California, 1948.
[22] Chao K. C. and J. D. Seader (1961). A General Correlation of Vapor-Liquid
Equilibria in Hydrocarbon Mixtures. A.I.Ch.E. Journal, 1961.
[23] Shah M. K. and Bishnoi P. R. (1978). Multistage Multicomponent Separation
Calculations using Thermodynamic Properties Evaluated by the SRK/PR Equation
of State. The Canadian Journal of Chemical Engineering, Vol. 56, 1978.
[24] Guevara Ureña Alejandro (2012). Prediccion y validacion del Equilibrio Liquido –
Vapor para una mezcla multicomponente No Ideal mediante el modelo UNIFAC.
[Tesis para optar el título]. Trujillo: Universidad Nacional de Trujillo. Facultad de
Ingeniería Química.
[25] Yeong Koo Yeo (2017) Chemical Engineering Computation with MATLAB.. Taylor
& Francis Group. 2017
[26] W. Fred Ramirez (1998). Computational Methods in Process Simulation. Elsevier
Science & Technology Books. 1998.
[27] Gilat, A. and Subramaniam, V. 2014. Numerical methods for engineers and
scientists. Third Edition. John Wiley & Sons, Inc. U.S.A.
[28] Chapra, S. 2011. Métodos numéricos para ingenieros. Sexta Edición. McGraw-
Hill/Interamericana Editores, S.A. de C.V. México.

Robert Eduardo Graus Peláez Jorge Pablo Escobar Esquivel

Tesista 1 Tesista 2

Ms. Walter Moreno Eustaquio

Asesor