Destilacion Ponchon y Savarit y Mccabe y Thiele FINAL
Destilacion Ponchon y Savarit y Mccabe y Thiele FINAL
Destilacion Ponchon y Savarit y Mccabe y Thiele FINAL
& THIELE.
2010
Problema: a continuacin se dan los datos del equilibrio Vapor Lquido y entalpa
para el sistema heptano - etilbenceno, a una (1) atm abs de presin. Utilice los dos
mtodos referenciados arriba; recuerde que ambos son complementarios. Establecer
diferencias.
t, C xH yH
136.2 0 0
129.5 0.08 0.233
122.9 0.185 0.428
119.7 0,251 0.514
116.0 0.335 0.608
110.8 0.487 0.729
106.2 0.651 0.834
103.0 0.788 0.904
100.2 0.914 0.963
98.5 1.000 1.000
Destilacin segn los mtodos de Ponchon & Savarit y McCabe & Thiele.
1. Cmo cree que se obtuvieron los datos de las entalpas HL y HG? De cules
ecuaciones se obtuvieron? Cmo se les denominan a estas entalpas?
Observe bien los datos del enunciado del problema propuesto # 9.4 del texto
de Treybal.
2
Esta ecuacin consta de dos partes, la primera parte corresponde al calor sensible de
la mezcla y la segunda corresponde al calor de mezclado de los componentes.
Donde:
CpL= capacidad calorfica de la solucin a TL.
TL = temperatura del liquido en el punto de formacin de burbuja correspondiente a la
concentracin del liquido a la presin predominante.
T0 = temperatura de referencia.
MA= peso molecular promedio de la mezcla.
Hdis= calor de disolucin a la temperatura TL (punto de formacin de burbuja).
Por otro lado tenemos q HG representa la entalpa del vapor saturado y se halla por
medio de la siguiente ecuacin
Donde:
Y= fraccin mol del vapor en el equilibrio del componente mas voltil en el vapor q en
nuestro caso es el heptano.
CLA, CLB= capacidades calorficas del liquido de cada componente.
TG= temperatura del vapor en el punto de formacin de roco.
T0= temperatura de referencia.
A, B = calores latentes de evaporacin de las sustancias puras a TG.
NOTA: cabe recordar que A representa el componente mas voltil es decir el heptano
y B el menos voltil es decir el etilbenceno.
Hdis= Lo llaman calor integral de solucin esta hallado a 15C y las unidades son [KJ/
Kmol sol.], a medida que aumenta la fraccin mol en el liquido del componente mas
voltil (acetona) aumenta el calor desprendido, llega a un mximo en 0.3 fraccin mol
y luego empieza a disminuir.
CL= capacidad calorfica del liquido a 17.2C, las unidades son KJ/ (Kg. sol* C),
disminuye su valor al aumentar la temperatura del liquido.
3
El procedimiento grafico inicia teniendo muy claro las entalpas del liquido-vapor y
representndolas en un grafico H vs. x, y, despus procedemos a encontrar las lneas
de unin en este grafico esto se logra haciendo coincidir la entalpa del liquido en el
equilibrio HL con la entalpa del vapor en el equilibrio HG a la misma temperatura.
Para hallar la curva de equilibrio en el grafico x vs. y, se traza en este diagrama una
lnea de 45 y=x, que tenga como intercepto el punto (0,0) y luego se dispone los dos
diagramas seguido el grafico H vs. x, y en la parte superior y el grafico x, y en la parte
inferior. Despus proyectamos los puntos de saturacin, correspondientes a las lneas
de reparto, con un lnea vertical del primer diagrama hasta que corte la lnea y=x,
desde el primer HL y luego se traza otra vertical desde el HG hasta la recta de 45
cuando se tienen estas dos lneas ubicados en la recta de 45, se procede a tomar el
punto resultante de la proyeccin de HG y se traza una lnea horizontal que corte la
vertical que resulto como proyeccin de HL, el punto resultante, es el primer punto de
la curva de equilibrio. Este procedimiento se sigue realizando para los valores de HG y
HL que correspondan a las lneas de unin, y as se obtienen los puntos
correspondientes a la curva de equilibrio en X vs. Y.
Se escoge un punto aleatoriamente entre las dos fases este representa q hay presente
liquido-vapor, por ejemplo el punto propuesto (0.675, 37500 KJ/ Kmol), la lnea de
unin para ese punto se determina de la siguiente manera: se toma la lnea de unin
mas cercana al punto y se traza una paralela a esta lnea que toque el punto, el punto
correspondiente a esta lnea HL, lo proyectamos hasta que corte la curva de equilibrio
y luego trazamos una lnea horizontal hasta la recta Y= X, este punto lo proyectamos
subiendo verticalmente hacia la grafica H vs. X, Y, en donde corte la regin de vapor
saturado obteniendo un HG, despus unimos los dos puntos HG y HL, si esta lnea de
reparto pasa por el punto propuesto, encontramos la lnea de unin. En caso contrario
repetimos el procedimiento para un nuevo punto de HL, cercano a la lnea de reparto
que ms se aproxima al punto propuesto, hasta que al unir los puntos HG y HL la lnea
formada pase por dicho punto.
4
L M V
Balance de masa:
6000 0.5
ZF 100 0.5147
58.3019
6000 Kg
( PM AV ) MEZCLA 102.9126
58.3019 Kmol
5
1KmolH
95Kg
ZD 100 Kg 0.95 0.9526
1KmolH 1KmolEB 0.9971
95Kg 5Kg
100 Kg 106 Kg
100 Kg
( PM AV ) destilado 100.2908
0.9971 Kmol
1KmolEB
5Kg
0.0467 0.05
106 Kg
XW
1KmolEB 1KmolH 0.9434
5Kg 95Kg
106 Kg 100 Kg
100 Kg
( PM AV ) colas 105.6829
0.9462 Kmol
BALANCE DE MASA
F ZF D ZD W XW
58.3019(0.5147) D(0.9526) W (0.05)
6
El flujo molar de colas: W= 28.2854 Kmol/h.
Relacin de Reflujo
L0 G1 D
R 1.5
D L0 G1
G1 D 1.5L0 G1
G1 L0 D
G1 RD D D( R 1) 30.0165(1.3978 1)
Kmol
G0 71.9745
h
=
Donde;
Lo: es el flujo de regreso al fraccionador.
D: Es el flujo del destilado que sale de la torre.
Para reflujo total no se entrega destilado, sino que todo el gas que sale de la torre se
condensa y se devuelve como reflujo, por tanto D tiende a cero y R tiende a infinito,
adems la relacin Ln/Gn+1 es igual a la relacin ****** que a su vez es igual 1.Por
tanto las curvas de operacin coinciden con la linea de 45 en la grfica X vs Y, y los
puntos D Y W en la grfica X,Y vrs H estan en el infinito, en consecuencia el
nmero de platos requeridos es el valor mnimo Nm.
7
punto donde XD y la recta de 45 se unen, esa horizontal se traza hasta donde corte la
curva de equilibrio y en ese punto se lanza una vertical hasta que corte la recta de 45,
all el triangulo que se forma de esta vertical,esta horizontal y la linea de 45 implica la
primera etapa. El resto de etapas se pueden hallar de la misma manera que e hallo la
primera pero empezando donde se termin el triangulo anterior.
Este proceso anterior se puede observar en las hojas anexas donde se observa que
grficamente hay 5.2 etapas para el mtodo de Ponchon y Savarit, que en total son 4
platos. Por el mtodo de Mccbe y Thiele el nmero de platos y de etapas es igual que
en el mtodo de Ponchon. Analticamente se puede determinar el nmero de platos
mnimos usando la siguiente ecuacin.
1
log ( 1 ) ( )
+1 =
Donde;
(1
= (1 )
(1 )
=
(1 )
8
0.98(1 0.9547)
= = 2.325
0.9547(1 0.98)
0.145(10.05)
= = 3.22
0.05(10.45)
0.9547 1 0.9547
log ( )( )
+1 = 1 0.9547 0.05
log(2.7376)
2.6025
+1 = = 5.9534
0.4371
1
=
1
9
Para MccCabey Thiele
10
9. Hallar el nmero de platos tericos de la torre.
# Etapas = 9.4
# Platos = 8.4 9
Ver grafica 3
# Etapas= 10
# Platos= 9
Ver grafica 4.
11
Ahora obtenemos el perfil de temperaturas.
160
140
120
100
T (C) 80
60
40
20
0
0 2 4 6 8 10 12
Nmero de plato.
12. Determine los requerimientos energticos; o sea, las cargas trmicas del
condensador y del rehervidor.
Los siguientes datos fueron obtenidos del diagrama de Ponchon & Savarit.
R 1.5Rm 1.3978
Q' H G1
R
H G1 H L 0
Q'53700
1.3978
53700 21400
12
KJ
Q' 98850
Kmol
QC QC
Q' H D 21400
D 30.0165
KJ
QC 1873029.6
h
= 58200 /
Ahora,
= = 24200
28,2854
= 2330716.96 /
14. Disee un plato para esta torre. Recuerde el problema propuesto # 6.8 del
Treybal. Tambin recuerde el ejemplo # 6.3 del mismo texto.
n RT (0,0199)(0,08205)(382,56)
Qgas = = = 0,627 3 /
P 1
Perforaciones
Ao area del orificio do 2
Aa
= area activa =0,907 ( p' ) (6.1)
13
Entonces,
p' =(2,5)(5)=12,5mm
Y el espesor de la hoja es un medio del dimetro del orificio, este espesor equivale a
2,25mm.
Ahora la ecuacin (16)
Ao 5 2
Aa
=0,907 (12,5
) =0,14512
Dimetro de la torre
La velocidad superficial del gas durante la inundacin est dado por la ecuacin 6,29
del Treybal
12
= (
) (6.2)
es la tensin superficial
La relacin de velocidades superficial es igual a:
0,5 0,5 0,0032 659 0,5
( ) ( ) = ( )( ) = ( )( ) = 0,073
0,627 3,20
Por ende los valores son =0,0489 y =0,0302; reemplazando estos valores en la
ecuacin (6.3)
La constante de inundacin CF = 0,0419 y ahora en la ecuacin (6.2) la velocidad
superficial del gas es VF =0,599 m/s
Valores pequeos de V se usan para el diseo real y para lquidos que no hacen
espuma es normal un valor del 80%- 85% de VF
As el valor de V =(0,8)( 0,599)=0,48 m/s
Para el dimetro o longitud del derramadero suele ser 60-80% del dimetro de la torre,
escogiendo tentativamente 70%, W=0,70T. De la tabla 6,1, se tiene que el rea del
plato utilizado por un vertedero cuando W=0,70T es 8,808%.
De modo que el rea total
14
1,31
= 10.08808 = 10,08808 = 1,43 2 (6.7)
(1,4)2
= = 1,4343 2
4
Derramaderos
Cresta del derramadero h1 y altura del derramadero hw
Weff se representa como una cuerda del crculo de dimetro T, una distancia h1, ms
alejada del centro que el derramadero real, que la geometra de la figura 6.16 del
Treybal sugiere a partir de la ecuacin 6.33, la ecuacin 6.34
0,5 2
2 2 2
(
) = () {[() 1] + 2 1 } (6.9)
Reemplazando los valores la relacin ( )=0,9988, se cumple la condicin Wef W,
por tanto h1=8,75x10-6
Fijamos hw, altura del derramadero hw=50mm=0,05m
La distancia h3, la diferencia en el nivel del liquido dentro e inmediatamente afuera del
vertedero, ser la suma de las prdidas de presin que resultan del flujo del lquido y
del gas en el plato superior (ecuacin 6,44):
3 = + 2 (6.10)
= + + (6.11)
15
hD: cada de presin en el plato seco (ecuacin 6.36, treybal)
2 4 2
= 2
[0,40 (1,25
) + (1
) ]
(6.12)
CO es un coeficiente de orificio que desprende del espesor del plato/ dimetro del
orificio (ecuacin 6.37)
0,25
= 1,09 ( ) (6.13)
()
De la tabla 6.2 se toma el valor de = 0,39, por ende el valor de
l=1,95mm=0,00195m.
Reemplazando en (6.13), = 1,370.
0,627
= = = 4,50 /
0,1391
El factor de friccin de Fanning, f, se toma de una tabla estndar, para esto debemos
calcular el Reynolds del gas y como consecuencia de esto la viscosidad del gas.
= 0,00029/
= 0,00049/
= 0,00039/
(0,005)(4,50)(3,20)
= = = 185048,535
0,00039
= 0,0178
4.428
=
= 2,4224 = 0,6539/ (6.17)
Sustituyendo en (6.15)
= 0,0318
16
6 6(0,04)
= = 659(0,005)(9,8) = 0,0074 (6.18)
0
es la menor de dos reas, la seccin transversal del vertedero o el rea libre entre
el vertedero y el plato.
El rea para el flujo del lquido debajo del faldn 0,025x0,9459=0,0236m2 y esta es la
menor rea
En (6.19)
3 0,0032 2
2 = ( ) = 0,0028m
2(9,8) 0,0236
Finalmente h3
Verificacin de inundacin
+ 1 + 3 < 2 (6.20)
0,50
0,05 + 0,0005 + 0,0525 = 0,1030 < = 0,25
2
Verificar lloriqueo
Si la velocidad del gas a travs de los orificios es muy pequea, el lquido gotear a
travs
de ellos y se perder el contacto sobre el plato para el lquido.
Para = 0,70; el derramadero se coloca en 0,3296 = 0,3296(1,35) = 0,445
Por lo tanto
= 2(0,445) = 0.8908
17
Segn este valor el plato no presenta lloriqueo excesivo hasta que la velocidad del gas
a travs de los orificios se reduzca cerca de este valor
Arrastre:
0,5
( ) = 0,07
= 0,8
De la grfica 6.17 (p. 196) se puede calcular el arrastre en platos perforados (E)
= 0,043
18. Si hace una destilacin flash para la condicin del alimento inicialmente
propuesta Qu esperara como resultado en las composiciones de cabeza y de
cola para una relacin de residuos a destilado de 50%?
Cuando se utiliza vapor vivo la relacin de reflujo y la composicin del destilado estn
fijos el nmero de platos aumenta, pero al eliminar el rehervidor se ahorra mas dinero
que el que se gasta al poner mas platos.
18
20. Analice su problema inicial si se decide a usar un condensador parcial
Cmo y en qu cambia?
Cuando se decide usar un condensador parcial se hace con frecuencia cuando resulta
muy elevada la presin que se requiere para la condensacin total del vapor, en el
diagrama de Ponchon y Savarit la condicin del destilado y el reflujo no es la misma,
entonces la lnea DL0 no es recta, entonces se debe construir una nueva lnea de
enriquecimiento.
21. Analice su problema inicial si adiciona otra alimentacin ms; es decir doble
alimentacin Cmo y en qu cambia?
Segn el mtodo de Ponchon & Savarit, la construccin sobre el diagrama Hxy para
las secciones de la columna arriba de F1 y debajo de F2 es igual que para la columna
de una sola alimentacin, con los puntos D y W ubicados de igual forma. En la
seccin media entre las alimentaciones, el punto de diferencia W puede localizarse
considerando los balances de materia y entalpia hacia la parte superior o hacia el
fondo; el resultado neto es el mismo.
Por ejemplo para una mezcla de doble alimentacin liquida donde M (corriente ficticia
igual al flujo neto ascendente y saliente) cae sobre la lnea DF1y sobre ala lnea WF2.
Una solucin que representa la mezcla de alimentacin compuesta debe caer sobre la
lnea DW, con:
F1 z F1 F2 z F2 F1 H F1 F2 H F2
Z F , pr , H F , pr
F1 F2 F1 F2
Considerando el libro Separation process principles de Seader J.D and Henley Ernest.
19
Para la tercera parte la lnea de operacin se emplea el mismo procedimiento utilizado
en la parte dos de la torre con lo que se tiene el balance:
20
23. Analice los cambios generados en la relacin mnima de reflujo cuando se
tiene una curva de equilibrio cncava en diferentes sectores de la misma.
Sabemos que la relacin de reflujo mnimo Rmin va a estar dada por la siguiente
relacin:
+ 1
En el intercepto con el eje Y dicha lnea me da una relacin que permite calcular Rmin.
Ahora analizamos la situacin cuando se presentan concavidades como las que se
pueden ver en la siguiente figura;
21
Para determinar el Rmin. La lnea d enriquecimiento es uno de los puntos determinantes
de dicha relacin y est se encuentra en una posicin tangente a la lnea de equilibrio
en la zona de enriquecimiento; como lo podemos observar en el punto P, si la lnea
pasara por el punto K, el reflujo que se obtendra sera muy pequeo, comparado con
el obtenido por otras rectas, concluimos a medida que el intercepto va subiendo, la
relacin de reflujo disminuye, hasta el punto en que la lnea tiende a una lnea
horizontal y por tanto el reflujo mnimo es igual a cero.
22