Université Mohammed V

Ecole Mohammadia d’Ingénieurs

Département Génie des Procédés

Cours Opérations Unitaires II :

Absorption & Stripping

Enseignant :
Y. ZALIM
 Table des matières

I- Introduction générale ..........................................................................................................2

II- Solubilité à l’équilibre des gaz dans les liquides .................................................................4
II-1- Systèmes à deux composants .............................................................................................4
II-2- Système multi-composants .................................................................................................5
II-3- Solutions liquides idéales ....................................................................................................5
II-4- Solutions non-idéales ..........................................................................................................7
II-5- Choix de solvant pour l’absorption .....................................................................................7
III- Transfert de matière d’un constituant ; bilan de matière ..................................................9
III-1- Bilan de matière pour un écoulement à contre-courant ...................................................9
III-2- Rapport minimale liquide-gaz pour un absorbeur .......................................................... 13

III-3- Bilan de matière pour un écoulement en co-courant ..................................................... 15

IV- Colonnes à garnissage ....................................................................................................... 16

IV-1- Généralités ..................................................................................................................... 16
IV-2- Le garnissage .................................................................................................................. 18
IV-3- Support du garnissage .................................................................................................... 22
IV-4- Distributeur du liquide .................................................................................................... 22
V- Hydrodynamique des colonnes à garnissage ................................................................... 24
V-1- Pertes de charge pour un écoulement monophasé ........................................................ 24
V-2- Pertes de charge pour un écoulement biphasé ............................................................... 25
V-3- Ecoulement à contre-courant du liquide et du gaz à travers les garnissages ................. 25
V-4- Choix du garnissage et calcul du diamètre ...................................................................... 27

VI- Transfert de matière dans les colonnes à garnissage ...................................................... 33

VI-1- Rétention liquide ............................................................................................................ 34
VI-2- Coefficients de transfert de matière............................................................................... 34
VI-3- Cinétique de transfert de matière et calcul de la hauteur d’une colonne à garnissage 40

1
 I- Introduction générale

La purification des matières premières, la séparation des produits désirables des produits
indésirables, le traitement des effluents gazeux et liquides pour l'élimination des composants
toxiques sont certains des problèmes courants fréquemment rencontrés par les ingénieurs
procédés travaillant dans les industries chimiques et les industries connexes.
Certains mélanges peuvent être séparés par des moyens purement mécaniques, par
exemple : séparation des particules en suspension des gaz par des cyclones ou des filtres à
manches ; la séparation des suspensions solides provenant de liquides par filtration,
décantation ou centrifugation ; etc. Mais il existe de nombreux autres mélanges et solutions
dont les composants ne peuvent pas être séparés par des méthodes purement mécaniques,
par exemple : l'élimination de l'ammoniac d'un mélange air-ammoniac, la séparation de
l'hydrogène et de l'azote dans une usine d'ammoniac, l'élimination du dioxyde de carbone et
de l'eau à partir du gaz naturel et ainsi de suite. Dans de tels cas, un groupe d’opérations
basées sur le transfert de matière, c'est-à-dire la capacité de certains composants à se
déplacer à l'échelle moléculaire dans une phase ou d'une phase à l'autre sous l'influence d'un
gradient de concentration, ont été jugés très utiles. L’absorption est un exemple de ces
opérations.
L’absorption, appelé aussi lavage du gaz par solution liquide, est une opération dans laquelle
un liquide set mis en contact avec un mélange gazeux pour dissoudre préférentiellement un
ou plusieurs constituants de mélange gazeux dans le liquide. Parfois, un soluté est récupéré à
partir d’un liquide en mettant un gaz inerte en contact avec le liquide ; cette opération
constitue l’inverse de l’absorption, et elle est appelée Stripping du gaz.

Exemple d’utilisations industrielles de l’absorption

• Le sous-produit gazeux d’une chaudière est lavé par l’eau pour absorber l’ammoniaque
(NH3) et par une huile pour récupérer le Benzène et le Toluène
• Le sulfure d'hydrogène H2S présent dans un gaz de chaudière ou des hydrocarbures
naturelles peuvent être récupéré par absorption dans des solutions alcalines
• Benzène et Toluène sont récupérés à partir de l’huile d’absorption par un stripping à
la vapeur d’eau

