Chapitre 4

Distillation
Méthode McCabe et Thiele
Pour dimensionner une colonne à distiller, il est nécessaire de déterminer le
nombre de plateau théorique de la colonne.

Il y a différents méthodes.

Dans le cas de mélanges binaires, on peut utiliser des techniques simples de

résolution graphique.
La simple est la méthode de McCabe-Thiele introduit dans le 1925.

Cette méthode est basée sur l'hypothèse de Lewis, selon laquelle le débit
molaire de chaque phase (liquide ou vapeur) reste constant d'un plateau à l'autre,
à l'intérieur de chaque section de la colonne.
Nous appellerons L et V ces débits dans la section de rectification, L' et V' ces
débits dans la section d'épuisement.
La colonne à distiller:

•Alimentation
•Rebouilleur
•Condenseur
•Reflux
•Distillat
•Résidu
•Plateaux
Reflux
L
reflux interne R 
V
L
reflux externe R 
D

1 1
L L  LD  D
    1  
V LD  L   L

1 1
L L V  L  V 
      1 
D V L  L   L 
 1. Bilans matière
On définit 3 zones:

1) zone de rectification
2) zone d'épuisement
3) toute l'unité
En régime permanent:

Ce qui entre  Ce qui sort

Sur toute l'unité:

F  DB
z F F  xD D  x B B
z F F  xD D  x B  F  D 

F  z F  xB   D  xD  xB 
 D  z F  xB  B  z F  xD 
 
F  xD  xB  F  xB  xD 

Dans un problème, zF et F étant connus, il faudra calculer les 4 grandeurs

B, D, xD et xB.

Les bilans matière (total et partiel) donnent 2 équations, l'énoncé du problème

doit permettre d'en écrire encore 2 autres.

On pourra rencontrer les spécifications suivantes:

•on spécifie xD et xB

•on spécifie les taux de récupérations d'un composé en tête et en pied:

(xDD) et (xBB)

• on fixe une concentration xD (ou xB) et un taux de récupération DxD (ou B xB).
1.1 Bilan dans la zone de rectification (zone 1)

Vn 1  Ln  D

yn1Vn1  xn Ln  xD D

(1)

L'équation (1) relie la composition du liquide xn avec celle de la vapeur

yn+1: il s'agit donc des compositions de deux courants qui se croisent.

Cette équation définit la ligne opératoire de la zone de rectification. Dans

un diagramme y vs x, on trouvera sur la ligne opératoire les points de
coordonnées (xn, yn+1) qui correspondent à deux courants qui se croisent.
 1.2 Bilan dans la zone d'épuisement (zone 2)

Lm  V m1  B

xm Lm  ym1V m1  xB B 

(2)

On obtient ainsi l'équation de la ligne opératoire de la zone d'épuisement

(points qui correspondent à deux courants qui se croisent)
Remarque: Les lignes opératoires ne sont pas des droites (pas encore),
car à priori, les débits des courants liquide et vapeur peuvent changer
dans la colonne..
2. Hypothèses simplificatrices de Lewis, méthode de Mc Cabe et Thiele.

Les hypothèses de Lewis ( constant molal overflow assumptions) qui

vont permettre de simplifier les expressions des équations des
lignes opératoires, sont les suivantes :

1) la colonne est adiabatique

2) les chaleurs de mélange des deux composés du

binaire sont négligeables

3) les chaleurs molaires de vaporisation de deux

composés sont égales

Ces 3 hypothèses impliquent: donc que lorsqu'une partie de la vapeur se

condense sur un plateau, la chaleur ainsi récupérée n'est pas dissipée
vers l'exterieur puisque le système est adiabatique, elle est donc
disponible pour entraîner la vaporisation d'une partie du liquide. Comme
les chaleurs de mélange sont négligeables et que les chaleurs de
vaporisation sont égales, la condensation d'une mole de vapeur
conduit à la vaporisation d'une mole du liquide.
Hypothèses simplificatrices
condensation d'une mole de vapeur  vaporisation d'une mole du liquide

