Distillation

Guy Gauthier ing. Ph.D.

SYS-823 : Été 2010
Équation d’Antoine
 La relation entre la pression de
vapeur et la température est
représentée par:

B
log10 P  A 
S
0

T C
 B
log10 P  A 
S
0

T C
Équation d’Antoine:
 La pression P est exprimée en
Pascal, la température T en degré
Celcius.

 A, B et C sont des constantes.

 Elles peuvent être obtenues avec trois
expériences ou l’on mesure la
température et la pression.
 B
log10 P  A 
S
0

T C
Équation d’Antoine:
 Pour l’eau (Pa,°C) :
 A = 10.23, B=1750 et C = 235.

 Pour du 1,2,4 Trichlorobenzène

(mmHg, °C)
 A = 7.42844, B=2063.53 et C =
235.161.
 Valide de 94.8835 à 461.097 °C.
Volatilité relative
 La volatilité relative est une mesure
de la difficulté avec laquelle deux
composantes chimiques peuvent
être séparées.

 Volatilité de la composante j par

rapport à la composante k est:
yj xj
 jk 
yk xk
Volatilité relative des systèmes binaires

 Se définit par:
y x

 1 y  1 x
 Diagramme du point d’ébullition
115
Mixture de
Vapeur surchauffée
110
deux produits
A et B.
105

Courbe
100 du point
Temperature (°C)

de rosée
Courbe
du point
95 de bulle
d’ébullition
(c)
90
(b)

85

(a)
80

Liquide sous-refroidi
75
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fraction molaire de A
 Diagramme du point d’ébullition
115
Mixture de
Vapeur surchauffée
110
deux produits
A et B.
105
Produit B pur
Courbe
100 du point
Temperature (°C)

de rosée
Courbe
du point
95 de bulle
d’ébullition
(c)
90
(b)

85

(a)
80

Liquide sous-refroidi Produit A pur

75
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fraction molaire de A
 Diagramme du point d’ébullition
115
Mixture de
Vapeur surchauffée
110
deux produits
A et B.
105

Courbe
100 du point
Temperature (°C)

de rosée
Courbe
95
du point
de bulle
Produit à 82.5 %
d’ébullition
de A à 90°C
(c)
90
(b)

Produit
85 à 40%
de A à 90°C
(a)
80

Produit à 40%
Liquide sous-refroidi
de
75
0 A à0.1
80°C 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Fraction molaire de A
Ainsi,…
 … la concentration de produit A à
doublé simplement en augmentant
la température.

 Cela parce que le point d’ébullition de

la composante A est inférieur à celui de
la composante B.
 Si on transfère la vapeur dans un
condenseur, le liquide sera à plus de
80% constitué de produit A.
Courbe de l’équilibre vapeur-liquide
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Courbe de l’équilibre vapeur-liquide
1 Proportion de produit
A dans la vapeur
0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2

Proportion de produit
0.1
A dans le liquide
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Courbe de l’équilibre vapeur-liquide
1
77.76% de produit A
0.9
dans la vapeur
Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2
50% de produit A
0.1 dans le liquide
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Colonne de distillation
 Le produit est
acheminé à un
plateau de la
colonne.
 Identifié « Feed ».
 Par gravité le liquide
s’écoule vers le bas.
 La vapeur monte
vers le haut.
La section de rectification
 Cette section est
localisée au dessus
du plateau recevant
l’alimentation
(Feed).
 Plus on monte plus la
concentration de la
composante légère
augmente.
La section de revaporisation
 Cette section est
localisée au dessous
du plateau recevant
l’alimentation
(Feed).
 Plus on descend plus
la concentration de la
composante légère
diminue.
Le condenseur
 Situé au haut de la
section de
rectification, le
condenseur sert à
condenser la vapeur
arrivant au haut du
distillateur.
Le condenseur
 Une partie du liquide
condensé est
retourné dans le
distillateur (Reflux).
 Le distillat est le
produit léger que
l’on a séparé du
produit lourd.
La bouilloire
 La bouilloire au bas
de la colonne de
distillation sert à
remettre le liquide
sous sa phase
vapeur.
 La vapeur est
retournée au
distillateur (Boilup).
La fraction de queue
 Une partie du liquide
au bas n’est pas
envoyé à la
bouilloire, c’est le
résidu (Bottom).
 Ce liquide contient la
composante lourde
qui a été séparée du
produit léger.
Bilan matière de la composante légère
au plateau i.

 La dynamique est:
d  M i xi 
 Li 1 xi 1  Vi 1 yi 1  Li xi  Vi yi
dt
Li-1xi-1 Viyi

Plateau i

Lixi Vi+1yi+1
Bilan matière de la composante légère
au plateau recevant l’alimentation.

 La dynamique est:
d  M nf xnf 
 Lnf 1 xnf 1  Vnf 1 ynf 1  Lnf xnf
dt
Vnf ynf  Fz F Lnf-1xnf-1 Vnfynf

FzF Plateau nf

Lnfxnf Vnf+1ynf+1
Bilan matière de la composante légère
au condenseur.

 La dynamique est:
d  M 1 x1 
 V2 y2  L1 x1  Dx1
dt
 V2 y2  LD x1  Dx1
D
Distillat
(liquide quittant
la colonne)
Plateau 1
Reflux
(liquide retournant
à la colonne)
L1=LD L1x1 V 2y 2
Bilan matière de la composante légère
à la bouilloire.

 La dynamique est:
