UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

FES-ZARAGOZA
CARRERA DE INGENIERÍA QUÍMICA
DGAPA-PAPIME: PE104312
ISBN: en trámite

Destilación Flash

Pérez Rodríguez Juan José

José Antonio Zamora Plata

FES ZARAGOZA
UNAM
2012
 DESTILACIÓN FLASH

El separador flash es el equipo de proceso con una etapa de equilibrio con el que se
puede encontrar un diseñador. No obstante, aun cuando interviene una sola etapa, el
cálculo de las composiciones y cantidades relativas de las fases de líquido y vapor para
una temperatura y presión dadas requiere una tediosa solución por tanteo. [Cap. 7 Henley]

OBJETIVOS

 Desarrollar un algoritmo para el cálculo de la destilación flash isotérmica.

 Generar los programas correspondientes en lenguaje FORTRAN, para
el cálculo de constantes de equilibrio usando ecuaciones de estado
cúbicas (Redlich-Kwong-Soave-Mathias y Peng-Robinson-Stryjek-Vera)
para el cálculo de la destilación flash con el método de Rachford-Rice.

1
INTRODUCCIÓN.
Un flash es una sola etapa de destilación en la que la alimentación se vaporiza
parcialmente para producir un vapor más rico en el componente más volátil. En la figura
1.a se presenta una alimentación líquida que se calienta a presión y se somete a una
operación adiabática de un flash mediante el descenso de presión a través de una
válvula, separándose el vapor del líquido residual en tanque flash. Si se suprime la
válvula, en el calentador se puede vaporizar parcialmente por un líquido de baja presión y
separar posteriormente las dos fases. Alternativamente, se puede enfriar una alimentación
de vapor y condensarla parcialmente, separando las fases en un tanque, tal como se
indica en la figura 1.B, para dar lugar a un líquido que es más rico en el componente
menos volátil. En ambos casos, si el equipo esta adecuadamente diseñado, el líquido y el
vapor que salen de la cámara están en equilibrio.

V, yi, HV
P V, T V
Calentador
Alimentación F, zi
líquida HF
TF, PF
Alimentación
líquida Q
L, xi, HL
Fig. 1.A PL, TL

V, yi, HV
Condensador P V, T V
parcial
Alimentación F, zi
vapor HF
TF, PF
Alimentación Q
líquida
Fig. 1.B L, xi, HL
PL, TL

Figura 1. Separadores continuos de una sola etapa (1.A Destilación flash

y 1.B Condensación parcial).

Excepto que la volatilidad relativa sea muy grande, el grado de separación de los
componentes que se puede alcanzar en una sola etapa es bajo. Por esta razón, tanto la
operación flash como la condensación parcial son generalmente operaciones auxiliares

2
para la preparación de corrientes de alimentación que sufrirán un posterior tratamiento.
Sin embargo, los métodos utilizados en el cálculo de una sola etapa son de gran
importancia. Los cálculos desarrollados para una sola etapa de flash son también
ampliamente utilizados para determinar la condición de las fases de una corriente de
composición, temperatura y presión conocidas.

Para las operaciones de una sola etapa de equilibrio que se presentan en la figura
1, con una corriente de alimentación y dos corrientes de producto, las (2C+6) ecuaciones
que se incluyen en la tabla 1 relacionan las (3C+10) variables (F, V, L, z, y, x, T F, TV, TL,
PF, TV, TL, Q), quedando (3C+4) grados de libertad. Suponiendo que (C+2) variables de la
alimentación F, TF, PF y (C-1) valores de zi se desarrollan procedimientos de cálculo para
varios conjuntos diferentes de especificaciones, siendo los más importantes:

1. PL (o PV) y V/F (o L/F): porcentaje de vaporización (porcentaje de

condensación).
2. PL (o PV) y TL (o TV): flash isotérmico.
3. PL (o PV) y Q=0: flash adiabático.

Ecuación No. De ecuaciones Ecuación No.

PV=PL (Equilibrio mecánico) 1 1
TV=TL (Equilibrio térmico) 1 2
Yi=Kixi (Equilibrio entre fases) C 3
Fzi=Vyi+Lxi (Balance de materia, de C-1 4
componentes)
F=V+L (Balance de materia global) 1 5
HFF+Q=HVV+HLL (Balance de energía) 1 6
zi=1 1 7

yi=1 1 8

xi=1 1 9

NE=2C+6

Tabla 1. Ecuaciones para operaciones de una sola etapa de flash y de condensación

parcial.

3
 RELACIONES DE EQUILIBRIO ENTRE FASES.
(Equilibrio Vapor – Líquido).

Las relaciones de equilibrio vapor-líquido son necesarias para resolver muchos

problemas de ingeniería. Los datos requeridos pueden ser (y frecuentemente lo son)
medidos en forma directa mediante la experimentación. Estas mediciones pueden ser
difíciles, incluso para sistemas binarios, y llegan a resultar progresivamente tediosas
cuando aumenta el número de componentes. Éste es el incentivo para aplicar la
termodinámica al cálculo de las relaciones de equilibrio de fases. La termodinámica no
genera datos; proporciona sólo el cuadro matemático que permite el máximo
aprovechamiento de cualquier dato disponible; por tanto, forma la base para lo correlación
sistemática, la extensión, la generalización, la evaluación y la interpretación de datos. Más
aún, proporciona los medios por los cuales pueden aplicarse para propósitos prácticos las
predicciones de varias teorías de física molecular y de mecánica estadística.

El problema fundamental del equilibrio vapor-líquido aparece en un sistema de

