Distillation

Paternité - Pas d'Utilisation Commerciale - Partage des Conditions Initiales

à l'Identique : [Link]
 Table des matières

Table des matières 3

I - Principe de la distillation 7

A. Schéma d'une distillation simple..........................................................7

B. Terminologie...................................................................................11

C. Amplification de la séparation : relation de Fenske..............................12

II - Calcul d'une colonne simple 17

A. Dénombrement des degrés de liberté.................................................17

B. Bilan de matière..............................................................................19

C. Plateaux adiabatiques : hypothèse de Lewis.......................................20

III - Méthode de Mc Cabe et Thiele pour le calcul d'une

distillation 23

A. Le problème...................................................................................23

B. La méthode de McCabe et Thiele.......................................................24

C. Droite opératoire de rectification.......................................................25

D. Droite opératoire d'épuisement.........................................................27

E. La construction graphique de McCabe et Thiele...................................28

F. Essayez vous-même........................................................................31

IV - Variations autour de la construction de McCabe et Thiele

33

Jacques Schwartzentruber (EMAC)

3
 A. Taux de reflux et nombre d'étages....................................................33

B. Reflux minimum..............................................................................33

C. Condenseur total ou partiel, consommation énergétique.......................34

D. Optimisation économique de la distillation..........................................35

E. Azéotropie......................................................................................39

F. Azéotropie......................................................................................39

G. Efficacité des plateaux.....................................................................39

H. Murphree.......................................................................................41

I. Influence de l'état de la charge..........................................................41

J. Etat de la charge.............................................................................43

Jacques Schwartzentruber (EMAC)

4
 Introduction

Test d'une référence multi : mes références

Jacques Schwartzentruber (EMAC)

5
 Principe de la
I -

I
distillation
Schéma d'une distillation simple 7
Terminologie 11
Amplification de la séparation : relation de Fenske 12

A. Schéma d'une distillation simple

Soit un mélange binaire A , B à séparer ( A étant le constituant le plus volatil).

On peut imaginer de placer ce mélange dans des conditions de pressions et de

température telles qu'il apparaisse deux phases (liquide et vapeur), la vapeur étant
naturellement enrichie en A, le liquide en B.

Jacques Schwartzentruber (EMAC)

7
Principe de la distillation

Image 1 : Distillation éclair (flash)

On peut admettre que le liquide et la vapeur qui quittent le ballon sont à l'équilibre,
et on sait donc déterminer les compositions des deux phases en fonction des
conditions du procédé (température, pression).

Jacques Schwartzentruber (EMAC)

8
 Principe de la distillation

On parle pour un tel procédé de "distillation éclair" ou de "flash".

Sauf si les deux constituants A et B sont de volatilités très différentes, cette

séparation sera insuffisante : la vapeur sera en général loin de contenir du
constituant A pur. Mais il suffit de refroidir la vapeur obtenue pour la condenser
partiellement, ce qui conduit à une nouvelle phase vapeur encore enrichie en A .
De même, en réchauffant le liquide issu du premier flash, on le vaporise
partiellement, et le liquide résiduel est encore enrichi en B ...

On peut ainsi "empiler" des opérations de flash jusqu'à obtenir les produits A et
B avec la pureté voulue.

Image 2 : Purification par flashes successifs

Naturellement, dans cet agencement, la quantité de produit A récupéré en sortie

du dernier flash est très inférieure à celle introduite ; il en va de même pour B .

Jacques Schwartzentruber (EMAC)

9
Principe de la distillation

D'où la nécessité de recycler vers les étages inférieurs le liquide condensé dans les
flashes situés au dessus de l'alimentation, et vers les étages supérieurs la vapeur
produite dans les flashes situés en dessous de l'alimentation. On obtient ainsi le
schéma de principe de la distillation.

Image 3 : Principe de la distillation étagée à contre-courant

Dans cet exemple, on voit que la séparation sélective des constituants d'un
mélange repose sur l'existence d'un gradient de température entre le bas et le
sommet de la colonne. Le plus souvent, ce gradient de température est maintenu
simplement en refroidissant le sommet de la colonne de distillation (condenseur), et
en réchauffant le fond (rebouilleur) : il n'est pas nécessaire d'apporter ou de
soutirer de la chaleur sur chaque flash, l'étagement des températures étant assuré
par le mélange des fluides (vapeur et liquide qui se croisent à chaque étage).

Le refroidissement du condenseur permet de créer une phase liquide, qui descend

dans la colonne, alors que le chauffage du rebouilleur produit une phase vapeur qui
remonte dans la colonne : on parle donc d'opération de séparation à contre-

Jacques Schwartzentruber (EMAC)

10
 Principe de la distillation

courant. La figure ci-dessous montre un exemple (simplifié) de réalisation pratique

d'une colonne de distillation : les flashes sont simplement constitués par des
plateaux (ou étages) où le liquide et la vapeur sont en contact.

Image 4 : Réalisation pratique d'une colonne de distillation

B. Terminologie

Méthode : Par convention

Nous conviendrons de numéroter les étages à partir du sommet de la colonne (le
condenseur étant l'étage 1) vers le bas de la colonne (le rebouilleur est le dernier
étage). Nous notons L i  le débit de liquide quittant le plateau i (qui descend vers
le plateau i1 ) et V  i  le débit de vapeur quittant la plateau i pour remonter au
plateau i−1 .

Les fractions molaires du constituant k sur le plateau i en phases liquide et

i i 
vapeur seront notées respectivement x k et y k .

On appelle distillat le débit de matière d soutiré au condenseur, et résidu celui

soutiré au rebouilleur (noté b ).

Jacques Schwartzentruber (EMAC)

11
Principe de la distillation

On appelle reflux le débit de liquide renvoyé par le condenseur vers le deuxième

étage de la colonne ( r ). Le taux de reflux est le rapport de ce débit au débit de
distillat :

r f =r /d
De même, le taux de rebouillage est le rapport du débit de vapeur rebouillie
(renvoyée vers l'avant dernier étage) au débit de résidu ( = V n /b ).

L'étage sur lequel est introduite la charge est appelé étage d'alimentation. La
portion de colonne située au-dessus de l'étage d'alimentation est la section de
rectification ; la portion de colonne située en dessous de cet étage est la section
d'épuisement . Las fractions molaires globales dans la charge (qui peut être liquide,
vapeur, ou un mélange de liquide et de vapeur) sont notées z k .

Dans une colonne, le rebouilleur est toujours à une pression sensiblement

supérieure au condenseur (pour permettre le bullage de vapeur à travers le
liquide). On fait ainsi apparaître une perte de charge le long de la colonne (que l'on
exprime aussi souvent comme une perte de charge par étage, en supposant un
profil de pression linéaire).

C. Amplification de la séparation : relation de Fenske

Considérons une colonne fonctionnant à reflux total, c'est-à-dire que les débits
d'alimentation et de soutirage sont négligeables par rapport aux débits internes
(reflux et rebouillage).

Jacques Schwartzentruber (EMAC)

12
 Principe de la distillation

Image 5 : Colonne à reflux total

Fondamental
i
Soit α AB la sélectivité de l'étage i (relative au couple de constituants A et B).
La relation de Fenske stipule que la sélectivité de l'ensemble de la colonne est le
produit des sélectivités de tous les étages :

colonne i
AB = AB
i

Jacques Schwartzentruber (EMAC)

13
Principe de la distillation

Démonstration
Soit en effet une colonne à reflux total (figure ci-dessus), et considérons le sous-
système formé des i plateaux supérieurs de la colonne. La sélectivité de ce sous-
système, pour un couple ( AB ) de constituants, est :

1 i 
1 , , i  yA / xA
AB = 1 i 
y B / xB

Ce sous-système n'échange de matière avec l'extérieur que par les débits Li (qui
le quitte) et V i 1 (qui y pénètre). En régime stationnaire, ces deux débits sont
donc égaux, et par conséquent :

L i = V i1 
x iA = yAi1
x iB  = yBi1
La sélectivité du sous-système formé des i1 plateaux supérieurs est :

y 1 i 1
1,... ,i 1 A /xA
 AB =
y 1 i 1
B /xB
i   i1
1,... ,i  x A / x A
=  AB ⋅ i   i1
xB / xB
 i1
1, ..., i y / x  i1
=  AB ⋅ Ai1 Ai1
yB / x B
Cette équation s'obtient à partir de la ligne précédente en introduisant les égalités
des débits. On reconnaît dans la fraction du deuxième membre la sélectivité de
l'étage i1 , et par conséquent :

,i 1 ..., i i1

1,...
AB =1,
AB ⋅AB
d'où, par récurrence, la relation :

colonne i
AB = AB
i

Remarque
On remarquera que la démonstration est valable pour n'importe quel type de
séparation à contre-courant et à reflux total. En particulier, on n'a pas fait
l'hypothèse que les débits de matière quittant chaque étage sont à l'équilibre
thermodynamique.

Cette relation, bien que valable dans un cas particulier seulement (taux de reflux
très grand) permet d'estimer le gain en sélectivité obtenu par une colonne de
distillation par rapport à un seul flash.

Jacques Schwartzentruber (EMAC)

14
 Calcul d'une
II -

II
colonne simple
Dénombrement des degrés de liberté 17
Bilan de matière 19
Plateaux adiabatiques : hypothèse de Lewis 20

A. Dénombrement des degrés de liberté