Les bases théoriques de l’absorption et du stripping sont les mêmes et évoquent le transfert
d’un ou plusieurs constituant d’une phase (gaz ou liquide) à une autre phase (liquide ou gaz.
Alors, les deux opérations seront traitées ensemble dans ce cours.
Une opération d’absorption gaz-liquide industrielle est, en général, effectuée en continu
dans une colonne verticale où les fluides circulent à contre-courant. Les principaux types de
contacteurs gaz liquide sont :

2
 ✓ Colonnes à plateau et à garnissage
✓ Colonne à bulles, cuve agitée, jet immergé
✓ Colonnes à film et à pulvérisation, venturi, éjecteur
La colonne à garnissage représente le plus simple et le plus traditionnel des contacteurs
gaz-liquide et ainsi le plus utilisé dans les opérations d’absorption. Dans ce cours seules les
colonnes à garnissage seront traitées.

Sortie du gaz

Entrée du liquide

Garnissage

Entrée du gaz

Sortie du liquide

Figure 1.1 : Colonne à garnissage pour

absorption ou stripping

Le diamètre de la colonne à garnissage dépend des quantités de gaz et du liquide mis en

jeu, leurs propriétés et le rapport de leurs débits. La hauteur de la colonne, et par la suite le
volume du garnissage, dépend des variations de concentrations et du taux de transfert de
matière par unité de volume du garnissage à réaliser. Le calcul de la hauteur du garnissage
dépend donc du calcul des bilans de matière, du facteur de potentialité et des coefficients de
transfert de matière.

3
 L’absorption se base sur la solubilité des constituants à transférer, d’une part sa solubilité
dans le mélange gazeux initial et d’autre part dans le liquide ou il va être transférer. Le chapitre
suivant traite la solubilité des gaz dans les liquides.

II- Solubilité à l'équilibre des gaz dans les liquides

La vitesse à laquelle un constituant gazeux d'un mélange va se dissoudre (être absorbé)
dans un liquide absorbant dépend de l'écart d'équilibre existant entre les deux phases. Par
conséquence, il est nécessaire de prendre en compte les caractéristiques d'équilibre des
systèmes gaz-liquide.

II-1. Systèmes à deux composants

Si une quantité d’un gaz pur est mis en contact avec un liquide relativement non volatil
pendant suffisamment de temps, un équilibre sera établi. La concentration résultante de gaz
dissous dans le liquide est la solubilité du gaz à la température et à la pression du travail. La
figure 2.1 présente la solubilité de certains gaz dans l’eau :
𝑝* : Pression partielle du soluté dans le gaz, mmHg

𝑥 : Fraction molaire du soluté dans le liquide

Figure 2.1 : solubilité des gaz dans l’eau

4
Remarques :

• A température constante, la solubilité augmente avec la pression. Exemple courbe A

• Différents gaz donnent différentes solubilités dans le même liquide. Courbe C et F
• Si la pression d'équilibre d'un gaz à une concentration de liquide donnée est élevée,
comme dans le cas de la courbe B, le gaz est dit relativement insoluble dans le liquide,
tandis que si elle est basse, comme pour la courbe C, la solubilité est dite élevée.
• Dans la plupart des cas, la solubilité d’un gaz diminue quand la température augmente.
Exemple courbe A et F
• Au point d’ébullition du solvant, la solubilité est presque nulle à condition que sa
pression de vapeur soit inférieure à celle du gaz ou de la vapeur du soluté.

II-2. Systèmes multi-composants

Si un mélange de gaz est mis en contact avec un liquide, dans certaines conditions la
solubilité à l'équilibre de chaque gaz sera indépendante des autres, à condition toutefois que
l'équilibre soit décrit en termes de pressions partielles dans le mélange de gaz. Si tous les
composants du gaz sauf un, sont pratiquement insolubles, leurs concentrations dans le liquide
seront si faibles qu'elles ne pourront pas influer sur la solubilité du composant relativement
soluble, et la généralisation s'applique.
Par exemple, la courbe A (Fig. 2.1) décrira également la solubilité de l'ammoniac dans l'eau
lorsque l'ammoniac est dilué avec de l'air, l'air étant insoluble dans l'eau, à condition que
l'ordonnée de la courbe soit considérée comme la pression partielle de l'ammoniac dans le
mélange de gaz.
Si plusieurs composants du mélange sont sensiblement solubles, la généralisation ne sera
applicable que si les gaz solutés sont indifférents à la nature du liquide, ce qui ne sera le cas
que pour les solutions idéales. Par exemple, un mélange de gaz de propane et de butane se
dissoudra dans une huile de paraffine non volatile indépendamment, car les solutions qui en
résultent sont essentiellement idéales. D'autre part, la présence de méthylamine peut influer
sur la solubilité de l'ammoniac dans l'eau, car les solutions résultantes de ces gaz ne sont pas
idéales. La solubilité d'un gaz sera également influencée par la présence d'un soluté non volatil
dans le liquide, tel qu'une solution de sel dans l'eau, lorsque ces solutions sont non idéales.

II-3. Solutions liquides idéales

Lorsque la phase liquide peut être considérée comme idéale, la pression partielle du gaz à
l'équilibre peut être déterminer sans avoir recours à l’expérimentation.
Une solution idéale est définie par quatre caractéristiques significatives essentielles,
toutes interdépendantes :

5
 1- Les forces moyennes d'attraction et de répulsion intermoléculaires dans la solution
restent constantes lors du mélange des constituants
2- Le volume de la solution varie linéairement avec la composition
3- Il n'y a ni absorption ni dégagement de chaleur lors du mélange des constituants. Pour
les gaz se dissolvant dans des liquides, toutefois, ce critère ne devrait pas inclure la
chaleur de condensation du gaz à l'état liquide.
4- La pression de vapeur totale de la solution varie linéairement avec la composition
exprimée en fractions molaires.

En réalité, il n'y a pas de solutions idéales et les mélanges actuels ne considèrent l'idéalité
que comme une limite. L'idéalité exigerait que les molécules des constituants soient similaires
en taille, structure et nature chimique, et l'approche la plus proche d'une telle condition est
peut-être illustrée par les solutions d'isomères optiques de composés organiques. Toutefois,
dans la pratique, de nombreuses solutions sont très proches de solutions idéales que pour des
fins d'ingénierie elles peuvent être considérées idéales. Les membres adjacents ou presque
adjacents d'une série homologue de composés organiques entrent particulièrement dans
cette catégorie. Ainsi, par exemple, des solutions de benzène dans du toluène, d'alcool
éthylique et d'alcool propylique, ou des gaz d'hydrocarbures paraffiniques dans des huiles de
paraffine peuvent normalement être considérées comme des solutions idéales.

Si un mélange gazeux, en équilibre avec une solution liquide idéale, suit la loi des gaz
parfaits ; la pression partielle 𝒑
̅ * d'un soluté gaz A est égale au produit de sa pression de
vapeur p à la même température et de sa fraction molaire dans la solution liquide 𝒙. C'est la
loi de Raoult. (* c’est pour indiquer l’équilibre)

̅∗ = 𝒑. 𝒙
𝒑 (2.1)

La nature du solvant liquide n’est considérée que dans la mesure où elle établit l’idéalité de
la solution ; la solubilité d’un gaz dans une solution idéale est donc toujours la même.
Pour des pressions totales supérieures à celles auxquelles s'applique la loi des gaz parfaits,
la loi de Raoult peut souvent être utilisée en substituant les pressions par les fugacités.

6
 II-4. Solutions non-idéales
Pour les solutions liquides non idéales, la loi de Raoult donnera des résultats hautement
incorrects. La ligne D (Fig. 2.1), par exemple, est la pression partielle calculée de l'ammoniac
en équilibre avec des solutions d’eau à 10 ° C en supposant que la loi de Raoult soit applicable,
mais on voit clairement que le résultat est incorrect. D'autre part, on voit que la droite E
représente très bien les valeurs de la solubilité de l'ammoniac dans l'eau à 10°C jusqu'à des
fractions molaires de 0,06 dans le liquide. Une telle droite est donnée par une loi appelée la
loi d’HENRY :

̅∗
𝒑
𝒚* = = m. 𝒙 (2.2)
𝒑𝒕
Avec :
m : est une constante
𝑝̅* : la pression partielle du soluté à l’équilibre
pt : pression totale du mélange gazeux

La non applicabilité de la loi de HENRY à des concentrations élevées peut être le résultat
d'interactions chimiques avec le liquide ou de dissociation électrolytique, comme c'est le cas
avec l'ammoniac-eau, ou la non-idéalité en phase gazeuse.

On peut s’attendre à ce que la plupart des gaz respectent la loi d’HENRY pour atteindre des
pressions d’équilibre d’environ 5 x 105 N/m2 (environ 5 atm), bien que, si la solubilité est faible,
comme dans le cas hydrogène-eau, elle peut être respectée à des pressions pouvant atteindre
34 N/m2. Les gaz du type vapeur (qui sont en dessous de leur température critique) suivront
généralement la loi jusqu'à des pressions d'environ 50% de la valeur de saturation à la
température de travail à condition qu'aucune action chimique ne se produise dans le liquide.
Dans tous les cas, m doit être déterminée expérimentalement.

II-5. Choix de solvant pour l’absorption

Si le but principal de l'opération d'absorption est de produire une solution spécifique,
comme dans la fabrication de l'acide chlorhydrique par exemple, le solvant est spécifié par la
nature du produit final. Si le but principal est d'éliminer certains composants du gaz, un choix
est souvent possible. L’eau est bien sûr le solvant le moins cher et le plus abondant, mais les
propriétés suivantes sont des considérations importantes :

7
1. Solubilité du gaz : La solubilité du gaz doit être élevée, pour augmenter le taux
d'absorption et diminuer la quantité de solvant requise. Généralement, les solvants de
nature chimique similaire à celle du soluté à absorber procurent une bonne solubilité.
Ainsi, les huiles d’hydrocarbures, et non pas l'eau, sont utilisées pour éliminer le
benzène du gaz sous-produit d’une chaudière.
Si une réaction chimique a lieu entre le soluté et le solvant, le taux d'absorption est
généralement très élevé. Mais si le solvant doit être récupéré pour être réutilisé, la
réaction doit être réversible. Par exemple, le sulfure d'hydrogène peut être éliminé des
mélanges de gaz à l'aide de solutions d'éthanolamine, car le sulfure est facilement
absorbé à basse température et facilement éliminé à haute température. La soude
caustique absorbe le sulfure d'hydrogène de manière excellente mais ne le libère pas
lors d'une opération de stripping.

2. Volatilité : Le solvant doit avoir une faible pression de vapeur car le gaz sortant une
opération d'absorption est habituellement saturé de solvant et une grande partie peut
en être perdue. Si nécessaire, un second solvant moins volatil peut être utilisé pour
récupérer la partie évaporée du premier, figure 2.2.

Gaz purifié
Récupération

Solvant non volatile

du solvant

Solvant récupéré + Absorbant

Solvant
Absorption

Entrée du Gaz

Solution

Figure 2.2 : Absorption avec section de récupération des solvants volatils

Par exemple, dans le cas d'absorbeurs d'hydrocarbures, où une huile solvant

relativement volatile est utilisée dans la partie principale de l'absorbeur, en raison des
caractéristiques de solubilité supérieures et où le solvant volatilisé est récupéré du gaz
par une huile non volatile. De même, le sulfure d’hydrogène peut être absorbé par une
solution aqueuse de phénolate de sodium, mais le gaz désulfuré est lavé avec de l’eau
pour récupérer le phénol évaporé.

8
 3. Corrosion : Le solvant non corrosif ou moins corrosif réduit les coûts de construction
l'équipement ainsi que les coûts d'entretien
4. Coût : Le solvant devrait être bon marché, pas cher, de sorte que les pertes seront
insignifiantes et devrait être facilement disponible.
5. Viscosité : Une viscosité assez basse est préférable. Elle donne de grande vitesse
d’absorption, améliore les caractéristiques d’engorgement des colonnes d’absorption,
résulte en des pertes de charge basse et de bonnes caractéristiques de transfert de
matière et de chaleur
6. Autres : Si possible, le solvant doit être non toxique, non-inflammable, chimiquement
stable et avoir un point de congélation bas

III- Transfert de matière d’un constituant ; Bilan de matière

Dans ce chapitre nous allons voir les bilans de transfert de matière appliqués aux opérations
d’absorption et de stripping. L’élaboration de ces bilans va nous permettre de trouver les
équations des droites opératoires qui régissent ces opérations. Ensuite, la résolution de ces
équations de bilans va nous permettre de dimensionner la colonne (colonne à garnissage) à
utiliser pour une opération d’absorption.

III-1. Bilan de matière pour un écoulement à contre-courant

La figure 3.1 montre le schéma d’une tour à contre-courant qui peut être soit une tour à
garnissage, soit une tour remplie de tout autres matériaux pour établir un contact liquide-gaz.
Le flux gazeux à n’importe quel point P de la colonne est défini par un débit molaire G
(mol/surface*temps), constitué d’un soluté qui diffuse A d’une fraction y, d’une pression
partielle p et d’un rapport molaire Y, et d’un gaz inerte non transféré de débit molaire Gs
(mol/surface*temps).

9
 moles/(temps)(surface)

moles/(temps)(surface)

Figure 3.1 : Représentation des flux d’une colonne d’absorption ou de stripping

Flux Gazeux Flux liquide

Symbole Désignation Symbole Désignation
Débit molaire gazeux à un point P de la Débit molaire liquide à un point P de
G L
colonne (mole/temps*surface) la colonne (mole/temps*surface)
Débit molaire gazeux à l’entrée de la Débit molaire liquide à la sortie de la
G1 L1
colonne (mole/temps*surface) colonne (mole/temps*surface)
Débit molaire gazeux à la sortie de la Débit molaire liquide à l’entrée de la
G2 L2
colonne (mole/temps*surface) colonne (mole/temps*surface)
Débit molaire gazeux du solvant sans Débit molaire liquide du solvant sans
Gs Ls
soluté (mole/temps*surface) soluté (mole/temps*surface)
Fraction molaire du soluté A dans le gaz à Fraction molaire du soluté A dans le
y1 l’entrée de la colonne (mole de A/mole x1 liquide à la sortie de la colonne (mole
du gaz total) de A/mole du liquide total)
Fraction molaire du soluté A dans le gaz à Fraction molaire du soluté A dans le
y2 la sortie de la colonne (mole de A/mole x2 liquide à l’entrée de la colonne (mole
du gaz total) de A/mole du liquide total)
Fraction molaire du soluté A dans le gaz Fraction molaire du soluté A dans le
y au point P dans la colonne (mole de x liquide au point P dans la colonne
A/mole du gaz total) (mole de A/mole du liquide total)
Rapport molaire du soluté A dans le gaz à Rapport molaire du soluté A dans le
Y1 l’entrée de la colonne (mole de A/mole X1 liquide à la sortie de la colonne (mole
du solvant gazeux) de A/mole du solvant liquide)

10
 Rapport molaire du soluté A dans le gaz à Rapport molaire du soluté A dans le
Y2 la sortie de la colonne (mole de A/mole X2 liquide à l’entrée de la colonne (mole
du solvant gazeux) de A/mole du solvant liquide)
Rapport molaire du soluté A dans le gaz Rapport molaire du soluté A dans le
Y au point P dans la colonne (mole de X liquide au point P dans la colonne
A/mole du solvant gazeux) (mole de A/mole du solvant liquide)
Pression partielle du soluté A à l’entrée
p1
de la colonne
Pression partielle du soluté A à la sortie
p2
de la colonne
pt Pression totale du mélange gazeux

Les relations qui relient les différents paramètres sont les suivantes :
𝑦 𝑝
𝑌= = (3.1)
1 − 𝑦 𝑝𝑡 − 𝑝

𝐺
𝐺𝑠 = 𝐺(1 − 𝑦) = (3.2)
1+𝑌

𝑥
𝑋=
1−𝑥 (3.3)

𝐿
𝐿𝑠 = 𝐿(1 − 𝑥) = (3.4)
1+𝑋

Remarque :
Les débits des solvants gazeux et liquide Gs et Ls reste essentiellement constant à travers
toute la colonne. Alors, il est commode d’exprimer les bilans de matières en ces termes.

Un bilan de matière sur le soluté A dans la partie inférieure de la colonne (enveloppe I figure
3.1) nous donne :

𝐺𝑠 (𝑌1 − 𝑌) = 𝐿𝑠 (𝑋1 − 𝑋) (3.5)

𝐿𝑠 𝐿𝑠
𝑌= 𝑋 + (𝑌1 − 𝑋1 ) (3.6)
𝐺𝑠 𝐺𝑠

11
 L’équation (3.5) et (3.6) représente l’équation d’une droite de coordonnées (X, Y) et de
𝐿𝑠
pente et qui passe par le point de coordonnées (X1, Y1). Cette droite est appelée la droite
𝐺𝑠
opératoire (figure 3.2).
La substitution de X par X2 et Y par Y2 montre que cette droite passe par le point (X2, Y2)
pour une colonne d’absorption.
Cette ligne droite indique la relation entre la concentration dans le gaz Y et la
concentration dans le liquide X à n’importe quel point de la colonne.

(bas)
Courbe
Droite opératoire, d’équilibre
Moles soluté /moles solvant dans le

pente Ls/Gs
Courbe
d’équilibre (Haut)

Droite
opératoire,
pente Ls/Gs
gaz

(bas)

Moles soluté/moles solvant dans le liquide

(b) Stripping
(a) Absorption
Figure 3.2 : Droite Opératoire pour une Absorption et un Stripping

Les données de solubilité à l'équilibre pour le soluté gaz dans le solvant liquide peuvent
également être tracées dans le même diagramme, avec les mêmes unités, comme la courbe
MN de la figure 3.2, par exemple. Chaque point de cette courbe représente la concentration
de gaz en équilibre avec le liquide à sa concentration correspondante et à une température
définie.

Remarque :
Pour un absorbeur (transfert de matière du gaz au liquide), la droite opératoire se trouve
toujours au-dessus de la courbe d'équilibre, tandis que pour le cas du stripping (transfert de
matière du liquide vers le gaz), la droite opératoire est toujours en dessous, comme illustré à
la figure 3.2.

12
 La courbe opératoire est une ligne droite uniquement lorsqu'elle est tracée en fonction des
rapports molaires. En fonction des fractions molaires ou des pressions partielles cette ligne
est courbée (non droite) comme sur la figure 3.3 pour un absorbeur. L’équation de cette
courbe est donnée par :
𝑦1 𝑦 𝑝1 𝑝 𝑥1 𝑥
𝐺𝑠 ( − ) = 𝐺𝑠 ( − ) = 𝐿𝑠 ( − ) (3.7)
1 − 𝑦1 1 − 𝑦 𝑝𝑡 − 𝑝1 𝑝𝑡 − 𝑝 1 − 𝑥1 1 − 𝑥

Courbe d’équilibre

y* = f(x)
Fraction molaire du soluté dans le gaz

Courbe opératoire

Fraction molaire du soluté dans le liquide

Figure 3.3 : Courbe Opératoire en fraction molaire

Remarque :
La pression total pt est généralement considérée constante à n’importe quel point de la colonne.

