Articulos Mezclas

Red de Revistas Cientficas de Amrica Latina, el Caribe, Espaa y Portugal
Sistema de Informacin Cientfica
Francisco Snchez Mares, Pedro Alonso Dvila

Prediccin de Puntos Crticos en Mezclas de Gas Natural
Conciencia Tecnolgica, nm. 34, julio-diciembre, 2007, pp. 32-36,
Instituto Tecnolgico de Aguascalientes
Mxico
Disponible en: http://www.redalyc.org/articulo.oa?id=94403407
Conciencia Tecnolgica,
ISSN (Versin impresa): 1405-5597
contec@correo.ita.mx
Mxico
Cmo citar?
Fascculo completo
Ms informacin del artculo
Pgina de la revista
www.redalyc.org
Proyecto acadmico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
Conciencia Tecnolgica No. 34, Julio - Diciembre 2007
Prediccin de Puntos Crticos en Mezclas de Gas Natural

Investigacin
M. C. Francisco Snchez Mares y Dr. Pedro Alonso Dvila

Doctorado en Ciencias en Ingeniera Qumica
Universidad Autnoma de San Lus Potos
Av. Dr. Manuel Nava No. 6, Zona Universitaria, San Lus Potos, S. L. P.
fsanchez96@hotmail.com
Resumen
Actualmente es el combustible ms limpio para la

generacin de energa elctrica, ofrece las mejores
opciones en trminos econmicos, de rendimiento e
impacto ambiental. Es utilizado en diversos procesos
industriales como la produccin de cermica, cemento,
vidrio, textiles, en los sectores comercial y domstico
y como combustible en vehculos.
Debido a que el gas natural se ha convertido en un
combustible importante, para procesarlo es necesario
conocer su comportamiento de fases, el cual es
descrito por su envolvente y el punto crtico.
La ubicacin del punto crtico sobre la envolvente
de fases proporciona informacin importante; por
ejemplo, si la mezcla est principalmente formada de
componentes ligeros o pesados o si existe condensacin
retrgrada. Conociendo lo anterior se pueden
determinar las alternativas adecuadas de diseo que
debern emplearse en la extraccin y procesamiento
del gas del yacimiento.
Como se puede observar el clculo de puntos crticos
en mezclas de gas natural es una actividad importante.
No obstante, el software comercial disponible para
realizar estos clculos tiene un costo elevado. Por
lo tanto, el objetivo de este trabajo es implementar
un mtodo eficiente para calcular puntos crticos en
mezclas de gas natural resolviendo la formulacin de
Heidemann y Khalil (1980) con un mtodo de Newton
multivariable.
La prediccin confiable de puntos crticos en mezclas

de gas natural permite conocer su comportamiento de
fases. En este trabajo se presenta el clculo de puntos
crticos de dos mezclas de gas natural resolviendo
la formulacin de Heidemann y Khalil (1980) con
un mtodo de Newton multivariable. Los resultados
muestran que el mtodo es eficiente y puede ser
aplicado a mezclas multicomponentes.
Palabras Clave
Puntos crticos, Gas
Heidemann y Khalil.
natural,
Formulacin
de
Abstract
The reliable calculation of critical points in
natural gas mixtures allows to predict the behavior
of phases. The calculation of critical points of two
mixtures of natural gas solving the Heidemann and
Khalils formulation (1980) with a multivariable
Newtons method is presented in this work. The results
show that method is efficient and can be applied to
multicomponent mixtures.
Keywords
Critical point, Natural gas, Heidemann and Khalils
formulation.
Trabajos previos
Peng y Robinson (1977) calcularon las propiedades
crticas de mezclas multicomponentes empleando su
famosa ecuacin de estado. El mtodo propuesto se
basa en la solucin de un conjunto de determinantes
derivados de la transformacin de Legendre de la
energa libre de Gibbs.
Hicks y Young (1977) calcularon los puntos crticos
de diversas mezclas binarias utilizando un mtodo de
bsqueda exhaustiva del punto donde coinciden las
derivadas de las condiciones crticas en trminos de la
energa libre de Helmholtz.
El mtodo ms empleado actualmente para calcular
puntos crticos en mezclas multicomponentes es el
propuesto por Heidemann y Khalil (1980) que emplea
una formulacin propia y que es aplicable a ecuaciones
Introduccin
El gas natural es una mezcla de hidrocarburos
parafnicos ligeros que tiene como principal componente
al metano CH4. Tambin contiene pequeas cantidades
de etano C2H6, propano C3H8 y en proporciones
variables nitrgeno N2, dixido de carbono CO2 y
cido sulfhdrico H2S (PEMEX, 2005).
Se puede encontrar en dos formas distintas asociado
y no asociado. Se denomina asociado cuando se
encuentra acompaado de petrleo, mientras que
no asociado es cuando se encuentra en yacimientos
exclusivos de gas natural.
32
PREDICCIN DE PUNTOS CRTICOS EN MEZCLAS DE GAS NATURAL

M.C. Francisco Snchez Mares y Dr. Pedro Alonso Dvila
de estado cbicas. La convergencia al punto crtico se

logra mediante un mtodo de Newton con dos lazos
de iteracin.
Michelsen (1984) propone una metodologa para
calcular puntos crticos de mezclas multicomponentes
empleando una modificacin del plano tangente de
Gibbs, mientras que Stradi y col. (2001) resolvieron la
formulacin de Heidemann y Khalil (1980) empleando
un mtodo de Newton con intervalos definidos
para calcular puntos crticos de mezclas binarias y
ternarias.
Henderson y col. (2004), Freitas y col. (2004) y JustoGarca y Garca-Snchez (2005) abordaron el clculo
de puntos crticos en mezclas multicomponentes como
un problema de optimizacin global donde la funcin a
minimizar era el criterio del plano tangente de Gibbs.
Recientemente, Snchez-Mares y Bonilla-Petriciolet
(2006) resolvieron la formulacin de Heidemann
y Khalil (1980) con Simulated Annealing (Goffe y
col., 1994) y determinaron las propiedades crticas de
diversas mezclas multicomponentes.
Solucin de las condiciones crticas de Heidemann

y Khalil (1980)
Las Ecuaciones (1) (4) fueron resueltas empleando
un mtodo de Newton multivariable, que para
converger a la solucin depende de una estimacin
inicial de las incgnitas adecuadas. Los valores
iniciales de Tc y Vc se calcularon de la siguiente
forma
donde zi es la composicin en la alimentacin del

componente i.
Para el caso de
la inicializacin fue
Fundamentos tericos
Heidemann y Khalil (1980) propusieron que para
determinar el punto crtico de una mezcla se deben
cumplir las siguientes condiciones
Materiales y mtodos
a) Las ecuaciones de estado cbicas de PengRobinson (PR) y Soave-Redlich-Kwong (SRK).
Las ecuaciones de estado cbicas PR y SRK en su
forma estndar son las siguientes
Ecuacin SRK
donde C es el nmero de componentes, A es la

energa libre de Helmholtz, ni representa el nmero de
moles del componente i,
representa
una perturbacin en el nmero de moles diferente
de cero y Q es una matriz cuadrada formada por los
elementos
donde P es la presin del sistema, T es la temperatura,

V es el volumen y R es la constante universal de los
gases. Los parmetros restantes se definen de la
siguiente forma
Las expresiones (1) (4) conforman un sistema de

C + 2 ecuaciones no lineales donde la temperatura
crtica Tc, el volumen crtico Vc y
sern las
incgnitas a determinar. Para calcular la presin crtica
Pc se sustituyen Tc y Vc calculados en la ecuacin de
estado cbica bajo anlisis.
33

donde Tc es la temperatura crtica de componente puro,

Pc es la presin crtica de componente puro, es el
factor acntrico del componente puro y Tr = T / Tc .
d) Propiedades crticas de componente puro.

Las propiedades crticas de componente puro
empleadas en el clculo de los puntos crticos de las
mezclas de gas natural se presentan en la Tabla 2 y
fueron obtenidas de Michelsen y Mollerup (2004).
Componente
CH4
C2H6
C3H8
i-C4H10
n-C4H10
i-C5H12
n-C5H12
n-C6H14
n-C7H16
N2
CO2
Tc (K)
190.6
305.4
369.8
408.1
425.2
460.4
469.6
507.4
540.2
126.2
304.2
Pc (MPa)
4.599
4.883
4.244
3.647
3.799
3.383
3.373
2.968
2.735
3.394
7.375
0.008
0.098
0.152
0.176
0.193
0.227
0.251
0.296
0.351
0.040
0.225
Tabla 2. Propiedades de componente puro

empleadas en el clculo de puntos crticos de mezclas
de gas natural.
e) Programa de cmputo empleado para calcular
puntos crticos en mezclas de gas natural.
b) Reglas de mezclado utilizadas.
Las reglas de mezclado empleadas para las dos
ecuaciones de estado fueron las convencionales, es
decir, las de Van der Waals y se definen de la siguiente
forma
La Figura 1 presenta el diagrama de flujo del mtodo

empleado para calcular puntos crticos en mezclas de
gas natural.
Los parmetros de interaccin binaria utilizados

en este trabajo se tomaron como cero para las mezclas
entre alcano-alcano. Para las mezclas entre no alcano alcano se utilizaron los valores presentados en la Tabla
1. Todos los parmetros de interaccin binaria fueron
tomados de Michelsen y Mollerup (2004).
N2
CO2
N2
CO2
CH4
C2H6
C3H8
n-C4H10
0.02
0.12
0.06
0.15
0.08
0.15
0.08
0.15
n-C5H12
n-C6H14
N2
CO2
0.08
0.15
0.08
0.15
0.00
0.00
0.00
0.00
Tabla 1. Parmetros de interaccin binaria kij

empleados para mezclas binarias no alcano-alcano.
Tc,Vc y
c) Parmetros de interaccin binaria empleados.
Proponer nuevos
valores de
donde kij es el parmetro de interaccin binaria.
1.0E-6
Figura 1. Diagrama de flujo para calcular las

condiciones crticas de mezclas de gas natural
empleando el mtodo de Newton multivariable.

34

El programa de cmputo empleado en este trabajo

para obtener las condiciones crticas de mezclas de
gas natural fue codificado en lenguaje FORTRAN y se
compone de un programa principal y dos subrutinas.
La primer subrutina contiene la codificacin de
la formulacin de Heidemann y Khalil (1980), es
decir las Ecuaciones (1) (4), mientras que en la
segunda subrutina se encuentra el mtodo de Newton
multivariable.
El criterio de paro margen de error establecido
en el programa fue 1.0E-6 y los clculos se realizaron
en un equipo de cmputo porttil con procesador Intel
Celeron M 1.50 GHz y 240 MB de RAM.
donde el superndice cal es el valor de la propiedad

crtica calculada con la ecuacin de estado y el
subndice exp es el valor experimental.
En general se puede observar que las dos ecuaciones
de estado predicen de manera adecuada las propiedades
crticas.
Resultados y discusin
Mezcla de gas natural con ocho componentes
En esta seccin slo se presentan los resultados

obtenidos para dos mezclas diferentes. A continuacin
se presentan los sistemas evaluados.
La composicin de esta mezcla de gas se presenta en la

Tabla 5 y tiene como caracterstica principal la adicin
de dixido de carbono en pequea cantidad. Para esta
mezcla no existen reportados datos experimentales en
la literatura.
Gas de Michelsen (1980)

Esta mezcla de gas natural fue estudiada por
Michelsen (1980) y tiene la composicin mostrada en
la Tabla 3.
Componente
CH4
Composicin (%)
94.30
C2H6
2.70
C3H8
0.74
n-C4H10
0.49
n-C5H12
0.10
n-C6H14
0.27
N2
1.40
La Tabla 4 presenta el punto crtico calculado

para esta mezcla de gas con las ecuaciones de estado
cbicas de PR y SRK y la comparacin con el valor
experimental reportado.
SRK
Tc
(K)
Pc
(MPa)
Composicin (%)
95.08
2.14
0.29
0.11
0.04
0.01
1.94
0.39
Tabla 5. Componentes y composiciones de una

mezcla de gas natural con adicin de CO2.
Tabla 3. Componentes y composiciones del gas de

Michelsen (1980).
Experimental
Tc
Pc
(K)
(MPa)
Componente
CH4
C2H6
C3H8
n-C4H10
n-C5H12
n-C6H14
N2
CO2
PR
Tc
(K)
Pc
(MPa)
207.09 5.578
202.39 5.830
202.47 5.850
DA
2.27
4.52
2.23
4.88
Tabla 4. Punto crtico calculado para el gas de
Michelsen (1980) con las ecuaciones de estado cbicas
PR y SRK y su comparacin con el valor experimental.
La Desviacin Absoluta (DA) para la temperatura y
presin crtica se calcul de la siguiente forma,
El punto crtico calculado para esta mezcla se

reporta en la Tabla 6.
Tc
(K)
197.53
SRK
Pc
(MPa)
Tc
(K)
5.24
197.36
PR
Pc
(MPa)
5.23
Tabla 6. Punto crtico predicho para la mezcla de gas

natural de ocho componentes.
En general, las ecuaciones de estado cbicas SRK y
PR son confiables en la prediccin de puntos crticos.
Por lo tanto, el software comercial empleado en la
industria petroqumica est basado en estas ecuaciones.
Sin embargo, en algunos casos se pueden presentar
errores en la prediccin. Los cuales vienen dados por
la formulacin terica de la ecuacin de estado cbica,
no por el mtodo de solucin de las condiciones
crticas.

35

Debido a lo anterior es de esperar que se presenten

errores en la prediccin de las propiedades crticas de
mezclas que contengan alcoholes o agua y en sistemas
con una elevada cantidad de nitrgeno. El error entre
los sistemas hidrocarburos-nitrgeno es ocasionado
por la alta asimetra molecular entre estos compuestos,
mientras que para las mezclas hidrocarburos-agua o
hidrocarburos-alcoholes es debido a la polaridad de
las molculas.
Para el caso de los sistemas evaluados en este trabajo
el punto crtico calculado es confiable, debido a que
existe poca cantidad de nitrgeno y no se incluyen en
la mezcla los componentes agua o alcohol. Adems
las ecuaciones de estado cbicas PR y SRK han sido
probadas con este tipo de sistemas y han mostrado
una gran capacidad para reproducir los resultados
experimentales. Por lo tanto, es adecuado modelar este
tipo de mezclas con esta alternativa.
Conclusiones
En este trabajo se han utilizado las ecuaciones
de estado cbicas de SRK y PR combinadas con el
mtodo de Newton multivariable con el fin de calcular
los puntos crticos presentes en mezclas de gas natural.
El mtodo propuesto es eficiente y puede ser aplicado
para determinar las propiedades crticas de diversos
sistemas multicomponentes, es decir, no est limitado
simplemente a mezclas de gas natural.
La metodologa utilizada en este trabajo es
dependiente a las estimaciones iniciales de las
propiedades crticas. Por lo tanto, para converger al
punto crtico de la mezcla bajo anlisis es necesario
proponer diferentes formas de inicializar el mtodo.
Para los dos casos evaluados las Ecuaciones (5) (7)
resultaron adecuadas. No obstante, para mezclas
compuestas por ms componentes, con una cantidad
elevada de nitrgeno o constituida por agua o algn
alcohol se debern buscar formas alternativas
de inicializar el mtodo, debido a que en estas
condiciones la tcnica propuesta presenta problemas
de convergencia.
Para casos de mezclas de gas natural que contengan
agua o metanol ser necesario derivar las condiciones
crticas de Heidemann y Khalil (1980) para ecuaciones
de estado ms complejas con el fin de predecir de
manera adecuada las condiciones crticas, o utilizar
ecuaciones de estado cbicas con parmetros de
interaccin binaria en funcin de la temperatura.
Referencias
[1] Henderson, N., Freitas, L., Platt, M.G., (2004),
Prediction of critical points: a new methodology
using global optimization, A. I. Ch. E. Journal,
Vol. 50, p. 1300-1314.
[2] Freitas,
L., Platt, G., Henderson, N., (2004), Novel

Approach for the Calculation of Critical Points
in Binary Mixtures using Global Optimization,

Fluid Phase Equilibria, volumen 225, p.
29
-37.
[3] Goffe, W.L., Ferrier, G.D. y Rogers J., (1994),
Global optimization of statistical functions with
Simulated Annealing, J. Econometrics, volumen
60, p. 65-100.
[4]
Heidemann, R.A., Khalil, A.M., (1980), The
Calculation of Critical Points, A. I. Ch. E.
Journal, volumen 15, p.
769
-779.
[5]
Hicks, C.P., Young, C.L., (1977), Theoretical
Prediction of Phase Behavior at High Temperatures
and Pressures for Non-Polar Mixtures: I. Computer
Solution Techniques and Stability Tests, J. Chem.
Soc. Faraday II , volumen 73, p. 597-612.
[6] Justo-Garca, D.N., Garca-Snchez, F., (2005),
Clculo de puntos crticos de sistemas
multicomponentes utilizando optimizacin global,
XX Congreso Nacional de Termodinmica;
(Apizaco, Tlaxcala, 5-9 de Septiembre de 2005),
p. 342-366.
[7] Michelsen, M.L., (1980), Calculation of Phase
Envelopes and Critical Points for Multicomponent
mixtures, Fluid Phase Equilibria, volumen 4, p.
1-
10.
[8] Michelsen,
M.L., (1984), Calculation of Critical

Points and Phase Boundaries in the Critical
Region, Fluid Phase Equilibria, volumen
16,
p.
57
-76.
[9]
Michelsen, M.L., y Mollerup J., (1984), Calculation
of Critical Points and Phase Boundaries in the
Critical Region, Fluid Phase Equilibria, volumen
16,
p. 57
-76.
[10]
Michelsen, M.L., y Mollerup J., (2004),
Thermodynamic Modelling: Fundamentals and

Computational Aspects, Tie-Line Publications.
[11] PEMEX Anuario estadstico
(2005).
[12]
Peng, D.Y., Robinson, D.B., (1977), A Rigorous
method for predicting the critical propierties of
multicomponent systems from an equation of
state, A. I. Ch. E. Journal, volumen 23,
p. 137144.
[13] Snchez-Mares, F., Bonilla-Petriciolet, A. (2006),
Clculo de puntos crticos empleando una
estrategia de optimizacin global estocstica,
Afinidad, volumen, No. 525, p. 396-403.
[14] Stradi,
A.B., Brennecke, F.J., Khon, P.J., Stadtherr,

A.M. (2001), Reliable Computation of Mixture
Critical Points. A. I. Ch. E. Journal, volumen
47,
p. 212- 220.
Artculo recibido: 9 agosto 2007
Aceptado para publicacin: 9 noviembre 2007

36
Marco Antonio Ponce Guerrero, Daniel Salazar Sotelo

Mtodo integral para mejorar el diseo de columnas de destilacin, a partir de conceptos de ahorro de
energa
Revista del Centro de Investigacin. Universidad La Salle, vol. 4, nm. 16, enero-junio, 2001, pp. 69-76,
Universidad La Salle
Mxico
Revista del Centro de Investigacin. Universidad

La Salle,
evm@ulsa.mx
Universidad La Salle
Mxico
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
Mtodo integral para mejorar el

diseo de columnas de destilacin,
a partir de conceptos de ahorro de energa
MarcoAntonio Ponce Guerrero y Daniel Salazar Sotelo
Escuela de Ciencias Qumicas, ULSA.
<dsalazar @www.imp.mx>
RESUMEN
Se presenta una metodologa para analizar esquemas de destilacin y proponer modificaciones que signifiquen ahorros de energa, mediante el uso de curvas de entalpa-temperatura que relacionan la
operacin simulada con el caso de una operacin ideal reversible, tal como lo proponen Dhole y Linhoff.
Se hace una aplicacin a una columna de destilacin atmosfrica de crudo, que presenta un arreglo complejo con pumparounds y extracciones laterales. El resultado de adicionar un pumparound muestra una
reduccin en los consumos de energa, del orden de 400 millones de BTU/da, comprobando la utilidad
de los mtodos de anlisis de procesos que utilizan la simulacin.
Palabras clave: esquemas de destilacin, ahorros de energa, curvas de entalpa-temperatura.
ABSTRACT
A methodology for analysis of distillation schemes is presented. Improvements meaning energy savings
can be proposed, by using enthalpy-temperature curves relating simulated operations with ideal
reversible operations, according to Dhole and Linhoff.
The methodology is applied to an atmospheric distillation column, presenting a complex arrangement with
pumparounds and lateral extractions. The results show that by adding a pumparound is possible to save
an amount of energy around 400 millions BTU/day, confirming the advantages of analysis methods based
on simulation.
Keywords: distillation schemes, energy savings, enthalpy-temperature curves.
INTRODUCCiN
Los incrementos de precios de petrleo crudo

en el sector petrolero provocan aumentos en los
costos de los energticos y, en general, en los
costos de operacin de los procesos de transformacin industrial. Para abatir estos costos se
busca una reduccin en el consumo de energa
en la industria; esta tendencia se inici a finales
de la dcada de los aos setenta, forzando a los
investigadores a buscar mtodos innovadores
para lograr el ahorro y uso eficiente de la
energa.
Rev. Centro Inv. (Mx) Vol. 4, Nm. 16, Ene-Jun. 2001
El uso eficiente de la energa ha sido tomado

seriamente en nuestro pas donde se elabor
un Plan Nacional de Desarrollo 1995-2000, el
cual pretende que los bienes y servicios producidos por el sector industrial satisfagan los
niveles de calidad exigidos internacionalmente,
requirindose al mismo tiempo que sean compatibles con los objetivos de preservacin y
mejoramiento ambiental.
La industria de la refinacin
El desarrollo de la industria de la refinacin est
69
Artculo
ntimamente ligado al de la economa en general, debido a la importancia y variedad de productos que genera. En Mxico se cuenta con 6
grandes refineras que representan una capacidad total de procesamiento de 1.3 millones de
barriles por da (BPD) de crudo. PEMEX-Refinacin actualmente atraviesa un proceso de
modernizacin para incrementar la calidad de
los productos de refinacin y mejorar la eficiencia operativa de sus instalaciones.
Una refinera consta de diferentes plantas
que se utilizan para la transformacin del
petrleo crudo. La primera fase consiste fundamentalmente en el desalado del petrleo crudo,
el calentamiento en un horno y su envo a una
planta de destilacin primaria que opera a presin atmosfrica, donde se separa en: gas
hmedo, gasolina despuntada, gasolina primaria, nafta pesada, kerosina, gasleo ligero
primario (GOLP), gasleo pesado primario
(GOPP) y residuo primario. El residuo primario
se enva a la planta de destilacin de alto vaco,
'donde se separa en: gasleo ligero de alto
vaco (GOLAV), gasleo pesado de alto vaco

(GOPAV)y residuo de vaco.
La planta de destilacin combinada es una
planta de destilacin atmosfrica y de destilacin al vaco integradas trmicamente (ver
Figura 1). Existen dos trenes paralelos de precalentamiento de la carga para intercambiar
calor con los productos tanto de la destilacin
atmosfrica como de la de vaco; tambin se
intercambia calor con corrientes de extraccin
de calor de las torres.
En una refinera, estas unidades son las que
ms capacidad tienen y adems presentan
mayor consumo energtico. Es por ello que,
para este trabajo, se seleccion una columna
de destilacin atmosfrica como caso de estudio, para la aplicacin de un mtodo de anlisis
mediante el cual se pueden proponer modificaciones en el esquema de la columna que conlleven un ahorro de energa sin detrimento de la
calidad de los productos.
es,duo pnmano
Gasleo pesado
Turbosina
Queroseno
ni Me'
Gasolina
Torre de
Crudo
Residuo al!
vaco
Figura 1. Diagrama simplificado de una unidad de destilacin combinada.

70
Articulo
METODOLOGA
En la dcada de los ochenta las grandes empresas de ingeniera se abocaron al anlisis de
procesos en la industria, obteniendo finalmente
una tcnica til para el ahorro de energa; esta
tcnica de integracin energtica es llamada
"Mtodo de punto de pliegue" y es utilizada para
el anlisis de procesos nuevos (diseo), as
como para procesos existentes (rediseos o
revamps).
La tecnologa de punto de pliegue fue introducida por Linnhoff, y es una tcnica de anlisis
de la integracin trmica de un proceso qumico, en la que se demuestra que es posible el
establecimiento de diseo de procesos de intercambio trmico con un nmero mnimo de intercambiadores de calor y una mnima rea de
transferencia de calor.
Con el transcurso del tiempo, la aplicacin de
la tecnologa de punto de pliegue se ha extendido para analizar la integracin trmica de un
proceso qumico completo, incluso un complejo
de plantas. Linnhoff y Flower en 1978 propusieron el uso de correlaciones TemperaturaEntalpa e introdujeron la idea de fijar una diferencia de temperatura mnima DTmin para un
problema dado. En 1979, Linnhoff adicion el
concepto de flujo de calor que cruza el punto de
pliegue cuando desarroll lo que ahora se llaman las reglas de oro de la tecnologa de
pliegue. Esto fue seguido por la colocacin
apropiada del concepto de la integracin correcta de bombas de calor y columnas de destilacin.
Proceso de separacin por Destilacin
La destilacin es el proceso de separacin ms
extensamente usado en la industria qumica. Es
adems la operacin unitaria que ms consume
energa comparada con otros procesos de separacin. Esto hace a la destilacin un objetivo
primario para estudios de conservacin de la
energa.
Normalmente existe solamente un punto de
pliegue en el flujo de calor en cascada que
existe en las columnas de destilacin, y es un
punto con exactamente cero de flujo de calor. El
calor aceptado en el rehervidor es requerido a
una temperatura ms alta que la temperatura
del punto de pliegue y el calor cedido en el condensador est por debajo de la temperatura de
pliegue. El calor nunca debe de ser transferido
a travs del Punto de Pliegue.
Liebmann y Dohle aplicaron la tecnologa de
punto de pliegue a una columna de destilacin
fijando la mnima condicin termodinmica
(MTC). Para que una columna opere en la mnima condicin termodinmica, es necesario que
la curva de equilibrio coincida con la curva de
operacin, esto se podra lograr manteniendo el
mnimo flujo neto de calor en todo el equipo, por
medio de la adicin de una cantidad de calor
dQa en cada plato de la seccin de agotamiento
por medio de un rehervidor parcial y remocin
de una cantidad de calor dQr en cada plato de
la seccin de rectificacin. La condicin ideal de
una columna operando en la MTC requiere, por
lo tanto, de un nmero de etapas infinito y condensadores o rehervidores parciales por cada
etapa de equilibrio.
De una columna operando en esta condicin
se puede obtener el perfil de entalpas en funcin de la temperatura, para comparar los flujos
mnimos de calor con los de una torre real. Con
el objeto de obtener los flujos mnimos de calor
de una torre es necesario conocer los flujos
mnimos de vapor y lquido en cada punto y
asociarles su entalpa especfica toda vez que
esta propiedad slo es funcin de la temperatura, composicin y presin del sistema.
Para obtener la MTC de un sistema multicomponente, es necesario resolver "n" ecuaciones de balances por componente y "n" ecuaciones de equilibrio que son complicadas y
requieren de un procedimiento iterativo (donde
n es el nmero de etapas de equilibrio).
Generalmente los mtodos propuestos para
determinar la condicin mnima termodinmica
de sistemas multicomponentes son procedimientos tediosos, consumidores de tiempo, y
por lo mismo, poco prcticos. Dhole presenta un
mtodo para obtener lo que llaman la condicin
mnima termodinmica prctica (practical near
minimum thermodynamic condition o PNMTC),
basado en los resultados de una simulacin
convergida de una columna de destilacin, considerando un componente clave pesado y un
componente clave ligero.
71
Artculo
En el PNMTC, se resuelven simultneamente
las ecuaciones de equilibrio y operacin de los
componentes clave. Para ello, se incorporan las
composiciones de equilibrio de las corrientes
lquido y vapor que emergen de la misma etapa
en las ecuaciones de balance de materia:
Como se mencion anteriormente, el vapor y

las corrientes lquidas emergen desde la misma
etapa donde est en equilibrio una con otra. Las
entalpas para los flujos de lquido y vapor mnimos se obtienen mediante la proporcionalidad
molar directa:
Vmin'y~ - 4un 'x~ = DL
Vmin.y~ - 4nin . x~ =
Tr
H
Vmin
DH
= H*
Lmm
Estas ecuaciones establecen el flujo mnimo

de vapor (Vmin.) y el flujo mnimo del lquido
(Lmin.) a la temperatura de la etapa. L corresponde al clave ligero y H al clave pesado. D corresponde al flujo molar de los destilados. Estos
flujos, junto con las fracciones al equilibrio, nos
generan un sistema de dos ecuaciones con dos
incognitas; las incognitas son los flujos molares
de vapor y de lquido mnimos por cada etapa, y
se resuelven mediante mtodos numricos. Los
flujos mnimos de calor as obtenidos corresponden a la operacin reversible de una columna y sirven como criterio de comparacin contra
una columna real. La diferencia de valores entre
las. dos condiciones representan el alcance
mximo para mejorar las fuerzas impulsoras
dentro de la columna.
donde:* significa
Vmin
= flv ( V *)
L
Lmin
( L
*)
que estn al equilibrio.
L hace referencia a la etapa lquida.

V hace referencia a la etapa de vapor.
Al calcular los dficits de entalpa para cada
plato de la columna real y graficarlos en un diagrama de temperatura entalpa, se obtiene la
Gran Curva Compuesta (GCC) de la columna. A
partir de la informacin obtenida de los dficits
de entalpa, se obtienen las curvas compuestas
fra y caliente, que representan las fuerzas
impulsoras de transferencia de masa y calor
entre las etapas. (Ver Figura 2).
De acuerdo a Dhole el rea dentro de la Gran
Curva Compuesta representa las fuerzas impulsoras de transferencia de materia y de calor
dentro de la columna, mientras que las curvas
compuestas fra y caliente proyectan estas
fuerzas impulsoras entre las curvas lquido y
vapor.
Para obtener la representacin temperatura Entalpa para PNMTC, es necesario expresar

los flujos de lquido y vapor mnimos en funcin
de entalpas. La simulacin de la columna proporciona resultados etapa por etapa de los valores de entalpa de lquido y de vapor.
T
T,,,,
5'
Ti+2
Ti+]
Ti
T~d
Datos de la columna real y de la condicin

