Ejercicios Resueltos de Destilación Completo
Ejercicios Resueltos de Destilación Completo
Ejercicios Resueltos de Destilación Completo
LEY DE RAOULT Para soluciones concentradas podemos aplicar la Ley de Raoult que es
una ley ideal y define para fases vapor – líquido en equilibrio:
Donde 𝑃̅𝑖 es la presión parcial del componente más volátil en el vapor [atm], 𝑃𝑖 es la presión
de vapor del componente i puro [atm] y x es la fracción mol del componente más volátil en
el líquido.
Esta ley sólo es válida para soluciones ideales
Fracción parcial del componente gaseoso
Si:
Volatilidad Relativa:
Para 98.40 ºC
760 𝑚𝑚𝐻𝑔 − 333 𝑚𝑚𝐻𝑔
𝑥𝐴 = =1
760 𝑚𝑚𝐻𝑔 − 333 𝑚𝑚𝐻𝑔
760 760 𝑚𝑚𝐻𝑔 − 333 𝑚𝑚𝐻𝑔
𝑦𝐴 = ∗ =1
760 760 𝑚𝑚𝐻𝑔 − 333 𝑚𝑚𝐻𝑔
Datos de Equilibrio a partir de las presiones de vapor de ambos componentes
P n-heptano P n- octano
T ºC Xa Ya α
[mmHg] [mmHg]
98,40 760 333 1,0000 1,00 2,2823
105,00 940 417 0,6558 0,8112 2,2542
110,00 1050 484 0,4876 0,6737 2,1694
115,00 1200 561 0,3114 0,4917 2,1390
120,00 1350 650 0,1571 0,2791 2,0769
125,60 1540 760 0,0000 0,0000 2,0263
2,1580
b) Tomando en cuenta la ecuación:
D = (40 % F)
YD
𝑆𝑖 𝐷 = 0.4 ∗ 100 = 40 𝑚𝑜𝑙𝑒𝑠
XB
Si tenemos en cuenta en la gráfica de x vs y
De la gráfica:
𝑥𝐵 = 0.28; 𝑦𝐷 = 0.454
Interpolando XB = 0.28 T = 116.017 ºC
115 0,3114
116,017498 0,28
120 0,1571
Inciso e)
ESQUEMA DEL PROBLEMA
D = (20 % F)
YD
𝑆𝑖 𝐷 = 0.2 ∗ 100 = 20 𝑚𝑜𝑙𝑒𝑠
XB
Xn-Heptano= 0.35
Xn-Octano = 0.65
𝑚𝑛−𝐻𝑒𝑝𝑡𝑎𝑛𝑜 = 100 ∗ 0.35 = 35 𝑔
𝑚𝑛−𝑂𝑐𝑡𝑎𝑛𝑜 = 100 ∗ 0.65 = 65 𝑔
35
𝑥𝑛−ℎ𝑒𝑝𝑡𝑎𝑛𝑜 = 100.204 = 0.38
35 65
100.204 + 114.231
110 0,4876
114.18842 0,34
115 0,3114
Obteniendo los datos de las fracciones molares del componente más volátil en las fases líquido y
vapor
T Xa Ya
YD
F=100 Kmol
XF = 0.2
XB
Balance global en la destilación cerrada:
𝐹 =𝐷+𝐵
100 = 𝐷 + 𝐵
Balance por el componente más volátil:
𝐹 ∗ 𝑥𝐹 = 𝐷 ∗ 𝑦𝐷 + 𝐵 ∗ 𝑥𝐵
(𝐷 + 𝐵) ∗ 𝑥𝐹 = 𝐷 ∗ 𝑦𝐷 + 𝐵 ∗ 𝑥𝐵
𝐷 ∗ 𝑥𝐹 + 𝐵 ∗ 𝑥𝐹 = 𝐷 ∗ 𝑦𝐷 + 𝐵 ∗ 𝑥𝐵
𝐵 𝑦𝐷 − 𝑥𝐹
=
𝐷 𝑥𝐹 − 𝑥𝐵
Hallando los valores de 𝑦𝐷 𝑦 𝑥𝐵
Interpolando el valor de la x e y en la tabla con la temperatura de operación 72 ºC
73 0,10
72 0,125
71 0,15
𝑥𝐵 = 0.125
73 0,262
72 0,305
71 0,348
𝑦𝐷 = 0.305
En la ecuación:
𝐵 0.305 − 0.20
= = 1.4
𝐷 0.20 − 0.125
𝐵 = 1.4 ∗ 𝐷
100 = 𝐷 + 1.4 ∗ 𝐷
𝐷 = 41.67 𝐾𝑚𝑜𝑙
𝐵 = 58.33 𝐾𝑚𝑜𝑙
41.67 58.33
𝒇= = 0.4167 𝝋= = 0.5833
100 100
𝐷 ∗ 𝑦𝐷 41.67 ∗ 0.305
𝑅= ∗ 100% = = 63.55 %
𝐹 ∗ 𝑥𝐹 100 ∗ 0.2
X Y
0,05 0,107
0,1 0,202
0,15 0,29
0,2 0,369
0,25 0,442
0,3 0,506
0,35 0,565
0,4 0,618
0,45 0,665
0,5 0,708
0,55 0,748
0,6 0,784
0,65 0,817
0,7 0,847
0,75 0,875
0,8 0,9
0,85 0,926
0,9 0,95
1 1
BASE DE CALCULO DE 100 MOLES
Y = 0.4 fm
F = 100 moles
Xf=0.5 (𝑦1 − 𝑦0 )
𝑦 = 𝑦0 + ∗ (𝑥 − 𝑥0 )
(𝑥1 − 𝑥0 )
0,2 0,369
0.22 0.4
0,25 0,442
𝑥𝐵 = 0.22
Para la Grafica
Grafica Ej 3
25
15
10
0 xF
0 0,1 0,2 xB 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
1
𝑠𝑖 → 𝑟𝑒𝑝𝑟𝑒𝑠𝑒𝑛𝑡𝑎 𝑙𝑎 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑟𝑣𝑎 𝑞𝑢𝑒 𝑠𝑖𝑔𝑢𝑒 𝑢𝑛𝑎 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑐𝑎
𝑦−𝑥
1
= 𝑓(𝑥) = 𝑎𝑥 5 + 𝑏𝑥 4 + 𝑐𝑥 3 + 𝑑𝑥 2 + 𝑒𝑥 + 𝑓
𝑦−𝑥
Si tomamos una ecuación de cuarto grado para representar el área entre los limites definidos:
𝑥𝐹
𝐹
𝐿𝑛 = ∫ (𝑎𝑥 4 + 𝑏𝑥 3 + 𝑐𝑥 2 + 𝑑𝑥 + 𝑒)𝑑𝑥
𝐵
𝑥𝐵
𝛼𝐴𝐵 = 2.16
a) Resolver con la Ecuación de Rayleigh
b) Resolver tomando la volatilidad del sistema como constante en el intervalo de
operación:
F= F=D+B
XF =0.4 Si tomamos una base de calculo
De cálculo de F =100 moles
D = 60 moles
B = 100 – 60 = 40 moles
2.16 ∗ 𝑥
𝑦=
1 + 1.16 ∗ 𝑥
x y 1/(y-x)
0,1 0,194 10,690
0,2 0,351 6,638
0,3 0,481 5,534
0,4 0,590 5,259
0,5 0,684 5,448
0,6 0,764 6,092
0,7 0,834 7,438
0,8 0,896 10,388
0,9 0,951 19,579
1 1,0000
En la Gráfica:
GRAFICO
25,000
15,000
10,000
5,000
0,000
0 0,1 0,2 0,3 0,4 x 0,5 0,6 0,7 0,8 0,9 1
XB=? F
𝑥𝐹
𝐹 𝑑𝑥
ln = ∫
𝐵 𝑥𝐵 𝑦 − 𝑥
a 360.56
b -665.61
c 459.06
d -138.11
e 20.582
0.4
100
𝐿𝑛 = ∫ (360.56𝑥 4 − 665.61𝑥 3 + 459.06𝑥 2 − 138.11𝑥 + 20.582)𝑑𝑥
40
𝑋𝐵
0.4
F= F=D+B
XF =0.4 Si tomamos una base de calculo
De cálculo de F =100 moles
D = 60 moles
B = 100 – 60 = 40 moles
Ecuación:
Entonces:
1.160.4 1 − 𝑥𝐵 2.16
2.5 = √ ∗[ ]
𝑥𝐵 1 − 0.4
Resolviendo la ecuación no Lineal:
𝑥𝐵 = 0.2341
Sección de empobrecimiento:
Columnas de Agotamiento:
Para el cálculo de estas columnas se emplean las ecuaciones deducidas de los balances de
materia y la ecuación relativa de la recta inferior de operación. Para trazar esta recta en el
diagrama x – y tenemos en cuenta que ha de pasar por el punto (xW, xW) y podemos
determinar otro punto de esta recta para y = 1, ya que para este punto el valor de x vendrá
dado por:
x 0.00 0.039 0.065 0.09 0.140 0.215 0.430 0.601 0.782 0.853 0.932 0.985
y 0.00 0.085 0.153 0.225 0.316 0.449 0.678 0.802 0.910 0.948 0.970 0.993
Datos:
D =?
XD = 0.975 RD = 1.4 Rm
Entrada de líquido en la alimentación:
F = 40 Kmol/h
XF = 0.55 Temperatura de punto de ebullición
B =?
XB = 0.02
Pendiente = Infinito
a) Flujos molares B y D
Balance Global en el sistema
𝐹 =𝐵+𝐷
40 = 𝐵 + 𝐷 ………………… (1)
Balance global por componente volátil
𝐹 ∗ 𝑥𝐹 = 𝐵 ∗ 𝑥𝐵 + 𝐷 ∗ 𝑥𝐷
40 ∗ 0.55 = 𝐵 ∗ 0.02 + 𝐷 ∗ 0.975 … … … … (2)
sistema de ecuaciones con la ecuación 1 y 2
40 = 𝐵 + 𝐷
22 = 0.02𝐵 + 0.975𝐷
Resolviendo el sistema de ecuaciones
𝐾𝑚𝑜𝑙 𝐾𝑚𝑜𝑙
𝐵 = 17.80 𝐷 = 22.20
ℎ ℎ
Puntos en la Grafica
𝑅𝐷 = 1.4 𝑅𝑚
La pendiente minina en relación al reflujo mínimo viene por la ecuación:
𝑅𝑚 𝑋𝐷 − 𝑦 ′
𝑚𝑚 = =
𝑅𝑚 + 1 𝑋𝐷 − 𝑥 ′
Mediante la Grafica
𝑥 ′ = 0.55 ; 𝑦 ′ = 0.77
En la ecuación
𝑅𝑚 0.975 − 0.77
= → 𝑅𝑚 = 0.93
𝑅𝑚 + 1 0.975 − 0.55
𝑅𝐷 = 1.4 ∗ 0.93 = 1.302
En la zona de rectificación, la nueva ordenada en el punto de origen para el reflujo real está dada por:
𝑥𝐷 0.975
0.00 = = = 0.423 → 𝑛𝑢𝑒𝑣𝑎 𝑐𝑜𝑜𝑟𝑑𝑒𝑛𝑎𝑑𝑎 𝑑𝑒 𝑜𝑟𝑖𝑔𝑒𝑛 𝑒𝑛 𝑒𝑙 𝑒𝑗𝑒 𝑑𝑒 𝑦
1 + 𝑅𝐷 1 + 1.302
Si unimos el nuevo punto de origen con XD se tendrá la recta en la parte de enriquecimiento:
x 0.00 0.04 0.08 0.10 0.15 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
y 0.00 0.23 0.365 0.418 0.517 0.579 0.665 0.729 0.779 0.825 0.870 0.915 0.958 1.00
D=
XD = 0.95
F = 100 Kmol/h
XF = 0.4
La alimentación ingresara en la columna como
liquido saturado:
𝑙𝑖𝑞𝑢𝑖𝑑𝑜 𝑎 𝑙𝑎 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑒𝑏𝑢𝑙𝑙𝑖𝑐𝑖ó𝑛
B=
XB =0.02 𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = 𝑖𝑛𝑓𝑖𝑛𝑖𝑡𝑜
𝑅𝐷 = ?
Puntos en la Grafica
Método de solución:
Es necesario usar un procedimiento de intento y error (tantear gráficamente) sobre el diagrama x vs
y, Se trazan las líneas de operación pertenecientes a las rectas de enriquecimiento y agotamiento que
partan de los puntos de composición xD y xB sobre la diagonal a 45º y que intersecten sobre la línea
de alimentación. El par de líneas de operación correcto es el que da 7 etapas de equilibrio (incluyendo
el caldero) en la construcción del escalonamiento.
Se puede observar con el tanteo de las rectas se obtuvo 6.8 etapas que serían 7 etapas exactas
incluyendo el rehervidor
De la recta de enriquecimiento se puede obtener el valor de la intersección en el eje Y que seria:
𝑥𝐷 0.95
0.00 = = = 0.35
1 + 𝑅𝐷 1 + 𝑅𝐷
Despejando el valor de reflujo del destilado
𝑹𝑫 = 𝟏. 𝟕𝟏𝟒
𝐹 =𝐷+𝐵
D =50 Kmol/h
XD = 0.95 𝐵 = 100 − 50 = 50 𝐾𝑚𝑜𝑙/ℎ
F = 100 Kmol/h
XF = 0.5
Balance por el componente más Volátil:
𝐹 ∗ 𝑥𝐹 = 𝐵 ∗ 𝑥𝐵 + 𝐷 ∗ 𝑦𝐷
Puntos en la Grafica
Según la gráfica el número de platos mínimos es igual a
El número mínimo de etapas es 3 etapas más 1 caldero
Si la pendiente de la recta de alimentación tiene un porcentaje de vapor del 60%
(𝟏 − 𝒇) (𝟏 − 𝟎. 𝟔)
𝑚=− =− = −𝟎. 𝟔𝟔𝟕
𝒇 𝟎. 𝟔
Angulo = -33.7º
En la Grafica ubicamos la recta de alimentación
a) Aplicando la relación de reflujo mínimos y las siguientes ecuaciones:
𝑅𝐷 = 1.5 𝑅𝑚
La pendiente minina en relación al reflujo mínimo viene por la ecuación:
𝑅𝑚 𝑋𝐷 − 𝑦 ′
𝑚𝑚 = =
𝑅𝑚 + 1 𝑋𝐷 − 𝑥 ′
En la ecuación
𝑅𝑚 0.95 − 0.652
= → 𝑅𝑚 = 0.77
𝑅𝑚 + 1 095 − 0.265
𝑅𝐷 = 1.5 ∗ 0.77 = 1.155
En la zona de rectificación, la nueva ordenada en el punto de origen para el reflujo real está
dada por:
𝑥𝐷 0.95
0.00 = = = 0.441 → 𝑛𝑢𝑒𝑣𝑎 𝑐𝑜𝑜𝑟𝑑𝑒𝑛𝑎𝑑𝑎 𝑑𝑒 𝑜𝑟𝑖𝑔𝑒𝑛 𝑒𝑛 𝑒𝑙 𝑒𝑗𝑒 𝑑𝑒 𝑦
1 + 𝑅𝐷 1 + 1.155
F=
79 % en Vol. de N
D=?
YD= 0.9
Balance global
𝐹 =𝐷+𝐵
El flujo B será representado por W
𝐹 =𝐷+𝑊