EJEMPLO DE CÁLCULOS DE PRÁCTICA DE RECTIFICACIÓN POR LOTES DE UNA

MEZCLA BINARIA DE ETANOL-AGUA

Versión corta
Secuencia de cálculos.
I.- Tabla de datos experimentales.
I.1.- Datos iniciales obtenidos de la muestra del hervidor antes de operar:

Hervidor °C ∆h (cm) Vol (lts.)  (g/cm3) % X peso

Alimentación Temperatura Lect. nivel De gráfica Densímetro De gráfica
Datos 25 27.5 38.5 0.972 15

I.2.- Datos de operación de Rectificación a Reflujo Total, en régimen permanente:

Lectura rotámetro a Reflujo Total:

L0 RT = VRT 0.12 l/min

Nota: este dato, nos sirve para establecer la relación de reflujo constante en la operación
de la extracción de los cortes de Destilado (Re = Lo/D).
I.3.- Datos de operación de Rectificación a Relación de Reflujo Constante durante la
extracción del corte 1 de 10 minutos:
Lecturas de la Relación de flujos constantes en los rotámetros del destilado y del
reflujo (Re = L0/D1 constante):

Destilado (D1) 0.04 (l/min)

Reflujo (L01) 0.08 (l/min)

0.08 l/min
ℜ= =2
0.04 l/min
Tamaño del Destilado del corte 1:

Volumen (litros) 0.4

Tiempo (minutos) 10
Gvc1 (litros/minuto) 0.04

Datos obtenidos del Destilado del corte 1 y de la muestra final del hervidor:

Corte 1 Temp. °C  (g/cm3) % Peso (gráfica)

Destilado (XD1 ) 25 0.814 90
Residuo (XW1) 25 0.9804 10

I.3.- Previamente a todos los cálculos, se deberá desarrollar las gráficas no ideales
de la mezcla Etanol-Agua a una presión absoluta de 585 mm de Hg. En este punto
deberá mostrar la tabla de cálculos y los gráficos correspondientes T vs X-Y y X vs
Y.
II.- Conversión de datos experimentales iniciales a fracciones, masas y flujos
molares.

