Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Manuel Absorption FIVA 2015 PDF

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 12

ABSORPTION

GAZ-LIQUIDE

T.P. 2me Anne

I Labsorption du CO2
I.1 Rappel :
L'absorption est une opration unitaire prsente dans un grand nombre de procds,
notamment trs utilise pour la purification des gaz. Dans ce cas, l'objectif de
l'opration est d'effectuer le transfert d'un compos, gnralement faiblement
concentr, d'une phase gazeuse vers la phase liquide. Le compos absorb peut se
dissoudre physiquement dans le liquide et ventuellement y ragir chimiquement.
Labsorption est donc un phnomne de dissolution dun produit ltat gazeux
dans un liquide dans le solvant.

I.2 Cas tudi


L'tude propose ici concerne l'absorption du CO2 contenu dans un flux d'air par de
l'eau (absorption "physique"), ou dans une solution de soude (absorption "chimique").
Laugmentation du CO2 atmosphrique par la pollution urbaine et industrielle est un
rel problme cologique. Elle contribue en effet lamplification de leffet de serre
atmosphrique.
La capture du CO2 est donc actuellement un thme fort de lenvironnement.
Labsorption du CO2 par un solvant est une technologie envisage actuellement pour
capter le CO2 la sortie des centrales lectriques, grosses mettrices de CO2 ; cest le
captage en post combustion.

I.3 Dispositif exprimental


Cette opration est ralise dans un contacteur gaz-liquide de type "colonne
garnissage" o le gaz et le liquide circulent contre courant. La colonne de diamtre 80
mm est constitue de deux portions de garnissages de 1 m chacune (voir schma de
linstallation annexe 1). Le garnissage est de type vrac, compos danneaux Raschig
6mm, en verre, vitant ainsi tout problme de corrosion par rapport au solvant. Les
caractristiques du garnissage sont fournies en tableau 1.

___________________________________________________________________________
0,8 mm

6 mm

Type :

Anneaux Raschig

Matriaux :

verre

Aire spcifique :

895 m2/m3

Nombre dlment/m3 :

3 023 000

Fraction de vide :

80%.

__________________________________________________________________________
6 mm
Tableau 1 : Caractristiques du garnissage

Mesures :
Diffrentes mesures peuvent tre ralises sur ce dispositif exprimental :
- les tempratures dentre et de sortie du gaz et du liquide.
- la teneur en CO2 en phase gaz par un appareil de mesure bas sur l'absorption du
rayonnement infrarouge par le CO2 contenu dans l'air.
- perte de charge dans la colonne par mesure de hauteur deau.

I.4 Dmarrage et scurit


Dmarrage :
La colonne se dmarre prfrentiellement en mettant seulement le dbit liquide. Celui-ci
doit tre tout dabord tre important, pendant quelques minutes, afin de bien mouiller le
garnissage, pour ensuite tre ramen la valeur du dbit liquide choisi pour ltude.
Ensuite le dbit gaz est amen sa valeur, progressivement.
Arrt :
La colonne doit tre rince aprs son utilisation avec de la soude. Mettre seulement un
dbit deau, durant environ 15 minutes.
Scurit :
Lors de ltude avec le solvant soude, le port des lunettes et des gants est ncessaire
lors de la prise dchantillons de la phase liquide.

II Etudes
La colonne peut tre schmatise ainsi o le dbit de gaz G1 et le dbit de liquide Ln
circulent contre-courant :
Sortie gaz
Gn, yn

Entre liquide
Ln, xn

Entre gaz
G1, y1

Sortie liquide
L1, x1

Deux grands types de phnomnes interviennent dans cette opration et permettent de


dimensionner la colonne et de la caractriser en termes d'efficacit :
- l'hydrodynamique de la colonne,
- les phnomnes d'quilibre thermodynamique et de transfert du CO2 entre la
phase gaz et la phase liquide.
Dans le cas o l'on considre une raction dans la phase liquide, lefficacit de
labsorbeur est augmente.
3

Lobjectif du TP est donc dtudier en premire partie lhydrodynamique afin de


dterminer les dbits de fonctionnements, et en deuxime partie de quantifier la
quantit de CO2 absorbe avec et sans raction.

II.1. Etude de lhydrodynamique


