CÁLCULOS

Tabla de datos.
2
Presión de Trabajo: 0.25 Kg/cm Presión de Vapor: 1.0 Kg/cm2
Temperatura de alimentación TF = 390 C P atm = 585 mm de Hg
 TABLA DE DATOS EXPERIMENTALES (METANOL-AGUA)
Flujo
Temperatura Densidad
%peso %mol %Rotámetro 100%
(°C) (g/cm3)
(l/min)

Alimentación 25 0.975 0.14 0.084 29 4.5

Destilado 25 0.818 0.89 0.82 10 1.01
Plato n-1 (4) 0.834 0.84 0.75

Plato n(5) 0.846 0.79 0.68

Plato n+1 (6) 0.855 0.75 0.63

Residuo 25 0.998 0.02 0.0113 72 1.64

Reflujo 25 0.818 0.89 0.82 29.5 1.64

Temperaturas de los platos

TPLATO1 73 68
TPLATO2 76 68
TPLATO3 78 68
TPLATO4 82 70
TPLATO5 91 72
TPLATO6 100 73
TPLATO7 102 78
TPLATO8 104 84
TPLATO9 104 87
TPLATO10 103 88
TPLATO11 104 86
TPLATO12 104 90
TPLATO13 103 92
TPLATO14 104 95
TPLATO15 102 95
TPLATO16 - -
TPLATO17 93 86
THERVIDOR 104 99
 1. Cálculo de la presión absoluta para desarrollar las gráficas de equilibrio NO
IDEALES:
Ecuación: Pabs= Pman + Patm
760 mmHg
Pabs=0.25 kg/cm2 +585 mmHg
1.033 kg /cm2
Pabs=768.9 mmHg=769 mmHg
El desarrollo de las gráficas se aprendió en la sesión de Introducción por lo
cual se da por hecho que ya el alumno, sabe desarrollar cualquier gráfica a
determinada presión absoluta de operación.
2. Cálculo de los flujos molares
a) Cálculo de fracciones mol.

% peso A / PM A
x=
Ecuación: % peso A % peso B
+
PM A PM B

14
32
xf = =0.084
14 86
+
32 18

89
32
xD= =0.8198
89 11
+
32 18

b) Cálculo de los flujos como agua

l
FG f H O =4.5 ∗0.29=1.305l /min
2
min
l
DG D H O =1.01 ∗0.10=0.101 l/min
2
min
l
WGW H O=1.64 ∗0.72=1.1808l /min
2
min
l
LG W H O =1.64 ∗0.295=0.4838l/min
2
min

c) Cálculo de los flujos como mezcla.

GV H O
GVmezcla = 2

Ecuación: 7.02 ρx
√ 8.02− ρx
 Considerando que 1g/cm3 es igual a 1 kg/l se tiene que:
1.305 l/min
GVf = =1.3239l /min
7.02(0.975)
√ 8.02−(0.975)
0.102 l/min
GVD= =0.1142 l/min
7.02( 0.818)
√ 8.02−( 0.818)
1.1808 l/min
GVW = =1.1821l/min
7.02( 0.998)
√ 8.02−( 0.998)
0.4838 l/min
GVLD = =0.5418 l/min
7.02(0.818)
√ 8.02−(0.818)

3. Cálculo de los flujos molares

ECUACIÓN : F=
G Vf ( minl )∗ρ( kgl )
kg
( PM ¿ ¿ A X A + PM B ( 1−X A ) )( )¿
kmol

F=
1.3239 ( minl )∗0.975( kgl ) =0.0673 kmol/min
kg
32 ( 0.084 )+ 18(0.916)( )
kmol

kmol
∗60 min
min kmol
0.0673 =4.0380
1h h

D=
0.1142 ( minl )∗0.818( kgl ) =0.00317 kmol
min
∗60 min
=0.1902
kmol
kg 1h h
32 ( 0.82 )+ 18(0.18)( )
kmol
 W=
0.1821 ( minl )∗0.998( kgl ) =0.0649
kmol
min
∗60 min
=3.8940
kmol
kg 1h h
32 ( 0.012 ) +18(0.998)( )
kmol
l kg kmol
L D=
0.5418 ( )min
∗0.818( )
l
=0.9024
min
∗60 min
=0. 9024
kmol
kg 1h h
32 ( 0.82 ) +18 (0.18)( )
kmol