II.1.- Cálculo de la fracción mol de la masa inicial del hervidor (X ¿ ¿WA 0)¿ y la masa
inicial en moles (W 0 ¿:

Ecuaciones:
Fracción mol inicial del metanol:
%W A 15
PM A 46
X WA 0= = =0. 06459
%W A %W B 15 85
+ +
PM A PM B 46 18
X WA 0=0.06459
Masa inicial molar de la mezcla etanol-agua en el hervidor:

V 0ρ ( 3 8.5 l )( 0.9 72 Kg/l )

W 0= 0
= =1.8892 Kmol
PM prom [ 46 ( 0. 06459 ) +18 ( 0. 9354 ) ] Kg/ Kmol
W 0 =1.8892 Kmol

II.1.1.- Cálculo de las fracciones mol del Destilado del corte 1 (X Dt1) y de la masa final
del hervidor (X ¿ ¿Wt 1)¿ :

Fracción mol del destilado, del corte 1:

%W A 90
PM A 46
X Dt 1= = =0. 77885
%W A %W B 90 10
+ +
PM A PM B 46 18
X Dt 1=0.77885

Fracción mol del residuo en el hervidor, después del corte 1:

%W A 10
PM A 46
X Wt 1= = =0 .04167
%W A %W B 10 90
+ +
PM A PM B 46 18
 X Wt 1=0. 04167

II.2.- Cálculo del flujo molar de la masa evaporada a Reflujo Total (V 1RT =
L0RT):
Nota: En los rotámetros, las lecturas que nos dan directamente, son como flujo de agua,
es decir, tenemos que corregir el flujo que en realidad es una mezcla de etanol-agua, para
ello, el mismo proveedor de los rotámetros nos proporciona la siguiente ecuación:

litros /min de agua=litros /min de mezcla 2

√ 7.02 ρ
ρf −ρ
Adonde:

f =Densidad del flotador = 8.02 g/cm3 para acero inoxidable 316.

 =Densidad del fluido en condiciones de operación, g/cm3.

7.02 = Factor de corrección constante.

Nota: Aquí se usa la densidad (0.8292 g/ml), de la mezcla del Destilado del
corte 1, porque es la misma masa que del reflujo.
Flujo como agua a flujo de mezcla metanol-agua.

L0 RT
litros
de mezcla=
litros
min de agua❑
=¿
0.12 ( minl ) =0 .1348 l/min ¿

√ √
min 2 7.02 ρ 2 7.02 ( 0.8 14 )
ρf −ρ ( 8.02−0.8 14 )
V1RT = L0RT = 0.1348 l/min
Flujo de mezcla a flujo molar:

GV ( minl )∗ρ ( kgl )

¿(
( l
)∗0.8 14 ( )
0. 1348
kg
L0 RT =
kg
) ¿=
min l

( 46 ( 0.77885 )+18 ( 0.22115 ) ) (

=
Kmol
〖 PM A X + PM ( 1− X ) ¿ kmol
A B A kg
kmol ) min

−3 Kmol
L0 t 1=2.7564∗10
min
II.3.- Cálculo de los Flujos molares de los rotámetros del destilado y reflujo
que permanecieron constantes durante la extracción del corte 1:
La relación de reflujo utilizada en la extracción del Corte 1, fue la siguiente:
 L0 0.08 lpm
Re = = =2
D 0.04 lpm
II.3.1.- Conversión del Destilado: Considerando la ecuación anterior podemos
hacer la conversión de flujo de agua a flujo de mezcla.

D
litros
de mezcla=
litros
min de agua❑
=¿
0. 04 ( minl ) =0. 0449l /min ¿

√ √
min 2 7.02 ρ 2 7.02 ( 0.8 14 )
ρf −ρ ( 8.02−0. 814 )
Por lo tanto, el flujo molar del destilado será:

GVD 1 ¿ ( minl )∗ρ ( kgl )

¿(
( l
) ∗0. 814 ( )
kg
0. 0449
D=
kg
) ¿=
min l

( 46 ( 0. 77885 ) +18 ( 0. 22115 ) ) (

=9.1813 ¿ 10
−4 Kmol

kmol )
〖 PM A X + PM ( 1−X ) ¿ kmol
A B A kg min

−4 Kmol
D=9.1813 ¿10
min

II.3.2.- Reflujo: Conversión del Reflujo a flujo de mezcla.

litros
de mezcla=
litros de agua
min
=¿
0.08 ( minl ) =0.08983 l/min ¿

√ √
min 2 7.02 ρ 2 7.02 ( 0.8 14 )
ρf −ρ ( 8.02−0.8 14 )

GV ( minl )∗ρ ( kgl )

¿(
( l
)∗0.8 14 ( )
kg
0. 08983
L0 t 1 =
kg
) ¿=
min

( 46 ( 0. 77885 ) +18 ( 0.22115 ) ) (

l
=1.8369 ¿ 10
−3 Kmo

kmol )
〖 PM A X + PM ( 1−X ) ¿ kmol
A B A kg min

−3 Kmol
L0 t 1=1.8369 ¿ 10
min

II.3.3.- Masa evaporada (Vt1):

−3 Kmol −4 Kmol Kmol
V t 1 =L0 t 1+ D t 1=1.8369 ¿ 10 +9.1813 ¿ 10 =
min min min
 −3 Kmol
V t 1 =2.75503 ¿10
min

II.3.4.-Cálculo de la relación de reflujo constante, durante la extracción del

corte 1, utilizando los flujos molares:
−3 Kmol
1.8369 ¿10
L min
ℜt 1= 0 t 1 = =2Por lo tanto, se comprueba que: ℜt 1=2
D0 t 1 −4 Kmol
9.1813 ¿10
min
Nota: El resultado es exactamente el mismo aplicado con los datos de los
rotámetros como flujo de agua.
III.- Aplicación del método gráfico de Mc Cabe – Thiele
III.1.- Cálculo del número de platos teóricos y eficiencia global de la columna,
durante la extracción del corte #1 con rectificación a una relación de Reflujo
constante, utilizando el método gráfico de Mc Cabe Thiele.
1.- Con los datos obtenidos en el corte 1 de la fracción mol del destilado XDt1 = 0.77885, la
fracción mol del residuo XWt1 = 0.04167 y la relación de reflujo constantes de Ret1 = 2,
podremos trazar la recta de operación del domo, calculando la intersección “b” en el eje
de las ordenadas. Sabemos que:

R x
y= x+ D =mX+b
La ec. del domo es: R+1 R+1
X D 0.77885
Por lo tanto: b= = =0.25962
ℜ+1 2+1
b=0.2 5962
Por lo que tendremos los resultados en la gráfica # 3:
1.1.- Se localiza XD, en la línea de 45° y se tira la recta de operación del domo hacia el eje
de las ordenadas (Y), en el valor de b = 0.25962
1.2.- Una vez tirada la recta de operación del domo, se forman escalones, iniciando con
una horizontal desde la línea de 45° en XD, hasta tocar la curva de equilibrio y de aquí, se
traza una vertical hasta tocar la línea de operación del domo y así sucesivamente hasta
rebasar XW. El último escalón será el que rebase a XW.
1.3.- En donde el número de escalones, representa el Número de Platos Teóricos (NPT).
Por lo que tendremos los resultados en la gráfica # 3:
2.- Cálculo de la eficiencia global de la columna, durante la extracción del corte #1 y
manteniendo una relación de reflujo constante. Los datos encontrados en la gráfica # 3,
fue de NPT= 10, y los platos reales que tiene la columna (NPR), son 16 platos, por lo que:
NPT −1 10−1
ηGC = ( 100 )= ( 100 )=56.25 % ηGC =56.25 %
NPR 16
Nota: Esta eficiencia, tiene dos interpretaciones teóricas, una es que la columna esta
sobrada de platos ya que con solo 9 platos daría la misma eficiencia, y la otra es que
todos los platos trabajaron en promedio global a un 56.25% de eficiencia.
IV.- Cálculos de la masa final (Wt1), después de la extracción del corte #1, cálculo de
la fracción mol promedio teórica del destilado (X Dt1) y cálculo del tiempo de
operación teórico (t1), utilizando los datos experimentales.
1.- Cálculo de la masa final, para esto se utilizará la siguiente ecuación:

W t 1=W A 0
( )
X Dt 1−X WA 0
X Dt 1−X Wt 1 (
=1. 8892 Kmol
0. 77885−0. 06459
0. 77885−0. 04167 )
=1. 8305 Kmol

W t 1=1. 8305 Kmol

2.- Cálculo de la fracción mol promedio teórica del Destilado (XDt1prom. teorica), ecuación:
 W A 0 X WA0 −W t 1 X Wt 1 ( 1.8892 Kmol ) ( 0. 06459 )−(1.8305 Kmol)(0.04167)
X D prom . teo. = =
W A 0−W t 1 (1 .8892 Kmol−1.8305 kmol)
X D prom . teo. =0. 77 885
3.- Cálculo del tiempo de operación teórico:
Ecuación 1:
W A 0 −W t 1 W A 0 −W t 1 1.8892 kmol−1.8305 kmol
t teo= = = =63.9343 min
L Dt 1 −4 kmol
V (1− ) 9.1813 ¿ 10
V min
t teo=63.9343 min
Ecuación 2:
ℜ+1 2+ 1
t teo=
Vt1
( W A 0−W t 1 )= −3
( 1. 8892 kmol−1. 8305 kmol ) =63.9343 min
2.75503 ¿ 10 kmol /min
t teo=63.9343 min

Nota: Como podemos observar con ambas ecuaciones nos da el mismo tiempo de
operación teórico de 63.9343 min, y este tiempo teórico está suponiendo que toda la
masa que se evaporó, incluyendo el líquido y vapor que se encuentra dentro de la
columna se evaporó a la misma velocidad en que se extrajo la masa del destilado del
corte #1, que fue de .4 l/10 min = 0.04 l/min = 2.75503 ¿ 10−3 kmol/min.

V.- Método para “Predecir gráfica y matemáticamente” una concentración promedio

de Destilado (XD prom) y el tiempo de operación teórico (t teo), utilizando los datos
experimentales de esta práctica, desarrollando la curva de operación de la
Destilación y un método de integración de la ecuación de Rayleigh.

Lo primero que tenemos que hacer es desarrollar los datos para trazar la curva de
operación del proceso de destilación a reflujo constante, mediante el siguiente
procedimiento, mostrado en la siguiente tabla:
 Paso 1 Paso 2 Paso 3 Paso 4
Suponer XD Calcular b Graficar Calcular f(x)

Dato XD

1 0.8 0.26667
2 0.77885 0.2596
3 0.7 0.2333
4 0.6 0.2
5 0.5 .16667
6 0.4 0.1333
7 0.3 0.1

Paso 1.- Suponer valores en una columna de X D, por ejemplo, para este caso, podemos
iniciar con XD=0.8, posteriormente con 0.77885 (que es dato de la práctica),0.7, 0.6, 0.5,
0.4, y 0.3.
Paso 2.- En otra columna, para cada valor supuesto de X D, obtener su valor de la
ordenada al origen “b”, con el valor obtenido del Reflujo constante del corte 1, a saber:
Para XD=0.8, tenemos que:
XD 0.8
b= = =0.2667
ℜ+1 2+ 1
Para XD=0.7:
XD 0.7
b= = =0.2333
ℜ+1 2+ 1
Para XD=0.77885:
X D 0.77885
b= = =0.2596
ℜ+1 2+1

Paso 3.- Una vez obtenidos los valores correspondientes de XD y b, trazar las rectas de
operación en la gráfica de equilibrio. Resultarán líneas paralelas, sí las traza en una sola
gráfica, ya que tienen la misma pendiente. Con el número de etapas teóricas de operación
del corte 1 igual a 10, trace los escalones en cada recta de operación y obtenga los
valores correspondientes de XW, utilizando los 10 escalones exactos. Ver gráficas
siguientes.
Nota: Para fines didácticos, se trazaron las líneas de operación, una en cada gráfica, ya
que en una sola gráfica no se visualizan bien los datos obtenidos.
4.- En una última columna, una vez obtenido el valor correspondiente de X W, para cada
valor de XD, se calcula la función de (X) a cada valor obtenido, que es igual a la inversa de
la diferencial de XD menos XW, es decir:
Ejemplo para XD= 0.8 y correspondiente valor gráfico de XW= 0.7, tendremos:
1 1
f ( x ) 1= = =10
X D− X w 0.8−0.7
para XD= 0.7 y correspondiente valor gráfico de XW= 0.18:
1 1
f ( x ) 2= = =1.92
X D− X w 0.7−0.18
Por lo que, deduciendo:
1 1
f ( x ) 3= = =1.7135
X D− X w 0.6736−0.09

1 1
f ( x ) 4= = =1.7699
X D− X w 0.6−0.035

1 1
f ( x ) 5= = =2.0833
X D− X w 0.5−0.02

1 1
f ( x ) 6= = =2.5974
X D −X w 0.4−0.015

1 1
f ( x )= = =3.4483
X D −X w 0.3−0.01
Finalmente, todo se resume en la siguiente tabla de resultados:

Paso 1 Paso 2 Paso 3 Paso 4

Dato XD

1 0.8 0.2992 0.700 10.0000

2 0.7 0.2618 0.180 1.9200
3 0.6736 0.2516 0.090 1.7135
4 0.6 0.2244 0.035 1.7699
5 0.5 0.1870 0.020 2.0833
6 0.4 0.1496 0.015 2.5974
7 0.3 0.1122 0.010 3.4483

Paso 5.- Ahora sí graficamos las dos últimas columnas de X W contra la f(x), tendremos el
siguiente diagrama, el cual nos facilitará el hacer una predicción, calculando el área bajo
la curva por medio de una integración gráfica utilizando el método del trapezoide para
calcular la fracción mol promedio (XDprom. teo) y el tiempo de operación teórico (t teo). Ver
siguiente gráfico:
5.- Con esta gráfica de operación podremos hacer una predicción, suponiendo que
realizaremos otra destilación con los mismos datos originales de la mezcla: X WA0=0.1753 y
WA0=1.52 Kmol y que la llevaremos a una concentración final supuesta de X WAt=0.05,
¿Cuál sería la masa final (Wt), la fracción mol promedio (XDprom teo) y el tiempo de
operación teórico (t teo)?
a).- Lo primero que hacemos es agrandar el área de integración para tener una mayor
exactitud en las lecturas gráficas.
b).- Trazar los límites del área bajo la curva, con los datos de X WA0=0.1753 y XWAt=0.05 y
obtener los datos de la f(x) correspondientes:

c).- Dividir el área bajo la curva en partes iguales nones, de acuerdo al método de
integración del trapezoide, en este caso la dividiremos en 3 partes iguales y tomaremos
sus datos de la función de X (f(x)) y de sus XW correspondientes:
Por lo que finalmente tendremos estos datos, para desarrollar la integración de la
ecuación de Rayleigh:

Área XW f (x)
P#1 0.1753 1.90
P#2 0.1330 1.82
P#3 0.0910 1.70
P#4 0.0500 1.75

Ecuación de Integración de Rayleigh, por el método de integración del trapezoide:

( ) [ ]
X Wt
dX Wt 1
∫ Y −X =ln W0
=∆ X
Y −X
prom .
XW 0

De donde:
∆ X=X Wt−X W 0=0.05−0.1753=−0.1253 y

[ 1
Y −X ]
prom .=
[
1.9+1.82+ 1.82+ 1.70+1.70+1.75
6
=1.78166
]
Por lo tanto:

ln
( )
Wt
W0
=(−0.1253 ) [ 1.78166 ] =−0.223242

Despejando a Wt, tendremos:

−0.223242
W t =e ( W 0) =0.79992 ( 1.52 Kmol )=1.21588 Kmol
Con estos datos, hacemos los cálculos de la posible concentración del Destilado y el
tiempo estimado de operación:
W A 0 X WA0 −W t 1 X Wt 1 ( 1.52 Kmol )( 0.1753 )−(1.21588 Kmol)(0.05)
X D prom . teo. = =
W A 0−W t 1 (1.52 Kmol−1.21588 kmol)
X D prom . teo. =0.676253

W A 0 −W t 1 W A 0 −W t 1 1.52 kmol−1.21588 kmol

t teo= = = =56.528 min
L Dt 1 kmol
V (1− ) 0.00538
V min
t teo=56.528 min
Conclusión: Lo que podemos ver es que el promedio teórico de la fracción mol del
destilado (XD prom. teo.=0.676253) prácticamente permanece constante desde una
concentración del residuo inicial de XW0=0.1753 hasta una concentración residual final de
XWt=0.05, esto es siempre y cuando se mantenga la relación de reflujo constante de
Re=1.67435, durante la extracción del Destilado. También vemos que el tiempo de
operación teórica aumenta casi al doble (56.528 min).

NOTA: FALTA TABLA DE RESULTADOS