I) Définitions

1
La distillation est une méthode permettant de séparer
les différents constituants d’un mélange en fonction de
leur température d’ébullition. Elle peut être effectuée
à pression atmosphérique ou sous pression réduite.

•Distillation simple: distillation en une seule étage.

•Distillation fractionnée: colonne de distillation de

plusieurs étages où s’établit un équilibre
thermodynamique liquide-vapeur. Du bas au haut de la
colonne, la vapeur s’enrichit en composé le plus volatil.
 Approche générale pour la résolution des
problèmes de distillation.
2
Les problèmes associés à n’importe quelle opération de
séparations sont de deux types:

1. Problème de conception (design problem): le degré de

séparation désiré est donné et l’équipement nécessaire pour
effectuer cette séparation est conçu.

2. Problème d’évaluation des performances d’un équipement

existant (rating problem); l’équipement dans ce cas est
connu et on désire de calculer ses performances c’est-à-dire
le degré de séparation que l’on peut obtenir pour des
conditions de fonctionnement données. Ce genre de calcul
peut par exemple être effectué pour la détermination des
conditions optimales de fonctionnement d’un appareil de
séparation donné.
 Etapes à suivre pour la conception ou
l’évaluation de performance de l’opération de
distillation:
3

1. Définition du mélange
2. Etablissement des critères de séparations
3. Obtention des données relatives aux substances à
traiter
4. Sélection d’un modèle pour le calcul du nombre
d’étages ou d’unités de transfert
5. Calcul du nombre d’étages ou d’unités de transfert
6. Dimensionnement de la colonne de distillation
 II) Rappels thermodynamiques
II.1) Variables et notation
4 utilisées:

T: température
P: pression

Mélange binaire de deux constituants chimiques A1 et

A2.
x et y fraction molaire, pourcentage en moles.
x1: fraction molaire globale de A1 dans l’ensemble des phases
P1 : la pression de vapeur en A émise par le liquide à la
température T.
P1°(T) : la pression de vapeur saturante de A à la
température du mélange T.
II.2) Définitions:
Mi
Nombres de moles: N i = Wi 5

Fractions molaires:

Mi
Fractions massiques: zi =
k M k

Concentration volumique:

Molalité:
II.3) Quelquesconversions

6
Wi xi
xi → zi  zi =
 k Wk x k
zi
Wi
zi → xi  xi =
k Wk
zk
II.4) Variance d’un
système
La variance d’un système est donnée par :
7

v = c + n −
où v est la variance du système, c est le nombre
d’espèces indépendant constituant le système, n est
le nombre de facteurs physiques (usuellement la
température et la pression), φ est le nombre de
phases.
 II.5) Obtention de diagramme de phases
Diagramme isobare P = constante
x2 = 0
x2 = 0,3 x2 = 0,7
T T

VAPEUR
x2 = 1 8
T (K)

L=V

LIQUIDE

t (s) 1
COURBES ANALYSE THERMIQUE 0
x2
 T II.6) Signification des courbes frontières
T
*
T 1 Vapeur

9
Liq = vap

Liquide *
T
2
x2
0 1
courbe d'ébullition : Cette courbe est la séparation entre
le liquide et le mélange liquide -vapeur.
courbe de rosée : Cette courbe est la séparation entre la
vapeur et le mélange liquide-vapeur
II.7) Visualisation des plateaux théoriques : colonnes à plateaux
0 1 xB
x2l
x1l x3l x4l x5 l
Vapeur ascendante
* Liquide descendant
M T B x5 l y 5
T5
T4
10 y4 x4l
T3
x3l y3
LIQUIDE

x2l
T2
y2
T1 y1
VAPEUR
Aliq + B liq miscibles
*
T A y1 y2 y3 y4 y
T
x1l
T 5 CHAUFFAGE
II.8) Calcul d’equilibre liquide - vapeur

L’équilibre entre phase se traduit11par l’égalité de fugacité:

f i L (x1 , x2 ,..., xk , P, T ) = f iV ( y1 , y2 ,..., yk , P, T ) Où

 P V 
f i ( x1 , x2 ,..., xk , P, T ) = xi   i ( x1 , x2 ,..., xk , T )  Pi (T )  i  exp  sat
L Sat
dP 
 Pi RT

coefficient d’activité
facteur de Poynting tend
vers 1 à des pressions modérées
Et
f iV ( y1 , y2 ,..., yk , P, T ) = yi  i  P
II.8.1) Volatilité et volatilité relative:
La volatilité  j y du
j constituant j est égale au rapport de la
fraction molaire de ce constituant dans la phase vapeur à la
fraction molaire x j de ce même constituant dans la phase
liquide : yj
j =
xj
Cette grandeur est en fait incommode. En effet, elle varie très
rapidement avec la température et ce, à pression totale fixée. Il
est préférable d’utiliser la volatilité relative du constituant j par
rapport au constituant k :  j y j xk
 jk = =
 k x j yk
Pour un solution binaire formée de deux constituants a et b, on
aura l’expression suivante :  ab xa
ya =
1 + ( ab − 1) xa
On choisit en général les constituants de telle manière que :
12
II.8.2) Calcul du point de bulle et du point
de rosée:
➢ Le point de bulle ou point d’ébullition d’un mélange peut être
calculé à partir de l’équation:
n n

 y = x 
i =1
i
i =1
i i =1

où n est le nombre de constituants présents dans le mélange

➢ le point de rosée ou point d’initiation de condensation d’un

mélange, peut être calculé à partir de l’équation:

n n
yi
 xi =
i =1 i =1 i
=1

13
II.8.3) Équilibre liquide-vapeur d'une solution
idéale
Une solution est dite idéale lorsque deux conditions sont
remplies :
➢ Si elle obéit à la loi de
Raoult : l
f j = p j = x j p sat
j
pj est la pression partielle de j dans la solution,

p sat
j est la pression de vapeur saturante de j.

➢ la vapeur peut être considérée comme gaz parfait :

fiv = p j = y j P

P est la pression totale du système.

14
Considérons un mélange binaire, qui forme une solution idéale en
phase liquide, et un gaz parfait en phase vapeur. On cherche à
déterminer les courbes de bulle et de rosée isothermes de ce
mélange. Il suffit pour cela d'écrire l'égalité des potentiels
chimiques ou des fugacités des constituants du mélange entre les
deux phases :

L'équation de la courbe de rosée s'en déduit, il suffit d'éliminer

les compositions en phase liquide entre les relations d'équilibre :

15
Exemple 1:

Les valeurs des pressions de vapeur de n-hexane

(H) et n- octane (O) sont données dans le tableau
ci-dessous :

a) En supposant que le mélange de ces deux

constituants vérifie les lois de Raoult et de
Dalton ; construire les diagrammes (Txy) et
(xy) pour ce système à la pression totale de 1
atm (101 KPa).

b) Lorsqu’un mélange liquide qui contient 30% de n-hexane est chauffé, quelle est la
composition la vapeur qui est en équilibre à la température d’ébullition.

c) La vapeur du mélange obtenu en b) est condensée de la température de rosée à la température

d’ébullition et on sépare le liquide obtenu. Si ce mélange liquide est vaporisé, quelle est la
composition de la vapeur en équilibre qui est formée. Représenter ces transformations b) et c)
dans les diagrammes x-y et Txy.

16
III. Colonne de distillation

17
III.1) Plateau théorique

18
III.1.1) Bilan Matière et enthalpique:

B.M.T
:
B.M.P:

B.E:

Eq. L V:

19
Ces équations sont appliquées pour chaque
étage de la colonne de distillation.

Exemple 2: 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.
Calculer les débits de distillat et de résidu,
ainsi que le taux de récupération de chaque
constituant dans le distillat.
20
III.1.2) HYPOTHESES DE LEWIS
•La chaleur de mélange des constituants à séparer est très
négligeable.
•Les chaleurs molaires de vaporisation sont pratiquement égale.
• la variation des chaleurs sensibles pour le liquide et la vapeur de
plateau à plateau sont négligeable. H n = 0
•La colonne est bien calorifugée de sorte que
Ces hypothèses conduisent aux :
Ln +1 = Ln = ... = L Vn +1 = Vn = ... = V
H L ,n +1 = H L , n = ... = H L H V , n +1 = H V , n = ... = H V

21
II.2) Méthode de Mc Cabe et Thiele pour le
calcul d'une distillation
III.2.1) 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
➢ la composition des produits (distillat et résidu ) sont
imposées dans le cahier des charges.
22
Dans un premier temps, nous choisirons arbitrairement un taux de
reflux ‫ז‬. 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. 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

Pour nous fixer les idées, nous travaillerons sur l'exemple suivant
:
Exemple-3: 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éthanol pur à 98%. On souhaite que le résidu ne
contienne pas plus de 3% De methanol (en fraction molaire).
Nous prendrons pour commencer un taux de reflux de 1,5. La
colonne fonctionne à pression atmosphérique. La volatilité relative
23
est 2,5. calculer le nombre de plateaux théorique.
III.2.2) La méthode de McCabe et Thiele

Nous admettons l'hypothèse de Lewis, selon laquelle le débit

molaire de chaque phase 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.

24
 a) Droite opératoire de rectification

Bilan global pour les deux constituants :

et pour le constituant A :

On obtient donc :

25
 b) Droite opératoire d'épuisement

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

L' = V '+ R
et le bilan sur le constituant le plus volatil donne :

L' x j = y j +1  V ' + xR  R

Ce qui nous aussi à l'équation d'une droite opératoire dans la

section d'épuisement de la colonne :

L' R
y j +1 =  x j -  xR
V' V'
26
 c) Droite d’alimentation

Considérons le plateau d’alimentation et faisons autour de

celui-ci le bilan matière et un bilan thermique

L, HL V, HV

Plateau
F , HF d’alimentation

L’, H’L V’, H’V

Fig. : Plateau d’alimentation

27
Bilan Matière: F + V '+ L = L'+V

Bilan thermique:F  H F + V ' H '

V + L  H L = L ' H L + V  HV
'

Hypothèse de Lewis: H L = H L' = hL et H V = H V' = hV

Enthalpie molaire Enthalpie molaire

liquide Vapeur
hV − H F  hV − H F 
L'− L = F et V − V ' = F 1 − 
hV − hL  hV − hL 

H V = hV − hL = Chaleur latente de vaporisation

28
 hV − H F hV − H F
On pose : q= =
hV − hL H V

On obtient : L'− L = q  F et V − V ' = (1 − q )  F

q=1: liquide d’alimentation bouillant ou saturé

q=0: l’alimentation sous forme vapeur saturé
q < 0: vapeur surchauffée
q >1 : liquide sous bouillant
0 < q < 1: liquide - vapeur

29
Equation de la droite d’alimentation:

(xi,yi) le point d’intersection de deux droites opératoires, on aura

donc:

V  yi = L  xi + D  xD et V ' yi = L' xi - R  xR
La différence membre à membre donne:

(V − V ')  yi = (L − L')  xi + D  xD + R  xR

q z
yi =  xi −
q −1 q −1

30
 1

q =1 q >1
0 <q <1 Fig. : Les q-lines et les
conditions thermiques de
q =0
l’alimentation
q <0

0 Z x 1

31
 d) Problème de conception
Pour fixer les idées, nous travaillerons sur l'exemple suivant

Exemple 4:

On désire distiller un mélange binaire contenant 40% en

moles en constituant le plus volatile. Le distillat et le résidu
doivent contenir des pourcentages molaires en ce constituant
de 90% et 5%, resp. l’alimentation est un liquide saturé de
débit 100 kmol/h. Le taux de reflux à utiliser est 2,625. La
volatilité relative du mélange est de 2,5. En supposant le
condenseur total, le liquide de reflux saturé et le rebouilleur
partiel, déterminer par la méthode de McCabe et Thiele le
nombre de plateaux théoriques de la colonne et la position du
plateau d’alimentation.

32
 d) Problème d’évaluation des performances

Évaluer les performances d’une colonne existante, à savoir

la détermination de:
Il➢existe une relation
la consommation entre le pour
énergétique nombre de un degré
effectuer
plateaux théoriques et le taux de reflux.
de séparation spécifié,
➢ le taux de reflux
2 situations
➢ la position optimale limites de
du plateau fonctionnement
d’alimentation, …

33
 ❑ Taux de reflux minimum

Exemple
Le taux5:
de reflux peut être déterminé à partir de la relation
suivante:
Déterminer le taux de reflux Pente
minimum pour la colonne de
 =
1 − Pente
l’exemple 4 min

Pente = Pente de la droite d’enrichissement.

34
❑ Nombre minimum de plateaux

La colonne fonctionne en circuit fermé, c-à-dire il

n’y a ni courant d’alimentation, ni celui de
soutirage et on a un reflux total.

Les 2 droites opératoires se superposent avec la

droite y=x

35
EQUATION DE FENSKE:
Le nombre de plateaux théoriques peut être déterminé
aussi analytiquement:
 xLD xHD 
ln  
N min + 1 =  xLB xHB 
ln ( LH )

➢Méthode rapide pour déterminer Nmin.

➢Le nombre de plateaux d’une colonne est le double de

Nmin.
36
Exemple 6:

Utiliser l’équation de Fenske pour estimer le nombre de

plateaux théoriques minimum pour effectuer la
séparation indiquée dans l’exemple 4.

37
38
f) Consommation énergétique

✓ Cas de rebouilleur

Le rebouilleur est utilisé pour générer un débit de vapeur égal à

(LR+D). Le constituant le moins volatil (B) étant très largement
majoritaire, une bonne estimation de la puissance consommée
par le rebouilleur est :

 = (L + D )  Δh V
Q B R B

où Δh V
B est l'enthalpie de vaporisation du constituant B.

39
 ✓Cas de condenseur
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 de débit LR:

 = L  Δh V
Q C R A

➢ Dans le cas d'un condenseur total qui condense Le débit

LR+D:
 = (L + D)  Δh V
Q C R A

40
III.3) Efficacité des colonnes à plateaux

L’efficacité de Murphree vapeur, définie par :

yn − yn +1
EMV = * 100
yn − yn +1

41