Método McCabe−Thiele para la separación de
Metanol/Agua mediante programación en MATLAB
Andrés Gallego, Yulieth Daza, Stiven Cruz, Hary Correa
Programa de Ingeniería Química, Facultad de Ingenierías y Arquitectura, Universidad de
Pamplona.
Resumen: El método McCabe−Thiele se basa en la destilación y en utilizar un diagrama de
equilibrio para determinar el número de etapas necesarias y la eficiencia de la columna de
destilación. La metodología utilizada se describe en detalle, incluyendo los cálculos y las
aplicaciones utilizadas. Los resultados obtenidos demuestran que el método proporciona una
buena aproximación matemática para la separación de metanol/agua y puede ser utilizado
como una herramienta de aprendizaje efectiva para los estudiantes en el área de la destilación
continua, además se demuestra la eficacia del método McCabe−Thiele en la separación de
metanol/agua y su utilidad en la aplicabilidad educativa.
Palabras claves: Destilación, Separación, Metanol, Agua, Método McCabe−Thiele,
programación en MATLAB.
Introducción
La destilación es un proceso utilizado para separar componentes o sustancias de una mezcla
líquida a través de la ebullición y la condensación selectivas. Este método se basa en que
cada uno de los componentes de la mezcla posee un punto de ebullición diferente, lo que
permite que se evapore y condense a diferentes temperaturas. Para lograr la separación, es
necesario crear una segunda fase, en la que tanto el líquido como el vapor estén presentes y
puedan entrar en contacto mientras fluyen en contracorriente en una columna de platos o
empaquetada. De igual manera, es de gran importancia que los componentes tengan
diferentes volatilidades para que se separen en diferentes niveles de fase. Por último, las
dos fases se deben poder separar mediante la gravedad o por medios mecánicos. Es
importante considerar la naturaleza de la mezcla que se desea separar y la eficiencia del
proceso a la hora de seleccionar el método de destilación adecuado [1].
El sistema de McCabe−Thiele es una técnica empleada en la destilación y permite evaluar
el comportamiento de una mezcla de líquidos binaria. Esta técnica se fundamenta en la
elaboración de un diagrama que representa la correlación entre la composición del líquido
destilado y la composición del vapor que se encuentra en equilibrio con él. El método de
McCabe−Thiele se utiliza para calcular la cantidad de etapas que se necesitan para separar
los componentes de una mezcla binaria y para determinar la eficacia de una columna de
destilación [1].
El equilibrio de metanol−agua es fundamental en el método de McCabe−Thiele para separar
los componentes de una mezcla binaria mediante destilación. El diagrama de
McCabe−Thiele se construye a partir de los datos de equilibrio de la mezcla binaria, que
muestra la relación entre la composición del líquido destilado y la composición del vapor en
equilibrio con él. Este diagrama se utiliza para calcular el número de etapas necesarias para
separar los componentes de la mezcla y para determinar la eficiencia de la columna de
destilación.
1. Metodología
Especificaciones para la separación
Se propuso separar 100 kmol/h de una mezcla de metanol/agua con composición equimolar
en la alimentación cuyas especificaciones se muestran en la Tabla. 1.
Tabla 1. Datos de entrada para el método de McCabe−Thiele [1].
Especificaciones Valor
zF 0.5
xD 0.85
xB 0.082
q 0.5
i, j 3.27
R / Rmin 1.31
Construcción de la curva de equilibrio
Para la elaboración de la curva de equilibrio, se utilizó la siguiente ecuación que utilizó los
datos de la Tabla. 2:
𝛼𝑖,𝑗 𝑥𝑖
𝑦𝑖 = (1)
1 + 𝑥𝑖 (𝛼𝑖,𝑗 − 1)
Esta ecuación depende de la composición en líquido del compuesto más volátil de la mezcla
y su volatilidad relativa media para poder construirse.
Tabla 2. Operaciones comerciales representativas de destilación binaria [1].
Volatilidad Relación entre reflujo
Mezcla binaria
relativa media y reflujo mínimo
Metanol / Agua 3.27 1.31
Construcción de q−line
Dado que el valor de q representa la calidad del alimento en cinco categorías, Tabla. 3. Esto
se tuvo en cuenta más adelante en la construcción del método. Tal y como lo ilustra la Tabla.
1, para este trabajo el alimento es una mezcla parcialmente vaporizada.
Tabla 3. Valores de q para las cinco condiciones de alimentación.
Condición de la alimentación q
Líquido subenfriado >1
Líquido en punto de burbuja 1
Parcialmente vaporizada [0 , 1]
Vapor en punto de rocío 0
Vapor sobrecalentado <0
Una vez identificada la condición de la alimentación, se calculó la línea q mediante:
𝑞 𝑧𝐹
𝑦=( )𝑥 − ( ) (2)
𝑞−1 𝑞−1
Determinación de Rmin y R
En esta sección, se utilizó la línea de la zona de rectificación cuando esta intersecta con la
curva de equilibrio:
𝑅min 1
𝑦=( )𝑥 + ( )𝑥 (3)
𝑅min + 1 𝑅min + 1 𝐷
𝑉𝐵 + 1 1
𝑦=( ) 𝑥 − ( ) 𝑥𝐵 (4)
𝑉𝐵 𝑉𝐵
La ecuación (3) corresponde a la línea de operación de la zona de agotamiento o
empobrecimiento.
Procedimiento en MATLAB
A continuación, se ingresó todos los parámetros de entrada como la composición en la mezcla
del componente más volátil, la composición del mismo en el destilado y fondos para la
resolución del método. Así mismo, se ingresó el valor de q y la volatilidad relativa promedio
para la mezcla problema.
En la línea 17 se inició el trazo de la línea de 45° cuyo título del diagrama se ingresó en la
línea 20. Los valores de xD, zF y xB se ingresaron como líneas no sólidas según sus coordenadas
correspondientes hasta la línea 24. Al tenerse ya la línea de 45° y las composiciones ubicadas en el
diagrama, se construye la curva de equilibrio a través de la ecuación (1) ingresada en la línea 25 como una
función de manejo para declarar implícitamente su variable independiente. Finalmente, se grafica la curva
mediante el uso de “fplot” con un ajuste en ambos ejes de un valor mínimo 0 y un valor máximo 1.
Cuando se observa la ecuación (2), puede apreciarse que su pendiente e intercepto se ingresó
en la línea 31 y 32, respectivamente. Como la ecuación de la alimentación corresponde a la
definición y = mx + b, se agregó de la misma forma en la línea 33. A partir de la línea 34 hasta la línea
47 se ingresaron los condicionales para los casos de la Tabla. 3
Al saberse que el valor de Rmin se calcula cuando se ubica el “pinch point”, se crea una
variable con su nombre en la línea 51 cuyo valor que tomará será el valor de la composición
en el alimento cuando se tenga un líquido saturado. En los demás casos, el valor de x de ese
punto pinch será calculado mediante la intersección de la ecuación (1) con la (2) como se
indicó en la línea 53.
Después de realizarse ello, el valor de x será conocido y por lo tanto, evaluándolo en la
ecuación (1) como se muestra en la línea 55, se podrá calcular su par ordenado. Como ya se
conocen las coordenadas del punto pinch y el punto de destilado, se pude calcular la
pendiente de dicha línea como se observa en la línea 56. De la ecuación de la recta, se despeja
el intercepto ya que se conoce la pendiente y se ingresa en la línea 57. Mediante la ecuación
(3) y resolviendo para el intercepto, se determina el valor de Rmin. Finalmente, se determina
el valor del reflujo real.
La ecuación (3) correspondiente la zona de rectificación se ingresó en la línea 64.
Nuevamente se determina la intersección entre la línea de rectificación y la línea de la
alimentación, por esa razón, en la línea 68 se ingresó “fsolve” con el fin de determinar el
valor de “x_intersect” que convierte en cero la función objetivo. La línea de rectificación y
agotamiento se grafica gracias a los códigos de las líneas 72 y 73, respectivamente. Como se
requiere determinar y trazar las etapas en la rectificación, x se resuelve de la ecuación (1)
para encontrar su valor en el equilibrio (línea 76).
Como ahora se requiere trazar líneas para el conteo de etapas, se crea una nueva variable
ya que dichas líneas poseen coordenadas distintas. Dado que el conteo es de derecha a
izquierda, la primera coordenada (xD, yD) = (xD, xD) se conoce y por eso toman los valores de
las líneas 78 y 79 del código. Dado que, el programa necesita saber cuándo detenerse en el
conteo, se ingresa un “while” hasta la línea 95. Finalmente, en la línea 96 se ingresa el
comando para imprimir la cantidad de etapas en la zona de rectificación, teniendo en cuenta
que la función “ceil” aproxima el valor a la cantidad superior (Ej. 4.03 5).
En este apartado, el código posee la misma lógica que se realizó en la zona de rectificación.
Finalmente, el número de etapas totales se imprime en la línea 124.
2. Resultados
A continuación, se presentan los resultados obtenidos:
Fig. 1. Cálculo del Rmin.
Fig. 2. Diagrama de McCabe−Thiele para el sistema metanol/agua.
Tabla 4. Resultados obtenidos a partir del método de McCabe−Thiele.
Resultados Valor
Rmin 0.7160
R 0.9380
NR 4
NS 5
N 9
3. Análisis de resultados
Al observarse la Fig. 1, puede evidenciarse la correcta definición de la ubicación del
punto pinch, permitiendo así calcular el resultado obtenido del reflujo mínimo de la
Tabla. 4. En la Fig. 2, cuando se trazan las líneas de rectificación y agotamiento, estas
muestran un comportamiento congruente con las definiciones de McCabe y Thiele en su
trabajo. Es claro que las etapas de rectificación y agotamiento poseen un desfase en el
cambio de zona de la columna de destilación al pasar por la línea de alimentación, esto
es debido a que el programa realizado presentó un modelo matemático básico para la
elaboración de las etapas sin buscar una forma de comparación en el plano de las
mismas, pues lo anterior requeriría niveles más grandes de programación y
pensamiento computacional. No obstante, para la separación propuesta se tiene una
buena aproximación matemática y entendimiento del procedimiento McCabe−Thiele con
el fin de ilustrar y mejorar el aprendizaje de los estudiantes de ingeniería química en el
área de la destilación continua.
Es relevante mencionar el comportamiento de la volatilidad relativa en estos cálculos,
pues cuanto más pequeña sea esta, más difícil será la separación, requiriendo así un
mayor número de etapas a una relación de reflujo R/Rmin constante.
4. Conclusiones
El método McCabe−Thiele implementado en MATLAB fue una herramienta efectiva
para la separación de metanol−agua mediante destilación continua. Los resultados
obtenidos mostraron una buena aproximación matemática y un entendimiento claro del
procedimiento McCabe-Thiele. Esto es especialmente relevante para los estudiantes de
ingeniería química, ya que les permite comprender y mejorar su aprendizaje en el área
de la destilación.
A pesar de que el programa utilizado (MATLAB) presentó un modelo matemático
básico, los resultados obtenidos fueron congruentes con las definiciones establecidas por
McCabe y Thiele en su trabajo original. Aunque se demostró un desfase en el cambio de
zona de la columna de destilación al pasar por la línea de alimentación, esto no afectó
significativamente la calidad de la separación propuesta. Sin embargo, sería
interesante explorar en futuras investigaciones la posibilidad de desarrollar un modelo
más avanzado que permita una comparación más precisa en el plano de las etapas de
la columna de destilación.
5. Referencias
[1] J. D. Seader, E. J. Henley, and K. D. Roper, “Distillation of Binary Mixtures,” in
Separation Process Principles, 3rd ed.John Wiley & Sons, 2013, p. 260.
