REPUBLICA BOLIVARIANA DE VENEZUELA

UNIVERSIDAD NACIONAL EXPERIMENTAL

“RAFAEL MARÍA BARALT”
MINISTERIO DEL PODER POPULAR PARA LA EDUCACIÓN
VICERRECTORADO ACADÉMICO

EQUILIBRIO LÍQUIDO-VAPOR DE SISTEMAS MULTICOMPONENTES

Profesor: Autora:
Ing. Yohana Madriz Br. Maria Ybañez C.I. 25945430

Santa Ana de Coro, mayo de 2022

 INTRODUCCIÓN

El equilibrio de las fases líquido y vapor (ELV) es de gran utilidad en muchos

procesos de separación y purificación, por lo que se ha dirigido un extenso número
de investigaciones teóricas y experimentales para su mejor entendimiento. A pesar
de que las mezclas binarias son un poco comunes, un amplio conocimiento de
ellas resulta esencial, ya que la extrapolación de su comportamiento es necesaria
para caracterizar a las mezclas multicomponentes.

En tal sentido, en el aprendizaje del equilibrio líquido-vapor, es necesario

conocer y entender conceptos útiles para su interpretación y explicación, como lo
son: la fugacidad, el coeficiente de fugacidad y otros derivados de los criterios
para el equilibrio de fases.

Por otra parte, el cálculo de las condiciones de equilibrio de fases líquido vapor
en mezclas multicomponentes es un tema de interés general para la ingeniería, ya
que sirve de base para el diseño de equipos de separación asociados con dichas
fases.

Finalmente, en la presente investigación se establecen los cálculos de los

puntos de Burbuja y Rocío que pueden ser estimados a través de correlaciones
matemáticas relativamente sencillas en aproximaciones para sistemas de
comportamiento ideal y no ideal.
 Equilibrio líquido-vapor de sistemas multicomponentes

[Link] líquido vapor

El equilibrio líquido vapor es el estado en que se encuentran las fases liquida
y vapor de un sistema de componentes al que se ha dado suficiente tiempo
para que alcance un estado estacionario, en el cual la composición de las
fases difiere entre sí, y se tiene una energía mínima de Gibbs.

La determinación del equilibrio permite el estudio del comportamiento de

sistemas, en los cuales coexisten dos o más componentes. Bajo este
fundamento se tratan los procesos de separación, en los cuales se busca
enriquecer o purificar uno de los componentes del sistema para obtener un
beneficio, ya sea para su aplicación en un proceso posterior o para lograr un
producto final.

Cuando se presenta coexistencia entre las fases líquido y vapor

necesariamente se debe considerar la existencia de un equilibrio entre las
mismas. La primera aproximación a este sistema es el modelo ideal; el modelo
asume idealidad tanto en la fase liquida como en la fase vapor, dicho sistema
puede ser descrito por la Ley de Raoult; que es aplicable a presiones bajas o
moderadas y si las especies involucradas son semejantes.

Sat
y i ∙ P= xi ∙ Pi (1)

La Ley de Raoult asegura que la presión que se presenta en la fase de vapor

es debida a las partículas que se evaporan provenientes del líquido, en otras
palabras, la presión ejercida sobre la fase de vapor corresponde a la suma de
partículas de cada especie que está saliendo de la fase líquida. Esta fue una
buena aproximación para un grupo moderado de sistemas de equilibrio, pero se
decidió abandonar la suposición de líquido ideal para desarrollar un modelo
matemático que describiera mejor un mayor grupo de sistemas. En otras
palabras, para sistemas altamente alejados de la idealidad (como los
azeotrópicos) es necesario utilizar un parámetro de corrección que correlacione
 las fases.
Sat
y i ∙ P=γ i ∙ xi ∙ Pi (2)

[Link] líquido vapor en sistemas multicomponentes

1.
2.
Los sistemas no ideales de tres componentes pueden visualizarse
gráficamente, utilizando coordenadas triangulares para expresar las
composiciones; no obstante, para más de tres componentes, la evaluación
gráfica se vuelve muy compleja. En la actualidad, el conocimiento acerca de los
sistemas no ideales de este tipo es extremadamente limitado; más aún, son
relativamente pocos los datos que se han reunido. Muchos de los sistemas de
multicomponentes de importancia industrial pueden considerarse casi ideales
en la fase líquida para todos los fines prácticos. Esto es particularmente cierto
para mezclas de hidrocarburos de la misma serie homóloga, como las mezclas
de la serie parafinica o los hidrocarburos aromáticos de bajo punto de
ebullición. En estos casos, puede aplicarse la Ley de Raoult o su equivalente
en función de las fugacidades y calcular el equilibrio a partir de las propiedades
de los componentes puros. Sin embargo, generalmente no es seguro predecir
el comportamiento detallado de un sistema de multicomponentes a partir,
únicamente, de la consideración de los componentes puros, ni siquiera del
conocimiento de los sistemas binarios simples que pueden formarse a partir de
los componentes. Por ejemplo, algunas veces los sistemas de tres
componentes forman azeotropos ternarios, cuyo vapor y liquido en el equilibrio
tienen composiciones idénticas. Sin embargo, el hecho de que existan uno, dos
o tres azeotropos binarios entre los componentes no asegura la formación de
un azeotropo ternario; además, un azeotropo ternario no coincide
necesariamente con la composición a la temperatura del punto mínimo o
máximo de burbuja para el sistema a presión constante. La complejidad de los
sistemas hace que las computadoras digitales sean útiles para manejar los
datos.
 Un sistema azeotrópico se presenta cuando bajo condiciones específicas