III-2. Le rapport minimal liquide-gaz pour un absorbeur

Lors du dimensionnement des colonnes d’absorption, la quantité du gaz à traiter G ou Gs,
les concentrations extrêmes Y1 et Y2 et la composition du liquide entrant X2 sont généralement
connues et fixées par les exigences du procédé. Par contre la quantité du liquide à utiliser doit
être déterminer.

13
Moles soluté/moles solvant dans le gaz

Pente =
Ls /Gs

Pente =
Ls(min) /Gs

Moles soluté/moles solvant dans le liquide

Figure 3.4 : Rapport liquide-gaz minimal pour une colonne d’absorption

Pour la figure 3.4.a :

• La droite opératoire doit passer par le point D et doit se terminer à l’ordonnée Y1

• Si la quantité du liquide utilisée est telle sorte qu’on obtient la ligne DE, la
concentration du liquide à la sortie serait X1
• Si moins de liquide est utilisé, la composition du liquide de sortie sera nettement plus
grande, comme au point F, mais étant donné que les forces motrices de diffusion
seront plus petit, l’absorption sera plus difficile. Le temps de contact entre le gaz et le
liquide doit alors être plus long, et l'absorbeur correspondant doit être plus grand
• Le débit minimal du liquide pouvant être utilisé correspond à la droite opératoire DM,
qui présente la plus grande pente pour toute droite touchant la courbe d’équilibre et
elle est tangente à la courbe en P. En ce point, la force motrice par diffusion est nulle,
le temps de contact requis pour le changement de concentration souhaité est infini et
une colonne infiniment grande en résulte. Ceci représente alors le rapport liquide / gaz
minimal Ls(min) /Gs.
Pour la figure 3.4.b :

• La courbe d'équilibre est le plus souvent concave vers le haut, et le rapport minimal
liquide-gaz correspond alors à une concentration de liquide de sortie en équilibre avec
le gaz entrant.

Ces principes mentionnés ci-dessus s’appliquent également au stripping, où la droite opératoire

qui touche la courbe d’équilibre à un point quelconque représente le rapport maximal liquide- gaz
et une concentration maximale dans le gaz de sortie.

14
 III-3. Bilan de matière pour un écoulement en co-courant

moles/(temps)(surface
moles/(temps)(surface
)
)

Droite opératoire,
pente -Ls/Gs

P
L G
Ls Gs Courbe
X Y d’équilibre

Figure 3.5 : Absorbeur à co-courant

Lorsque le gaz et le liquide s'écoulent à co-courant, comme dans la Fig. 3.5, la droite
opératoire présente une pente négative - Ls / Gs (voir démonstration ci-dessous).
Il n'y a pas de limite au rapport -Ls/Gs, mais une tour infiniment longue produirait une sortie
liquide et gaz en équilibre, comme pour le point (Xe, Ye) de la courbe figure 3.5.
Un bilan de matière sur le soluté A dans la partie supérieure de la colonne (enveloppe J
figure 3.5) nous donne :

𝐿𝑠 (𝑋1 − 𝑋) = 𝐺𝑠 (𝑌 − 𝑌1 ) (3.8)
𝐿𝑠 𝐿𝑠
𝑌=− 𝑋 + (𝑌1 + 𝑋1 ) (3.9)
𝐺𝑠 𝐺𝑠
C’est l’équation de la droite opératoire de pente -Ls/Gs et qui passe par les deux points
(X1,Y1) et (X2,Y2).
Comme déjà vu, si on travaille avec les fractions molaires au lieu des rapports molaires, la
ligne opératoire sera une courbe.

15
 L’écoulement à co-courant peut être utilisé lorsqu'une colonne suffisamment longue est
construite en deux sections, comme dans la Fig.3.6. La deuxième section opère généralement en
écoulement à co-courant pour économiser sur la conduite de gaz de grand diamètre reliant les deux
colonnes.
Il peut également être utilisé si le gaz à dissoudre dans le liquide est une substance pure, si le
fonctionnement à contre-courant ne présente aucun avantage, ou si une réaction chimique rapide
et irréversible avec le soluté dissous se produit dans le liquide.

Figure 3.6 : Arrangement à contre-courant et à co-courant pour les

très longues tours

IV- Colonnes à garnissage

IV.1 Généralités
Une colonne à garnissage est une colonne qui est remplie d’éléments permettant
d’augmenter la surface de contact entre la phase gazeuse et la phase liquide. Ce sont ces
éléments qui constituent le garnissage (figure 4.1). La colonne contient également :
▪ Dans sa partie inférieure, un plateau pour supporter le garnissage et un injecteur de gaz
▪ Dans sa partie supérieure, un distributeur de liquide et parfois un dévésiculeur
▪ Parfois un redistributeur
En considérant la façon d’alimenter la colonne en gaz et en liquide, on peut distinguer trois
types de colonnes : à contre-courant, à co-courant et à courant croisé
Dans les colonnes à garnissage, le liquide est pulvérisé sur le garnissage placé entre des
grilles. La phase liquide qui contient l’absorbant, forme un film sur les éléments de
garnissage (zone de mouillage). Le garnissage qui créé une série d’obstacles sur le
parcours du liquide et du gaz au sein de la colonne permet :

16
✓ D’augmenter la surface de contact entre les deux phases, et par là, la surface
d’échange par unité de volume de la colonne
✓ D’augmenter le temps de séjour des deux phases dans la colonne et donc leur
durée de contact
✓ De créer un régime turbulent de la phase gaz favorisant ainsi le brassage et
l’homogénéisation local des concentrations

Sortie du gaz

Entrée du liquide
Distributeur du liquide

Colonne

Garnissage en vrac

Redistributeur du liquide

Support du garnissage

Entrée du gaz

Sortie du liquide

Figure 4.1 : Colonne à garnissage à contre-courant

17
 Les colonnes garnies peuvent avoir une grande efficacité. L’efficacité pour une colonne
donnée, dépend d’un grand nombre de facteurs, à savoir :

• La vitesse du gaz
• La solubilité du soluté dans la solution
• La concentration de l’absorbant dans la solution
• La réactivité éventuelle du gaz avec les constituants du liquide
• La diffusivité du soluté dans la phase gazeuse et de l’absorbant dans la phase
liquide
• La teneur en soluté de la phase gazeuse
• La hauteur de la colonne

Les colonnes à garnissage présentent un grand nombre d’avantage dont on peut citer :
✓ Elles sont moins chères que les colonnes à plateaux, surtout si le diamètre de la
colonne n’est pas trop grand.
✓ Elles sont bien adaptées pour fonctionner dans un environnement acide et corrosif
✓ Comme le liquide est peu agité, il est possible d’utiliser les absorbants qui ont tendance
à mousser
✓ Elles peuvent atteindre une grande efficacité d’absorption pour beaucoup de gaz
✓ La faible rétention du liquide peut être aussi un avantage, principalement si
l’absorbant est sensible à la température

IV.2 Le garnissage
Pour répondre aux exigences de l’opération d’absorption, le garnissage doit avoir certaines
propriétés dont on cite :
o Il doit être chimiquement inerte vis-à-vis des fluides utilisés
o Il doit être assez résistant, mais pas trop lourd
o Il doit contenir des passages adéquats pour les fluides en circulation sans répartition
de liquide ou de pertes excessives
o Il doit permettre un bon contact entre le gaz et le liquide
o Il ne doit pas être trop cher

Les garnissages sont généralement caractérisés par deux propriétés intrinsèques :

- La surface du garnissage par unité de volume garni ap (m2/m3)

- Fraction du vide dans le garnissage ou porosité ε0

18
Ainsi, un bon garnissage doit présenter :
֍ Une surface par unité de volume suffisamment grande pour offrir une surface de
contact importante entres les deux fluides
֍ Un vide élevé pour permettre le passage à des débits de fluide également importants
sans créer de pertes de charge élevées
Il existe deux types principaux de garnissage, garnissage rempli en vrac « Random » et
garnissage régulier ou systématique « Regular packing ».
a) Garnissage en vrac
Les garnissages en vrac sont simplement déversés dans la colonne lors de l’installation et
laissés tomber d’une façon aléatoire. Dans le passé, des matériaux facilement disponibles tels
que les pierres brisées, du gravier ou des morceaux de coke étaient utilisés comme garnissage.
Mais même s’ils sont très peu coûteux, ils ne sont pas bien appréciés pour des raisons telles
que leur surface réduite et les mauvaises caractéristiques d'écoulement de fluide qu’il
présentent. Les garnissages en vrac les plus utilisés de nos jours sont fabriqués. Les types les
plus courants sont illustrés à la figure 4.2.

Figure 4.2 : différents types du garnissage en vrac

19
Anneaux Raschig (figure 4.2a)
Ce sont des manchons cylindriques de hauteur égale au diamètre extérieur (dimension
nominale). Le choix du matériau est extrêmement vaste : céramique (grès ou porcelaine),
graphite, métaux et alliages, plastiques et verres. Les dimensions nominales courantes vont
de 6,4 à 102 mm (1/4 à 4 in).
Anneaux Pall (figure 4.2b)
La surface intérieure d’un anneau Raschig se révèle beaucoup moins accessible au
mouillage que sa surface extérieure, d’où l’idée de ménager des fentes transversales dans les
parois en repoussant la matière à l’intérieure. Les anneaux Pall ainsi constitués permettent de
réduire le volume de garnissage à mettre en œuvre, pour une opération donnée, d’environ 30
à 40 %. Ces anneaux, de plus en plus répandus, existent en céramique, métaux et plastiques ;
les dimensions nominales courantes vont de 15,9 à 76 mm (5/8 à 3 in).
Selles de Berl (figure 4.2c)
Ces corps ont une forme dérivée du paraboloïde hyperbolique avec renforcement des
arêtes et adjonction de nervures. Cette forme assure au lit une structure enchevêtrée
provoquant un autoblocage qui diminue la poussée exercée par la charge sur les parois de la
colonne. On constante également une amélioration de l’homogénéité de la fraction de vide.
Le volume à mettre en œuvre pour une opération déterminée est sensiblement de 20 %
inférieur à ce qu’il serait avec des anneaux Raschig. La production, essentiellement en
céramique ou en plastique, s’effectue par moulage. Les dimensions nominales courantes vont
de 6,4 à 76 mm (1/4 à 3 in).
Selles Intalox (figure 1d)
La forme générale de ces corps est celle d’un demi-tore engendré par un demi-cercle dont
la concavité est tournée vers l’extérieur. Leur production, essentiellement en céramique ou
en plastique, peut être obtenue par extrusion. L’emboîtement des selles Intalox est moins
marqué que celui des selles de Berl; c’est à la fois un avantage (meilleure porosité de la charge,
diminution des points de contact générateurs d’angles morts) et un inconvénient (moins
bonne stabilité d’où poussée latérale plus importante). En première approximation, le volume
à mettre en œuvre pour une opération déterminée est le même qu’avec des selles de Berl.

b) Garnissages systématiques
Anneaux ordonnés
Les empilages d’anneaux ordonnés sont souvent utilisés comme premières couches
interposées entre une grille formée de poutres et le garnissage en vrac (figure 4.3). Cette
disposition est particulièrement recommandée pour les colonnes de grand diamètre lorsque
la charge de garnissage est importante ou lorsque les problèmes de corrosion requièrent une
grille briquetée.

20
 Figure 4.3 : Interposition d’anneaux ordonnés entre une grille