Ln  Ln1  Constante  L Vn  Vn1  Constante  V

Lm  Lm1  Constante  L V m  V m1  Constante  V

Dans ces conditions, les lignes opératoires sont alors simplement des
droites opératoires dont les équations sont les suivantes:

(3)

(4)

Sur un diagramme y-x, les courants qui se croisent entre les plateaux
sont des points de coordonnées [xn,yn+1] ou [xm,ym+1] qui appartiennent
aux droites opératoires.
Pour tracer une droite il faut soit connaitre deux points,
soit connaître un point et la pente.
Nous allons chercher un point particulier de chacune des droites opératoires
en trouvant l'intersection de ces droites avec la droite y=x.
Zone de rectification:

L D   L D 
y  x    x   xD   x  1     xD 
V  V   V V 
V  L  D D 
 x   x     xD   x  x D
 V  V  V 
La droite opératoire de la zone de rectification:

- passe donc par le point de coordonnée  xD , xD 

- et a une pente égale à de pente L
 
V 
Zone d’épuisement

L B   L B 
y  x    x   xB   x 1      xB 
V  V   V V 

V  L   B  B 
 x   x     xB   x  xB
 V  V  V 

La droite opératoire de la zone d’épuisement:

- passe donc par le point de coordonnée  xB , xB 

- et a une pente égale à de pente  L
 
V 
 3. Principe de la méthode de McCabe et Thiele
Pour résoudre un problème de distillation, il faut satisfaire des équations
d'équilibre liquide-vapeur sur chaque étage et de bilan matière (total et
partiel) entre les étages.

a) les courants qui sortent

d'un étage sont à
l'équilibre et
correspondent donc à
un point de la courbe
d'équilibre

b) les courants qui se

croisent entre les
étages correspondent à
des points appartenant
aux droites opératoires.

La construction
consiste en un tracé en
escalier qui passe
successivement de la
droite opératoire à la
courbe d'équilibre.
 4. Bilan sur le plateau d'alimentation
Sur le plateau d'alimentation, la redistribution du courant d’alimentation
dépend de ses conditions thermodynamiques par rapport à celles du
plateau.

Si l'alimentation est sous forme d'un liquide saturé, tout le débit introduit va
redescendre dans la zone d'épuisement.

Pour une vapeur saturée, elle va monter dans la zone de rectification.

Bilan sur le plateau d'alimentation

Dans le cas d’une alimentation sous forme d’un mélange gaz-liquide, celui-
ci va se répartir dans chacune des zones.
Bilan sur le plateau d'alimentation
Un liquide sous refroidi, non seulement va descendre dans la zone
d'épuisement mais il va aussi entrainer la condensation d'une partie de la
vapeur présente sur le plateau.

Une vapeur surchauffée entraînera la vaporisation d'une partie du liquide.

 Considérons les équations des deux droites opératoires et cherchons
l'intersection de ces deux droites. Il faut résoudre:

(eq1) (eq2)

On multiplie (eq1) par V (eq2) par V, et on retranche les 2 relations

Vy  Lx  DxD V y  Lx  BxB    
 y V  V  x L  L   DxD  BxB 

zF F  xD D  xB B   
 y V  V  x L  L  zF F 
(5)

Ainsi l'intersection des deux droites vérifie l'équation (5) qui est l'équation
d'une troisième droite qu'on appelle la droite de l'alimentation.
On introduit alors le coefficient, q, comme
caractérisant les variations du débit liquide entre les
zones de rectification et d'épuisement.
Cette variation va être quantifiée par rapport au débit
F. On appelle q, la fraction de F qui descend sous
forme liquide

(6)

Bilan sur l’étage d’alimentation F  L V  V  L  L  L  F V V

(7)

(8)
 Avec les équations (6) et (7) l'équation (5) de la droite d'alimentation s'écrit:

Un point particulier de la droite d'alimentation s'obtient en cherchant

l'intersection avec la droite y=x:

 q   zF   q   zF 
y  x  x 
    x   1  
.
 q  1   1  q   q  1   q  1 

 x   zF 
    x  zF
 q 1   q 1 