Nous limiterons notre étude au cas d'une colonne simple : colonne avec une seule
alimentation, pas de soutirage latéral (on ne sort que le distillat en tête et le résidu
en pied), et où seuls le condenseur (étage 1) et le rebouilleur (étage N) sont les
sièges de transferts thermiques avec l'extérieur.

La pression est fixée, et on la suppose à peu près homogène dans toute la colonne.

L'alimentation (débit, composition, température, pression) est imposée.

Jacques Schwartzentruber (EMAC)

15
Calcul d'une colonne simple

Image 6 : Schéma d'une colonne simple

Soit c le nombre de constituants présents dans la charge. Chaque plateau de cette

colonne est décrit par :

Jacques Schwartzentruber (EMAC)

16
 Calcul d'une colonne simple

 c équations d'équilibre (une égalité entre potentiel chimique en phase

liquide et en phase vapeur pour chaque constituant)
 c équations de bilan matière (une par constituant)
 1 équation de bilan enthalpique
Il y a donc N  2 c1 équations descriptives du problème.

Les variables qui interviennent dans ces équations sont :

 les débits molaires partiels de chaque constituant quittant chaque plateau,
en phase liquide et vapeur, soit 2 Nc variables
 les températures de chaque plateau, soit N variables
 la puissance thermique reçue par rebouilleur (dans le bilan enthalpique sur
le rebouilleur) et la puissance thermique cédée par le condenseur (dans le
bilan enthalpique sur le condenseur)
On a donc N  2 c12 variables décrivant le fonctionnement de la colonne. A ces
variables, il convient aussi d'ajouter le nombre de plateaux de la colonne, a priori
indéterminé, et qui, s'il n'apparaît pas explicitement dans l'une des équations,
détermine le nombre d'équations et d'inconnues du système.

On a donc, en tout N  2 c1 équations, reliant N  2 c13 variables : le

nombre de degrés de libertés (variables que l'on peut s'imposer à peu près
arbitrairement) est donc de 3.

Les variables que l'on peut s'imposer dans la conception d'une colonne pourront par
exemple être :
 Une fraction molaire dans le distillat (pureté), une fraction molaire dans le
résidu (pureté), taux de reflux
 taux de récupération de deux constituants dans l'un ou l'autre des débits de
sortie, taux de reflux
 nombre d'étages, taux de reflux et débit de distillat
 quantités de chaleur échangées au condenseur et au rebouilleur, nombre
d'étages...

On remarquera aussi que certaines variables ne peuvent pas être spécifiées

simultanément, par exemple les débits de distillat et de résidu (puisque leur somme
est égale au débit total d'alimentation, connu).

Lors des étapes de conception d'une distillation, on part habituellement de

performances de séparation, par exemple pureté du distillat et du résidu, ou taux
de récupération de deux constituants-clés dans le distillat ou le résidu, ce qui fixe
deux variables. Par le biais des équations de bilan matière sur la colonne, les débits
de distillat et de résidu sont eux aussi déterminés. On choisira souvent comme
troisième variable le taux de reflux, que l'on pourra faire varier de façon à choisir la
configuration de colonne la plus économique possible.

B. Bilan de matière
On souhaite séparer un mélange méthanol-eau, contenant 40% de méthanol en
fraction molaire, pour obtenir au distillat du méthanol pur à 98% et au résidu de
l'eau pure à 97%. Le débit molaire de la charge est de 1 kmol/h.

Jacques Schwartzentruber (EMAC)

17
Calcul d'une colonne simple

Question
Calculer les débits de distillat et de résidu, ainsi que le taux de récupération de
chaque constituant dans le distillat.

C. Plateaux adiabatiques : hypothèse de Lewis

Cette hypothèse permet une simplification considérable du calcul d'une colonne de
distillation. Elle repose sur la constatation que les enthalpies de vaporisation
molaires des différents constituants que l'on sépare dans une colonne sont souvent
assez proches (ou du moins du même ordre de grandeur).

Considérons alors un plateau interne de la colonne (adiabatique), dans une

séparation de deux constituants A et B ( A étant le plus volatil).

Jacques Schwartzentruber (EMAC)

18
 Calcul d'une colonne simple

Image 7 : Début de la construction de Mc Cabe et Thiele pour un condenseur

Sur le plateau  i , le liquide qui descend du plateau  i−1  croise la vapeur

montant du plateau  i1  . Un certain flux de constituant A passe du liquide à la
vapeur sur le plateau; de même, un certain flux de constituant B passe de la

Jacques Schwartzentruber (EMAC)

19
Calcul d'une colonne simple

vapeur au liquide sur ce même plateau.

Or, la vaporisation d'un flux de A nécessite de l'énergie, alors que la condensation

d'un flux de B fournit de l'énergie. Les énergies mises en jeu par ces phénomènes
de vaporisation/condensation sont de loin les termes dominants dans le bilan
énergétique du plateau (les termes de chaleur sensible ou de chaleur de mélange
sont négligeables devant les énergies mises en jeu dans un changement de phase).

Le fait que le plateau soit adiabatique impose donc que l'énergie nécessaire pour
la vaporisation de A est égale à l'énergie cédée par la condensation de B
sur le plateau.

Si nous admettons que les enthalpies de vaporisation molaires des deux

constituants sont proches, on en déduit que le flux de A évaporé est égal au
flux de B condensé sur le plateau.

Si l'on regarde le débit de liquide traversant le plateau, cela signifie que chaque
molécule de A qu'il va perdre sera compensée par une molécule de B . Le débit
molaire total de liquide n'est donc pas modifié à la traversée du plateau (mais sa
composition change). Il en va de même pour le débit de vapeur.
Il en résulte donc :