mnima termodiumica (MTC).
Construccin de la gran curva compuesta

GCC de la columna.
Figura 2. Diagramas de la Gran Curva Compuesta (GCC) para una columna de destilacin.
72
Artculo
ETAPA
T
REHERVIDOR
.......
Gran Curva Compuesta (T-H)
::
REHERVIDOR
............................................
Las modificaciones posibles en el diseo de

una columna de destilacin para reducir los
consumos de servicios auxiliares, como pueden
ser cambios en la presin de operacin, en la
relacin de reflujo, la adicin de cambiadores
de calor o bien de "pumparounds" (corrientes
de recirculacin), puede ser analizada con la
construccin de las curvas compuestas, donde
se vern fcilmente los niveles en la columna
que requieren modificaciones y el efecto de las
mismas sobre el balance de calor global. La
metodologa aplicada se puede resumir en los
siguientes pasos:
Simulacin rigurosa de la columna de destilacin, obteniendo perfiles de temperatura,

presin, flujo y entalpas, as como las composiciones de lquido y vapor en cada una de
las etapas de equilibrio.
Obtencin del flujo mnimo de vapor y lquido
con su Entalpa especfica asociada; elaborando un balance de materia en la columna
de destilacin para obtener el factor de separacin, el factor de enriquecimiento, el flujo
del destilado, el flujo del producto y obtener
as los flujos mnimos requeridos para las
secciones de rectificacin y agotamiento.
Elaboracin de la grfica del nmero de
platos reales contra el flujo de vapor real,
ideal y mnimo.
. Obtencin de la gran curva compuesta (diagrama entalpa-temperatura), y sus respectivas curvas caliente y fra.
Proposicin de nuevos casos de arreglos de
la columna de destilacin.
. Determinacin de los costosde inversinpor
::
Gran Curva Compuesta (Etapa-H)
concepto de las modificaciones y la disminucin de los costos de operacin de la columna, debido al ahorro de energa.
RESULTADOSY DISCUSiN
Para este estudio se tom como ejemplo la
planta combinada de la refinera de Tula, en
particular la columna atmosfrica que tiene una
capacidad de 150,000 BPD Y que consta de 34
platos tericos, tal como se muestra en la Figura 3.
El crudo desalado llega a las torres de
despunte D01A y D01B, donde se realiza la
primera separacin, obteniendo nafta de
despunte en el domo, la cual se manda a
tratamiento custico, y crudo despuntado en el
fondo. Por la parte superior de las torres de
despunte se alimenta un reflujo de nafta ligera,
con objeto de recuperar la nafta y fracciones
ms pesadas que van en el vapor de la alimentacin de la torre de despunte.
El crudo se manda a un horno donde se
calienta hasta 370C. En la columna de destilacin primaria se alimenta el crudo calentado y
vapor de agua, para disminuir la presin parcial
de los hidrocarburos y mantener la temperatura
y presin adecuadas en la zona de vaporizacin
de la torre, asi como para contribuir como
agente de arrastre de los productos ligeros. Se
obtienen diversos cortes como son los incondensables, las naftas ligeras y pesadas, la
kerosina, el GOLP, el GOPP y el residuo pri-
marioa diferentespresionesy temperaturas.

73
Articulo
La columna en su diseo original tiene tres agotadores laterales que salen de los platos 35, 39
Y 43; estos agotadores se identifican como
D03A, D038 y D03C obteniendo Nafta pesada,
Querosina y GOLP como productos.
rencia de masa entre estos platos aumentando

la eficacia en la separacin.
Al aplicar la metodologa y construir la Gran
Curva Compuesta (ver Figura nm. 4) , se
observa que en los platos 31 a 33 existe un
cambio abrupto en el consumo de energa, siendo que en estos puntos se genera energa que
puede ser aprovechada.
La torre de destilacin primaria tambien tiene

dos pumparounds que salen del plato 17 y del
plato 27, cuyo objetivo es aumentar la transfe-
8
SA14
3r.
~~
SAl5
6C
DOIA
DOIB
L----
9C
4A
ROlA
OJ3A
4C
-,
9A
ROIR
(JJ3B
4B
4D
~-.
9B
M1?T
--uJ
D02A I
IOD:
-,
m,
nA
,
jnR
731 MMBTU/d
13AG:
174MMBTU/d
12A:
495 MMBTU/d
85!C
I1A
14A
o1
(JJ7
- -
?lR
248 MMB1u/d
T;=121!C
T.=285!C
193 MMB-IU/d
T,,=180X:
T,=318!C I
ZiAiDO*_I
ISA
mR
~
L
21266 Mfv1BTU/d
355X:
1612 MMBTU/d
450!C
L--8170MMBTUld
347!C
Figura 3. Diagrama de la columna tomada como caso de estudio.

74
Rev. Centro Inv. (Mx) Vol. 4. Nm. 16, Ene-Jun. 2001
nnnn
__m-
Artculo
Etapa
-1000000
-500000
500000
10aoooo
1500000
200000O
2500000
H (MMBtuIh)
Figura 4. Gran Curva Compuesta para la condicin de diseo de la columna
En esta Figura 4, se observa la existencia de

2 puntos de pliegue por arriba de la alimentacin (platos 15 y 21) que no pueden ser
modificados debido a que son alimentacin a
los agotadores laterales; y un cruce del punto
de pliegue (plato 31) en la alimentacin donde
la entalpa de vapor es menor a la de lquido
debido a que la temperatura de alimentacin es
mucho mayor a la generada en el rehervidor.
Con estas observaciones se decidi generar
y simular una nueva alternativa como se menciona a continuacin:
Tratando de aplicar el mtodo de punto de
pliegue, se encontr que no existe posibilidad
para modificar la relacin de reflujo, ya que al
hacerlo
cambian
significativamente
las
propiedades y los rendimientos de los productos.
Modificar el precalentamiento, para extraer
calor de la alimentacin obliga a incrementar la
cantidad de vapor de agotamiento para obtener
una buena separacin.
Si se incorpora un tren de intercambio de
calor, a manera de condensador parcial, se

cambian fuertemente
productos destilados.
las propiedades
de los
Finalmente se opt por agregar un

pumparound; (ver Figura 5) esta corriente de
recirculacin se agreg del plato 28 al 30,
extrayendo un flujo de la corriente liquida y
pasndolo a una etapa superior, donde existe
una temperatura menor, recuperando calor y
bajando considerablemente el deficit de
entalpa en los platos inferiores a sta.
Con esta modificacin obtuvimos un deficit
de entalpa en las etapas donde se hizo la
modificacin de un orden de 3400 MMBTU/da,
comparado con el diseo orginal que tiene una
entalpa de orden de 3800 MMBTU/da, es
decir, un ahorro de energa de 400 MMBTU/da,
sin detrimento de la calidad y cantidad de los
productos destilados.
Con estos cambios se demuestra el potencial
del mtodo para estudiar modificaciones a
columnas de destilacin que bajan considerablemente el consumo de energa, que puede
aprovecharse en otros procesos, significando
un incremento
en la rentabilidadde las instala-
ciones.
75
Artculo
SA14
3B
13C
9C
D02A
21A
22A
(r.+
~
11~511b,jd
TE-121X:
T ,~2631C
72~5Ibm/d
TE-lSOX:
T,~3021C
:'-""""-"""'~"':
97l1lbm/d
::
:
TE=l92X:
T,=336X:
<t.
..d...J
23.B+
23A
DO:
::
..
Figura
5. Esquema modificado de la Columna
BIBLlOGRAFIA
(1) Liebmann K., Dohle R., "Integrated Crude

(2)
(3)
(4)
(5)
76
Distillation Design", Computers Chem.

Engng., vol. 19, pp. 8119-8124,1995.
Linhoff B., "Heat Integration of Distillation
Columns into Overall Processes", Chemical
engineering Science, vol. 38, nm. 8, pp.
1175-1188,1983.
Linhoff B., "Pinch Analysis - A State of the
Art Overview", Document of University of
Manchester Institute of 8cience & technology.503-522.
Chavez l., "Ahorro de Energa en columnas
de destilacin aplicando la Metodologa de
Punto de Pliegue", XXXV Convencin Nal.
delIMIQ, Morelia, Mich.
Dhole V. R., "Distillation Column Targets",
Computers Chem. Engng, vol. 17, nm.
5/6, pp. 549-560, 1993.
(6) Linhoff B., "General Proeess improvements

through Pinch Teehnology", Chemical engineering progress, pp 51-58, junio, 1998.
(7) Linhoff B., "Pinch Teehnology has come of
Age", Chemical engineering progress, pp.
34-40, julio, 1984.
(8) Fony Z., Thermodynamie Analysis Rectification 1. Reversible Model of rectification".
International chemical engineering, vol. 14,
nm. 1, pp. 18-27, enero, 1974.
(9) Fony Z., "ThermodynamicAnalysis Rectifieation 11.Finite aseade Models", International chemical engineering, vol. 14, nm.
2, pp. 203-210, abril,1974.
(10) Trivedi K. K., "Optimize a Licensor's Design
using Pinch Technology", Hydrocarbon processing, pp. 113-126, mayo, 1996.
Enrique Eduardo Tarifa, Eleonora Erdmann, Demetrio Humana, Samuel Franco Domnguez, Lorgio Mercado
Fuentes
Determinacin de la curva de destilacin flash para cortes de petrleo utilizando simuladores comerciales
Ingeniera e Investigacin, vol. 28, nm. 1, abril, 2008, pp. 50-56,
Universidad Nacional de Colombia
Colombia
Ingeniera e Investigacin,
revii_bog@unal.edu.co
Universidad Nacional de Colombia
Colombia
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
REVISTA INGENIERA E INVESTIGACIN VOL. 28 No. 1, ABRIL DE 2008 (50-56)
Determinacin de la curva de destilacin flash para cortes de

petrleo utilizando simuladores comerciales
Using commercial simulators for determining flash distillation curves for
petroleum fractions
Enrique Eduardo Tarifa1, Eleonora Erdmann2, Demetrio Humana3,
Samuel Franco Domnguez4 y Lorgio Mercado Fuentes5
RESUMEN
El presente trabajo describe un nuevo mtodo para estimar la curva de destilacin flash EFV (equilibrium flash vaporization) para cortes de petrleo utilizando simuladores comerciales. Para ello se implementa un modelo estacionario de la destilacin flash en un simulador comercial, y se ajusta el modelo utilizando una curva de destilacin obtenida con procedimientos estndares de laboratorio. Dicha curva puede ser del tipo TBP, ASTM D86, D1160 o
D2887, entre otras; e involucra un procedimiento experimental ms simple que el requerido para obtener la curva
EFV. Para la simulacin se puede utilizar cualquier simulador comercial que sea capaz de modelar petrleo, en este
trabajo se emplearon los simuladores HYSYS y CHEMCAD. Para evaluar el mtodo propuesto se analizaron experimentalmente varios tipos de petrleos y cortes de petrleo. Luego, de acuerdo con el mtodo propuesto, los datos
obtenidos fueron ingresados a un simulador para estimar las correspondientes curvas EFV. Las curvas estimadas
utilizando HYSYS y CHEMCAD fueron comparadas con las producidas por dos mtodos tradicionales de estimacin: el de Edmister y el mtodo de Maxwell. En todos los casos, las curvas estimadas por simulacin se aproximaron a la curva promedio de las de Edmister y Maxwell. El mtodo propuesto tiene varias ventajas: 1) evita la necesidad de obtener la curva EFV en forma experimental; 2) es independiente del tipo de curva experimental a utilizar
para ajustar el modelo; 3) permite realizar estimaciones para diferentes presiones utilizando una nica curva experimental como dato.
Palabras clave: petrleo, EVF, TBP, simulacin, HYSYS, CHEMCAD.

ABSTRACT
This work describes a new method for estimating the equilibrium flash vaporisation (EFV) distillation curve for petroleum fractions by using commercial simulators. A commercial simulator was used for implementing a stationary model for flash distillation; this model was adjusted by using a distillation curve obtained from standard laboratory analytical assays. Such curve can be one of many types (eg ASTM D86, D1160 or D2887) and involves an experimental
procedure simpler than that required for obtaining an EFV curve. Any commercial simulator able to model petroleum
can be used for the simulation (HYSYS and CHEMCAD simulators were used here). Several types of petroleum and
fractions were experimentally analysed for evaluating the proposed method; this data was then put into a process simulator (according to the proposed method) to estimate the corresponding EFV curves. HYSYS- and CHEMCADestimated curves were compared to those produced by two traditional estimation methods (Edmisters and Maswells
methods). Simulation-estimated curves were close to average Edmister and Maxwell curves in all cases. The proposed method has several advantages; it avoids the need for experimentally obtaining an EFV curve, it does not depend on the type of experimental curve used to fit the model and it enables estimating several pressures by using just
one experimental curve as data.
Keywords: petroleum, EVF, TBP, simulation, HYSYS, CHEMCAD.

Recibido: diciembre 7 de 2006
Aceptado: diciembre 18 de 2007
1
Dr., Ingeniera Qumica. Profesor asociado, Dedicacin exclusiva, Facultad de Ingeniera, Universidad Nacional de Jujuy (UNJu). Argentina.
Investigador adjunto, Consejo Nacional de Investigaciones Cientficas y Tcnicas (CONICET), Argentina. eetarifa@arnet.com.ar
2
Ingeniera qumica. Dra., Ciencia de Materiales. Profesor adjunto, Dedicacin exclusiva, Facultad de Ingeniera, Universidad Nacional de Salta
(UNSa), Argentina. Consejo de Investigaciones de la Universidad Nacional de Salta (CIUNSa) - Instituto de Investigaciones para la Industria Qumica
(INIQUI)CONICET, Argentina. eleonora@unsa.edu.ar
3
Ingeniero qumico. Jefe de trabajos prcticos, Dedicacin exclusiva, Facultad de Ingeniera, Universidad Nacional de Jujuy (UNJu), Argentina.
dhumana@arnet.com.ar
4
Ingeniero qumico. Profesor adjunto, Dedicacin exclusiva, Facultad de Ingeniera, Universidad Nacional de Jujuy (UNJu), Argentina. sfdominguez
@yahoo.com
5
Ingeniero qumico. Profesor titular, Dedicacin exclusiva, Facultad de Ingeniera, Universidad Nacional de Salta (UNSa), Argentina. Consejo de
Investigaciones de la Universidad Nacional de Salta (CIUNSa) - Instituto de Investigaciones para la Industria Qumica (INIQUI) CONICET,
Argentina. lmercad@unsa.edu.ar
50
TARIFA, ERDMANN, HUMANA, FRANCO, MERCADO
Introduccin
Para modelar un proceso que involucre petrleo es necesario
conocer o al menos estimar las propiedades del mismo. La
determinacin exacta de la composicin del crudo no es posible debido a la enorme cantidad de componentes que integran la mezcla. En su lugar, es prctica generalizada caracterizar el petrleo por medio de una curva de destilacin (temperatura frente a porcentaje acumulado de lquido destilado)
realizada de acuerdo a procedimientos estandarizados por
normas. Luego, se emplean correlaciones existentes para estimar las propiedades fsicas del petrleo a partir de las temperaturas que conforman la curva de destilacin (Behrenbruch
y Dedigama, 2007). En base a lo anterior, el crudo puede ser
modelado como un grupo de componentes hipotticos cuyas
propiedades dan forma a la curva de destilacin. La incorporacin de mediciones experimentales adicionales, como el
peso molecular, la densidad, o la viscosidad en funcin del
porcentaje acumulado de lquido destilado, generan nuevas
curvas que permiten refinar las estimaciones para el crudo
(Riazi et al., 2004). Por lo general, el peso molecular y la
densidad se miden para el volumen total de la muestra.
Por otra parte, los componentes livianos (light ends) de la
muestra son componentes orgnicos de bajo peso molecular,
como por ejemplo: metano, etano, propano, etc. Los componentes livianos pueden ser estimados a partir de la curva
de destilacin, o de un anlisis por separado, o pueden ser
ignorados, dependiendo ello de los objetivos del estudio. Si
se incorporan al estudio los componentes livianos, la cantidad de estos representa la cantidad de crudo que tiene un
punto de ebullicin promedio igual al punto especificado para dichos componentes. Esto es, la base elegida para los componentes livianos es la misma que la tomada para la muestra
de crudo total.
La conveniencia de estimar un tipo de curva de destilacin a
partir de otra ha sido reconocida desde hace tiempo (Hoffman, 1969). Lo mismo ocurre con la estimacin de las curvas
de destilacin a partir de otros tipos de datos experimentales
(Mondragn y Ouchi, 1984; Falla et al., 2006; Pasquini y Ferreira Bueno, 2007). En este trabajo se presenta un nuevo
mtodo para estimar la curva de destilacin EFV (equilibrium
flash vaporization) a partir de curvas obtenidas por procedimientos estndares de laboratorio, tales como: TBP (true boiling point distillation), ASTM D86, la destilacin D1160 (separada o combinada), y la destilacin D2887 simulada a partir de datos de cromatografa (Watkins, 1981). Por definicin,
la curva EFV es generada por una serie de experimentos realizados a presin constante (por lo general, 1 atm). El resultado vincula la temperatura de equilibrio con el porcentaje
de volumen de lquido destilado, donde el vapor total est
en equilibrio con el lquido que no se evapor. Esta curva es
ampliamente utilizada en el diseo de equipos para la industria petroqumica (Cerutti, 2002). Sin embargo, debido a que
el procedimiento experimental requerido para obtener la
curva EFV es complicado, esta es generalmente estimada a
partir de otras curvas que involucran procedimientos experi-
mentales ms simples. Para tal fin, se desarrollaron mtodos

para estimar la curva EFV, siendo el de Edmister y el de Maxwell los ms utilizados (Maxwell, 1950; Nelson, 1958; Edmister and Okamoto, 1959). Estos mtodos realizan la estimacin empleando un procedimiento grfico que requiere datos obtenidos de diagramas y tablas desarrollados para dicho
efecto. Este empleo intensivo de grficos y tablas, hace que
los mtodos sean engorrosos, pocos exactos y difciles de seguir sin cometer errores.
Para resolver los problemas que afectan a los mtodos anteriormente descritos, en este trabajo se presenta un nuevo
mtodo para estimar la curva EFV utilizando simuladores de
procesos qumicos. El mtodo se ejemplifica utilizando los
simuladores comerciales HYSYS y CHEMCAD; pero igualmente puede utilizarse cualquier otro simulador, siempre y
cuando sea capaz de manejar petrleo. A fin de comparar el
mtodo propuesto con los de Edmister y Maxwell, se analizaron muestras de varios tipos de petrleo. En todos los casos
estudiados, la curva EFV estimada por simulacin se ubic
entre las estimadas por los otros dos mtodos. Una gran ventaja del nuevo mtodo es que no depende del tipo de curva
experimental empleada; esto es, la curva EFV puede ser estimada siguiendo el mismo procedimiento sin importar el tipo
de curva de destilacin experimental ingresada al simulador
para caracterizar al petrleo.
Caracterizacin del petrleo

La destilacin de petrleo en laboratorio es de gran utilidad
en el refinamiento y formulacin de combustibles. Tres de
los ensayos ms importantes son (Wuithier, 1971; Speight,
1998; Cohen, 2003): la curva ASTM (la ms difundida), la
TBP y la EFV. En las secciones siguientes se describen los citados ensayos.
La curva ASTM
Para este ensayo se utiliza el baln Engler con dimensiones
estndar. Se le debe suministrar calor a una velocidad tal que
pueda recogerse 5 cm3 por minuto de destilado. Durante el
experimento se registra para la curva la temperatura a la cual
se obtiene la primera gota de destilado; lo mismo se hace
cada vez que se recoge un nuevo 10% de destilado. La temperatura mxima es 370 C, por lo cual se trata de una destilacin suave sin divisin.
La curva TBP
La destilacin TBP no est estandarizada. El destilado se obtiene de la divisin de la muestra en las fracciones que la
componen. Se emplea el equipo Geters con una carga de
1.000 a 5.000 cm3 de muestra. La columna de fraccionamiento est rodeada por un calentador que compensa las
prdidas de calor, de esta forma se pretende emular una columna adiabtica. En el tope de la columna se instala un condensador que condensa el vapor que llega a l. El condensado desciende y entra en contacto con el vapor que saciende, all se produce la separacin. Nuevamente, la curva se
51
DETERMINACIN DE LA CURVA DE DESTILACIN FLASH PARA CORTES DE PETRLEO UTILIZANDO SIMULADORES COMERCIALES
obtiene registrando la temperatura en el momento en que se

produce la primera gota de destilado, y cada vez que se recoge un nuevo 10% de este.
La curva EFV
Esta curva grafica la temperatura de equilibrio lquido-vapor
frente al porcentaje de volumen de lquido destilado. La
Figura 1 muestra el equipo de laboratorio empleado para la
destilacin flash, el mismo est simplificado para facilitar la simulacin; sin embargo, la modificacin realizada no afecta la
exactitud de los resultados. El tanque T-1 contiene el petrleo a ser analizado en condiciones estndar (15 C, 1 atm). El
lquido a destilar pasa por el calentador H-1 para alcanzar la
temperatura de equilibrio. Cuando el lquido calentado alcanza la torre F-1 se produce la evaporacin parcial a la presin presente en dicho equipo (generalmente 1 atm). El lquido que no evapor es recogido en el tanque T-2, mientras
que el vapor es condensado y enfriado a 15 C en el intercambiador C-1, y finalmente es recogido en el tanque T-3. El
porcentaje de lquido destilado es: 100(volumen del lquido
en T-3 al final)/(volumen del lquido en T-1 en el inicio); o lo
que es mejor, de la siguiente manera: 100(flujo volumtrico
de FVL)/(flujo volumtrico de F0). De esta forma, repitiendo
el experimento para distintas temperaturas de equilibrio, y
graficndolas frente a los correspondientes porcentajes de
volumen destilado, se obtiene la curva EFV.
dos y aproximados. En las secciones siguientes se describirn

dos de los mtodos ms difundidos.
Mtodo de Maxwell
El mtodo de Maxwell utiliza una lnea recta auxiliar que corta la curva TBP en 10% y 70% del porcentaje de volumen
destilado; esta lnea se llama DRL (distillation reference line).
Luego, se define la lnea FRL (flash reference line) en funcin
de la lnea anterior y de datos extrados de diagramas desarrollados para tal fin. Finalmente, la curva EFV se obtiene a
partir de la FRL y datos provenientes de diagramas adicionales (Maxwell, 1950). La intensa utilizacin de diagramas hace
que este mtodo sea complejo e inexacto.
Mtodo de Edmister
El mtodo de Edmister localiza un punto clave (la temperatura correspondiente al 50% de volumen destilado). Luego, estima los incrementos de la curva EFV en funcin de los aumentos de la TBP utilizando datos extrados de diagramas
(Edmister y Okamoto, 1959). Finalmente, la curva EFV se obtiene procesando todos los incrementos. Una vez ms, el intenso uso de diagramas es la mayor debilidad del mtodo.
Correccin por presin

La curva EFV estimada para 1 atm por cualquiera de los dos
mtodos descritos anteriormente puede ser utilizada para estimar una nueva curva EFV para una presin diferente. Esto
puede hacerse mediante el trazado de una curva paralela cuya posicin est determinada por el diagrama de Cox y la
nueva presin.
Estimacin de la curva EFV por simulacin
Figura 1. Equipo para la destilacin flash.
Estimacin de la curva EFV

Conocer las condiciones del equilibrio de fases vapor-lquido
de las fracciones del petrleo es esencial para el diseo de la
mayora de los equipos empleados en la industria petroqumica. Esta informacin puede ser obtenida experimentalmente mediante la determinacin de la curva EFV. Sin embargo, el procedimiento es laborioso y costoso. Por estos motivos, generalmente se emplean correlaciones empricas para
estimar la curva EFV a partir de curvas ASTM o TBP, que son
ms fciles de obtener experimentalmente.
En la literatura se presentan varias correlaciones que involucran las curvas ASTM, TBP y EFV a presin atmosfrica. No
obstante, ninguno de esos mtodos es aplicable para todo el
rango de presin y temperatura existente en la industria del
petrleo. Ms an, los mtodos de estimacin son complica-
52
Los simuladores de procesos estn cada vez ms difundidos,

y habiendo discutido los inconvenientes de los mtodos tradicionales de estimacin de la curva EFV, es conveniente desarrollar un nuevo mtodo de estimacin que emplee simulacin con todas las ventajas que ello implica: exactitud, flexibilidad, velocidad, costo reducido, etc. A continuacin se
describe cmo construir un modelo para el equipo de la
Figura 1. Ese modelo puede ser implementado en cualquier
simulador que sea capaz de manejar petrleo, por ejemplo,
CHEMCAD. Luego, un segundo modelo es presentado, el
cual es ms simple de simular pero ms complejo conceptualmente. El simulador HYSYS puede manejar ambos modelos.
Simulacin del equipo de destilacin con HYSYS

La Figura 2 muestra el PDF (Process Flow Diagram) del modelo implementado en HYSYS del equipo experimental presentado en la Figura 1. Se simula slo la parte que comienza
en la corriente F0. El paquete fsico-qumico (fluid package)
recomendado es Peng Robinson. La Tabla 1 detalla los objetos y datos requeridos para la simulacin. La variacin de la
temperatura se realiza utilizando la herramienta Case Studies.
Los pasos para introducir la composicin del petrleo se
describen en la siguiente seccin, donde todo el PDF de la

Figura 2 es reemplazado por una nica corriente llamada FT.
3.
4.
Figura 2. Equipo flash de laboratorio simulado con HYSYS.
Tabla 1. Unidades de operacin y corrientes empleadas en la simulacin
con HYSYS.
Nombre
Tipo
Datos
Heater
H-1
Delta P = 0
Separator
F-1
Delta P = 0
Cooler
C-1
Delta P = 0
Energy Stream
EH
Energy Stream
EC
Material Stream
F0
Composicin del petrleo,
T = 15 C, P = 1 atm,
Liq Vol Flow @Std Cond (m3/h) =
100
Material Stream
F1
T = temperatura de equilibrio
deseada
Material Stream
FV
Material Stream
FL
Material Stream
FVL
T = 15 C
5.
Simulacin de una corriente con HYSYS

Aunque el modelo de la seccin anterior es correcto, los mismos resultados pueden ser obtenidos utilizando un modelo
ms simple que emplea caractersticas especiales de HYSYS.
El nuevo modelo est formado por una nica corriente. El
procedimiento para estimar la curva EFV con este nuevo modelo es:
1. Cree un nuevo caso en HYSYS.
2. En el PFD, inserte una corriente material llamada FT que
contenga la composicin del petrleo. Para ello, siga
los siguientes pasos:
2.1. En la solapa Fluid Pkgs, seleccione el paquete Pen
Robinson y agregue los componentes pertenecientes a los Light Ends.
2.2. En la solapa Oil Manager, presione el botn Enter
Oil Environment.
2.3. En la solapa Assay, defina un nuevo Assay introduciendo la curva experimental que posee (TBP,
ASTM, o cualquiera de las que acepta HYSYS).
Ingrese all cualquier dato adicional que disponga
(propiedades globales, curvas de densidad, etc.).
Presione el botn Calculate para determinar las
correspondientes Working Curves. Cierre la ventana Assay.
2.4. Para el recin definido Assay, en la solapa Cut/Blend
defina un nuevo Blend. Cierre la ventana Blend.
6.
2.5. Finalmente, en la solapa Install Oil instale una corriente material, y nmbrela FT. Presione el botn
Return Basic Environment.
2.6. Presione el botn Return to Simulation Environment.
En la corriente FT, especifique las siguientes condiciones:
P=1 atm, T= temperatura inicial de la curva empleada
en Assay, Liq Vol Flor@Std Cond = 100 m3/h (base de
clculo). En este momento, HYSYS puede resolver la
corriente, y esta se vuelve azul.
Para generar todos los puntos de la curva automticamente, utilice la herramienta Databook. Para ello, siga
los pasos que se detallan a continuacin:
4.1. En la solapa Variables, agregue las siguientes variables pertenecientes a FT: Temperature y Phase Liquid Volume Flow STD (Vapour Phase).
4.2. En la solapa Case Studies, agregue un nuevo caso
llamado Curva EFV. Seleccione Temperature como variable independiente, y Phase Liquid Volume Flow STD (Vapour Phase) como variable dependiente.
4.3. Presione el botn View y especifique el rango deseado para Temperature.
4.4. Presione el botn Start.
4.5. Una vez finalizado el clculo, presione el botn Results, seleccione Table. Copie la tabla (seleccinela y presione ctrl.+C), o seleccione File/Print para generar un archivo de texto.
Pegue o importe los datos en una hoja de clculo (por
ejemplo, Excel).
Grafique la curva EFV.
Simulacin con CHEMCAD