Balance de materia en la columna (Entradas = Salidas)

F= D + W
Sustituyendo datos, tenemos:

4.0380 kmol/h= 0.1902 kmol/h + 3.8940 kmol/h

4.0380 kmol/h ≠ 4.0842 kmol/h

Conclusión: Como la salida (D+W = 4.0842), es un poco mayor que la entrada

(F=4.0380), se dice que se estaba operando con tendencia a secar la columna,
si los valores hubieran estado invertidos, es decir, que la entrada fuera mayor
que la salida se dice que la columna se estaba operando con tendencia a
inundar la columna. Una vez conocido el balance se hace el ajuste, ya sea
aumentando F o disminuyendo W y/o D. Según sea lo más conveniente para la
operación, hay que hacer un análisis previo.

4. Cálculo de reflujos externo e interno.

Cálculo del Reflujo externo:

L D 0.9024 kmol /h
ℜ= = =4.7444
D 0.1902 kmol /h
Calculo de Reflujo interno:
Ri=m( pendiente de larecta del domo )

ℜ 4.7444
Ri= = =0.8259−−−−− pendiente
ℜ+1 4.7444 +1

Ecuación de la recta del domo:

ℜ XD
Y n= ∗X n−1 +
ℜ+1 ℜ+1
0.82
Y n=0.8259∗X n−1+
4.7444 +1
Y n=0.8259∗X n−1+ 0.14274−−−−Ecuación del domo

Nota: Con esta ecuación se traza la recta de operación en la gráfica de

equilibrio X vs Y, usando XD=0.82 y b=0.14274 (intersección al eje de las
ordenadas (Y)), ver gráfica 2.

Confirmación del valor de Re con datos de cualquier plato.

Para el caso del plato #5: Sí X n−1 = 0.68 y sustituyendo en la ecuación del
domo, tenemos:
Y n=0.8259∗(0.68)+ 0.14274
Y n=0.704352

X D −Y n 0.82−0.704352
De la ecuación: ℜ= = =4.749
Y n−X n −1 0.704352−0.68

Conclusión: El resultado es casi exacto, sí cumple contra:

Re = Lo/D = 4.744 (Hay una pequeña diferencia por los decimales).

Eficiencia de la columna.

Para calcular la eficiencia de la columna, es necesario obtener la ecuación

de la recta de operación del área de alimentación (F) y para ello se requiere
calcular (q), que es la relación del calor necesario para evaporar una mol de
la alimentación contra el calor necesario para evaporar una mol del líquido
en el plato de alimentación.
Ecuaciones:
H −h
q= v f
H v −h L
Sí consideramos a Treferencia.= TF, las ecuaciones de las entalpías quedan
como sigue:

H v =(C pA X A +C pB X B )( T R−T F ) + λ A X A + λ B X B
h L =(C pA X A +C pB X B ) ( T B−T F )
h F =( C pA X A +C pB X B ) ( T F −T F ) =0

Los valores de Cp´s y calores latentes, se obtienen de tablas, y T R y TB de

la gráfica T vs XY (ver gráfica 1). Donde:

CpA = 19.38 Kcal/Kmol 0C CpB = 18 Kcal/Kmol 0C

 λ A = 8415 Kcal/Kmol λ B = 9718 Kcal/Kmol.

De gráfica 1, tenemos: TR = 98 0C y TB = 88.6 0C; sabemos que: TF = 39 0C

Sustituyendo datos, tenemos:

H v = (19.38∗0.084+ 18∗0.916 )( 98−39 ) + [ 8415(0.084)+9718(0.916) ]

H v =10677.39 Kcal/Kmol

h L =(19.38∗0.084+18∗0.916) ( 88.6−39 )=923.91 Kcal/Kmol

10677.39−0.0
q=
10677.39−923.91

q=1.095

Conclusión: Tenemos un valor de q>1.0, lo cual indica que se alimentó un

líquido subenfriado y por lo tanto su pendiente es positiva, es decir, tendencia
hacia la derecha.

Conociendo el valor de q, se procede a calcular la ecuación de la zona de

alimentación.

q XF
Ecuación: Y q= X q−
q−1 q−1
Sustituyendo datos tenemos:

1.095 0.084
Y q= X −
0.095 q 0.095

Y q =11.5263 X q −0.8842

Esta ecuación se utiliza para trazar la línea de operación en la gráfica y se