Matelas
Ils sont obtenus par juxtaposition de bandes de tissus métalliques ou plastiques gaufrées.
Les ondes de gaufrage forment des canaux inclinés à 45o parfois simples, plus souvent en
chevrons. Les éléments, de 0,1 à 0,15 m (4 à 6 in) de hauteur, sont cylindriques, d’une seule
pièce pour un diamètre inférieur ou égal à 0,3 m (1 ft), en plusieurs pièces pour des
diamètres supérieurs.
Garnissages structurés
Ce sont des éléments autoporteurs formés de canaux inclinés à 30 ou 45o sur la verticale,
le diamètre hydraulique de ces canaux est variable depuis environ 0,01 m suivant les
produits ; ces canaux peuvent être éventuellement munis de perforations. Ces éléments
sont fournis en acier, en plastique ou en céramique.
La figure 4.4 montre d’autres exemples de garnissage systématique.

a) Anneaux Rashing ordonnés

b) Anneaux double spirales

c) Garnissage en métal déployé

d) Grilles de bois

Figure 4.4 : Garnissage systématique

21
 IV.3 Support du garnissage
Un espace libre au bas de la tour est nécessaire pour assurer une bonne distribution du gaz
dans le garnissage. Par conséquence, le garnissage doit être supportée au-dessus de l’espace
libre. Le support doit, bien entendu, être suffisamment solide pour supporter le poids d’un
garnissage d’une hauteur raisonnable, et il doit avoir suffisamment de surface libre pour
permettre le passage du liquide et du gaz. Une grille de barres, du type montré à la Figure 4.1,
peut être utilisée. Mais des supports spécialement conçus qui fournissent des passages
séparés pour le gaz et le liquide sont préférés. La figure 4.5 montre un exemple dont la zone
libre pour l'écoulement est de l'ordre de 85% et qui peut être réalisée avec diverses
modifications et avec de nombreux matériaux, notamment des métaux, des céramiques et
des plastiques.

Gaz Liquide

Figure 4.5 : Exemple de support de garnissage

IV.4 Distributeur du liquide

L'importance d'une distribution initiale adéquate du liquide au sommet du garnissage est
indiquée à la figure 4.6. Divers dispositifs sont utilisés pour la distribution de liquide. Les buses
de pulvérisation généralement conduits à un grand entrainement du liquide dans le gaz pour
être utiles. La disposition illustrée à la Figure 4.1 peut être utilisé dans de petites colonnes.
Pour les grands diamètres, un distributeur comme celui illustré à la Fig. 4.7 peut être utilisé,
et de nombreuses autres dispositions sont disponibles. Il est généralement jugé nécessaire de
prévoir au moins cinq points d’introduction de liquide pour chaque 0,1 m² de section
transversale de la colonne pour les grandes tours (d > 1,2 m) et un nombre supérieur pour les
plus petits diamètres.

22
 Entrée Liquide Entrée Liquide

Garnissage
sec

Garnissage
mouillé

Figure 4.6 : distribution du liquide et mouillage du garnissage Figure 4.7 : distributeur du liquide

En particulier à des vitesses de gaz élevées, le gaz sortant du sommet du garnissage peut
entraîner des gouttelettes de liquide. Les gouttes sont généralement arrêtées par un matelas
dévésiculeur tissé en métal ou en plastique de 0,10 à 0,15 m d’épaisseur solidement maintenu
sur une grille. Les filets gazeux contournent les fils tandis que les gouttes, d’inertie plus élevée,
les heurtent et s’y accrochent. Les gouttes s’agglomèrent sur les fils et s’en détachent dès que
la pesanteur compense les forces capillaires.
Le dévésiculeur est généralement placé à une hauteur égale à un diamètre de colonne, au-
dessus du lit de garnissage.
La détermination de la section d’un dévésiculeur découle conventionnellement du choix
d’une vitesse superficielle de la phase gazeuse G’/ρG fondée sur la formule de Souders et
Brown :
(4.1)
𝐺′
= 𝐾 √(ρ𝐿 − ρ𝐺 )/ρ𝐺
ρ𝐺

Avec : G’(kg/m².s) débit massique du gaz

ρG, ρL masse volumique respectivement du gaz et du liquide
K (m/s) constante de Souders et Brown
Les constructeurs préconisent K=0.08 m/s, le fonctionnement étant toutefois satisfaisant
dans la fourchette 0,032 à 0,107 m/s. Le calcul du diamètre d’un dévésiculeur fondé sur la
valeur K=0.08 m/s conduit à un diamètre de matelas significativement inférieur (de 70 à 80 %
en règle générale), à celui de la colonne.

23
 V- Hydrodynamiques des colonnes à garnissage

V.1 Perte de charge pour un écoulement monophasé

La perte de charge subie par un seul fluide lors de son écoulement à travers un lit de
solides compacts est bien raisonnablement représentée par l’équation d’ERGUN :

∆𝑝 𝑔𝑐 𝜀 3 𝑑𝑝 𝜌𝐺 150(1 − 𝜀) (5.1)
= + 1.75
𝑍 (1 − 𝜀)𝐺′2 𝑅𝑒
Cette équation s'applique aussi bien aux flux gazeux qu’aux flux liquides.
Le terme à gauche de l’équation est un facteur de friction. Ceux de droite représentent les
contributions au facteur de friction.
𝟏𝟓𝟎(𝟏−𝜺)
Le premier terme [ ] est pour un écoulement purement laminaire.
𝑹𝒆

Le second terme [1.75] est pour un écoulement complètement turbulent.

L’importance relative de chacun de ces deux termes change à mesure que le débit change.
𝑑𝑝 𝐺′
𝑅𝑒 = : le nombre de Reynolds.
𝜇

𝑑𝑝 : le diamètre effectif des grains solides (garnissage). Il correspond au diamètre d’une

sphère ayant le même rapport surface/volume que le solide en question.

Si ap est la surface spécifique du garnissage (m²/m3 lit fixe). La surface par unité du volume
𝒂𝒑
des grains est .
𝟏−𝜺
Alors :
6(1 − 𝜀)
𝑑𝑝 = (5.2)
𝑎𝑝

Normalement, Ce ne sera pas la même chose que la taille nominale des particules.
Pour des débits de gaz G’ supérieur un peut prêt à 0.7 kg/m².s (500lb/ft².h), le premier terme
à droite de l’équation (5.1), celui du régime laminaire, peut être négligé. Pour un type de
garnissage et une taille de particule donnée, l’équation (5.1) peut être simplifiée par
l’expression empirique suivante :
∆𝑝 𝐺′
= 𝐶𝐷 (5.3)
𝑍 𝜌𝐺
Les valeurs du coefficient CD sont listées par le tableau 5.1.

24
 Paramètre Désignation SI Unités anglo-saxonnes
∆𝑝 Perte de charge linaire 𝑁/𝑚² 𝑙𝑏/𝑓𝑡²
𝑍 𝑚 𝑓𝑡
Diamètre effectif des grains 𝑚 𝑓𝑡
𝑑𝑝
solides (garnissage)
𝜀 Porosité du lit - -
𝜌𝐺 Masse volumique du gaz 𝑘𝑔/𝑚3 𝑙𝑏/𝑓𝑡 3
𝜌𝐿 Masse volumique du liquide 𝑘𝑔/𝑚3 𝑙𝑏/𝑓𝑡 3
𝐺′ Débit massique du gaz 𝑘𝑔/𝑚2 . 𝑠 𝑙𝑏/𝑓𝑡 2 . ℎ
𝑔𝑐 Facteur de correction des unités 𝑔𝑐 = 1 𝑔𝑐 = 4.18 ∗ 108
𝑅𝑒 Nombre de Reynolds - -
𝜇 Viscosité dynamique 𝑘𝑔/𝑚. 𝑠 𝑐𝑃

𝑎𝑝 Surface spécifique du garnissage 𝑚2 /𝑚3 𝑓𝑡 2 /𝑓𝑡 3

𝐽 Facteur de correction 𝐽=1 𝐽 = 1.502

𝐶𝑓
Facteur caractéristique du Valeur 𝐶𝑓 à partir du tableau 5.1
garnissage
𝑑 Diamètre nominal du garnissage 𝑚 𝑓𝑡

V.2 Perte de charge pour un écoulement biphasé (deux phases)

Pour un écoulement simultané de liquide et de gaz à contre-courant, les données de perte

de charge de divers chercheurs font apparaître des écarts importants dus à des différences de
densité et de fabrication du garnissage. Les estimations ne peuvent donc pas être très
précises.
Dans la plupart des cas, la corrélation généralisée de la figure 5.1 sera utile. Les valeurs de Cf
figurent dans le tableau 5.1.

V.3 Ecoulement à contre-courant du liquide et du gaz à travers les garnissages

Pour la plupart des garnissages en vrac, la perte de charge subie par le gaz est influencée à
la fois par les débits de gaz et de liquide de manière similaire à celle illustrée à la Figure5.2.

La pente de la droite représentant la perte de charge dans le garnissage sec est généralement
comprise entre 1.8 et 2, indiquant un écoulement turbulent pour la plupart des vitesses de gaz.

25
Figure 5.2 : Perte de charge linéique et rétention liquide opératoire dans un garnissage en
fonction de la vitesse massique du gaz 26
 ✓ A une vitesse de gaz fixe, la perte de charge augmente avec l’augmentation du débit
liquide car, la section de passage libre du gaz est réduite par la présence du liquide.
✓ Dans la région au-dessous de la courbe II, la rétention liquide h0, la quantité d’eau
retenue dans le garnissage, reste constante avec le changement de la vitesse du gaz et
elle augmente avec l’augmentation de la vitesse du liquide
✓ Dans la région entre II et I la rétention liquide augmente rapidement avec la vitesse
du gaz. La section libre du passage du gaz devient plus petite et la perte de charge
augmente rapidement. C’est la zone de charge
✓ Au fur et à mesure que le débit du gaz augmente pour atteindre la courbe I, à un débit
de liquide fixe, le liquide est entrainé par le flux gazeux et le phénomène
d’engorgement se produit. Ainsi la perte de charge augmente d’une façon très rapide
✓ Il n’est pas pratique qu’une colonne à garnissage opère dans les conditions
d’engorgement.

La zone de charge est moins bien connue que la limite d’engorgement et les colonnes à
garnissage industrielles sont généralement dimensionnées pour fonctionner entre 60 et 70 % du
débit gazeux à l’engorgement (ce qui pratiquement correspond à la zone de charge).

V.4 Choix du garnissage et calcul du diamètre

a) Taux de mouillage
Le débit-volume de liquide doit être suffisant pour assurer un film qui ruisselle sur le
maximum de la surface du garnissage ; si le débit de liquide est insuffisant, une partie de la
surface disponible ne sera pas utilisée, soit parce qu’elle reste sèche, soit encore parce que le
film ne se renouvelle pas.

• Morris et Jackson appellent taux de mouillage (wetting rate) le quotient du débit-

volume de liquide par le périmètre maximal mouillable dans une section droite (A*ap)
et prescrivent :
𝐿′ (5.4)
≥ 0.000025 𝑚2 /𝑠
𝑎𝑝 𝜌𝐿

27
 • Par une approche un peu différente, Mersmann arrive à :
𝑎𝑝 𝑑 2
∗ 𝑅𝑒𝐿 (𝐺𝑎𝐿 )−3 ≥ 0.0002 (5.5)
6𝜀

𝑅𝑒𝐿 est le nombre de Reynolds et 𝐺𝑎𝐿 le nombre de Galillée :

𝑑 𝐿′
𝑅𝑒𝐿 = (5.6
𝜇𝐿 (5.5)

𝑔 𝑑 3 𝜌𝐿2 )
𝐺𝑎𝐿 = (5.7)
𝜇𝐿 ²
g : accélération de la pesanteur : g = 9,81 m/s2

Ces réquisitions établies pour des anneaux Raschig en acier ou en céramique ne semblent pas avoir
été actualisées pour les garnissages actuels en plastique ; en l’absence d’autres informations, il est
prudent de reconduire le minimum de Morris et Jackson (0,000025 m2/s).
Il ne faut pas non plus tomber dans l’excès contraire. Il ne faut pas avoir un taux de mouillage
supérieur à la valeur estimé par Bugarel qui est de 0,000 2 m2/s.