El procedimiento completo para crear la curva EFV utilizando
el simulador CHEMCAD es el siguiente:
1. Comience un nuevo trabajo.
2. Cree el flowsheet mostrado en la Figura 3 y emplee los
da-tos de la Tabla 3.
3. Seleccione los compuestos livianos utilizando Component
List.
4. Seleccione Peng-Robinson para K-value y Enthalpy.
5. Defina la corriente 1 ingresando la curva de destilacin
como sigue:
5.1. En la opcin ThermoPhysical, seleccione Distillation
Curve. La primera vez que ejecute el comando, el
programa mostrar una pantalla preguntando por el
ID de la corriente a ser caracterizada. Ingrese 1 para
definir la corriente 1.
5.2. Seleccione la correlacin a usar para la caracterizacin de la curva de destilacin y para la estimacin de las propiedades de los pseudo componentes. Seleccione Chemstations para Molecular
Weight Method y para ASTM D86-TBP Interconversion Method, y seleccione Lee Kesler para Crit.cal Properties Method.
53
5.3. En el siguiente formulario (Curve Temperature Cut

Range), en la primera fila ingrese la temperatura
inicial y la final de la curva experimental, y 50 para No of_points.
5.4. Seleccione Edit Curve Data para ingresar las propiedades globales. Seleccione 0) User specified for
Distillation Curve, y seleccione el tipo apropiado
para la curva experimental (TBP, ASTM, etc.). Si
no conoce la curva de densidad, ingrese la correspondiente Bulk gravity. Seleccione 4) m3/h para
Selected total flow units, e ingrese 100 para Total
flow rate (base de clculo).
5.5. Ingrese la curva experimental. Opcionalmente, puede ingresar tambin la curva de densidad.
5.6. Presione Save and Exit para iniciar el clculo.
6. En la corriente 1, especifique las siguientes condiciones:
P=1 atm y T=15C. En este punto, CHEMCAD puede
resolver la corriente, para ello presione el botn Flash.
7. En el equipo heater 1, ingrese la temperatura inicial de la
curva experimental en Temperature of stream 2.
8. En el equipo cooler 2, ingrese 15 C en Temperature of
stream 4.
9. Resuelva el modelo presionando Run all.
10. Para generar automticamente todos los puntos de la
curva EFV, utilice la herramienta Sensitivity Study. Para
ello, proceda de la siguiente manera:
10.1. En el men Run, seleccione la opcin Sensitivity Study.
10.2. Seleccin la opcin New Analysis.
10.3. Ingrese el nombre Curva EFV para el anlisis de
sensibilidad.
10.4. Presione Edit Independent Variable, seleccione la
opcin Equipment, especifique 1 para ID, seleccione 4 T Out Srt 1 for Variable, escriba Teq en
Variable name, e ingrese la temperatura inicial y final de la curva experimental en Vary this variable
from/to, y 50 en equal steps.
10.5. Presione Edit Recorded Variables, seleccione Stream,
especifique 4 en ID, seleccione 7 Total std liq. as
Variable, escriba Vol % en Variable name y seleccione 29 Actual liq vol rate para Variable units.
10.6. Presione Run All. CHEMCAD genera la curva EFV.
10.7. Presione Plot Results, seleccione Teq para el eje Y,
y Vol % para el eje X.
11. Pegue o importe los datos a una planilla de clculo
(por ejemplo, Excel).
12. Nuevamente grafique la curva EFV.
Tabla 2. Unidades de operacin empleadas en la simulacin con
CHEMCAD.
Nombre
1
2
3
54
Tipo
Heater
Cooler
Separator
Figura 3. Equipo flash de laboratorio simulado con CHEMCAD.

Tabla 3. Corrientes utilizadas en la simulacin con CHEMCAD.
Nombre
Tipo
Datos
Material Stream
1
Composicin del petrleo,
T = 15 C, P = 1 atm,
Liq Vol Flow @Std Cond (m3/h) = 100
Material Stream
2
T = temperatura de equilibrio
deseada
Material Stream
3
Material Stream
4
T = 15 C
Material Stream
5
Casos de estudio
Para evaluar el mtodo propuesto, se estudiaron varios tipos
de petrleo. La Tabla 4 caracteriza las muestras estudiadas,
donde puede verse que se emple una amplia variedad de
muestras. En los primeros tres ejemplos, se utiliz la curva
TBP para caracterizar la muestra, mientras que en los dos ltimos ejemplos se emple la curva ASTM D86. Todas ellas
fueron obtenidas experimentalmente. Para todos los ejemplos, se estim la curva EFV utilizando los mtodos de Edmister y Maxwell, como as tambin los simuladores HYSYS y
CHEMCAD. Las curvas obtenidas para los ejemplos son
representadas desde la Figura 4 a la Figura 8. Para todos los
ejemplos, las curvas EFV obtenida por simulacin se ubicaron
entre las obtenidas por los mtodos de Edmister y Maxwell.
Sin embargo, en el ejemplo 2, donde se emple un petrleo
pesado, la curva obtenida con CHEMCAD muestra una desviacin hacia abajo para porcentaje de destilados superiores
a 50%. Finalmente, la Figura 9 muestra las curvas EFV para el
ejemplo 1 pero a una presin distinta de 1 atm. Para ello se
realiz la correspondiente correccin por presin en los mtodos tradicionales; mientras que para obtener las curvas
simuladas, simplemente se ingres la nueva presin P=50
mmHg en el modelo, y ningn otro cambio fue necesario.
Nuevamente, las curvas producidas por HYSYS y CHEMCAD
estn entre las producidas por los mtodos de Edmister y
Maxwell.
Tabla 4. Muestras analizadas.
Ejemplo
Curva
experimental
1
2
3
4
5
TBP
TBP
TBP
ASTM D86
ASTM D86
Datos
Delta P = 0
Delta P = 0
Rango de
temperatura
(C)
25-289
61-461
58-88
71-256
56-358
15/15
Comentario
0.80
0.84
0.76
0.75
0.75
Crudo
Crudo
Nafta
Crudo
Crudo
Conclusiones
350
300
T (C)
250
200
150
100
50
0
0
10
20
30
40
50
60
70
80
90
100
V%
TBP
Maxwell
Edmister
Hysys
Chemcad
Figura 4. Ejemplo 1
450
400
350
T (C)
300
250
200
150
100
50
0
0
10
20
30
40
50
60
70
80
V%
TBP
Maxwell
Edmister
Hysys
Chemcad
Figura 5. Ejemplo 2
100
Un nuevo mtodo para estimar la curva EFV

por simulacin fue presentado. Para evaluar el
mtodo se emplearon los simuladores HYSYS
y CHEMCAD, pero cualquier otro simulador
con capacidad para manejar petrleo puede
ser utilizado. Si bien el mtodo presentado exige que se realice un esfuerzo para construir
el modelo, este esfuerzo se realiza una nica
vez; y una vez disponible, el modelo permite
estimar las curvas EFV que sean necesarias, ahorrando todo el trabajo que implica obtener
en forma experimental dichas curvas. Por otra
parte, es cierto que an es necesario realizar
experiencias para obtener las curvas que se utilizan para ajustar el modelo, pero las experiencias que se deben llevar a cabo son mucho ms simples que las requeridas para producir la curva EFV.
Las curvas estimadas utilizando HYSYS y
CHEMCAD fueron comparadas con aquellas
producidas por dos mtodos tradicionales de
estimacin: el mtodo de Edmister y el mtodo de Maxwell. En todos los casos, las curvas
estimadas por simulacin se aproximaron a la
curva promedio de las curvas de Edmister y
Maxwell. Aunque los resultados son comparables, el mtodo propuesto tiene las ventajas
que provienen de automatizar el clculo, como por ejemplo: mayor velocidad, ms exactitud, menor probabilidad de cometer errores,
etc. A esto se suma la gran variedad de curvas
que pueden ser utilizadas por el modelo dependiendo del simulador utilizado, siendo los
pasos a seguir independientes del tipo de curva utilizada. Finalmente, el mtodo permite
estimar fcilmente la curva EFV a diferentes
presiones utilizando una nica curva como
dato; para ello, basta con indicar en el modelo cul es la presin deseada.
90
Bibliografa
80
70
Behrenbruch, P., Dedigama, T., Classification

and characterisation of crude oils based on
distillation properties., Journal of Petroleum
Science and Engineering., 57, pp. 166180,
2007.
T (C)
60
50
40
30
Cerutti, A. A., La Refinacin del petrleo., Tomo I, Ed. IAPG, pp. 9-25, 2002.
20
10
0
0
10
20
30
40
50
60
70
80
90
V%
TBP
Figura 6. Ejemplo 3
Maxwell
Edmister
Hysys
Chemcad
100
Cohen, L., Diseo y Simulacin de Procesos

Qumicos., Segunda Ed., Algeciras-Espaa,
Editor Len Cohen Mesonero, 2003.
Edmister, W. C., Okamoto, K. K., Applied Hydrocarbon
Thermodynamics-Part
12:
55
Falla, F. S., Larini, C., Le Roux, G. A. C., Quina, F. H., Moro,

L. F. L., Nascimento, C. A. O., Characterization of crude
petroleum by NIR., Journal of Petroleum
Science and Engi-neering, 51, pp. 127137,
2006.
Equilibrium Flash Vapo-rization Correlations for Petroleum

Fractions, August- Petro-leum Refiner, 1959.
300
Hoffman, E. J., Relations between true boiling

point and ASTM distillation curves., Chemical
Engineering Science, 24 (1), January, pp. 113
-117, 1969.
250
T (C)
200
150
Maxwell, J. B., Data Book on Hydrocarbons.

London, D. van Nostrand Company, 1950.
100
Mondragon, F., Ouchi, K., New method for

obtaining the distillation curves of petroleum
products and coal-derived liquids using a
small amount of sample., Fuel, 63 (1), January, pp. 61-65, 1984.
50
0
0
10
20
30
40
50
60
70
80
90
100
V%
ASTM D86
Maxwell
Edmister
Hysys
Chemcad
Figura 7. Ejemplo 4
400
350
T (C)
250
200
150
Speight, J. G., The Chemistry and Technology

of Petroleum., Third Ed, New York, 1998.
100
Watkins, R. N., Petroleum Refinery Distillation.,

Second Ed., Gulf Publishing, 1981.
50
0
0
10
20
30
40
50
60
70
80
90
100
V%
ASTM D86
Maxwell
Edmister
Hysys
Chemcad
Figura 8. Ejemplo 5
350
300
250
T (C)
Pasquini, C., Ferreira Bueno, A., Characterization of petroleum using near-infrared spectroscopy: Quantitative modeling for the true
boiling point curve and specific gravity., Fuel,
86, pp. 19271934, 2007.
Riazi, M. R., Al-Adwani, H. A., Bishara, A., The
impact of characterization methods on properties of reservoir., fluids and crude oils: options
and restrictions, Journal of Petroleum Science
and Engineering, 42, pp. 195 207, 2004.
300
200
150
100
50
0
0
10
20
30
40
50
60
70
80
90
V%
TBP
Maxwell
Edmister
Hysys
Chemcad
Figura 9. Correccin por presin, P = 50 mmHg
56
Nelson, W. L., Petroleum Refinery Engineering., Ed. Mc Graw Hill, 1958.
100
Wuithier, P., Refino y Tratamiento Qumico: el

petrleo., Madrid, Cepsa S. A, Marcel Dekker
Inc., 1971.
A. Castro Agero, A. Jimnez Gutirrez

Mtodo corto para el diseo de columnas de destilacin trmicamente acoplada petlyuk
Revista Mexicana de Ingeniera Qumica, vol. 1, nm. 3, 2002, pp. 85-96,
Universidad Autnoma Metropolitana Unidad Iztapalapa
Mxico
Revista Mexicana de Ingeniera Qumica,

amidiq@xanum.uam.mx
Universidad Autnoma Metropolitana Unidad
Iztapalapa
Mxico
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
REVISTA MEXICANA DE INGENIERA QUMICA Vol. 1 (2002) 85-96
AMIDIQ
METODO CORTO PARA EL DISEO DE COLUMNAS DE DESTILACION

TERMICAMENTE ACOPLADA PETLYUK
A SHORT CUT METHOD FOR THE DESIGN OF PETLYUK THERMALLY COUPLED
DISTILLATION COLUMNS
A. Castro-Agero* y A. Jimnez-Gutirrez.
Depto. Ingeniera Qumica, Instituto Tecnolgico de Celaya.
Av. Tecnolgico y Garca Cubas S/N, 38010, Celaya, Gto.
Resumen
Se propone un mtodo de diseo para el sistema de destilacin trmicamente acoplado, mejor conocido como
sistema Petlyuk. El mtodo se basa en las ecuaciones de Fenske, Underwood y Gilliland, con las cuales se puede
tener un buen estimado de los flujos internos, nmero de etapas tericas, as como las zonas o etapas de
interconexin entre el prefraccionador y la columna principal. Para mostrar la aplicacin del mtodo de diseo, se
incluye el caso de estudio de separacin de una mezcla ternaria.
Palabras clave: destilacin, petlyuk, destilacin trmicamente acoplada.
Abstract
A method for the design of the Fully Thermally Coupled Distillation system, or Petlyuk column, is presented. The
method is based on the Fenske, Underwood and Gilliland equations, which provide the internal flowrates, number of
ideal stages, as well as the stages where the prefractionator and the main column must be interlinked. A case study
for the separation of a ternary mixture is presented to show the application of the design method.
Key words: distillation, petlyuk column, thermally coupled distillation
1. Introduccin
Los sistemas convencionales de
destilacin son relativamente fciles de
disear y operar, pero tienen un consumo de
energa significativo. En los ltimos 20 aos,
se han hecho estudios sobre los sistemas de
destilacin para mejorar su consumo de
energa (Finn, 1993). Se han desarrollado
secuencias no convencionales como las
columnas con corrientes laterales y los
sistemas trmicamente acoplados, entre
otros. Dentro de los sistemas acoplados, se ha
demostrado que el sistema denominado
Petlyuk puede ahorrar hasta un 30% del costo
de energa comparado con los sistemas
convencionales (Tedder y Rudd, 1978;
Glinos y Malone, 1985). Una de las razones
por la que se ha limitado el uso de las
columnas trmicamente acopladas es la falta
de mtodos confiables de diseo.
En este campo los esfuerzos se han

concentrado en utilizar la ecuacin de
Underwood (1948) para comparar los
diferentes arreglos desde el punto de vista de
consumo de energa.
Tedder y Rudd (1978) hicieron un
estudio paramtrico de ocho sistemas de
destilacin para separar mezclas ternarias.
Dentro de estos sistemas se encuentran
incluidos las secuencias convencionales
(Directa e Indirecta) y tres sistemas acoplados
(rectificador lateral, agotador lateral y sistema
con prefraccionador). Para este estudio
utilizaron mtodos de diseo basados en
Nodos de Composicin, los cuales hacen un
leve sobrediseo de la columna para
separaciones perfectas. Este sobrediseo se
incrementa cuando la separacin es ms
imperfecta.
Autor para la correspondencia. E-mail: arturo@iqcelaya.itc.mx

Tel. (46161) 17801 Fax: 17575
Publicado por la Academia Mexicana de Investigacin y Docencia en Ingeniera Qumica, A. C.
85
Castro-Agero y Jimnez-Gutirrez / Revista Mexicana de Ingeniera Qumica Vol. 1 (2002) 85-96
Fidkowski y Krlikowski (1986, 1987

y 1990) desarrollaron un procedimiento para
obtener el flujo de vapor ptimo para
sistemas trmicamente acoplados, incluido el
sistema Petlyuk. Para calcular los flujos
internos utilizan la ecuacin de Underwood
(1948). El procedimiento se presenta para
una solucin ideal ternaria, es decir, valores
constantes de volatilidades relativas y
velocidades de flujos internos equimolares, y
separaciones perfectas. El estudio se ve
limitado a sistemas con alimentacin como
lquido saturado.
Glinos y Malone (1985) desarrollaron
un procedimiento simplificado para calcular
los flujos mnimos de vapor y lquido dentro
de una columna de destilacin con
rectificador lateral. En este estudio, la
columna con rectificador lateral se
descompone en un arreglo equivalente (en
cuanto a flujos) de dos columnas, con el
destilado de la primera sirviendo como
alimentacin a la segunda columna. Glinos y
Malone muestran que bajo ciertas
condiciones de operacin, los flujos que
sirven de interconexin entre la primera y la
segunda columna se pueden considerar como
una seudo-alimentacin con una condicin
trmica desplazada.
El sistema Petlyuk puede ser
reproducido en una sola torre de destilacin
con una divisin interna. Este sistema se
denomina de pared dividida. Para que el
sistema Petlyuk y el de pared dividida sean
termodinmicamente equivalentes, no debe
de haber transferencia de energa a travs de
la pared divisora.
El presente trabajo propone un mtodo
simplificado para el diseo del sistema
Petlyuk o de pared dividida. La importancia
de este mtodo es que proporciona una
primera aproximacin de manera rpida, la
cual puede luego validarse mediante
simulaciones rigurosas.
86
2. Teora
Cuando se trata de separar una mezcla
ternaria o de multicomponentes en tres
productos, normalmente se usan dos torres de
destilacin
(secuencias
convencionales
directa e indirecta).
En el caso del sistema Petlyuk (Fig. 1a)
o su equivalente de pared dividida (Fig. 1b),
la forma de separar los componentes es
diferente a los sistemas convencionales. La
alimentacin
se
introduce
a
un
prefraccionador o seccin I, el cual hace una
separacin inicial entre el componente ms
voltil y el ms pesado. La composicin del
componente intermedio se distribuye entre el
domo y el fondo de la seccin I.
La columna principal se puede dividir
en dos secciones, la seccin superior a la
corriente de producto S o seccin II de las
Fig. 1a y 1b, y la seccin inferior a la
corriente S, o seccin III. La seccin II tiene
una alimentacin de vapor saturado y una
corriente lateral lquida que es el reflujo de la
seccin I, mientras que la seccin III tiene
una alimentacin de lquido saturado y una
corriente lateral en forma de vapor.
Cada una de las secciones de la
columna principal realiza una separacin
entre compuestos adyacentes. Adems, si la
separacin que se realiza en la seccin I es lo
ms perfecta posible, se puede considerar
como una separacin binaria. Finalmente lo
que se obtiene son tres corrientes de
productos: el destilado D, la corriente lateral
S y el fondo B. Tericamente se puede
alcanzar cualquier pureza deseada en cada
uno de los productos.
Algunos investigadores han llegado a la
conclusin que este sistema ahorra energa
debido a que elimina los retromezclados
que se presentan en las columnas
convencionales (Triantafyllou y Smith, 1990;
Finn, 1993).
Fig. 1. (a) Sistema Petlyuk, (b) Sistema de pared dividida.
3. Desarrollo del modelo

El sistema Petlyuk (Fig. 2a) se
simplifica como se muestra en la Fig. 2b. En
la seccin I, adems de la alimentacin, se
definen los productos netos; el destilado neto
es la diferencia entre el vapor que sale por el
domo y el lquido que entra, y el fondo neto
es la diferencia entre el lquido que sale y el
vapor que entra en la parte baja. Las
funciones que hacen las etapas de
interconexin ahora las hacen el condensador
y rehervidor que se presentan en la seccin I.
Finalmente, la columna principal puede verse
como la suma de dos secciones o columnas
convencionales, tal que la corriente lateral, S,
es la suma de la corriente del fondo de la
seccin II ( B II ) y la corriente de destilado de
la seccin III ( D III ).
Para nuestro estudio tomamos una

alimentacin de tres componentes (F1, F2, F3).
Se utilizan las ecuaciones de Underwood
(1948), Fenske (Henley y Seader, 1990) y
Gilliland (1940) con sus respectivas
consideraciones.
3.1 Diseo del prefraccionador (Seccin I)
En el caso de la seccin I, se deben
especificar inicialmente las recuperaciones de
los componentes clave ligero (f1) y clave
pesado (f3) en el domo. De acuerdo con los
grados de libertad de una columna
convencional, es necesario especificar una
variable adicional; en este caso se usa la
recuperacin del componente intermedio (f2).
87
Fig. 2. Resultados del diseo para el caso de estudio.
De acuerdo con el procedimiento de

Underwood (1948), el valor del flujo del
reflujo mnimo se obtiene sustituyendo las
races obtenidas de (1) en (2) y es el mayor
de las dos soluciones de (2). Fidkowski y
Krlowkoski (1986) encontraron que el valor
ptimo del reflujo se obtiene cuando (2) se
resuelve simultneamente con ambos valores
de . De esta forma se pueden obtener los
valores de LImin y de f2.
3
XF
i i = 1 q,
i =1 i
(1)
1 > 1 > 2 > 2 > 3

3 Fi j
LImin = max
fi
j i =1 i j
88
(2)
Con los resultados del procedimiento de

Underwood (1948) se pueden calcular los
productos netos de la seccin I ( D I y B I ),
las composiciones ( XD I y XB I ), el reflujo
I
mnimo ( Rmin
), y si se considera que los
flujos molares internos son constantes, se
pueden estimar los flujos internos.
Conociendo las composiciones de los
productos netos de la seccin I, se puede
estimar el nmero mnimo de etapas tericas
mediante la ecuacin de Fenske (Henley y
Seader, 1990). Adems, si se ha establecido
una relacin de reflujo para esta seccin, se
puede estimar el nmero de etapas tericas
por medio de la Ecuacin del Gilliland
(1940). Por ltimo, con la correlacin de
Kirkbride (Henley y Seader, 1990) se puede
estimar la localizacin de la alimentacin
para esta seccin. En el Apndice A.1 se
puede encontrar el procedimiento completo
para el diseo de la seccin I.
3.2. Diseo de la columna principal

Las dos secciones de la columna
principal se consideran como la suma de dos
columnas convencionales. La seccin II se
disea separando el componente clave ligero
(f1) del intermedio (f2); el clave pesado se va
completamente al fondo de esta seccin. La
seccin III separa el clave intermedio (f2) del
clave pesado (f3); el clave ligero se va
completamente al domo de esta seccin.
Del diseo de la seccin I, se conocen
los flujos de las alimentaciones a las
secciones II y III pero no se conocen sus
composiciones. Tampoco se conocen las
composiciones de las corrientes laterales que
van hacia la seccin I, y no estn en
equilibrio con las alimentaciones. Por lo
tanto el procedimiento de Underwood (1948)
no se puede utilizar directamente para
determinar el reflujo mnimo de las secciones
II y III.
Glinos y Malone (1985) desarrollaron
un procedimiento para estimar el reflujo
mnimo en sistemas con rectificador lateral.
Estos sistemas presentan un problema similar
al que se tiene para la seccin II. La
conclusin a la que llegaron fue que una
alimentacin a una columna de destilacin
con una corriente lateral a una o unas pocas
etapas por arriba o debajo de dicha
alimentacin puede representarse por una
seudoalimentacin con la condicin trmica
desplazada (3) y (4).
q II =
3
i XD I i
i =1 i
II
L
I
V L
= 1 q II ,
DI
1 > II > 2
(3)
(4)
Si aplicamos el mismo procedimiento

para la seccin III se obtienen las expresiones
equivalentes para esta seccin, ecuaciones
(5) y (6). De esta forma en ambas secciones
se pueden reducir la alimentacin y corriente

lateral a una seudo-alimentacin con flujos y
composiciones conocidas.
q III =
___
I
___
I
___
I
L V
i XBIi
i =1i
III
= 1 q III ,
____
I
BI
2 > III > 3
(5)
(6)
Las Ec. (3) y (4) proporcionan el valor

de la raz de para la seccin II, de la misma
forma las que las Ec. (5) y (6) lo hacen para la
seccin III. Con estas races se puede estimar
el reflujo mnimo necesario para cada
seccin.
Una vez calculados los flujos mnimos
internos, se procede a establecer cual de las
dos secciones es la dominante. Para poder
comparar los flujos internos es necesario que
el flujo de reflujo mnimo de la seccin III
sea llevado hasta el domo de la seccin II,
para comparar directamente el flujo mnimo
de la seccin II con el flujo mnimo necesario
para la seccin III (7). De acuerdo con
Nikolaides y Malone (1987), para sistemas
con mltiples alimentaciones la seccin que
requiera el mayor flujo interno es la
dominante.
LIImin
D
R min = max
III
I II
L min + S D q
(7)
Se pueden en seguida calcular los flujos

internos a la condicion del reflujo de
operacin ( R ):
LII = R*D
LIII = LII + DIq II-S
(8)
89
Finalmente, habiendo calculado los

flujos a condiciones mnimas y a condiciones
de operacin, se pueden terminar de disear
ambas secciones de la columna principal. En
el Apndice A.2 se describe el procedimiento
para disear completamente la columna
principal, lo cual implica la determinacin
del nmero de etapas tericas y las etapas
donde se debe de interconectar la columna
principal con el prefraccionador.
3.3 Ajuste del nmero de etapas tericas de

la columna principal con el prefraccionador
Normalmente el nmero de etapas para
el prefraccionador no coincide con el nmero
de etapas tericas que suman la parte baja de
la seccin II y la parte de arriba de la seccin
III. Esto no representa un problema crtico
para el sistema Petlyuk, en donde el
prefraccionador y la columna principal se
pueden construir independientes. En el caso
de la columna de pared dividida s es
necesario que coincidan, aunque algunos
autores (Triantafyllou y Smith, 1990) han
sugerido que se puede utilizar empaques con
diferentes alturas equivalentes. Para evitar
este problema, se puede hacer el diseo tal
que el nmero de etapas tericas entre en
prefraccionador y las secciones mencionadas
de la columna principal coincidan.
Una forma de hacer coincidir el
nmero de etapas tericas es ajustar las
recuperaciones
especificadas
para
el
prefraccionador. Sin embargo, esta opcin
podra presentar el problema de que en
ocasiones las recuperaciones especificadas
no se alcanzaran.
Otra forma de manipular el nmero de
etapas tericas es ajustar la relacin de
reflujo que se especific para el
prefraccionador, lo cual es mas conveniente,
ya que el flujo neto y las composiciones del
destilado y fondos del prefraccionador
permanecen constantes.
90
4. Resultados y discusin
4.1 Caso de estudio

Para hacer una validacin del modelo
de diseo aqu presentado, se tom un
problema base (Tablas 1 y 2). La relacin de
reflujo ( R ) se fij en 1.2Rmin
y las
recuperaciones
en
el
domo
del
prefraccionador fueron 0.95 y 0.05 para f1 y f3
respectivamente. Se dise la columna
necesaria para llevar a cabo la separacin, y
los resultados se utilizaron como una base
para una simulacin rigurosa.
Tabla 1. Especificacin de la alimentacin para el
caso de estudio.
Compuestos
F
q
i
(kg-mol/hr)
n-Pentano
14.966
5.366
n-Hexano
14.966
2.281
n-Heptano
15.419
1.000
Total
45.351
0
Tabla 2. Especificacin de los productos.
Compuestos
n-Pentano
n-Hexano
n-Heptano
Total
XD
0.98
0.02
0.00
1.00
XS
0.09
0.82
0.09
1.00
XB
0.00
0.04
0.96
1.00
La simulacin rigurosa se realiz

mediante el simulador comercial ASPEN
PLUS 10.1. Las propiedades termodinmicas
se estimaron utilizando la ecuacin de estado
de Redlich-Kwong.
4.2 Resultado
La Fig. 3 presenta el diseo preeliminar
de la columna de destilacin. En este caso se
recurri a modificar la relacin de reflujo del
prefraccionador para hacer coincidir el
nmero de etapas tericas en la zona del
prefraccionador.
Fig. 3. Perfil de la fraccin mol del lquido en la columna principal (diseo original).
La Fig. 4 presenta los resultados

obtenidos por el simulador ASPEN PLUS. Se
observa que, aunque la simulacin rigurosa
no alcanza las composiciones establecidas, la
aproximacin inicial es buena. En el
destilado de la columna Petlyuk casi se logr
la separacin deseada, ya que las
composiciones del destilado coinciden con
las especificadas. En el caso de la corriente
lateral y de fondo, sin embargo, se presenta
una mayor desviacin. Esto se debe
principalmente a que la seccin III es la
seccin dominante, es decir, que necesitaba el

mayor flujo interno.
La Fig. 5 presenta los resultados de la
simulacin rigurosa, en la cual se tuvo que
incrementar el reflujo de la columna principal
desde 3.06 hasta 3.6 para lograr que las tres
composiciones de los productos se
satisfacieran, sin sobrepurificarse. Debe
resaltarse que el nmero de etapas ideales del
prefraccionador y de la columna principal, as
como las etapas de alimentacin e
interconexin entre el prefraccionador y
columna principal, no se modificaron.
91
Fracc. Mol Lquido
Etapas Ideales
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.80
0.90
1.00
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
n-Pentano
n-Hexano
n-Heptano
n-Pentano
n-Hexano
n-Heptano
Fig. 4. Perfil de la fraccin mol del lquido en la columna principal (diseo original).
Fracc. Mol Lquido
Etapas Ideales
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.80
0.90
n-Pentano
n-Hexano
n-Heptano
n-Pentano
n-Hexano
n-Heptano
Fig. 5. Perfil de la fraccin mol del lquido en la columna principal (reflujo a 3.6).
92
1.00
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Conclusiones
En este trabajo se propone un mtodo
de diseo para la secuencia de destilacin
trmicamente acoplada Petlyuk y/o pared
dividida sencillo de aplicar. El mtodo
propuesto se puede utilizar como una buena
base para diseos detallados mediante
simulaciones rigurosas. Adems el mtodo
calcula flujos de interconexin entre el
prefraccionador y columna principal, los
cuales se pueden utilizar como primera
aproximacin para la optimizacin del
consumo de energa del sistema Petlyuk. Para
el caso de estudio analizado, por ejemplo, se
requiri slo de leves modificaciones para
obtener los requerimientos de diseo
establecidos.
El mtodo se utiliz para disear una
columna Petlyuk que separa una mezcla de
tres componentes, pero tambin puede ser
utilizado para el diseo de columnas Petlyuk
que
necesiten
separar
mezclas
multicomponentes en tres productos.
Como parte del desarrollo del
mtodo, se obtuvieron relaciones directas
para estimar el reflujo mnimo y la
recuperacin del componente intermedio en
el prefraccionador, cuando la alimentacin es
como lquido saturado o vapor saturado.
Estas relaciones simplifican en gran medida
el clculo del diseo, considerando que en
muchos casos, la alimentacin a una columna
de destilacin cumple o est cerca de esas
condiciones de saturacin.
Nomenclatura
f
Fraccin de recuperacin de un
componente en el domo de la columna, (0 a
1)
q
Condicin trmica de la alimentacin
B
Flujo del producto en el fondo de la
columna, Kg-mol/hr
D
Flujo del producto en el domo de la
columna, Kg-mol/hr
F
Flujo de alimentacin al sistema, Kgmol/hr
L
Flujo lquido en la parte superior a la
alimentacin de la columna
__
L
Flujo lquido en la parte inferior a la
alimentacin de la columna, Kg-mol/hr
N
Nmero de etapas o platos ideales
P
Presin en el domo de la columna, atm
R
Relacin de reflujo
V
Flujo de vapor en la parte superior a la
__
V
Flujo de vapor en la parte inferior a la
XB
Fraccin mol del destilado
XD
Fraccin mol del destilado
XF
Fraccin mol de la alimentacin al
sistema
smbolos griegos
Volatilidades relativas
Races de la Ecuacin de Underwood

Subndices
i
Referencia al nmero de componente
en la mezcla
cl
Componente clave ligero
cp
Componente clave pesado
min Mnimo
Superndices
I, II, III
Seccin del sistema Petlyuk
(ver Fig. 2)
Agradecimientos
Este trabajo fue apoyado por una Beca
de CONACYT a uno de los autores. Adems
la investigacin fue apoyada por el SIHGO,
bajo el No. del proyecto 19980205002.
Apndice A
Descripcin del mtodo de diseo
A.1 Diseo del prefraccionador (Seccin I)

Para disear el prefraccionador o
seccin I, es necesario especificar f1 y f3 en el
domo de esta seccin, as como especificar la
relacin de reflujo ( R I ).
La Ecuacin de Underwood (1948) se
pueden escribir de la siguiente forma:
93
c
XF
i i = 1 q,
i =1 i
LImin
LImin
XD I XB I
cp
cl
log
p
p
XB XD cp
cl
I
N min =

log cl
cp
1 > 1 > 2 > 2 > 3

(9)
c
F
= i 1 fi
i =1 i 1
c
F
= i 2 fi
i =1 i 2
D I = Fi * f i
i =1
F *f
XD iI = i i
DI
B I = Fi * (1 f i )
i =1
R Imin =
F *(1 f i )
XBiI = i
BI
LImin
D
El nmero de etapas ideales N I para el

prefraccioandor se estima por medio de la
correlacin de Gilliland (1940).
(11)
x=
(12)
(13)
Para calcular el flujo mnimo en esta

seccin se deben estimar las dos races de (9)
y resolver el sistema de ecuaciones (10). De
esta forma se obtienen el flujo mnimo y la
recuperacin del componente intermedio.
Con las relaciones (11-13) se calculan los
productos netos del prefraccionador y sus
composiciones.
Dos casos especiales se presentan
cuando la condicin de alimentacin es como
lquido saturado o vapor saturado. En estos
casos la solucin simultnea del sistema de
ecuaciones
(10)
presenta
soluciones
analticas. En el Apndice B.1 se presenta el
desarrollo
de
las
ecuaciones
para
alimentacin como lquido saturado y en el
Apndice B.2 se presenta cuando la
alimentacin es vapor saturado.
El nmero mnimo de etapas ideales
para el prefraccionador se calcula con la
ecuacin de Fenske (Henley y Seader, 1990):
94
(14)
(10)
R I R Imin
RI +1
Donde:
c
y=
N I N Imin
NI + 1
(15)
1 + 54.4x x 1
y = 1 exp
11 + 117.2 x x 0.5
La etapa de alimentacin se puede

estimar por medio de la ecuacin de
Kirkbride:
I
N sup
log
NI
inf
I
I
= 0.206 log B XFcp XB cl

I XF
I
cl XD cl
D
(16)
A.2 Diseo de la seccin II y III

A.2.1 Estimacin del reflujo mnimo y
balance de masa entre las dos secciones
Por balance de masa en cada una de las
secciones se obtienen el flujo y
composiciones de las corrientes auxiliares
B II y D III . Se considera que el componente
clave pesado que contiene D I se va
completamente a B II , y que el componente
ligero que sale por B I se va completamente a
D III .
El destilado del prefraccionador, que es
la alimentacin a la seccin II, va a tener una
condicin trmica de acuerdo con la Ec. (3).
Asimismo, para la seccin III la

condicin trmica de la alimentacin ser de
acuerdo con (5). Con estas condiciones
trmicas para cada seccin, se calculan las
races de con las Ecs. (4) y (6). Finalmente
el reflujo mnimo para cada seccin se
obtiene mediante:
3
R II
min =
i XDi
i =1i
II
1,
(17)
II
LIII
min = R min D
3
i XDIIIi
R III
=
1,
min
III
i =1 i
(18)
Ahora se compara cual de las dos

secciones es la dominante mediante (7). Se
obtienen luego los flujos a condiciones de
operacin con la ecuacin (8). Finalmente el
reflujo al que van a operar cada una de las
secciones por separado es el siguiente:
R III =
LII
D
LIII
D
III
Apndice B
Flujos mnimos para la seccin I cuando
q=0 y q=1
III
III
LIII
min = R min D
R II =
correspondientes corrientes de destilados y

fondos. As tambin, el nmero de etapas
ideales se obtiene aplicando la correlacin de
Gilliland (15) en cada seccin. Finalmente,
las etapas donde se van a conectar las
secciones II y III con la seccin I se estiman
utilizando la ecuacin de Kirkbride (16),
considerado a D I y a B I como las
alimentaciones.
(19)
(20)
A.2.2 Diseo final de la seccin II y seccin

III
El componente calve intermedio es
ahora el componente clave pesado para la
seccin II y el clave ligero para la seccin III.
El destilado de la seccin II es el destilado
final del sistema D, y el fondo de esta
seccin es la corriente auxiliar B II ; en la
seccin III, el destilado es la corriente D III y
el fondo es el fondo final del sistema B. El
nmero mnimo de etapas ideales de ambas
secciones se obtiene por medio de la
ecuacin de Fenske (14), utilizando las
B.1 Alimentacin como lquido saturado

(q=0)
Cuando la alimentacin a la columna se
tiene como lquido saturado, el clculo de
en la ecuacin (9) se convierte en una
solucin de una ecuacin algebraica de
segundo orden. Por lo tanto se pueden
encontrar las races analticamente.
1XF1 2 XF2 3 XF3
+
+
=0
1 2 3
(21)
donde:
1 =
b + b 2 4ac
2a
2 =
b b 2 4ac
2a
(22)
a = F11 + F22 + F33
F11(2 + 3) + F22(1 + 3)
b =
+ F313(1 + 2)
c = (F1 + F2 + F3)123
Como consecuencia, el flujo mnimo de

reflujo y la recuperacin del componente
intermedio en el prefraccionador se obtienen
en este caso mediante las siguientes
expresiones.
95
f2 =
2 3
2
f1 1
f3
1 3
1 3
LImin =
F 3
F1
f1 +
f3
1
1 3
Referencias
(23)
(24)
B.2 Alimentacin como vapor saturado

(q=1)
Cuando la alimentacin a la columna es
en forma de vapor saturado, la ecuacin (9)
se convierte en una ecuacin cbica. Una de
las races de esta ecuacin es igual a cero. Al
factorizar dicha raz se identifica una
ecuacin cuadrtica como en el caso de
lquido saturado.
1XF1 2 XF2 3 XF3
+
+
=1
1 2 3
(25)
donde:
1 =
2 =
b + b 2 4ac
2a
b b 2 4ac
y 3 = 0
2a
(26)
a = F1 + F2 + F3
F1( 2 + 3 ) + F2 (1 + 3 )
b =
+ F31(1 + 2 )
c = F1 23 + F213 + F31 2
Para este caso, se obtienen las siguientes

expresiones analticas.
( 3 )
( 2 )
f2 = 1 2
f1 + 3 1
f3
2 (1 3 )
2 (1 3 )
(27)
F + F2 1 3 + F31 2
(f1 f 3 ) (28)
LImin = 1 2 3
2 (1 3 )
96
Fidkowski, Z. y Krlikowski, L. (1986).

Thermally coupled system of distillation
columns: optimization procedure. AIChE
Journl 32, 537-546.
Fidkowski, Z. y Krlikowski, L. (1987).
Minimum requirements of thermally
coupled distillation systems. AIChE
Journal 33, 643-653.
Fidkowski, Z. y Krlikowski, L. (1990). Energy
requirements
of
nonconventional
distillation systems. AIChE Journal 36,
1275-1277.
Finn, A.J. (1993). Consider thermally coupled
distillation,
Chemical
Engineering
Progress, October, 41-45.
Gilliland, E. R. (1940). Multicomponent
Rectification;
Industrial
Engineering
Chemistry 32, 1220-1223.
Glinos, K. y Malone, M.F. (1985) Minimum
vapor flows in a distillation column with a
sidestream stripper. Industrial Engineering
Chemical Process Design and Development
24, 1087-1090.
Henley, E.J. y Seader, J.D. (1990), Operaciones
de Separacin por Etapas de Equilibrio en
Ingeniera Qumica. Revert, S. A.
Nikolaides, I.P
y Malone, M.F. (1987).
Approximate design for multiple-feed/sidestream distillation systems, Industrial
Engineering Chemical Research 26, 18391845.
Tedder, D.W. y Rudd, D.F. (1978). Parametric
studies in industrial distillation. AIChE
Journal 24, 303-334.
Triantafyllou, C. y Smith, R. (1990). The design
and optimization of fully thermally coupled
distillation columns. Transactions of the
Institution of Chemical Engineers, 70.
Underwood, A. J. V. (1948). Fractional
distillation of multicomponent mixtures.
Chemical Engineering Progress 44.
Narvez-Garca, A.; Zavala-Lora, J.C.; Rocha-Uribe, A.; Rubio-Atoche, C.

MTODO CORTO PARA LA DESTILACIN DISCONTINUA MULTICOMPONENTE CONSIDERANDO UNA
POLTICA DE REFLUJO VARIABLE
Revista Mexicana de Ingeniera Qumica, vol. 12, nm. 3, diciembre, 2013, pp. 621-637
Distrito Federal, Mxico

amidiq@xanum.uam.mx
Iztapalapa
Mxico
Cmo citar?
Nmero completo
Pgina de la revista
www.redalyc.org
Revista Mexicana de Ingeniera Qumica

Vol. 12, No.CONTENIDO
3 (2013) 621-637
Volumen 8, nmero 3, 2009 / Volume 8, number 3, 2009
DISCONTINUA MULTICOMPONENTE
METODO
CORTO PARA LA DESTILACION
CONSIDERANDO
UNA
POL
ITICA
DE REFLUJO VARIABLE
213 Derivation and application of the Stefan-Maxwell equations
SHORTCUT METHOD
FORde MULTICOMPONENT
BATCH DISTILLATION
(Desarrollo y aplicacin
las ecuaciones de Stefan-Maxwell)
CONSIDERING
A
VARIABLE
REFLUX
POLICY
Stephen Whitaker
A. Narvaez-Garca1 , J.C. Zavala-Lora2 , A. Rocha-Uribe1 y C. Rubio-Atoche1
Biotecnologa / Biotechnology
Universidad Autonoma de Yucatan, Facultad de Ingeniera Qumica, Campus de Ingenieras y Ciencias Exactas,
245 Modelado de la biodegradacin en biorreactores de lodos de hidrocarburos totales del petrleo
Periferico Norte Kilometro 33.5, Tablaje Catastral 13615, Colonia Chuburna de Hidalgo Inn, C.P. 97203, Merida,
intemperizados
en suelos
sedimentos
Yucat
an,y M
exico. Tel: +52 (999) 9460956.
2
Universidad Autonoma(Biodegradation
del Carmen,
Dependencia
Academica
Ciencias
Qumicas
y Petrolera
(DACQyP),
modeling
of sludge bioreactors
of total de
petroleum
hydrocarbons
weathering
in soil
Facultad de Qumica, Cuerpo
Academico Consolidado de Ingeniera Qumica Aplicada, Campus Principal, Calle
and sediments)
56 No. 4 Esquina Avenida
Colonia
BenitoC.A.
JuaLucho-Constantino,
rez, C.P. 24180,
Cd. del Carmen,
Campeche, Mexico.
S.A.Concordia,
Medina-Moreno,
S. Huerta-Ochoa,
L. Aguilera-Vzquez,
A. JimnezTel: +52 (938) 3828484.
Gonzlez y M. Gutirrez-Rojas
259 Crecimiento, sobrevivencia y adaptacin de Bifidobacterium infantis a condiciones cidas
Recibido 26 de Octubre de 2012; Aceptado 14 de Julio de 2013

(Growth, survival and adaptation of Bifidobacterium infantis to acidic conditions)
Resumen
L. Mayorga-Reyes, P. Bustamante-Camilo, A. Gutirrez-Nava, E. Barranco-Florido y A. AzaolaSe ha propuesto un metodo
corto para destilacion discontinua considerando una poltica de reflujo variable. El
Espinosa
metodo esta basado en265elStatistical
trabajoapproach
de Sundaram
y Evans (1993) y considera un balance global de materia, as
to optimization of ethanol fermentation by Saccharomyces cerevisiae in the
como un balance parcial de materia. La relacion funcional entre las concentraciones del fondo y del domo es
presence of Valfor zeolite NaA
calculada mediante la utilizaci
on de las ecuaciones de Fenske, Underwood y Gilliland (FUG). La solucion al modelo
(Optimizacin
estadstica
fermentacin
etanlica de
Saccharomyces
cerevisiae
en presencia
de
matematico propuesto se ha
comparado
con de
unlam
etodo riguroso
utilizando
varias
mezclas
del tipo
Clase I y Clase
zeolitaes
Valfor
zeolite NaA)
II. El metodo corto propuesto
de simple
aplicacion y puede ser utilizado para trabajos preliminares de diseno,
simulacion y optimizacion.G. Inei-Shizukawa, H. A. Velasco-Bedrn, G. F. Gutirrez-Lpez and H. Hernndez-Snchez
Palabras clave: metodo corto, FUG, destilacion discontinua, poltica de reflujo variable.
Ingeniera de procesos / Process engineering
Abstract
271 Localizacin de una planta industrial: Revisin crtica y adecuacin de los criterios empleados en
A Shortcut method for batch
has proposed considering a variable reflux policy. The method is based on
esta distillation
decisin
the work of Sundaram and(Plant
Evans
(1993)
and considered
an overall
material
balance, as well as a partial material
site selection: Critical
review and adequation
criteria used
in this decision)
balance. The functional relationship
between
the
concentrations
of
the
bottom
and the dome is calculated using
J.R. Medina, R.L. Romero y G.A. Prez
the Fenske, Underwood and Gilliland equations (FUG). The solution to the proposed mathematical model has been
compared with the rigorous method using Class I and Class II mixtures. The proposed shortcut method is very
simple and it can be used for preliminary design, simulation and optimization work.
Keywords: shortcut method, FUG, batch distillation, variable reflux policy.
Autor
para la correspondencia. E-mail: narvaez1410@hotmail.com Tel.: +52 (999) 9460956, Fax +52 (999) 9460956
Publicado por la Academia Mexicana de Investigacion y Docencia en Ingeniera Qumica A.C.
621
Narvaez-Garca y col./ Revista Mexicana de Ingeniera Qumica Vol. 12, No. 3 (2013) 621-637
Introduccion
La destilacion batch o discontinua es un proceso

ampliamente utilizado para la separacion de pequenas
cantidades de mezclas, la recuperacion de pequenas
cantidades de materiales peligrosos en corrientes de
desechos, la recuperacion de solventes, as como
para la obtencion de productos farmaceuticos y
biotecnologicos de alto valor agregado entre otros,
por tanto, el desarrollo de modelos matematicos para
la prediccion del comportamiento del proceso es un
trabajo que ha ganado interes a lo largo de los u ltimos
tiempos.
Una columna de destilacion batch o discontinua
puede operar utilizando cualquiera de las siguientes
polticas: reflujo constante, reflujo variable y
reflujo o ptimo. El comportamiento del proceso
puede predecirse mediante la utilizacion de modelos
matematicos basados en los balances de materia y
energa. Los modelos matematicos que se obtienen
pueden clasificarse como modelos simplificados
(metodos cortos), modelos semirigurosos, modelos
rigurosos y modelos de reduccion de orden.
En la actualidad la obtencion de modelos rigurosos
es un campo de mucho interes sobre todo porque las
computadoras actuales cuentan con un alta precision y
capacidades de procesamiento mayores, sin embargo,
el uso generalizado de equipos tales como las tabletas
y/o computadoras portatiles de menor capacidad de
procesamiento de datos, hacen factible la busqueda
de metodos simplificados que puedan predecir el
comportamiento de los procesos sin tener que recurrir
a mayores capacidades de procesamiento. Aunado a
esto, siempre sera necesario contar con este tipo de
metodos para obtener datos iniciales para el proceso
de optimizacion matematica.
Cabe mencionar que el modelado y la simulacion
del proceso de destilacion discontinua, utilizando una
poltica de reflujo variable, ha sido abordado de
forma exitosa utilizando teora de control (Monroy
Loperena, R. y Alvarez-Ram
rez, J., 2003) y redes
neuronales (Cressy, D.C. y col., 1993). Tambien, se
han realizado trabajos tendientes a encontrar el efecto
de la derivacion del reflujo (Zavala y col., 2006).
Entre los trabajos reportados en la literatura se
considera que el primer modelo de la dinamica
completa de la columna fue resuelto por Distefano
(1968) y es conocido como un modelo o metodo
riguroso. Domenech y Enjalbert (1981) presentaron
reducciones al modelo riguroso considerando entre
otras restricciones que la columna puede trabajar sin
acumulacion de lquido en cada uno de los platos y en
622
el condensador, este tipo de modelo es conocido como

semiriguroso.
Diwekar (1988) presento un metodo simplificado
o corto basado en un balance global de materia, un
balance parcial respecto al componente (i) y, para la
relacion funcional entre la concentracion del fondo y
el domo, las ecuaciones de Hensgestebeck-Geddes,
Underwood y Gilliland. La solucion presentada por
Diwekar considero el uso de las polticas de reflujo
constante y reflujo variable.
Sundaram y Evans (1993) tambien presentaron un
metodo corto basado en un balance global de materia
y un balance parcial respecto al componente (i). En
ese caso el trabajo u nicamente considero una poltica
de reflujo constante. Para la relacion funcional entre
la concentracion del fondo y el domo, utilizaron las
ecuaciones de Fenske, Underwood y Gilliland (FUG).
Otros autores que han presentados trabajos al respecto
en tiempos mas recientes son Barolo y Guarise (1996)
y Ehsani (2002).
A diferencia de los metodos rigurosos que
consideran la dinamica completa de la columna, los
metodos reducidos o cortos son modelos matematicos
que predicen el comportamiento del proceso
considerando la menor cantidad de ecuaciones,
generalmente realizando un balance global de materia
y los balances parciales en el fondo de la columna
considerando un componente cualesquiera i. Aun
cuando los equipos de computo actuales tienen
capacidades de procesamiento y almacenamiento
mayores, los metodos cortos continuan justificandose
a partir de considerar la reduccion en los tiempos
de computo, as como una precision aceptable en
los resultados obtenidos respecto al metodo riguroso,
ademas, son una herramienta adecuada para encontrar
valores iniciales para el proceso de optimizacion
matematica. Tambien, son utilizados para el diseno
de la columna al permitir la obtencion de condiciones
lmite del proceso. Los metodos cortos son muy
sencillos de aplicar y por lo consiguiente son de
utilidad en los procesos de ensenanza-aprendizaje.
Ademas, los metodos cortos pueden proporcionar
condiciones lmites del proceso (numero mnimo de
etapas y relacion de reflujo mnima) y son metodos
que permiten obtener de forma rapida valores iniciales
para los procesos de optimizacion matematica donde
la complejidad de los metodos requieren datos muy
cercanos a la solucion. Por otra parte, el uso
extensivo de equipos de computo pequenos y de menor
capacidad tales como las tabletas, los Smartphones
y las computadoras mini, cuyas capacidades son
limitadas o reducidas, requieren de programas con
www.rmiq.org
menos requerimientos, de ah que el desarrollo de

un metodo corto sea de mucha utilidad en este
sentido. Por tanto, el presente trabajo retoma la
metodologa propuesta por Sundaram y Evans (1993)
y propone un metodo corto considerando una poltica
de reflujo variable. El modelo matematico obtenido
considera volatilidades relativas constantes, flujo
molar constante as como acumulacion despreciable
de vapor y lquido en los platos y el condensador.
La diferencia de esta propuesta con los metodos
cortos reportados en la literatura, es que el presente
trabajo utiliza directamente el sistema de ecuaciones
de Fenske-Underwood-Gilliland, mientras que otras
propuestas consideran la ecuacion de HengestebeckGeddes, la cual, es una derivacion de la ecuacion
de la Fenske. Tambien, esta propuesta no realiza la
comparacion entre las relaciones de reflujo mnimo de
Underwood y de Gilliland, lo que tiende a reducir el
tiempo de calculo.
Al igual que otros metodos reportados, esta
propuesta puede ser utilizada para obtener condiciones
de diseno y/o parametros iniciales para la optimizacion
matematica del proceso de destilacion discontinua,
tambien conocida como solucion de problemas de
control o ptimo. Las limitaciones que impone el uso
de volatilidades relativas constantes a lo largo del
tiempo y de la columna pueden ser superadas al
calcularse las volatilidades relativas del fondo en cada
paso y obteniendo una volatilidad relativa promedio.
Para cada corte, la volatilidad relativa del domo
permanece sin cambios significativos debido a que
la concentracion del producto deseado permanece sin
variaciones o sus variaciones son muy pequenas.
Este trabajo solo considera el uso de volatilidades
relativas constantes, sin embargo, si existieran
cambios significativos en la volatilidad relativa se
puede considerar lo siguiente: para mezclas ideales, el
equilibrio lquido-vapor se calcula utilizando la Ley de
Raoult y se aplican un punto de burbuja para conocer
la temperatura en el fondo y en el domo. Para el
caso en que el coeficiente de actividad () no pueda
considerarse cercano a la unidad se puede utilizar la
Ley de Raoult modificada y calcular dicho coeficiente
utilizando algun modelo de solucion (Wilson, NRTL,
UNIQUAC, UNIFAC, etc.). En los casos en que el
vapor de la mezcla no tenga un comportamiento ideal
sera necesario utilizar una ecuacion de estado para
obtener el coeficiente de fugacidad ().
Para la utilizacion de las ecuaciones de
Underwood se han considerado separaciones Clase
I y tambien separaciones Clase II, definidas por
Shiras y col., (1950) de la siguiente manera: Clase
I son aquellas separaciones que al efectuarse en una

columna, de platos infinitos, todos los componentes de
la alimentacion se encuentran presentes en el producto
del domo y en el producto del fondo, es decir, todos los
componentes se encuentran distribuidos (presentes)
a lo largo de la columna y, Clase II son aquellas
separaciones que al efectuarse en una columna, de
platos infinitos, uno o varios de los componentes de
la alimentacion solo se encuentran presentes en el
producto del domo o en el producto del fondo, es
decir, uno o mas de los componentes no se distribuye
a lo largo de la columna. De igual forma, un concepto
importante en el desarrollo del modelo es el concepto
de componente clave ligero y componente clave
pesado definidos de la siguiente manera: componente
clave ligero (lk) es aquel componente ligero de una
mezcla que se encuentra presente en el residuo en
cantidades importantes y el componente clave pesado
(hk) es aquel componente pesado que se encuentra
presente en cantidades importantes en el destilado.
A continuacion se presentan las ecuaciones que
conforman el metodo corto desarrollado por Sundaram
y Evans:
Balance global:
V
t
(1)
Bnueva = Banterior
1 + R anterior
Balance parcial respecto al componente i
(i)
x(i)
B,nueva =x B,anterior +
(i)
[x(i)
D x B ]anterior
B
nueva
Banterior
Banterior
Ecuacion de Fenske y modificaciones:

nh x(i) ih x(k) io
B
log x(i)D (k)
xD
B
Nmin =
log ai,k
x(k)
(i) D
Nmin
x(i)
ai,k
D = xB
x(k)
B
x(k)
D = Pc
i=1
x(k)
B
Nmin
x(i)
B ai,k
Ecuaciones de Gilliland:
(1 + 54.4X)(X 1)
Y = 1 exp
(11 + 117.2X) X
N Nmin
Y=
N+1
R Rmin
X=
R+1
Ecuacion de Eduljee:
www.rmiq.org
Y = 0.75 X 0.5668
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
623
Ecuacion de Underwood (Clase I)

Rmin =
Nmin
l,k
l,k
Pc
Nmin
(l,k 1) i=1 xBi l,k
(10)
esta ecuacion (14) requiere el calculo de y la

concentracion del domo. Para se utiliza la ecuacion
(11) [Ecuacion de Underwood]:
c
X
i,k x(i)
B
=0
i=1 i,k
Ecuaciones de Underwood (Clase II):

c
X
i,k x(i)
B
=0
i,k
i=1
Rmin + 1 =
c
X
i,k x(i)
D
i,k
i=1
(11)
(12)
"
x(i)
D
Cabe senalar que el trabajo de Sundaram y Evans

(1993) solo considera la separacion de mezclas Clase
I. Las ecuaciones (11) y (12) se presentan en este
trabajo porque la propuesta que se desarrolla hace uso
de ellas en la separacion de mezclas Clase II. Las
ecuaciones (11) y (12) son ecuaciones no lineales que
requieren de algun metodo numerico para resolverlas,
tal y como el metodo de Newton-Raphson.
y, para las concentraciones en el domo se utiliza la

ecuacion (4) [Ecuacion de Fenske] para el componente
(i):
=
x(i)
B
x(k)
D
x(k)
B
#
i = 1, 2, ..., n;
Nmin
i,k
i,k
siendo k el componente de referencia. La composicion

del componente de referencia (k), puede obtenerse
utilizando una modificacion de la ecuacion (5)
de forma que quede expresada en funcion del
componente clave ligero (lk). La modificacion se
obtiene utilizando la restriccion de la suma de las
composiciones en el domo
c
X
Metodo simplificado o corto
x(i)
D =1
(15)
i=1
Hasta la actualidad, no se han reportado en la

literatura del a rea, propuestas tendientes a obtener un
metodo corto para destilacion discontinua utilizando
el trabajo de Sundaram y Evans (1993) considerando
la poltica de reflujo variable, por tanto, el presente
trabajo, retoma el metodo corto mencionado y
propone una metodologa que considera la poltica
de reflujo variable. Para el caso de reflujo variable,
la concentracion del componente deseado permanece
constante durante toda la operacion, por tanto, en cada
tiempo es necesario calcular la relacion de reflujo,
calculo que puede realizarse utilizando un despeje de
la ecuacion (8):
R=
X + Rmin
1X
Rmin
624
c
X
(lk)
(lk)
x(i)
D xD = 1 xD
(16)
i=1
y, sustituyendo la ecuacion (4) en esta ecuacion (16)

se tiene que:
)
c (" (k) #
X
xD
(lk)
Nmin (i)
(i,k )
xB x(lk)
D = 1 xD
(k)
x
i=1
B
(17)
La propia ecuacion (4) para la concentracion del

componente lk queda expresada como:
(13)
"
Esta ecuacion (13) requiere el calculo de Rmin

por lo que se puede utilizar la ecuacion (10) o
la ecuacion (12). La ecuacion (10) es aplicable
cuando se considera que todos los componentes estan
distribuidos desde el fondo hasta el domo. La ecuacion
(12) se aplica cuando alguno de los componentes no
se distribuye a traves de la columna. Para el caso de
utilizar la ecuacion (12) se tiene que:
c
X
i,k x(i)
D
= 1 +
i,k
i=1
expresion a la cual puede restarse la concentracion del

componente ligero (lk):
(14)
x(lk)
D =
x(k)
D
x(k)
B
#
Nmin
x(lk)
B lk,k
(18)
Ahora bien, sustituyendo esta ecuacion (18) en la

ecuacion (17) y despejando para obtener x(k)
D , se tiene:
"
x(k)
D
#(" X
c
x(k)
B
i=1
x(k)
D
www.rmiq.org
#
)
Nmin (i)
Nmin (lk)
i,k
xB lk,k
xB = 1 x(lk)
(19)
D
(k)
(1 x(lk)
D )x B
= hP
i
(lk) Nmin
c
Nmin (i)
i=1 i,k x B x B lk,k
(20)
Entonces, la ecuacion (4) puede quedar expresada

como:
Nmin (i)
[1 x(lk)
D ]i,k x B
=
i
(lk) Nmin
Nmin ( j)
j=1 i,k x B x B lk,k
x(i)
D = h Pc
etapas, por tanto, se puede utilizar la ecuacion

(21) o la ecuacion (22) para obtener la funcion de
Nmin necesaria para el metodo de Newton-Raphson.
Entonces, cuando i = lk = 1, la ecuacion (17) puede
expresarse de la siguiente manera:
(1 x(lk)
D )
n P j,k Nmin h x( j) i Nmin h x(lk) io
c
B
B
lk,k
j=1 i,k
i,k
x(i)
x(i)
B
c
x(1) X
Nmin
D
(1)
xB
(21)
f (Nmin ) =x(1)
D
Si se considera que el componente clave ligero (lk)

es el componente 1, entonces, esta ecuacion (21) se
reduce a:
(1 x(1)
D )
j,k Nmin h x( j) io ;
c
j=2
i = 1, 2, ...n
(23)
Si se utiliza la ecuacion (6) el proceso es un poco mas

complejo y se requiere un proceso iterativo. En ambos
casos es necesario el valor del numero mnimo de
etapas, por tanto, obtener Nmin es de maxima prioridad.
La ecuacion de Fenske [Ecuacion (3)] permite
calcular el numero mnimo de etapas considerando
los componentes clave ligero (lk) y el componente de
referencia (k), entonces:
nh x(lk) ih x(k) io
B
log D(k) (lk)
xD
xB
Nmin =
(24)
log(lk,k )
(hk)
donde x(lk)
a dada
D es una cantidad conocida y xD est
por la ecuacion (20), entonces, la ecuacion (4) se
puede expresar como:
nh x(lk) ih Pcj=2 Nj,kmin x(Bj) io

D
Nmin =
log
=
nh
1x(lk)
D
log lk,k
ih P
x(lk)
D
1x(lk)
D
c
j=2
Nj,kmin
log lk,k
x( j) io
B
x(lk)
B
(25)
Ahora bien, esta ecuacion (25) requiere un proceso

iterativo para la solucion del numero mnimo de
j,k
B
(1)
xB
(26)
1,k

+1
f 0 (Nmin ) =x(1)
D
(27)
c ( j)
X
x
Nmin
ln
j,k
j=2

4 N Nmin 1.7643
X = 1
3 N+1
x(lk)
B
c ( j)
X
x
Nmin
1=0
B
x(i)
B
i,k
1,k
j=2
(22)
Para las ecuaciones (20), (21) y (22) se requiere

la obtencion del numero mnimo de etapas, ademas
de que es necesario considerar que la ecuacion (12)
requiere el valor de X, que se puede obtener utilizando
las ecuaciones (6) o (9). Si se utiliza la correlacion de
Eduljee [ecuacion (9)] el despeje es:
log
j=2

x(Bj) = 1 x(1)
D
i = 1, 2, ...n
x(i)
D = nP
j,k
x(1)
B
1,k

j,k
1,k
(28)
f (Nmin )
(29)
f 0 (Nmin )
Con el valor del numero mnimo de etapas ya se
puede obtener la relacion de reflujo y los demas
valores que involucran a la variable Nmin . El equilibrio
lquido-vapor de las mezclas es obtenido utilizando la
ecuacion (30):
Nmin,nueva = Nmin,anterior
i,k x(i)
( j)
j=i [ j,k x ]
y(i) = Pc
(30)
siendo k el componente de referencia.
Procedimiento de solucion
A continuacion se detalla el procedimiento para