requiere tener dos puntos de referencia y uno es X F localizado en la línea
de 450, y el otro punto se obtiene con la ecuación calculada anteriormente,
como sigue:

Ejemplo: Sí consideramos que Y q =1 y despejando a Xq, tendremos:

1+ 0.8842
X q= =0.1634
11.5263
 Con esto obtenemos el otro punto de referencia, siendo P (0.1634, 1.0),
para trazar la línea de operación de la zona de alimentación (ver gráfica 2).

Para trazar la línea de operación del fondo, no se requiere hacer cálculos,

únicamente se localiza XW en la línea de 450 y se traza la recta al punto de
intersección de las líneas de operación del domo y alimentación, (Ver
gráfica 2).

Finalmente se trazan los escalones con rectas paralelas a los ejes de la

gráfica, iniciando desde el punto XD localizado en la línea de 450 y tomando
como referencia la curva de equilibrio. Tomando en cuenta que durante el
trazado de los escalones se hace antes del punto de intersección de las
tres rectas de operación, sobre la línea del domo y posteriormente a la
intersección sobre la línea del fondo y hasta rebasar X w, se detiene el
trazado de escalones, (Ver gráfica 2). De la gráfica 2 obtenemos los
siguientes datos:

Numero de platos teóricos = 5 (hasta rebasar X W)

Plato de alimentación teórico = # 3.

De la siguiente ecuación de la eficiencia tenemos:

NET −1
η= ∗100
NER
Donde:
NET= Número de platos teóricos = 5
NER= Número de platos reales (En este caso nuestra columna de
destilación tiene 17 platos).
El restarle 1 al NET, significa que este procedimiento considera que el
hervidor funciona como un plato.

5−1
ηGc = x 100=23.53%
17

Conclusión:
5. Cálculo Lade
de la masa teoría
vapordice que con 4 platos es suficiente, por lo tanto,
utilizada
solo se está usando el 23.53% de la capacidad de la columna, o que
en promedio todos los platos trabajaron al 23.53%.

Cálculo de la Masa de vapor utilizada en la experimentación.

 Diagrama general de una columna de destilación para hacer el balance de
calor:
V1
QC = D(Re+1)

LD
D = DhD

QF F
FhF

= Mv

L QH
W Qw= W
Ecuación del balance general de calor de la columna:
Entradas = Salidas
Q F +Q H =Q c + Q D +Q w
Despejando:
Q H =QC +QD + QW −Q F
Sabemos que:
Q H =M V λV ; QC = D ( ℜ+1 ) λm; Q D=Dh D ; Q W =WhW ; Q F=Fh F , , y
λ m=λ A X A + λ B X B
Sustituyendo tenemos:
M v λv =D ( ℜ+1 ) λ m+ D hD +W hW −F hF

Despejando de la ecuaci ó n anterior la masade vapor

D ( ℜ+1 ) λm + DhD +Wh W −Fh F

M v=
λv

El calor latente del vapor, se obtiene con su presión absoluta.

λ v con Pabs vap=Pman+ Patm
kg mmHg∗1.033 kg/cm 2
Pabs vap=1.0 2 +585
cm 760 mmHg
Pabs vap=1.795138=1.8 kg/cm2

Para obtener el λ v se usan tablas de propiedades del agua líquido vapor,

por lo tanto;
1.8 kg /cm2 corresponde a λv =528.5 kcal /kg y una TW = 116.5 0C
Por otro lado, con la gráfica #1, y con XD = 0.82, tenemos: TD = 67.5°C
 Nuevamente sí consideramos a T referencia.= TF, las ecuaciones de las
entalpías quedan como sigue:

h D =(C pA X A +C pB X B ) ( T D−T F )
hW =(C pA X A + C pB X B) ( T W −T F )
h F =( C pA X A +C pB X B ) ( T F −T F ) =0

Los valores de Cp´s y calores latentes, ya los tenemos y son:

CpA = 19.38 Kcal/Kmol 0C CpB = 18 Kcal/Kmol 0C

λ A = 8415 Kcal/Kmol λ B = 9718 Kcal/Kmol.

Por lo que haciendo las sustituciones respectivas, procedemos a hacer los

cálculos correspondientes:
λ m=8415 ( 0.82 ) +9718 ( 0.18 ) =8649.54 Kcal /kmol

h D = [19.38 (0.82) + 18 (0.18)] (67.5 – 39) = 545.25 Kcal/Kmol