b) Surface mouillée
Même lorsque la distribution du liquide au-dessus du lit de garnissage est assurée de façon
uniforme et le taux de mouillage compris entre 0,000025 et 0,0002 m2/s, la surface volumique
mouillée aW reste inférieure à la surface volumique géométrique ap (figure 5.3). Avant la zone
de charge, la surface mouillée est :
- Sensiblement indépendante de la vitesse du gaz
- Une fonction croissante du débit de liquide
Onda et coll proposent :
𝑎𝑤 𝛾𝑐 0.75 −0.35
= 1 − exp (−1.45 ( ) (𝑎𝑝 𝑑) (𝐺𝑎𝐿 )0.05 (𝑊𝑒𝐿 )0.2 ) (5.8)
𝑎𝑝 𝛾

𝐺𝑎𝐿 est le nombre de Galilée et 𝑊𝑒𝐿 le nombre de Weber :

𝑑 𝐿′² (5.9)
𝑊𝑒𝐿 =
𝜌𝐿2 𝛾
Avec : 𝛾 (kg.s-2 ou N.m-1 ) tension superficielle du liquide.
𝛾𝑐 (kg.s-2 ou N.m-1 ) valeur caractéristique de la tension superficielle dépendant du
matériau de garnissage (voir tableau ci-dessous)

28
 Figure 5.3 : Aire mouillée et aire effective en fonction de la vitesse superficielle du liquide

La surface volumique effective ae est la surface effectivement utilisée pour le transfert de

matière et peut être décomposée en deux termes :
o Une fraction seulement de la surface mouillée (en effet, il y a lieu d’exclure les
zones mortes, localisations principales de la rétention statique)
o L’aire des gouttes qui tombent d’un élément sur le suivant
Tant que la zone de charge n’est pas atteinte, la surface effective est sensiblement
indépendante de la vitesse du gaz. Aux valeurs modérées de la vitesse superficielle du liquide
(L’/ρL ≈ 1 m/s), la contribution de la surface des gouttes est faible et la surface effectivement
utilisée est inférieure à la surface mouillée ; à mesure que la vitesse superficielle du liquide
croît, la contribution des gouttes compense l’influence des zones mortes et la surface effective
devient égale puis supérieure à la surface mouillée (figure 5.3).

29
 Tableau 5.1 : Caractéristiques des garnissages en vrac

Taille nominale d, mm (inch)

Garnissage 6 9.5 13 16 19 25 32 38 50 76 89
(1/4) (3/8) (1/2) (5/8) (3/4) (1) (11/4) (11/2) (2) (3) (31/2)
Anneaux Rashig
Céramique :
Epaisseur, mm 0.8 1.6 2.4 2.4 2.4 3 4.8 4.8 6 9.5
Cf 1600 1000 580 380 255 155 125 95 65 37
CD 909 749 457 301 181.8 135.6
𝜺 0.73 0.68 0.63 0.68 0.73 0.73 0.74 0.71 0.74 0.78
ap, m²/m3 (ft²/ft3) 787(240) 508(155) 364(111) 328(100) 262(80) 190(58) 148(45) 125(38) 92(28) 62(19)

Métal :
0.8 mm epaisseur
Cf 700 390 300 170 155 115
𝜺 0.69 0.84 0.88 0.92
ap, m²/m3 (ft²/ft3) 774(236) 420(128) 274(83.5) 206(62.7)

1.6 mm epaisseur
Cf 410 290 220 137 110 83 57 32
CD 688 431 485 304 172.9 133.5
𝜺 0.73 0.78 0.85 0.87 0.9 0.92 0.95
ap, m²/m3 (ft²/ft3) 387(118) 236(71.8) 186(56.7) 162(49.3) 135(41.2) 103(31.4) 68(20.6)

Anneaux Pall
Plastique :
Cf 97 52 40 25 16
CD 207 105.2 61.8 47.5 23.9
𝜺 0.87 0.9 0.91 0.92 0.92
ap, m²/m3 (ft²/ft3) 341(104) 206(63) 128(39) 102(31) 85(26)

Métal :
Cf 70 48 28 20 16
CD 133.4 95.5 56.6 36.5
𝜺 0.93 0.94 0.95 0.96
ap, m²/m3 (ft²/ft3) 341(104) 206(63) 128(39) 102(31)

Flexirings :
Cf 78 45 28 22 18
𝜺 0.92 0.94 0.96 0.96 0.97
ap, m²/m3 (ft²/ft3) 345(105) 213(65) 131(40) 115(35) 92(28)

Hy-pak :
Cf 45 18 15
CD 88.1 28.7 26.6
𝜺 0.96 0.97 0.97

Selles de Berl
Céramique :
Cf 900 240 170 110 65 45
CD 508 295 184
𝜺 0.6 0.63 0.66 0.69 0.75 0.72
ap, m²/m3 (ft²/ft3) 899(274) 466(142) 269(82) 249(76) 144(44) 105(32)

Selles Intalox
Céramique :
Cf 725 330 200 145 98 52 40 22
CD 399 256 241.5 96.2 71.3 40.6
𝜺 0.75 0.78 0.77 0.775 0.81 0.79
ap, m²/m3 (ft²/ft3) 984(300) 623(190) 335(102) 256(78) 195(59.5) 118(36)

Palastique :
Cf 33 21 16
CD 96.7 56.5 30.1
𝜺 0.91 0.93 0.94
ap, m²/m3 (ft²/ft3) 207(63) 108(33) 89(27)

30
 )1/2
Engorgement
(flooding)

𝐺 𝜌𝐿 − 𝜌𝐺
𝜌𝐺
Perte de charge dans le gaz

𝑋𝑒 = ′ (
𝐿′
(N/m²)/m

𝜌𝐺 (𝜌𝐿 − 𝜌𝐺 ) 𝑔𝑐
𝑌𝑒 =
𝐺′2 𝐶𝑓 𝜇𝐿0.1 𝐽

Figure 5.1 : Engorgement et perte de charge dans les colonnes à garnissage en vrac 31
 c) Calcul du diamètre de la colonne
Le mécanisme du calcul du diamètre de la colonne est présenté par les étapes suivantes :

Calcul
𝐿′ 𝜌𝐺
𝑋𝑒 = ′
( )1/2
𝐺 𝜌𝐿 − 𝜌𝐺

A partir de la figure 5.1 :

𝐺𝑒′2 𝐶𝑓 𝜇𝐿0.1 𝐽
𝑌𝑒 =
𝜌𝐺 (𝜌𝐿 − 𝜌𝐺 ) 𝑔𝑐

Calcul du flux gazeux à

l’engorgement 𝐺𝑒′

Calcul du flux gazeux

𝐺 ′ = 𝛼𝐺𝑒′
𝛼 = 0.6 à 0.7

Section droite de la colonne : Perte de charge

𝐷é𝑏𝑖𝑡 𝑚𝑎𝑠𝑠𝑖𝑞𝑢𝑒 𝑑𝑢 𝑔𝑎𝑧 𝑘𝑔/𝑠
par unité de
𝑆= ∆𝑝
𝐺′ longueur
𝑍

Calcul du diamètre : Flux massique

Diamètre corrigé Section corrigée
corrigé
4𝑆 Dcorrigé Scorrigé G’corrigé
𝐷=ඨ
𝜋

Remarque :
𝐷/𝑑 ≥ 12 entre le diamètre de la colonne et le diamètre nominal du garnissage doit être satisfaite
pour éviter la présence d’important vide dans la périphérique de la colonne et ainsi limiter les
écoulement préférentiels du liquide le long des parois.

32
 VI- Transfert de matière dans les colonnes garnies

Lorsqu'une tour à garnissage est utilisée de la manière habituelle pour une absorption ou
un stripping à contre-courant, pour le transfert d’un soluté entre le gaz et le liquide, le taux
de transfert du soluté peut être calculé à partir des valeurs mesurées du débit du gaz et du
liquide et des concentrations de soluté dans le flux entrant et sortant.
En raison des difficultés de mesure des concentrations en soluté à l’interface gaz-liquide et
des difficultés de mesure de la surface d’échange entre le gaz et le liquide, le flux de transfert
de matière ne peut pas être mesuré directement. On se contente donc de mesurer la vitesse
de transfert de matière, comme étant le produit du flux de transfert de matière et la surface
interfacial. Cette vitesse de transfert de matière, une fois divisée par le volume du garnissage,
donne le coefficient volumique globaux de transfert de matière, Kx a, Ky a, KG a, KL a, FOG a,FOL
a, etc. avec a est la surface interfaciale par unité de volume de garnissage (m2/m3).
Les coefficients partiels de transfert de matière (kx , ky, FG , FL) et la surface interfaciale a
qui forment ces coefficients volumétriques globaux dépendent différemment des propriétés
du fluide, des débits et du type de garnissage. Les coefficients volumétriques globaux ne sont
donc utiles que pour le dimensionnement de colonnes remplis du même garnissage et opèrent
avec les mêmes produits chimiques aux mêmes débits et concentrations que ceux utilisés lors
des mesures de ces coefficients.
Pour obtenir des coefficients partiels, l’approche générale est de choisir des conditions
expérimentales telles que la résistance au transfert de matière en phase gazeuse soit
négligeable par rapport à celle du liquide. C'est le cas pour l'absorption de gaz très insolubles,
par exemple l'oxygène ou l'hydrogène dans l'eau. Les mesures dans de tels systèmes
conduisent alors à des valeurs de kx a, kL a et FL a, qui peuvent être corrélés en termes de
variables du système. Il n’existe évidemment aucun système d’absorption dans lequel le
soluté est si soluble dans le liquide que la résistance en phase liquide est tout à fait
négligeable. Mais en soustrayant la résistance au liquide connue des résistances globales, il
est possible d’obtenir les coefficients en phase gazeuse ky a, kG a et FG a et ensuite de les
corréler en termes de variables du système.

33
 VI.1 Rétention liquide

La rétention liquide désigne le liquide retenu dans la colonne sous forme de films mouillant
le garnissage et de flaques emprisonnées dans les crevasses entre les particules du garnissage.
On constate que la rétention liquide totale Φ𝐿𝑡 est composé de deux parties :
Φ𝐿𝑡 = Φ𝐿𝑜 + Φ𝐿𝑠 (6.1)
Avec :
Φ𝐿𝑡 : Rétention totale
Φ𝐿𝑜 : Rétention opératoire. Constitue le liquide en mouvement continu sur le garnissage. Lors
de l’arrêt du flux de gaz et de liquide, ce liquide s’écoule du garnissage. Cette rétention est
fonction du taux de mouillage
Φ𝐿𝑠 : Rétention statique. Elle est surtout localisée aux points de contact entre les éléments du
garnissage et dans les pores du solide. Cette quantité du liquide est quasi-stagnante, ne quitte
pas la colonne lors de l’arrêt et indépendante du taux de mouillage.

Remarque :
Le liquide de la retenue statique s'équilibre rapidement avec le gaz adjacent, sa surface interfaciale ne
contribue donc pas au transfert de matière.

VI.2 Coefficients de transfert de matière

Pour les anneaux de Raschig et les selles de Berl, Shulman et ses associés ont déterminé les
caractéristiques des coefficients de transfert de matière sans surface « kG » en faisant passer
des gaz à travers des lits remplis de garnissage en naphtalène qui se subliment dans le gaz. En
comparant ceux-ci avec les « kG a » à partir de l'absorption liquide, les surfaces interfaciales
d'absorption et de vaporisation ont été obtenues. Les données sur les « kL a » ont ensuite
fourni la corrélation pour le coefficient de transfert de la phase liquide kL.
Pour les anneaux de Raschig et les selles de Berl le coefficient de transfert de matière coté
gazeux est donné par la corrélation suivante :