las composiciones del líquido y el vapor para un componente i son iguales;
este es un comportamiento muy alejado de la idealidad. Si las desviaciones
son grandes y positivas respecto de Raoult, el sistema genera un azeotropo de
mínimo punto de ebullición en un diagrama T-xy. De otra parte, si el sistema
presenta desviaciones grandes y negativas respecto de Raoult, se genera un
azeotropo de máximo punto de ebullición en un diagrama T-xy

Por otra parte, para los sistemas de multicomponentes se acostumbra

describir los datos en el equilibrio mediante el coeficiente de distribución K.
Para el componente j,
yj
K j= (3)
xj

en donde K j depende por lo general de la temperatura, presión y composición de

la mezcla. La volatilidad relativa α ij del componente 1 con respecto a j es
y i / xi K i
α ij = = (4)
y j/ x j K j

Para soluciones ideales a presión moderada, K, es independiente de la

composición y sólo depende de la temperatura (puesto que ésta afecta a la
presión de vapor) y la presión total
Sat
Pj
K j= (5)
Ptot
Y
Ptot
α ij = Sat
(6)
Pj

Punto de burbuja para el vapor en el punto de burbuja,

∑ yi =1(7)
K A ∙ x A + K B ∙ x B + K C ∙ x C +⋯=1(8)

Con el componente j como componente de referencia

 K A ∙ x A K B ∙ x B K C ∙ xC 1
+ + + ⋯= (9)
Kj Kj Kj Kj
1
α Aj ∙ x A +α Bj ∙ x B + α Cj ∙ x C +⋯=∑ α i , j ∙ x i= (10)
Kj

La composición del vapor en el punto de burbuja en el equilibrio está dada por

α i , j ∙ xi
y j= (11)
∑ αi , j ∙ xi
Habiendo fijado la composición del líquido y la presión total, el cálculo de la
temperatura se hace por ensayo hasta satisfacer la ecuación (10). La
convergencia es rápida, puesto que las α sólo varían lentamente al cambiar la
temperatura.

Punto de rocío para el líquido en el punto de rocío,

∑ x i=1(12)
y A yB yC
+ + + ⋯=1(13)
KA KB KB

Con el componente J como componente de referencia

K j ∙ y A K j ∙ yB K j ∙ yC
+ + + ⋯=K j (14 )
KA KB KC
y A y B yC yi
+ + +⋯=∑ =¿ K j (15)¿
α Aj α Bj α Cj αi,j

La composición del líquido en el punto de rocío está dada por

y i /α i , j
x i= (16)
∑ ( y i /αi , j )
2.1. Termodinámica del equilibrio líquido vapor (Métodos de predicción)
El equilibrio termodinámico entre las fases vapor y líquida de un
sistema multicomponentes requiere tres condiciones:

Equilibrio térmico, para lo cual la temperatura debe ser igual en

ambas fases;
 Equilibrio mecánico, que implica igualdad de la presión en ambas
fases y
Equilibrio químico, que exige la igualdad de la fugacidad en la
mezcla de cada componente en cada fase.

Esta última condición se expresa de distintas maneras, según el grado

de idealidad que se admita para el sistema. Así:

a) Si ambas fases se comportan en forma ideal, el equilibrio se

plantea a través de la ley de Raoult
Sat
y i ∙ P= xi ∙ Pi (17)

b) Si existe una ligera o moderada no idealidad en ambas fases se

puede utilizar una formulación del equilibrio que incluye sendos
coeficientes de fugacidad ∅ , los que se calculan a partir de una
ecuación de estado para sistemas reales: ∅ iV ∙ y i ∙ P=∅ iL ∙ x i ∙ PiSat (18)

La fugacidad, ∅ , mide la energía de Gibbs de un gas real