hW = [19.38 (0.0113) + 18 (0.9887)] (116.5 – 39) = 1396.21 Kcal/Kmol

λ m=8415 ( 0.82 ) +9718 ( 0.18 ) = 8649.54 Kcal/Kmol

Por lo que sustituyendo los resultados anteriores en la ecuación del

balance general, tendremos:

( 0.1902 ) ( 4.744+1 ) ( 8649.54 )+ ( 0.1902 )( 545.25 ) +3.894 ( 1396.21 )−0

M v=
528.5

M v =28.364 Kg /h

Conclusión: El consumo de vapor en la experimentación

efectuada en esta columna de Destilación fue de 28.364 Kg/h.

6. Eficiencia en platos de la columna de Rectificación.

A cada plato se le puede calcular tanto la eficiencia líquida, como la

eficiencia evaporativa.

Ecuación general.
 datos experimentales reales
ηp=
datosequilibrio ideal o teórico

Eficiencia líquida de un plato:

X n−1−X n
η plato = x 100
X n−1−X ¿n
nlíquido

Nota: Xn* en el equilibrio es obtenida en la gráfica de equilibrio.

Cálculo de la eficiencia del Plato 5,

P 4−P5 0.75−0.68
η plato 5 = = x 100=21.875 %
P 4−P 5¿ 0.75−0.43
líquido

Con P4 se obtiene P5*. Ver resultados en la gráfica #3.

Para la eficiencia de vapor, Murphree establece la siguiente regla:

“La concentración del vapor que emite un plato es igual a la concentración
del líquido del plato superior”, es decir:

Y n= X n−1; Y n +1=X n ; etc.

Por lo tanto, sí la eficiencia evaporativa de los platos, se puede poner en

función de las concentraciones líquidas de los platos, como sigue:

Y n−Y n+1 X n−1−X n

η plato 5 = ¿ = ¿ x 100
vapor
Y n−Y n+1 Y n−X n

Cálculo de la eficiencia evaporativa del plato 5:

P 4−P 5 0.75−0.68
η plato 5 = = x 100=36.84 %
vapor
P 4 ¿−P 5 0.87−0.68

Ahora con P5, se obtiene P4* en la gráfica de equilibrio. Ver gráfica #3.
Una vez comprendido lo anterior, se calculan rápidamente las eficiencias del
PLATO #6:
0.68−0.63
η plato 6 = x 100=13.51%
líquido
0.68−0.31
0.68−0.63
η plato 6 = x 100=23.26 %
vapor
0.845−0.63
 Conclusión: Los valores de los porcentajes en la eficiencia de los platos 5
(21.875%) y 6 (13.51%), nos indica que al ir bajando el líquido en la columna
va perdiendo concentración del más volátil, se va agotando, y para el vapor
que va subiendo en la columna la eficiencia del plato 5 (36.84%) y del plato 6
(23.26%).nos indican que el vapor se va enriqueciendo en los más volátiles.
La eficiencia de los platos es directamente proporcional a la concentración de
los volátiles.

NOTA: FALTA TABLA DE RESULTADOS

6. GRÁFICAS

TR=98

TB=88.6

XF =0.084
 Gráfica # 1: Diagrama T vs X-Y (mezcla metanol-agua) @ 769 mm de Hg
NO IDEAL

105

100

liquido Vapor
95

90
Temperatura, T (°C)

85

80

75

70

65

60
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Fracción molar x, y
 P (0.16, 1.0)
𝝶P5VAP= (0.68, 0.87)

𝝶P5LIQ= (0.43, 0.75) XD=0.82

1 = (0.63, 0.845)
𝝶P6VAP

𝝶P6LIQ= (0.31, 0.68)

b=0.14

4
Xf =0.084

XW=0.011
 Gráfica # 2: Diagrama X vs Y (mezcla Metanol-agua) @ 769 mm de Hg
NO IDEAL
1.0

0.9

0.8

0.7
Fracción molar del vapor, y

0.6

0.5

0.4

0.3

0.2

0.1

0.0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Fracción molar del líquido, x

 Gráfica # 3: Diagrama X vs Y (mezcla Metanol-agua) @ 769 mm de Hg
NO IDEAL
1.0

0.9

0.8

0.7
Fracción molar del vapor, y

0.6

0.5

0.4

0.3

0.2

0.1

0.0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

Fracción molar del líquido, x