𝐹𝐺 𝑆𝑐𝐺 2/3 𝑘𝐺 𝑝𝐵,𝑀 𝑆𝑐𝐺 2/3 𝑑𝑠 𝐺 ′ (6.2)

= = 1.195 [ ]−0.36
𝐺 𝐺 𝜇𝐺 (1 − 𝜀𝐿𝑜 )
Avec 𝜀𝐿𝑜 est la porosité opératoire du garnissage :
𝜀𝐿𝑜 = 𝜀 − Φ𝐿𝑡 (6.3)

34
𝜀 : la porosité du lit de garnissage (tableau 5.1)
Φ𝐿𝑡 : rétention liquide total (Equation 6.1)
𝑑𝑠 : diamètre d’une sphère ayant la même surface qu’un élément de garnissage (m)
(tableau 6.1) (noter ds ≠ dp )
𝜇𝐺 : viscosité dynamique du gaz ( 𝑘𝑔/𝑚. 𝑠)
G : débit molaire spécifique du gaz (mol/m².s)
G’ : débit massique spécifique du gaz (kg/m².s)
𝑘𝐺 : coefficient de transfert de matière gazeux (mol/s.m².Pa) (utilisé pour les solutions
diluées)
𝑝𝐵,𝑀 : moyenne logarithmique des pressions partielles du soluté gazeux aux deux extrémités
de la colonne (Pa)
𝐹𝐺 : coefficient de transfert de matière gazeux (mol/s.m²)
𝜇𝐺
𝑆𝑐𝐺 : nombre de Schmidt gazeux 𝑆𝑐𝐺 = 𝜌
𝐺 𝐷𝐺

𝜌𝐺 : masse volumique du gaz (kg/m3)

𝐷𝐺 : diffusivité gazeuse (m²/s)

Les propriétés du fluide doivent être évaluées dans les conditions moyennes du travail entre
l’interface et le gaz.
Le coefficient de transfert de matière coté liquide est donné par :
𝑘𝐿 𝑑𝑠 𝑑𝑠 𝐿′ 0.45
= 25.1 ( ) 𝑆𝑐𝐿 1/2 (6.4)
𝐷𝐿 𝜇𝐿
Vu que les données sur les liquides ont été obtenues à de très faibles concentrations de
soluté, il est possible de convertir 𝑘𝐿 en 𝐹𝐿 par 𝑭𝑳 = 𝒌𝑳 ∗ 𝒄 avec 𝒄 est la masse volumique
molaire du solvant liquide.
Avec :
𝐿′ : débit massique spécifique du liquide (kg/m².s)
𝜇𝐿 : viscosité dynamique du liquide ( 𝑘𝑔/𝑚. 𝑠)
𝑘𝐿 : coefficient de transfert de matière liquide ( 𝑚/𝑠)
𝐷𝐿 : diffusivité liquide (m²/s)
𝜇𝐿
𝑆𝑐𝐿 : nombre de Schmidt gazeux 𝑆𝑐𝐿 = 𝜌
𝐿 𝐷𝐿

35
 • Calcul de la porosité opératoire
Pour calculer 𝜀𝐿𝑜 qui figure dans les corrélations de calcul des coefficients de transfert de
matière la rétention liquide totale Φ𝐿𝑡 doit être calculer :
D’après l’équation 6.3 :

𝜀𝐿𝑜 = 𝜀 − Φ𝐿𝑡
𝜀 est la porosité du lit de garnissage donnée par le tableau 5.1
On a aussi d’après l’équation 6.1 :
Φ𝐿𝑡 = Φ𝐿𝑜 + Φ𝐿𝑠

Φ𝐿𝑠 rétention liquide statique est donnée directement par le tableau 6.1.
Φ𝐿𝑜 rétention liquide opératoire est calculé comme suite :
(6.5)
Φ𝐿𝑜 = Φ𝐿𝑜𝑊 ∗ 𝐻
(6.6)
Φ𝐿𝑜𝑊 = Φ𝐿𝑡𝑊 − Φ𝐿𝑠𝑊
Φ𝐿𝑜𝑊 : rétention opératoire en eau (w=water)
Φ𝐿𝑠𝑊 : rétention statique en eau
Φ𝐿𝑡𝑊 : rétention totale en eau
𝐻 : facteur de correction
Les expressions de Φ𝐿𝑜𝑊 , Φ𝐿𝑠𝑊 et 𝐻 sont données par le tableau 5.1.
Dans le tableau 6.1:
𝜇𝐿 : viscosité du liquide kg/m.s
𝜎 : tension superficielle du liquide N/m
𝐿′ : flux massique du liquide kg/m².s
𝜌𝐿 : la masse volumique du liquide kg/m3

36
 Tableau 6.1 : la rétention liquide dans le garnissage

Taille
nominale 𝝁𝑳 ,
Garnissage ds, m 𝚽𝑳𝒔 Eau, (températures ordinaires)
kg/m.s
𝑯
mm in

𝛽 = 1.508𝑑𝑠0.376 0.57 0.13

975.7𝐿′ 𝜇𝐿 𝜎 0.1737−0.262 log(𝐿′ )
2.47 ∗ 10 −4 < 0.012 0.430 ( )
Anneaux 13 0.5 0.01774 Φ𝐿𝑠𝑊 = 𝜌𝐿0.84 (2.024𝐿′ − 1) 0.073
0.0486𝜇𝐿0.02 𝜎 0.99 𝑑𝑠1.21
Raschig au 25 1 0.0356
céramique 38 1.5 0.0530 𝑑𝑠1.21 𝜌𝐿0.37 (2.09 ∗ 10−6 )(737.5𝐿′ )𝛽 0.57 0.31
Φ𝐿𝑡𝑊 = 2168𝐿′ 𝜇𝐿 𝜎 0.1737−0.262 log(𝐿′ )
50 2 0.0725 𝑑𝑠2 > 0.012 ()
𝜌𝐿0.84 (2.024𝐿 ′ 0.430− 1) 0.073

𝛽 = 1.104𝑑𝑠0.376 0.57 0.13

407.9𝐿′ 𝜇𝐿 𝜎 0.1737−0.262 log(𝐿′ )
−4 < 0.012 ()
Anneaux 5.94 ∗ 10 𝜌𝐿0.84 (1.393𝐿 ′ 0.315
− 1) 0.073
25 1 0.01301 0.0237𝜇𝐿0.02 𝜎 0.23 Φ𝐿𝑠𝑊 =
Raschig au 𝑑𝑠1.21
carbon
38 1.5 0.0543 𝑑𝑠1.21 𝜌𝐿0.37
50 2 0.0716 (7.34 ∗ 10−6 )(737.5𝐿′ )𝛽 0.57 0.31
Φ𝐿𝑡𝑊 = 901𝐿′ 𝜇𝐿 𝜎 0.1737−0.262 log(𝐿′ )
𝑑𝑠2 > 0.012 ()
𝜌𝐿0.84 (1.393𝐿 ′ 0.315
− 1) 0.073

𝛽 = 1.508𝑑𝑠0.376
0.57 0.13
1404𝐿′ 𝜇𝐿 𝜎 0.2817−0.262 log(𝐿′ )
< 0.020 0.413 ( )
Selles de 13 0.5 0.31622 4.23 ∗ 10−3 𝜇 0.04 𝜎 0.55 5.014 ∗ 10−5 𝜌𝐿0.84 (3.24𝐿′ − 1) 0.073
Berl au 25 1 0.0320 𝐿 Φ𝐿𝑠𝑊 =
1.56 0.37
𝑑𝑠 𝜌𝐿 𝑑𝑠1.56
céramique 38 1.5 0.0472
(2.32 ∗ 10−6 )(737.5𝐿′ )𝛽 0.57 0.31
Φ𝐿𝑡𝑊 = 2830𝐿′ 𝜇𝐿 𝜎 0.28177−0.262 log(𝐿′ )
𝑑𝑠2 > 0.020 ( )
𝜌𝐿0.84 (3.24𝐿 ′ 0.413− 1) 0.073

37
 • Calcul de la surface interfaciale

Pour le cas d’une absorption ou d’un stripping avec l’eau ou des solutions aqueuses très
diluées, les surfaces interfaciales 𝒂𝑨𝑾 (m2/m3) sont regroupées sous forme de graphe ou
donnée par des expressions empiriques. Tableau 6.2
808 𝐺 ′ 𝑛 𝑝
𝑎𝐴𝑊 = 𝑚( ) 𝐿′ (6.7)
𝜌𝐺 0.5
𝐿′ : flux massique du liquide (kg/m².s)
𝐺 ′ : flux massique du gaz (kg/m².s)
𝜌𝐺 : massa volumique du gaz (kg/m3)
m, n et p sont par le tableau 6.2

Pour le cas d’une absorption ou d’un stripping avec des solutions non aqueuses on a :
Φ𝐿𝑜
𝑎𝐴 = 𝑎𝐴𝑊 (6.8)
Φ𝐿𝑜𝑊
Pour le cas du contact entre un gaz et un liquide pur, l’évaporation par exemple, on a :
Φ𝐿𝑡
𝑎𝑉 = 0.85𝑎𝐴 (6.9)
Φ𝐿𝑜
Φ𝐿𝑡𝑊
𝑎𝑉𝑊 = 0.85𝑎𝐴𝑊 (6.10)
Φ𝐿𝑜𝑊

𝑎𝑝 : surface spécifique du garnissage (m2/m3 de garnissage)

𝑎𝐴 : surface interfaciale spécifique d’absorption, stripping ou distillation (m2/m3 de

garnissage)
𝑎𝑉 : surface interfaciale spécifique de contact entre un gaz et liquide pur (m2/m3 de
garnissage)
𝑎 : surface interfaciale spécifique de transfert de matière (m2/m3 de garnissage)

38
 Tableau 6.2 : Surface interfaciale d’absorption et stripping pour les solutions aqueuses

Taille nominale Intervalle de L’

Garnissage
mm in Kg/m².s Lb/ft².h m n p
0.68-2.0 500-1500 28.01 0.2323L’- 0.3 -1.04
13 0.5
2.0-6.1 1500-4500 14.69 0.01114L’+ 0.148 -0.111
0.68-2.0 500-1500 34.42 0 0.552
25 1
Anneaux 2.0-6.1 1500-4500 68.2 0.0389L’- 0.0793 -0.47
Raschig 0.68-2.0 500-1500 36.5 0.0498L’- 0.1013 0.274
38 1.5
2.0-6.1 1500-4500 40.11 0.01091L’- 0.022 0.14
0.68-2.0 500-1500 31.52 0 0.481
50 2
2.0-6.1 1500-4500 34.03 0 0.362

0.68-2.0 500-1500 16.28 0.0529 0.761

13* 0.5*
2.0-6.1 1500-4500 25.61 0.0529 0.170
0.68-2.0 500-1500 52.14 0.0506L’- 0.1029 0
Selles de Berl 25 1
2.0-6.1 1500-4500 73 0.0310L’- 0.0630 -0.359
0.68-2.0 500-1500 40.6 -0.0508 0.455
38 1.5
2.0-6.1 1500-4500 62.4 0.0240L’- 0.0996 -0.1355

* Seulement pour des débit de gaz G’ < 1.08 kg/m².s (800 lb/ft².h)

39
 VI.3 Cinétique de transfert de matière et calcul de la hauteur d’une colonne à
garnissage