d  M NS xNS 
 LNS 1 xNS 1  VNS y NS  BxNS
dt
 Les débits dans la colonne
Distillat
D  Comme rien ne se
perd et rien ne se
Section de
LR=LD VR rectification crée, les débits
doivent balancer,
F(1-qF)
Alimentation comme un bilan
FqF comptable !
Section de
LS VS revaposisation

B
Fraction
de queue
 Les débits dans la colonne
Distillat
D  Au haut de la colonne:
Section de
VR  LR  D
LR=LD VR rectification

F(1-qF)  Au bas de la colonne:

Alimentation

FqF LS  VS  B
Section de
LS VS revaposisation

B
Fraction
de queue
 Les débits dans la colonne
Distillat
D  A l’alimentation:
 Liquide:
Section de
LR=LD VR rectification
LS  LR  FqF
F(1-qF)
Alimentation
 Gaz:
FqF
Section de
VR  VS  F  1  qF 
LS VS revaposisation

B
Fraction
de queue
Design d’un distillateur binaire
 Tout d’abord, il faut déterminer la
fraction molaire de la composante
légère au distillat, à la fraction de
queue ainsi que dans l’alimentation.

 Par exemple:
 Respectivement 0.9, 0.05 et 0.5.
On positionne ces points sur le
diagramme
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Ensuite…
 On trace la ligne « q » de
l’alimentation.
 q est la fraction molaire sous forme
liquide.

 La pente de cette ligne est q/(q-1).

 Supposons que l’alimentation est

totalement liquide (q = 1).
On positionne la droite sur le
diagramme
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Pente q/(q-1)
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Ensuite…
 On trace ensuite la ligne d’opération
du distillateur pour la section de
rectification.
 Au dessus de l’étage recevant
l’alimentation.

 Tracée à partir du point

correspondant au distillat:
 Pente = L/(L+D).
On positionne la droite sur le
diagramme
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Pente L/(D+L)
Vapeur (y)

0.5

0.4

0.3

0.2
xDD/(D+L)
0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Extrêmes
 Si le débit du distillat est nul (D=0),
tout est envoyé au reflux R (Total
reflux condition) et la ligne
d’opération est superposée à la
ligne x=y. Le nombre d’étage est
alors minimal.
Extrêmes
 Le rapport L/D définit le ratio de
reflux R.
 Au ratio de reflux minimal Rmin, la
ligne croise la ligne d’équilibre au
même endroit que la ligne « q » de
l’alimentation.

 Le nombre d’étages est alors infini !

Solution mitoyenne
 … et minimisant les coûts
d’opération:
 C’est d’avoir un ratio de reflux (R =
L/D) égal à environ 1.2 à 1.5 fois le
ratio de reflux minimal Rmin.
On remarque l’intersection entre les
deux lignes
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Intersection
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Ensuite
 On trace une dernière ligne, entre le
point correspondant à la fraction de
queue et l’intersection identifiée
dans la figure précédente.

 La pente de cette droite est: LS/VS.

On trace ainsi la dernière ligne
d’opération
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

Pente = LS/VS
0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Ensuite…
 En partant du point correspondant
au distillat, on trace une ligne
horizontale jusqu’à la ligne
d’équilibre.

 Puis une ligne verticale jusqu’à l’une

des deux lignes récemment
ajoutées.
On remarque l’intersection entre les
deux lignes
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Répétons l’opération autant de fois que
nécessaire
1

0.9

Ligne d'équilibre
0.8

0.7

0.6
Vapeur (y)

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Nombre d’étages
 On dénombre ensuite le nombre
d’étages du distillateur.
 On ajoute généralement des étages
(jusqu’à 10 %) pour prendre en
compte que la colonne ait été
conçue trop petite.
 Sans compter le rendement des étages
qui est de 50 à 70%.
Comptons le nombre d’étages
(il y en a 6)
1
1
0.9

Ligne d'équilibre 1 2
0.8

2 3
0.7

0.6
3 4
Vapeur (y)

0.5

0.4
4 5

0.3

0.2 5 6
0.1

6 7
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Liquide (x)
Renumérotation à cause de
l’environnement de simulation…

 Pour la modélisation, le condenseur

au haut du distillateur est l’étage
#1.
 Ainsi, l’alimentation entre au 3e
étage du distillateur et est considéré
le 4e étage du modèle.
 Dû au fait que le premier élément d’un
vecteur en Matlab est numéroté 1
plutôt que 0.
Exemple avec alpha = 3
 Courbe de l’équilibre vapeur-liquide:
1

0.9

0.8 Distillat
0.7

0.6

0.5
y

0.4 Alimentation
0.3
Résidu
0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Exemple (alpha = 3)

0.9
q=1
0.8

0.7 Pente = 0.375

0.6 Rmin = 0.6
0.5
y

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Exemple
  L = 0.8xD
1

0.9

0.8

0.7
R = 0.8
0.6
De 0.4737
0.5
à 0.5233
y

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Exemple
 Ls/Vs = 1.4938, car croisement à
0.7222; 1

0.9

0.8

0.7

0.6

0.5
y

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Les étages (il y en a 10)
1
1
0.9
2
3
0.8
4
0.7
5
6
0.6

0.5 7
y

0.4
8
0.3

0.2 9
10
0.1

0 Démo sur
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x MATLAB®
Sur MATLAB®

On visait 0.9

On visait 0.05