Cours Et TD
Cours Et TD
/ Université de Médéa
CHAPITRE 01 :
Généralités et rappels de génie
de la réaction I ‘’Réacteurs 1’’
I. Introduction :
Le génie de la réaction s’intéresse à l’étude de la réaction chimique dans son aspect micro
(cinétique) et macro (réacteur) par l’élaboration des modèles mathématiques (lois cinétiques,
bilan de masse et d’énergie).
II. Types de réacteur :
Selon le mode de fonctionnement on distingue:
Réacteurs continus (ouvert) :
Dans ce cours, nous allons voir ‘’les réacteurs hétérogènes’’ catalytique et non catalytique
largement rencontrés dans l’industrie : pétrochimique, pharmaceutique, minéralogique....
1
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑛𝑛0 : Nombre de modes des réactifs actifs, il est utilisé surtout pour les systèmes complexes
(plusieurs réactions).
b) Taux de conversion d’un constituant clés X A soit A le réactif limitant, d’où la conversion
est X A .
𝑛𝑛𝐴𝐴 = 𝑛𝑛𝐴𝐴0 + 𝜈𝜈𝑖𝑖 𝑋𝑋𝐴𝐴 𝑛𝑛𝐴𝐴0 = 𝑛𝑛𝐴𝐴0 − 𝑋𝑋𝐴𝐴 𝑛𝑛𝐴𝐴0 (soit :𝜈𝜈𝐴𝐴 = −1)
2
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Hypothèses :
-mélange homogène (la concentration et la température sont le même dans tous le réacteur).
Bilan :
𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼 − 𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜 + (𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 − 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐) =Accumulation
*Pour un seul réacteur :
𝑑𝑑𝑑𝑑𝑖𝑖
0 − 0 + 𝑉𝑉. 𝑟𝑟𝑖𝑖 = Où : 𝑉𝑉. 𝑟𝑟𝑖𝑖 : temps de séjours (temps de réacteur)
𝑑𝑑𝑑𝑑
*pour plusieurs réacteurs :
𝑑𝑑𝑑𝑑𝑖𝑖
0 − 0 + 𝑉𝑉. ∑𝜈𝜈𝑖𝑖𝑖𝑖 𝑟𝑟𝑗𝑗 =
𝑑𝑑𝑑𝑑
Hypothèses :
-Mélange parfait
-Volume constant V
-Les conditions qui règnent sont de sortie (C S , T° S )
-Débit constant = Q S =Q 0.
𝑑𝑑𝑑𝑑𝑖𝑖
Bilan : En état stationnaire =0.
𝑑𝑑𝑑𝑑
𝑉𝑉
𝐹𝐹𝑖𝑖𝑖𝑖 − 𝐹𝐹𝑖𝑖𝑖𝑖 + 𝑉𝑉. 𝑟𝑟𝑖𝑖 = 0 , 𝐹𝐹𝑖𝑖 = 𝑄𝑄0 . 𝐶𝐶𝑖𝑖 𝑒𝑒𝑒𝑒 𝜏𝜏 = :temps de passage
𝑄𝑄0
3
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Hypothèses :
Bilan :
4
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
CHAPITRE 02 :
Partie 01 : Réacteurs fluides-fluides : Réactions
Gaz-liquide
I. Introduction :
Les réactions gaz-liquide sont multiples, et ils sont présentés dans plusieurs industries :
Chloration, Hydrogénation, Oxydation. 𝜈𝜈𝐴𝐴 𝐴𝐴(𝑔𝑔) + ⋯ → 𝐴𝐴(𝑙𝑙)
5
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Hypothèses :
-Ecoulement laminaire dans chaque filme (pas de turbulence filme stagnant).
-Transfert de masse gouverné par la diffusion unidirectionnelle et perpendiculaire à l’interface.
-Au-delà de chaque filme, la composition de chaque belk est considéré comme uniforme.
-A l’interface L-G les deux phases sont en équilibre thermodynamique.
*Comme résultat de cette théorie ou modèle on a :
a /- Les coefficients de transfert de matière dans chaque filme est :
𝐷𝐷𝐴𝐴𝐴𝐴 𝐷𝐷𝐴𝐴𝐴𝐴
𝐾𝐾𝐶𝐶 = 𝑒𝑒𝑒𝑒 𝐾𝐾𝐿𝐿 =
𝜕𝜕𝐺𝐺 𝜕𝜕𝐿𝐿
Avec :
𝐷𝐷𝐴𝐴𝐴𝐴 , 𝐷𝐷𝐴𝐴𝐴𝐴 : Diffusivité de A dans le film gazeux (G) et liquide (L) respectivement.
𝜕𝜕𝐺𝐺 , 𝜕𝜕𝐿𝐿 :Epaisseur du filme G et du filme L respectivement.
𝐾𝐾𝐺𝐺 , 𝐾𝐾𝐿𝐿 :Coefficient de T.M de A dans la phase G et la phase L respectivement.
La vitesse de transfert du composé A de la phase G vers la phase L est donné par r= Aire
spécifique x Flux surfacique
6
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Avec :
1 1 𝐾𝐾𝐻𝐻
= + : Résistance globale aux T.M dans tout le système Gaz-liquide.
𝑎𝑎.𝐾𝐾𝐺𝐺𝐺𝐺 𝑎𝑎.𝐾𝐾𝐺𝐺 𝑎𝑎.𝐾𝐾𝐿𝐿
*Remarque : le modèle du double filmes est le plus simple parmi d’autre tels le : modèle de
pénétration de Higbie, ou du renouvellement de surface de Dankwert.
a /- Facteur de Hatta :
*Selon l’équation (1), les conditions aux limites sont :
𝑑𝑑2 𝐶𝐶𝐴𝐴
𝐷𝐷𝐴𝐴 . = −𝑟𝑟𝐴𝐴 = 𝑟𝑟 ; 𝑍𝑍 = 0 → 𝐶𝐶𝐴𝐴 = 𝐶𝐶𝐴𝐴𝐴𝐴 𝑒𝑒𝑒𝑒 à 𝑍𝑍 = 𝛿𝛿 → 𝐶𝐶𝐴𝐴 = 𝐶𝐶𝐴𝐴𝐴𝐴
𝑑𝑑𝑑𝑑 2
𝑑𝑑2 𝐶𝐶𝐵𝐵 𝑑𝑑𝑑𝑑𝐵𝐵
𝐷𝐷𝐵𝐵 . = 𝛾𝛾. 𝑟𝑟 ; 𝑍𝑍 = 0 → = 0 𝑒𝑒𝑒𝑒 à 𝑍𝑍 = 𝛿𝛿 → 𝐶𝐶𝐵𝐵 = 𝐶𝐶𝐵𝐵𝐵𝐵
𝑑𝑑𝑑𝑑 2 𝑑𝑑𝑑𝑑
*La résolution de ces équations différentielles dépend de la forme de r et elle est facile
seulement pour des réactions d’ordre 1 par rapport à A.
*Dans ce contexte, on introduit le nombre ou le facteur de Hatta dont la forme générale est :
1 2 𝑚𝑚−1
𝐻𝐻𝐻𝐻 = � . 𝑘𝑘. 𝐷𝐷𝐴𝐴𝐴𝐴 . 𝐶𝐶𝐴𝐴𝐴𝐴 . 𝐶𝐶𝐵𝐵𝐵𝐵
𝐾𝐾𝐿𝐿 𝑚𝑚 + 1
7
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Avec :
*Dans ce qui suit, on traite les réactions chimiques avec l’ordre m=1 par rapport à A. D’où le
𝑛𝑛 .𝐷𝐷
�𝑘𝑘.𝐶𝐶𝐵𝐵𝐵𝐵 𝐴𝐴𝐴𝐴 𝑛𝑛
𝑘𝑘.𝐶𝐶𝐵𝐵𝐵𝐵
facteur Hatta est : 𝐻𝐻𝐻𝐻 = = 𝛿𝛿. �
𝐾𝐾𝐴𝐴𝐴𝐴 𝐷𝐷𝐴𝐴𝐴𝐴
b/-Facteur d’accélération :
Nommé aussi facteur d’intensification E a exprime l’amélioration du flux dans le film liquide :
𝜑𝜑 = 𝐸𝐸𝑎𝑎 . 𝐾𝐾𝐿𝐿 . (𝐶𝐶𝐴𝐴𝐴𝐴 − 𝐶𝐶𝐴𝐴𝐴𝐴 )
𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹 𝑑𝑑𝑑𝑑 𝐴𝐴 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓ℎ𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑙𝑙 ′ 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑒𝑒𝑒𝑒 𝑝𝑝𝑝𝑝é𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑟𝑟é𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎
𝐸𝐸𝑎𝑎 =
𝐹𝐹𝐹𝐹𝐹𝐹𝐹𝐹 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑑𝑑𝑑𝑑 𝐴𝐴 𝑒𝑒𝑒𝑒 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑑𝑑𝑑𝑑 𝑙𝑙𝑙𝑙 𝑟𝑟é𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 𝑐𝑐ℎ𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖
*Le calcul de 𝐸𝐸𝑎𝑎 est très difficile (𝐶𝐶𝐴𝐴𝐴𝐴 ??) , d’où l’utilisation d’un nombre 𝐸𝐸𝐴𝐴𝐴𝐴 ‘’Facteur
d’accélération limite’’.
𝐷𝐷𝐵𝐵𝐵𝐵 . 𝐶𝐶𝐵𝐵𝐵𝐵
𝐸𝐸𝐴𝐴𝐴𝐴 = 1 +
𝜈𝜈. 𝐷𝐷𝐴𝐴𝐴𝐴 . 𝐶𝐶𝐴𝐴𝐴𝐴
*En utilisant le diagramme ‘’Van Krevelen-Hoftyzer’’ 𝐸𝐸𝑎𝑎 = 𝑓𝑓(𝐻𝐻𝐻𝐻) dans lequel l’assymtote
𝐻𝐻𝐻𝐻 → ∞ donne la valeur de 𝐸𝐸𝐴𝐴𝐴𝐴 .
8
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑎𝑎.𝐾𝐾𝐵𝐵𝐵𝐵 .𝐶𝐶𝐵𝐵𝐵𝐵
*Si 𝐶𝐶𝐵𝐵𝐵𝐵 est élevée ⟹ > 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 . 𝑃𝑃𝐴𝐴𝐴𝐴 … . . (∗∗).
𝜈𝜈
Remarque :
*Si les valeurs de 𝐶𝐶𝐵𝐵𝐵𝐵 donnent les cas (∗) et (∗∗) en même temps ⟹Il y a un point 𝐶𝐶′𝐵𝐵 , 𝑃𝑃′𝐴𝐴 )
point de changement de régime (ou de transition) qu’il faut déterminer en
𝑎𝑎.𝐾𝐾𝐵𝐵𝐵𝐵 .𝐶𝐶’𝐵𝐵𝐵𝐵
utilisant : = 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 . 𝑃𝑃𝐴𝐴𝐴𝐴
𝜈𝜈
9
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
10
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
11
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
*Il faut calculer𝐻𝐻𝑎𝑎 ; ensuite 𝐸𝐸𝐴𝐴𝐴𝐴 (en utilisant l’hypothèse 𝑃𝑃𝑖𝑖𝑖𝑖𝑖𝑖è𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 = 𝑃𝑃𝐴𝐴𝐴𝐴 = 𝐾𝐾𝐻𝐻 . 𝐶𝐶𝐴𝐴𝐴𝐴 )⟹
comparer 𝐻𝐻𝐻𝐻 avec les valeurs de 𝐸𝐸𝐴𝐴𝐴𝐴 ⟹ lire le régime qui correspond à la valeur de Hatta
calculée.
𝐷𝐷𝐵𝐵𝐵𝐵
* La valeur de 𝐷𝐷𝐵𝐵𝐵𝐵 peut etre déduite de la formule : 𝐾𝐾𝐵𝐵𝐵𝐵 = aussi 𝐷𝐷𝐴𝐴𝐴𝐴 = 𝐾𝐾𝐵𝐵𝐵𝐵 . 𝛿𝛿
𝛿𝛿
*Cas particulier :
Souvent les 08 régimes de tableau sont classés selon les valeurs référentielles de Ha 2 et
0.02 comme suivant :
• Si 𝐻𝐻𝐻𝐻 > 2 : Réaction dans le film liquide ou à l’interface (les 04 premier régimes).
• Si 0.02< 𝐻𝐻𝐻𝐻 < 2 : Réaction de vitesse moyenne à lente dans le film liquide et dans le
bulk liquide (les régimes de la ligne 5.6 et 7 dans le tableau).
• Si 𝐻𝐻𝐻𝐻 < 0.02 : réaction très lente (infiniment lente) dans le bulk liquide seulement
(dernier régime).
12
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
EXERCICE 1 :
Soit la réaction gaz-liquide suivante : 𝐴𝐴(𝑔𝑔) + 2𝐵𝐵(𝑙𝑙) → 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝(𝑙𝑙)
𝑟𝑟 = 𝐾𝐾. 𝐶𝐶𝐴𝐴 . 𝐶𝐶𝐵𝐵2 et 𝐾𝐾 = 106 m6/mol2.h
Avec : 𝑃𝑃𝐴𝐴 = 100𝑃𝑃𝑃𝑃 ; 𝐶𝐶𝐵𝐵 =100mol/m3 ; 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 =10-2 mole/m3. h. Pa ; 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 =20 h-1
𝑎𝑎=20m-1 ; 𝐾𝐾𝐻𝐻 =105 Pa.m3/mol ; 𝑓𝑓𝐿𝐿 =0.98 ; 𝐷𝐷𝐴𝐴𝐴𝐴 = 𝐷𝐷𝐵𝐵𝐵𝐵 =10-6 m2/h.
1. Déterminer le régime qui contrôle la réaction.
2. Calculer la vitesse de la réaction.
3. Calculer en pourcentage toutes les résistances au transfert de masse.
4. Refaite les mêmes questions si 𝐶𝐶𝐵𝐵 =1mol/m3
EXERCICE 2 :
Refaite le même EXERCICE 1 pour les données suivantes :
𝑟𝑟 = 𝐾𝐾. 𝐶𝐶𝐴𝐴 . 𝐶𝐶𝐵𝐵 Et 𝑃𝑃𝐴𝐴 = 100𝑃𝑃𝑃𝑃 , 𝐶𝐶𝐵𝐵 =100mol/m3
1. Avec k= 10m3/mol.h ; 𝐾𝐾𝐻𝐻 =103 Pa.m3
2. Avec k= 10-2m3/mol.h ; 𝐾𝐾𝐻𝐻 =1 Pa.m3
On donne : 𝑃𝑃𝐴𝐴 = 105 𝑃𝑃𝑃𝑃 , 𝐶𝐶𝐵𝐵 =1 mol/m3, 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 =10-1 mole/m3. h. Pa, 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 102 h-1
13
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
CHAPITRE 02 :
Partie 02 : Réacteurs G-L
I. Introduction :
Les réacteurs G-L appelés aussi contacteurs ou absorbeurs sont nombreux dans plusieurs
industries. Ils sont classés en deux types : colonnes (trous) et cuves (Tanks).
14
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
15
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑑𝑑𝑑𝑑𝐴𝐴 𝐹𝐹𝐿𝐿
𝐹𝐹𝐺𝐺 . = . 𝑑𝑑𝑑𝑑𝐴𝐴 = 𝑎𝑎. 𝐾𝐾𝐺𝐺𝐺𝐺 (𝑃𝑃𝐴𝐴 − 𝑃𝑃𝐴𝐴∗ ). 𝑑𝑑𝑑𝑑𝑅𝑅
𝜋𝜋 𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
16
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
D’où :
𝐹𝐹𝐺𝐺 𝑃𝑃 𝑑𝑑𝑑𝑑𝐴𝐴 𝐹𝐹𝐿𝐿 𝐶𝐶𝐴𝐴2 𝑑𝑑𝑑𝑑𝐴𝐴
𝑉𝑉𝑅𝑅 = 𝑆𝑆. ℎ = . ∫𝑃𝑃 𝐴𝐴2 ∗ = ∫𝐶𝐶
𝜋𝜋.𝑎𝑎.𝐾𝐾𝐺𝐺𝐺𝐺 𝐴𝐴1 𝑃𝑃𝐴𝐴 −𝑃𝑃𝐴𝐴 𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇 .𝑎𝑎.𝐾𝐾𝐿𝐿𝐿𝐿 𝐴𝐴1 𝐶𝐶𝐴𝐴∗ −𝐶𝐶𝐴𝐴
Avec :
Remarque : Dans le cas d’une absorption physique (sans réaction chimique) le calcul des
performance ℎ ou 𝑉𝑉𝑅𝑅 nécessite l’utilisation des équations (3) et celle du bilan global.
17
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑑𝑑𝑑𝑑𝐴𝐴 𝑑𝑑𝑑𝑑𝐵𝐵
*L’équation bilan devient en substituant 𝑑𝑑𝑑𝑑𝐴𝐴 par ( ) et 𝑑𝑑𝑑𝑑𝐵𝐵 par ( ):
𝜋𝜋 𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
𝐹𝐹𝐺𝐺 −𝐹𝐹𝐿𝐿
𝑟𝑟. 𝑑𝑑𝑑𝑑𝑅𝑅 = 𝑑𝑑𝑑𝑑𝐴𝐴 = 𝑑𝑑𝑑𝑑𝐵𝐵 ……… (4)
𝜋𝜋 𝜈𝜈.𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
*L’équation (5) est dite équation de performance qui permit de calculer 𝑉𝑉𝑅𝑅 ou ℎ de la colonne.
De plus, l’équation bilan de masse sur toute la colonne est :
𝐹𝐹𝐺𝐺 −𝐹𝐹𝐿𝐿
(𝑃𝑃𝐴𝐴2 − 𝑃𝑃𝐴𝐴1 ) = (𝐶𝐶𝐴𝐴2 − 𝐶𝐶𝐴𝐴1 ) …………… (6)
𝜋𝜋 𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
Remarque :
• Le volume du réacteur 𝑉𝑉𝑅𝑅 peut etre obtenu graphiquement par l’aire sous la courbe :
1 1
𝑓𝑓(𝑃𝑃𝐴𝐴 ) Ou bien 𝑓𝑓(𝐶𝐶𝐵𝐵 ).
𝑟𝑟 𝑟𝑟
Où ; aire est donnée par :
1 𝑉𝑉𝑅𝑅 .𝜋𝜋
𝑓𝑓(𝑃𝑃𝐴𝐴 )→Aire=
𝑟𝑟 𝐹𝐹𝐺𝐺
1 −𝜈𝜈.𝑉𝑉𝑅𝑅 .𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
𝑓𝑓(𝐶𝐶𝐵𝐵 )→Aire=
𝑟𝑟 𝐹𝐹𝐿𝐿
Et en utilisant les conversions 𝑦𝑦𝐴𝐴 et 𝑥𝑥𝐵𝐵 :
1 𝑉𝑉𝑅𝑅
𝑓𝑓(𝑦𝑦𝐴𝐴 )→Aire=
𝑟𝑟 𝐹𝐹𝐺𝐺
1 −𝑉𝑉𝑅𝑅 .𝜈𝜈
𝑓𝑓(𝑥𝑥𝐵𝐵 )→Aire=
𝑟𝑟 𝐹𝐹𝐿𝐿
𝐹𝐹𝐿𝐿
𝐹𝐹𝐺𝐺 . 𝑑𝑑𝑑𝑑𝐴𝐴 = 𝑑𝑑𝑑𝑑𝐵𝐵 = 𝑟𝑟. 𝑑𝑑𝑑𝑑𝑅𝑅
𝜈𝜈
18
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞é 𝑑𝑑𝑑𝑑
𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞é 𝑑𝑑𝑑𝑑
𝐵𝐵 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝
• Phase liquide : � � = �𝐵𝐵 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡é𝑒𝑒�
𝑝𝑝𝑝𝑝𝑝𝑝 𝑝𝑝ℎ𝑎𝑎𝑎𝑎𝑎𝑎
𝑝𝑝𝑝𝑝𝑝𝑝 𝑟𝑟é𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎
𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿
1 𝑟𝑟
𝐹𝐹𝐿𝐿 (𝑥𝑥𝐵𝐵𝐵𝐵 − 𝑥𝑥𝐵𝐵𝐵𝐵 ) = 𝐿𝐿�𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 . 𝑉𝑉𝑅𝑅 ≠ 𝑟𝑟�𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 . 𝑉𝑉𝑅𝑅
𝜈𝜈
𝑟𝑟𝐿𝐿 : vitesse de transformation de B dans la phase liquide seulement,( 𝑟𝑟𝐿𝐿 ≠ 𝑟𝑟). Elle ne tient pas
compte de ce qui se passe dans G.
∗ Pour tout le réacteur (phase G + phase L) on a :
𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞é 𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 𝑑𝑑𝑑𝑑 𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞𝑞é 𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔𝑔 𝑑𝑑𝑑𝑑
1
� 𝐴𝐴 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑝𝑝𝑝𝑝𝑝𝑝 � = � 𝐵𝐵 𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝𝑝 𝑝𝑝𝑝𝑝𝑝𝑝 �
𝜈𝜈
𝐺𝐺 𝐿𝐿
𝐹𝐹𝐿𝐿
𝐹𝐹𝐺𝐺 . (𝑦𝑦𝐴𝐴𝐴𝐴 − 𝑦𝑦𝐴𝐴𝐴𝐴 ) = (𝑥𝑥 − 𝑥𝑥𝐵𝐵𝐵𝐵 )
𝜈𝜈 𝐵𝐵𝐵𝐵
𝐹𝐹𝐺𝐺 𝐹𝐹𝐿𝐿
Pour un système dilué : (𝑃𝑃𝐴𝐴𝐴𝐴 − 𝑃𝑃𝐴𝐴𝐴𝐴 ) = (𝐶𝐶𝐵𝐵𝐵𝐵 − 𝐶𝐶𝐵𝐵𝐵𝐵 ) …. (*)
𝜋𝜋 𝜈𝜈𝐶𝐶𝑇𝑇𝑇𝑇𝑇𝑇
19
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
2/ calculer les valeurs de Hatta (Ha), 𝐸𝐸𝐴𝐴𝐴𝐴 et ensuite 𝐸𝐸𝑎𝑎 et 𝑟𝑟 (parfois il faut utiliser la méthode
essaie et erreur.)
𝐹𝐹𝐺𝐺
3/ comparer le terme ( (𝑃𝑃𝐴𝐴𝐴𝐴 − 𝑃𝑃𝐴𝐴𝐴𝐴 )) avec le terme (𝑟𝑟�𝑃𝑃 . 𝑉𝑉𝑅𝑅 ) s’il sont égaux ou non, et
𝜋𝜋 𝐴𝐴𝐴𝐴
réessayer en modifiant la valeur de 𝑃𝑃𝐴𝐴𝐴𝐴 jusqu’à ce que l’égalité se réalise entre les deux termes
et noter 𝐶𝐶𝐵𝐵0 et 𝐶𝐶𝐵𝐵𝐵𝐵 .
−VL dCB
4/ en utilisant 𝐶𝐶𝐵𝐵0 et 𝐶𝐶𝐵𝐵𝐵𝐵 et . = 𝑟𝑟�𝑃𝑃 . 𝑉𝑉𝑅𝑅 , vous obtenez :
ν dt 𝐴𝐴𝐴𝐴
20
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
1
L’équation (∗) est généralement calculée graphiquement par = f(CB ) pour déterminer l’air
r
t.ν
sous la courbe = .
fL
Remarques :
- Le temps minimal pour que toute la quantité de A réagie est noté 𝑡𝑡𝑚𝑚𝑚𝑚𝑚𝑚 pour un système dilué
−VL .π(CBF −CB0 )
(𝑃𝑃𝐴𝐴 << 𝜋𝜋) ,on a t min = .
ν.FG .PAe
tmin
- L’efficacité de l’utilisation de la réaction chimique avec A est exprimée par : ɱ =
tutile
21
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
1/ le coefficient de transfert de matière dans la phase liquide est 5.10-5 m/s , quel type de
réacteur suggérer vous pour un tel système ?
Exercice 02 :
La concentration d’une impureté A dans l’air (105 Pa) doit être abaissé de 100 Pas vers 20 Pas
par absorption dans une colonne à contre-courant en utilisant l’eau pure.
Déterminer la hauteur de la colonne nécessaire.
On donne : 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 0.32𝑚𝑚𝑚𝑚𝑚𝑚/𝑚𝑚3 . ℎ. 𝑃𝑃𝑃𝑃 ; 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 0.1ℎ−1 ; 𝐾𝐾𝐻𝐻 = 12.5𝑃𝑃𝑃𝑃. 𝑚𝑚3 /𝑚𝑚𝑚𝑚𝑚𝑚
; 𝐹𝐹𝐿𝐿/𝑆𝑆 = 7. 105 𝑚𝑚𝑚𝑚𝑚𝑚/𝑙𝑙. 𝑚𝑚2 ; 𝐹𝐹𝐺𝐺/𝑆𝑆 = 105 𝑚𝑚𝑚𝑚𝑚𝑚 ⁄ 𝑚𝑚2 . ℎ ; 𝐶𝐶𝑡𝑡𝑡𝑡𝑡𝑡 = 56. 103 𝑚𝑚𝑚𝑚𝑚𝑚 ⁄𝑚𝑚3
22
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑘𝑘 = 10 𝑚𝑚3 /𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚. 𝑠𝑠
Exercice 03 :
Un polluant gazeux A de pression P A égale à 100pa dans l’air est récupéré à l’aide d’une
réaction chimique instantanée réalisée dans une colonne à contre-courant en utilisant un
réactif B de concentration C B. 𝐴𝐴(𝑔𝑔) + 𝐵𝐵(𝑙𝑙) → 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃(𝑙𝑙)
Les critères de sécurité exigent que la pression de ce gaz soit égale à 20 pa, Déterminer la
hauteur de la colonne nécessaire pour respecter ce critère.
1/ si la concentration initiale de réactif B est 0.8 mol/l.
2/ si la concentration initiale de réactif B est 0.032mol/l.
3/ si la concentration initiale de réactif B est 0.128 mol/l.
On donne :
𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 0.32𝑚𝑚𝑚𝑚𝑚𝑚/𝑚𝑚3 . ℎ. 𝑃𝑃𝑃𝑃 ; 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 0.1ℎ−1 ; 𝐾𝐾𝐻𝐻 = 12.5𝑃𝑃𝑃𝑃. 𝑚𝑚3 /𝑚𝑚𝑚𝑚𝑚𝑚 ; 𝑘𝑘 = 10 𝑚𝑚3 /
𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚. 𝑠𝑠 ;𝐷𝐷𝐴𝐴𝐴𝐴 = 𝐷𝐷𝐵𝐵𝐵𝐵 ;
𝐹𝐹𝐿𝐿/𝑆𝑆 = 7. 105 𝑚𝑚𝑚𝑚𝑚𝑚/ℎ. 𝑚𝑚2 ; 𝐹𝐹𝐺𝐺/𝑆𝑆 = 105 𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚2 . ℎ ; 𝐶𝐶𝑡𝑡𝑡𝑡𝑡𝑡 = 56. 103 𝑚𝑚𝑚𝑚𝑚𝑚 ⁄𝑚𝑚3 ; 𝜋𝜋 = 1𝑏𝑏𝑏𝑏𝑏𝑏
Exercice 04 :
Un composé gazeux A est introduit à une fraction 0.01 dans une cuve agitée et subit une
réaction extrêmement rapide sous les conditions suivantes : 𝐴𝐴(𝑔𝑔) + 𝐵𝐵(𝑙𝑙) → 𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃(𝑙𝑙)
𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 0.72𝑚𝑚𝑜𝑜𝑜𝑜/𝑚𝑚3 . ℎ. 𝑃𝑃𝑃𝑃 ; 𝑎𝑎. 𝐾𝐾𝐴𝐴𝐴𝐴 = 144ℎ−1 ; 𝐾𝐾𝐻𝐻 = 0.1𝑃𝑃𝑃𝑃. 𝑚𝑚3 /𝑚𝑚𝑚𝑚𝑚𝑚 ; 𝜋𝜋 = 1𝑎𝑎𝑎𝑎𝑎𝑎 ;
𝑎𝑎 = 200𝑚𝑚−1 ; 𝐹𝐹𝐿𝐿 = 9. 103 𝑚𝑚𝑚𝑚𝑚𝑚/ℎ ; 𝐹𝐹𝐺𝐺 = 9.104 𝑚𝑚𝑚𝑚𝑚𝑚/ℎ ; 𝐶𝐶𝐵𝐵 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒é = 5.56. 103 𝑚𝑚𝑚𝑚𝑚𝑚/𝑚𝑚3 ; 𝐶𝐶𝑡𝑡𝑡𝑡𝑡𝑡 =
56. 103 𝑚𝑚𝑚𝑚𝑚𝑚 ⁄𝑚𝑚3 ; 𝐷𝐷𝐴𝐴𝐴𝐴 = 𝐷𝐷𝐵𝐵𝐵𝐵 .
1/ montrer que la réaction est contrôlée par la résistance dans le film gazeux.
23
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
CHAPITRE 03 :
Réacteurs G-S (Résumé)
−1 𝑑𝑑𝑑𝑑𝐴𝐴
𝑟𝑟𝐴𝐴 = . = 𝑘𝑘𝑆𝑆 . 𝐶𝐶𝐴𝐴𝐴𝐴
𝑆𝑆𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑
Bilan de masse :
24
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
−1 𝑑𝑑𝑑𝑑𝐴𝐴
𝑄𝑄𝑔𝑔 = = 𝐾𝐾𝑔𝑔 (𝐶𝐶𝐴𝐴𝐴𝐴 − 𝐶𝐶𝐴𝐴𝐴𝐴 )
𝛿𝛿𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑
III. Etude d’un réaction G-S : solide de taille constante (sans changement de la taille) :
Le réactif A (gazeux) réagit avec le solide B sur sa surface en laissant accumuler le produit solide
sur le cœur de réactif B le produit est appelé « cendres ».
Soit la réaction d’ordre 1 : 𝐴𝐴(𝑔𝑔) + 𝑏𝑏𝑏𝑏(𝑠𝑠) → 𝑃𝑃(𝑠𝑠)
On distingue les 03 cas ou régimes :
𝑉𝑉 4
On a 𝑑𝑑𝑑𝑑𝐵𝐵 = 𝑏𝑏𝑑𝑑𝑑𝑑𝐴𝐴 𝑒𝑒𝑒𝑒 𝑁𝑁𝐵𝐵 = 𝜌𝜌𝐵𝐵 . = 𝜋𝜋𝑟𝑟𝑐𝑐3
𝑀𝑀𝐵𝐵 3
𝑁𝑁𝐴𝐴 𝑒𝑒𝑒𝑒 𝑁𝑁𝐵𝐵 : nombre de mole de B et A à un temps t>0.
𝜌𝜌𝐵𝐵 ; Masse volumique de B.
𝑀𝑀𝐵𝐵 : Masse molaire de B.
−1 𝑑𝑑𝑑𝑑𝐵𝐵 −𝑏𝑏 𝑑𝑑𝑑𝑑𝐴𝐴
Aussi : = = 𝑏𝑏. 𝐾𝐾𝑔𝑔 . 𝐶𝐶𝐴𝐴𝐴𝐴 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐
𝛿𝛿𝑒𝑒𝑒𝑒𝑒𝑒 𝑑𝑑𝑑𝑑 4𝜋𝜋𝜋𝜋 2 𝑑𝑑𝑑𝑑
25
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
−𝜌𝜌𝐵𝐵 𝑑𝑑𝑑𝑑𝑐𝑐
= . 𝑟𝑟𝑐𝑐2 . = 𝑏𝑏. 𝐾𝐾𝑔𝑔 . 𝐶𝐶𝐴𝐴𝐴𝐴
𝑀𝑀𝐵𝐵 .𝑅𝑅 2 𝑑𝑑𝑑𝑑
𝑡𝑡 −𝜌𝜌𝐵𝐵 𝑟𝑟
∫0 𝑑𝑑𝑑𝑑 = 𝑀𝑀 .∫𝑅𝑅 𝑐𝑐 𝑟𝑟𝑐𝑐2 𝑑𝑑𝑑𝑑𝑐𝑐
𝑏𝑏 .𝑏𝑏.𝐾𝐾𝑔𝑔 .𝐶𝐶𝐴𝐴𝐴𝐴 .𝑅𝑅 2
−𝜌𝜌𝐵𝐵
𝑡𝑡 = . ( 𝑟𝑟𝑐𝑐3 − 𝑅𝑅3 )
3.𝑏𝑏.𝑀𝑀𝐵𝐵 .𝐾𝐾𝑔𝑔 .𝐶𝐶𝐴𝐴𝐴𝐴 .𝑅𝑅 2
𝜌𝜌𝐵𝐵 .𝑅𝑅 𝑟𝑟 3
= . �1 − � 𝑐𝑐� �
3.𝑏𝑏.𝑀𝑀𝐵𝐵 .𝐾𝐾𝑔𝑔 .𝐶𝐶𝐴𝐴𝐴𝐴 . 𝑅𝑅
-Pour une réaction complète (𝑟𝑟𝑐𝑐 = 0) le temps de réaction totale est noté 𝜏𝜏.
𝜌𝜌𝐵𝐵
𝜏𝜏 = . 𝑅𝑅
3. 𝑏𝑏. 𝑀𝑀𝐵𝐵 . 𝐾𝐾𝑔𝑔 . 𝐶𝐶𝐴𝐴𝐴𝐴
𝑡𝑡 𝑟𝑟 3
-on a : = 1 − � 𝑐𝑐�
𝜏𝜏 𝑅𝑅
4� 𝜋𝜋𝑟𝑟 3 𝑟𝑟 3
3 𝑐𝑐
Et : 1 − 𝑋𝑋𝐵𝐵 = 4� 𝑅𝑅 3 = � 𝑐𝑐 � C’est la fraction de B qui n’a pas réagi.
3 𝑅𝑅
𝑟𝑟 3 𝑡𝑡
Donc : 𝑋𝑋𝐵𝐵 = 1 − � 𝑐𝑐 � =
𝑅𝑅 𝜏𝜏
−1 𝑑𝑑𝑑𝑑𝐴𝐴 𝑑𝑑𝑑𝑑𝐴𝐴
. = 𝐷𝐷 . ; D : coefficient de diffusion
4.𝜋𝜋.𝑟𝑟 2 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
−𝑑𝑑𝑑𝑑𝐴𝐴 𝑅𝑅 𝑑𝑑𝑑𝑑 𝐶𝐶
D’où : . ∫𝑟𝑟 = 4𝜋𝜋. 𝐷𝐷=∫𝐶𝐶 𝐴𝐴𝐴𝐴=0
=𝐶𝐶
𝑑𝑑𝑑𝑑𝐴𝐴
𝑑𝑑𝑑𝑑 𝑐𝑐 𝑟𝑟 2 𝐴𝐴𝐴𝐴 𝐴𝐴𝐴𝐴
−𝑑𝑑𝑑𝑑𝐴𝐴 1 1
� − � = 4𝜋𝜋. 𝐷𝐷. 𝐶𝐶𝐴𝐴𝐴𝐴 … (∗)
𝑑𝑑𝑑𝑑 𝑟𝑟𝑐𝑐 𝑅𝑅
𝑑𝑑𝑑𝑑𝐵𝐵 4𝜋𝜋𝜌𝜌𝐵𝐵
Substituants 𝑑𝑑𝑑𝑑𝐴𝐴 par = . 𝑟𝑟𝑐𝑐2 . 𝑑𝑑𝑑𝑑𝑐𝑐
𝑏𝑏 3.𝑏𝑏.𝑀𝑀𝐵𝐵
26
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝜌𝜌𝐵𝐵 .𝑅𝑅 2 𝑟𝑟 2 𝑟𝑟 3
L’équation (∗) ⟹ = . �1 − 3 � 𝑐𝑐� +2 � 𝑐𝑐� �
6.𝑏𝑏.𝑀𝑀𝐵𝐵 .𝐷𝐷.𝐶𝐶𝐴𝐴𝐴𝐴 𝑅𝑅 𝑅𝑅
𝜌𝜌𝐵𝐵
∗ 𝑟𝑟𝑐𝑐 = 0 ⟹ 𝜏𝜏 = . 𝑅𝑅2
6. 𝑏𝑏. 𝑀𝑀𝐵𝐵 . 𝐷𝐷. 𝐶𝐶𝐴𝐴𝐴𝐴
𝑟𝑟 3 𝑡𝑡 2�
∗ 1 − 𝑋𝑋𝐵𝐵 = � 𝑐𝑐� 𝑒𝑒𝑒𝑒 ∶ = 1 − 3(1 − 𝑋𝑋𝐵𝐵 ) 3 + 2(1 − 𝑋𝑋𝐵𝐵 )
𝑅𝑅 𝜏𝜏
𝜌𝜌𝐵𝐵
𝑡𝑡 = . (𝑅𝑅 − 𝑟𝑟𝑐𝑐 )
𝑏𝑏.𝑀𝑀𝐵𝐵 .𝐾𝐾𝑆𝑆 .𝐶𝐶𝐴𝐴𝐴𝐴
𝜌𝜌𝐵𝐵
𝜏𝜏 = . 𝑅𝑅
𝑏𝑏.𝑀𝑀𝐵𝐵 .𝐾𝐾𝑆𝑆 .𝐶𝐶𝐴𝐴𝐴𝐴
𝑟𝑟 3
ET : 1 − 𝑋𝑋𝐵𝐵 = � 𝑐𝑐�
𝑅𝑅
𝑡𝑡 1�
= 1 − (1 − 𝑋𝑋𝐵𝐵 ) 3
𝜏𝜏
27
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
𝑅𝑅𝑚𝑚𝑚𝑚𝑚𝑚 𝐹𝐹(𝑅𝑅𝑖𝑖 )
Et : ����
1 − 𝑋𝑋𝐵𝐵 = ∑𝑅𝑅(𝜏𝜏>𝐸𝐸𝑆𝑆 )[1 − 𝑋𝑋𝐵𝐵 (𝑅𝑅𝑖𝑖 )]. pour un piston.
𝐹𝐹
Où : 𝐹𝐹 = ∑ 𝐹𝐹(𝑅𝑅𝑖𝑖 )
𝐹𝐹(𝑅𝑅𝑖𝑖 ): Débit de chaque fraction de rayon 𝑅𝑅𝑖𝑖 .
𝑋𝑋𝐵𝐵 (𝑅𝑅𝑖𝑖 ): Taux de conversion de la fraction de rayon 𝑅𝑅𝑖𝑖 .
Série N° :02
Chapitre 03 : Réacteur G-S
Exercice 01 :
Soit la réaction : 2𝑍𝑍𝑍𝑍𝑍𝑍(𝑠𝑠) + 3𝑂𝑂2 (𝑔𝑔) → 2𝑍𝑍𝑍𝑍𝑍𝑍 + 2𝑆𝑆𝑆𝑆2 réalisée dans l’air à 900°C et 1 atm, la
réaction est du premier ordre (k=2.10-2m/s) si la diffusion externe est négligée.
1/ quelle est l’étape limitante du processus si R i =10-4 m et 𝜌𝜌 =4130 kg/m3 et la diffusivité de O 2
est 8.10-6 m2/s.
2/ calculer le temps de réaction complète.
3/ calculer X B pour que les particules auront seulement la moitié de leur rayon initial.
On donne : M(ZnS)=97.4
Exercice 02 :
Des particules solides de rayon 4mm sont convertis à 87.5% pendant 5 min sans changement de
����
leur taille, si l’étape limitante est la diffusion interne, Calculer la conversion( 𝑋𝑋 𝐵𝐵 )moyenne que
on obtiendrait dans un réacteur à lit fluidisé opérant dans les mêmes conditions précédentes,
mais en utilisant une phase solide constituée de 50% de particules R 1 =1mm et 50% de particules
R 2 =2mm. Le débit d’alimentation est 10kg/min avec une charge de 300kg de chaque type de
particules (R 1 et R 2 ).
28
Conception des réacteurs/ Dr. Khalladi R./ Université de Médéa
Exercice 03 :
Des particules solides de rayon 4mm sont convertis à 87.5% pendant 5 min sans changement de
leur taille ; si la vitesse de la réaction est contrôlée par la diffusion dans les cendres, calculer la
conversion( 𝑋𝑋 ����
𝐵𝐵 ) moyenne que on obtiendrait dans un réacteur à lit fluidisé opérant dans les
mêmes conditions, mais en utilisant des solides constitués de 50% de particules R 1 =1mm et 50%
de particules R 2 =2mm. Le réacteur est alimenté avec un débit 10kg/min ayant une même charge
solide 300kg pour toutes les particules.
29