(mientras la presión mide la energía de Gibbs de un gas ideal), y
depende de la temperatura y la presión. En gases no ideales, la
fugacidad debe ser uniforme en todo el recipiente, y como es una
función de la temperatura y la presión, a una temperatura dada,
valores iguales de fugacidad implican valores iguales de presión.
Para un componente de una mezcla ideal se tiene la fugacidad
como una presión parcial corregida y el comportamiento de una
solución no electrolítica se aproxima asintóticamente al
comportamiento lineal.

c) Si la fase líquida presenta una marcada no linealidad, por ejemplo,

debido a la existencia de puentes de hidrógeno, acudiremos al
modelo de soluciones, a través del coeficiente de actividad γ :
V Sat Sat
y i ∙ ∅i ∙ P=γ i ∙ x i ∙ ∅ i ∙ Pi ∙ POY i (19)
 En la expresión anterior se agrega, además, el factor de Poyting
POY i , que permite tener en cuenta la influencia de la variación de la
presión entre la presión de vapor PSat y la del sistema P.

La actividad fue caracterizada por Lewis como una cantidad que

define qué tan activa es una solución comparada con su
comportamiento ideal. En una solución ideal el número de moles de
las especies se define por medio de una escala de concentración,
sea molaridad, molalidad o fracción molar; así, para una solución
ideal la actividad es igual a la concentración, y para una solución
real, se tiene, en donde γ es un coeficiente de actividad, que para
las soluciones ideales es igual a la unidad.

a)

b)

c)

d) Debemos tener presente que el apartamiento de la idealidad se

produce en la fase líquida en condiciones de presión y temperatura
menos severas que para la gaseosa; esto es, la fugacidad en el
vapor se sigue expresando por y i P pero en el líquido se requiere
hacer uso del modelo de soluciones. Tendremos, en este caso, una
ley de Raoult modificada como expresión de la igualdad entre
fugacidades:
Sat
y i ∙ P=γ i ∙ xi ∙ Pi (20)

La consideración del ELV en el caso de gases como el N 2

merece un párrafo aparte. En este caso, se trata de un gas que se
encuentra disuelto en el líquido y la fugacidad no tiene relación,
como antes, con un cambio de fase, por lo que corresponde aplicar
una extensión de la ley de Henry más que la de Raoult. Así, para
esos casos, la igualdad de fugacidades quedará expresada como:
 V mix
∅ i ∙ y i ∙ P=H i ∙ x i (21)

Entonces los coeficientes de fugacidad se calculan a partir de

ecuaciones de estado para sistemas reales. Las más comunes son:
Soave-Redlich-Kwong (SRK) y sus variantes
Peng-Robinson (PR) y sus variantes

Todas estas ecuaciones requieren el uso de coeficientes de

interacción binarias para considerar adecuadamente las mezclas
multicomponentes.

Para el cálculo de coeficientes de actividad se cuenta entre las

más comunes con:
Ecuación de Wilson
Ecuaciones de Margules
Modelo Non Random Two Liquids (NRTL)
Modelo UNIversal QUAsi Chemical (UNIQUAC)

Nota: Las ecuaciones de estado y coeficientes de actividad son

modelos termodinámicos ajustables para diferentes sistemas de
componentes, pero complejos y extensos. Es por ello que, no son
definidas cada una de las ecuaciones matemáticas representativas
en esta investigación.

Ecuación de Antoine
La forma más usual de expresar la función de Antoine es:

B
ln ( P iSat )= A− (22)
C+T

B
A−
Sat C+T
Pi =e (23)

En esta fórmula P representa la presión de vapor de saturación expresada en

milímetros de mercurio (mmHg), T es la temperatura que viene siendo la variable
independiente y se expresa en ℃. Las letras A, B y C son las constantes o
parámetros de la fórmula de Antoine.

La importancia de esta fórmula que, a pesar de ser empírica, da una expresión

analítica sencilla que puede ser usada con facilidad en los cálculos
termodinámicos.

Ejemplo 1. (Apartado 2: Equilibrio líquido vapor en sistemas multicomponentes)

Una solución de hidrocarburos a una presión total de 350 KN /m2 ( 508 lbf /¿2 ) tiene el
siguiente análisis
n−C3 H 8=5.0 % , n−C 4 H 10=30.0 % , n−C5 H 12=40.0 % , n−C 6 H 14 =25.0 % en mol.
Calcular los puntos de burbuja y de rocío.

