CÓDIGO DE LA PRÁCTICA TM 01 y TM 02
PRÁCTICA
TM 01 – DETERMINACIÓN EXPERIMENTAL DE DATOS DE
EQUILIBRIO
TM 02 - RECTIFICACIÓN AEPT DE RELLENOS Y ET DE
PISOS
NOMBRE DEL PROFESOR AYUDANTE:
Marcos Larriba
FECHA REALIZACIÓN FECHA DE ENTREGA
08/02/2021 22/02/2021
NOMBRE DEL ALUMNO DNI
Manuel Liébana Díaz 02749119K
Beatriz López Cabezas 51555051E
Pablo López Piñeiro 02303988D
CÓDIGO DE LA PRÁCTICA
PRÁCTICA
CALIFICACIÓN
PRESENTACIÓN
CONCEPTO MEMORIA
EXPTAL. INTRODUCCIÓN
Presentación 10%
Introducción 10% PROCEDIMIENTO
Instalación
10% RESULTADOS
Procedimiento
Resultados
50% DISCUSIÓN
Discusión
Conclusiones 20%
CONCLUSIONES
NOMBRE DEL PROFESOR AYUDANTE:
FECHA RECOGIDA CONFORME
NOMBRE DEL ALUMNO DNI
2
1. INTRODUCCIÓN Y OBJETIVOS:
La destilación es un proceso que consiste en separar los componentes o sustancias de una mezcla
líquida mediante el uso de la ebullición selectiva y la condensación. La destilación puede resultar
en una separación esencialmente completa (componentes casi puros), o puede ser una separación
parcial que aumenta la concentración de los componentes seleccionados en la mezcla.
Los diagramas de equilibrio de mezclas binarias entre dos fases suelen ser isóbaros o isotermos y
en menor medida diagramas a composición constante de un componente en una de las fases. Otra
forma de representar las composiciones de las dos fases en equilibrio es en gráficos en los que se
muestra la composición de un componente de una de las fases frente a la composición del mismo
componente en la otra fase en equilibrio. Estos gráficos se denominan diagramas de distribución,
y pueden ser isóbaros o isotermos.
Figura 1. Diagrama isóbaro Figura 2. Diagrama de distribución
El objetivo de esta práctica es determinar experimentalmente datos de equilibrio 1íquido vapor de
una mezcla binaria real y utilizar estos datos experimentales para determinar el número mínimo
de pisos de tres columnas de rectificación, mediante el método gráfico de McCabe-Thiele y el
método analítico de Fenske, así como la altura equivalente de piso teórico (AETP) de los rellenos
y la eficacia global de los pisos.
2. INSTALACIÓN Y PROCEDIMIENTO EXPERIMENTAL:
La práctica se realiza en tres instalaciones diferentes: una columna de rectificación, un ebullómetro
y un cromatógrafo de gases.
3
En primer lugar, se introduce la mezcla tolueno-ciclohexano con la ayuda de un embudo en el
ebullómetro (figura 1). La mezcla empieza a hervir por la acción de un calefactor eléctrico, el
vapor se dirige hacia la parte superior del ebullómetro y con la ayuda de un refrigerante condensa
mientras que el líquido se recoge por la parte inferior. El equilibrio se alcanza cuando la
temperatura del vapor y la del líquido son constantes e iguales. Una vez se alcanza el equilibrio,
se recoge una muestra del vapor y del líquido para analizarla posteriormente en el cromatógrafo
de gases. Este procedimiento se repite nueve veces y cada grupo debía estudiar tres mezclas de
composición variable.
Figura 1. Ebullómetro
Rectificación en torres de relleno
En esta parte de la práctica se va a llevar a cabo la determinación de las composiciones del piso 1
y destilado de la columna de relleno (figura 2) en cada uno de los experimentos para
posteriormente calcular el número de pisos.
4
Figura 2. Columna de rectificación
Primero encendemos el calentador situado debajo de la alimentación a la columna que hará ebullir
al compuesto más volátil (ciclohexano) y esperamos hasta que el termómetro de las corrientes de
destilado y de piso 1 marquen una temperatura constante que nos refleja que se ha alcanzado el
equilibrio líquido-vapor.
Tomamos el valor de ambas temperaturas y recogemos las muestras de destilado y de piso 1 que
se llevarán finalmente al cromatógrafo de gases.
Repetimos este procedimiento en los tres experimentos.
3. RESULTADOS Y DISCUSIÓN
En primer lugar, se usaron los datos bibliográficos de Myers para calcular los parámetros de
iteración binaria de la ecuación de Van Laar.
Tabla 1. Datos bibliográficos de temperatura y composición
T/K x1 y1 x2 y2
383,80 0,000 0,000 1,000 1,000
381,40 0,041 0,102 0,959 0,898
378,60 0,091 0,212 0,909 0,788
5
377,00 0,118 0,264 0,882 0,736
376,00 0,143 0,308 0,857 0,692
374,90 0,164 0,348 0,836 0,652
373,70 0,192 0,386 0,808 0,614
372,65 0,217 0,422 0,783 0,578
371,50 0,245 0,457 0,755 0,543
371,50 0,243 0,450 0,757 0,550
370,55 0,273 0,492 0,727 0,508
370,10 0,283 0,504 0,717 0,496
369,50 0,304 0,532 0,696 0,468
368,65 0,323 0,547 0,677 0,453
368,40 0,335 0,560 0,665 0,440
367,35 0,368 0,596 0,632 0,404
366,95 0,379 0,599 0,621 0,401
365,90 0,416 0,633 0,584 0,367
365,00 0,452 0,662 0,548 0,338
363,70 0,504 0,702 0,496 0,298
362,90 0,533 0,724 0,467 0,276
362,00 0,559 0,749 0,441 0,251
361,10 0,599 0,774 0,401 0,226
361,15 0,602 0,777 0,398 0,223
360,50 0,634 0,794 0,366 0,206
359,70 0,672 0,811 0,328 0,189
358,40 0,727 0,852 0,273 0,148
357,95 0,763 0,864 0,237 0,136
357,60 0,780 0,877 0,220 0,123
356,90 0,814 0,895 0,186 0,105
355,85 0,874 0,926 0,126 0,074
354,25 0,964 0,973 0,036 0,027
353,85 1,000 1,000 0,000 0,000
Con estos datos, se representan los valores de y1 frente a x1 para obtener el diagrama de
distribución.
6
Diagrama de distribución
1,0
0,8
0,6
y1
0,4
0,2
0,0
0,0 0,2 0,4 0,6 0,8 1,0
x1
Figura 3. Diagrama de distribución
Se representó a su vez el diagrama isóbaro del equilibrio líquido – vapor de la mezcla. Los puntos
experimentales se ajustaron a una ecuación de segundo grado para obtener una ecuación teórica
que nos permita relacionar la temperatura con la composición.
Diagrama T vs x1,y1
115,00
110,00
y = -6,3589x2 - 24,445x + 110,84
105,00 R² = 0,9995
100,00
T ºC
95,00
90,00
85,00 y = 20,637x2 - 48,707x + 109,5
R² = 0,9972
80,00
75,00
0,0 0,2 0,4 0,6 0,8 1,0
x1,y1
Figura 4. Diagrama isóbaro del equilibrio líquido – vapor
A continuación, se calcularon los coeficientes de actividad de las mezclas estudiadas. Para ello,
con las temperaturas experimentales obtenidas en el laboratorio se calcularon las presiones de
saturación en el equilibrio con la ecuación de Antoine (ecuación 1). Una vez obtenidas las
7
presiones de saturación, se obtuvo el coeficiente de actividad despejando la ecuación 2. Hay que
tener en cuenta que los datos bibliográficos fueron tomados a una presión de 760 mmHg.
(1)
(2)
Tabla 2. Presiones de saturación y coeficientes de actividad
T(K) Psat1(mmHg) Psat2(mmHg) ϒ1 ϒ2
383,80 1723,27 760,52 - 1,00
381,40 1622,58 710,03 1,17 1,00
378,60 1510,84 654,51 1,17 1,01
377,00 1449,68 624,37 1,17 1,02
376,00 1412,43 606,09 1,16 1,01
374,90 1372,31 586,47 1,18 1,01
373,70 1329,53 565,65 1,15 1,02
372,65 1292,95 547,92 1,14 1,02
371,50 1253,78 529,01 1,13 1,03
371,50 1253,78 529,01 1,12 1,04
370,55 1222,11 513,78 1,12 1,03
370,10 1207,32 506,69 1,12 1,04
369,50 1187,82 497,36 1,12 1,03
368,65 1160,62 484,37 1,11 1,05
368,40 1152,71 480,60 1,10 1,05
367,35 1119,94 465,03 1,10 1,04
366,95 1107,65 459,21 1,08 1,07
365,90 1075,87 444,20 1,07 1,08
365,00 1049,20 431,64 1,06 1,09
363,70 1011,58 414,01 1,05 1,10
362,90 988,96 403,45 1,04 1,11
362,00 963,97 391,83 1,06 1,10
361,10 939,48 380,48 1,05 1,13
361,15 940,83 381,10 1,04 1,12
360,50 923,42 373,06 1,03 1,15
359,70 902,34 363,34 1,02 1,21
8
358,40 868,89 347,99 1,03 1,18
357,95 857,54 342,80 1,00 1,27
357,60 848,79 338,81 1,01 1,25
356,90 831,50 330,93 1,00 1,30
355,85 806,08 319,40 1,00 1,40
354,25 768,51 302,44 1,00 1,88
353,85 759,33 298,32 1,00 -
Finalmente, se determinaron los coeficientes de actividad a dilución infinita mediante
extrapolación de las siguientes curvas:
ϒi vs xi (bibliográfica)
2,00
1,80
1,60
1,40
1,20
1,00
ϒi
gamma1
0,80
0,60 gamma2
0,40
0,20
0,00
0,0 0,2 0,4 0,6 0,8 1,0
xi
Figura 5. Coeficientes de actividad frente a composiciones
Extrapolando de la figura, podemos estimar el valor de γ1ꚙ , γ2ꚙ . Con estos valores de coeficiente
de actividad a dilución infinita se pueden estimar A12 y A21, que son los parámetros de iteración
binaria de la ecuación de Van Laar mediante la ecuación (3).
(3)
Tabla 3. Coeficientes de actividad a dilución infinita y parámetros de Van Laar
ϒ1ꚙ 1,167
ϒ2 ꚙ
2
A12 0,1544
A21 0,6931
9
A continuación, se tomaron datos experimentales de temperatura y composición y se realizaron
los mismos cálculos para determinar los parámetros de iteración binaria de la ecuación de Van
Laar. Para determinar las composiciones de las distintas mezclas, se analizaron las muestras en un
cromatógrafo de gases.
Los resultados de las áreas obtenidas en el cromatógrafo de gases se recogen en la siguiente tabla:
Tabla 4. Áreas obtenidas para las mezclas ciclohexano-tolueno
TM01 Fase Área
Mezcla 1 Vapor 5792
Mezcla 2 Vapor 14312
Mezcla 3 Vapor 15124
Mezcla 4 Vapor 19373
Mezcla 5 Vapor 24222
Mezcla 6 Vapor 23380
Mezcla 7 Vapor 26664
Mezcla 8 Vapor 27600
Mezcla 9 Vapor 28292
Mezcla 1 Líquida 3251
Mezcla 2 Líquida 9787
Mezcla 3 Líquida 8888
Mezcla 4 Líquida 14517
Mezcla 5 Líquida 21049
Mezcla 6 Líquida 16960
Mezcla 7 Líquida 24629
Mezcla 8 Líquida 23609
Mezcla 9 Líquida 27452
Con estos valores de área, se obtuvieron las composiciones de vapor y líquido para cada mezcla.
Para ello, se sustituyeron las áreas en la siguiente ecuación de la recta de calibrado:
Á𝑟𝑒𝑎 = 31135 × 𝐶𝑜𝑚𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛(𝑥𝑖 𝑜 𝑦𝑖)
Tabla 5. Composiciones de líquido y vapor para las mezclas ciclohexano-tolueno
x1 y1 x2 y2
Mezcla 1 0,104 0,186 0,896 0,814
Mezcla 2 0,314 0,460 0,686 0,540
Mezcla 3 0,285 0,486 0,715 0,514
Mezcla 4 0,466 0,622 0,534 0,378
10
Mezcla 5 0,676 0,778 0,324 0,222
Mezcla 6 0,545 0,751 0,455 0,249
Mezcla 7 0,791 0,856 0,209 0,144
Mezcla 8 0,758 0,886 0,242 0,114
Mezcla 9 0,882 0,909 0,118 0,091
Se representan los valores de y1 frente a x1 para obtener el diagrama de distribución.
Diagrama de distribución experimental
1,200
y = -0,8158x2 + 1,7186x + 0,0231
1,000 R² = 0,987
0,800
0,600
y1
0,400
0,200
0,000
0,000 0,200 0,400 0,600 0,800 1,000
x1
Figura 6. Diagrama de distribución
Como se puede observar en la figura 6, se han ajustado los datos bibliográficos a un polinomio de
grado 2. Sustituyendo en el polinomio obtenemos los ybibl y con la siguiente ecuación, calculamos
el error relativo:
|𝑦𝑒𝑥𝑝 − 𝑦𝑏𝑖𝑏𝑙 | × 100
𝐸𝑅 (%) =
𝑦𝑏𝑖𝑏𝑙
Tabla 6. Errores experimentales de las composiciones de líquido y vapor
y1 ybib Error y
0,18 0,23 21,86
0,45 0,52 11,73
0,48 0,48 0,073
0,62 0,68 8,65
0,77 0,84 7,56
0,75 0,75 0,19
0,85 0,89 4,80
0,88 0,88 0,13
11
0,90 0,93 2,32
A continuación, se calcularon los coeficientes de actividad de las mezclas estudiadas. Para ello,
con las temperaturas experimentales obtenidas en el laboratorio se calcularon las presiones de
saturación en el equilibrio con la ecuación de Antoine (ecuación 1). Una vez obtenidas las
presiones de saturación, se obtuvo el coeficiente de actividad despejando la ecuación 2. Hay que
tener en cuenta que los datos experimentales fueron tomados a una presión de 696,9 mmHg.
Tabla 7. Presiones de saturación y coeficientes de actividad
Temp (ºC) Psat1(mmHg) Psat2(mmHg) ϒ1 ϒ2
104 1455,33 627,14 0,85 1,00
99,7 1299,85 551,26 0,78 0,99
95 1144,83 476,86 1,03 1,05
92 1053,61 433,71 0,88 1,13
88,4 951,66 386,11 0,84 1,23
88,9 965,34 392,46 0,99 0,97
83,5 825,38 328,15 0,91 1,45
82,3 796,55 315,08 1,02 1,03
78,8 716,91 279,34 1,00 1,92
Se representó el diagrama isóbaro del equilibrio líquido – vapor de la mezcla. Para comparar
gráficamente los datos experimentales y los bibliográficos, se representaron los puntos
experimentales y el ajuste a una ecuación de segundo grado obtenido anteriormente.
Diagrama T vs x1,y1
115,00
110,00
105,00
100,00
T ºC
95,00
90,00
85,00
80,00
75,00
0,0 0,2 0,4 0,6 0,8 1,0
x1,y1
12
Figura 7. Diagrama isóbaro del equilibrio líquido – vapor
Se calculó a su vez el error relativo cometido en el experimento mediante la siguiente ecuación.
Para obtener los valores de temperatura bibliográficos, sustituimos las composiciones de vapor de
la tabla en la ecuación de segundo grado que obtuvimos anteriormente.
|𝑇𝑒𝑥𝑝 − 𝑇𝑏𝑖𝑏𝑙 | × 100
𝐸𝑅 (%) =
𝑇𝑏𝑖𝑏𝑙
Tabla 8. Errores relativos de las temperaturas de la mezcla ciclohexano-tolueno
T exp (ºC) Tx (ºC) Ty bib (ºC) Error Tx Error Ty
105 104,64 106,07 0,34 1,01
99,7 96,23 98,26 3,61 1,47
95 97,28 97,47 2,34 2,53
92 91,28 93,17 0,79 1,25
88,4 86,00 87,97 2,79 0,48
88,9 89,09 88,90 0,22 0,00
83,5 83,88 85,24 0,46 2,04
82,3 84,43 84,17 2,53 2,23
79,8 82,60 83,38 3,39 4,29
Finalmente, se determinaron los coeficientes de actividad a dilución infinita mediante
extrapolación de las siguientes curvas:
ϒi vs xi(experimental)
2,5
1,5
Yi
gamma 2
1
gamma 1
0,5
0
0,000 0,200 0,400 0,600 0,800 1,000
Xi
13
Figura 8. Coeficientes de actividad frente a composiciones
Extrapolando de la figura, podemos estimar el valor de γ1ꚙ , γ2ꚙ . Con estos valores de coeficiente
de actividad a dilución infinita se pueden estimar A12 y A21, que son los parámetros de iteración
binaria de la ecuación de Van Laar mediante la ecuación (3).
Tabla 9. Coeficientes de actividad a dilución infinita y parámetros de Van Laar
ϒ1ꚙ 1,3
ϒ2ꚙ 2,5
A12 0,2623
A21 0,9162
Rectificación en torres de relleno
Los valores experimentales de temperaturas de equilibrio de las corrientes de destilado y del piso
1 son:
Tabla 10. Temperaturas de las corrientes de piso 1 y destilado.
E1 E2 E3
T (ºC) piso 1 91,3 92,4 94
T (ºC) destilado 79,2 78,9 79
Como cabe esperar, las temperaturas de la corriente de destilado son menores que las del piso 1
ya que, viendo el diagrama de la figura 7, la T es menor para una composición de volátil mayor.
Al igual pasa en la ascendencia de las temperaturas del piso 1 ya que la composición del más
volátil va disminuyendo en cada experimento. Sin embargo, esta tendencia no ocurre en las
temperaturas del destilado y esto se puede deber a una mala calibración del termómetro de esa
corriente.
-Método de McCabe-Thiele
Este método consiste en hallar gráficamente el número de pisos teóricos.
14
A partir de los valores de área que se han obtenido para las muestras de destilado y del piso 1 en
cada uno de los experimentos se pueden calcular las composiciones de las corrientes dividiendo
el área entre la pendiente de la recta de calibrado: área=31135*x o y.
Tabla 11. Composiciones de las muestras de destilado y del piso 1 en cada experimento.
TM02 Piso Area Composición
Columna 1 1-muestra 1 17139 0,550
Columna 1 1-muestra 2 17319 0,556
Columna 1 1-muestra 3 14494 0,466
Columna 1 D-muestra 1 30502 0,980
Columna 1 D-muestra 2 29929 0,961
Columna 1 D-muestra 3 30068 0,966
El valor de composición del piso 1 del tercer experimento no es del todo representativo ya que se
perdió muestra al traspasarla de un frasco a otro.
A continuación, representamos los pisos teóricos en el diagrama de distribución (de los datos
bibliográficos) entre las composiciones de piso 1 y de destilado de cada experimento:
(M+N)min=6.6
Figura 9.1. Número mínimo de pisos teóricos del experimento 1.
15
(M+N)min=4.7
Figura 9.2. Número mínimo de pisos teóricos del experimento 2.
(M+N)min=5.65
Figura 9.3. Número mínimo de pisos teóricos del experimento 3.
16
-Método analítico de Fenske
Este método consiste en hallar analíticamente el número mínimo de pisos teóricos.
El procedimiento de cálculo es el siguiente:
En primer lugar, se obtienen los coeficientes de actividad mediante la ecuación de Van Laar:
Las presiones de saturación se obtienen de la ecuación de Antoine (ecuación 1) para cada
temperatura.
Las volatilidades relativas se calculan con la siguiente ecuación:
Como se puede observar, las volatilidades relativas del destilado son menores que las respectivas
del piso 1 ya que la composición del volátil es mayor en el destilado por lo que la facilidad de
separación de ambos compuestos (lo que representa la volatilidad) es menor.
Además, la composición del volátil en el destilado aumenta para cada experimento por lo que la
tendencia de la volatilidad debería ser decreciente, sin embargo, del experimento 1 al 2 esto no
ocurre ya sea por un error experimental.
En el caso de la tendencia de la volatilidad relativa del piso 1 es la contraria al destilado, esta
aumenta ya que la composición del volátil en cada experimento es menor.
Tabla 12. Procedimiento de cálculo del método de Fenske
Experimentos
E1 E2 E3
Destilado
xD1 0,98 0,961 0,966
xD2 0,02 0,039 0,034
1,00116 1,00405 1,00315
2,22182 2,01944 2,06797
TD (K) 352,35 352,05 352,15
PS1 (mmHg) 725,689 719,102 721,292
PS2 (mmHg) 283,256 280,319 281,295
12 1,1544 1,2738 1,2428
Piso 1
17
xP1 0,55 0,556 0,446
xP2 0,45 0,444 0,554
1,15482 1,15272 1,18915
1,06349 1,06589 1,03267
TP1 (K) 364,45 365,55 367,15
PS1 (mmHg) 1033,156 1065,438 1113,779
PS2 (mmHg) 424,111 439,280 462,115
2,645 2,635 2,805
Tabla 13. Volatilidad media, número de pisos y AEPT para cada experimento.
E1 E2 E3
12m 1,747 1,832 1,867
(M+N)min 6,613 4,922 5,707
AEPT (cm) 7,712 10,362 8,936
La volatilidad relativa media se calcula:
El número mínimo de pisos teóricos se calcula:
Teniendo en cuenta que la altura de relleno de la columna son 51 cm, la altura equivalente de
piso teórico es:
4. CONCLUSIONES:
Se han obtenido experimentalmente unos parámetros de iteración binaria de la ecuación de Van
Laar (A12=0,2623 ; A21=0,9162) semejantes a los valores bibliográficos de Myers (A12=0,1544 ;
A21=0,6931). Los errores relativos de temperatura y composición son bajos para todos los
experimentos excepto para la mezcla 1. Los errores cometidos en esta mezcla se pueden deber a
errores experimentales cometidos durante el transcurso de la práctica, tales como no haber
esperado el suficiente tiempo para alcanzar el equilibrio en el ebullómetro, errores en la medida
de las temperaturas, etc.
Si comparamos el número mínimo de pisos teóricos del experimento 1 para el método de McCabe-
Thiele y el de Fenske obtenemos valores de 6.6 y 6.613 respectivamente, en el experimento 2 se
18
obtienen valores de 4.7 y 4.922 y en el experimento 3 los valores son de 5.65 y 5.707, como se
puede observar, son valores muy parecidos entre sí en los tres experimentos por lo que podemos
concluir que ambos métodos son válidos para hallar el número mínimo de pisos.
19