obtener el perfil de las relaciones de reflujo requeridas
para mantener constante la concentracion del
componente clave ligero (lk) en el producto. De inicio
se debe contar con los siguientes datos conocidos:
on (F),
composicion [x(i)
B ] y cantidad de alimentaci
flujo de vapor (V) a traves de la columna, numero de
etapas teoricas (N), volatilidades relativas (i ) de los
componentes de la mezcla, el tiempo de produccion
(tprod ) y el componente de referencia (k). En la
metodologa presentada por Sundaram y Evans (1993)
el componente clave pesado (hk) es el componente de
referencia.
La solucion inicia obteniendo el numero de etapas
teoricas mnima (Nmin ) que es un valor necesario para
calcular las composiciones en el domo [x(i)
D ] utilizando
la ecuacion de Fenske. Despues, es necesario calcular
www.rmiq.org
625
la cantidad de reflujo mnima (Rmin ) cuyo valor es

utilizado para calcular la relacion de reflujo (R).
Conociendo el valor de R, se puede calcular el valor de
la cantidad de destilado (D) instantaneo. Se considera
el incremento del tiempo (t) y mediante los balances
total y parcial se calcula el remanente en los fondos (B)
y su composicion [x(i)
en se puede conocer la
B ]. Tambi
cantidad de destilado acumulado. La Figura 1 muestra
el diagrama de flujo para la solucion del metodo corto
propuesto.
Casos de estudio
El modelo matematico del metodo corto propuesto ha

sido resuelto considerando varias mezclas (binaria y
multicomponentes) cuyas condiciones de entrada al
proceso se muestran en la Tabla 1. En cada uno
de los casos de estudio se considera que la mezcla
es alimentada a la temperatura de ebullicion. La
tolerancia al error es de 1.0 104 , el paso de
integracion es t = 1.0 101 h y el tiempo de
produccion es el necesario para agotar al componente
clave ligero. Se ha considerado una alimentacion de
200 kmol y un flujo de vapor de 110 kmol/h.
Los resultados de la solucion a los problemas
presentados en la Tabla 1 se comparan con un metodo
riguroso despreciando la acumulacion de lquido en
la columna y en el condensador-tanque de reflujo. En
este trabajo se desprecio la acumulacion de vapor.
El comparativo se realiza utilizando el programa
R (1991) para el modelo riguroso.
MultiBatchDS
4.1
Discusion de resultados
En cada caso, los perfiles de concentracion en el

fondo de la columna (rehervidor) no alcanzaron una
desviacion mayor al 7% para las concentraciones y un
7% para la relacion de reflujo permitiendo considerar
que los valores obtenidos en la solucion del metodo
propuesto son apropiados. La comparacion de los
resultados fue realizada de la siguiente manera: los
datos que se utilizaron para la solucion del modelo
matematico propuesto tambien fueron utilizados en
R
el software MultiBatchDS
obteniendose los perfiles
de composicion y relacion de reflujo para un tiempo
amplio y adecuado para agotar al componente clave
ligero. En cada caso, se tomaron mezclas consideradas
en la literatura con la intencion de validar los datos de
la solucion.
Por ejemplo, el caso 5 fue tomado de Diwekar
y Madhavan (1991). Ellos reportaron una desviacion
626
para los valores de la concentracion del componente

clave (componente 1) del 8.1% maximo. En este
trabajo la desviacion maxima obtenida fue del 6.35%
considerando hasta el 51.8% de destilado, ya que a
partir de ese punto las cantidades de reflujo que se
requieren son muy grandes y la obtencion de producto
es mnima, lo cual implica gastos energeticos mayores
que pueden influir en la economa del proceso, sobre
todo porque gran parte de la energa se retira por el
condensador. La Figura 7 muestra un comparativo del
comportamiento de las concentraciones en el fondo
de la columna, as como del comportamiento de las
relaciones de reflujo, observandose que el metodo
propuesto tiene un mejor ajuste con el modelo riguroso
utilizado.
Los perfiles de las relaciones de reflujo muestran
discrepancia cuando el componente clave ligero (lk)
se encuentra en cantidades pequenas y, se considera
que esto es debido a que existen procesos iterativos
en el metodo corto que pueden introducir pequenos
errores de redondeo que se pueden incrementar en
cada iteracion. La influencia del proceso iterativo en
el calculo de la relacion de reflujo R para el metodo
corto se observa en la ecuacion (13), ya que el valor de
R depende de los valores de Rmin y X. Para el metodo
riguroso el valor de R depende de las concentraciones
del condensador y los dos u ltimos platos de la columna
(Diwekar, 2012).
Tambien, un factor de importancia es el numero
de etapas utilizadas, ya que en los casos resueltos, la
mayor discrepancia se encuentra cuando se utiliza un
numero mayor de etapas (15 y 30 platos) debido a que
la solucion al sistema de ecuaciones de Underwood y
Gilliland se ve afectada con la precision del proceso
iterativo para resolverlo. Otro factor que influye en
los resultados es la desviacion de la idealidad de la
mezcla, y para ese caso, sera necesario considerar
la influencia de la temperatura y la concentracion
para predecir el comportamiento de las volatilidades
relativas, factor que no ha sido considerado en este
trabajo.
Una de las limitantes que se ha observado para
la adecuada solucion del metodo es la eleccion
de los componentes claves, ya que si se hace de
forma inadecuada puede causar una sobrestimacion o
una subestimacion en los balances como mencionan
Zamar y col., (1998). Una forma sencilla de elegir
dichos componentes es predecir la distribucion de
los mismos en la columna de destilacion discontinua
utilizando la ecuacion (31) obtenida por los autores
antes mencionados.
www.rmiq.org
Datos iniciales:
Alimentacin, Composicin de los Fondos, Flujo de Vapor, Nmero de platos,
Tiempo de Produccin, Volatilidades Relativas
Suponer N
min
Calcular N
min
utilizando el mtodo de Newton-Raphson

Ecs. (26)-(28)
Calcular X utilizando las correlaciones de Gilliland o Eduljee

Ec. (14)
Calcular las composiciones del domo utilizando la ecuacin de Fenske Ecs. (4) y
(21)
Calcular utilizando las ecuaciones de Underwood

Ecs. (10) y (13)
Calcular R utilizando la correlacin de Gilliland o la de Eduljee

Ec. (12)
Calcular B y la composicin del fondo utilizando el balance global y el balance

parcial. Ecs. (1) y (2)
Incrementar el tiempo
No
Tiempo tprod
Terminar
530
531
532
Si
Figura 1. Diagrama de flujo para la solucin del mtodo corto considerando una poltica de
reflujo variable.
533
Manuscrito
sometido
lam
Revista
Mexicana
de Ingeniera
Qumica
Fig. 1. Diagrama de
flujo para la
solucion adel
etodo corto
considerando
una poltica
de reflujo variable.19
www.rmiq.org
627
Tabla 1. Condiciones de entrada para los casos de estudio a reflujo variable.

Caso
1
2
3
4
5
Fraccion molar de alimentacion

1
2
3
4
0.40
0.50
0.33
0.41
0.50
0.20
0.30
0.33
0.39
0.50
0.30
0.10
0.34
0.20
-
0.10
0.10
Nmin
(1 i ) = (1 r )ir
Volatilidades relativas
1
2
3
4
1.67
1.67
1.33
1.76
2.40
1.25
1.25
1.00
1.00
1.00
(31)
siendo i el numero de moles del componente i en el

fondo (i = di / fi ) y r el componente de referencia.
Los datos calculados con esta ecuacion (31) se
trazan y ajustan a una lnea recta. De los componentes
que se encuentran sobre la lnea recta se obtienen los
componentes clave, el componente mas ligero sera el
componente clave ligero (lk) y el componente mas
pesado sera el componente clave pesado (hk). Una
descripcion detallada de este proceso se encuentra
en el trabajo de Zamar y col., (1998). La Figura
(8) muestra el comportamiento de la concentracion
y la relacion de reflujo considerando dos diferentes
componentes claves pesados (hk). En este caso, se
obtiene una subestimacion de la concentracion del
componente clave ligero (lk) de hasta un 15.5% y de
un 24.3% en la relacion de reflujo.
Aunque este trabajo es una derivacion del trabajo
de Sundaram y Evans (1993), aqu se considera que el
componente clave ligero (lk) es el componente ligero
que se encuentra en mayor cantidad en los fondos
y el componente clave pesado (hk) el componente
pesado que se encuentra en mayor cantidad en el
domo de acuerdo con la definicion dada anteriormente.
Otra limitante para una adecuada solucion son los
procesos iterativos involucrados en las ecuaciones
de Fenske, Underwood y la correlacion de Gilliland
ya que requieren de un adecuado estimado inicial
incrementando los tiempos de calculo y/o procesos de
ciclados si no se alcanzan valores dentro del lmite
de error deseado y, es que, de acuerdo con Barolo y
Guarise (1996) el sistema formado por las ecuaciones
de Underwood y Gilliland puede no tener solucion,
sin embargo, cuando el numero de etapas es pequena
generalmente si la tiene.
Como el perfil de reflujo es una funcion continua
(una funcion es continua cuando para puntos cercanos
en el dominio se producen pequenas variaciones en
la funcion) la dificultad de solucion generalmente se
produce en el primer paso, siendo los demas pasos
mas sencillos de resolver porque ya se conoce el punto
628
1.00
1.00
0.67
0.68
-
0.83
0.83
-
N+
x(lk)
D
Clase
5
30
10
10
9
3
3
2
2
2
0.70
0.95
0.80
0.99
0.95
I
II
I
I
I
anterior, sin embargo, las propias restricciones en el

error aceptado pueden derivar en procesos ciclados o
desbordes en la memoria de la computadora.
Otra dificultad que puede existir en la solucion
de la ecuacion de Underwood es que en su forma
original es una funcion discontinua cuando i ,
sin embargo, esta dificultad puede ser superada
considerando que la ecuacion de Underwood (en su
forma original) puede ser representada adecuadamente
por un polinomio de orden N=Numero de
componentes (Monroy-Loperena y Vargas-Villamil,
2001) cuyas races se obtienen utilizando cualquier
metodo numerico convencional. Para el caso de la
destilacion discontinua [ecuacion (11)] el polinomio
es de orden N 1 y esta representado por la ecuacion
(32)
c
c
Y
X

(i)
i,k x
j,k = 0
(32)
B
j=1
i=1
j,i
donde k es el componente de referencia y c es el

numero de componentes de la mezcla.
Conclusiones
En este trabajo se ha desarrollado un metodo corto
considerando una poltica de reflujo variable para
el proceso de destilacion discontinua. El modelo
matematico propuesto ha sido desarrollado a partir
del trabajo de Evans y Sundaram (1993) y pretende
ser una alternativa viable dentro del contexto de los
metodos cortos utilizados para el diseno, simulacion,
optimizacion y problemas de control del proceso de
destilacion discontinua. Cabe decir, que a diferencia
del metodo original de Evans y Sundaram para
reflujo constante (donde solo se presentan mezclas con
distribucion Clase I), este trabajo tambien considera
la separacion de mezclas con distribucion Clase II.
Ademas, se consideran las definiciones de componente
clave ligero (lk) y clave pesado (hk) dadas por Shiras
y col., (1950) a diferencia del trabajo de Sundaram y
www.rmiq.org
534
534
535
536
(a) Perfil de concentraciones en el fondo.
535
536
537
538
537
539
538
540
(b) Perfil de relaciones de reflujo.

R

Figura
2. Caso
1: Mezcla
utilizando
MultibatchDS
. cuaternaria, clase I. Comparacin del mtodo corto propuesto
(b) Perfil
de relaciones de reflujo.
con un modelo riguroso utilizando
MultibatchDS.
539
541
540
Figura 2. Caso 1: Mezcla cuaternaria, clase I. Comparacin del mtodo corto propuesto
con un modelo riguroso utilizando MultibatchDS.
541
Manuscrito sometido a la Revista Mexicana de Ingeniera Qumica
Fig. 2. Caso 1: Mezcla cuaternaria, clase I. Comparacion del metodo corto propuesto con un modelo riguroso
www.rmiq.org
629
20
542
543
544
542
543
544
545
546
547
548
545

Fig. 3. Caso 2: Mezcla cuaternaria, clase II. Comparacion del metodo corto propuesto con un modelo riguroso
R

Figura
3. Caso
2: Mezcla
utilizando
MultibatchDS
. cuaternaria, clase II. Comparacin del mtodo corto propuesto
con un modelo riguroso utilizando MultibatchDS.
630
www.rmiq.org
549
546
547

21
Figura 3. Caso 2: Mezcla cuaternaria, clase II. Comparacin del mtodo corto propuesto
550
551
552
550
551
552
553
554
555
556
553
Figura
4. Caso R 3:
MultibatchDS
. Mezcla ternaria, clase I. Comparacin del mtodo corto propuesto con un
modelo riguroso utilizando MultibatchDS.
557
554
555
556
Manuscrito
sometido
a la
Revista
Mexicana dedel
Ingeniera
Qumica
22
Figura 4. Caso
3: Mezcla
ternaria,
clase
I. Comparacin
mtodo corto
propuesto con un
Fig. 4. Caso 3: Mezcla ternaria, clase I. Comparacion del metodo corto propuesto con un modelo riguroso utilizando
www.rmiq.org
631
558
559
560
558
559
560
561
562
563
564
561
565
562

Figura 5. Caso 4: Mezcla ternaria, clase II. Comparacin del mtodo corto propuesto con
unFig.
modelo
riguroso utilizando MultibatchDS.
5. Caso 4: Mezcla ternaria, clase II. Comparacion del metodo corto propuesto con un modelo riguroso
R
utilizando MultibatchDS
.

632
563
564
565
www.rmiq.org
Manuscrito
sometido
a laclase
Revista
Mexicana de del
Ingeniera
Figura 5. Caso
4: Mezcla
ternaria,
II. Comparacin
mtodo Qumica
corto propuesto con23
un modelo riguroso utilizando MultibatchDS.
566
567
568
566
567
568
569
570
571
572
569
Figura 6. Caso 5: Mezcla binaria, Clase I. Comparacin del mtodo corto propuesto con un
573
570
571
572
573
Fig. 6. Caso 5: Mezcla binaria, Clase I. Comparacion del metodo corto propuesto con un modelo riguroso utilizando
R
MultibatchDS
.
www.rmiq.org
633
Manuscrito
sometido
a la
Revista
Mexicana dedel
Ingeniera
Qumica
24
Figura 6. Caso
5: Mezcla
binaria,
Clase
I. Comparacin
mtodo corto
propuesto con un
574
575
576
574
575
576
577
578
579
580
577
581
578
mtodo
un modelo
riguroso
utilizando
MultibatchDS.
Fig. 7. corto
Caso 5:yMezcla
binaria, Clase
I. Comparaci
on del m
etodo corto propuesto con un metodo corto y un modelo
579
580
R
riguroso utilizando MultibatchDS
.
(b) Perfil de
relaciones de reflujo.
www.rmiq.org
25
mtodo corto y un modelo riguroso utilizando MultibatchDS.
634
582
583
584
585
582
583
584
585
586
587
588
589
586
590
587
588
589
590
Fig. 8. Caso 1: Mezcla cuaternaria, Clase I. Solucion del metodo corto propuesto eligiendo como componente
(b)
Perfil(4)dey relaciones
reflujo.
clave pesado (hk) al componente m
as pesado
al componentede
pesado
(3) que se encuentra en mayor cantidad
inicialmente.
Figura 8. Caso 1: Mezcla cuaternaria, Clase I. Solucin del mtodo corto propuesto
eligiendo como componente clave pesado (hk) al componente ms pesado (4) y al
componente pesado (3) que se encuentra en mayor cantidad inicialmente.
www.rmiq.org
635
26
Figura 8. Caso 1: Mezcla cuaternaria, Clase I. Solucin del mtodo corto propuesto
eligiendo como componente clave pesado (hk) al componente ms pesado (4) y al
componente pesado (3) que se encuentra en mayor cantidad inicialmente.
Evans que consideran como componente clave ligero

(lk) al componente mas volatil y como componente
clave pesado (hk) al componente mas pesado.
La eleccion de los componentes claves es
parte fundamental en la solucion del modelo
matematico propuesto ya que puede obtenerse una
sobreestimacion o subestimacion de los balances de
materia. En este trabajo, se presento el impacto de
no elegir correctamente a los componentes claves.
La solucion del modelo matematico se comparo con
los resultados obtenidos mediante un metodo riguroso
considerando una acumulacion de lquido y vapor
igual a CERO y, que se encuentra implementado en
R
el simulador MultiBatchDS
. En todos los casos,
la diferencia de resultados respecto a los perfiles
de concentracion, tanto en el fondo como en el
domo, no excedio al 7%, lo mismo se observa
para los perfiles de las relaciones de reflujo, por lo
que se considera que los resultados obtenidos son
satisfactorios. En cuanto a los perfiles de la relacion
de reflujo, el comportamiento permite observar que se
puede terminar el proceso de separacion cuando una
pequena variacion en los fondos de la concentracion
del componente deseado requiere de una relacion de
reflujo muy grande.
Agradecimientos
Agradecemos a la Universidad Autonoma del Carmen
y a la Universidad Autonoma de Yucatan, en Mexico,
por todas las facilidades brindadas para el desarrollo
de este trabajo, de la misma manera agradecemos al
PROMEP por aportar lo recursos financieros mediante
el convenio PROMEP/103.5/10/5126.
Nomenclatura
B
D
hk
N
c
Nmin
lk
k
R
Rmin
t
V
x
Indices
636
fondos (kmol)
producto (kmol)
componente clave pesado
numero de etapas
numero de componentes
numero mnimo de etapas
componente clave ligero
componente de referencia
relacion de reflujo
relacion de reflujo mnima
tiempo (h)
flujo de vapor (kmol/h)
composicion del lquido
B
fondos
D producto
hk componente clave pesado
lk componente clave ligero
i
componente i
Smbolos griegos
Volatilidad relativa
Parametro de la ecuacion de Underwood
Referencias
Barolo, M. y Guarise, G. B. (1996).
Batch
distillation of multicomponent systems with
constant relative volatilities. Transactions of the
Institution of Chemical Engineers 74, 863-871.
Cressy, D.C., Nabney, Ian T. y Simper, A. (1993).
Neural control of a batch distillation. Neural
Computing and Applications 1, 115-123.
Distefano, G.P. (1968). Mathematical modeling and
numerical integration of multicomponent batch
distillation equations. AIChE Journal 14, 190199.
Diwekar, U.M. (1988).
Simulation, Design,
Optimization and Optimal Control of
Multicomponent Batch Distillation Columns.
Ph. D. Thesis. Indian Institute of Technology,
Bombay, India.
Diwekar, U.M. y Madhavan, K.P. (1991). BatchDist: A comprehensive package for simulation,
design, optimization and optimal control
of multicomponent,
multifraction batch
distillation columns. Computers and Chemical
Engineering 15, 833-842.
Diwekar, U. (2012). Batch Distillation: Simulation,
Optimal Design, and Control. Second Edition.
CRC Press. U.S.A.
Domenech, S. y Enjalbert, M. (1981). Program
for simulating batch rectification as a
unit operation.
Computers and Chemical
Engineering 5, 181-184.
Ehsani, Mohammad R. (2002).
Prediction of
the operating conditions in a batch distillation
column using a shortcut method. Iran Journal
of Chemistry & Chemical Engineering 21, 118125.
Monroy-Loperena, R. y Alvarez-Ram
rez, J. (2003).
A note on the identification and control of batch
www.rmiq.org
distillation columns. Chemical Engineering

Science 58, 4729-4737.
operations. Industrial & Engineering Chemistry

Research 32, 511-517.
Monroy-Loperena, R. y Vargas-Villamil, F.D.

(2001).
Determination of the polynomial
defining Underwoods equations in short-cut
distillation design. Industrial and Engineering
Chemistry Research 40, 5810-5814.
Zamar, Silvina D., Salomone, Enrique y Iribarren,

Oscar A. (1998).
Shortcut method for
multiple task batch distillations. Industrial
& Engineering Chemistry Research 37, 48014807.
Shiras, R.N., Hanson, D.N. y Gibson, C.H. (1950).

Calculation of minimum reflux in distillation
columns. Industrial & Engineering Chemistry
42, 871-876.
Sundaram, S. y Evans, L.B. (1993). Shortcut
procedure for simulating batch distillation
Zavala-Lora, J.C., Cordova-Quiroz, A.V., RoblesHeredia, J.C., Anguebes-Franseschi, F.,

Narvaez-Garca, A. y Elvira-Antonio, N.
(2006). Efecto de la derivacion del reflujo en
destilacion discontinua. Revista Mexicana de
Ingeniera Qumica 5, 109-113.
www.rmiq.org
637
Revista Mexicana de Ingeniera Qumica

ISSN: 1665-2738
amidiq@xanum.uam.mx
Iztapalapa
Mxico
Barroso-Muoz, F. O.; Hernndez -Castro, S.; Segovia-Hernndez, J. G.; Aguilera-Alvarado, A. F.

Estudio de sistemas de destilacin trmicamente acoplados con reaccin qumica
Revista Mexicana de Ingeniera Qumica, vol. 5, nm. Su1, noviembre, 2006, pp. 137-142
Distrito Federal, Mxico
Cmo citar el artculo

Nmero completo
Pgina de la revista en redalyc.org

REVISTA MEXICANA DE INGENIERA QUMICA Vol. 5, Supl. 1 (2006) 137-142
AMIDIQ
ESTUDIO DE SISTEMAS DE DESTILACIN TRMICAMENTE ACOPLADOS

CON REACCIN QUMICA
STUDY OF THERMALLY COUPLED DISTILLATION SYSTEMS
WITH CHEMICAL REACTION
F. O. Barroso-Muoz, S. Hernndez -Castro*, J. G. Segovia-Hernndez, y A. F. Aguilera-Alvarado
Facultad de Qumica, Universidad de Guanajuato, Noria alta s/n C.P. 36050 Guanajuato, Gto, Mxico.
Recibido 1 de Julio 2005; Aceptado 4 de Diciembre 2005
Resumen
Las secuencias de destilacin con acoplamiento trmico han resultado ser una buena opcin para reducir tanto los
consumos de energa como los costos de capital, para la separacin de algunas mezclas multicomponentes.
Recientemente se ha publicado informacin referente al diseo, operacin y control de este tipo de esquemas, que
ha permitido su implementacin prctica en la industria de procesos qumicos. En este trabajo se estudia el proceso
de produccin de acetato de etilo, usando secuencias de destilacin con acoplamientos trmicos a columnas
auxiliares (rectificador o agotador) para llevar a cabo la reaccin-separacin. Se encuentra que la secuencia con
acoplamiento trmico a un agotador lateral es la que presenta ms ventajas, ya que reduce el consumo de energa
del proceso tradicional y adems mediante el agotador remueve el agua que se forma en la reaccin, con lo cual se
favorece la produccin de acetato de etilo y se evita la posterior formacin de azetropos. Este mismo resultado ha
sido obtenido en la columna tipo Petlyuk.
Palabras clave: sistemas trmicamente acoplados, reaccin qumica, acetato de etilo.
Abstract
Thermally coupled distillation systems have been excellent options to save energy and capital investment for the
separation of multicomponent mixtures. Latest information about design, operation and control of these systems
has been published and it has allowed a practical implementation in the chemical industry. In this paper the ethyl
acetate process is studied, using distillation systems with side columns (either stripper or rectifier) to carry out the
reaction-separation. Thermally coupled sequence with side stripper shows more advantages because it reduces
traditional process energy demand and also eliminates the water produced by the chemical reaction which allows
higher ethyl acetate production and avoid azeotrope formation. A similar result has been obtained for the case of
the Petlyuk distillation column.
Keywords: thermally coupled systems, chemical reaction, ethyl acetate.
1. Introduccin
A pesar del alto consumo de energa y baja
eficiencia termodinmica (Tedder y Rudd, 1978;
Flores y col., 2003), la destilacin contina siendo la
tcnica ms utilizada para separar mezclas
multicomponentes, de hecho cada nuevo proceso de
separacin debe ser comparado con sta como una
forma estndar de separacin. Una tendencia actual
en los procesos de separacin es la bsqueda de
nuevos esquemas de destilacin, que puedan reducir
sustancialmente el costo anual de operacin. Una
forma de lograr el objetivo anterior es el empleo de
sistemas de destilacin trmicamente acoplados, los
cuales pueden tener una reduccin en el consumo de
energa de hasta 50% en comparacin con los trenes
de destilacin convencionales de amplio uso en la
industria de procesos qumicos (Hernndez y
Jimnez, 1999a; Triantafyllou y Smith, 1992; Rong y

col., 2003). Este ahorro de energa se logra debido a
que los esquemas con acoplamiento trmico no
presentan remezclado en comparacin con los
esquemas de destilacin convencionales. En los
trenes de columnas de destilacin convencionales,
por ejemplo, la secuencia directa de dos columnas
para la separacin de una mezcla ternaria, el perfil de
composicin del componente intermedio en la
primera columna aumenta desde el domo hasta un
mximo por debajo de la alimentacin y luego
disminuye. Esta disminucin de composicin se
conoce como remezclado y es causa de una baja
eficiencia termodinmica, debido a que se deber de
suministrar ms energa en la segunda columna para
lograr las purezas deseadas. Este problema se evita
en las columnas con acoplamiento trmico, pues en
las secuencias con columnas laterales se extrae una
*Autor para la correspondencia: E-mail: hernasa@quijote.ugto.mx

Tel: 473-7320006 x 8142, Fax: 473-7320006 x 8108
Publicado por la Academia Mexicana de Investigacin y Docencia en Ingeniera Qumica A.C.
137
F. O. Barroso-Muoz y col./ Revista Mexicana de Ingeniera Qumica Vol. 5, Supl. 1 (2006) 137-142
corriente lateral del plato en donde se alcanza el

mximo de la concentracin del componente
intermedio para suministrar la alimentacin a la
segunda columna (Triantafyllou y Smith, 1992;
Hernndez y col., 2003). En general, en un sistema
con acoplamiento trmico, un par de corrientes en
fase lquida y vapor sustituyen ya sea a un
condensador o rehervidor, por lo que la transferencia
de calor se logra mediante contacto directo de
corrientes de materia.
El desarrollo de mtodos de diseo y
optimizacin
de
esquemas
de
destilacin
trmicamente
acoplados
ha
permitido
su
implementacin prctica, reportndose ahorros en el
consumo de energa y capital de 30%. BASF se ha
convertido en la industria lder en la construccin de
columnas de destilacin trmicamente acopladas,
mediante implementacin de paredes divisorias en
las columnas, creando con esto una columna que es
termodinmicamente equivalente a la columna
Petlyuk (Kaibel y Schoenmarkers, 2002).
La implementacin real de sistemas de
destilacin trmicamente acoplados ha sido posible
debido a que se conocen suficientes aspectos sobre el
control y la operacin. Se pensaba que la
implementacin prctica poda presentar problemas
debido a la naturaleza integrada de las columnas.
Estudios reportados (Wolff y Skogestad, 1995;
Hernndez y Jimnez, 1999b; Segovia-Hernndez y
col., 2004) en el rea de control han puesto de
manifiesto que las propiedades tericas de control de
los esquemas integrados son tan buenas como las de
los esquemas convencionales. En este sentido una
seleccin adecuada de estrategias de control y
sintonizacin de controladores, permitirn que la
reduccin en el costo anual de operacin de los
esquemas trmicamente acoplados sea alcanzable en
la prctica, sin introducir problemas potenciales de
control.
La tendencias actuales en el campo de
esquemas de destilacin trmicamente acoplados
son: estudios de modos avanzados de control que
permitirn mejorar las propiedades de control,
separacin de sistemas azeotrpicos, destilacin
reactiva y estudios sobre el diseo y la optimizacin
mediante modelos basados en transferencia de masa
(Grossmann, 2005).
En el contexto anterior, el propsito de este
trabajo es el estudio de sistemas de reaccinseparacin, en esquemas de destilacin trmicamente
acoplados (TCDS) con columnas laterales como los
indicados en la Fig. 1.
En el diseo y optimizacin de los esquemas
de destilacin con acoplamiento trmico es muy
importante minimizar el consumo de energa con
respecto al flujo de interconexin. Un mtodo de
diseo para TCDS usa las secciones de platos de las
secuencias de destilacin convencionales directa e
indirecta. Por ejemplo, en el caso de la TCDS con
rectificador lateral de la Fig. 1a, el diseo se muestra
de manera conceptual en la Fig. 2a. En el caso de la

secuencia convencional directa, la primera columna
tiene secciones de rectificacin y agotamiento
denotadas por 1 y 2, para la segunda columna se
tienen las secciones de rectificacin y agotamiento
llamadas 3 y 4.
Fig. 1. Esquemas de destilacin trmicamente

acoplados para separaciones ternarias.
Fig. 2. Diseo de TCDS partiendo de una secuencia

