Tr9_3-V2.e-McWill.

xmcd, 18/03/2023 IQF1017 Procesos de Separación III, Tema 1 de 5

Problema 9.3 (Treybal)

Una solución tiene la siguiente composición, expresada como fracción mol: etano, 0.0025; propano, 0.25;
i-butano, 0.185;n-butano, 0.560; i-pentano 0.0025. A continuación, la presión es de 10 bars (145 lbf/in2) abs. Utilícense los
coeficientes de distribución en el equilibrio del nomograma de Depriester.
a) Calcule el punto de burbuja.
b) Calcule el punto de rocío.
c) La solución se va a evaporar instantáneamente con el fin de evaporar 40% en mol de la alimentación. Calcule la
composición de los productos.
d) La solución se destila diferencialmente para producir un residuo que contenga 0.80 fracción mol de n-butano. Calcule la
composición total del residuo y el porcentaje de la alimentación que se evapora.

Datos: Presión: p := 10 bar

Número de c := 5
componentes:

Nota: en los siguientes vectores introducir los nombres de los componentes, su fracción molar y las constantes de la ec.
de McWilliams, sin modificar el primer elemento de cada vector

 "Componente"   0 
 "Etano"   0.0025 
   
 "Propano"   0.25 
C := z :=
 "i-butano"   0.185 
 "n-butano"   0.560 
   
 "i-pentano"   0.0025 

Ecuación de McWilliams: La ecuación se ha adaptado parta usar temperaturas en °K o °R y presión en cualquier unidad

( )
Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p := "Expresamos T en °R y p en psia"
T
T
R
p
p
psi
 aT1 aT2 ap2 ap3 
Kt  exp + + aT6 + ap1 ln( p ) + + 
 T2 T 2 p 
 p 

De Bibliografía (Wankat) tenemos las constantes de la ec de McWilliams, para lo s componentes partic ulares:

 "aT1"   "aT2"   "aT6"   "aT2"   "aT2"   "aT2" 

 -687248.25   0   7.90694   -0.88600   49.02654   0 
           
 -970688.5625   0   7.15059   -0.76984  a :=  0   6.90224 
aT1 := a := aT6 := a := ap3 :=
 -1166846  T2  0   7.72668  p1  -0.92213  p2  0   0 
 -1280557   0   7.94986   -0.96455   0   0 
           
 -1481583   0   7.58071   -0.93159   0   0 
Autor: M.A.M.N.
Tr9_3-V2.e-McWill.xmcd, 18/03/2023 IQF1017 Procesos de Separación III, Tema 2 de 5

SOLUCIÓN:
a) Punto de burbuja
Dado que se proporciona la presión, calculamos los coeficientes de distribución de equilibrio proponiendo una
temperatura y verificando que la suma de las concentraciones de la fase V sea igual con uno

Tomando x := z

Suponemos T := 331.2 K

i := 1 , 2 .. c Ke := Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p 

i i  i i i i i  (1)

y := Ke  x
i i i (2)

Componente Coeficiente de distribución

 "Componente"   0   0   0 
 "Etano"   4.787   0.0025   0.012 
       
 "Propano"   1.888   0.25   0.472 
C= Ke = x= y=
 "i-butano"   0.864   0.185   0.1599 
 "n-butano"   0.635   0.56   0.3557 
       
 "i-pentano"   0.294   0.0025   0.0007 
Verificamos la suma de las
fracciones molares de la fase V
y = 1 (3)

Temperatura de burbuja: T = 331.2 K (iterar hasta converger)

Una forma más sistemática se obtiene analizando el procedimiento anterior. Esto conduce a un problema de
raíces de ecuaciones donde la función a resolver se obtiene sustituyendo ec. (1) en la ec. (2) y el resultado a su vez en la
ec. (3). Este ultimo resultado se iguala a cero obteniendo finalmente:

c
F( T) :=
  Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p  xi - 1
  i i i i i i  
(4)
i =1
Graficamos la función para ubicar la raíz que sería la temperatura de burbuja
T := 290  K , 291  K .. 410  K

1 La raiz de la ecuación (4) se

ubica en la intersección con
F ( T) el eje horizontal.

300 350 400

-1

Autor: M.A.M.N.
Tr9_3-V2.e-McWill.xmcd, 18/03/2023 IQF1017 Procesos de Separación III, Tema 3 de 5

Escogemos un valor cercano a la raiz y resolvemos la ecuación (4)

Suponemos Tini := 300  K

(( )
Tb := root F Tini , Tini )
Temperatura de burbuja: Tb = 331.184 K

a) Punto de rocío
Dado que se proporciona la presión, calculamos los coeficientes de distribución de equilibrio proponiendo una
temperatura y verificando que la suma de las concentraciones de la fase L sea ig ual con uno

Tomando: y := z

Suponemos T := 341.9 K

i := 1 , 2 .. c Ke := Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p 

i i i  i i i i 
y
i
x :=
i Ke
i

Componente Coeficiente de distribución

 "Componente"   0   0   0 
 "Etano"   5.393   0.0005   0.0025 
       
 "Propano"   2.234   0.1119   0.25 
C= Ke = x= y=
 "i-butano"   1.058   0.1749   0.185 
 "n-butano"   0.793   0.7061   0.56 
       
 "i-pentano"   0.38   0.0066   0.0025 
Verificamos la suma de las
fracciones molares de la fase V
x = 1 y = 1
Temperatura de rocío: T = 341.9 K (iterar hasta converger)

También puede implementarse un procedimiento similar al del cálculo del punto de burbuja. En este caso, la función a
resolver sería:

c y
i 
F( T) :=
  -1
 Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p  
(5)
i =1   i i i i i i 

Graficamos la función para ubicar la raíz que sería la temperatura de rocío

T := 290  K , 291  K .. 410  K

Autor: M.A.M.N.
Tr9_3-V2.e-McWill.xmcd, 18/03/2023 IQF1017 Procesos de Separación III, Tema 4 de 5

La raiz de la ecuación (5) se

1 ubica en la intersección con
el eje horizontal.
F ( T)

300 350 400

-1

Escogemos un valor cercano a la raiz y resolvemos la ecuación (5). El valor de temperatura deberá estar entre la de
burbuja y la de rocío
Suponemos Tini := 300  K

(( )
Td := root F Tini , Tini )
Temperatura de rocío: Td = 341.894 K

El cálculo de los coeficientes de distribución de equilibrio puede efectuarse por algún otro modelo termodinámico que sea
aplicable a la mezcla

c) La solución se va a evaporar instantáneamente con el fin de evaporar 40% en mol de la

alimentación. Calcule la composición de los productos.

En este caso la fracción vaporizada es ψ := 0.40

Haremos uso de la ecuación de Rachford-Rice, sin embargo la incógnita no es la fracción de vapor, ψ, sino la temperatura
de la destilación instantánea.

Procedemos por tanteo suponiendo una temperatura del flash e iterando hasta que la fracción de vapor converja al valor del
enunciado.

Suponemos T := 336.71 K

Calculamos los coeficientes de distribución de equilibrio:

i := 1 , 2 .. c Ke := Keq aT1 , aT2 , aT6 , ap1 , ap2 , ap3 , T , p 

i 
i i i i i i 

Resolvemos la ecuación de Rachford-Rice para determinar la fracción de vapor ψ

c  zi  1 - Ke  
  i 
f ( ψ) :=
  1 +  Ke - 1  ψ
(6)
i =1   i  

El valor de la fracción de vapor debe estar en el inervalo [0,1] y deberá coincidir con el dado en el enunciado. En caso
contrario, deberemos cambiar la suposicipon de la temperarura
Autor: M.A.M.N.
Tr9_3-V2.e-McWill.xmcd, 18/03/2023 IQF1017 Procesos de Separación III, Tema 5 de 5

Graficamos la ecuación de Rachford-Rice

0.2

0.1 La raiz de la ecuación (6) se

ubica en la intersección con
el eje horizontal.
f ( ψ)
0 0.2 0.4 0.6 0.8

- 0.1

- 0.2

Tomamos como valor iniciañ ψini := 0.5

(( )
ψ := root f ψini , ψini )
(Verificamos que este valor coincida con el del
ψ = 0.4
enunciado e iteramos de ser necesario
cambiando el valor propuesto de la temperatura)

Composición de los productos:

Una vez lograda la convergencia, tanto el valor de ψ como los valores de las constantes de disrtribución de equilibrio, Ke,
son los correctos y podemos usarlos para el cálculo de las composiciones de los productos

z  Ke z
i i i (7) y (8)
i := 1 , 2 .. c y := x :=
i 1 +  Ke - 1   ψ i 1 +  Ke - 1   ψ
 i   i 

x = y =
 "Componente"  i i
 "Etano"  0.0009 0.0048
 
"Propano"  0.1754 0.3618
C= 
 "i-butano"  0.1879 0.1807
 "n-butano"  0.6323 0.4515
 
0.0034 0.0011
 "i-pentano" 

x = 1 y = 1

Como ejercicio para el estudiante, se deja el elaborar un procedimiento que realice automáticamente la solución de este
inciso de igual forma que se hizo para los incisos (a) y (b)

También puede resolverse el problema usando otros modelos para el cálculo de los coeficientes de distribución de
equilibrio, Ke, tal como el modelo ideal (Ley de Raoult y Ec de Antoine para las presiones de vapor).

Autor: M.A.M.N.