L i  = L i −1
V i  = V  i 1

Les débits molaires de liquide et de vapeur sont constants dans chaque section de
la colonne (le raisonnement ci-dessus n'est pas valide pour le plateau
d'alimentation).

Jacques Schwartzentruber (EMAC)

20
 Méthode de Mc
III -

III
Cabe et Thiele pour
le calcul d'une
distillation
Le problème 23
La méthode de McCabe et Thiele 24
Droite opératoire de rectification 25
Droite opératoire d'épuisement 27
La construction graphique de McCabe et Thiele 28
Essayez vous-même 31

A. Le problème
Nous nous intéressons à une colonne simple, c'est-à-dire qui possède une seule
alimentation (que nous supposerons pour le moment à l'état de liquide saturé), et
uniquement deux sorties :
 un distillat vapeur, soutiré au condenseur
 un résidu liquide, soutiré au rebouilleur

Nous cherchons à dimensionner une telle colonne, sachant :

 qu'elle doit traiter une charge binaire, composée d'un mélange deux espèces
A et B, de composition connue z A
d b
 la composition des produits (distillat yA et résidu x A ) sont imposées
dans le cahier des charges.

Dans un premier temps, nous choisirons arbitrairement un taux de reflux r f (nous

verrons par la suite comment le choisir de façon raisonnée).

La colonne est supposée isobare : cela revient à négliger la perte de charge entre
plateaux.
Tous les plateaux interne de la colonne sont supposé adiabatiques, e.

L'objectif est de déterminer :

 le nombre de plateaux nécessaire à cette séparation
 la consommation énergétique de cette colonne
 une estimation (grossière) de l'investissement et du coût global de la
séparation

Jacques Schwartzentruber (EMAC)

21
Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Pour nous fixer les idées, nous travaillerons sur l'exemple suivant :
Nous cherchons à distiller un mélange de méthanol (40% en fraction molaire) et
d'eau (60%), de façon à obtenir au distillat du méthol pur à 98%.On souhaite que
le résidu ne contienne pas plus de 3% d'eau (en fraction molaire). Nous prendrons
pour commencer un taux de reflux de 1,5.
La colonne fonctionne à pression atmosphérique.

Simulateur : Construction de McCabe et Thiele

B. La méthode de McCabe et Thiele

Nous admettons 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.

Jacques Schwartzentruber (EMAC)

22
 Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Image 8 : Colonne de distillation

C. Droite opératoire de rectification

Soit (i) un plateau dans la section de rectification.

Jacques Schwartzentruber (EMAC)

23
Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Considérons le système composé des plateaux 1 (condenseur) à i (système entouré

dans le graphique ce dessus). Ce système reçoit de la vapeur à partir du plateau (i
+1) situé juste en dessous de (i), et cède de la matière à travers le distillat et le
liquide quittant le plateau (i). Nous pouvons écrire des bilans matière en régime
permanent sur ce système :

Bilan global pour les deux constituants :

V =Ld
et pour le constituant A :

V y Ai1 =d y Ad L x Ai 

Or, le débit L de liquide est égal au débit de reflux, r, provenant du condenseur

(hypothèse de Lewis). Nous avons donc :
r
V =r d =d  1r f  où r f = est le taux de reflux.
d

On obtient donc :

 i1 rf i  y dA 
y A = x 
1r f A 1r f

relation entre la compositions de la vapeur ( y Ai1 ) et la composition du liquide (

i 
x A ) qui se croisent entre les plateaux (i) et (i+1) : dans un plan ( x A , y A ), sa
représentation est une droite, que nous appellerons droite opératoire de
rectification. Cette droite a les deux propriétés suivantes, qui permettent de la
tracer facilement :
 elle passe par le point de coordonnées ( y Ad  , y Ad  )
d
yA
 son ordonnée à l'origine est
1r f

Simulateur : Construction de McCabe et Thiele

D. Droite opératoire d'épuisement

Nous pouvons effectuer un calcul analogue dans la section d'épuisement de la
colonne (en dessous de l'alimentation) en délimitant un système formé par les
plateaux j à N (rebouilleur) de la colonne.

Le bilan matière global sur ce système donne :

L′ =V ′ b
et le bilan sur le constituant A donne :

Jacques Schwartzentruber (EMAC)

24
 Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

′  j ′  j1   b
L x A =V y A b x A

Ce qui nous conduit aussi à l'équation d'une droite opératoire dans la section
d'épuisement de la colonne :

L′  j  b  b 
y Aj1 = x A − ′ xA
V′ V

Les conditions d'alimentation permettent d'exprimer les débits molaires de liquide

et vapeur dans la section d'épuisement (L' et V') en fonction des débits de liquide et
vapeur dans la section de rectification (L et V). En effet, si l'alimentation est un
liquide saturé, nous pouvons admettre que le liquide amené par la charge rejoint
simplement le débit de liquide descendant dans la colonne, et donc :
L'=L+F
V'=V

Ce qui conduit à une nouvelle forme de l'équation de la droite opératoire

d'épuisement :

d  1r f  b b
y Aj1 = x Aj − x b 
d  1r f  d  1r f  A

On vérifie sans mal que cette droite opératoire :

 passe par le point de coordonnées ( x Ab  , x Ab  )
 croise la droite opératoire de rectification en un point d'abscisse zA

E. La construction graphique de McCabe et Thiele

Nous disposons à ce stade de trois éléments graphiques, dans le plan  x A , y A

 une courbe d'équilibre, qui représente le comportement thermodynamique
du mélange binaire (A, B) à la pression de la colonne. Les plateaux étant des
lieux où, en première approximation, l'équilibre entre phases est atteint,
chaque plateau de la colonne peut être représenté par un point de cette
courbe
 deux droites opératoires, valables chacune dans l'un des sections de la
colonne, qui relient la composition de la vapeur quittant le plateau i+1 en
fonction de celle du liquide quittant le plateau i : les droites opératoires
permettent de passer d'un plateau à l'autre de la colonne.

Jacques Schwartzentruber (EMAC)

25
Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Image 9 : Courbe d'équilibre et droites opératoires

Méthode
La construction est dès lors très facile :
d  1
on part de la composition du distillat, en ordonnée ( y A ou y A , puisqu'il s'agit
bien de la composition de la vapeur au condenseur, plateau 1)
1
la courbe d'équilibre permet d'accéder à x A (composition du liquide dans le
condenseur, ie composition du reflux qui retourne au plateau 2 de la colonne)
2
la droite opératoire de rectification permet alors de déterminer y A en fonction de
1 2
x A connaissant y A , la courbe d'équilibre permet de déterminer la composition
2 
du liquide x A

Jacques Schwartzentruber (EMAC)

26
 Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Image 10 : Construction graphique de McCabe et Thiele

On voit que la construction revient à tracer des « marches » entre la courbe

d'équilibre et la droite opératoire, chaque arête (numérotées 1, 2, sur le schéma)
correspondant à un plateau d'équilibre.

Jacques Schwartzentruber (EMAC)

27
Méthode de Mc Cabe et Thiele pour le calcul d'une distillation

Pour limiter le nombre de plateaux, on va toujours utiliser la droite opératoire la

plus éloignée de la courbe d'équilibre : on change donc de droite opératoire juste
après leur intersection. Ce changement de droite opératoire permet de positionner
le plateau d'alimentation

On arrête la construction lorsque la composition du liquide atteint ou passe en

b
dessous de la composition souhaitée pour le résidu, x A : on a alors déterminé une
configuration de colonne qui permet de répondre au cahier des charges.

F. Essayez vous-même
On souhaite distiller une charge contenant 40% (en fraction molaire) de méthanol
et 60% d'eau. On souhaite un distillat pur à 98% (en méthanol). Le résidu ne doit
pas contenir plus de 3% de méthanol.
Faites la construction pour différents taux de reflux (entre 0,5 et 2).
Vous constaterez :
 qu'il y a un taux de reflux minimal : il faut en effet que les droites
opératoires soient, tout au long de la construction, du même côté de la
courbe d'équilibre, sinon la construction devient impossible. Cela signifie que
si on choisit un taux de reflux inférieur au taux de reflux minimal, il sera
impossible de faire la séparation avec les puretés souhaitées. Lorsqu'on est
juste au taux de reflux minimal, le nombre d'étages devient très grand (en
toute rigueur, il est infini)
 qu'il y a un nombre d'étages minimal : lorsque le taux de reflux devient très
grand, les droites opératoires sont confondues avec la première bissectrice.
Estimez le taux de reflux minimal et le nombre d'étages minimal requis pour cette
séparation. Entre ces deux configurations extrêmes, il existe une infinité de
solutions possibles au problème, qui se distinguent par leur coût.

Simulateur : Construction de McCabe et Thiele

Jacques Schwartzentruber (EMAC)

28
 Variations autour
IV -

IV
de la construction
de McCabe et Thiele
Taux de reflux et nombre d'étages 33
Reflux minimum 33
Condenseur total ou partiel, consommation énergétique 34
Optimisation économique de la distillation 35
Azéotropie 39
Azéotropie 39
Efficacité des plateaux 39
Murphree 41
Influence de l'état de la charge 41
Etat de la charge 43

A. Taux de reflux et nombre d'étages

La relation entre le taux de reflux et le nombre d'étages nécessaires à une
séparation donnée peut être facilement visualisée grâce à la construction de
MacCabe et Thiele.

En effet, augmenter le taux de reflux revient à éloigner les droites opératoires de la

courbe d'équilibre, jusqu'à les confondre avec la première bissectrice. On trouve
alors le nombre minimum d'étages, correspondant à un taux de reflux infini
(colonne à reflux total).

Le taux de reflux minimal, correspondant à un nombre infini d'étages, est obtenu

lorsque la partie utile des droites opératoires vient "toucher" la courbe d'équilibre,
sans la traverser. Cela peut correspondre, selon le mélange traité et les
spécifications de séparation, à deux situations :
 les droites opératoires ont leur intersection sur la courbe d'équilibre
 l'une des droites opératoires est tangente à la courbe d'équilibre

B. Reflux minimum
Déterminez le taux de reflux minimal et le nombre minimal d'étages pour les deux
séparations spécifiées.
Les compositions sont toujours données en fractions molaires.

Jacques Schwartzentruber (EMAC)

29
Variations autour de la construction de McCabe et Thiele

Question 1
On souhaite séparer un mélange méthanol-eau à 40% de méthanol pour obtenir,
au distillat, du méthanol pur à 98% et au résidu de l'eau pure à 97%.

Question 2
On souhaite séparer un mélange acétone-eau à 20% d'acétone pour obtenir au
distillat, de l'acétone à 95% et au résidu de l'eau à 99%.

C. Condenseur total ou partiel, consommation

énergétique
Tout le développement que nous avons fait jusqu'à présent s'applique au cas où le
distillat est vapeur. On dit alors que le condenseur est partiel, puisqu'il ne condense
que le reflux, et pas le distillat. Un condenseur partiel est un étage d'équilibre qui
effectue une séparation (le distillat et le reflux n'ont pas la même composition).

Très souvent, on préfère obtenir un distillat liquide, ce qui facilite son transport et
son stockage. On demande donc au condenseur de liquéfier distillat et reflux, on
parle donc de condenseur total.

Un condenseur total n'est plus un étage de séparation. Par rapport à une colonne
conçue sur la base d'un condenseur partiel, il faut donc rajouter un étage interne à
la colonne si le condenseur est doit être total. En fait, la construction de McCabe et
Thiele est exactement la même, que le condenseur soit total ou partiel, mais elle
commence au deuxième plateau (au lieu du premier) dans le cas d'un condenseur
total.

Que le condenseur soit total ou partiel, le rebouilleur doit toujours générer un débit
r d de vapeur. Le constituant le moins volatil ( B ) étant très largement
majoritaire, une bonne estimation de la puissance consommée par le rebouilleur
est :

 v−l 
Q̇b = rd   h B
 v−l 
où  h B est l'enthalpie de vaporisation du constituant B .

Pour le condenseur, dans lequel le constituant le plus volatil ( A ) est majoritaire,

le puissance thermique à évacuer est :
 dans le cas d'un condenseur partiel, qui ne condense donc que le reflux r :
Q̇c =r  hAv−l 
 dans le cas d'un condenseur total, qui ne condense le débit r d :
 v−l 
Q̇c =r d  h A

D'après l'hypothèse de Lewis, les enthalpies de vaporisation des deux constituants

A et B sont très proches. Cela signifie qu'une colonne à condenseur total
consomme une certaine puissance thermique au rebouilleur (à haut niveau de
température) et rend à peu près la même puissance thermique au condenseur
(mais à température plus basse) : une telle colonne dégrade en fait la puissance
thermique qui la traverse.

Jacques Schwartzentruber (EMAC)

30
 Variations autour de la construction de McCabe et Thiele

D. Optimisation économique de la distillation

La relation entre taux de reflux et nombre d'étages nécessaires est représentée
graphiquement ci-dessous.

Jacques Schwartzentruber (EMAC)

31
Variations autour de la construction de McCabe et Thiele

Image 11 : Relation entre taux de reflux et nombre d'étages

Tous les points de la courbe correspondent à des solutions vérifiant les

spécifications de séparation données. Mais ces solutions diffèrent par leurs coûts.

Jacques Schwartzentruber (EMAC)

32
 Variations autour de la construction de McCabe et Thiele

Le choix d'un taux de reflux augmente la consommation énergétique ; en effet, la

puissance thermique nécessaire au rebouilleur, qui génère un débit r d de
vapeur, est approximativement égale à :

 v−l 
Q̇b = rd   h
où  h v−l  est l'enthalpie de vaporisation du constituant majoritaire dans le
rebouilleur. Par contre, on diminue le nombre d'étages (colonne moins haute), mais
on a besoin de colonnes plus larges (pour permettre à un grand débit de vapeur de
passer).

Il y a donc un compromis à trouver entre coûts d'investissements et coûts de

production. Afin de pouvoir calculer un coût global, combinant investissement (en
EUR) et fonctionnement (en EUR/an), on utilise des coûts annualisés (sans rentrer
dans des détails comptables, il suffit de supposer que l'investissement est payé par
un emprunt remboursé annuellement).

La structure des coûts de séparation a typiquement l'allure suivante :

Jacques Schwartzentruber (EMAC)

33
Variations autour de la construction de McCabe et Thiele

Image 12 : Structure des coûts dans un procédé de distillation

Assez habituellement, on trouve qu'il existe un optimum au voisinage d'un taux de

reflux égal à 1,2 r f , min . La position de cet optimum dépend néanmoins fortement
de la conjoncture économique, puisqu'il s'agit d'un compromis entre des coûts
d'investissement (liés en particulier au prix des matériaux, à la main-d'œuvre) et
des coûts énergétiques (cours du pétrole).

Jacques Schwartzentruber (EMAC)

34
 Variations autour de la construction de McCabe et Thiele

E. Azéotropie
Lorsqu'un mélange binaire présente un azéotrope, la distillation simple ne permet
pas de séparer les corps purs avec n'importe quelle pureté ; tout au plus peut-on
récupérer l'un des corps purs à une extrémité de la colonne, le mélange
azéotropique sortant à l'autre extrémité : la construction de Mc Cabe et Thiele ne
permet pas de "traverser" un azéotrope.

Dans le cas le plus fréquent d'un azéotrope positif, le mélange azéotropique est
plus volatil que chacun des corps purs constitutifs du mélange. Si on distille un tel
mélange, on obtiendra toujours en tête un mélange proche de l'azéotrope, et en
pied l'un ou l'autre des corps purs, selon la composition du mélange d'origine.

F. Azéotropie

Question 1
On dispose d'un mélange 2-propanol-eau à 80% de 2-propanol (en fraction
molaire).
Si on distille ce mélange, quelles compositions peut-on attendre en pied et en tête
de colonne ?

Question 2
Même question pour un mélange contenant 20% de 2-propanol.

G. Efficacité des plateaux

Les plateaux réels d'une colonne de distillation ne permettent jamais d'obtenir
l'équilibre thermodynamique, les temps de séjour étant limités. Il faut donc
disposer de plus de plateaux réels que le nombre de plateaux théoriques estimé par
la méthode de McCabe et Thiele.

Pour tenir compte du fait que les fluides quittant un plateau ne sont pas à
l'équilibre, on affecte à chaque plateau de la colonne une grandeur, appelée
efficacité de Murphree vapeur, définie par :

y Ai − y Ai1 
E MV =  i  , E
y A − y Ai1 
i , E
où y A est la fraction molaire du constituant A dans la vapeur qui serait à
l'équilibre avec le liquide quittant le plateau. On voit que si l'efficacité du plateau
est égal à l'unité, le liquide et la vapeur sortent à l'équilibre (plateau théorique). Si
l'efficacité vaut 0, la concentration de la vapeur en constituant A ne varie pas au
passge du plateau : ce constituant n'est pas transféré d'une phase à l'autre.

On peut en général admettre que l'efficacité de Murphree vapeur est à peu près
constante pour tous les plateaux de la colonne. Elle a souvent une valeur de l'ordre
de 0,7.

La prise en compte d'une efficacité de Murphree dans la construction de McCabe et

Jacques Schwartzentruber (EMAC)

35
Variations autour de la construction de McCabe et Thiele

Thiele revient à utiliser, au lieu de la courbe d'équilibre thermodynamique, une

courbe d'équilibre corrigée, rapprochée de la droite opératoire d'un facteur égal à
l'efficacité de Murphree.

Image 13 : Prise en compte de l'efficacité de Murphree

Jacques Schwartzentruber (EMAC)

36
 Variations autour de la construction de McCabe et Thiele

Simulateur : Construction de McCabe et Thiele

La version du simulateur ci-dessous permet de prendre en compte une efficacité de
Murphree identique pour tous les plateaux de la colonne.

H. Murphree
On sépare un mélange methanol-eau (40% en fraction molaire de méthanol) pour
obtenir du méthanol pur à 98% et de l'eau pure à 97%. En supposant des plateaux
d'équilibre, nous avons déjà déterminé le taux de reflux minimal ( r f , min =0 ,713 )
et le nombre minimum d'étages nécessaires à cette séparation.

Question
Que deviennent ces grandeurs si l'efficacité de Murphree des plateaux est de 0,6 ?

I. Influence de l'état de la charge

Jusqu'à présent, nous avons toujours considéré que la charge était à l'état de
liquide saturé. Supposons maintenant qu'elle soit partiellement vaporisée, avec une
fraction de vapeur (en nombre de moles)  .

Changer l'état de la charge n'a pas d'effet sur les débits dans la section de
rectification de la colonne, qui doit de toute façon être parcourue par un débit de
vapeur égal à r d et un débit de liquide égal à r : la droite opératoire de
rectification n'est donc pas affectée par l'état de la charge, et son équation est :

 i1  i  d 
 rd  y A =r x 1 dx A

Par contre, les débits dans la section d'épuisement seront modifiés.

Ecrivons donc l'équation de la droite opératoire d'épuisement, si L ' et V ' sont les
débits de liquide et de vapeur dans la section d'épuisement :

V ' yAj 1 = L' x Aj −b x Ab 

On admet que la fraction vaporisée de la charge, soit  F , rejoint la vapeur qui

monte dans la colonne, et que la fraction liquide s'ajoute au débit de liquide
s'écoulant vers la section d'épuisement, comme le montre la figure ci-dessous :

Jacques Schwartzentruber (EMAC)

37
Variations autour de la construction de McCabe et Thiele

Image 14 : Répartition de la charge au plateau d'alimentation

On a donc :

L' = r F 1−
r d = V 'F 

Les équations des deux droites opératoires deviennent :

 r d  y Ai1 = r x i1 dx Ad 

 j1 
 rd −F  y A = [ rF  1− ] x1 j  −bx 1b 

Jacques Schwartzentruber (EMAC)

38
 Variations autour de la construction de McCabe et Thiele

et l'intersection de ces deux droites,  x int , y int  vérifie :

 y int = −1 x int z A

Cette dernière équation est celle d'une droite, passant par  z A , z A  et de pente
−1
. On l'appelle "droite d'état thermique de la charge".


Le tracé des droites opératoires lorsque la charge est partiellement vaporisée se fait
donc de la façon suivante :
 d d 
 tracer la droite opératoire de rectification (passe par  xA ,xA  et
 d
x A
0, 
1r f
 déterminer le lieu de l'intersection avec la droite opératoire d'épuisement :
droite passant par  z A , z A  et de pente −1 .

 b  b
 tracer la droite opératoire d'épuisement, issue de  x A , x A  et coupant la
droite opératoire de rectification au point déterminé plus haut.

Simulateur : Construction de McCabe et Thiele

Ce calcul peut facilement être étendu à un liquide sous-refroidi, ou une vapeur

surchauffée : il suffit de remplacer le concept de fraction vaporisée (limité à la
situation où on a deux phases à l'équilibre) par celui de titre vapeur, défini comme
suit :

h−h L ,sat 
=
hV , sat −h L ,sat 
où h est l'enthalpie du mélange considéré, h L ,sat  l'enthalpie de ce fluide porté à
sa température de bulle (liquide saturé) et h V ,sat  l'enthalpie de ce fluide porté à sa
température de rosée (vapeur saturée).

Cette notion de titre vapeur recoupe exactement la notion de fraction vaporisée

pour un mélange liquide-vapeur. Le titre vapeur d'un liquide sous-refroidi est un
nombre négatif et une vapeur surchauffée a un titre vapeur supérieur à l'unité.
Attention, des titres vapeur inférieurs à -0,1 ou supérieurs à 1,1 sont en général
peu réalistes !

J. Etat de la charge
On reprend l'exemple initial de la distillation d'un mélange méthanol-eau (fraction
molaire de méthanol=40%) pour obtenir du méthanol pur à 98% et de l'eau pure à
97%. Le taux de reflux est de 1,5. On envisage deux possibilités :
 charger la colonne avec un mélange liquide saturé
 charger la colonne avec un mélange vaporisé à 50% (en nombres de moles)

Question
Comparez les nombres d'étages nécessaires dans les deux cas, ainsi que la

Jacques Schwartzentruber (EMAC)

39
Variations autour de la construction de McCabe et Thiele

consommation énergétique du rebouilleur. On suppose une efficacité de Murphree

des plateaux de 0,7. Le condenseur est total (on récupère le distillat liquide).

Jacques Schwartzentruber (EMAC)

40