Solución. Los valores de K se van a tomar de la monografía de DePriester [Chem.

Eng. Progr. Symp. Ser., 49(T), 1 (1953), y de Chem. Engrs. Handbook, 5a., ed.,
Fig. 13-661, que supone soluciones líquidas ideales.

Determinación del Punto de Burbuja

Se establece una tabla con los datos y resultados para la temperatura en el Punto
de Burbuja. Los cálculos para dicha tabla se muestran a continuación:
 Paso 1
En la columna (1) se listan las composiciones en fracción mol.
ni
x i=
100 %
5.0 % 30.0 %
x n−C H = =0.05 , x n−C H 10 = =0.30
3 8
100 % 4
100 %
40.0 % 25.0 %
x n−C H = =0.40 , x n−C H 14 = =0.25
5 12
100 % 6
100 %
 Paso 2
Como primera estimación del punto de burbuja, se escoge de la
monografía 60 ℃ ; en la columna (2) se listan las K correspondientes a
2
60 ℃ , 350 KN /m .
K n−C 3 H8 =4.70 , K n−C 4 H 10 =1.70 , K n−C 5 H 12 =0.62
 K n−C 6 H14 =0.25

 Paso 3
Se escoge el pentano como componente de referencia por ser el de mayor
K i ∙ xi
composición en la mezcla; en la columna (3) se listan las de la
Kj
ecuación 9.
K n−C H ∙ x n−C H 4.70 ∙ 0.05
3 8 3 8
= =0.379
K n−C H 5 12
0.62
K n−C H ∙ x n−C H 10 1.70 ∙0.30
4 10 4
= =0.822
K n−C H 5 12
0.62
K n−C H ∙ x n−C H 0.62 ∙ 0.40
5 12 5 12
= =0.400
K n−C H 5 12
0.62
K n−C H ∙ x n−C H 0.25 ∙0.25
6 14 6 14
= =0.1008
K n−C H 5 12
0.62
K i ∙ xi
∑ Kj
=¿ 0.379+0.822+0.400+0.1008=1.702 ¿

1
K n−C H12 = =0.588
5
1.702
K n−C 5 H12 teórica ≠ K n−C H estimada ;0.62 ≠0.588 5 12

El coeficiente de equilibrio del pentano teórico a la temperatura de 60 ℃

seleccionado difiere del coeficiente de equilibrio del pentano estimado.
Desde este punto se repiten los pasos anteriores hasta que dichos
coeficientes de equilibrio coincidan.

Repetición del proceso desde el paso 2

 Paso 2
El punto de partida es K n−C 5 H12 =0.588 correspondiente a una temperatura
de 56.8 ℃, a partir de dicha temperatura se toman de la monografía los
coeficiente de equilibrio K del resto de los componentes y se listan en la
columna (4).
K n−C 3 H8 =4.60 , K n−C 4 H 10 =1.60 , K n−C 5 H 12 =0.588
 K n−C 6 H14 =0.235

 Paso 3
Se mantiene el pentano como componente de referencia y los resultados
se listan en la columna (5).
K n−C H ∙ x n−C H 4.60 ∙ 0.05
3 8 3 8
= =0.391
K n−C H 5 12
0.588
K n−C H ∙ x n−C H 10 1.60 ∙0.30
4 10 4
= =0.816
K n−C H 5 12
0.588
K n−C H ∙ x n−C H 0.588 ∙ 0.40
5 12 5 12
= =0.400
K n−C H 5 12
0.588
K n−C H ∙ x n−C H 0.235 ∙0.25
6 14 6 14
= =0.100
K n−C H 5 12
0.588
K i ∙ xi
∑ Kj
=¿ 0.391+0.816+0.400+ 0.100=1.707 ¿

1
K n−C H12 = =0.586≈ 0.588
5
1.707
K n−C 5 H12 teórica=K n−C H estimada; 0.588=0.588 5 12

El coeficiente de equilibrio del pentano teórico a la temperatura de 56.8 ℃

seleccionado coincide con el coeficiente de equilibrio del pentano
estimado. Es decir que la temperatura en el punto de burbuja es de
56.8 ℃ .

 Paso 4
Este es el punto de burbuja; la composición correspondiente del vapor en
el punto de burbuja está dada en la columna (6). Los cálculos son
realizados con la ecuación 11.
α i , j ∙ xi
y j=
∑ αi , j ∙ x i
0.391 0.816
y n−C H = =0.229 , y n−C H 10 = =0.478
3 8
1.707 4
1.707
 0.400 0.100
y n−C H = =0.234 , y n−C H 14 = =0.0586
5 12
1.707 6
1.707
∑ y j=¿ 0.229+ 0.478+ 0.234+0.0586=1 ¿
Tabla 1: Punto de Burbuja
K i ∙ xi K i ∙ xi
xi K i @60 ℃ K i @56.8 ℃ yi
Componente Kj Kj
(1) (2) (4) (6)
(3) (5)
n−C3 H 8 0.05 4.70 0.379 4.60 0.391 0.229
n−C 4 H 10 0.30 1.70 0.822 1.60 0.816 0.478
n−C5 H 12 0.40 0.62 0.400 0.588 0.400 0.234
n−C6 H 14 0.25 0.25 0.1008 0.235 0.100 0.0586
1 1.702 1.707 1

Determinación del Punto de Rocío

Se establece una tabla con los datos y resultados para la temperatura en el
Punto de Rocío. El procedimiento es similar al punto de burbuja y muestran a
continuación:
 Paso 1
En la columna (1) se listan las composiciones en fracción mol.
ni
y i=
100 %
5.0 % 30.0 %
y n−C H = =0.05 , y n−C H = =0.30
3 8
100 % 4
100 %
10

40.0 % 25.0 %
y n−C H = =0.40 , y n−C H = =0.25
5 12
100 % 6
100 % 14

 Paso 2
 Como primera estimación del punto de rocío, se escoge de la monografía
80 ℃ ; en la columna (2) se listan las K correspondientes a
2
80 ℃ , 350 KN /m .
K n−C 3 H8 =6.30 , K n−C H =2.50 , K n−C
4 10 5 H12 =0.96
K n−C 6 H14 =0.43

 Paso 3
Se escoge el pentano como componente de referencia por ser el de mayor
composición en la mezcla; en la columna (3) se listan las volatilidades
relativa α ij de la ecuación 4.
Ki
α ij =
Kj
K n−C H8 6.30 K n−C H 10 2.50
α n−C H = 3
= =6.56 , α n−C H = 4
= =2.60
3 8
K n−C 5 H12 0.96 K n−C 4 10
5 H 12 0.96
K n−C H 0.96 K n−C H 0.43
α n−C H = 5 12
= =1.0 , α n−C H 14 = 6 14
= =0.448
5 12
K n−C H 5 12
0.96 6
K n−C H 5 12
0.96

 Paso 4
yi
Se listan en la columna (4) las de la ecuación 15.
α ij
y n−C H 0.05 y n−C H 0.30
3 8
= =0.00762 , = =0.1154
4 10

α n−C 3 H8 6.56 α n−C H 2.60 4 10

y n−C H 0.40 y n−C H 0.25

5 12
= =0.400 , = =0.558
6 14

α n−C 5 H 12 1.00 α n−C H 0.448 4 10

y
∑ α i =¿ 0.00762+0.1154 +0.400+0.558=1.081¿
ij

1.081≠ 0.96=K n−C 5 H12

El coeficiente de equilibrio del pentano a la temperatura de 80 ℃ difiere

del coeficiente de equilibrio del pentano estimado. Desde este punto se
repiten los pasos anteriores hasta que dichos coeficientes de equilibrio
coincidan.
 Repetición del proceso desde el paso 2
 Paso 2
El punto de partida es K n−C 5 H12 =1.08 correspondiente a una temperatura
de 83.7 ℃ , a partir de dicha temperatura se toman de la monografía los
coeficiente de equilibrio K del resto de los componentes y se listan en la
columna (5).
K n−C 3 H8 =6.60 , K n−C H =2.70 , K n−C
4 10 5 H12 =1.08
K n−C 6 H14 =0.47

 Paso 3
Se mantiene el pentano como componente de referencia; en la columna
(6) se listan las volatilidades relativa α ij de la ecuación 4.
Ki
α ij =
Kj
K n−C H8 6.60 K n−C H 2.70
α n−C H = 3
= =6.11 , α n−C H 10 = 4 10
= =2.50
3 8
K n−C 5 H12 1.08 4
K n−C H5 12
1.08
K n−C H 1.08 K n−C H 14 0.47
α n−C H = 5 12
= =1.0 ,α n−C H = 6
= =0.435
5 12
K n−C H 5 12
1.08 K n−C 6 14
5 H 12 1.08

 Paso 4
yi
Se listan en la columna (7) las de la ecuación 15.
α ij
y n−C H 0.05 y n−C H 10 0.30
3 8
= =0.00818 , 4
= =0.120
α n−C 3 H8 6.11 α n−C 4 H10 2.50
y n−C H 0.40 y n−C H 0.25
5 12
= =0.400 , = =0.575 6 14

α n−C 5 H 12 1.00 α n−C H 0.435 4 10

y
∑ α i =¿ 0.00818+ 0.120+0.400+0.575=1.103 ¿
ij

1.08 ≠1.103=K n−C H 5 12

El coeficiente de equilibrio del pentano a la temperatura de 83.7 ℃

nuevamente difiere del coeficiente de equilibrio del pentano estimado.
 Para no plantear el proceso repetitivo de estimación y hacer el ejercicio
más extenso se resume que coeficiente de equilibrio del pentano a partir
de 1.103 coincide con el próximo estimado, es decir, 1.103=1.103=K n−C 5 H12

y que representa una temperatura de 84 ℃ .

En consecuencia la temperatura y formación del punto de rocío es de

84 ℃ .

 Paso 5
Este es el punto de rocío; la composición correspondiente del líquido en el
punto de rocío está dada en la columna (8). Los cálculos son realizados
con la ecuación 11.
y i / αi , j
x j=
∑ y i /α i , j
0.00818 0.120
x n−C H = =0.0074 , x n−C H 10 = =0.1088
3 8
1.103 4
1.103
0.400 0.575
x n−C H = =0.3626 , x n−C H = =0.5213
5 12
1.103 6
1.707
14

∑ y j=¿ 0.0074 +0.1088+0.3626+ 0.5213=1 ¿

Tabla 2: Punto de Rocío

yi yi
Component yi K i @80 ℃ α ij K i @83.7 ℃ α ij xi
α ij α ij
e (1) (2) (3) (5) (6) (8)
(4) (7)
0.0076 0.0081 0.007
n−C3 H 8 0.05 6.30 6.56 6.60 6.11
2 8 4
 0.108
n−C 4 H 10 0.30 2.50 2.60 0.1154 2.70 2.50 0.120
8
0.362
n−C5 H 12 0.40 0.96 1.00 0.400 1.08 1.00 0.400
6
0.44 0.43 0.521
n−C6 H 14 0.25 0.43 0.558 0.47 0.575
8 5 3
1 1.081 1.103 1

Ejemplo 2. (Mezclas Ideales)

Se tiene una mezcla ideal con la siguiente composición: 20% mol n-Hexano, 50%
mol Benceno y 30% mol Tolueno.
Determinar:
a) La presión de Punto de Rocío de la mezcla a 150 ºF.
b) La temperatura de Punto de Rocío de la mezcla a 85 psia.
c) El Punto de Burbuja de la mezcla a 70 psig (Presión Barométrica 780
mmHg).
d) La presión en el Punto de Burbuja si la mezcla está a 150 ºF.

Solución:
a) Presión de Punto de Rocío de la mezcla a 150 ºF
La presión en el Punto de Rocío es la presión a la cual el primer vapor
comienza a condensar a 150 ºF. La mezcla se asume toda vapor y la
composición del condensado es determinada por:
n
1 yi
=∑ Sat
P T i=1 Pi

Al asumir que se tiene una mezcla ideal, se puede usar la Ecuación de

Antoine para determinar la presión de vapor de los componentes puros ( Pi ) :
Sat

B
ln ( P iSat )= A− T =150 ℉=338,6 K
C+T

Tabla 3
 yi
Comp. A B C Pi
Sat
yi Sat
Pi
−4
n-Hexano 15.8366 2697.55 -48.78 684.67 0.20 2.921 ×10
−3
Benceno 15.9008 2788.51 -52.36 472.93 0.50 1.057 ×10
−3
Tolueno 16.0137 3096.52 -53.67 171.75 0.30 1.747 ×10
−3
3.096 ×10

n
1 yi 1
=∑ Sat =3.096 ×10−3 → =3.096× 10−3
P T i=1 Pi PT
PT =333 mmHg

b) La temperatura en el Punto de Rocío de la mezcla a 85 psia

La Temperatura del Punto de Rocío es la temperatura a la cual el primer vapor
comienza a condensar a 85 psia. Para una solución ideal a 85 psia (4935
mmHg) la relación a usar es:
1
n
yi
=∑ Sat
P T i=1 Pi

Sat
P =exp A−
[ B
C +T ]

[ ]
0.20 0.50
+
1 [ 15.8366−
2697.55
−48.78+T ] e[ 15.9008−
2788.51
−52.36+T ]
= e
4935 +0.30

e
[ 16.0137−
3096.52
−53.67 +T ]

Se resolverá para T mediante el método de ensayo y error para obtener:

T=500K PT =13105 mmHg
T=425K PT =3842 mmHg
T=435K PT =4321 mmHg

T 435 K
c) La temperatura en el Punto de Burbuja de la mezcla a 70 psig (Presión
Barométrica 780 mmHg)
La Temperatura del Punto de Burbuja a una presión dada es la temperatura a
la cual la presión de vapor de la mezcla iguala la presión total, y en este
momento la primera parte del líquido se convierte en vapor. La mezcla se
asume toda líquida y la composición de vapor se asume en equilibrio con el
líquido. La relación usada bajo estas suposiciones para un líquido ideal es:

PT =∑ Pi ∙ x i
Sat

760 mmHg
70 psi∙ =3619 mmHg
14.696 psi

PT =( 780+3619 ) mmHg=4399mmHg|¿|

[ [ ]+0.50 ∙ e[ ]
]
2697.55 2788.51
15.8366− 15.9008−
− 48.78+T −52.36+T
4399= 0.20 ∙ e