multicomponentes, formado por N especies químicas no reactivas, para el cual las
variables de la regla de las fases son: T, P, N-1 fracciones molares en el líquido y N-1
fracciones molares en el vapor. Por ende, se tiene 2N variables. La aplicación de la regla
de las fases establece que F=N, lo cual significa que para un estado en equilibrio sólo N
variables de las 2N son independientes. Una vez que las variables de las regla de las
fases se han especificado, las N variables que quedan pueden determinarse, en principio,
por la solución simultánea de las N relaciones de equilibrio de la forma de la ecuación:

fiV  fiL (10)

donde los superíndices v y l denotan las fases vapor y líquida respectivamente.

En la práctica comúnmente se especifica ya sea T o P y cualquier composición de

la fase líquida o de la fase vapor, fijando 1+(N-1)=N variables de la regla de las fases. Las
N variables restantes constituyen el objetivo de los cálculos, dado que se tienen
suficientes datos disponibles para determinar todas las propiedades termodinámicas
necesarias. La solución matemática formal de este problema fundamental se presenta a
continuación:

4
 La expresión que define a la fugacidad se escribe para cada fase como sigue:

Vapor:

fˆiV  ˆiV yi P (11)

Líquido

fˆi l  ˆil xi P (12)

Sustituyendo estas ecuaciones en la ecuación (10) resulta:

yiˆil  xiˆiv (13)

Las fracciones molares del vapor y del líquido, yi y xi, las cantidades físicas reales,
aparecen ahora en las relaciones de equilibrio, pero ninguna es explícita porque las i son
funciones de la composición. Además, las i son funciones de T y P; por tanto, la
ecuación anterior representa N relaciones complejas vinculando T, P, las xi y las yi.
Incluso para resolver por computación estas ecuaciones, es necesario que las i se
expresen analíticamente como funciones de T, P y la composición. Esto requiere una
ecuación de estado que represente con exactitud la propiedad volumétrica de ambas
fases, líquida y de vapor, para intervalos de temperaturas, presiones y composiciones de
interés. Al utilizar una ecuación de estado se deben establecer los límites de aplicación de
ecuaciones de estado.

Si se van a determinar valores precisos de las constantes de una ecuación de

estado tan compleja, se requiere una gran cantidad de datos volumétricos, incluso para un
componente puro. Estos datos se encuentran solo disponibles para algunas sustancias.
Los datos para mezclas son aun más raros y no hay reglas desarrolladas, teórica o
empíricamente, que puedan aplicarse a mezclas y que proporcionen la dependencia de
las constantes respecto a la composición. La principal dificultad se presenta en el cálculo
del coeficiente de fugacidad de la fase líquida; no se conoce ninguna ecuación de estado
que describa adecuadamente el comportamiento general de la fase líquida. Como
resultado, esta solución general para el problema fundamental del equilibrio vapor-líquido
ha tenido un uso limitado.

5
 ECUACIONES DE ESTADO.

Los procedimientos de diseño del equipo para las operaciones de separación

requieren del conocimiento de entalpías y densidades además de relaciones de equilibrio
entre fases. La termodinámica clásica proporciona un medio de obtener todas estas
magnitudes a partir de P, V y T, que se conocen generalmente como ecuaciones de
estado. Si bien se han propuesto un gran número de ecuaciones P, V y T, relativamente
pocas resultan adecuadas para los cálculos de diseño. Algunas de ellas se relacionan en
la tabla 2. En todas las ecuaciones de la tabla 2 interviene la constante universal de los
gases R y, en todos los casos excepto en 2, intervienen otras constantes que son únicas
para una especie particular. Todas las ecuaciones de estado se pueden aplicar a mezclas
utilizando reglas de mezclado para combinar constantes de las especies puras.

La ecuación del gas ideal, es ampliamente utilizada tanto para gases puros como
para mezclas gaseosas. Esta ecuación no tiene en cuenta el tamaño molecular ni el
potencial de energía de las interacciones moleculares. Cuando cada una de las especies
de una mezcla, así como la mezcla misma, siguen la ley de los gases ideales, se cumplen
la ley de Dalton de la aditividad de las presiones parciales y la ley de Amagat de aditividad
de los volúmenes de las especies puras. La ecuación de mezcla en función de la
densidad molar /M es:

1 P  ni
   i 1 (14)
M v RT V

La ley de los gases ideales es generalmente bastante exacta para presiones bajas
hasta una atmósfera. A 50 psia la ecuación anterior puede presentar desviaciones con
respecto a los datos experimentales de hasta un 10%.

No existe ninguna ecuación de estado sencilla para la fase líquida que no este
basada en el uso de una densidad conocida de las especies líquidas puras y en las
suposiciones de incomprensibilidad y de volúmenes aditivos. Cuando un vapor no se
comporta como un gas ideal resulta difícil formular una ecuación de estado exacta debido
a la necesidad de tener en cuenta las interacciones moleculares.

6
 El principio de los estados correspondientes, que está basado en la similitud de
comportamiento molecular, especialmente en el punto crítico, se puede utilizar para
deducir correlaciones generalizadas, gráficas o tabuladas, del factor de compresibilidad Z
de la ecuación generalizada como una función de la temperatura (absoluta) reducida
Tr=T/Tc, de la presión (absoluta) reducida Pr=P/Pc y de un conveniente tercer parámetro.
Con frecuencia las ecuaciones generalizadas de estado se ajustan a ecuaciones
empíricas de Tr y Pr para su empleo en los métodos de diseño con ordenadores.

La ecuación de estado virial proporciona una buena base teórica para el cálculo de
las relaciones P, V y T de especies puras, tanto polares como no polares, y de mezclas en
fase de vapor. Los coeficientes del virial B, C y superiores, pueden en principio obtenerse
a partir de la mecánica estadística. Sin embargo, el estado actual de los conocimientos es
tal que en la mayor parte de los casos la ecuación se trunca en B, el segundo coeficiente
virial, que se estima a partir de una correlación generalizada. En esta forma, la ecuación
del virial es exacta para densidades de aproximadamente la mitad de la densidad crítica.

La ecuación de cinco constantes de Beattie y Bridgeman, la ecuación de ocho

constantes de Benedict, Webb y Rubin (B-W-R) y la ecuación de dos constantes de
Redlich y Kwong (R-K) son relaciones empíricas aplicables en un amplio intervalo de
presión. La ecuación de R-K es particularmente atractiva debido a que contiene
solamente dos constantes que pueden determinarse directamente a partir de la
temperatura crítica Tc y de la presión crítica Pc. Además, la ecuación de R-K tiene una
exactitud que se compara favorablemente con otras ecuaciones de estado más
complejas, y es capaz de aproximar la región líquida. La ecuación de van der Waals
puede fallar considerablemente en este aspecto.

Cuando la ecuación R-K se aplica a mezclas se utilizan las siguientes formas

equivalentes del factor de compresibilidad:

1 A2 h
Z 
1  h B(1  h)
(15)
 A2  A2
Z 3  Z 2  BP   BP  1 Z  ( BP ) 2  0
 B  B

7
donde las reglas de mezcla son:

C C
A   Ai yi
i 1
o bien Ax
i 1
i i (16)

C C
B   Bi yi o bien B x i i (17)
i 1 i 1

siendo

1/ 2
 0.4278 
1/ 2
 a 
Ai   2 i 2.5    
2.5 
(18)
R T   PcTr 

bi 0.0867
Bi   (19)
RT PcTr
BP
h (20)
Z

Si se utiliza la ecuación R-K directamente, las equivalentes reglas de mezcla para

la fase de vapor son:
C  C 
a     yi y j aij  (21)
i 1  j 1 
donde:
aij  (ai a j )1 / 2 (22)

y
C
b   yibi (23)
i 1

La ecuación (15) es cúbica en Z. Según demostró Edmister, solamente se obtiene

una raíz real positiva para temperaturas supercríticas cuando sólo existe una única fase.
En caso contrario, se obtienen tres raíces reales, correspondiendo el mayor valor de Z a
una fase de vapor y el menor valor de Z a una fase líquida.

8
 PROPIEDADES TERMODINÁMICAS DEDUCIDAS PARA LAS FASES DE
LÍQUIDO Y VAPOR.

Si se dispone de ecuaciones para el calor específico o la entalpía del gas ideal

(presión nula), así como de una ecuación de estado, se pueden deducir en forma
consistente propiedades termodinámicas por aplicación de las ecuaciones de la
termodinámica clásica.
Entalpía.

La capacidad calorífica molar de los gases se expresa convenientemente como un

polinomio de la temperatura

cP0 v  a1  a2T  a3T 2  a4T 3  a5T 4 (24)

La integración de la ecuación anterior conduce a una ecuación de la entalpía molar

del gas ideal a la temperatura T con respecto a la temperatura de referencia T 0,

ak (T k  T0k )
T 5
H   C dT  
0
V
0
PV (25)
T0 k 1 k

Cuando no es válida la suposición de la ley del gas ideal, se utiliza la ecuación 26 para
corregir el efecto de la presión sobre la entalpía. Para una especie pura o para mezclas a
la temperatura T y presión P, la entalpía del vapor es:

 
HV   ( yi H iV0 )  ( HV  HV0 ) (26)
 

La ecuación (26) es especialmente adecuada para utilizar con ecuaciones de

estado que son funciones explícitas de la presión. Las mismas dos ecuaciones pueden
utilizarse para determinar la entalpía de la fase líquida. La aplicación se facilita si la

9
 ecuación de estado es una función continua que pasa entre las regiones de líquido y
vapor. Por tanto,

 
H L   ( xi H iV0 )  ( H L  HV0 ) (27)
 

para especies puras a temperaturas inferiores de la crítica, la ecuación anterior se puede

escindir analíticamente por medio de la siguiente ecuación, donde el subíndice s se refiere
a las condiciones de la presión de saturación.

  P    P    P  
H L  HV0  ( Pv )Vs  RT    P  T    dv  T   (vVS  vLS )  ( Pv ) L  ( Pv ) Ls    P  T   dv
  T v   T  s   T v 

....(28)

La ecuación (27) es la forma preferida, sobre todo si la ecuación de estado es una

función continua.

Si la ecuación de estado R-K se sustituye en la ecuación (28), la integración que

se requiere se realiza en la forma detallada por Edmister, y el resultado se sustituye en la
ecuación (26) y (27), las ecuaciones para la entalpía molar de la mezcla se transforman
en

C  3 A2  BP 
HV   ( yi H iV0 )  RT ZV  1  ln 1   (29)
i 1  2 B  ZV 

C  3 A2  BP 
H L   ( xi H )  RT Z L  1 
0
ln 1   (30)
2 B  Z L 
iV
i 1 

10
 COEFICIENTE DE FUGACIDAD DE UNA ESPECIE PURA.

El coeficiente de fugacidad 0 de una especie pura a la temperatura T y la presión

P se puede obtener directamente a partir de la ecuación de estado. Si P<Psi, 0 es el
coeficiente de fugacidad del vapor y si P>Psi, 0 es el coeficiente de fugacidad del líquido.
La presión de saturación corresponde a la condición 0L=0V. Integrando con la ecuación
de estado R-K se obtiene:

 A2  BP 
V0  exp ZV  1  ln( ZV  BP )  ln 1   (31)
 B  ZV 

 A2  BP 
  exp ZV  1  ln( ZV  BP )  ln 1 
0
V
 (32)
 B  ZV 

COEFICIENTES DE FUGACIDAD DE MEZCLAS.

Los coeficientes de fugacidad de mezclas de especies en fases líquida o de vapor

se pueden obtener a partir de las ecuaciones de estado. Si se aplica la ecuación de
estado de R-K, siguiendo un procedimiento bastante tedioso, como el utilizado por Redlich
y Kwong, se obtienen las siguientes ecuaciones:

 Bi A2  2 Ai Bi   BP 
iV  exp ( ZV  1)  ln( ZV  BP )     ln 1   (33)
 B B A B  ZV 

 Bi A2  2 Ai Bi   BP 
iL  exp ( Z L  1)  ln( Z L  BP )     ln 1   (34)
 B B A B  Z L 

11
Nombre Ecuación Evaluación de funciones o
Constantes de especies
Ley del gas ideal RT Ninguna
P
v
Generalizada Z=(Pr, Tr, Zc, o ) tal como se obtiene a
ZRT
P partir de datos experimentales.
v
Virial (Onnes) A partir de datos volumétricos,
RT  B C  correlaciones generalizadas o
P 1   2  ... mecánica estadística, donde B, C, etc.,
v  v v  dependen de la temperatura.

Van der Waals RT a 27 R 2Tc2

P  2 a
RT
,b  C
v b v 64 PC 8PC

Beattie-Bridgeman Constantes, a, A0, b, B0 y c a partir de

RT  c   b  A  a 
P  2 1  3 v  B0 1    20 1   datos experimentales.
v  vT   v  v  v 

Benedict-Webb-Rubin
(B-W-R) RT  C 1 1 Constantes, a, A0, b, B0, c, C0, ,  a
P   B0 RT  A0  02  2  (bRT  a) 3 partir de datos experimentales.
v  T v v
     
c1  2  exp   2 
a
 6  
v   v 
3 2
v vT
Redlich-Kwong R 2TC2.5 RT
(R-K) RT a a  0.4278 , b  0.0867 C
P  0.5 PC PC
v  b T v (v  b )
Tabla 2. Ecuaciones de estado.
 CALCULO DE FLASH ISOTERMICO.

Si la temperatura de equilibrio TV (o TL) y la presión de equilibrio PV (o PL) de una

mezcla multicomponente están especificadas, los valores de las restantes variables 2C+6
se determinan a partir de la siguiente tabla por medio de un cálculo de flash isotérmico.
Sin embargo, el procedimiento de cálculo no es directo debido a que che no es una
ecuación lineal de las incógnitas V, L yi, y xi. Se han desarrollado numerosas estrategias
de solución, pero el procedimiento generalmente preferido es el de Rachford y Rice
cuando los valores de K son independientes de la composición.

Etapas Numero de la ecuación

TL=TV 35

PL=PV 36

C 37
zC  1   zi
i 1

Resolver:

C
zi (1  Ki )
f     0 38
i 1 1   ( ki  1)

para =V/F, donde Ki=Ki(TV, PV).

39
V=F

zi
xi  40
1   ( Ki  1)

zi Ki
yi   xi Ki 41
1   ( Ki )

L=F-V 42
Q=HVV+HLL-HFF 43

Tabla 3. Procedimiento Rachford-Rice para el cálculo de flash isotérmico.

13
 Primeramente se resuelven las ecuaciones que solamente contienen una
incógnita. Así, se obtienen PL, TL y el valor no especificado de zi a partir de 35 y 36.
Después, puesto que la incógnita Q aparece solamente en 43, dicha incógnita se calcula
solamente una vez que todas las demás ecuaciones han sido resueltas. Queda entonces
3, 4, 5, 8 y 9 para obtener V, L y, y x. Estas ecuaciones pueden agruparse para obtener
las incógnitas de una forma secuencial sustituyendo 5 en 4 para eliminar L; Combinando
el resultado con 3 para obtener dos ecuaciones 40 y 41, una con xi pero no con yi,
mientras que la otra con yi pero no con xi, y sumando estas dos ecuaciones y
combinándolas con 8 y 9 en la forma yi-xi=0 para eliminar yi y xi y obtener una sola
ecuación no lineal solamente con V (o ). Resolviendo V de esta ecuación mediante un
procedimiento iterativo, se obtienen directamente las demás incógnitas a partir de 39 a 43.
Cuando TF y/o PF no están especificadas no se obtiene Q a partir de 43. En este caso
está determinada la condición de equilibrio entre fases para una temperatura (T V=TL) y
una presión (PV=PL) conocidas.

La ecuación 38 se puede resolver por tanteo probando valores de  comprendidos

entre 0 y 1 hasta que f()=0.

Los métodos de cálculo mas ampliamente utilizados para la resolución de che se

basan en el método de Newton. Se calcula un valor de predicción de la raíz  para la
iteración k+1 a partir de la relación recursiva:

 ( k 1)   ( k ) 
 
f  (k )
 
f '  (k )
(44)

donde la derivada de 44 es:

zi (1  Ki ) 2
 
C
f '  (k )   (45)
i 1 [1   ( ki  1)]
k 2

14
 CÁLCULO DE FLASH ADIABÁTICO PARA MULTICOMPONENTES.

Cuando la presión de una corriente líquida de composición, velocidad de flujo y

temperatura (o entalpía) conocidas, se reduce adiabáticamente a través de una válvula
como se muestra en la figura 1.A, se puede realizar un cálculo de flash adiabático para
determinar la temperatura, composiciones y velocidades de flujo que resultan para las
corrientes de vapor y líquido en equilibrio para una presión especificada de la cámara
flash. Para este caso se aplica el procedimiento de flash isotérmico en forma iterativa y se
supone una temperatura TV como variable de iteración, cuyo valor se supone. Después se
determinan , x, y, y L como para un flash isotérmico. El valor supuesto de T V (igual a TL)
se contrasta mediante un balance de entalpía obtenido combinando 43 para Q=0 con 42,
con el fin de obtener:

HV  (1   ) H L  H F
f (TV )  0 (46)
1000

dónde la división por 1000 da lugar a 46 del orden de la unidad. El método de falsa
posición se puede utilizar para iterar sobre la temperatura TV hasta que |TV(k+1)-TV(k)|
/TV(k)<0.0001. el procedimiento se facilita cuando las mezclas forman soluciones ideales, o
cuando los valores de K pueden estimarse a partir de Ki=0iL.

Este método resulta muy conveniente cuando 38 no es sensible a T V. Este es el

caso de mezclas con amplio intervalo de ebullición. Para mezclas con estrecho intervalo
de ebullición (p.e. isómeros) el método puede fracasar si 38 es extremadamente sensible
al valor de TV. En este caso es preferible seleccionar  como variable de tanteo y obtener
TV por resolución iterativa de 38

C
zi (1  Ki )
f TV    0 (47)
i 1 1   ( ki  1)

después se obtienen x e y a partir de 39 y 40, respectivamente, para luego obtener

directamente  de 46 teniendo en cuenta que

15
 HV  (1   ) H L  H F
f ( )  0 (48)
1000

a partir de la cual

HF  HL
 (49)
HV  H L

Si  obtenido a partir de 49 no es igual al valor de  supuesto para resolver 47, se

utiliza el nuevo valor de  para comenzar nuevamente el lazo con 47.

16
 APENDICE 1

CÁLCULO DEL EQUILIBRIO LIQUIDO - VAPOR USANDO LA ECUACIÓN DE

ESTADO REDLICH-KWONG-SOAVE.

Dos versiones modificadas de la ecuación Redlich-Kwong se describen aquí. Se

describe propiamente el equilibrio líquido-vapor de especies saturadas no polares.

Calculo de los coeficientes de fugacidad y actividad con la ecuación de estado

Redlich-Kwong-Soave.

Experimentos han mostrado que la ecuación original de la ecuación de estado

Redlich-Kwong no es muy satisfactoria para cálculos de equilibrio líquido-vapor. Por
consiguiente Soave modifico la forma original reemplazando la constante a de la ecuación
con una función de temperatura envolviendo también el factor omega de Pitzer. Así, la
ecuación de estado modificada queda:

RT aT 
P 
  b  (  b)

Basándonos en consideraciones fundamentales, es requerido que a sea dada en

el punto crítico por:

R 2TC2
aC  0.42747
PC

Así, a cualquier temperatura:

a  aC T 

Dónde [T] representa un factor de corrección que es una función solamente de la

temperatura con el factor omega como un parámetro. ¿Cómo encontraremos esta
función?

17
 La ecuación de estado, después de rearreglarse, puede ser escrita para z:

z 3  z 2  z( A  B  B2 )  AB  0

dónde:

aT P
A
R 2T 2

bP
B
RT

Esta ecuación también es válida para la presión de vapor en equilibrio, y las

soluciones para z, z(l) que son los factores de compresibilidad del vapor y líquido
respectivamente. Es más conocido que:

f P  1
ln     dP
P 0
 RT P 

Usando la ecuación de estado RKS, para el coeficiente de fugacidad tenemos:

f A Z  B
ln ( Z  1)  ln( Z  B)  ln  
P B  Z 

Para el equilibrio, el coeficiente de fugacidad del vapor saturado es igual que el del
líquido saturado:

f f (l )
 i  ln  ln  i
P0 P0

Lo que significa por otro lado, que a cualquier presión P0, correspondiente a una
temperatura dada T, hay un simple valor para A, el cual hará que la diferencia

18
desaparezca entre el lado izquierdo y derecho de la ecuación anterior. Los valores
resultantes pueden ser efectivamente linearizados como:

  1  mi (1  Tri )

donde: i = factor de corrección a ser determinado.

mi = parámetro específico de la sustancia para el componente i.

El tercer parámetro en la ecuación de estado es el factor omega. Los parámetros

wi y mi son funcionalmente relacionados, y la dependencia es no linear. Asumiendo una
relación cuadrática, una estimación del parámetro resulta en:

mi  0.480  1.57wi  17wi2

Así es como la situación parece para un sistema de un solo componente. Si hay

mas de un componente en el sistema, la pregunta es como usar la ecuación de estado
RKS para el cálculo del coeficiente de fugacidad del equilibrio entre las fases. Para el
coeficiente de fugacidad la siguiente relación termodinámica fue derivada:

1 
fi   1  P  d  ln Z
Pxi 
ln 
  RT  n  
  i T , P , n j 

Después de sustituir la ecuación de estado y las siguientes reglas de mezclado

que son usadas:

a   ( xi ai0.5 ) 2

b   xibi
i
Tenemos:

fi b A  a 0.5 b   B 
ln  i ( Z  1)  ln( Z  B)   2 i0.5  i  ln 1  
Pxi b B a b  Z

19
 Entonces, si la ecuación de estado es aplicada a cálculos de equilibrio, entonces
de a(T) a la temperatura T puede ser obtenido por el conocimiento de wi, de la siguiente
manera:

i Tr   (1  mi (1  Tri ))2

R 2TC2
ai [Tr ]  i [Tr ]0.42747
PC

RTCi
bi  0.08664
PCi

Para los parámetros envueltos en la ecuación para determinar el coeficiente de

fugacidad obtenemos:

ai0.5  i0.5TCi / PCi0.5


a 0.5  xii0.5TCi / PCi0.5 
i

bi TCi / PCi

b  ( xiTCi / PCi )
i

2
P  x  0.5T 
A  0.42747 2   i i 0.5 Ci 
T  i PCi 

P xiTCi
B  0.08664 
T i PCi

Si los valores de A y B obtenidos con las fórmulas anteriores son sustituidos en la

ecuación de factor de compresibilidad, tenemos:

Z 3  Z 2  Z ( A  B  B 2 )  AB  0

Resolviendo para ambas fases vapor y líquido, los valores para determinar el
coeficiente de fugacidad son obtenidos.

20
 METODOLOGIA.

ALGORITMO PARA ÉL CÁLCULO DE UN FLASH ADIABÁTICO E

ISOTÉRMICO MULTICOMPONENTE CON LA ECUACIÓN DE ESTADO DE
REDLICH-KWONG-SOAVE.

1. Leer los componentes y su fracción mol respectivos en la mezcla.

2. Leer la temperatura y la presión.
3. Calcular i para cada componente en la mezcla.
4. Calcular ai y bi para cada componente en la mezcla.
5. Calcular a y b con las reglas de mezclado correspondientes para la mezcla.
6. Calcular los coeficientes A y B para la ecuación de estado.
7. Resolver la ecuación de estado por iteración directa para el cálculo del factor de
compresibilidad de las fases líquida y vapor.
8. Calcular los coeficientes de fugacidad de la fase líquida y vapor para cada
componente iL y iV.
9. Calcular las constantes de equilibrio de cada componente en la mezcla Ki.
10. Calcular iterativamente el porcentaje de vaporización de la mezcla .
11. Calcular las fracciones mol del líquido y el vapor ya separados de cada componente
en la mezcla xi y yi. (*)
12. Verificar el porcentaje de error para , en caso de no ser el deseado regresar al paso
número 4.

21
 PROBLEMA RESUELTO Y RESULTADOS.
Para una mezcla de 26.54% mol de etano y 73.46% mol de n-heptano a 10 bar y
430 K, calcular la fracción del líquido y la composición de las dos fases líquido y vapor.
Usando la ecuación de estado Redlich-Kwong-Soave.
Solución:
Tc (K) Pc (bar) w
Etano 305.4 48.8 .099
n-heptano 540.3 27.4 .349

Con la metodología que se dio y las formulas dadas en el apéndice 1, se obtienen

los siguientes resultados de forma manual y después se dan los resultados generados
con el programa.
Para el caso de la solución que se hizo en forma manual se pondrán la primera y
ultima iteración y con el programa solo se pondrán los resultados.

1ª. Iteración:

Tr m  a b A B L v K x y  error
Etano 1.40 .5246 .8137 4.595 .045 - - - - 15.83 .0242 .715 .2476 .5046
n-heptano .79 .880 1.198 37.74 .142 - - - - .1588 .9757 .287 .2478 .5048

5ª. Iteración:

Tr m  a b A B L v K x y  error
Etano 1.40 .5246 .8137 4.595 .045 .2860 .0390 9.04 1.014 8.914 .0339 .354 .419 .0006
n-heptano .79 .880 1.198 37.74 .142 .0848 .0204 .4411 .773 .5701 .9660 .645 .419 .006

22
Con el programa se obtienen los siguientes resultados:

L v K x y 
Etano 9.5001 .4011 9.3790 0.0593 .5369 .5769
n-heptano 1.0129 .7871 .5096 .9407 .4631 .5769

CONCLUSIONES.

De lo anterior podemos concluir que el uso de un programa para cálculos en la

ingeniería es de gran utilidad y confianza, en este trabajo se cumple con los objetivos
mencionados, ya que se llegó a la generación de un programa para el cálculo de un flash
isotérmico, de manera que este pueda ser accesado por los alumnos de la carrera de
ingeniería química, para que estos se puedan involucrar a la práctica de ejercicios del
modulo de procesos de separación, ya sea de manera manual o por medio de un
programa.

23
 CODIGO DE PROGRAMACIÓN EN LENGUAJE FORTRAN
PROGRAM FLASH
DIMENSION PC(10), TC(10), NCOMP(10), NUCOMP(10), W(10), PM(10), TR(10), FRACC(10), FUG(10)
DIMENSION ALFA(10), AA(10), BB(10), M(10), K(10), ACT(10), AM(10,10), BM(10), Q(10)
DIMENSION X(10), YX(10), Y(10), YD(10)
REAL :: K, VAP
INTEGER :: N, NUCOMP
CHARACTER (LEN=15) :: NCOMP
PRINT *,"INTRODUCIR TEMPERATURA ( K)"
READ (*,*) T
PRINT *,"INTRODUCIR PRESION (kPa)"
READ (*,*) P
PRINT *,"INTRODUCIR EL NUMERO DE COMPUESTOS (DE 2 a 10)"
READ (*,*) N
PRINT *,"INTRODUCIR EL NUMERO DEL COMPUESTO Y SU FRACCION MOL"
DO I=1,N
CALL COMPUESTOS(NUCOMP(I), NCOMP(I), PM(I), PC(I), TC(I), W(I), Q(I))
READ (*,*) FRACC(I)
END DO
PRINT *,"________________________________"
PRINT *," No. Compuesto Composicion"
PRINT *,"________________________________"
DO I=1,N
PRINT 150, NUCOMP(I), NCOMP(I), FRACC(I)
END DO
! CALCULO DE LAS CONSTANTES DE EQUILIBRIO
P=P*1000
R=831430
DO I=1,N
TR(I)=T/TC(I)
END DO
DO I=1,N
M(I)=0.48508+1.55171*W(I)-0.15613*(W(I)**2)
END DO
IF (T<=TC(I)) THEN
DO I=1,N
ALFA(I)=(1+M(I)*(1-SQRT(TR(I)))-Q(I)*(1-TR(I))*(0.7-TR(I)))**2
END DO
ELSE
DO I=1,N
ALFA(I)=EXP(2*((1+.005+0.3*Q(I)*M(I))-1)*(1-TR(I)**(1+.005+0.3*Q(I)*M(I)))/(1+.005+0.3*Q(I)*M(I)))

24
 END DO
END IF
DO I=1,N
AA(I)=(0.42747*(R*TC(I))**2*ALFA(I))/PC(I)
BB(I)=(0.08664*R*TC(I))/PC(I)
END DO
DO I=1,N
DO J=1,N
AM(I,J)=FRACC(I)*FRACC(J)*SQRT(AA(I)*AA(J))
END DO
END DO
DO I=1,N
BM(I)=FRACC(I)*BB(I)
END DO
SA=AM(1,1)+AM(1,2)+AM(1,3)+AM(1,4)+AM(1,5)+AM(1,6)+AM(1,7)+AM(1,8)+AM(1,9)+AM(1,10)+&
AM(2,1)+AM(2,2)+AM(2,3)+AM(2,4)+AM(2,5)+AM(2,6)+AM(2,7)+AM(2,8)+AM(2,9)+AM(2,10)+&
AM(3,1)+AM(3,2)+AM(3,3)+AM(3,4)+AM(3,5)+AM(3,6)+AM(3,7)+AM(3,8)+AM(3,9)+AM(3,10)+&
AM(4,1)+AM(4,2)+AM(4,3)+AM(4,4)+AM(4,5)+AM(4,6)+AM(4,7)+AM(4,8)+AM(4,9)+AM(4,10)+&
AM(5,1)+AM(5,2)+AM(5,3)+AM(5,4)+AM(5,5)+AM(5,6)+AM(5,7)+AM(5,8)+AM(5,9)+AM(5,10)+&
AM(6,1)+AM(6,2)+AM(6,3)+AM(6,4)+AM(6,5)+AM(6,6)+AM(6,7)+AM(6,8)+AM(6,9)+AM(6,10)+&
AM(7,1)+AM(7,2)+AM(7,3)+AM(7,4)+AM(7,5)+AM(7,6)+AM(7,7)+AM(7,8)+AM(7,9)+AM(7,10)+&
AM(8,1)+AM(8,2)+AM(8,3)+AM(8,4)+AM(8,5)+AM(8,6)+AM(8,7)+AM(8,8)+AM(8,9)+AM(8,10)+&
AM(9,1)+AM(9,2)+AM(9,3)+AM(9,4)+AM(9,5)+AM(9,6)+AM(9,7)+AM(9,8)+AM(9,9)+AM(9,10)+&

AM(10,1)+AM(10,2)+AM(10,3)+AM(10,4)+AM(10,5)+AM(10,6)+AM(10,7)+AM(10,8)+AM(10,9)+AM(10,10)
SB=BM(1)+BM(2)+BM(3)+BM(4)+BM(5)+BM(6)+BM(7)+BM(8)+BM(9)+BM(10)
A=(SA*P)/(R**2*T**2)
B=(SB*P)/(R*T)
CALL COMPRESIBILIDAD(A, B, ZL, ZV)
DO I=1,N
DO J=1,N
ACT(I)=EXP(BB(I)/SB*(ZL-1)-LOG(ZL-B)-A/B*(2*AA(I)/(SA)-BB(I)/SB)*LOG(1+B/ZL))
FUG(I)=EXP(BB(I)/SB*(ZV-1)-LOG(ZV-B)-A/B*(2*AA(I)/(SA)-BB(I)/SB)*LOG(1+B/ZV))
END DO
END DO
DO I=1,N
K(I)=ACT(I)/FUG(I)
END DO
VAP=0.5
VAPB=VAP
J=0
40 J=J+1

25
 DO I=1,N