Dans les colonnes à garnissage, les fluides sont en contact continu toute au long de leur
passage dans la colonne. Ainsi, la composition du liquide et du gaz change en continu avec la
hauteur du garnissage. Chaque point de la courbe opératoire représente donc des conditions
trouvées à un niveau donné dans la colonne.
a- Coefficients et unité de transfert partiels
Considérons une colonne à garnissage de section droite d'unité (A= 1 m²) Figure 6.1. La
surface interfaciale effective totale pour le transfert de matière résultante de la dispersion du
liquide sur le garnissage par unité de section droite de la colonne est S.
S est définit comme le produit de la surface interfaciale spécifique (a) et la hauteur du
garnissage, ou volume spécifique (volume/unité de section droite), (Z). Pour un volume
différentiel dZ, la surface interfaciale est :
dS= a dZ (6.11)

La quantité de soluté A dans le gaz passant par la section différentielle considérée de la

colonne est G.y (mol/(surface)(temps)) et le taux de transfert de matière est alors d(G.y)
(molA/(volume différentiel)(temps)).

Figure 6.1 : Colonne à garnissage

40
 • Cas de solutions diluées et le transfert d’un seul constituant
Le flux de transfert de matière est :
𝑁𝐴 = 𝑘𝑦 (𝑦 − 𝑦𝑖 ) (6.12)
𝑁𝐴 : flux de transfert de matière mol/(surface)(temps)
𝑘𝑦 : coefficient partiel de transfert de matière coté gaz mol/(surface)(temps)

𝑦 : fraction molaire de A dans le gaz

𝑦𝑖 : fraction molaire de A à l’interface gaz-liquide

On a aussi la quantité de A transféré dans dZ est d(G.y) :

𝐺𝑠 . 𝑦 𝑦 𝑑𝑦 𝑑𝑦
𝑑(𝐺. 𝑦) = 𝑑 ( ) = 𝐺𝑠 𝑑 ( ) = 𝐺𝑠 =𝐺 (6.13)
1−𝑦 1−𝑦 (1 − 𝑦)² (1 − 𝑦)
Un bilan de matière sur dZ nous donne :
𝑑(𝐺. 𝑦) 𝑑(𝐺. 𝑦)
𝑁𝐴 = = (6.14)
𝑑𝑆 𝑎 𝑑𝑍
D’après (6.12), (6.13) et (6.14) on a :
𝑑𝑦
𝐺 = 𝑘𝑦 (𝑦 − 𝑦𝑖 )
(1 − 𝑦) 𝑎 𝑑𝑍
𝐺 𝑑𝑦
𝑑𝑍 =
𝑘𝑦 𝑎(1 − 𝑦)(𝑦 − 𝑦𝑖 ) (6.15)

La hauteur du garnissage est donnée donc par :

𝑍 𝑦1
𝐺 𝑑𝑦
𝑍 = ∫ 𝑑𝑍 = ∫ (6.16)
0 𝑦2 𝑘𝑦 𝑎(1 − 𝑦)(𝑦 − 𝑦𝑖 )

De même si on travaille coté liquide on a :

𝑍 𝑥1
𝐿 𝑑𝑥
𝑍 = ∫ 𝑑𝑍 = ∫ (6.17)
0 𝑥2 𝑘𝑥 𝑎(1 − 𝑥)(𝑥𝑖 − 𝑥)

Les valeurs de 𝑦𝑖 et 𝑥𝑖 sont obtenues par :

𝑁𝐴 = 𝑘𝑦 (𝑦 − 𝑦𝑖 ) = 𝑘𝑥 (𝑥𝑖 − 𝑥)
(𝑦 − 𝑦𝑖 ) −𝑘𝑥
=
(𝑥 − 𝑥𝑖 ) 𝑘𝑦 (6.18)

L’équation (6.18) indique le point (𝑥𝑖 , 𝑦𝑖 ) appartient à l’intersection entre la courbe

−𝑘𝑥
d’équilibre et la droite de pente passant par le point (𝑥, 𝑦). Figure 6.2
𝑘𝑦

41
 Courbe d’équilibre yi=f(xi)

y Pente

yi

x xi
Figure 6.2 : concentration à l’interface à partir de la courbe d’équilibre

• Cas de solutions non diluées ou transfert de plus qu’un seul constituant

Lorsque on traite des situations qui n'impliquent ni la diffusion d'une seule substance, ni
la contre-diffusion équimolaire ou si les taux de transfert de matière sont élevés, les
coefficients de type F doivent être utilisés.
L’approche générale reste la même, bien que les expressions résultantes soient plus
complexes que celles développées ci-dessus.
La forme générale du flux dans ce cas est :
𝑁𝐴 (𝑁𝐴 / ∑ 𝑁) − 𝑦𝑖 𝑁𝐴 (𝑁𝐴 / ∑ 𝑁) − 𝑥
𝑁𝐴 = 𝐹𝐺 𝑙𝑛 = 𝐹𝐿 𝑙𝑛 (6.19)
∑𝑁 (𝑁𝐴 / ∑ 𝑁) − 𝑦 ∑ 𝑁 (𝑁𝐴 / ∑ 𝑁) − 𝑥𝑖

Si on considère un mélange de deux constituant, soluté A et inerte B (𝑁𝐵 = 0) on a :

𝑁
𝑁𝐴 (𝑁 +𝐴 𝑁 ) − 𝑦𝑖 1 − 𝑦𝑖
𝐴 𝐵
𝑁𝐴 = 𝐹𝐺 𝑙𝑛 = 𝐹𝐺 𝑙𝑛 (6.20)
𝑁𝐴 + 𝑁𝐵 𝑁 1−𝑦
(𝑁 +𝐴 𝑁 ) − 𝑦
𝐴 𝐵

D’après (6.13), (6.14) et (6.20) on a :

𝑍 𝑦1
𝐺 𝑑𝑦 (6.21)
𝑍 = ∫ 𝑑𝑍 = ∫
1−𝑦
0 𝑦2 𝐹𝐺 𝑎(1 − 𝑦) 𝑙𝑛 1 − 𝑦𝑖

42
 𝑍 𝑥1
𝐿 𝑑𝑥 (6.22)
𝑍 = ∫ 𝑑𝑍 = ∫
1−𝑥
0 𝑥2 𝐹𝐿 𝑎(1 − 𝑥) 𝑙𝑛 1 − 𝑥
𝑖

Les valeurs de 𝑦𝑖 et 𝑥𝑖 sont obtenues par :

1 − 𝑦𝑖 1−𝑥
𝑁𝐴 = 𝐹𝐺 𝑙𝑛 = 𝐹𝐿 𝑙𝑛
1−𝑦 1 − 𝑥𝑖
1 − 𝑦𝑖 1 − 𝑥 𝐹 /𝐹 1 − 𝑥 𝐹 𝑎/𝐹 𝑎
=( ) 𝐿 𝐺 =( ) 𝐿 𝐺 (6.23)
1−𝑦 1 − 𝑥𝑖 1 − 𝑥𝑖

Pour chaque valeur de (𝑥, 𝑦) de la courbe opératoire, une courbe représentant 𝑦𝑖 en fonction de
𝑥𝑖 peut-être tracée. L’intersection de cette courbe avec la courbe d’équilibre fournit le 𝑦𝑖 qui va être
utilisé avec 𝑦 dans la résolution de l’équation (6.21). L’équation (6.21) peut être donc intégrée
graphiquement en représentant la fonction à intégrer en fonction de 𝑦.

On a aussi :
𝑦 − 𝑦𝑖 = (1 − 𝑦𝑖 ) − (1 − 𝑦)
Donc on peut écrire l’équation (6.21) comme suit :
𝑍 𝑦1
𝐺(1 − 𝑦)𝑖𝑀 𝑑𝑦 (6.24)
𝑍 = ∫ 𝑑𝑍 = ∫
0 𝑦2 𝐹𝐺 𝑎(1 − 𝑦) (𝑦 − 𝑦𝑖 )

Avec (1 − 𝑦)𝑖𝑀 est la moyenne logarithmique de (1 − 𝑦𝑖 ) et (1 − 𝑦).

(1 − 𝑦𝑖 ) − (1 − 𝑦)
(1 − 𝑦)𝑖𝑀 = (6.25)
1−𝑦
𝑙𝑛 1 − 𝑦𝑖

D’après les équations (6.16) et (6.24) on note que :

𝐹𝐺 = 𝑘𝑦 (1 − 𝑦)𝑖𝑀 = 𝑘𝐺 𝑃𝑡 (1 − 𝑦)𝑖𝑀

Pour les solutions diluées :

𝐹𝐺 ≈ 𝑘𝑦 ≈ 𝑘𝐺 𝑃𝑡

On définit la hauteur d’unité de transfert gazeux 𝐻𝑡𝐺 :

𝐺 𝐺 𝐺
𝐻𝑡𝐺 = = = (6.26)
𝐹𝐺 𝑎 𝑘𝑦 𝑎 (1 − 𝑦)𝑖𝑀 𝑘𝐺 𝑎 𝑃𝑡 (1 − 𝑦)𝑖𝑀

On définit la hauteur d’unité de transfert gazeux 𝑁𝑡𝐺 :

𝑦1
(1 − 𝑦)𝑖𝑀 𝑑𝑦
𝑁𝑡𝐺 = ∫ (6.27)
𝑦2 (1 − 𝑦) (𝑦 − 𝑦𝑖 )

43
 L’équation (6.24) devient donc :
𝑍 𝑦1
𝐺(1 − 𝑦)𝑖𝑀 𝑑𝑦 (6.28)
𝑍 = ∫ 𝑑𝑍 = ∫ = 𝐻𝑡𝐺 ∗ 𝑁𝑡𝐺
0 𝑦2 𝐹𝐺 𝑎(1 − 𝑦) (𝑦 − 𝑦𝑖 )

Remarque :
𝐺
Dans l’équation (6.28) on a pris l’avantage que le rapport 𝐹 est presque constant, ce qui n’est pas
𝐺𝑎
le cas pour 𝐺 et 𝐹𝐺 𝑎, pour le faire sortir de l’intégrale.

𝑁𝑡𝐺 : représente le nombre de fois qu’on a la moyenne (𝑦 − 𝑦𝑖 ) /(𝑦1 − 𝑦2 ). Il mesure la difficulté

d’absorption
𝐻𝑡𝐺 : est la hauteur du garnissage offrant une unité de transfert gazeux

L’équation (6.24) peut être simplifiée d’avantage par la substitution de la moyenne

logarithmique (1 − 𝑦)𝑖𝑀 par une moyenne arithmétique. Ceci peut introduire une petite
erreur :
(1 − 𝑦𝑖 ) − (1 − 𝑦) (1 − 𝑦𝑖 ) + (1 − 𝑦)
(1 − 𝑦)𝑖𝑀 = ≈ (6.29)
1 − 𝑦𝑖 2
𝑙𝑛 1 − 𝑦

𝑁𝑡𝐺 devient donc :

𝑦1
𝑑𝑦 1 1 − 𝑦2
𝑁𝑡𝐺 = ∫ + ln ( ) (6.30)
𝑦2 (𝑦 − 𝑦𝑖 ) 2 1 − 𝑦1

L’équation (6.30) est plus simple à intégrer graphiquement

Aussi :
log (𝑦1)
𝑑𝑙𝑜𝑔(𝑦) 1 1 − 𝑦2
𝑁𝑡𝐺 = 2.303 ∫ + log ( ) (6.31)
log (𝑦2 ) (𝑦 − 𝑦𝑖 ) 2 1 − 𝑦1

44
 Les relations ci-dessus ont toutes leurs contreparties en termes de concentrations liquides,
obtenues exactement de la même manière :
𝑍 𝑥1 𝑥1
𝐿 𝑑𝑥 𝐿(1 − 𝑥)𝑖𝑀 𝑑𝑥
𝑍 = ∫ 𝑑𝑍 = ∫ =∫ (6.32)
1−𝑥 𝐹𝐿 𝑎(1 − 𝑥) (𝑥𝑖 − 𝑥)
0 𝑥2 𝐹𝐿 𝑎(1 − 𝑥) 𝑙𝑛 1 − 𝑥 𝑥2
𝑖