II.1.2 Rappel thorique
L'hydrodynamique de la colonne est lie aux fluides et la colonne utilise. Les
principaux paramtres qui interviennent sont :
- la nature des fluides (tat, viscosit, densit),
- le dbit des fluides,
- les grandeurs gomtriques de la colonne (hauteur, diamtre),
- les caractristiques du garnissage (type, taille, densit, etc.).
L'coulement du liquide se fait sous forme d'un film sur la surface du garnissage.
L'paisseur de ce film, ainsi que la surface recouverte, augmentent avec le dbit. Pour
un dbit de liquide modr donn, avec un dbit de gaz, contre courant, il est possible
de distinguer deux zones de fonctionnement qui sont identifiables partir de la mesure
de la perte de charge du gaz la traverse du garnissage :
- faible valeur de dbit de gaz, la loi de variation de perte de charge du gaz au travers
du garnissage est voisine de celle obtenue pour un garnissage sec.
- des valeurs plus leves, il y a interaction entre les deux coulements et la rtention
de liquide augmente. La perte de charge du gaz augmente alors plus rapidement avec
le dbit de gaz. Dans ce type de fonctionnement, les turbulences cres au voisinage
de l'interface gaz-liquide sont les plus favorables au transfert de matire. Au del d'une
certaine valeur du dbit de gaz, la perte de charge crot extrmement rapidement et
provoque l'engorgement de la colonne.
De l'hydrodynamique va dpendre le temps de contact entre les deux phases, la
surface d'change "ae" et donc la quantit de matire change entre les deux phases.
On distingue ainsi trois surfaces :
- la surface gomtrique du garnissage (ou aire spcifique) : a
- la surface mouille : aw
- la surface dchange : ae
La surface dchange est proche de la surface mouille ; en revanche, la surface
mouille peut tre beaucoup plus faible que la surface gomtrique si le dbit liquide
est trop faible pour bien arros le garnissage.
Le taux de mouillage peut tre valu par la corrlation de Onda :

aw
= 1 exp 1,45 c
a

L
Avec Ga L =

g d3 L

L 2

WeL =

0.75

d L2
L L

(a d )

0 , 35

Ga L

0.05

0.2
We L

Et
L: viscosit dynamique (Pa.s)

L: masse volumique (kg/m3)


L : dbit spcifique liquide dans les conditions de fonctionnement (kg/m.s)
aw: surface spcifique mouille du garnissage (m/m3)
a : surface spcifique du garnissage (m/m3)
d : diamtre du garnissage (m)
pour raschig a.d = 4,7
L: tension de surface du liquide (0,072 N /m pour leau, idem pour la soude)
c: tension de surface critique du matriau du garnissage (verre 0,073 N /m)

D'autre part, pour que la colonne fonctionne normalement, il faut viter l'engorgement
de la colonne pour le gaz et pour le liquide.
Les points d'engorgement peuvent tre prdits par la corrlation de LOBO qui relie le
"flow factor" FL Ye avec

FL =

L G
G L

et

1 1
( l ) 0.2 eau Fp
g G L eau
L

masse volumique du gaz (kg/m3), pour l'air 20C, 1.29 kg/m3


masse volumique du liquide (kg/m3), pour l'eau, 1000 kg/m3
viscosits (Pa s)
Facteur de garnissage (m2/m3)

Avec G
L

FP
Avec Fp =

Ye = G e

( 0 )3

comme valeur thorique

0
fraction de vide du garnissage.
g
acclration de la pesanteur = 9.81 ms-2
Ge , L ou G densits de dbit massique
= (dbit massique/section droite de la colonne vide) en kg s-1 m-2
Le lissage de la courbe exprimentale de LOBO est donn par
ln Ye = -[1.023 ln FL + 0.1065(ln FL)2 + 3.89]

II.1.2 Etudes
L'tude de l'hydrodynamique est ralise avec de l'eau pour la phase liquide et de l'air
pour la phase gaz. Pour cette tude on utilisera de leau du rseau en alimentation
directe sur la colonne.
Faire varier le dbit de gaz G dbit de liquide L constant pour 2 dbits de liquide L et
relever la perte de charge dans la colonne.
Dans chaque cas :
5

- Tracer la courbe log P = f (log G) et observer visuellement lengorgement de


la colonne. Attention de ne pas recommencer une exprience avant que la
colonne ne soit dsengorge.
- Dterminer les points de charge et d'engorgement (Gc et Ge).
- Calculer Fp le facteur de garnissage en utilisant corrlation de LOBO.
- Calculer le taux de mouillage par chaque dbit liquide.
- Comparer et commenter vos rsultats.

II.2. Etude du transfert de matire en absorption physique


II.2.2 Rappel thorique
Le transfert de matire est souvent reprsent par la thorie du double film. Aux abords
de l'interface gaz-liquide, il y a un gradient de concentration de CO2 ct gaz et un
gradient de concentration ct liquide.
Direction du transfert
Cg
C

Film
gazeux

Phase
gazeuse