+0.30 ∙ e
[ 16.0137−
3096.52
−53.67+ T ]

Al resolver por ensayo y error se obtiene:

T=400K PT =2372 mmHg
T=425K PT =4070 mmHg
T=430K PT =4496 mmHg

T 429 K

d) Presión de Punto de Burbuja de la mezcla a 150 ºF

La presión en el Punto de Rocío a 150 ºF (338.6K) es la presión a la cual el
primer líquido comienza a vaporizar cuando el vapor y el líquido están en
equilibrio. Para un líquido ideal:
PT =∑ Pi ∙ x i
Sat
 Tabla 4
Comp. A B C Pi
Sat
yi Sat
Pi ∙ x i
n-Hexano 15.8366 2697.55 -48.78 684.67 0.20 136.93
Benceno 15.9008 2788.51 -52.36 472.93 0.50 236.47
Tolueno 16.0137 3096.52 -53.67 171.75 0.30 51.53
424.9

PT =425 mmHg
Ejemplo 3.
Un depósito de vaporización instantánea se alimenta con una mezcla líquida
de 42 % molar de propano, y 58 % molar de etano. En el depósito se reduce
la presión y se separan corrientes de L y V. El sistema está a 305 K y ha de
vaporizar 40 % molar de líquido suministrado. ¿Cuál es la presión a la que
debe operar el depósito? ¿Cuáles son las composiciones en el equilibrio?

Solución:
T =305 K
P=? B . C .=100 Kgmol de Alimentación
z 1=0,58
z 2=0,42
Asumiendo
V =40 Kgmol

Podemos establecer los límites de la P, con el cálculo de Presión de burbuja y

cálculo de Presión de rocío.

Presión de Burbuja: se aplica la condición de la ecuación (3). Asumiendo la

alimentación como la composición de la fase líquida nos queda:
z Etano =x Etano =0,58
z Propano =x Propano =0,42

Asumimos valores de Presión de burbuja hasta que se cumpla la condición de la

ecuación (3).
Asumiendo P = 2800 kPa. Con la presión asumida y la T dada, se traza una línea
entre la P y T en el nomograma, que corte cada una de las líneas de las
constantes del etano y propano. Las constantes de equilibrio K se leen
directamente en el nomograma.

Tabla 5: P = 2800 kPa. Condición Burbuja

Componente Ki zi K i ∙ zi

Etano 1.5 0.58 0.87

Propano 0.52 0.42 0.2184

yi
K i= → y i=K i ∙ x i → y i=K i ∙ z i
xi
y Etano=K Etano ∙ z Etano =1.5 ∙0.58=0.87

y Propano=K Propano ∙ z Propano =0.52∙ 0.42=0.2184

∑ K i ∙ zi =0.87+0.2184=1.0884
La condición (3) no se cumple, se debe asumir un nuevo valor de P, que
disminuya la constante de equilibrio, para así disminuir la sumatoria. Observando
el nomograma, para que disminuya la constante de equilibrio debe aumentar la
presión. Asumiendo otro valor de P:

Tabla 6: P = 3000 kPa. Condición Burbuja

Componente Ki zi K i ∙ zi

Etano 1.38 0.58 0.8004

Propano 0.49 0.42 0.2058

yi
K i= → y i=K i ∙ x i → y i=K i ∙ z i
xi

y Etano=K Etano ∙ z Etano =1.38 ∙0.58=0.8004

y Propano=K Propano ∙ z Propano =0.49∙ 0.42=0.2058

∑ K i ∙ zi =0.8004+0.2058=1.0062
Como se cumple la condición, la presión de burbuja es 3000 kPa.
Presión de Rocío
Asumiendo las z i= y i tenemos que:
z Etano = y Etano =0,58
z Propano = y Propano=0,42

Asumimos valor de P, hasta que la condición de sumatoria de la ecuación (3) se

cumpla.
Tabla 7: P = 2500 kPa. Condición de Rocío

Componente Ki zi zi / K i

Etano 1.7 0.58 0.34

Propano 0.58 0.42 0.72

yi y z
K i= → xi = i → x i = i
xi Ki Ki
z Etano 0.58
x Etano = = =0.34
K Etano 1,7
z Propano 0.42
x Propano = = =0.72
K Propano 0.58
z
∑ Ki =0.34+ 0.72=1.06
i

La condición de sumatoria de la ecuación (3) no se cumple, se asume un nuevo

valor de P. Para disminuir la sumatoria, se debe incrementar el valor de Ki, para
que esto ocurra, se observa en nomograma, se debe disminuir la P.

Asumiendo un nuevo valor de 2300 kPa se tiene:

Tabla 8: P = 2300 kPa. Condición de Rocío

Componente Ki zi zi / K i

