Práctica No.
5
DESTILACIÓN BACTH
1. MODELO DE CÁLCULO Y RESULTADOS:
1.1 Fracción experimental de soluto
Para la práctica de laboratorio se modeló una destilación tipo Bach, donde se
preparó la solución de Etanol-agua con 25ml y 150 ml de cada sustancia de manera
respectiva, posteriormente se inició el calentamiento y se registraron cambios en la
temperatura del destilado, peso y concentración. Para poder seguir el cambio de
concentración con el tiempo se tomaron alícuotas cada 30 segundos del destilado y
se llevaron a un refractómetro.
Los datos encontrados se evidencian en las gráficas 1, 2, y 3 respectivamente.
Grafica 1. Gradiente de temperatura
Grafica 2. Peso del destilado
Grafico 3. Índice de refracción
Los graficos 1, 2 y 3 representan los datos experimentales tomados durante el
ensayo de la destilación batch.
Entonces con el objetivo de calcular las concentraciones del destilado se empleó el
método de diferencias finitas el cual define la concentración en el destilado como:
(Ecuación 1). (Seader, 2011)
𝑊𝑜𝑥𝑜 − 𝑊𝑥
(𝑋𝐷)𝑎𝑣𝑔 = (𝑌𝐷)𝑎𝑣𝑔 = (𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 1)
𝑊𝑜 − 𝑊
𝐷𝑜𝑛𝑑𝑒:
𝐷 = 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑑𝑒𝑠𝑡𝑖𝑙𝑎𝑑𝑜 𝑚𝑜𝑙/ℎ
𝑌 = 𝑌𝐷 = 𝑋𝐷 = 𝑓𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑚𝑜𝑙𝑎𝑟 𝑒𝑛 𝑒𝑙 𝑑𝑒𝑠𝑡𝑖𝑙𝑎𝑑𝑜
𝑊 = 𝑚𝑜𝑙𝑒𝑠 𝑑𝑒 𝑙í𝑞𝑢𝑖𝑑𝑜 𝑟𝑒𝑠𝑡𝑎𝑛𝑡𝑒𝑠
𝑋 = 𝑋𝑤 = 𝑓𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑚𝑜𝑙𝑎𝑟 𝑒𝑛 𝑒𝑙 𝑙𝑖𝑞𝑢𝑖𝑑𝑜
0 = 𝑖𝑛𝑑𝑖𝑐𝑎 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑖𝑛𝑖𝑐𝑖𝑎𝑙
Así mismo, para adaptar el modelo a los datos experimentales es necesario recurrir
a propiedades físicas de la mezcla Etanol-agua empleada, para lo cual se
encontraron modelos que describen la fracción molar liquida de etanol con respecto
al índice de refracción experimental. Teniendo en cuanta la literatura se encontró
que , entre los modelos considerados y condiciones de operación, el modelo que
mejor describe los datos experimentales se muestra a continuación. (Osney Pérez
Ones*, 2009) .
Para R entre 1.3634-1.3658
0.9937 − 0.7254%R
𝑥 = 𝑥𝑤 = (𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 2)
1 + 0.3267%R − 0.7733%R2
𝐷𝑜𝑛𝑑𝑒:
𝑅 = 𝑖𝑛𝑑𝑖𝑐𝑒 𝑑𝑒 𝑟𝑒𝑓𝑟𝑎𝑐𝑐𝑖ó𝑛 𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙
Y para darle valores a 𝑊 .(Ecuación 3)
𝑊𝐷
𝑊= (𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 3)
𝑋𝑤 ∗ 𝑀𝑤𝐸𝑡𝑎𝑛𝑜𝑙 + 𝑋𝐻2𝑂 ∗ 𝑀𝑤𝐻2𝑂
Donde
𝑊𝐷 = 𝑃𝑒𝑠𝑜 𝑒𝑛 𝑒𝑙 𝑑𝑒𝑠𝑡𝑖𝑙𝑎𝑑𝑜(𝑘𝑔)
Reemplazando en la ecuación 1 de diferencias finitas para calcular la fracción en el
destilado, se encontraron las fracciones en el destilado con un ajuste de 0.97, lo que
indica que el modelo de diferencias finitas es confiable y puede detallar los datos
experimentales.(Grafica 4)
Fracción molar en el tiempo con diferencias finitas
0.8
y = 9E-09x2 - 0.0002x + 0.6749
0.6
Xliquido
R² = 0.978
0.4
0.2
0
0 500 1000 1500 2000 2500
Tiempo (s)
Series1 Poly. (Series1)
Grafica 4. Cambio de la fracción en el tiempo con modelo de propiedades fiscas para mezcla Etanol-
agua.
En la gráfica 4 se observan la disminución de la fracción molar con el tiempo y el
ajuste al modelo.
1.2 Fracción teórica de soluto
Para la fracción teórica se planteó la simulación de una torre de destilación con 33
etapas, las cuales se determinaron con las condiciones de los intervalos de tiempo
tenidos en cuenta en el laboratorio a la hora de tomar las mediciones (mediciones
cada 30s), se identificaron entonces tres corrientes en cada etapa, una para el
líquido que entra a la etapa en el momento t y dos corrientes de salida en equilibrio,
una liquida y otra vapor.
Para la modelación de las etapas se plantearon los balances de energía que se
definen igualados a cero como flujos de entrada menos flujos de salida. (Ecuación
4)
𝑁𝑡 = 𝑁𝑣 t+ 𝑁𝐿 t+t (𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 4)
Los equilibrios en la corriente de vapor y la corriente liquidan se salida (Ecuación
5)
𝐿 ∗Il - 𝑉 ∗iV (𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 5)
Se especificó la alimentación y para los cálculos de entalpias, volumen especifico,
fugacidades se utilización los paquetes termodinámicos TermoPR-NRT, finalmente
el sistema se corre teniendo en cuenta el cálculo del calor suministrado, el cual fue
determinado al inicio de la práctica de laboratorio tomando el tiempo en el que la
solución llegó a una temperatura t.
dQ/dt(W) 68.72
t(s) 30
Tabla1. Calor y cambio de temperatura para cada etapa
Se conectaron las corrientes y se corrió el sistema con las 33 etapas. (Grafico 5)
Fracción molar en delta de t para simulacion
de 33 etapas
y = -1E-09x2 - 3E-05x + 0.302
0.4 R² = 1
Xetanol simulacion
0.3
0.2
0.1
0
0 500 1000 1500 2000 2500
Tiempo (s)
Series1 Poly. (Series1)
Grafica 5. Cambios en la Fracción molar en delta de t para simulación de 33 etapas
En el grafico se evidencia entonces el avance de la destilación evidenciado en el
cambio de concentración respecto a deltas de t de 30 segundos, se adaptó
igualmente un polinomio cuadrado que, según el R2, describe en su totalidad el
comportamiento de los datos.
Finalmente, para realizar un análisis comparativo se encontraron los índices de error
entre el modelo teórico y el experimental, los resultados arrojados se muestran a
continuación. (Tabla 2).
Tabla.2 Datos de discrepancia entre modelo experimental con diferencias finitas y modelo teorico a partir de simulacion por etapas
X etanol,DF X etanol, modelo %e X etanol,DF xetanol,modelo %e X etanol,DF xetanol,modelo %e
0.6457 0.3010 53.39 0.5603 0.2759 50.77 0.4389 0.2516 39.64
0.7001 0.2990 57.30 0.5605 0.2738 48.23 0.4020 0.2493 35.54
0.6318 0.2969 53.00 0.5289 0.2716 48.67 0.4168 0.2470 37.72
0.6548 0.2949 54.97 0.5291 0.2694 47.85 0.3868 0.2447 38.32
0.5624 0.2887 48.67 0.5166 0.2672 45.33 0.3966 0.2424 38.93
0.5907 0.2866 51.49 0.5166 0.2628 45.47 0.3850 0.2401 37.62
0.5791 0.2845 50.88 0.4888 0.2606 44.29 0.3710 0.2378 35.90
0.5694 0.2824 50.41 0.4819 0.2583 41.11 0.3576 0.2355 34.15
0.5555 0.2802 49.56 0.4677 0.2561 41.65 0.3384 0.2332 31.09
0.5603 0.2781 50.37 0.4387 0.2538 36.85 0.3323 0.2309 30.53
Tabla 2. Datos de discrepancia entre modelo experimental con diferencias finitas y modelo teórico a
partir de simulación por etapas
En la tabla 2 se muestran los cálculos de error, que ayudan a encontrar la
discrepancia entre los cálculos teóricos y experimentales.
2. ANÁLISIS DE RESULTADOS.
Como se describió en el apartado del modelo de cálculo, se halló la concentración de etanol
en el destilado, siguiendo dos métodos diferentes: mediante diferencias finitas (Seader,
2011), y realizando una simulación por etapas. Con el método de diferencias finitas, el cual
parte de las propiedades medidas experimentalmente, se halló que la concentración del
destilado disminuyo a través del tiempo empezando con xetanol=0,64 y con un valor de
xetanol=0,33 al final del proceso. Mientras que realizando la simulación se reporta una
xetanol=0,30 a los 30 segundos de que el sistema empiece a destilar, y xetanol=0,23 en la etapa
final (etapa 33). Que la fracción final del método de diferencias finitas corresponda con la
fracción molar inicial de la simulación, nos da una idea de la diferencia entre los resultados
obtenidos por ambos métodos (Tabla 2).
Comparando los resultados de ambos métodos se obtuvieron errores relativos porcentuales
por debajo del 53% y el menor error obtenido fue de 30%. De los resultados obtenidos se
observa que la concentración de etanol en el destilado disminuye al avanzar el tiempo, esto
se explica ya que cuando la temperatura de la mezcla aumenta y comienza la ebullición, la
primera sustancia que se evapora y se dirige al condensador, es en mayor proporción el
etanol, debido a que su temperatura crítica es menor que la del agua; por lo tanto, las
primeras cantidades de destilado tendrán una composición mayor en etanol. Pasado cierto
tiempo de iniciada la ebullición, el agua también empieza a evaporarse y al poseer la mezcla
una cantidad de agua mayor que la del etanol, el flujo de vapor de agua condensado que
irá a parar al destilado seguirá aumentando haciendo que la fracción molar del etanol en el
destilado vaya disminuyendo conforme avanza el tiempo.
Aspectos que han podido influir en la desviación de los resultados respecto al
comportamiento real son: imprecisión al medir los índices de refracción en el refractómetro,
ya que durante la práctica se presentaron problemas al hacer esta medición; también cabe
la posibilidad de que el modelo no describa rigurosamente el proceso de destilación batch
que se está desarrollando.
3. CONCLUSIONES
Para una mezcla etanol-agua no es posible hacer una separación total de los dos
componentes por medio de una destilación en batch, ya que además de que parte
del agua se evapora y se condensa para llegar al destilado, la mezcla etanol-agua
forma un azeótropo antes de separarse completamente, por lo que habría que
recurrir a otros métodos de separación más sofisticados.
En la práctica realizada se logró obtener un destilado con una fracción molar de
alrededor de 0.3 .Para aumentar la concentración de etanol en el destilado se podría
aplicar control a la temperatura de la mezcla para que esta se mantuviera por encima
de la temperatura crítica del etanol, y por debajo de la temperatura crítica del agua.
4. REFENCIAS
Geankoplis, J. (1998). Capítulo 9. Secado. En J. Geankoplis, Procesos de transporte y
operaciones unitarias (págs. 595-605). México: CECSA.
Osney Pérez Ones*, J. D. (2009). Evaluación de propiedades físicas de mezclas etanol-agua
(II). Medellin.
Seader, H. y. (2011). Operaciones de separación por etapas de equilibrio en ingeniería
química.