UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II

FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo

MÉTODO DE MC CABE – THIELE

El método desarrollado por Mc Cabe y Thiele está basado en la

representación gráfica de las ecuaciones de los balances de materia
en un diagrama de distribución.

El método se desarrollará para una columna provista de un

condensador total y considerando despreciable las pérdidas de calor
fuera del sistema, estableciendo los balances de materia en estado
estacionario para cada sección de la columna. Además se considerará
que al interior de la columna, las corrientes de líquido y de vapor son
saturadas.

Consideraciones adicionales

 Las etapas se numeran en orden creciente y en forma

descendente. QC

1
2 L0 D
3 x0 xD
F 4
zF

B
QR xB

-1-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
 En la Zona de Enriquecimiento:
QC
L flujo molar de líquido
V
V flujo molar del vapor y1

 En la Zona de Empobrecimiento:

L flujo molar de líquido

V flujo molar del vapor 1 Lo D

xo xD
 Las composiciones, referidas al componente más volátil, se

expresan en fracción molar:

n
x = fracción molar de A en el líquido

y = fracción molar de A en el vapor L V

xn yn+1
 Relación de reflujo o relación de reflujo externa se define F
zF
como R = L0 / D

 Las entalpías molares se representan por H, con subíndice L ó Figura Nº 1 : Zona de Enriquecimiento

V, según se trate de un líquido o de un vapor.

Balance de materia total:
HL = entalpía molar de un líquido
V = L + D (1)
HV = entalpía molar de un vapor

Zona de Enriquecimiento
Balance de materia para el componente más volátil:
Balances de materia, en estado estacionario, entre una etapa “n”
cualquiera ubicada en la zona de enriquecimiento, y el extremo
V . yn+1 = L . xn + D . xD (2)
superior de la columna:

-2-
 UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo

Dividiendo entre el flujo de vapor V:

L D
yn+1 = ---- . xn + ---- . xD (3) Línea de operación de la
V V Zona de enriquecimiento

Esta ecuación corresponde a la de una línea recta ( y = m . x + b) de

pendiente L / V y con intercepto en el eje de ordenadas ( D / V ) . x D. xD/(R+1)
y que se conoce como la Ecuación de la Línea de Operación para la Zona
de Enriquecimiento.

Realizando un balance de materia total, alrededor del sistema de xB xD

condensación:
Figura Nº.2: Línea de operación de la Zona de enriquecimiento

V = Lo + D (4)

Conociendo los datos de equilibrio para un sistema dado,

Utilizando la ecuación (4) y la definición de relación de reflujo R = representados de un diagrama de distribución, la línea de operación se
L0/D, para las suposiciones del método, se tiene: puede representar fácilmente dando valores a xn para obtener
valores de yn+1 , tal como se muestra en la figura Nº 2.
L = Lo = R . D (5)

V = (R + 1) . D (6) Así se tienen dos puntos característicos:

xn = xD yn+1 = xD

xn = 0 yn+1 = xD / (R+1)
Sustituyendo las ecuaciones (5) y (6) en la ecuación (4) :
R xD
yn+1 =------ . xn + ------ (7)
R + 1 R + 1 Zona de Empobrecimiento

Se obtiene una ecuación para la línea de operación de la Zona de Balances de materia, en estado estacionario, entre una etapa “m”
Enriquecimiento, análoga a la anterior, y que también corresponde cualquiera ubicada en la zona de empobrecimiento y el extremo
a la de una línea recta. inferior de la columna:

-3-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo

F
zF Esta ecuación corresponde a la de una línea recta (y = m . x + b) de
m
pendiente L / V y con intercepto (-B/V) . xB en el eje de ordenadas y
L V
que se conoce como la Ecuación de la Línea de Operación para la Zona
xm ym+1
de Empobrecimiento.
m+1
Sustituyendo la ecuación (8) en (10) se obtiene una ecuación análoga
a la anterior y que también corresponde a la de una línea recta:
L B
ym+1 = -------- . xm - ------- . xB (11)
L - B L - B
B
QR xB
En este caso se tiene un punto característico:

xm = xB ym+1 = xB

Figura Nº 3: Sección de Empobrecimiento Balances de materia y energía en la bandeja de alimentación

Balance de materia total:
L V
HL f-1 HVf
L = V + B
(8)
F
zF f
Balance de materia para el componente más volátil: HF

L . xm = V . ym+1 + B . xB L V
(9) HLf HV f+1

Dividiendo entre el flujo de vapor V: Balance de materia total, en estado estacionario:

L B F + L + V = L + V (12)
ym+1 = ---- . xm - ---- . xB (10)
V V
(V - V) = (L - L) - F (13)

-4-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo

Balance de energía, considerando despreciables las pérdidas de De la ecuación (19) y relacionando con la ecuación (18) se obtienen
calor al exterior: relaciones para evaluar los flujos de líquido y vapor para zonas
adyacentes:
F . HF + L . HLf-1 + V . HV f+1 = L . HLf + V . HV f (14)
L - L = F . q (20)
Bajo las suposiciones del método de que no existe variación marcada V - V = F . (q - 1) (21)
de entalpías para corrientes de líquido y de vapor que proceden de
etapas adyacentes:
Intersección de las líneas de Operación
HVf = HVf+1  HV (15)

Si se encuentra el punto de intersección de las dos líneas de operación,

HLf-1 = HLf  HL (16)
esto permitirá que se puede graficar la línea de operación de la Zona
de Empobrecimiento en el diagrama de distribución.
Luego el balance de energía se puede expresar de la siguiente
manera:
De las ecuaciones de los balances de materia en cada zona de la
F . HF + L . HL + V . HV = L . HL + V . HV (17) columna:
Zona de Enriquecimiento: V . yn+1 = L . xn + D . xD (22)
(V - V) . HV = (L - L) . HL - F . HF (18)
Zona de Empobrecimiento: V . ym+1 = L . xm - B . xB (23)

Sustituyendo (13) en (18): En el punto de intersección de las líneas de operación:

yn+1 = ym+1 = y
L - L HV - HF xn = xm = x
------- = ----------- = q (19)
F HV - HL Introduciendo estas coordenadas en las ecuaciones de las líneas de
operación:
Esta ecuación representa el Factor q definido como la relación
entre la cantidad de calor necesario por mol para llevar la V . y = L . x + D . xD (24)
alimentación desde su condición hasta la de vapor saturado y el
calor latente molar de vaporización. V . y = L . x - B . xB (25)

-5-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
Restando ambas ecuaciones: de operación de la Zona de Enriquecimiento, permite encontrar el
punto de la línea de operación de la Zona de Empobrecimiento.
(V - V) . y = (L - L) . x - (B . xB + D . xD) (26)

Introduciendo el factor q en estas ecuaciones y de un balance de Condición térmica HF q q/q-1

materia total alrededor de toda la columna: Líquido por debajo del punto de
HF < H L q>1 >1
burbuja
F . (q - 1) . y = F . q . x - F . zF (27) Líquido en su punto de burbuja
HF = HL q=1 
q zF Vapor en su punto de rocío
HF = HV q=0 0
y = ------ . x - ------ (28)
q - 1 q - 1 Vapor por encima de su punto de
HF > HV q<0 1 > q/q-1 > 0
rocío
Mezcla líquido – vapor
HL < HF < HV 1>q>0
Esta ecuación se conoce como Ecuación de la línea q, que
representa el lugar geométrico de los puntos donde es posible la
intersección de las líneas de operación. Desde le punto de vista Líquido saturado
geométrico, corresponde a la ecuación de una línea recta con Líquido subenfriado
pendiente (q / q – 1) y con intercepto (zF / q - 1).
Mezcla líquido - vapor

Un punto característico de esta línea tiene como coordenadas:

Vapor saturado
x = zF y = zF,
Vapor sobrecalentado
mientras que la pendiente dependerá de la condición térmica de la
alimentación: líquido saturado, líquido sub-enfriado, mezcla líquido
y vapor, vapor saturado o vapor sobrecalentado, tal como se
muestra en la figura Nº 6 . zF

Figura Nº 6 : Posición de la línea q en función de la condición

Conociendo la condición térmica de la alimentación, se pueden térmica de la alimentación

encontrar los datos necesarios para representar la línea q en el

diagrama de distribución. La intersección de la línea q con la línea

-6-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
a) Cuando existe libertad de diseño (diseño de una columna nueva),
el plato de alimentación o carga será aquel para el cual el número
total de etapas es el menor, en cuyo caso se denomina Plato óptimo
de alimentación o carga. Apenas se rodea el punto de intersección
de ambas líneas de operación, se cambia de una línea a otra, tal
como se muestra en la figura Nº 8.
xD/(R+1)
b) Cuando la columna se ha diseñado con un plato de alimentación en
una zona determinada. En este caso el cambio de una línea de
operación a otra se efectúa, de ser posible, en dicho plato de
alimentación.
xB zF xD 1.0

0.9

Figura Nº 7: Representación gráfica de las líneas de operación

0.8

Determinación del número de etapas de equilibrio

0.7

0.6

Habiéndose determinado las ecuaciones para las líneas de operación 0.5

y representadas en el diagrama de distribución, se procede a 0.4

determinar el número de etapas de equilibrio. 0.3

0.2

Para tal efecto, se trazan escalones entre las líneas de operación y 0.1

de equilibrio, comenzando por el punto de coordenadas ( x D , xD ) y

finalizando en el punto de coordenadas ( x B , xB ) . El número de
0.0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

escalones trazados corresponderá al número de etapas de equilibrio

incluyendo al recalentador. Figura Nº 8: Determinación del número de etapas de equilibrio – Posición óptima

El cambio de una línea de operación a otra dependerá de cada uno

de las siguientes situaciones:

-7-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
1.0

0.9

0.8

0.7

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

Figura Nº 9: Determinación del número de etapas de equilibrio – Posición adelantada

1.0

0.9

0.8

0.7

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

Figura Nº 10: Determinación del número de etapas de equilibrio – Posición retrasada

-8-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
Algunos autores brindan valores de la relación de reflujo en función
del Rmín, así tenemos:
VARIACIÓN DE LA RELACIÓN DE REFLUJO
R = 1.10 - 1.50 Rmín (Backhurst - Harker)
Para efectuar una separación determinada en una columna de
R = 1.20 - 1.50 Rmín (Treybal)
destilación continua con rectificación se puede elegir
arbitrariamente una relación de reflujo, dentro de ciertos límites R = 1.30 Rmín + 0.30 (Pavlov)
como son la relación de reflujo mínima y la relación de reflujo total.

En la práctica la relación de reflujo varía dentro de este rango:

Relación de reflujo mínimo

Rtotal > R > Rmín La relación de reflujo mínimo se define como la máxima relación que
requerirá un número infinito de etapas para obtener la separación
De esta relación de reflujo dependerán los costos totales, es decir deseada, por debajo de la cual no puede operar la columna de
los costos fijos (costos de equipos) y los costos de operación. Para destilación.
un nuevo proyecto la relación de reflujo que debiera utilizarse ha
de ser la más económica, es decir aquella para la cual los costos Para esta condición corresponden la mínima capacidad de
totales son los menores. Esta se conoce como la relación de reflujo calentamiento en el recalentador y de enfriamiento en el condensador.
óptima, que se determina basándose en un balance económico.
A medida que disminuye la relación de reflujo, la pendiente de la línea
de operación de la zona de enriquecimiento se hace menor y se acerca
hacia la curva de equilibrio; ocurriendo lo mismo con la línea de
Costo
Costo total
anual

operación de la zona de empobrecimiento; por lo que el número de

etapas de equilibrio se incrementa, a la vez que se hace mayor la
Costos de
operación

capacidad de la columna.
Costos fijos

Dependiendo de la forma de la curva de equilibrio, la relación de

reflujo mínima puede determinarse de la siguiente manera:
Rmínimo Róptimo R