Etano 1.8 0.58 0.32

Propano 0.62 0.42 0.68

yi y z
K i= → xi = i → x i= i
xi Ki Ki
 z Etano 0.58
x Etano = = =0.32
K Etano 1,8
z Propano 0.42
x Propano = = =0.68
K Propano 0.62
z
∑ Ki =0.32+0.68=1.00
i

La presión de rocío es 2300 kPa.

La presión del sistema, es la incógnita del problema, debe estar entre la presión
de burbuja y la de rocío para que se vaporice el 40% de la alimentación.

Para determinar la P del sistema, se deben asumir valores de P, hasta que se

cumpla la siguiente ecuación de vaporización instantánea.
zi ∙ K i
1=∑
1+V ( K i−1 )

Asumiendo P = 2700 kPa

Se leen las constantes de equilibrio en el nomograma, a la presión de 2700 kPa y

la temperatura del sistema 31,85 °C.
Tabla 9: P = 2700 kPa. Vaporización

Componente Ki zi yi

Etano 1.51 0.58 0.727

Propano 0.52 0.42 0.270

Ya que:
z ∙K
∑ 1+V i K i−1 =
( i ) ( z Etano ∙ K Etano
1+V ( K Etano −1 )
z
)(
∙K
+ Propano Propano
1+V ( K Propano−1 ) )
( 1+ 0.40∙ ) ( 1+ 0.40∙ ( 0.52−1 ) )
z ∙K 0.58 ∙1.51 0.42 ∙0.52
∑ 1+V i( K i−1 ) = ( 1.51−1 )
+
i

z ∙K
∑ 1+V i( K i−1 ) =0,998
i
La presión del sistema es de 2700 kpa.

Ejemplo 4. (Coeficiente de Fugacidad y Coeficiente de Actividad)

La fugacidad del componente 1 en una mezcla líquida binaria a 298 K y 20 bar
viene dada por:
2 3
f 1=50 x1 −80 x1 +40 x 1

Donde la fugacidad esta en bar y x 1 es la fracción molar. Determine:

a) La fugacidad del componente puro 1.
b) El coeficiente de fugacidad ∅ .
c) La constante de Henry K 1
d) El coeficiente de actividad γ .

Solución:
a) Cuando x 1=1 entonces:
2 3
f 1=50 ( 1 )−80 ( 1 ) + 40 ( 1 )
f 1=10 ¯
¿
b)
f1 ¯
∅= =20 ¯ ¿ ¿
P 10 ¿ =0.5 ¿
c) La constante de Henry será
^
f1
K 1=lim
x i→ 0 x1

K 1=lim ( 50−80 x 1+ 40 x1 ) =50 ¯¿

2

xi→ 0

d) .
2 3
f1 50 x 1−80 x 1 + 40 x 1 2
γ 1= = =5−8 x 1+ 4 x 1
x1 ∙ f 1 10 x 1

γ 1=5−8 x 1+ 4 x 12=5−8 ( 1 ) + 4 ( 1 )2
γ 1=1
 CONCLUSIÓN
El equilibrio de fases líquido vapor ELV es uno de los temas de interés en la
Ingeniería, puesto que en muchos procesos industriales se presenta el hecho de
tener que realizar operaciones unitarias de separación que dependen de las
características del equilibrio en el sistema.

En este sentido, se pueden hacer estimaciones del equilibrio, evaluar las

composiciones que se prevén para la fase gaseosa a partir de la información dada
para la fase líquida, utilizando correlaciones matemáticas. Para sistemas de
solución ideal, la ecuación de Raoult es una excelente herramienta para
establecer esos valores.

Sin embargo, las mezclas en solución no siempre presentan comportamiento

ideal. Para sistemas de mezcla con comportamiento no ideal es preciso introducir
factores de corrección, conocidos como coeficientes de actividades para la fase
líquida y coeficientes de fugacidad para la fase gaseosa, y que fueron
desarrollados en los apartados de esta investigación.

Para culminar, se plantearon problemas típicos del equilibrio liquido vapor para
sistemas ideales o no ideales, teniendo una evaluación de cada uno de los casos
planteados.
 BIBLIOGRAFÍA

 Goitia, K. (2010). Equilibrio de fases. Universidad Nacional Experimental

Francisco de Miranda. Punto Fijo. Venezuela.

 Pimente, I. (2013). Selección de modelos termodinámicos Liquido Vapor.

Universidad Nacional Experimental de las Fuerzas Armadas. Maracaibo.
Venezuela.

 Treybal, R. (1994). Operaciones de transferencia de masa. Ed. McGraw-

Hill. 2da Edición.

 Salazar, F. (2014). Determinación de puntos de burbuja y de rocío.

Universidad Rafael Landívar. Guatemala.