de destilacin convencional: (a) TCDS con
rectificador lateral, (b) TCDS con agotador lateral.
La TCDS con rectificador lateral, se genera al
mover la seccin de agotamiento 4 de la segunda
columna al fondo de la primera columna, con lo cual
se genera una columna con secciones 1, 2 y 4. La
seccin de rectificacin 3 se acopla mediante un par
de corrientes de vapor y lquido a la columna de tres
secciones. En el caso de TCDS con rectificador
lateral, el flujo de vapor (VF) se vara hasta que se
detecta el mnimo consumo de energa (Hernndez y
Jimnez, 1996). Para la TCDS con agotador lateral el
flujo de interconexin LF es en fase lquida y el
diseo se obtiene de manera similar.
138
2. Metodologa
Usando el mtodo de diseo descrito en
Hernndez y Jimnez (1996) se procedi al estudio
del proceso industrial de produccin de acetato de
etilo. El proceso convencional industrial mostrado en
la Fig. 3, consta de una columna-reactor en la cual se
introducen corrientes de alcohol etlico y cido
actico en el rehervidor de la columna, que contiene
cido p-toluensulfnico como catalizador. En la
columna- reactor se produce acetato de etilo y agua
que se destilan juntamente con el alcohol y el cido
que no reaccionaron. Se pasan luego a un tanque
colector de destilados, en el que se introduce una
corriente de agua con el propsito de formar una
capa orgnica y otra acuosa. Parte de la capa
orgnica se regresa a la columna-reactor como
reflujo, la parte restante se introduce a una columna
de purificacin, en la que se obtiene una corriente de
acetato de etilo con una pureza de 99.5 % en peso.
La capa acuosa se introduce a una columna de
destilacin con el propsito de recuperar una parte
del alcohol y regresarlo a la fase de reaccin. Dos
problemas importantes del proceso son: primero, se
trata de una reaccin que est gobernada por el
equilibrio y, segundo, se forman azotropos que
dificultan la posterior purificacin del acetato y
recuperacin de alcohol y cido.
Dentro del estudio realizado y motivo del
presente artculo, se consideraron esquemas TCDS
con rectificador o agotador laterales (Fig. 1), con el
fin de estudiar la etapa de reaccin y determinar el
efecto que tienen en la posterior purificacin de las
corrientes de producto y reciclo. Al final, se
presentan los resultados preliminares para el estudio
de este sistema reactivo en una columna tipo Petlyuk.
acoplados y reactivos no existen. En la mayora de

los trabajos reportados para procesos de
esterificacin en columnas de destilacin (Chien y
col., 2005; Taylor y Krishna, 2000), el diseo ptimo
se obtiene mediante un procedimiento de
optimizacin en el cual se vara el nmero de platos,
plato de alimentacin, relacin de reflujo y carga
trmica hasta lograr minimizar una funcin objetivo.
En el presente estudio se procedi a proponer una
estructura de platos de la secuencia con acoplamiento
trmico para estudiar los parmetros operacionales.
Esta propuesta se basa en las columnas reales usadas
en un proceso comercial de produccin de acetato de
etilo. Bsicamente se considera que los reactivos se
introducen en el rehervidor, el cual se puede modelar
mediante un reactor de equilibrio, intentando
desplazar el equilibrio hacia la formacin del acetato
de etilo, al extraer la mayor cantidad de agua,
mediante una corriente lateral. Para el modelado de
estos sistemas se emplean las ecuaciones MESH,
balances de materia, energa, relaciones de equilibrio
y restricciones como sumatorias, reportadas en el
apndice A. Los esquemas analizados en esta
modalidad son similares al mostrado en la Fig. 4,
para la TCDS con agotador lateral, de acuerdo a la
ecuacin:
A+ B C + D
A Alcohol Etlico
B cido Actico
C Acetato de Etilo
D Agua
Fig. 4. TCDS reactivo con agotador lateral.
Fig. 3. Esquema convencional de produccin de

acetato de etilo.
3. Resultados
Como se mencion anteriormente ya existen
mtodos de diseo y optimizacin de TCDS; sin
embargo, para el caso de sistemas trmicamente
El estudio se realiz usando el simulador de

procesos Aspen Plus 11.1TM, empleando para la
prediccin del equilibrio de fases, las ecuaciones
NRTL y Redlich-Kwong y considerando que la
reaccin se lleva a cabo en fase lquida.
En primera instancia se encontr que la
columna lateral juega un papel muy importante en la
eficiencia global del proceso de reaccin-separacin,
ya que extrae el agua que se forma en la reaccin,
con lo cual se incrementa la conversin en la
reaccin, al remover directamente uno de los
productos de la reaccin. La conversin lograda para
la secuencia acoplada directa fue del 68%, mientras
que para la secuencia indirecta, la conversin
alcanzada fue del 80%. Tambin, la remocin del
139
agua resulta de gran importancia debido a que se

evita la formacin de los azetropos que forma el
acetato de etilo y el alcohol etlico con el agua.
En el caso de la secuencia trmicamente
acoplada indirecta, la Fig. 5 muestra la remocin del
agua en la columna lateral al variar el flujo de
interconexin fase lquida y el producto de fondos
del agotador lateral. En esta grfica se nota que la
mxima cantidad de agua que se puede extraer,
corresponde a un flujo de fondos de 22.7 kg-mol/h.
puede observar la cantidad de agua extrada y la

carga trmica respectivamente, al variar el flujo de
interconexin de vapor hacia el rectificador lateral.
De acuerdo a la Fig. 8, la carga trmica
mnima se presenta con un flujo de vapor de 27.24
kg-mol/h, y un flujo de destilado de 22.7 kg-mol/h
con dichos valores se extrae la mxima cantidad de
agua segn la Fig. 7 y correspondiente a 15.9 kgmol/h.
17
16
BR
(Kmol/hr)
19
11.35
17
13.62
15
15.89
18.16
13
20.43
11
22.7
15
Agua Extrada(Kmol/hr)
Agua Extrada (Kmol/hr)
21
DR
(Kmol/hr)
14
13
13.62
15.89
12
18.16
22.7
11
10
10
20
30
40
50
20
60
22
24
26
30
32
34
En el desempeo energtico de las TCDS, el

flujo de interconexin tiene un efecto muy
importante. En este sentido, la Fig. 6 muestra la
variacin del consumo energtico en funcin del
flujo de interconexin para diferentes valores de
producto de fondos, para el caso de la TCDS con
agotador lateral. Es posible observar en la Fig. 6 que
la carga trmica menor corresponde a un producto de
fondos de 22.7 kg-mol/h y un flujo de interconexin
de 68.1 kg-mol/h. La cantidad mxima de agua se
obtuvo con los valores anteriores, lo cual indica que
los valores ptimos para operar la TCDS con
agotador lateral son los encontrados anteriormente
debido a que cumplen los objetivos de extraer la
mayor cantidad de agua y, a su vez, la mnima carga
trmica suministrada a los rehervidores.
1,900,000.00
1,800,000.00
BR
(Kmol/hr)
1,700,000.00
11.35
1,600,000.00
13.62
1,500,000.00
15.89
1,400,000.00
18.16
1,300,000.00
20.43
Fig. 7. Efecto del flujo de interconexin en el agua

extrada, para la TCDS con rectificador lateral.
2.20E+06
2.00E+06
Carga Trmica ( W )
Fig. 5. Efecto del flujo de interconexin en el agua

extrada, para la TCDS con agotador lateral.
DR
(Kmol/hr)
1.80E+06
13.62
15.89
18.16
22.7
1.60E+06
1.40E+06
1.20E+06
1.00E+06
20
22
24
26
28
30
32
34
FL (Kmol/hr)
Fig. 8. Efecto del flujo de interconexin en la carga

trmica, para la TCDS.
Cabe sealar que actualmente se est
analizando la columna Petlyuk y se han encontrado
resultados interesantes. Por ejemplo, respecto a la
carga trmica, lo cual puede observarse en la Fig. 9,
donde el consumo de energa es menor para este
esquema trmicamente acoplado. Lo anterior es
prometedor para posteriores anlisis.
698550
VF
(Kgmol/h)
22.7
1,200,000.00
698500
27.24
1,100,000.00
1,000,000.00
10
20
30
40
50
60
70
FL (Kmol/hr)
Fig. 6. Efecto del flujo de interconexin en el

consumo de energa, para la TCDS con agotador
lateral.
Carga Trmica ( W )
Carga Trmica (W)
28
FV (Kmol/hr)
FL (Kmol/hr)
698450
22.7
18.16
698400
15.89
13.62
698350
11.35
9.08
698300
4.54
698250
698200
0
10
15
20
25
30
FL (kgmol/hr)
De igual manera se explor la secuencia

trmicamente acoplada directa. En las Figs. 7 y 8 se
Fig. 9. Efecto del flujo de interconexin en la carga

trmica, para la columna Petlyuk.
140
Conclusiones
Al analizar las secuencias trmicamente
acopladas para la produccin de acetato de etilo, la
mayor cantidad de acetato fue producido por la
TCDS con agotador lateral, logrando un rendimiento
de la reaccin del 80%, mientras que la secuencia
trmicamente acoplada directa tuvo un rendimiento
del 68 %; lo anterior consecuencia directa de la
cantidad de agua extrada por la columna lateral en
cada configuracin.
De acuerdo a los resultados anteriores, se
puede concluir que los esquemas trmicamente
acoplados son factibles de aplicar en este proceso
debido a que aumentan el rendimiento de la reaccin,
al desplazar el equilibrio hacia la formacin del
acetato, mediante la extraccin del agua producida.
Adems con tal accin, se evita la formacin de dos
fases lquidas y de azetropos en el producto de la
columna de esterificacin.
Agradecimientos
Este proyecto se est desarrollando en parte
gracias al apoyo recibido de PROMEP y la
Universidad de Guanajuato.
Nomenclatura
BR Producto de fondos en agotador lateral.
c
Nmero de componentes.
DR Destilado retornado en rectificador lateral.
f
Flujo de alimentacin a la etapa de equilibrio.
h
Entalpa del lquido.
H
Entalpa del vapor.
i
Componente i.
j
Etapa de equilibrio.
k
Constante de equilibrio.
LF
Flujo de lquido.
l
Flujo de lquido que sale de la etapa.
M
Moles de lquido retenido en la etapa.
Q
Carga trmica.
TCDS Sistema de destilacin trmicamente
acoplados.
u
Corriente lateral de lquido que sale de la
etapa de equilibrio.
VF Flujo de vapor.
v
Flujo de vapor que sale de la etapa de
equilibrio.
w
Flujo lateral de vapor extrado de la etapa de
equilibrio.
x
Fraccin mol en la fase lquida.
y
Fraccin mol en la fase vapor.
H Entalpa de reaccin.
Referencias
Configuration. Journal of Process Control 15,

435.
Flores, O. A.; Crdenas, J.C.; Hernndez, S. y RicoRamrez, V. (2003). Thermodynamic Analysis
of Thermally Coupled Distillation Sequences.
Industrial and Engineering Chemistry
Research 42, 5940.
Grossmann, I. E. (2005). Comunicacin Personal.
Hernndez, S. y Jimnez, A. (1996). Design of
Optimal
Thermally-Coupled
Distillation
Systems
Using
a
Dynamic
Model.
Transactions of the Institution of Chemical
Engineers 74, 357.
Hernndez, S. y Jimnez, A. (1999a) Design of
Energy-Efficient Petlyuk Systems. Computers
and Chemical Engineering 23, 1005.
Hernndez, S. y Jimnez, A. (1999b). Controllability
Analysis of Thermally Coupled Distillation
Systems.
Industrial
and
Engineering
Chemistry Research 38, 3957.
Hernndez, S.; Pereira-Pech, S.; Jimnez, A. y RicoRamrez, V. (2003). Energy efficiency of an
indirect thermally coupled distillation
sequence. The Canadian Journal of Chemical
Engineering 81, 1087.
Kaibel, G. y Schoenmarkers, H. (2002). Process
Synthesis and Design in Industrial Practice.
Proc. ESCAPE-12 (Computer Aided Process
Engineering,10), Eds. J. Grievink and J.V.
Schijndel, Elsevier, Amsterdam 9.
Rong, B. G.; Kraslawski, A. y Turunen, I. (2003).
Synthesis of Functionally Distinct Thermally
Coupled Configurations for Quaternary
Distillation. Industrial and Engineering
Chemistry Research 42, 1204.
Segovia-Hernndez, J.G.; Hernndez, S.; RicoRamrez, V. y Jimnez, A. (2004). A
Comparison of the Feedback Control
Behavior between Thermally Coupled and
Conventional
Distillation
Schemes.
Computers and Chemical Engineering 28,
811.
Taylor, R. y Krishna, R. (2000). Modelling Reactive
Distillation. Chemical Engineering Science
55, 5183.
Tedder, D. W. y Rudd, D. F. (1978) Parametric
Studies in Industrial Distillation: Part I.
Design Comparisons. AIChE Journal 24, 303.
Triantafyllou, C. y Smith, R. (1992). The Design and
Optimisation of Fully Thermally Coupled
Distillation Columns. Transactions of the
Institution of Chemical Engineers 70, 118.
Wolff, E. A. y Skogestad, S. (1995). Operation of
Integrated
Three-Product
(Petlyuk)
Distillation
Columns
Industrial
and
Engineering Chemistry Research 34, 2094.
Chien, I.-L.; Teng, Y.-P.; Huang, H.-P. y Tang, Y. T.

(2005). Design and Control of an Ethyl
Acetate Process: Coupled Reactor/Column
141
Apndice A
El modelo del esquema reactivo se puede

obtener en base a la etapa de equilibrio mostrada en
la Fig. A.1.
Las ecuaciones que describen el modelo son
las siguientes:
Balances de materia por componente
f i , j + l i , j -1 + vi , j +1 l i , j (1 + S L J )
(A.1)
vi , j (1 + SVJ ) M j ri , j = 0
SVJ =
S LJ =
wj
Vj
uj
Lj
Equilibrio de fases.
yi , j = ki , j xi , j
y
i =1
=0
=0
(A.6)
Balance de energa.
f j h j + l j 1h j 1 + v j +1 H j +1 (v j + w j ) H j
(A.7)
(l j + u j ) h j + Q r H rxn M j = 0
vij
lij1
wij
(A.2)
fij
(A.3)
uij
(A.4)
Restricciones como sumatorias.

c
i =1
(A.5)
lij
Fig. A.1 Etapa de equilibrio.
vij+1
142
Francisco Snchez Mares

Uso de Arreglos Matriciales para obtener Secuencias ptimas de Separacin de Mezclas Multicomponentes
Conciencia Tecnolgica, nm. 27-30, 2005
Mxico
Conciencia Tecnolgica,
contec@correo.ita.mx
Mxico
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
Uso de Arreglos Matriciales para obtener Secuencias ptimas de Separacin

de Mezclas Multicomponentes
Divulgacin tecnolgica
Ing. Francisco Snchez Mares
Maestra en Ciencias en Ingeniera Qumica
fsanchez96@hotmail.com
Av. Lpez Mateos 1801 Ote., Esq. Av. Tecnolgico Aguascalientes, Ags.,
Tel. 9105002 Ext. 103 fax (0149) 700423 C.P.20256
Resumen
Introduccin
La separacin de sistemas multicomponentes ha

sido de gran importancia para la industria. Un ejemplo
es la petroqumica, la separacin del petrleo (mezcla
de hasta 277 compuestos de hidrocarburos) tiene como
resultado gran cantidad de subproductos con distinto
valor agregado. Una alternativa ampliamente utilizada
en la separacin de sistemas multicomponentes son los
trenes de separacin debido a que son eficientes y
econmicos, sin embargo, el encontrar cual es la
secuencia ptima es un problema que ha atrado la
atencin de diversos investigadores.
El presente trabajo muestra una alternativa
novedosa en la sntesis de secuencias de separacin de
mezclas multicomponentes, basndose en dos reglas
heursticas bsicas: tratar de separar los componentes
ms abundantes primero y hacer la separacin ms
difcil al ltimo.
Estas dos reglas se conjugan en un arreglo matricial
que las har interactuar entre si, de tal manera que la
respuesta del arreglo sea la ptima para los
componentes de la mezcla con base en las reglas antes
mencionadas, sin embargo, esto no garantiza que sea
la mejor en base a economa.
La secuencia propuesta puede ser analizada en base
a sus costos, y proponer una nueva secuencia a partir
del mtodo evolutivo, al tener las dos secuencias se
hace una comparacin en costos y se concluye si es la
mejor o todava existe otra ms econmica. Si existe
otra, entonces seguiremos evolucionando la secuencia
hasta obtener la ptima.
Al aplicar los arreglos matriciales en los ejemplos
presentados en este trabajo se obtiene en la mayora de
los casos la mejor secuencia, pero tambin hay algunos
donde se presenta otra alternativa muy cercana a la
ptima.
An cuando con frecuencia se puede disear un

separador complejo para obtener todos los productos
deseados de mezclas multicomponentes, generalmente
se utiliza una secuencia de separadores simples debido
a que frecuentemente resulta ms econmico que un
separador complejo.
Al aumentar al nmero de componentes de una
mezcla, el problema de sntesis de separacin se
complica por la naturaleza combinatorial del
problema. El nmero de subproblemas de separacin y
por lo tanto la cantidad de posibles secuencias que se
pueden generar dependen del nmero de componentes
de la mezcla original. La Tabla 1 muestra el nmero de
secuencias segn el nmero de componentes de la
mezcla.
Palabras clave
Reglas heursticas,
evolutivo.
arreglo
matricial,
mtodo
Nmero de
Nmero de
Nmero de
componentes
subgrupos
subproblemas
Nmero de
posibles
[N]
[N(N+1)/2] - N
de separacin
secuencias
[(N - 1)(N)(N + 1)/6]
[2(N-1)]!/N!(N - 1)!
10
10
20
14
15
35
42
21
56
132
28
84
429
36
120
1,430
10
45
165
4,862
Tabla 1. Nmero de posibles secuencias de separacin

en funcin del nmero de componentes.
Por lo tanto el problema se complica
significativamente al considerar ms de una posible
operacin unitaria para la separacin de los
componentes. Es necesario entonces tener un
procedimiento eficiente que detecte la secuencia de
separacin ptima (o cercana a la ptima), y de ser
posible, con un mnimo de esfuerzo computacional.
Aunque existen distintos enfoques para encontrar

una alternativa ptima como lo son mtodos
heursticos y los mtodos algortmicos, su principal
desventaja es que se tiene que analizar secuencia por
secuencia y determinar su costo para cada una, esto es
una actividad que puede ocupar algo de tiempo, aun
cuando con la ayuda de la computadora se reduce al
mximo este trabajo, es necesario buscar otras
alternativas distintas.
Fundamentos tericos
secuencia es preciso evaluar su coste con el fin de

determinar si es mejor. Las siguientes reglas
evolutivas de Stephanopoulos y Westerberg
proporcionan una aproximacin sistemtica.
1. Intercambiar las posiciones relativas de dos
separadores adyacentes.
2. Sustituir una determinada separacin usando el
mtodo I por la separacin del mtodo II.
Con frecuencia la sntesis evolutiva puede conducir
a la secuencia ptima sin necesidad de examinar todas
las secuencias. La eficacia depende de la estrategia
utilizada para dirigir las modificaciones.
1. Tcnicas heursticas y sntesis evolutiva.

2. El mtodo rpido.
Los mtodos heursticos, que abordan la solucin
de un problema por medio de reglas plausibles, pero
no infalibles, se utilizan ampliamente para evitar tener
que examinar todas las secuencias posibles con el fin
de encontrar las disposiciones ptimas y casi ptimas.
Utilizando mtodos heursticos se pueden
establecer rpidamente buenas secuencias sin disear y
estimar el coste del equipo aun en el caso de que haya
un gran nmero de componentes a separar.
Se han propuesto numerosas reglas heursticas para
secuencias sencillas compuestas de columnas de
destilacin ordinaria. Las reglas heursticas ms tiles
se pueden seleccionar, de acuerdo con Seader y
Westerberg[1] para desarrollar secuencias partiendo de
la alimentacin del proceso, en la forma que se indica
a continuacin.
1. Cuando los componentes adyacentes ordenados
de la alimentacin del proceso varan grandemente en
la volatilidad relativa, la secuencia de las escisiones se
establece en el orden de volatilidad relativa
decreciente.
2. La secuencia de las escisiones para la separacin
de componentes se establece en el orden del porcentaje
molar decreciente de la alimentacin del proceso,
cuando dicho porcentaje varia grandemente pero en
cambio la variacin de la volatilidad relativa es
pequea.
3. Cuando ni la volatilidad relativa ni el porcentaje
molar en la alimentacin varan grandemente, deben
separarse los componentes uno a uno como productos
de cabeza.
Estas tres reglas heursticas estn de acuerdo con
numerosas observaciones relacionadas con el efecto de
la presencia de componentes no clave sobre el coste de
escisin de dos componentes clave.
Una vez que se ha desarrollado una buena
secuencia de separacin a partir de reglas heursticas,
puede tratar de mejorarse mediante sntesis evolutiva
en la forma descrita por Seader y Westerberg. Este
mtodo consiste en desplazarse desde una secuencia de
partida hacia secuencias cada vez mejores a travs de
una sucesin de pequeas modificaciones. Para cada
Con la aplicacin de mtodos algortmicos se han

detectado secuencias de separacin ptimas. Al revisar
las reglas heursticas y su posible ingerencia para
explicar la estructura de secuencias ptimas, puede
encontrarse que en la mayora de los casos las
secuencias que minimizan el consumo de servicios
tienden a minimizar el costo anual del sistema.
Entonces, se puede utilizar una regla heurstica que
cumpla con este objetivo y es de esperarse que la
estructura que se genere sea ptima o una excelente
aproximacin al ptimo.
Un mtodo rpido[2] para la sntesis de trenes de
separacin consiste en buscar la alternativa que
minimice los flujos de alimentacin a la secuencia:
Min
Fi
[1]
donde Fi es la alimentacin a cada una de las

columnas que contiene la secuencia que separa n+1
componentes. Esta funcin objetivo es parte de un
problema de optimizacin bastante simple, que puede
resolverse manualmente por principios como
programacin dinmica, mtodos evolutivos o
inclusive mediante simple inspeccin en algunos
casos. Una vez que la estructura ptima se detecta,
cada columna deber disearse, para lo cual las
variables de presin y razn de reflujo a su vez pueden
establecerse por las heursticas ya tratadas.
3. Heursticas aplicadas.
Este trabajo se basa en las dos primeras reglas
heursticas propuestas por Seader y Westerberg, para
desarrollar secuencias partiendo de la alimentacin del
proceso.
El motivo de seleccionar estas heursticas:
a)
Provienen de la experiencia de diseadores.
b) Son fciles de medir, ya que la volatilidad y

el flujo manejado son conocidos desde un
principio.
c) El mtodo rpido es un mtodo basado
solamente en flujos y provee de buenas
secuencias al diseador.
4. Criterios para flujos y volatilidades.
En esta etapa se relaciona el problema de la

secuencia como si fuera un diseo de experimentos y
se llega a la conclusin de que esta comparacin puede
generar resultados interesantes, mismos que pueden
obtener una buena secuencia a partir de esta analoga.
En esta etapa del proyecto:
Al desear obtener una secuencia directa, se tienen

que implicar los cortes en cada arreglo, por ejemplo si
se desea evaluar el corte A/B se necesita tener como
datos la razn entre la volatilidad () de A y la de B, es
decir:
Corte ij = (
i/
j)
[2]
Volatilidad relativa ij = (
i/
j)
[3]
Para los flujos en el corte A/B se presentan dos
opciones para evaluarlo, la primera es el evaluar la
razn entre dos flujos, por ejemplo A/B.
Corte ij = Flujo i/ Flujo j
2. Diseo de los arreglos matriciales.
a). Se intenta aplicar arreglos ortogonales ya

existentes a la sntesis de procesos, relacionando su
nmero de componentes con los factores del arreglo y
el nmero de corridas con las alternativas, esto genera
un intento fallido pero presenta resultados interesantes.
b). Cuando se observa que existe la posibilidad de
generar una buena secuencia, pero que no existe un
arreglo adecuado para nuestro problema, se realiza el
diseo de los arreglos matriciales adecuados.
3. Pruebas con sistema de tres componentes en
arreglos matriciales.
[4]
Sin embargo, si obtenemos la razn de A/B y si A

tiene un caudal muy grande y B uno muy pequeo
entonces al dividirlos se obtendr un resultado muy
alto, mismo que engaara al arreglo matricial,
dicindole que posiblemente este corte presenta el
flujo ms grande, de tal manera que la propuesta
efectuada se descarta y se utiliza la siguiente:
Corte ij = (Flujo de i + Flujo de j)/suma de todos los
flujos involucrados
[5]
De esta manera el arreglo matricial tiene una buena
percepcin del orden de las magnitudes de flujos de
cada corte propuesto.
Materiales y mtodos
El proyecto se basa en el siguiente objetivo:
generar una alternativa novedosa, rpida y confiable
en la obtencin de la secuencia de separacin ptima
de una mezcla multicomponentes.
El desarrollo de este proyecto se lleva a cabo de la
siguiente forma:
1. Revisin de la teora de sntesis de secuencias de
separacin.
En esta parte se analizan todas las aportaciones a la
sntesis de procesos de distintos autores tanto en
mtodos algortmicos como heursticos.
a). Teniendo los arreglos adecuados, se procede a

evaluar un problema de cinco secuencias, pero por
separado, es decir, primero para volatilidades y se
observ que gener la secuencia ptima para esa
condicin. Posteriormente se evala la condicin de
flujos y tambin se genera la secuencia ptima.
b). Sabiendo que el arreglo funcionaba
perfectamente, se procede unir las dos condiciones en
nuevo arreglo matricial, mismo que vuelve a generar
una secuencia ptima idntica.
4. Desarrollo de un programa de cmputo para
evaluar sistemas de ms de tres componentes con
arreglos matriciales.
Al observar que el arreglo generaba secuencias
aceptables, se procede a la programacin de varios
arreglos, para facilitar el trabajo de experimentacin y
anlisis de resultados.
5. Resultados.
En esta etapa se evalan dos problemas
encontrados en la literatura y se comprueba la
efectividad de los arreglos generados para encontrar
buenas secuencias. Se genera un paquete de
programas, cada programa evala secuencias, dando
distinto peso a cada opcin (volatilidad y flujos) y se
obtienen opciones interesantes. Estos programas se
realizan en C++ Builder.
Resultados y discusin
A
Para probar la efectividad del mtodo se

seleccionaron varios problemas de la literatura, los
cuales han sido ampliamente estudiados por diversos
autores. Sin embargo debido a la brevedad del
documento solo se exponen dos ejemplos.
1. Separacin de parafinas ligeras por destilacin

ordinaria.
Considere la separacin de una mezcla de cinco
hidrocarburos ligeros en componentes puros por
destilacin ordinaria estudiada por Heaven[3] (1969),
Usando una base de clculo de 100 mol/hora. Ver la
Tabla 2.
Especie
Fraccin
Corte
Flujos
A: Propano
0.5
A/B
0.33334
2.0
1.33
A
B
A
B
C
A
B
C
D
E
B
C
D
D
E
E
Figura 1. Secuencia ptima. (50% volatilidad y 50%

flujos).
Volatilidad
relativa
B: Iso butano
0.15
B/C
0.6
C: N butano
0.25
C/D
1.25
2.4
D: Iso pentano
0.20
D/E
0.571428571
1.25
E: N pentano
0.35
Tabla 2. Datos para el problema de Heaven.
Y si se trabaja 10% volatilidad y 90% flujos se

obtiene:
A/B = 6.8433337(segundo en efectuar)
B/C = 6.3550004 (ltimo corte)
C/D = 7.4941401 (efectuar primero)
D/E = 6.6146802(tercero en efectuar)
Los datos se introducen en el arreglo matricial

mostrado en la Tabla 3.
A
Corte
Respuesta
A/B
Y1
B/C
Y2
C/D
Y3
D/E
Y4
A/B
Y5
B/C
Y6
C/D
Y7
D/E
Y8
Tabla 3. Arreglo matricial L8(45).