La droite d'alimentation passe donc par le point de coordonnées  zF , zF 
 q 
et sa pente vaut  
 q  1 
Bilan enthalpique sur l'étage d'alimentation:

(9)

Avec les hypothèses de Lewis, on suppose HL  HL et HV  HV

FH F  LH L  V HV  LH L  VHV

 FH F  H L L  L  HV V  V   
et avec (6) et (7) :

FH F  H L qF  HV (1  q) F
  HV  H F   q  HV  H L 

(10)
 Pour une alimentation qui est sous forme d'un liquide sous-refroidi, le terme
 HV  H F  correspond à l'énergie nécessaire pour réchauffer le liquide
jusqu'à sa température d'ébullition (point de bulle Tb ) puis pour le vaporiser.

 HV  H L  correspond à la chaleur latente de vaporisation du liquide,H Vap

On a donc l'expression

(11)

Pour une vapeur surchauffée, le terme  HV  H F  correspond à l'énergie

:

nécessaire pour réfroidir la vapeur jusqu'à sa température de condensation

(point de rosée Tr ) et on a donc

(12)
5. Construction graphique
Dans un Pb de distillation, on a les variables connues et d’autres que l’on doit calculer
analytiquement ou graphiquement :

Voila les différentes variables :

zF , F , xB , B, xD , D, L, V , L, V , q, nbre d ' étages
La méthode de McCabe et Thiele va nous permettre de chercher le nombre de
plateaux nécessaire à la séparation. Le nombre de plateau découle des
constructions suivantes:

a) on place le point  xD , xD  et on trace la droite opératoire de pente L

V

 xB , xB  et on trace la droite opératoire de pente

.
b) on place le point L
V
c) on place le point  zF , zF  et on trace la droite de l'alimentation  q 
de pente  
 q  1 
d) on fait la construction en escalier en partant de xB ou xD
 Résumé des constructions

Concernant la droite
d’alimentation, on a plusieurs
scénarios
 6. Position optimale du plateau de l'alimentation:

Il est important de connaitre

la position du plateau de
l’alimentation.

Au dessus du plateau de
l’alimentation, on utilise la
droite d’enrichissement

En dessous du plateau de
l’alimentation, on utilisera la
droite d’appauvrissement

La position optimale du
plateau de l’alimentation
correspond à la construction
qui donnera le nbre minimal
de plateau.
Ici c’est le n° 3. No 3: étage optimal
C’est le plateau sous la
marche duquel se trouve le
pt d’intersection des 3 droites
6. Positions non-optimales du plateau de l'alimentation : Deux Exemples

Pour ces constructions on aura un nombre plateaux supérieur au cas optimal

No 4: Trop bas No 2: Trop haut

7. Nombre minimum de plateaux
On peut apprécier la difficulté de faire la séparation, en calculant le nombre
minimum de plateaux. Il s'obtient dans le cas le plus favorable à la
séparation: fonctionnement à reflux total. On ferme les courants F,B,D et
toute la vapeur en tête est retournée dans la colonne:

L  L V V
Les pentes L/V et L /V sont
égales à 1.

Les droites opératoires se

superposent avec la
première bissectrice.
8. Reflux minimum
Reflux minimum dans le cas d'une colonne idéale munie d'un nombre infini
de plateaux.
Graphiquement, un nombre de plateau infini s'observe lorsque le point
intersection des droites opératoires et d'alimentation devient un point de la
courbe d'équilibre (la droite opératoire de la zone de rectification coupe ou
est tangent à la courbe d'équilibre)
 9. Situations complexes

On peut disposer sur une

installation de différents
courants du même
mélange binaire mais
avec des compositions
différentes.
Prévoir plusieurs zones
d’alimentation.

On alimente avec le
courant le plus riche en
produit le plus volatil dans
la partie haute de la
colonne.

Et le courant moins riche

dans la partie basse de la
colonne.
9.1 Deux alimentations
Si F1 est une vapeur saturée:

V  V   F1
L  L
Si F2 est un liquide saturé:

L  L  F2
V V
9.2 Soutirage de liquide
On va pouvoir tenir compte des
soutirages dans la colonne.
Supposons que l’on a un
soutirage à un étage particulier.

Bilan matière au niveau de l’étage.

V V
L  L  P
9.2 Soutirage de liquide
La construction est similaire au
cas avec plusieurs alimentation.

On a un soutirage d’un liquide

saturé donc la pente de la droite
d’alimentation est verticale .

La pente L’/V’ est plus faible que

L/V.
V V
L  L  P
9.3 Injection d'une vapeur d’eau en pied de colonne
(au lieu d’un rebouilleur) pour les mélanges aqueux

XB Y 0
rebouilleur
 9.3 Injection de vapeur d’eau en pied de colonne

Représentation schématique d’une

colonne que l’on opère avec
injection de vapeur d’eau en pied de
colonne.

Bilans matière :

F V  D  B
z F F  0  V  xD D  xB B
10. Rebouilleurs

Rebouilleur total =
toute l'alimentation liquide
du rebouilleur est vaporisée.
(1 entrée liq. et 1 sortie vap.)

Rebouilleur partiel =
formation d'un mélange
liquide-vapeur en équilibre=
équivalent à un étage théorique
(1 entrée et 2 sorties)
10. Condenseur

Le condenseur total condense

tout le courant vapeur sortant de
la colonne

Le condenseur partiel conduit à

la formation d'un mélange
liquide-vapeur en équilibre. Le
condenseur partiel est
équivalent à un étage
théorique.

Dans le décompte des étages à installer dans une colonne à distiller, il

faudra penser aux étages théoriques que représentent un rebouilleur partiel
et/ou un condenseur partiel
Le condenseur et le rebouilleur sont alimentés par des fluides caloporteurs.
Afin d'obtenir le transfert de chaleur, il faut une certaine différence de
températures entre les fluides.
Les tableaux ci-dessous donnent les différences de températures qui sont en
général utilisées en fonction de la nature des fluides caloporteurs utilisés.
Fluide caloporteur du gamme de ΔT (OC)
condenseur

système de réfrigération 3-10 Fluide caloporteur gamme de ΔT (OC)

du
rebouilleur
fluide de procédé 10-20
eau de refroidissement 6-20

vapeur 10-60

fluide de procédé 10-20

huile 20-60

eau bouillante 20-40

air 20-50
11. Reflux optimum : Considérations économiques

Frais totaux
Si R Rminimum

Frais de fonctionnement

Si R alors N Nminimum
Frais fixes

L/D

Valeur optimale Roptimal  1.2 Rminimum

12. Efficacité d'un plateau

Le plateau théorique correspond au cas idéal où

les courants qui en sortent, Vn et Ln, sont à
l'équilibre thermodynamique.

Dans des conditions réelles, les courants Vn et

yn  yn 1
Ln n'atteignent pas toujours cet équilibre.
EMV  *
yn  yn 1
On caractérise cet écart à l'équilibre en
définissant l'efficacité de Murphree en phase
vapeur par:

où yn* est la composition d'une vapeur

hypothétique qui serait en équilibre avec le
liquide Ln (xn)
12. Efficacité d'un plateau

De même, on définit l'efficacité de Murphree en phase liquide par:

xn  xn 1
EML  *
xn  xn 1

où xn* est la composition du liquide hypothétique qui

serait en équilibre avec la vapeur Vn.

L'efficacité de Murphree n'est pas une constante.

Constructions sur la base des
efficacités
- en phase vapeur ac xn  xn 1
- en phase liquide. EML   *
ab xn  xn 1
Le liquide qui sort réellement du
plateau avec la vapeur
correspond au point c tel que

ac
EML 
ab
df y  yn 1
EMV   n*
Une vapeur correspondant de yn  yn 1
au point a serait en
équilibre avec un liquide
correspondant au point b.

df
EMV 
de
 L'efficacité de Murphree est une information précise
mais difficile à calculer.

Le concept d'efficacité globale d'une colonne est

plus simple: