Unidad 3: Destilación fraccionada

Método de Ponchon y Savarit

Dra. María Azucena Herrera González

Métodos simples para calcular parámetros de
separación en sistemas binarios

• 1. Método de Mc-Cabe Thiele

• 2. Método de Ponchon y Savarit

Condensador
Destilado

Acumulador
enriquecimiento
Sector de

Alimento
agotamiento

Platos
Sector de

Residuo
Caldera
Características del método
• El método de Ponchon y Savarit es más riguroso que el método de Mc-
Cabe Thiele, ya que este método no admite hipótesis de simplificación.
Por lo tanto el método de Ponchon y Savarit toma en cuenta que:

▫ Los caudales molares de vapor y líquido son distintos en cada plato y en

cada sección de la columna.
▫ Variación de las entalpías específicas de ambas fases no es despreciable frente al
calor latente de vaporización
▫ Calor de mezcla no despreciable
▫ No igualdad de calores latentes molares de los componentes de la mezcla. Calor
latente molar de vaporización de la mezcla dependiente de la composición
▫ Perdidas despreciables de calor en el destilador

u El método de Ponchon y Savarit es más utilizado para diseñar torres en

sistema reales, en cambio el método de Mc-Cabe Thiele es una buena
aproximación del sistema y se utiliza como un primer paso en el diseño de
la columna.
Características del método

• Para sistemas reales, donde el calor latente

molar depende de la composición y donde se
presentan considerables efectos caloríficos de
mezcla, los cálculos resultan más complicados
no siendo aplicable el citado método.

• Para sistemas binarios, Ruheman Ponchon y

Savarit desarrollaron un método gráfico para
resolver este tipo de problemas en sistemas
reales, que se basa en el uso del diagrama
entalpía específica - composición.
Método de Ponchon y Savarit
• Datos de partida:
1. Conocer los datos de equilibrio líquido-vapor del sistema

Equilibrio

y=x
Linea auxiliar

x
Método de Ponchon y Savarit
• Datos de partida:
2. Datos entalpía vapor/líquido-composición vapor/líquido
• La figura muestra un diagrama
típico entalpía específica en
Curva de Vapor saturado donde:
HV/y
H L+V • La curva V corresponde a las
temperaturas de condensación
(vapor saturado) frente a la
composición y separa las zonas de
HL/x vapor (V) y líquido más vapor
Curva de Líquido saturado (V+L).

• La curva L corresponde a las

temperaturas de ebullición
(líquido saturado) y separa las
zonas, V+L y L.
y, x
Método de Ponchon y Savarit
• Datos de partida:
2. Datos entalpía vapor/líquido-composición vapor/líquido

Vapor
• La recta AB es una recta de
equilibrio.
B
HV/y
H L+V • Los puntos A y B representan las
composiciones y entalpías del
C líquido y vapor en equilibrio,
HL/x respectivamente.
A
• Un punto C en la zona V+L se
desdoblará en dos fases de
composición y entalpía dadas por
A para el líquido y B para el vapor.
Líquido

y, x
• Datos de partida:

2. Datos entalpía vapor/líquido-composición vapor/líquido

Este diagrama, como todos aquéllos en los
Vapor que se representan propiedades específicas
frente a la composición, cumple la regla
de la palanca, y los puntos representativos
B
HV/y de corrientes suma o diferencia están
H L+V alineados.
C
En efecto: si se mezclan las corrientes A y B
HL/x para dar C, el proceso se rige por el balance
A de materia total, el balance de materia para
el componente más volátil y el balance de
entalpía:
FA
FC
xA
HA xC
Líquido HC
FB
xB
y, x HB
 Método de Ponchon y Savarit
• Datos de partida:
2. Datos entalpía vapor/líquido-composición vapor/líquido

Vapor Balance de materia total:

mA + mB = mC
B
HV/y
H L+V Balance para el componente más volátil:
C mA x A + mB x B = mC x C
HL/x Balance de entalpía:
A
mAH A + mB HB = mC HC
Combinado las tres ecuaciones:

Líquido mA x B − xC HB − HC Ecuación de una

= =
mB xC − x A HC − H A recta que pasa por
los puntos A, B y C
y, x
Desarrollo del método de Ponchon y Savarit
Q’

V Condensador total

D
L xD
L

F
zF
HF
V
V

L Q’’
B Rehervidor total
xB
Método Ponchon y Savarit
Balance global de materia y energía
QC
F = D+B
Fz F = Dx D + Bx B

1
D FHF + QR = DHD + BHB + QC

L xD
FHF = DHD + BHB + QC − QR
HD
F ⎛ Qc ⎞ ⎛ QR ⎞
FHF = D ⎜ HD + ⎟ + B ⎜ HB + ⎟
zF ⎝ D⎠ ⎝ B⎠
HF Cargas térmicas
N V Qc QR
Q' = HD + Q'' = HB −
D B
QR B
xB FHF = DQ'+ BQ"
HB
Sección de enriquecimiento
Balance global de materia:

QC Vn+1 = Ln + D (1)

Balance del componente más volátil:

V1 , y 1
yn+1 Vn+1 = xn Ln + yD D (2)
Lo, x0
D,
! yD Balance de entalpía global
L1, x1
Vn+1 HVn+1 = Ln HL n + D HD + QC (3)
V2 , y 2
"
L2, x2 (Considerando pérdidas despreciables)
Vn , y n #
Ln , xn Donde Q’ es el calor eliminado del
condensador:
Vn+1, yn+1
Q’ = (D HD + QC) / D (4)
 Sección de enriquecimiento
QC A partir de las ecuaciones (3) y (4):

V1 , y 1 Vn+1 HVn+1 - Ln HL n = D Q’ (5) (D y Q’ =ctte)

Sustituyendo el valor de D de la ecuación (1)

Lo, x0 en (2) y (5).
D,
! yD
Despejando Ln/Vn+1 (relación de reflujo
L1, x1
interno) se tiene:
V2 , y 2
"
L2, x2
Ln y D − yn+1 Q'− HVn+1 (6)
Vn , y n # = =
Ln , xn Vn+1 yD − xn Q'− HLn
Vn+1, yn+1
Curva de operación para la zona de
enriquecimiento
 Sección de enriquecimiento
QC Evaluando la ec. (6) y la ec.(1) en n = 0:

V1 , y 1 L0 y D − y1 Q'− H V 1
Lo, x0 D,
= =
V1 y D − x0 Q'− HL0 (7)
yD
Balance de masa total en n = 0:

! V1 = L0 + D (8)
L1, x1
V2 , y 2 Sustituyendo (8) en (7):
"
L2, x2 L0 Q'− H V 1
Vn , y n = (9)
# D H V 1 −HL0
Ln , xn
Entoces, si el reflujo se define como R=L0/D
Vn+1, yn+1
Q'− H V 1
R= (10)
H V 1 −HL0
Representación gráfica de la sección de enriquecimiento

yD Q’ ΔD: Punto de diferencia (yD, Q’)

Hv vs..y
• La ecuación (6) representa una línea recta
en el diagrama H, x/y.
H 3
2 • La línea de operación pasará por los
1
puntos (xn, HLn), (yn+1, Hvn+1) y (yD, Q’).
HL vs. x
x, y • En el diagrama x,y la ec. (6) permite
graficar la curva de operación para la
zona de enriquecimiento
!"#$%"#&"'
y (!')'*!"#

x yD = x0
Sección de agotamiento
Balance global de materia:
LN-3 , xN-3
VN-2, yN- $%"
LN-3 = VN-2 + B (11)
2
LN-2, xN-2
Balance del componente más volátil:
VN-1, yN-1 $%!
xN-3 LN-3 = yN-2 VN-2 + xB B (12)
LN-1, xN-1
VN , y N $ Balance de entalpía global:

LN-3 HLN-3 + QB = VN-2 HVn-2 + BHB (13)

LN , xN
(Considerando pérdidas despreciables)
Vn+1, yn+1
Sea:
B, xB
Q’’ = (B HB - QB) / B (14)
QB
Sección de agotamiento

LN-3 , xN-3 A partir de las ecuaciones (13) y (14) se obtiene:

VN-2, yN- $%"
2
LN-3 HLN-3 - VN-2 HVn-2 = BQ’’ ' ' (15)
LN-2, xN-2
(donde B Q’’ ctte)
VN-1, yN-1 $%!
Sustituyendo el valor de B de la ecuación (11)
LN-1, xN-1
en (12) y (15), y despejando LN-3 / VN-2
VN , y N $
(relación de reflujo interno) se tiene:

LN , xN LN−3 y N−2 − x B HVN−2 − Q''

= = (16)
Vn+1, yn+1 VN−2 x N−3 − x B HLN−3 − Q''
B, xB Curva de operación de la zona de
QB agotamiento
Representación gráfica de la sección de agotamiento

ΔB: Punto de diferencia (xB, Q’’)

Hv vs..y • La ecuación (16) representa una línea
N-3 recta en el diagrama H, x/y.
H • La línea de operación pasará por los (xN-3,
N-2 HLN-3), (yN-2, HvN-2) y (xB, Q’’).

Q’’ HL vs. x • Sin embargo. en el diagrama x,y la ec.

(16) permite graficar la curva de
xB x, y
operación para la zona de
agotamiento

!"#$%"#&"'($%&')'*$%'

xB x
Zona de alimentación
Líquido saturado: Vapor saturado:
q=1 q=0
HF = HFL HF = HFV

F, z NF Mezcla Líquido-Vapor: Líquido subenfriado:

HF
0< q <1 q >1
(
HF = HFV − q HFV − HFL ) (
HF = HFV − q HFV − HFL )
Vapor sobrecalentado:
HFV − HF q <1
q=
H −H
( )
V L
F F HF = HFV − q HFV − HFL
Despejando la entalpía en el
alimento, tal que:

(
HF = HFV − q HFV − HFL )
 Zona de alimentación
18000
16000
14000
Entalpía (cal/mol)

12000
10000
8000
6000
4000
2000
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x/y
Representación gráfica de la columna de destilación
yD Q’
Hv vs. y

La ec. (20) representa una línea recta en el diagrama

xB Hxy la cual pasa por los puntos (zF, HF), (xB, Q’’) y
(yD, Q’)
zF HL vs. x
Q’’
La construcción para los platos de toda la columna
será la siguiente:
1. Ubicar yD en la curva de vapor saturado
2. Buscar la composición de equilibro del líquido
x en el diabrama de x vs y
3. Ubicar x en el diagrama H, x,y
4.Unir el punto yD con el x encontrado
5.Trazar una liena de Q’ a la x encontrada
xB zF yD = x 0
Número mínimo de etapas (Nmin) y relación de reflujo
mínimo (Rmin)
yD
Nmin—>Rtotal=∞ Hv vs. y Hv vs. y
DDm
Q’=H1V+R(H1V-H0L) XB

Q’=∞ Q’’=-∞
HL vs. x HL vs. x xD
zF

xW yD xW zF yD = x 0

Número mínimo de etapas a reflujo total Relación de reflujo mínimo