En el arreglo matricial se observan nmeros
uno y dos, el nmero uno se refiere a flujos y el dos a
valores de volatilidad.
Evaluando el arreglo matricial dndole el mismo
peso tanto a flujos como a volatilidades (50% flujos y
50% volatilidad) se obtienen los siguientes resultados:
A/B = 10.566699 (ltimo en realizar)
B/C = 12.710000 (segundo en efectuar)
C/D = 14.301428(efectuar primero)
D/E = 11.095714 (tercero en efectuar)
La Figura 1 muestra la secuencia para la relacin
50/50 y la Figura 2 presenta el tren de separacin
obtenido para una relacin 10% flujos y 90% flujos.
A
B
C
D
E
A
B
C
B
B
C
D
C
D
E
E
Figura 2. Segunda mejor secuencia. (10% volatilidad

y 90% flujos).
Cabe mencionar que la primer secuencia esta
reportada en la literatura como la mejor de las catorce
posibles con un costo de 858,780 USD y la segunda
como otra muy buena opcin solo 0.53% arriba de la
primera.
2. Separacin de olefinas ligeras y parafinas por
destilacin ordinaria.
Considere la separacin de una mezcla de olefinas
ligeras y parafinas por destilacin ordinaria. La mezcla
de alimentacin es presentada por Thompson y King
(1972[3], utilizando una base de clculo de 100

mol/hora. La Tabla 4 presenta los datos del problema.
Especie
Fraccin
Corte
Flujos
A: Etano
0.20
A/B
1.33334
3.5
B: Propileno
0.15
B/C
0.75
1.2
Volatilidad
relativa
C: Propano
0.20
C/D
1.33334
2.7
D: 1 buteno
0.15
D/E
1.21
E: N butano
0.15
E/F
3.0
F: N pentano
0.15
Tabla 4. Datos para el problema de

King.
A
A
B
C
A
B
C
D
E
F
Thompson y
B
C
C
D
E
F
D
E
Para este problema se utiliz el arreglo matricial

mostrado en la Tabla 5.
Figura 4. Secuencia (10% volatilidad y 90% flujos).

Conclusiones
Corte
Respuesta
A/B
Y1
B/C
Y2
C/D
Y3
D/E
Y4
E/F
Y5
A/B
Y6
B/C
Y7
C/D
Y8
D/E
Y9
E/F
Y10
El presente trabajo muestra una forma rpida de

obtener secuencias de separacin econmicas a partir
de dos reglas heursticas importantes, basadas en la
experiencia de diseadores de secuencias de
separacin. Las mejores secuencias generadas por este
mtodo son las que tienen igual porcentaje de
importancia tanto en flujos como en volatilidades
(50% y 50%), aunque tambin se obtienen excelentes
secuencias al darle mayor peso a los flujos (90% flujos
y 10% volatilidades).
Tabla 5. Arreglo matricial L10 (56).

La Figura 3 muestra la secuencia generada para la
relacin 50% volatilidad y 50% flujos.
A
A
B
C
D
E
F
B
B
C
B
C
D
E
F
El autor de este artculo agradece a los profesores Dr.

Omar Jair Purata Sifuentes y al Dr. Adrin Bonilla
Petriciolet, de la Maestra en Ciencias en Ingeniera
Qumica del Instituto Tecnolgico de Aguascalientes
por sus brillantes aportaciones a este proyecto.
Bibliografa
[1] Ernest J. Henley y J.D. Seader (2003);

Operaciones de separacin por etapas de equilibrio en
ingeniera qumica. Editorial Revert. Espaa
D
E
D
E
F
Agradecimientos
Figura 3. Secuencia (50% volatilidad y 50% flujos).

Esta secuencia viene reportada como una de las
cuatro mejores secuencias obtenidas para este
problema por Nadgir y Liu[3].
La secuencia generada para 10% volatilidad y 90%
flujos se presenta en la Figura 4.
[2] Arturo Jimnez Gutirrez (2003); Diseo de

procesos en ingeniera qumica. Editorial Revert.
Espaa.
[3] V. M. Nadgir y Y. A. Liu (1983); A simple
heuristic method for systematic sntesis of inicial
sequences for multicomponent separations.
AIChE Journal (Vol 29, No. 6).
TECNOLGICO DE ESTUDIOS SUPERIORES DE SAN FELIPE DEL PROGRESO

INGENIERA QUMICA
Aplicacin de procesos de separacin de mezclas multicomponentes
Determinacin de la definicin del polinomio de lquido-vapor

separado de mezclas multicomponentes
PRESENTAN:
Cortez Alvarado Christian Alhi
Garca Lpez Daniela Sarahy
Snchez Nava Dayra Margarita
Zubieta Aguirre Jessica
Catedrtico :
I.Q Guadalupe Lpez Garca
ABRIL 2015
INTRODUCCIN
Rosenberg (1963, 1977) mostr que la determinacin dividida lquidovapor en un problema

isotrmico-isobrico flash puede ser
representada por un polinomio de grado n - 1 con respecto a la
vaporizacin instantnea, donde N es el nmero de componentes
presentes en la mezcla.
En su trabajo, las expresiones para los coeficientes de este polinomio

se han presentado durante un mximo de 6 productos, que se puede
utilizar durante un mximo de 14 componentes
Las expresiones para el clculo de los coeficientes de este polinomio

se derivan para cualquier nmero de componentes utilizando la
ecuacin Rachford y Rice.
2
PLANTEAMIENTO DEL PROBLEMA
Considere el proceso de separacin flash vapor-lquido
Representacin esquemtica de una separacin flash.
Donde:
- relacin de la velocidad de flujo molar total del vapor a la de la

corriente de alimentacin, que generalmente se llama vaporizacin
instantnea (flash);
N es el nmero de especies en la mezcla
xi, yi, y zi son las fracciones molares de componente i (i = 1,,N); en

la fase lquida, en fase de vapor, y corriente de alimentacin,
respectivamente
Ki es el factor de separacin del componente i que se define como yi

= xi.
En general, los valores K se obtienen a partir de un modelo

termodinmico y son una funcin de la temperatura T, presin P, las
fracciones molares lquidos xi, y las fracciones molares de vapor yi se
puede escribir como:
Donde la Eq. (6) es la ecuacin Rachford-Rice
Ecuacin Rachford-Rice.
Una funcin objetivo que a menudo permite una convergencia ms

rpida para la determinacin de la fraccin vaporizada en una
corriente a temperatura y presin constante.
Un simple anlisis muestra que hay N + 2 grados de libertad, que por

lo general se fijan mediante la especificacin N fracciones molares de
alimentacin zi, la temperatura T, y la presin P, del proceso de
separacin flash.
Suponiendo que las especificaciones son elegidas por debajo del

estado crtico de una mezcla no azeotrpica (para evitar los efectos
retrgrados) para definir al menos una solucin de vapor-lquido real
(cantidades positivas de material en cada fase). Entonces tenemos un
problema isotrmico-isobrico.
Este estudio es importante porque est delimitado

entre cero y uno.
Los valores de fuera de la regin cero-uno implican
que o bien el lquido el vapor se alimenta en vez de ser
representado como un producto.
Por
lo tanto, si el lquido se suministra a la torre, adems

de la alimentacin, como consecuencia ser mayor que
uno, y si el vapor se alimenta, entonces ser negativo.
Estos
casos se han identificado como el flash negativo.
El
propsito principal de este concepto es ser capaz de

obtener la convergencia en la regin de una sola fase.
7
VAPORIZACIN FLASH (INSTANTANEA) DEL

POLINOMIO
De la ecuacin Rachford-Rice (Eq. 6), una ecuacin polinmica en se puede

derivar, multiplicando a travs de los denominadores de los trminos
individuales. Cuando Ki-1 se define como Ri, esta relacin es:
Se puede volver a escribir como un polinomio de grado N-1
Las races de este polinomio que se encuentran entre cero y uno van a ser los
valores deseados de .
La raz o races se pueden encontrar por cualquier mtodo convencional. El

mtodo de Newton podra ser utilizado.
CUESTIONES POLINOMICAS
Dado que, slo estamos interesados en la raz de [0,1]; si esta existe, se deben
cumplir las siguientes condiciones:
Estas condiciones son equivalentes a las observaciones para la existencia de una

fraccin de vapor-lquido. Slo se produce cuando el punto de operacin est entre
el roco y el punto de burbuja.
Que se puede obtener mediante la evaluacin de la ecuacin. (6) para

0, y = 1, respectivamente.
COEFICIENTES POLINOMIALES
Definiendo los siguientes vectores:
10
Donde:
Para k = 1,,N. Por lo tanto, los coeficientes

del polinomio son dados por
11
EJEMPLO:
Considere una mezcla de 6 componentes como se describe en la Tabla 1.
Tabla 1. Caractersticas de la mezcla a una presin de 2068 kPa y a una temperatura

de 400 K, para caso de estudio nmero 1.
12
Los valores Ki se asumen independientes de la composicin y se calculan

utilizando las expresiones polinmicas dadas.
De la ecuacin (6) y de la Tabla 1 se obtiene la siguiente expresin:
Y a partir de las ecuaciones. (13) a (17)
13
Las ecuaciones (18) y (19) se representan en la Fig. 2.
Las races del polinomio propuesto y la ecuacin Rachford-Rice.
14
Tenga en cuenta que si bien ambas ecuaciones tienen las mismas races, las
formas son diferentes. Eq. (18) es discontinua, mientras que la ecuacin. (19)
es continua y lisa. Las races son:
T = {- 5.00870, -1.68650, -0.38141, 0.72074, 2.41440}.
15
EJEMPLO 2
Ahora, consideremos la mezcla de 19 componentes. La composicin de la

alimentacin y los valores de Ki se dan en la Tabla 2.
16
En este caso se obtiene el siguiente polinomio:
17
Que tiene las siguientes races:
T = {- 11.093, -6.2422, -4.9901, -3.6911, -1.9426, -0.33238, -0.27747, 0.84614,

1.3847 5.7406 10-2 i, 1.6414 0.20527 i, 2.1747 0.21774 i, 2.8134, 3.5256,
4.8823, 6.9578}.
Como se mencion anteriormente, la nica raz que es de inters es

entre cero y uno.
En ambos casos de estudio, las desigualdades (9) y (10) se cumplen
18
Las figuras 3 y 4 muestran la ec. Rachford-Rice y las formas polinmicas para

el caso de estudio 1 y 2 respectivamente, en el intervalo de cero y uno.
Fig. 3 Forma del polinomio propuesto y la

ecuacin Rachford-Rice en el intervalo de
0 1 para Caso de estudio 1.
Fig. 4 Forma del polinomio propuesto y la

19 intervalo de
ecuacin Rachford-Rice en el
0 1 para Caso de estudio 2.
CONCLUSINES
Se ocupa del clculo de vaporizacin instantanea isobrica- isotrmica, que se

realiza extensamente en programas de simulacin de procesos asistidos por
ordenador.
Se han derivado expresiones para calcular los coeficientes necesarios para

definir un sistema de cualquier nmero de componentes para obtener el
polinomio que define el problema flash isotrmico- isobrica como una
funcin de .
Se encontr que la solucin propuesta es altamente estable y convergente

para cualquier tipo de vector de inicio y que el nmero de iteraciones es
afectado crticamente por ste vector y no necesariamente por la cantidad de
fases o nmero de componentes.
20
M. Vaca, A. Jimnez-Gutirrez, R. Monroy-Loperena

Aplicacin de tcnicas de colocacin en el diseo de columnas de destilacin petlyuk
Revista Mexicana de Ingeniera Qumica, vol. 5, nm. Su1, noviembre, 2006, pp. 121-127,
Mxico

amidiq@xanum.uam.mx
Iztapalapa
Mxico
Cmo citar?
Fascculo completo
Pgina de la revista
www.redalyc.org
REVISTA MEXICANA DE INGENIERA QUMICA Vol. 5, Supl. 1 (2006) 121-127
AMIDIQ
APLICACIN DE TCNICAS DE COLOCACIN EN EL DISEO

DE COLUMNAS DE DESTILACIN PETLYUK
APPLICATION OF COLLOCATION TECHNIQUES TO THE DESIGN
OF PETLYUK DISTILLATION COLUMNS
M. Vaca1,3, A. Jimnez-Gutirrez1* y R. Monroy-Loperena2
1
Departamento de Ingeniera Qumica, Instituto Tecnolgico de Celaya

Av. Tecnolgico y A. Garca Cubas, Celaya, 38010, Gto, Mxico.
2
ROMON, S.A. Paseo de los Pirules 124, Mxico, 04250, Mxico.
3
Divisin de Ingeniera Qumica y Bioqumica, Tecnolgico de Estudios Superiores de Ecatepec
Av. Tecnolgico S/N, Ecatepec, 55210, Mxico
Recibido 1 de Julio 2005; Aceptado 4 de Diciembre 2005
Resumen
Las tcnicas de colocacin ortogonal reducen la dimensionalidad de los modelos de procesos de separacin de
etapas; adems, transforman las variables discretas en continuas. Esta ltima caracterstica ha permitido la
aplicacin de mtodos estndar de programacin no lineal para disear columnas convencionales de destilacin. En
este trabajo se aplican tcnicas de colocacin para disear columnas de destilacin Petlyuk para la separacin de
mezclas ternarias. Utilizando como estimacin inicial el diseo preliminar obtenido mediante un mtodo corto, se
aplicaron tcnicas de colocacin para refinar la localizacin de las etapas de interconexin y de producto lateral as
como los flujos de interconexin que proporcionan purezas especificadas en los productos. El mtodo propuesto
proporciona una distribucin adecuada de etapas en el sistema Petlyuk, como se muestra para dos casos de estudio.
El diseo final se obtiene rpidamente mediante la metodologa presentada, con un esfuerzo notablemente menor
que el requerido cuando el ajuste del diseo preliminar se realiza mediante simulacin repetida y paquetes
comerciales. Los diseos fueron validados exitosamente mediante un simulador comercial.
Palabras clave: destilacin, columna Petlyuk, columnas trmicamente acopladas, mtodos de diseo.
Abstract
Orthogonal collocation techniques reduce the dimensionality of staged, separation process models; in addition, they
transform the discrete variables into continuous variables. This characteristic has allowed the application of
standard nonlinear-programming methods to the design of conventional distillation columns. In this work,
collocation techniques are applied to design Petlyuk distillation columns for the separation of ternary mixtures.
Using the preliminary design generated by a shortcut method, collocation techniques were applied to refine the
values for interconnecting and sidestream product stages as well as the interlinking flows that provide specified
product purities. The proposed method gives a proper distribution of stages for the Petlyuk system, as shown in two
cases of study. The final designs are easily obtained by means of the proposed methodology, with a significant
lower effort than that required when preliminary designs are corrected using recursive simulations. The designs
were successfully validated through the use of a commercial simulation package.
Keywords: distillation, Petlyuk column, thermally coupled columns, design methods.
1. Introduccin
El uso de acoplamiento trmico para
disminuir el consumo de energa en sistemas de
destilacin ha recibido mucha atencin en los
ltimos aos. La transferencia directa de energa
entre dos flujos a contracorriente provenientes de
columnas de destilacin adyacentes es el medio para
eliminar hervidores de las secuencias convencionales
de destilacin (SCD) y para reducir el consumo de
energa. Ejemplos de estos sistemas son las columnas

con rectificador lateral, con agotador lateral y las
columnas
Petlyuk,
llamadas
colectivamente
secuencias trmicamente acopladas (STA). Diversos
estudios han mostrado que los ahorros potenciales
pueden ir del 30 al 50% con respecto del consumo de
energa de las SCD (Glinos y Malone, 1988;
Hernndez y Jimnez, 1996, 1999; Jimnez y col.,
2003). Las STA no han sido empleadas ampliamente
a escala industrial; se ha dicho que el diseo y la
*Autor para la correspondencia: E-mail: arturo@iqcelaya.itc.mx

Fax: (01 461) 611 77 44
Publicado por la Academia Mexicana de Investigacin y Docencia en Ingeniera Qumica A.C.
121
M. Vaca y col./ Revista Mexicana de Ingeniera Qumica Vol. 5, Supl. 1 (2006) 121-127
operacin puede ser ms complejos que en el caso de

las SCD, debido al mayor nmero de grados de
libertad. Una columna Petlyuk (Fig. 1) tiene 11
grados de libertad, que deben especificarse explcita
o implcitamente a fin de determinar variables como
son: el nmero de etapas en cada una de sus 6
secciones ( N k , k = 1, , 6) , los flujos de las
corrientes de interconexin ( LS , VI ) , el reflujo de la
columna principal ( R) y el flujo de dos corrientes de

producto ( D, S , B) . Estas variables definen el
nmero de etapas en el prefraccionador, el nmero
de etapas en la columna principal, la etapa de
alimentacin, las etapas de interconexin superior e
inferior ( N L , NV ) y la etapa de extraccin de
producto lateral ( N S ) .
Fig. 1 Representacin esquemtica de una columna

Petlyuk
A la fecha se han empleado bsicamente dos
enfoques para disear columnas Petlyuk, que
Annakkou y Mizsey (1996) clasifican como enfoque
jerrquico y enfoque algortmico. En el enfoque
jerrquico, se efecta un diseo preliminar mediante
mtodos cortos, basados en las ecuaciones de
Fenske-Underwood-Gilliland-Kirkbride (FUGK). En
el siguiente paso, el sistema es validado empleando
simulacin rigurosa y se efectan cambios al diseo
de acuerdo al comportamiento deseado, optimizando
ya sea el consumo de energa en los rehervidores o el
costo total de la secuencia. Mtodos basados en el
enfoque jerrquico son los presentados por
Triantafyllou y Smith (1992), Annakkou y Mizsey
(1996), Hernndez (1998), Castro y Jimnez (2002),
y Muralikrishna y col. (2002). Por otro lado, en el
enfoque algortmico se desarrollan modelos de
optimizacin mixta entera no lineal con base en
superestructuras, en las cuales se incluyen todas las
alternativas de inters para el sistema (Dnnabier y
Pantelides, 1999; Yeomans y Grossman, 2000a;
2000b; Grossmann y col., 2005). Como resultado se
obtiene la mejor secuencia compleja, con el menor
costo, para una separacin especificada, que no es
necesariamente una columna Petlyuk. En este
enfoque se afronta la naturaleza discreta del
problema de sntesis y de diseo. Tiene la importante
caracterstica de involucrar modelos rigurosos pero

puede presentar dificultades computacionales
significativas, tales como la introduccin de
ecuaciones que llegan a ser singulares, la solucin de
muchas ecuaciones redundantes y el requerimiento
de excelentes estimaciones iniciales y lmites para las
variables para alcanzar convergencia. Su demanda
computacional puede ser excesiva. Grossmann y col.
(2005) proporcionan una revisin del estado del
enfoque, y concluyen que la combinacin de nuevas
representaciones de columnas individuales y de
superestructuras con tcnicas de programacin
matemtica disyuntiva y con esquemas robustos de
inicializacin basados en conceptos termodinmicos
hace posible resolver con eficiencia computacional
razonable estos problemas. En otro trabajo
relacionado, Shah y Kokossis (2002) proponen un
nuevo marco de referencia para la sntesis de
sistemas complejos de destilacin basado en el uso
de supertareas en lugar de superestructuras. En este
enfoque, el problema de optimizacin se formula
como un problema de programacin lineal mixta
entera que admite una solucin global ptima. Una
importante caracterstica del enfoque es que
aprovecha los esfuerzos desarrollados por otros
investigadores en la modelacin y la simulacin de
sistemas complejos de destilacin, tales como los
mtodos cortos, mtodos novedosos para clculo de
reflujo mnimo y anlisis de factibilidad.
Se han reportado tambin esfuerzos
recientes basados en clculos etapa a etapa.
Amminudin y col. (2001) desarrollan un complicado
mtodo de diseo semi-riguroso, que utiliza la
descomposicin en tres columnas convencionales
usada comnmente al aplicar mtodos cortos. Esta
representacin slo es vlida a reflujo mnimo
(Carlberg y Westerberg, 1989). Kim (2002) separa el
diseo estructural de la columna Petlyuk, que se
efecta primero a partir de clculos a reflujo total,
del diseo operativo, en el que se hallan el reflujo y
los flujos de interconexin requeridos; en ambos
casos, emplea clculos rigurosos de equilibrio. No
obstante, la estructura conseguida requiere tambin
correcciones.
Las correcciones que requieren los diseos
preliminares
basados
en
mtodos
cortos
normalmente se efectan sobre el reflujo y los flujos
de las corrientes de interconexin. En otros casos,
incluso debe modificarse el nmero de etapas en
algunas de las secciones de la columna Petlyuk
(Ramrez, 2003). La necesidad de efectuar estas
modificaciones puede atribuirse a las suposiciones de
derrame molar y volatilidades relativas constantes,
implcitas en el uso del mtodo FUGK. Amminudin
y col. (2001) mencionan tambin las limitaciones de
las correlaciones empricas del mtodo corto.
Las correcciones al reflujo y a los flujos de
interconexin puede efectuarse con relativa facilidad
debido a la naturaleza continua de esas variables. Sin
embargo, debido a la naturaleza discreta de las
122
variables involucradas, la correccin al nmero de

etapas en cada seccin slo puede efectuarse
mediante tediosa simulacin repetida cuando se
emplean simuladores comerciales (Kim, 2002;
Ramrez, 2002).
Otra alternativa para afrontar la naturaleza
discreta de los problemas de diseo de procesos de
separacin de etapas ha sido el uso de tcnicas de
colocacin. Cuando se aplican estas tcnicas, los
perfiles de las variables de proceso se aproximan
mediante polinomios, generalmente de bajo orden.
Las variables de proceso se tratan entonces como
funciones continuas de la posicin. Se consigue una
reduccin del tamao del problema, ya que un
polinomio de bajo orden puede aproximar los
perfiles en una columna de muchas etapas. Otra
consecuencia es que las variables discretas se pueden
tratar como variables continuas: el polinomio de
aproximacin permite conocer las variables an para
un nmero no entero de etapa. Esta caracterstica ha
permitido la aplicacin de tcnicas estndar de
programacin no lineal para resolver problemas de
diseo, donde deben definirse variables discretas
como el nmero de etapas y la localizacin de las
etapas de alimentacin o de productos para lograr
purezas o recuperaciones especificadas (Swartz y
Stewart, 1986; Seferlis y Hrymak, 1994). Debe
notarse que las tcnicas de colocacin no han sido
aplicadas a secuencias complejas de columnas de
destilacin, donde el nmero de variables discretas
es mayor.
En este trabajo, en el contexto del enfoque
jerrquico, se presenta una metodologa para
corregir, de una manera rpida y confiable, el diseo
preliminar de columnas Petlyuk. Se aplican tcnicas
de colocacin ortogonal y mtodos de programacin
no lineal estndar para modificar tanto la
localizacin de las etapas de interconexin superior e
inferior ( N L , NV ) y del producto lateral ( N S ) como
los flujos internos; con esto se busca obtener purezas
especificadas de los componentes en los productos
de domo, lateral y de fondo, respectivamente.
2. Simulacin mediante tcnicas de colocacin
Stewart y col. (1985) y Seferlis y Hrymak

(1994) han explicado con suficiente detalle la
aplicacin de tcnicas de colocacin para simular
columnas de destilacin convencionales. En esta
seccin se muestran slo las ideas centrales del
enfoque, considerando una cascada de N E etapas
que separa una mezcla de N C componentes. Una
etapa de equilibrio representativa de la cascada,
denotada s j , se muestra en la Fig. 2.
La mezcla a separar presenta volatilidades
relativas constantes y se supone que los flujos
internos tambin lo son. Bajo estas condiciones, el
modelo matemtico para la etapa representativa se
obtiene al combinar el balance molar para el
componente i en trminos de flujos molares de los

componentes:
f i ( s j ) + vi ( s j + 1) + li ( s j 1)
U (s j )
W (s j )
(1)
1 +
li ( s j ) 1 +
vi ( s j ) = 0
L( s j )
V ( s j )
con la relacin de equilibrio lquido-vapor para
sistemas de volatilidad relativa constante:
i li ( s j )
(2)
vi ( s j ) = NC
V (s j )
k lk ( s j )
k =1
Fig. 2 Etapa representativa

El modelo habitual del sistema, llamado
modelo de orden completo por Stewart y col. (1985),
se obtiene aplicando estas ecuaciones al total de
etapas discretas existentes, es decir haciendo
s j = j , j = 1, , N E . El modelo es entonces un
conjunto de N E N C ecuaciones, cuya solucin
proporciona los N E N C flujos molares desconocidos
de los componentes en la fase lquida, li , j .
Considrese ahora el modelo basado en colocacin,
tambin llamado modelo de orden reducido. Suponga
que los perfiles de flujos molares se representan
como funcin del nmero de etapa mediante un
polinomio de grado N PC 1 . Ya que el polinomio
contiene N PC coeficientes desconocidos slo se
requiere, en principio, aplicar el modelo del sistema
en N PC etapas o puntos de colocacin de la cascada.
En este caso, por lo tanto, j = 1, , N PC . En
principio, la ubicacin s j de los puntos de
colocacin puede ser arbitraria. Sin embargo,
estudiando columnas convencionales, Stewart y col.
(1985) mostraron la conveniencia de localizar dichos
puntos como las races de los polinomios ortogonales
discretos de Hahn de grado N PC 1 . En suma, el
modelo de colocacin es un conjunto de solamente
N PC N C ecuaciones. Su solucin proporciona los
N PC N C flujos molares desconocidas en los puntos
de colocacin. Los flujos molares en cualquier otra
etapa real o ficticia, s j [1, N E ] , se obtienen durante
y al final del procedimiento mediante interpolacin

de Lagrange a partir de las flujos molares en los
puntos de colocacin.
123
Con frecuencia es necesario considerar varios

elementos finitos, es decir, ms de un polinomio,
para representar los perfiles en una cascada o
seccin. Si este es el caso, al modelo matemtico
anterior deben aadirse ecuaciones de continuidad de
flujos en las fronteras de los elementos finitos
(Seferlis y Hrymak, 1994). Detalles adicionales
acerca de la implementacin de mtodos de
colocacin pueden encontrarse en las referencias
citadas.
Generalmente estas aproximaciones son adecuadas.

Para aplicar el mtodo de colocacin se defini
generalmente un elemento finito por cada 7 a 15
etapas discretas y se utilizaron cuatro puntos de
colocacin por elemento. Adems, para evitar la
adicin las ecuaciones de continuidad de los perfiles,
se localiz un punto de colocacin en cada frontera
de los elementos finitos.
DISEO
PRELIMINAR
3. Diseo refinado de columnas Petlyuk
La Fig. 3 muestra la metodologa empleada

para corregir los diseos preliminares. La columna
Petlyuk se dise mediante el mtodo corto de
Castro y Jimnez (2002), utilizando un programa de
cmputo de estos investigadores. Se especifican la
corriente de alimentacin, la presin de operacin,
las purezas requeridas, la relacin de reflujo
(respecto del reflujo mnimo) en ambas columnas y
las recuperaciones de los componentes ligero y
pesado en el domo de la columna prefraccionadora.
El programa proporciona las variables de diseo
D, S , R, LS , VI , N1 , N 2 , N L , N S , NV , N E . El valor de
un subconjunto X de estas variables de diseo se
considera sujeto a correccin, pero dejando constante
el total de etapas del diseo preliminar, tanto en la
columna prefraccionadora como en la columna
principal. El diseo se simula para conocer las
purezas de los componentes A, B y C en destilado,
producto lateral y residuo. Si la funcin objetivo
f ( X ) = ( x A, D x Aesp, D ) 2 + ( xB , S xBesp, S ) 2
(3)
+ ( xC , B xCesp, B ) 2
ha alcanzado su valor mnimo, el diseo se considera
satisfactorio y puede validarse mediante simulacin
rigurosa. En caso contrario, el subconjunto de
variables X se actualiza en el modulo de
optimizacin, que emplea el mtodo Simplex de
Nelder y Mead (Edgar y Himmelblau, 1988). El
diseo refinado que resulta se valid mediante el
simulador Aspen-Plus 10.2, empleando el mtodo de
propiedades IDEAL. En el modulo de simulacin se
suponen volatilidades relativas y derrame molar
constantes, es decir, es congruente con el mtodo
FUGK, pero utiliza tcnicas de colocacin para
permitir valores continuos de N L , N S , NV .
Para ilustrar la aplicacin del mtodo se
reportan los dos casos descritos en las Tablas 1 y 2,
aunque se analizaron todos los casos estudiados por
Jimnez y col. (2003). El flujo de alimentacin fue
de 100 kmol/h y la presin de diseo se fij para
asegurar el empleo de agua como medio de
enfriamiento. Las purezas de diseo se indican en
negritas en las tablas. Para completar la distribucin
de componentes en los productos se suponen
separaciones ntidas entre B y C en domo, y entre A y
B en fondo. Asimismo, se supone separacin
simtrica entre A y C en el producto lateral.
L S VI R D S
Ni
X
X
OPTIMIZACIN
SIMULACIN
X
NO
f(X)Min ?
SI
X
Fig. 3 Diagrama de flujo para obtener diseos

refinados.
Tabla 1 Caso de estudio I
COMPONENTE i
zi
xi , D xi , S
i-pentano
n-pentano
n-hexano
3.04
2.44
1.00
0.33 0.98 0.02 0.00

0.33 0.02 0.96 0.02
0.34 0.00 0.02 0.98
Tabla 2 Caso de estudio II

COMPONENTE i
zi
xi , D
xi , S
n-pentano
n-hexano
n-heptano
5.06
2.20
1.00
xi , B
0.15
0.70
0.15
0.91
0.09
0.00
0.005
0.990
0.005
xi , B
0.00
0.09
0.91
4. Resultados y discusin
Los resultados para el caso I se proporcionan

en la Tabla 3. Para este caso el diseo preliminar
indic que en la columna principal N L = 39 ,
N S = 72 y NV = 77 . La simulacin de la secuencia
con esta configuracin indic que los productos
lateral y de fondo no alcanzaban las purezas
requeridas (rengln 1). La aplicacin de la
metodologa de refinamiento del diseo, con como
variables
independientes
de
optimizacin,
N L = 30.89 ,
N S = 69.90
y
proporcion
NV = 76.98 . Al simular el diseo final, mediante el
modelo de orden completo, se encontr que se
lograban prcticamente las purezas de diseo
(rengln 2). La validacin del diseo obtenido indic
que se satisfacen las especificaciones (rengln 3). En
este caso de estudio debi modificarse tanto la
124
localizacin de la interconexin superior como la del

producto lateral sin modificar ningn flujo.
En la Tabla 4 se muestra que el diseo
preliminar no es satisfactorio para el caso de estudio
II: en el rengln 1 las purezas de los componentes,
obtenidas con el reflujo del diseo preliminar
R = 6.3 , estn lejos de la especificacin. Para este
caso, la correccin al conjunto de variables
independientes X = ( N L , N S , NV ) no fue suficiente
para lograr un diseo refinado adecuado. Se decidi,
por lo tanto, probar el conjunto de variables
independientes
de
optimizacin
X = ( N L , N S , NV , LS , VI ) a distintos valores de
reflujo en la columna principal. El diseo refinado
obtenido con R = 14 se da en el rengln 2 y su
simulacin indic la satisfaccin de las
especificaciones. La validacin mediante el
simulador comercial indic que es necesario un
ligero incremento adicional al reflujo para llevar a
especificacin la pureza del componente pesado
(rengln 3). Finalmente, debe resaltarse el
incremento tan importante que debi hacerse al
reflujo preliminar, de 6.3 a 14.0, para satisfacer las
especificaciones.
En ambos casos de estudio existen al menos
3 variables discretas por determinar, por lo que
deben examinarse un gran nmero de combinaciones
mediante simulacin repetida cuando el diseo se
refina usando simuladores comerciales. En contraste,
el diseo refinado se hall rpidamente empleando la
metodologa indicada. Por otra parte, es importante
enfatizar que el mdulo de simulacin utiliza un
modelo basado en volatilidades y derrame molar
constantes, implcitos tambin en el uso del mtodo
FUGK. La validacin de los diseos refinados indica
que ese modelo es efectivo para las mezclas de
hidrocarburos estudiadas. Por lo tanto, la necesidad
de correcciones a los diseos preliminares obtenidos
en el enfoque jerrquico se debe atribuir ms a las
limitaciones de las correlaciones empricas (Gilliland
o Kirkbride) que a las suposiciones de derrame molar
y volatilidades constantes, consistente con lo
DISEO
indicado por Amminudin y col. (2001). Actualmente

se desarrollan trabajos para incorporar modelos
semi-rigurosos o rigurosos tanto en el mdulo de
diseo como en el simulacin-optimizacin, a fin de
superar las limitaciones del enfoque jerrquico y de
considerar separaciones de mezclas no ideales.
Conclusiones
Se ha presentado una metodologa para

mejorar el diseo preliminar de columnas Petlyuk
generado mediante mtodos cortos. La aplicacin de
la metodologa permite hallar la distribucin de
etapas y/o los flujos de interconexin que
proporcionan purezas especificadas en sus productos,
aplicando tcnicas de colocacin ortogonal a un
modelo basado en volatilidades y flujos internos
constantes. Una caracterstica relevante del mtodo
es la posibilidad de usar tcnicas estndar de
optimizacin no lineal, tal como el mtodo Simplex
de Nelder y Mead. Dada la simplicidad de los
mtodos cortos y las caractersticas del mtodo de
colocacin presentado, la metodologa es
globalmente sencilla y, para mezclas ideales,
compite favorablemente con mtodos que utilizan
clculos etapa a etapa, como los de Kim(2002), y
Amminudin y col. (2001). Se confirm que la
necesidad de correcciones al diseo preliminar se
debe ms a las limitaciones de las correlaciones
empricas del mtodo FUGK que a las suposiciones
de flujos internos y volatilidades constantes. El
mtodo desarrollado disminuye notablemente el
esfuerzo para conseguir diseos mejorados, en
comparacin con el necesario cuando se utiliza
simulacin repetida mediante paquetes comerciales.
Modificando la representacin de las propiedades
termodinmicas, la metodologa puede extenderse
para efectuar estudios de diseo de separaciones
azeotrpicas y de destilacin reactiva.
Tabla 3 Resultados para caso de estudio I

SIMULACIN
NL
NS
NV
x A, D = 0.98 xB , S = 0.96 xC , B = 0.98
Preliminar 39
72
77
Refinado 30.89 69.90 76.98
Validacin 31
70
77
0.98
0.98
0.98
0.94
0.96
0.97
0.96
0.97
0.99
Tabla 4 Resultados para caso de estudio II

SIMULACIN
DISEO
NL
NS
NV
x A, D = 0.91 xB , S = 0.99 xC , B = 0.91
LS
VI
F
F
Preliminar 10
18
23 0.3885 0.7429
0.88
0.92
0.65
Refinado 10.94 15.94 23.44 0.5025 1.1053
0.92
0.99
0.91
Validacin 11
16
24 0.5025 1.1053
0.92
0.99
0.90
125
Superndices
esp
especificacin
Agradecimientos
Los autores agradecen el financiamiento a este
trabajo otorgado por CONACYT e IMP.
Nomenclatura
B
D
fi
li
L
LS
NC
NE
NL
N PC
NS
NV
R
sj
S
U
vi
V
VI
x
X
y
W
flujo molar de residuo

flujo molar de destilado
flujo molar del componente i en la
alimentacin
flujo molar del componente i en la fase
lquida
flujo molar de lquido que abandona la
etapa
flujo molar de lquido de interconexin
superior
nmero de componentes
nmero de etapas
etapa de interconexin superior en la
columna principal
nmero de puntos de colocacin
etapa de producto lateral en la columna
principal
etapa de interconexin inferior en la
columna principal
relacin de reflujo en la columna principal
etapa genrica
flujo molar de producto lateral
flujo molar de salida lateral de lquido
flujo molar del componente i en la fase
vapor
flujo molar de vapor que abandona una
etapa
flujo molar de vapor de interconexin
inferior
fraccin molar en fase lquida
vector de variables independientes del
diseo refinado
fraccin molar en fase vapor
flujo molar de salida lateral de vapor
Smbolos griegos
volatilidad relativa
Subndices
A
componente ligero
B
componente de volatilidad
residuo
C
componente pesado
D
destilado
i
componente
j
etapa o punto de colocacin
k
seccin de la columna Petlyuk
S
producto lateral
intermedia,
Referencias
Amminudin, K.A., Smith, R., Thong, D.Y.C. y
Towler, G.P. (2001). Design and Optimization
of Fully Thermally Coupled Distillation
Columns. Part 1: Preliminary Design and
Optimization Methodology. Transactions of
the Institution of Chemical Engineers 79, Part
A, 701-715.
Annakou, O. y Mizsey, P. (1996). Rigorous
Comparative Study of Energy-Integrated
Distillation Schemes. Industrial Engineering
Carlberg, N. A. y Westerberg, A. W. (1989).
Temperature-Heat Diagrams for Complex
Columns. 3. Underwoods Method for the
Petlyuk Configuration. Industrial Engineering
Castro-Agero, A. y Jimnez-Gutirrez, A. (2002).
Mtodo Corto para el Diseo de Columnas de
Destilacin Trmicamente Acoplada Petlyuk.
Revista Mexicana de Ingeniera Qumica 1,
85-96.
Dunnebier, G. y Pantelides, C. (1999). Optimal
Design of Thermally Coupled Distillation
Columns. Industrial Engineering Chemistry
Research 38, 162-176.
Edgar, T.F. y Himmelblau, D.M. (1988).
Optimization
of
Chemical
Processes.
McGraw-Hill, New York.
Finn, A. (1993). Consider Thermally Coupled
Distillations. Chemical Engineering Progress,
October, 41-45.
Glinos, K. y Malone, M. F. (1988). Optimality
Regions for Complex Column Alternatives in
Distillation Systems. Chemical Engineering
Research & Design 66, Part A, 229-240.
Grossmann, I.E., Aguirre, P. y Barttfeld, M. (2005).
Optimal Synthesis of Complex Distillation
Columns Using Rigorous Models. Computers
& Chemical Engineering 29, 1203-1215.
Hernndez, S. (1998). Estudio sobre la Optimizacin
y el Control de Esquemas de Destilacin
Trmicamente
Acoplados.
Tesis
de
Doctorado, Instituto Tecnolgico de Celaya,
Mxico.
Optimal
Thermally-Coupled
Distillation
System
Using
a
Dynamic
Model.
Transactions of the Institution of Chemical
Engineers 74, Part A, 357-362.
Energy-Efficient Petlyuk Systems. Computers
& Chemical Engineering 23, 1005-1010.
Jimnez, A., Ramrez, N., Castro, A. y Hernndez, S.
(2003). Design and Energy Performance of
Alternative Schemes to the Petlyuk
126
Distillation System. Transactions of the

Institution of Chemical Engineers 81, Part A,
518-524.
Kim, Y.H. (2002). Structural Design and Operation
of a Fully Thermally Coupled Distillation
Column. The Chemical Engineering Journal
85, 289-301.
Muralikrishna, K., Madhavan, K.P. y Shah, S.S.
(2002). Development of Dividing Wall
Distillation Column Design Space for a
Specified Separation. Transactions of the
155-166.
Ramrez Corona, N. (2002). Estudio Comparativo de
Secuencias de Destilacin no Convencionales
para Mezclas Ternarias. Tesis de Maestra,
Instituto Tecnolgico de Celaya, Mxico.
Ramrez Corona, N. (2003). Comunicacin personal.
Seader, J. D. y Henley, E. J. (1998). Separation
Process Principles. Wiley, New York.
Seferlis, P. y Hrymak, A.N. (1994). Optimization of
Distillation Units Using Collocation Models.
AIChE Journal 40, 813-825.
Shah, P.B.y Kokossis, A.C. (2002). New synthesis
framework for the optimization of complex
distillation systems. AIChE Journal 48, 527550.

Stewart, W.E., Levien, K.L. y Morari, M. (1985).
Simulation of Fractionation by Orthogonal
Collocation. Chemical Engineering Science
40, 409-421.
Swartz, C.L.E. y Stewart, W.E. (1986). A
Collocation Approach to Distillation Column
Design. AIChE Journal 32, 1832-1838.
Triantafyllou, C. y Smith, R. (1992). The Design and
Optimisation of Fully Thermally Coupled
Distillation Columns. Transactions of the
118-132.
Yeomans, H. y Grossmann, I.E. (2000a). Optimal
Design of Complex Distillation Columns
Using Rigorous Tray-by-tray Disjunctive
Programming Models. Industrial Engineering
Yeomans, H. y Grossmann, I. E. (2000b). Disjuntive
Programming Models for the Optimal Design
of Distillation Columns and Separation
Sequences. Industrial and Engineering
127
40
Revista Ingeniera e Investigacin No. 48 Diciembre de 2001
Mdulo Thtorial de Simulacin para el Clculo

de Torres de Destilacin por Etapas
Hermes A. Rangel Jara= , Gustavo Alberto Rodrguez Gmez**, Juan Carlos Rosas Molina**
REsUMEN
En este artculo se presenta un software desarrollado para

la enseanza del proceso de destilacin por etapas, enfocado
para ser utilizado por los estudiantes de pregrado de Ingeniera
Qumica. Se expone en l un marco general de los procesos
de separacin,
haciendo nfasis en el proceso de la
destilacin, sus orgenes, su importancia y sus perspectivas.
Posteriormente se desarrolla el tema del equilibrio de fases,
que,junto con el concepto de etapa de equilibrio, fundamentan
el clculo de la destilacin por etapas. En este ltimo aspecto
se realiza una explicacin de los modelos matemticos ms
representativos para las mezclas binarias, de una manera
relativamente
completa, mientras que se establece una
aproximacin al estudio de algunos modelos para mezclas
multicomponentes.
Para esto se implementan programas
tutores, que presentan paso a paso el procedimiento para
resolver un problema planteado por el usuario, a la vez que se
obtienen los resultados. Adems, con una intencin centrada
en solucionar un problema sin revisar el procedimiento de
clculo, se adicionan programas
simuladores
para los
diferentes modelos estudiados.
INTRODUCCIN
adestilacin es el procedimiento de separacin de especies

qumicas ms extensamente utilizado en las industrias de
proceso. Pese a ser una tcnica muy antigua, contina a la
vanguardia gracias a que posee una base tecnolgica relativamente
madura, siendo el mtodo de separacin clave contra el cual
otros mtodos alternativos pueden compararse [2].
Desde sus inicios, la ingeniera qumica ha estado ligada al
desarrollo de la destilacin, y son sus profesionales los que tienen
competencia ms directa con ella. Lo anterior pone de manifiesto
la necesidad de una adecuada formacin de los ingenieros
qumicos en esta importante tcnica. El trabajo que aqu se
presenta es un intento por proveer nuevas herramientas que
complementen la formacin acadmica tradicional en esta rea.
La preparacin del ingeniero qumico no debe limitarse al
conocimiento de los aspectos conceptuales de la operacin,
sino que debe complementarse con modelos matemticos que
sean representativos de las situaciones reales. En este sentido,
la simulacin de procesos brinda algunas herramientas que
permiten solucionar eficaz y eficientemente dichos modelos.
ABSTRACT
This article presents a computer program developed for

teaching of the stagewise distillation process, to be used by
undergraduate
chemical engineering students.Tt exposes a
separation processes overview, emphasizing in distillation,
its origins, importance and prospects. Phase equilibrium is
then studied, as well as the equilibrium stage concept, both
fundamental for stagewise distillation calculations. In the
latter aspect an explanation of representative mathematic
models is presented in a relatively complete way for binary
mixtures, whereas an approximation is made to the study
of multicomponent
mixture models. Tutor programs are
implemented presenting step-by-step procedures to solve
user specified problems at the time results are obtained.
Moreover, with a problem solving focused intention,
simulation programs are added for the different studied
models.
Sin lugar a dudas, la aparicin del computador y el desarrollo

de la informtica han influido de manera importante sobre la
sociedad de fin de siglo en diversas reas. En el campo de la
educacin, el computador se ha convertido en una herramienta
indispensable por diferentes razones: permite el manejo fcil y
eficiente de la informacin, brinda la capacidad de desarrollar
trabajos de diversa ndole hbilmente, permite mayor interaccin
entre el usuario y la informacin, etc. En la simulacin de
procesos, el computador tiene la mayor relevancia, pues es el
medio para poder implementar los modelos matemticos y hallar
soluciones satisfactoriamente. Lo anterior hace que sea el
instrumento esencial en el desarrollo de este trabajo, que en
esencia es un software, que adems de ser un simulador para el
clculo de torres de destilacin por etapas, es un tutor de la
operacin.
Ingeniero Qumico, M.Sc.I.Q., profesor titular, Universidad

.. Ingeniero Qumico, Universidad Nacional de Colombia.
Nacional de Colombia .
Mdulo Tutorial de Simulacin para el Clculo de Torres de Destilacin por Etapas
Algunos programas de esta ndole estn siendo desarrollados

en diferentes universidades de otros pases, y es la intencin
de este trabajo introducir la enseanza de la Ingeniera Qumica
de la Universidad Nacional en esta tendencia global de crear
alternativas
pedaggicas,
que utilicen las bondades
tecnolgicas y la gran versatilidad que representa la tecnologa
computacional. Adems, se busca promover la creacin de
trabajos orientados en la misma direccin en diferentes reas
de la ingeniera qumica, con la intencin de poder crear un
software acadmico integral.
l. LA GUA
TUTORlAL DE DESTILACIN
En el software desarrollado se presenta una recopilacin

terica de la destilacin, que recoge los aspectos conceptuales
necesarios para una razonable cobertura de las nociones bsicas
que permiten el entendimiento de la operacin, as como el
estudio de la misma.
41
modelar, descrito por los comportamientos ideal y no ideal,

para poder abocar finalmente el clculo de los parmetros
termodinmicos de un sistema en equilibrio. En este sentido, se
trata una ecuacin de estado cbica (Peng - Robinson) para el
clculo del equilibrio lquido-vapor y los trminos de energa,
tanto para sistemas binarios como para multicomponentes.
Dicha ecuacin trabaja con una base de datos de 50 hidrocarburos.
Adems, se incluye un modelo de clculos de actividad (Unifac)
para describir la fase lquida, que utiliza 75 compuestos como
base de datos.
En la segunda seccin se presentan los conocimientos
elementales de los procesos de separacin por etapas de
equilibrio, los diferentes tipos de configuraciones que pueden
presentar, los balances de materia y energa que se aplican y.
finalmente la representacin grfica que se puede hacer de un
proceso de separacin de etapa mltiple.
En un aparte adicional se hace una breve descripcin de la
estructura de los equipos utilizados para el contacto entre fases
en la destilacin y las unidades de transferencia de calor que
requieren. Adems, se tratan los aspectos ms notables de su
operacin.
Figura l. Pantalla de inicio de la Gua tutorial de destilacin.
Se expone primero un marco general de los procesos de

separacin, destacando la vital importancia que tienen en los
procesos industriales, y luego ejemplos de su aplicacin.
Despus de esta observacin global, se realiza un acercamiento
a los aspectos bsicos de la destilacin, incluyendo su
definicin, su desarrollo histrico y su relevancia presente y
futura. Aqu, se le compara con las tcnicas de separacin
alternativas y se hace un compendio de los principales campos
de estudio referentes a ella que estn desarrollndose en la
actualidad.
Posteriormente, se estudian los aspectos que dan fundamento
a la operacin: el equilibrio de fases y las operaciones de etapa
de equilibrio.
En la primera seccin se tratan inicialmente los conceptos
bsicos termodinmicos, dentro de los cuales est el estudio
del equilibrio lquido-vapor. Se sigue con el anlisis del estado
del sistema, el tratamiento matemtico con el cual se puede
Finalmente, tratados los aspectos fundamentales, se afronta

el estudio de los esquemas de clculo ms representativos de
cada uno de los tipos de operacin que se pueden presentar en
la destilacin por etapas y su descripcin. Primero se abarcan
los procesos que operan de manera discontinua: destilacin
diferencial y por lotes (con reflujo variable y constante), ambos
para mezclas binarias. Luego se trata la destilacin instantnea
oflash (isotrmica, adiabtica y con fraccin de vaporizacin)
para sistemas multicompuestos. Posteriormente se trabaja la
rectificacin continua de mezclas binarias, en donde se tratan
los conocidos e ilustrativos mtodos de Ponchon-Savarit, Me
Cabe- Thiele y Sorel. Finalmente se toca la rectificacin continua
de mezclas multicomponentes
en donde se ocupa primero
de los mtodos cortos para la determinacin
de la
distribucin
de los componentes,
la determinacin
del
nmero mnimo de etapas, el clculo del reflujo mnimo, el
clculo del nmero de etapas ideales y la localizacin del plato
de alimentacin. Finalmente se hace una aproximacin inicial a
dos mtodos rigurosos: el de Wang-Henke y el mtodo theta
de convergencia.
Dada la densidad de la informacin, se busc que la manera
en que se presentara fuera la ms amigable y fcilmente
accesible para el usuario. Para dirimir esta necesidad, se
utilizaron la versatilidad y el dinamismo que caracterizan a las
pginas Web de Internet, as que se elabor un total de 83
pginas para cobijar todo el marco terico (el programa con el
cual fueron realizadas fue Microsoft Frontl'age=). Este tipo de
formato permite navegar a travs de diferentes hojas por medio
de hipervnculos
conectando
trminos relacionados
o
referencias; tiene la capacidad de utilizar imgenes animadas,
facilita el acceso a informacin adicional en pantalla en forma
de cuadros de texto, etc., acciones que pueden realizarse con
una simple pulsacin del ratn.
42
Necesariamente
todas las hojas deban estar bajo una
estructura que las aunara. Entonces, se utiliz un programa
diseado para elaborar archivos de ayuda, llamado HTML Help
Workshop=, que permite la compilacin de hojas Web. Este
programa habilita la organizacin del contenido en forma
jerrquica en una estructura en forma de rbol o tabla de
contenido, en la cual se puede encontrar cualquiera de las hojas
del cuerpo del software para que sea visualizada (tal como se
presenta en la figura 2). Adems, posee herramientas de
bsqueda de palabras clave o de trminos aislados.
En la figura 3 se presentan irnoenes correspondientes a

algunas de las presentaciones explicativas incluidas en el
programa.
SI.. Ir .., pIIo1o 1'1'*10 C1JIIqJIor coro del polIocto, por eJon1>lo 111
pIIo1o
lien1>Orlh
- ~,
11propIedod exc~Ido, .......
CMO 11 presin, _i
c"' Ir'" de tDde IIIl4*fIde
del pIIno Irtado, Y al volardepondIri de
dnde corbI el pIIo1o su eje. Curdo pIonoI 1'1'_
111 ...
de ....,......
o
presin .....
conc:aWlCIn, .......
Ir'" del -'"
del 0CJJIbrI0 lo CIIftIn, cIIndo
~
de ..,.,,10
vopcr -1fqJIdo presin o lien1>Or...... -1196>'" el CMO. _10000._, CIClo ~ de _ ..... pueden _o ... _oodoI
.... .., pIIo1o bIdrnenIIonII.
~
.. "".,..,dtllhlRrlL;..
.... *"hoITO' ~~~tn
... aqcdDIdlbln_
.....p-*".IAi~~dI
__ ".)7......
.........
lC .....
Er ...
qL.e~~n'!lJCh::lporEXr'(lQlr'*-qAtH
tlglll'lGdllhotrtn.
en
t. d6rwtw ....... oprt:lCa)l. Ui
""''I
nohliHl) ....
1M
,.,.)7.
o.dI DemI:!trliDt.c.II~
tldlllhllUIIClInIIt"adllll'unrlldld_",*,port.N: ..........
.t,I"O.p.- ...
W1onW:lfttllnt)dllot"O.lil__...
dt.,.... NlDkDI q..-I'W\ ~W"
.-....-..:
.. Un 1I..1da
dIf~.,
.....
.,"~c..-..:ItG.o.ntodl_~IIn.~
,......t....-IfA~~
... QIn .. dIIIrroIIDdI.dllcMc:lI!n
l1li""'''''
a)
Figura 2. Tabla de contenido desplegada.
Como resultado
una gran variedad
de las herramientas
fue realizado. Este
destilacin.
se tiene un libro electrnico, que permite

de acciones con mucha facilidad en virtud
que exhiben los programas con los cuales
libro recibe el nombre de Gua tutorial de
Por otra parte, desde el libro electrnico es posible llamar

otros programas por medio de botones que son programados
para que ejecuten una aplicacin especfica. Esta habilidad se
utiliza en el programa para ejecutar las presentaciones
explicativas, los tutores y los simuladores, los cuales se detallan
a continuacin.
b)
Figura 3. Presentaciones explicativas. a) Equilibrio
b) Mtodo de Ponchon Savarit.
m.
11. LAS
PRESENTACIONES
Los
lquido-vapor
TUTORES
EXPLICATIVAS
Para algunos temas especficos se desarrollan presentaciones

que amplan la explicacin del tema en cuestin. stas son
exposiciones
con alto contenido de elementos visuales
animados que avanzan automticamente mientras presentan la
informacin pertinente. Dichas presentaciones
han sido
realizadas en el programa Microsoft Powerf'oint=, y se
encuentran identificadas en la gua tutorial por el siguiente
icono:
La presentacin arrancar simplemente pulsando sobre el

botn que se encuentra inmediatamente debajo de la figura.
Constituyen una parte fundamental del mdulo tutorial, ya

que permiten realizar clculos mientras se aprende el
procedimiento paso a paso. Fueron desarrollados sobre la hoja
electrnica Microsoft Excel".
Los programas tutores tienen como una de sus principales
ventajas la capacidad de desarrollar problemas diversos, debido
a que poseen una suficiente base de datos de compuestos,
pueden modificarse las variables de operacin (temperatura,
presin, composicin, etc.) y seleccionarse entre diferentes tipos
de clculo. Por medio de botones se avanza o se retrocede en
el clculo apareciendo los resultados pertinentes sobre la hoja
electrnica; tambin es posible reiniciar el clculo o abandonar
el ejemplo en cualquier momento. Los tutores se identifican
dentro de la gua por la siguiente imagen, que se encuentra
acompaada de un botn para ejecutar la aplicacin:
43
La figura 4 presenta algunas imgenes del tutor del mtodo

de Ponchon-Savarit.
d)
Figura 4. Imgenes del tutor del mtodo Ponchon-Savarit.
a)
Bienvenida. b) Seleccin de compuestos. e) Introduccin de datos.
d) Pasos del clculo.
IV. Los
a)
SIMULADORES
Estos programas tambin estn estructurados sobre la hoja

electrnica y se diferencian de los tutores en que no presentan
el procedimiento de clculo, sino que de la introduccin de
datos pasan directamente a la presentacin de resultados. La
intencin de estos programas es permitir el anlisis de los
resultados, as como observar la incidencia de las diferentes
variables sobre ellos.
En estos programas no se observa ninguna modificacin
sobre la hoja de Excel= ya que el clculo se realiza de manera
oculta. Se distinguen de las dems aplicaciones por el siguiente
icono:
b)
c)
La figura 5 presenta algunas imgenes del simulador de

destilacin instantnea (flash).
a)
44
En el siguiente cuadro se hace una relacin de los programas

tutores y simuladores incluidos en la gua:
Programa 1. Tutor del clculo de equilibrio de mezclas binarias segn la Ley
de Raoult.
Programa 2. Tutor del clculo de temperaturas de burbuja y roco de mezclas
multicomponentes utlizando la Lev de Raoult.
Programa 3. Tutor del clculo de la temperatura de burbuja utilizando la
ecuacin de estado de Peng - Robinson.
Programa 4. Tutor del clculo de la temperatura de roco utilizando la
ecuacin de estado de Pen - Robinson.
Programa 5. Tutor del clculo de la presin de burbuja utilizando la ecuacin
de estado de Pena - Robinson.
Programa 6. Tutor del clculo de la presin de roco utilizando la ecuacin de
estado de Peng - Robinson.
Programa 7. Simulador para el clculo de los puntos de burbuja y roco
utilizando la ecuacin de estado de Pen - Robinson.
Programa 8. Tutor para el clculo de los coeficientes de actividad por medio
del modelo UNIFAC.
Programa 9. Simulador para el clculo de los coeficientes de actividad por
medio del modelo UNIFAC.
Programa 10. Simulador para el clculo de los puntos de burbuja y roco
utilizando la ecuacin de estado de Peng - Robinson y el modelo de
actividad (UNIFAC).
Programa 11. Tutor para el clculo qrfico del nmero de etapas ideales.
Programa 12. Tutor para el clculo analtico del nmero de etapas ideales
I por el mtodo de Kremser.
Proarama 13. Tutor cara el clculo de la destilacin diferencial.
Programa 14. Tutor para el clculo de la destilacin discontinua con relacin
de reflujo_constante.
Programa 15. Tutor para el clculo de la destilacin discontinua con relacin
de reflujo variable.
Programa 16. Simulador para el clculo de la destilacin discontinua.
Programa 17. Tutor para el clculo de la destilacin instantnea isotrmica
con coeficientes de distribucin independientes de la composicin.
Programa 18. Simulador para el clculo de la destilacin instantnea
isotrmica con coeficientes de distribucin calculados por la ecuacin de
estado de Pena - Robinson.
Programa 19. Tutor para el clculo de la destilacin instantnea con fraccin
de vaporizacin con coeficientes de distribucin independientes de la
composicin.
Programa 20. Simulador para el clculo de la destilacin instantnea con
coeficientes de distribucin independientes de la composicin.
Programa 21. Tutor para el clculo de la destilacin continua de mezclas
binarias por medio del mtodo de Ponchon - Savarit.
binarias por medio del mtodo de McCabe - Thiele.
binarias por medio del mtodo de Sorel.
Programa 24. Tutor para el clculo de la distribucin de los componentes en
la destilacin continua de mezclas multicomponentes.
Programa 25. Tutor para el clculo del nmero mnimo de etapas por el
b)
Cc)
CONCLUSIONES
~---
......................
----,.
--"";"'WI~~x:::;::
_:---w.~c:n::.::D:r..:::=
El trabajo desarrollado puede ser de gran ayuda para

complementar la formacin conceptual de los estudiantes de
pregrado de Ingeniera Qumica en el rea de la destilacin.
Esto se debe a que ofrece variedad de opciones para el
acercamiento
a un tema determinado,
como acceso a
informacin terica, visualizacin de presentaciones explicativas
y ejecucin de tutores y de simuladores; todo organizado bajo
una misma estructura dinmica que permite la eleccin del tipo
de herramienta que se quiera utilizar .
r-.~----.c::;;:;:::i:.::!:i:=:~
.:-~~~-'~
.-__.
t:.:::.::
uc.:: ~
El beneficio que puede brindar este trabajo dentro de la

formacin acadmica de los ingenieros qumicos depende en
gran parte de que sea divulgado a travs de los cursos o
actividades en donde se realice el estudio de la destilacin o de
la simulacin de procesos.
.?"--'~~.!..:...::.
d)
Figura
5. Imgenes
problema.
del simulador
b) Introduccin
de destilacin
de datos.
flash.
e) Resultados.
a) Seleccin
d) Otros
del
resultados.
Contina latente la necesidad de realizar proyectos de este

tipo en diferentes campos de la ingeniera qumica, para avanzar
hacia el objetivo global de constituir una herramienta integral
de simulacin con un fin estrictamente didctico.

BIBLIOGRAFA
1. HUMPHREY, J. "Separation Processes: Playing a critical
Engineering Progress, October, 1995, Vol. 91(10): 31-41.
2. KUNESH, J. el al., "Distillation:
Engineering Progress, October,
role",
En: Chemical
Still towering over other options",

1995, Vol. 91(10): 43-54.
En: Chemical
3. RODRGUEZ, G. A. Y ROSAS 1. C .. "Mdulo Tutorial de Simulacin para el Clculo

de Torres de Destilacin por Etapas". Proyecto de grado presentado a la Universidad
Nacional de Colombia para optar al grado de Ingeniero Quntico, Santaf de Bogot,
D.C.,1999.
4. WEXLER, S., The official Microsoft HTML help authoring kit: Understanding,
crealing and migraling to Microsoft HTMLHelp.
United States of America, Microsoft
Press, 1998.
45

Articulos Mezclas

Cargado por

Copyright:

Formatos disponibles

Articulos Mezclas

Cargado por

Información del documento

Título 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

Articulos Mezclas

Cargado por

Copyright:

Formatos disponibles

Red de Revistas Cientficas de Amrica Latina, el Caribe, Espaa y Portugal

Sistema de Informacin Cientfica

Francisco Snchez Mares, Pedro Alonso Dvila

Ms informacin del artculo

Conciencia Tecnolgica No. 34, Julio - Diciembre 2007

Prediccin de Puntos Crticos en Mezclas de Gas Natural

M. C. Francisco Snchez Mares y Dr. Pedro Alonso Dvila

Actualmente es el combustible ms limpio para la

La prediccin confiable de puntos crticos en mezclas

PREDICCIN DE PUNTOS CRTICOS EN MEZCLAS DE GAS NATURAL

de estado cbicas. La convergencia al punto crtico se

Solucin de las condiciones crticas de Heidemann

donde zi es la composicin en la alimentacin del

donde C es el nmero de componentes, A es la

donde P es la presin del sistema, T es la temperatura,

Las expresiones (1) (4) conforman un sistema de

PREDICCIN DE PUNTOS CRTICOS EN MEZCLAS DE GAS NATURAL

donde Tc es la temperatura crtica de componente puro,

d) Propiedades crticas de componente puro.

Tabla 2. Propiedades de componente puro

La Figura 1 presenta el diagrama de flujo del mtodo

Los parmetros de interaccin binaria utilizados

Tabla 1. Parmetros de interaccin binaria kij

c) Parmetros de interaccin binaria empleados.

donde kij es el parmetro de interaccin binaria.

Figura 1. Diagrama de flujo para calcular las

Conciencia Tecnolgica No. 34, Julio - Diciembre 2007

PREDICCIN DE PUNTOS CRTICOS EN MEZCLAS DE GAS NATURAL

El programa de cmputo empleado en este trabajo

donde el superndice cal es el valor de la propiedad

Mezcla de gas natural con ocho componentes

En esta seccin slo se presentan los resultados

La composicin de esta mezcla de gas se presenta en la

Gas de Michelsen (1980)

La Tabla 4 presenta el punto crtico calculado

Tabla 5. Componentes y composiciones de una

Tabla 3. Componentes y composiciones del gas de

El punto crtico calculado para esta mezcla se

Tabla 6. Punto crtico predicho para la mezcla de gas

Conciencia Tecnolgica No. 34, Julio - Diciembre 2007

PREDICCIN DE PUNTOS CRTICOS EN MEZCLAS DE GAS NATURAL

Debido a lo anterior es de esperar que se presenten

L., Platt, G., Henderson, N., (2004), Novel

in Binary Mixtures using Global Optimization,

M.L., (1984), Calculation of Critical

Thermodynamic Modelling: Fundamentals and

A.B., Brennecke, F.J., Khon, P.J., Stadtherr,

Conciencia Tecnolgica No. 34, Julio - Diciembre 2007

Red de Revistas Cientficas de Amrica Latina, el Caribe, Espaa y Portugal

Sistema de Informacin Cientfica

Marco Antonio Ponce Guerrero, Daniel Salazar Sotelo

Revista del Centro de Investigacin. Universidad

Ms informacin del artculo

Mtodo integral para mejorar el

Los incrementos de precios de petrleo crudo

El uso eficiente de la energa ha sido tomado

vaco (GOLAV), gasleo pesado de alto vaco

Figura 1. Diagrama simplificado de una unidad de destilacin combinada.

Rev. Centro Inv. (Mx) Vol. 4, Nm. 16, Ene-Jun. 2001