Y(I)=FRACC(I)*(1-K(I))/(1+VAP*(K(I)-1))
YD(I)=(FRACC(I)*((1-K(I))**2))/((1+VAP*(K(I)-1))**2)
END DO
YY=Y(1)+Y(2)+Y(3)+Y(4)+Y(5)+Y(6)+Y(7)+Y(8)+Y(9)+Y(10)
YYD=YD(1)+YD(2)+YD(3)+YD(4)+YD(5)+YD(6)+YD(7)+YD(8)+YD(9)+YD(10)
VAP=VAP-YY/YYD
IF (ABS(VAPB-VAP).GE.1E-5) THEN
VAPB=VAP
GOTO 40
END IF
IF (VAP>=1.OR.VAP<=0) THEN
PRINT *,"PRESION FUERA DE RANGO"
END IF
DO I=1,N
X(I)=FRACC(I)/(1+VAP*(K(I)-1))
END DO
DO I=1,N
YX(I)=K(I)*X(I)
END DO
PRINT *,"_______________________________________________________________"
PRINT *,"Compuesto K(i) x(i) y(i) "
PRINT *,"_______________________________________________________________"
DO I=1,N
PRINT *, NCOMP(I), K(I), X(I), YX(I)
END DO
PRINT *,"_______________________________________________________________"
PRINT *,"Compuesto K(i) Act(i) Fug(i) "
PRINT *,"_______________________________________________________________"
DO I=1,N
PRINT *, NCOMP(I), K(I), ACT(I), FUG(I)
END DO
PRINT 160, VAP
PRINT 170, ZL
PRINT 175, ZV
150 FORMAT (1X,I3,5X,A15,2X,F5.3)
160 FORMAT (//," El porcentaje de vaporizacion es ",F8.6,//)
170 FORMAT (" El factor de comprersibilidad de la fase liquida es ", F8.6)
175 FORMAT (" El factor de compresibilidad de la fase vapor es ", F8.6)
PAUSE 1
END PROGRAM
SUBROUTINE COMPUESTOS(NUCOMP, NCOMP, PM, PC, TC, W, Q)

26
 INTEGER :: NUCOMP
CHARACTER (LEN=15) :: NCOMP
READ (*,*) N
OPEN (UNIT=61, FILE="datos.TXT", STATUS="OLD", ACTION="READ")
DO I=1,N
READ (UNIT=61, FMT=*) NUCOMP, NCOMP, PM, PC, TC, W, Q
END DO
REWIND (UNIT=61)
RETURN
END
SUBROUTINE COMPRESIBILIDAD(A, B, ZL, ZV)
Z=1
ZB=Z
I=0
30 I=I+1
Y=Z**3-Z**2+Z*(A-B-B**2)-(A*B)
YD=3*Z**2-2*Z+(A-B-B**2)
Z=Z-Y/YD
IF (ABS(Z-ZB)>=1E-6) THEN
ZB=Z
GOTO 30
END IF
Z1=Z
A1=Z1
A2=Z1
A3=(A-B-B**2)*Z1
Z2=(A2+SQRT(A2**2-4*A1*A3))/2*A1
Z3=(A2-SQRT(A2**2-4*A1*A3))/2*A1
IF (Z1>Z2.AND.Z1>Z3) THEN
ZV=Z1
ELSE
ZL=Z1
END IF
IF (Z2>Z1.AND.Z2>Z3) THEN
ZV=Z2
ELSE
ZL=Z2
END IF
IF (Z3>Z2.AND.Z3>Z1) THEN
ZV=Z3
ELSE
ZL=Z3

27
 END IF
RETURN
END SUBROUTINE

LISTA DE COMPUESTOS
1, Metano, 1.6430E-2, 4.6001E+6, 190.560, 4.06115E-03, 0
2, Etano, 3.0070E-2, 4.879E+6, 305.460, 9.40157E-02, 0
3, Propano, 4.4096E-2, 4.24552E+6, 369.860, 0.148716, 0
4, n-Butano, 5.8124E-2, 3.7997E+6, 425.2, 0.196726, 0
5, n-Pentano, 7.21509E-2, 3.37410E+6, 469.8, 0.248393, 0
6, n-Hexano, 8.61779E-2, 3.0296E+6, 507.9, 0.296379, 0
7, n-Heptano, 0.100206, 2.74E+6, 540.2, .348687, 0
8, n-Octano, 0.114233, 2.4926E+6, 569.4, 0.394001, 0
9, n-Nonano, 0.12826, 2.2798E+6, 595, 0.440465, 0
10, n-Decano, 0.142287, 2.1076E+6, 619, 0.479366, 0
11, n-Undecano, 0.156314, 1.9454E+6, 640, 0.521348, 0
12, n-Dodecano, 0.170341, 1.8137E+6, 659, 0.567165, 0
13, n-Hexadecano, 0.226449, 1.4185E+6, 725, 0.581241, 0
14, 2-Metilpropano, 5.8124E-2, 3.6477E+6, 408.1, 0.180685, 0
15, 2-Metilbutano, 7.21509E-2, 3.3336E+6, 461, 0.215465, 0
16, 2;2-Dimetilpropano, 7.21509E-2, 3.2019E+6, 433.8, 0.194894, 0
17, 2;3-Dimetilbutano, 8.61779E-2, 3.1309E+6, 500.3, 0.243007, 0
18, Hidrogeno, 2.016E-3, 1.2996E+6, 32.96, -.224469, 0
19, Benceno, 7.811E-2, 4.894E+6, 562.1, .210503, 0
20, Tolueno, 9.21E-2, 4.16446E+6, 591.96, .255292, 0
21, Etilbenceno, , 3.606e+6, 617.20, .3027,0
22, Acetona, , 4.696e+6,508.1, .30667,0
23, Metanol, , 8.095e+6, 512.58, .56533, -.16816
24, Etanol, , 6.148e+6, 513.92, .64439, -.03374
25, Nitrogeno, , 3.4e+6, 126.2, .03726, .01996
26, Oxigeno, , 5.09e+6, 154.77, .02128, .01512
27, Dio._de_carbono, , 7.383e+6, 304.21, .225, .04285
28, Amonio, , 11.289e+6, 405.55, .2517, .001
29, Agua, , 22.089e+6, 647.286, .3438, -.06635
30, Ac._clorhidrico, , 8.308e+6, 324.6, .12606, .01989

28