-9-
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
a) Curva de equilibrio cóncava xa , ya Coordenadas del punto de intersección de la
línea q con la curva de equilibrio (Pinch point)
En este caso la línea de operación que pasa por el punto de
coordenadas (xD , xD) y por el punto de intersección de la línea q con
la curva de equilibrio, corresponde a la mínima relación de reflujo.
Para un sistema con volatilidad relativa constante, estas
coordenadas se pueden obtener analíticamente por intersección de la
El valor de Rmín se puede evaluar a partir de la intersección de la
ecuación de la curva de equilibrio y de la línea q .
línea de operación con el eje de ordenadas o de su pendiente.

Se obtiene la siguiente ecuación cuadrática:

A partir de la intersección con el eje de ordenadas:

xD xD q ( - 1) xa2 + [ q -  (q – 1) - ( -1) zF ] xa - zF = 0
 = ------- Rmin = ---- - 1
Rmin + 1 
y También es factible determinar Rmín para sistemas con volatilidad
relativa constante por combinación de las ecuaciones de la línea de
operación, la línea q y de la línea de equilibrio, obteniéndose la
siguiente ecuación:
ya

Rmín . zF + q . xD  [ xD (q-1) + zF (Rmín +1) ]
------------------ = ------------------------
Rmín (1-zF) + q (1–xD) (Rmín+1) (1–zF) + (q–1) (1-xD)

xa zF xD
 Si la alimentación es un líquido saturado:
A partir de la pendiente: 1 xD  . (1 - xD)
xD - ya Rmín xD - ya Rmín = ----- ---- - ------------
m = ------ = -------- Rmín = -.------  - 1 zF (1 - zF)
xD - xa Rmín + 1 ya - x a

- 10 -
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
 Si la alimentación es un vapor saturado: xD xD
 = -------- Rmin = ---- - 1
1  . xD (1 - xD)
Rmin + 1 
Rmín = ----- ------ - ------- - 1
 - 1 zF (1 - zF)
Convexa en la Zona Inferior
b) Curva de equilibrio convexa

En este caso no se puede aplicar el método anterior, ya que la línea

de operación cortaría a la curva de equilibrio en otros puntos
intermedios indicando que el vapor procedente de una etapa sería
más rico que el vapor en equilibrio con el líquido proveniente de la
misma etapa, lo cual es imposible.

Convexa en la Zona Superior

En este caso la línea de operación que pasa por el punto de xB zF xD

coordenadas (xD , xD) y que es tangente a la curva de equilibrio,
corresponde a la mínima relación de reflujo. Relación de reflujo total

A medida que la relación de reflujo aumenta, la relación L/G también

aumenta y esto se traduce en un alejamiento de las líneas de operación
de la curva de equilibrio, por lo que el número de etapas de equilibrio
para una separación dada se hace cada vez menor.

Cuando se incrementa la relación de reflujo, de modo tal que tiende a

un valor infinito (R ), entonces la relación L/V = 1, por lo que las
líneas de operación de ambas zonas de la columna coinciden con la línea
diagonal de 45°. Para esta condición, el número de etapas de equilibrio
xD es mínimo.

- 11 -
UNIVERSIDAD NACIONAL DEL CALLAO TRANSFERENCIA DE MASA II
FACULTAD DE INGENIERÍA QUÍMICA Mg. Carlos Angeles Queirolo
1.0

0.9

0.8

0.7

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

xB zF xD
x

El número mínimo de etapas de equilibrio (para reflujo total) se

determina gráficamente en un diagrama de distribución, trazando
escalones entre la línea diagonal y la curva de equilibrio;
comenzando desde el punto de coordenadas (xD , xD) y finalizando
en el punto de coordenadas (xB , xB).

Si para un determinado sistema la volatilidad relativa es constante

o aproximadamente constante, se puede determinar analíticamente
el número de etapas de equilibrio mediante la Ecuación de Fenske.

xD 1 – xB
log
------ . -------
1 - xD xB
Nmín = ----------------------
log 

- 12 -