𝑍 = 𝐻𝑡𝐿 ∗ 𝑁𝑡𝐿 (6.33)

𝐿 𝐿 𝐿 (6.34)
𝐻𝑡𝐿 = = =
𝐹𝐿 𝑎 𝑘𝑥 𝑎 (1 − 𝑥)𝑖𝑀 𝑘𝐿 𝑎 𝑐 (1 − 𝑥)𝑖𝑀
𝑥1 𝑥1 (6.35)
(1 − 𝑥)𝑖𝑀 𝑑𝑥 𝑑𝑥 1 1 − 𝑥1
𝑁𝑡𝐿 = ∫ ≈∫ + ln ( )
𝑥2 (1 − 𝑥) (𝑥𝑖 − 𝑥) 𝑥2 (𝑥𝑖 − 𝑥) 2 1 − 𝑥2

Remarque :
Les mêmes valeurs de Z sont obtenues en raisonnant soit coté gaz soit coté liquide.

Pour le cas du stripping on obtient les mêmes formules que celles établies ci-dessus. La force motrice
(𝑦 − 𝑦𝑖 ) et (𝑥𝑖 − 𝑥) qui apparait dans ces formules est négative pour le cas du stripping, mais vu que
𝑥2 > 𝑥1 et 𝑦2 > 𝑦1la hauteur Z calculé est positive.

b- Coefficients et unité de transfert globaux

Pour les cas où la courbe de distribution d'équilibre est une droite et le rapport des coefficients
de transfert de matière est constant, les coefficients globaux de transfert de matière sont
convenables. Les expressions pour la hauteur du garnissage peuvent alors être écrites :
𝑁𝐴 = 𝐾𝑦 (𝑦 − 𝑦 ∗ ) = 𝐾𝑥 (𝑥 ∗ − 𝑥)

𝑍 = 𝐻𝑡𝑂𝐺 ∗ 𝑁𝑡𝑂𝐺 (6.36)

𝑦1
(1 − 𝑦)∗𝑀 𝑑𝑦
𝑁𝑡𝑂𝐺 =∫ ∗
(6.37)
𝑦2 (1 − 𝑦) (𝑦 − 𝑦 )
𝑦1
𝑑𝑦 1 1 − 𝑦2
𝑁𝑡𝑂𝐺 = ∫ + ln ( ) (6.38)
𝑦2 (𝑦 − 𝑦 ∗ ) 2 1 − 𝑦1
𝑌1
𝑑𝑌 1 1 + 𝑌2
𝑁𝑡𝑂𝐺 = ∫ + ln ( ) (6.39)
∗
(𝑌 − 𝑌 ) 2 1 + 𝑌1
𝑌2
𝐺 𝐺 𝐺
𝐻𝑡𝑂𝐺 = = = (6.40)
𝐹𝑂𝐺 𝑎 𝐾𝑦 𝑎 (1 − 𝑦)∗𝑀 𝐾𝐺 𝑎 𝑃𝑡 (1 − 𝑦)∗𝑀
(1 − 𝑦 ∗ ) − (1 − 𝑦)
(1 − 𝑦)∗𝑀 = (6.41)
1 − 𝑦∗
𝑙𝑛 1 − 𝑦

45
Avec 𝑦 ∗ (ou 𝑌 ∗ ) est la concentration du soluté dans le gaz s’il était en équilibre avec le
liquide de concentration 𝑥 (ou 𝑋), de sorte que (𝑦 − 𝑦 ∗ ) ou (𝑌 − 𝑌 ∗ ) est simplement la
distance verticale entre la droite opératoire et la courbe d’équilibre.
𝑁𝑡𝑂𝐺 est nombre global d’unité de transfert gazeux
𝐻𝑡𝑂𝐺 est la hauteur globale d’unité de transfert gazeux
𝐹𝑂𝐺 , 𝐾𝑦 et 𝐾𝐺 sont les coefficients globaux de transfert de matière coté gaz

Des équations similaires sont obtenues en cas du raisonnement coté liquide :

𝑥1
(1 − 𝑥)∗𝑀 𝑑𝑥
𝑁𝑡𝑂𝐿 = ∫ (6.42)
𝑥2 (1 − 𝑥) (𝑥 ∗ − 𝑥)
𝑥1
𝑑𝑥 1 1 − 𝑥1
𝑁𝑡𝑂𝐿 = ∫ + ln ( ) (6.43)
𝑥2 (𝑥 ∗
− 𝑥) 2 1 − 𝑥2
𝑋1
𝑑𝑋 1 1 + 𝑋1
𝑁𝑡𝑂𝐿 = ∫ + ln ( ) (6.44)
∗
𝑋2 (𝑋 − 𝑋) 2 1 + 𝑋2
𝐿 𝐿 𝐿
𝐻𝑡𝑂𝐿 = = = (6.45)
𝐹𝑂𝐿 𝑎 𝐾𝑥 𝑎 (1 − 𝑥)∗𝑀 𝐾𝐿 𝑎 𝑐 (1 − 𝑥)∗𝑀

(1 − 𝑥 ∗ ) − (1 − 𝑥)
(1 − 𝑥)∗𝑀 = (6.46)
1 − 𝑥∗
𝑙𝑛 1 − 𝑥

c- Nombre global d’unité de transfert : solutions diluées

L’expression du nombre d’unité de transfert peut être simplifiée d’avantage pour le cas de
solution diluées. Dans le cas d’un mélange gazeux dilué le deuxième terme de la définition de
𝑁𝑡𝑂𝐺 (équation (6.38)) peut être négligé et l’équation devient :
𝑦1
𝑑𝑦
𝑁𝑡𝑂𝐺 = ∫ (6.47)
𝑦2 (𝑦 − 𝑦 ∗ )

Si la courbe d’équilibre est une ligne droite d’équation :

𝑦 ∗ = 𝑚𝑥 + 𝑟 (6.48)

Pour des solutions diluées la courbe opératoire également peut être comme une droite
d’équation :
𝐿
𝑦= (𝑥 − 𝑥2 ) + 𝑦2 (6.49)
𝐺
𝐿1 ≃ 𝐿2 = 𝐿 & 𝐺1 ≃ 𝐺2 = 𝐺
L’équation (6.48) et (6.49) nous donne donc :
𝑦 − 𝑦 ∗ = 𝑞𝑥 + 𝑠 (6.50)

46
q et s sont des constantes d’expression :
𝐿 𝐿
𝑞 =𝐺−𝑚 et 𝑠 = (1 − 𝐺 ) 𝑦2 + 𝑟

L’équation (6.47) devient donc :

𝐿 𝑥1 𝑑𝑥 𝐿 𝑦1 − 𝑦1∗ 𝑦1 − 𝑦2
𝑁𝑡𝑂𝐺 = ∫ = ln ( ∗ ) = (𝑦 − 𝑦 ∗ ) − (𝑦 − 𝑦 ∗ ) (6.51)
𝐺 𝑥2 𝑞𝑥 + 𝑠 𝐺𝑞 𝑦2 − 𝑦2 1 1 2 2
𝑦1 − 𝑦1∗
ln ( )
𝑦2 − 𝑦2∗

𝑦1 − 𝑦2
𝑁𝑡𝑂𝐺 = (6.52)
(𝑦 − 𝑦 ∗ )𝑀
(𝑦 − 𝑦 ∗ )𝑀 est la moyenne logarithmique de la différence de concentrations aux extrémités
de la colonne.
Autre forme de cette équation :
𝐺(𝑦1 − 𝑦2 ) = 𝐾𝐺 𝑎𝑍𝑃𝑡 (𝑦 − 𝑦 ∗ )𝑀 (6.53)

Solutions diluées : Loi d’Henry

Si la loi d’Henry est applicable (r de l’équation (6.48) = 0) :
𝑦 ∗ = 𝑚𝑥 (6.54)
L’équation (6.47) nous donne :
𝑦 − 𝑚𝑥 1 1
ln ( 𝑦1 − 𝑚𝑥2 (1 − 𝐴) + 𝐴 )
2 2 (6.55)
𝑁𝑡𝑂𝐺 =
1
1−𝐴
𝐿
Avec 𝐴 = 𝑚𝐺 est le facteur d’absorption.

Pour le cas du stripping on a l’expression suivante :

𝑥2 − 𝑦1 /𝑚
ln ( (1 − 𝐴) + 𝐴 )
𝑥1 − 𝑦1 /𝑚 (6.56)
𝑁𝑡𝑂𝐿 =
1−𝐴
La représentation graphique des expressions précédente est donné par la figure 6.3 :

47
Figure 6.3 : Nombre d’unité de transfert pour absorption et stripping en
48
fonction du facteur d’absorption
 d- Hauteur globale d’unité de transfert
Lorsqu’on utilise les nombres globaux de transfert de matière, les hauteurs globales de transfert de
matière doivent être utilisées :

Dans le cours de transfert de matière, on a démonté que dans le cas ou la courbe d’équilibre est une
droite de pente m (loi d’Henry) la relation entre les coefficients globaux et les coefficients partiels est
la suivante :
1 1 𝑚
= + (6.57)
𝐾𝑦 𝑘𝑦 𝑘𝑥
On trouve donc :

𝐺 𝐺(1 − 𝑦)𝑖𝑀 𝑚𝐺 𝐿 (1 − 𝑥)𝑖𝑀

= + (6.58)
𝐹𝑂𝐺 𝑎 𝐹𝐺 𝑎(1 − 𝑦)∗𝑀 𝐿 𝐹𝐿 𝑎 (1 − 𝑦)∗𝑀

(1 − 𝑦)𝑖𝑀 𝑚𝐺 (1 − 𝑥)𝑖𝑀 (6.59)

𝐻𝑡𝑂𝐺 = 𝐻𝑡𝐺 + 𝐻𝑡𝐿
(1 − 𝑦)∗𝑀 𝐿 (1 − 𝑦)∗𝑀
Si la résistance au transfert de matière est surtout coté gazeux, 𝑦𝑖 ≃ 𝑦 ∗ on a :
𝑚𝐺 (1 − 𝑥)𝑖𝑀
𝐻𝑡𝑂𝐺 = 𝐻𝑡𝐺 + 𝐻𝑡𝐿 (6.60)
𝐿 (1 − 𝑦)∗𝑀
En plus pour les solutions diluées on a :
𝑚𝐺
𝐻𝑡𝑂𝐺 = 𝐻𝑡𝐺 + 𝐻 (6.61)
𝐿 𝑡𝐿
D’une manière similaire on a du côté liquide :
(1 − 𝑥)𝑖𝑀 𝐿 (1 − 𝑦)𝑖𝑀 (6.62)
𝐻𝑡𝑂𝐿 = 𝐻𝑡𝐿 + 𝐻𝑡𝐺
(1 − 𝑥)∗𝑀 𝑚𝐺 (1 − 𝑥)∗𝑀
Si la résistance au transfert de matière est surtout coté liquide, 𝑥𝑖 ≃ 𝑥 ∗ on a :
𝐿 (1 − 𝑦)𝑖𝑀
𝐻𝑡𝑂𝐿 = 𝐻𝑡𝐿 + 𝐻𝑡𝐺 (6.63)
𝑚𝐺 (1 − 𝑥)∗𝑀
En plus pour les solutions diluées on a :
𝐿
𝐻𝑡𝑂𝐿 = 𝐻𝑡𝐿 + 𝐻 (6.64)
𝑚𝐺 𝑡𝐺
Les données de calcul des coefficients partiels sont présentées dans le chapitre 5.

49
 Références

- Robert [Link], Mass-Transfer Operations, 3rd edition, Ed. McGraw-Hill, New-York,

1980
- J-C. CICILE & P. COPIGNEAUX, Distillation. Absorption, Technique de l’ingénieur, 2002

50