Film
liquide

Cgi

Phase
liquide

CLi

CL

Cg

CL

Interface gaz-liquide

Dans le cas de l'absorption de CO2 dans l'eau, on peut montrer que le "transfert est
limit par la phase liquide". On considre que l'on peut crire CLi = quilibre (Cgi)
quilibre (Cg)=C*. C* reprsente la solubilit du CO2 dans l'eau la temprature et la
pression partielle considre.
C* peut tre valu par la loi de Henry
PCO2= H CCO2*
Avec PCO2 pression partielle en CO2 (atm)
H
constante de Henry
CCO2* concentration molaire en CO2 dans la phase liquide lquilibre
La reprsentation du systme se simplifie donc de la manire suivante :

Direction du transfert
Film
liquide

Phase
gazeuse

Phase
liquide

Cg
Cgi

CLi
CL

CL

Interface gaz-liquide

Le transfert va donc tre caractris par le coefficient de transfert ct liquide kL et obit


l'quation qui exprime la densit de flux :
.
N CO = k (C* C
)
2
L
L CO
2
En intgrant cette expression sur la hauteur de colonne, on obtient lexpression du
coefficient de transfert :
C C Le
L
)
kL=
( LS
C
H.a w A C
Avec
H : hauteur de garnissage (m)
L : dbit liquide volumique (m3.s-1)
aw: aire mouille du garnissage (m2/m3)
AC : section droite de colonne (m2)
C : potentiel dchange (mol.m-3)
kL : coefficient de transfert (m.s-1)

C =

Avec

(C

*
e

C LS (C *S C Le )
C * C LS

ln *e
C S C Le

CLe concentration molaire en CO2 dans le liquide entrant (mol.m-3)


CLs concentration molaire en CO2 dans le liquide sortant (mol.m-3)
Ce*, la concentration en CO2 du liquide en quilibre avec CGe (mol.m-3)
Cs*, la concentration en CO2 du liquide en quilibre avec CGs (mol.m-3)
Avec galement,


k L L
L g
Avec

ReL =

1/ 3

= 0,0051 (Re L )

dL

Sc L =

2/3

(Sc L )

L
DL

Et
kL: coefficient de transfert de film ct liquide (m/s)

0.5

a
w
a

(a d )

0 , 27

L: viscosit dynamique (Pa.s)


L: viscosit cinmatique (m2.s-1)
L: masse volumique (kg/m3)
L : dbit spcifique liquide dans les conditions de fonctionnement (kg/m.s)
aw: surface spcifique mouille du garnissage (m/m3)
a : surface spcifique du garnissage (m/m3)
d : diamtre du garnissage (m)
pour les anneaux Raschig a.d = 4,67
DL: coefficient de diffusion liquide (m.s-1)

II.2.2 Etudes
Ltude se fera pour un dbit liquide et un dbit total de gaz correspondant au point de
charge du dbit liquide choisit, avec 4% de CO2.
Pour la phase gaz en entre et en sortie, relever la pression totale, la temprature, et la
composition molaire.
- Dterminer la PCO2 de la phase gaz en entre et en sortie de la colonne (attendre le
rgime permanent)
- Effectuer un bilan sur le CO2 sur lensemble de la colonne pour dterminer la
concentration en CO2 dans le liquide en sortie de colonne (CLs).
- Dterminer la valeur du coefficient de transfert liquide kL. Comparer avec la valeur
thorique.

II.3. Etude du transfert de matire en absorption chimique


II.3.1 Rappel thorique
Nous tudions labsorption du CO2 dun mlange CO2-air par une solution aqueuse de
soude. Cest une absorption avec raction chimique. En effet en solutions fortement
basiques (pH>12) le CO2 ragit avec les ions OH- :
CO2 + 2 NaOH-

Na2CO3 + H2O

Le solut qui transfre est alors consomm dans la phase liquide et sa concentration
est maintenue trs faible. Ceci permet dinfluencer la fois la vitesse dabsorption de
lespce (par modification du gradient de concentration dans le film liquide) et la
capacit dabsorption du solvant. Le flux travers linterface peut-tre exprim en
fonction du facteur dacclration E selon la relation :
.
N CO = E.k (C* C
)
2
L
L CO
2

k L'
est le rapport des coefficients de transfert de film avec raction chimique ( k L' )
kL
et sans raction chimique ( k L ). Ce facteur E caractrise l'importance de la raction
chimique sur le transfert vis--vis du processus de diffusion. Le facteur dacclration E
peut avoir une expression mathmatique pour des ractions simples.

O E =

Le facteur dacclration est accessible partir des nombres adimensionnels suivants :


le nombre de Hatta, Ha. Cest la racine carre du rapport de la quantit de
gaz qui peut au maximum ragir dans le film et de la quantit de gaz qui peut
tre transporte travers le film par diffusion molculaire ;
le facteur d'acclration limite Ei : il contient les rapports des diffusivits et
des concentrations et traduit le facteur dacclration obtenu pour une raction
instantane (la vitesse globale du processus est alors limite par la vitesse de
diffusion du solvant dans le film). Il reprsente la valeur asymptotique du facteur
dacclration E lorsque Ha augmente indfiniment.
Cas particulier de cintiques dordre 2 (ou de pseudo premier ordre) :
kR
Pour la raction dordre 2 entre A gazeux et B dissous : A Gaz + BLiq
C+D
La vitesse de raction sexprime selon : r = kR[A]L [B]L o kR est la constante de vitesse
de raction.
En supposant que [B]L est constant dans toute la colonne, on se ramne alors une
cintique du pseudo premier ordre : r = kR [A]L
Dans le cas de la soude en excs, on a alors une raction rapide du pseudo premier
ordre rapide. Le flux sexprime alors de la faon suivante :

.
N CO = E.k (C* )
2
L
Avec le nombre de Hatta dfini par :
Ha 2 =

k R [B ]L D A k R ' D A
=
k 2L
k 2L

Et le facteur dacclration limite scrit quant- lui :

Ei =

D B [B]L
+1
D A CA i

Le facteur dacclration peut alors tre dduit de labaque de Van Krevelen et Hoftijzer
(1948) figure 1 - si lon connat le nombre de Hatta Ha, et Ei.

1000

1000

100

100

Ei
10

10

1
0.1

10
Ha

100

1
1000

Figure 1 : Abaque de Van Krevelen et Hoftijzer pour une raction de pseudo premier ordre

II.2.2 Etudes
Pour les mmes dbits que labsorption sans raction :
- dterminer la quantit de CO2 transfre en mesurant la variation de la concentration
en CO2 de la phase gaz
- comparer les efficacits d'absorption obtenues avec celles de l'absorption sans
raction.
- calculer le nombre de Hatta
- calculer le facteur dacclration exprimentale et thorique.

10

Donnes physico chimiques


1- Diffusivit du CO2 dans leau.
Elle varie avec la temprature suivant la loi de Nernst :

D A L
= cte
T

On obtient les valeurs suivantes (Tableau 1) pour le CO2 :


Tableau 1 : influence de la temprature sur le coefficient de diffusion du CO2 dans leau

T (K)
288
L (cP)
1,22
DA105 (cm2.s- 1,58
1
)

289
1,20
1,62

290
1,17
1,66

291
1,15
1,70

292
1,12
1,75

293
1,10
1,79

294
1,09
1,81

295
1,06
1,87

296
1,04
1,91

297
1,02
1,95

298
1,00
2

2- Diffusivit des ions HO- dans leau


On prendra DB = 1,7 DA

3- Solubilit (constante de Henry) du CO2 dans leau


Linfluence de la temprature est donne par la relation :

1
1023,5
= 4,8945 +
log
T
He'0
Linfluence de la concentration en hydroxyde de sodium sexprime selon :

He' = He'0 10 h I
O

I est la force ionique de la solution en kmol.m-3

(Rappel :

I=

1
[Ci].z i2 )

2 ions i

T est en Kelvin
He, He0 sont les constantes de Henry en atm.m3.kmol-1
h est un facteur de solubilit en m3.kmol-1
Les valeurs de h ont t calcules en prenant en compte la contribution du gaz. Ces valeurs sont
reportes en fonction de la temprature dans le Tableau 2.
Tableau 2 : influence de la temprature sur le facteur de solubilit du CO2 dans leau
T (K)
hG (m3.kmol-1)
h (m3.kmol-1)

288
-0,01
0,147

289
-0,011
0,146

290
-0,012
0,145

291
-0,013
0,144

292
-0,014
0,143

293
-0,015
0,142

294
-0,016
0,141

295
-0,017
0,140

296
-0,018
0,139

297
-0,019
0,138

298
-0,019
0,138

4- Constante cintique kR de la raction entre le CO2 et la soude


On considre que la raction est dordre 1 par rapport chacun des ractifs.
La constante de vitesse dpend de la temprature et de la force ionique de la solution selon la relation :

k R = 10
O

2850

+ 0 ,133I
13 ,4
T

T est en Kelvin
I est en kmol.m-3
kR est en m3.kmol-1.s-1

11

ANNEXE 1
Schma de linstallation

12

Vous aimerez peut-être aussi