Questionnaire

Contrôle périodique
SOLUTION

GCH3110
Sigle du cours
Identification de l’étudiant(e)
Nom : Prénom :
Signature : Matricule : Groupe :

Sigle et titre du cours Groupe Trimestre

GCH3110 – Calcul des réacteurs
01 Été 2020
chimiques
Chargé de cours Local Téléphone
Non applicable (L’examen a
Alexandre Al-haiek lieu sur Moodle)
N/A (voir directives)
Jour Date Durée Heures
Lundi 15 juin 2020 3h 8h30 – 11h30
Documentation Calculatrice
Aucune Aucune
Les cellulaires, agendas
Toute Toutes électroniques ou
téléavertisseurs sont
Non programmable interdits.
Voir directives particulières
(AEP)
Directives particulières
-Toute la documentation est permise (notes et manuel du cours).
-Exprimez vos réponses en unités SI (kg, m, s, joule, etc.)
-Une partie de la pondération sera accordée à la clarté ET concision de votre démarche/réponse.
-Pour la démarche des questions chiffrées, utilisez uniquement l’espace disponible et ne perdez pas votre
temps dans l’affichage de votre équation. Allez-y directement avec une forme simpliste, ex. : des * pour les
multiplications, des / pour les divisions, des ^ pour les exposants, etc.
-Par souci d’équité, nous ne répondrons à aucune question lors de l’examen.
-Par contre, si vous avez des problèmes techniques avec Moodle, prenez une photo et envoyez la photo à
mon courriel ([Link]-haiek@[Link]). Je vous contacterai dans les minutes qui suivront.
Cet examen contient x9 questions sur un total de xx pages
Important

(excluant cette page)

La pondération de cet examen est de 35 %

Vous devez répondre sur : le questionnaire Moodle les deux

Vous devez remettre le questionnaire : oui non N/A

 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

RÉPONDEZ AUX 9 QUESTIONS SUIVANTES DANS LA PLATEFORME TEST DE

MOODLE.

DATE DE L’EXAMEN : LUNDI 15 JUIN 2020

DURÉE DE L’EXAMEN : 3 h (de 8 h 30 à 11 h 30)

Total : /27
Cet examen est composé de trois sections :

SECTION 0 : INFORMATIONS ET DÉCLARATION SUR L’HONNEUR (1 QUESTION)

Temps suggéré : 5 minutes

SECTION 1 : RÉSOLUTION DE PROBLÈMES (3 QUESTIONS)

Temps suggéré : 145 minutes

Total : /16
6,5
SECTION 2 : VRAI OU FAUX ET RÉPONSES COURTES (4 QUESTIONS)

Temps suggéré : 30 minutes

Total : /11
13,5
SECTION BONUS (1 QUESTION)

Total : /3
13,5

Page 2 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Section 0 : Informations et déclaration sur l’honneur

Question 0 (0 point)

Je déclare sur mon honneur que :

• ce travail est un travail original de ma part.

• j’ai répondu aux différentes questions par moi-même, sans l’aide d’une
autre personne par tout moyen de communication que ce soit.
• j’ai respecté toutes les conditions énoncées dans les consignes.
• personne n’a effectué ni rédigé ce travail ou une partie de ce travail à
ma place, gratuitement ou contre rémunération.
• je n’ai ni copié ni sauvegardé aucune partie du questionnaire ou du
matériel transmis pour cet examen.

i. Oui
ii. Non

Page 3 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Section 1 : Résolution de problèmes

Question 1 (6 points)
Soit la réaction réversible, isotherme et isobare en phase liquide suivante :

k
A↔ B

La réaction est élémentaire, la constante de vitesse (k) est une fonction de la

température selon la loi d’Arrhenius avec un facteur pré-exponentiel A =
46.5 𝑚𝑖𝑛−1 et une énergie d’activation de 14 800 J/mol. Considérez les conditions
suivantes :

➢ V = 8 m3
➢ T = 475 °C

a) Quelle constante d’équilibre permettrait d’obtenir une conversion d’au

moins 70% si la réaction est produite dans un réacteur PFR? Considérez
un débit volumétrique de 380 L/s. (2 points)

Étape 1 : Établir le bilan molaire (équation de design)

Utilisons le bilan molaire général :

𝑽 𝒅𝑵𝒋
𝑭𝒋𝟎 − 𝑭𝒋 + ∫ 𝒓𝒋 𝒅𝑽 =
𝒅𝒕

Pour un PFR, il faut maintenant différencier cette équation par rapport

au volume (V) :
∆𝑽
∆𝑮𝒋 = ∫ 𝒓𝒋 𝒅𝑽 = 𝒓𝒋 ∆𝑽

Le bilan molaire général s’écrit donc (en régime permanent) :

𝑭𝒋 | − 𝑭𝒋 | + 𝒓𝒋 ∆𝑽 = 𝟎
𝑽 𝑽+∆𝑽

𝑭𝒋 | − 𝑭𝒋 |
𝑽+∆𝑽 𝑽
− 𝒓𝒋 = 𝟎
∆𝑽

𝑭𝒋 | − 𝑭𝒋 |
𝑽+∆𝑽 𝑽
= 𝒓𝒋
∆𝑽

Par définition de la limite :

Page 4 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝑭𝒋 | − 𝑭𝒋 | 𝒅𝑭𝒋
𝑽+∆𝑽 𝑽
𝐥𝐢𝐦 =
∆𝑽→ 𝟎 ∆𝑽 𝒅𝑽

Qui donne,

𝒅𝑭𝒋
= 𝒓𝒋
𝒅𝑽

Avec 𝑭𝒋 = 𝑭𝒋𝟎 (𝟏 − 𝑿)

𝒅𝑭𝒋 𝒅[𝑭𝒋𝟎 (𝟏 − 𝑿)]

= = 𝒓𝒋
𝒅𝑽 𝒅𝑽
𝒅𝑭𝒋𝟎 𝒅𝑭𝒋𝟎 𝑿
− = 𝒓𝒋
𝒅𝑽 𝒅𝑽

𝑭𝒋𝟎 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆

L’équation de design d’un PFR est alors :

𝒅𝑿 𝒓𝒋
= −
𝒅𝑽 𝑭𝒋𝟎

Étape 2 : Exprimer la vitesse de réaction (-rA)

Soient les réactions suivantes :

𝒌𝟏
- 𝑨 → 𝑩 (R1) ;
𝒌−𝟏
- 𝑩 → 𝑨 (R2).

Selon R1 : −𝒓𝑨 = 𝒌𝟏 𝑪𝑨
Selon R2 : 𝒓𝑨 = 𝒌−𝟏 𝑪𝑩

Finalement,
𝒓𝑨,𝒏𝒆𝒕 = 𝒓𝑨,𝒇𝒐𝒘 + 𝒓𝑨,𝒓𝒆𝒗 = −𝒌𝟏 𝑪𝑨 + 𝒌−𝟏 𝑪𝑩

−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 𝑪𝑨 − 𝒌−𝟏 𝑪𝑩

𝒌
Avec 𝑲𝑪 = 𝒌 𝟏
−𝟏

𝑪𝑩
−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 ∗ [𝑪𝑨 − ]
𝑲𝑪
Page 5 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 3 : Détailler la stœchiométrie : les concentrations CA et CB en

fonction de la conversion X

Le réacteur PFR est un réacteur en écoulement et la réaction se

produit en phase liquide.

- 𝑪𝑨 = 𝑪𝑨𝟎 ∗ (𝟏 − 𝑿), A étant le réactif limitant ;

𝒚
- 𝑪𝑩 = 𝑪𝑨𝟎 ∗ (𝜽𝑩 + 𝑿), avec 𝜽𝑩 = 𝒚𝑩𝟎 . Puisque 𝒚𝑩𝟎 = 𝟎 → 𝜽𝑩 = 𝟎,
𝑨𝟎
𝑪𝑩 = 𝑪𝑨𝟎 ∗ 𝑿

Finalement,
𝑪𝑨𝟎 ∗ 𝑿
−𝒓𝑨 = 𝒌𝟏 ∗ [𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − ]
𝑲𝑪

Étape 4 : Combiner l’expression de la vitesse réactionnelle avec la

stœchiométrie

À l’équilibre, on sait que −𝒓𝑨 = 𝟎. En utilisant l’expression de la

vitesse de réaction trouvée à l’étape 2 et la stœchiométrie de l’étape
3, on obtient donc :

𝑪𝑨𝟎 ∗ 𝑿
−𝒓𝑨 = 𝒌𝟏 ∗ [𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − ]= 𝟎
𝑲𝑪

Étape 5 : Résoudre l’équation pour déterminer la valeur de la

constante d’équilibre Kc

𝑪𝑨𝟎 ∗ 𝑿
[𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − ]= 𝟎
𝑲𝑪

𝑿
𝑪𝑨𝟎 ∗ [(𝟏 − 𝑿) − ]= 𝟎
𝑲𝑪

𝑿
(𝟏 − 𝑿) − =𝟎
𝑲𝑪

𝑿
(𝟏 − 𝑿) =
𝑲𝑪

𝑿
𝑲𝑪 =
(𝟏 − 𝑿)

Page 6 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝟎. 𝟕
𝑲𝑪 =
(𝟏 − 𝟎. 𝟕)

Réponse : 𝑲𝑪 ∈ [𝟐 ; 𝟑]

b) Considérons une constante d’équilibre égale à 28. Quel débit volumétrique

permettrait d’obtenir une conversion de 60% si la réaction est produite dans
un réacteur CSTR? (2 points)

Étape 1 : Établir le bilan molaire (équation de design)

La réaction est produite en phase liquide dans un CSTR, l’équation de

design associée est (mêmes étapes qu’étape 1 question a) mais en
posant qu’il n’y a pas de variation spatiale de la vitesse de réaction
(parfaitement mélangé)) :

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
−𝒓𝑨,𝒏𝒆𝒕

Étape 2 : Exprimer la vitesse de réaction (-rA)

Soient les réactions suivantes :

𝒌𝟏
- 𝑨 → 𝑩 (R1) ;
𝒌−𝟏
- 𝑩 → 𝑨 (R2).

Selon R1 : −𝒓𝑨 = 𝒌𝟏 𝑪𝑨
Selon R2 : 𝒓𝑨 = 𝒌−𝟏 𝑪𝑩

Finalement,
𝒓𝑨,𝒏𝒆𝒕 = 𝒓𝑨,𝒇𝒐𝒘 + 𝒓𝑨,𝒓𝒆𝒗 = −𝒌𝟏 𝑪𝑨 + 𝒌−𝟏 𝑪𝑩

−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 𝑪𝑨 − 𝒌−𝟏 𝑪𝑩

𝒌
Avec 𝑲𝑪 = 𝒌 𝟏
−𝟏

𝑪𝑩
−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 ∗ [𝑪𝑨 −]
𝑲𝑪
Étape 3 : Détailler la stœchiométrie : les concentrations CA et CB en
fonction de la conversion X

Page 7 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Le réacteur PFR est un réacteur en écoulement et la réaction se

produit en phase liquide.

- 𝑪𝑨 = 𝑪𝑨𝟎 ∗ (𝟏 − 𝑿), A étant le réactif limitant ;

𝒚
- 𝑪𝑩 = 𝑪𝑨𝟎 ∗ (𝜽𝑩 + 𝑿), avec 𝜽𝑩 = 𝒚𝑩𝟎 . Puisque 𝒚𝑩𝟎 = 𝟎 → 𝜽𝑩 = 𝟎,
𝑨𝟎
𝑪𝑩 = 𝑪𝑨𝟎 ∗ 𝑿

Finalement,
𝑪𝑨𝟎 ∗ 𝑿
−𝒓𝑨 = 𝒌𝟏 ∗ [𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − ]
𝑲𝑪

Étape 4 : Insérer l’expression de la vitesse réactionnelle avec la

stœchiométrie dans l’équation de design

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
−𝒓𝑨,𝒏𝒆𝒕

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
𝑪
𝒌𝟏 ∗ [𝑪𝑨 − 𝑲𝑩 ]
𝑪

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
𝑪𝑨𝟎 ∗ 𝑿
𝒌𝟏 ∗ [𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − 𝑲𝑪 ]

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
𝑿
𝒌𝟏 ∗ 𝑪𝑨𝟎 ∗ [𝟏 − 𝑿 − 𝑲 ]
𝑪

𝑭
On sait que : 𝑪𝑨𝟎 = 𝒗𝟎 , l’expression du volume du CSTR devient alors :
𝑨𝟎
𝒗𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
𝑿
𝒌𝟏 ∗ [𝟏 − 𝑿 − 𝑲 ]
𝑪

Isolons 𝒗𝟎 :

𝑽𝑪𝑺𝑻𝑹 ∗ 𝒌𝟏 𝑿
𝒗𝟎 = ∗ [𝟏 − 𝑿 − ]
𝑿 𝑲𝑪

Page 8 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 5 : Résoudre l’équation pour déterminer la valeur du débit

volumétrique (𝒗𝟎 )

Avec :
- 𝑽𝑪𝑺𝑻𝑹 = 𝟖 𝒎𝟑 ;
- 𝒌𝟏 = 𝟒. 𝟑𝟏 𝒎𝒊𝒏−𝟏 ;
- 𝑲𝑪 = 𝟐𝟖 ;
- 𝑿 = 𝟎. 𝟔

𝟖 ∗ 𝟒. 𝟑𝟏 𝟎. 𝟔
𝒗𝟎 = ∗ [𝟏 − 𝟎. 𝟔 − ]
𝟎. 𝟔 𝟐𝟖

𝒎𝟑
𝒗𝟎 ∈ [𝟐𝟎 ; 𝟐𝟒]
𝒎𝒊𝒏

c) En gardant les conditions constantes mise à part la température, quelle doit

être la température d’opération d’un CSTR qui permet d’obtenir une
conversion identique à celle du PFR de la sous-question a)? Considérez le
débit volumétrique de la sous-question a). (2 points)

Étape 1 : Établir le bilan molaire (équation de design)

La réaction est produite en phase liquide dans un CSTR, l’équation de

design associée est (mêmes étapes qu’étape 1 question a) mais en
posant qu’il n’y a pas de variation spatiale de la vitesse de réaction
(parfaitement mélangé)) :

𝑭𝑨𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
−𝒓𝑨,𝒏𝒆𝒕

Étape 2 : Exprimer la vitesse de réaction (-rA)

Soient les réactions suivantes :

𝒌𝟏
- 𝑨 → 𝑩 (R1) ;
𝒌−𝟏
- 𝑩 → 𝑨 (R2).

Selon R1 : −𝒓𝑨 = 𝒌𝟏 (𝑻𝟏 )𝑪𝑨

Selon R2 : 𝒓𝑨 = 𝒌−𝟏 (𝑻𝟏 )𝑪𝑩

Finalement,
𝒓𝑨,𝒏𝒆𝒕 = 𝒓𝑨,𝒇𝒐𝒘 + 𝒓𝑨,𝒓𝒆𝒗 = −𝒌𝟏 (𝑻𝟏 )𝑪𝑨 + 𝒌−𝟏 (𝑻𝟏 )𝑪𝑩

−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 (𝑻𝟏 )𝑪𝑨 − 𝒌−𝟏 (𝑻𝟏 )𝑪𝑩

Page 9 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝒌 (𝑻 )
Avec 𝑲𝑪 = 𝒌 𝟏 (𝑻𝟏 )
−𝟏 𝟏

𝑪𝑩
−𝒓𝑨,𝒏𝒆𝒕 = 𝒌𝟏 (𝑻𝟏 ) ∗ [𝑪𝑨 − ]
𝑲𝑪

Étape 3 : Détailler la stœchiométrie : les concentrations CA et CB en

fonction de la conversion X

Le réacteur PFR est un réacteur en écoulement et la réaction se

produit en phase liquide.

- 𝑪𝑨 = 𝑪𝑨𝟎 ∗ (𝟏 − 𝑿), A étant le réactif limitant ;

𝒚
- 𝑪𝑩 = 𝑪𝑨𝟎 ∗ (𝜽𝑩 + 𝑿), avec 𝜽𝑩 = 𝒚𝑩𝟎 . Puisque 𝒚𝑩𝟎 = 𝟎 → 𝜽𝑩 = 𝟎,
𝑨𝟎
𝑪𝑩 = 𝑪𝑨𝟎 ∗ 𝑿

Finalement,
𝑪𝑨𝟎 ∗ 𝑿
−𝒓𝑨 = 𝒌𝟏 (𝑻𝟏 ) ∗ [𝑪𝑨𝟎 ∗ (𝟏 − 𝑿) − ]
𝑲𝑪

Étape 4 : Insérer l’expression de la vitesse réactionnelle avec la

stœchiométrie dans l’équation de design

𝒗𝟎 ∗ 𝑿
𝑽𝑪𝑺𝑻𝑹 =
𝑿
𝒌𝟏 (𝑻𝟏 ) ∗ [𝟏 − 𝑿 − 𝑲 ]
𝑪
Isolons 𝒌𝟏 :
𝒗𝟎 ∗ 𝑿
𝒌𝟏 (𝑻𝟏 ) =
𝑿
𝑽𝑪𝑺𝑻𝑹 ∗ [𝟏 − 𝑿 − 𝑲 ]
𝑪

Étape 5 : Résoudre l’équation pour déterminer la valeur de la

constante réactionnelle 𝒌𝟏 (𝑻𝟏 )

On désire atteindre une conversion de 70% dans un CSTR. En utilisant

l’équation de design du CSTR, on a alors que :

𝟐𝟐. 𝟖 ∗ 𝟎. 𝟕
𝒌𝟏 (𝑻𝟏 ) =
𝟎. 𝟕
𝟖 ∗ [𝟏 − 𝟎. 𝟕 − 𝟐𝟖 ]

𝒌𝟏 (𝑻𝟏 ) ∈ [𝟕 ; 𝟕. 𝟓] 𝒎𝒊𝒏−𝟏
Page 10 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 6 : Trouver la température 𝑻𝟏

𝑬
D’après l’équation d’Arrhénius : 𝒌𝟏 (𝑻𝟏 ) = 𝑨 ∗ 𝐞𝐱𝐩 (− 𝑹𝑻 ), avec 𝑻𝟏 la
𝟏
température recherchée.
Les données cinétiques sont fournies à 𝑻𝟐 = 𝟒𝟕𝟓°𝑪 (c.f énoncé). On
peut ainsi établir le système d’équation suivant :

𝑬
𝒌𝟏 (𝑻𝟏 ) = 𝑨 ∗ 𝐞𝐱𝐩 (− )
𝑹𝑻𝟏
𝑬
𝒌𝟐 (𝑻𝟐 ) = 𝑨 ∗ 𝐞𝐱𝐩 (− )
{ 𝑹𝑻𝟐

𝑬
𝒌𝟏 (𝑻𝟏 ) 𝐞𝐱𝐩 (− 𝑹𝑻𝟏 ) 𝑬 𝑬
= = 𝐞𝐱𝐩 (− + )
𝒌𝟐 (𝑻𝟐 ) 𝐞𝐱𝐩 (− 𝑬 ) 𝑹𝑻𝟏 𝑹𝑻𝟐
𝑹𝑻𝟐

𝒌𝟏 (𝑻𝟏 ) 𝑬 𝑬 𝑬 𝟏 𝟏
𝐥𝐧 ( )=− + = ∗( − )
𝒌𝟐 (𝑻𝟐 ) 𝑹𝑻𝟏 𝑹𝑻𝟐 𝑹 𝑻𝟐 𝑻𝟏

𝒌𝟏 (𝑻𝟏 ) 𝑹 𝟏 𝟏
𝐥𝐧 ( )∗ = −
𝒌𝟐 (𝑻𝟐 ) 𝑬 𝑻𝟐 𝑻𝟏
𝟏 𝟏 𝒌𝟏 𝟏 ) 𝑹
(𝑻
= − 𝐥𝐧 ( )∗
𝑻𝟏 𝑻𝟐 𝒌𝟐 (𝑻𝟐 ) 𝑬

−𝟏
𝟏 𝒌𝟏 (𝑻𝟏 ) 𝑹
𝑻𝟏 = ( − 𝐥𝐧 ( )∗ )
𝑻𝟐 𝒌𝟐 (𝑻𝟐 ) 𝑬

𝟏 𝟕. 𝟐𝟓 𝟖. 𝟑𝟏𝟒 −𝟏
𝑻𝟏 = ( − 𝐥𝐧 ( )∗ )
(𝟒𝟕𝟓 + 𝟐𝟕𝟑. 𝟏𝟓 𝟒. 𝟑𝟏 𝟏𝟒 𝟖𝟎𝟎

𝑻𝟏 ∈ [𝟔𝟔𝟓 ; 𝟔𝟗𝟎]°𝑪

Page 11 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Question 2 (6 points)
Soit la réaction isotherme en phase gazeuse suivante :

A + 2B → C + 2D

Le produit d’intérêt est l’espèce C. La réaction est produite dans un réacteur à lit
fixe. La réaction est d’ordre zéro en A et de premier ordre en B. Le débit
m3
volumétrique alimenté au réacteur est de 30 min, la pression initiale est de 10 bars,
la température est de 523.15 K et le paramètre de perte de charge, α, est de
0.01 kg −1 . L’alimentation est équimolaire en A et en B, et la constante de vitesse
de réaction est de 0.03 s −1.

a) Quelle serait la conversion si une masse de catalyseur de 80 kg est utilisée?

(2 points)

Étape 1 : Établir le bilan molaire (équation de design)

La réaction est réalisée en phase gazeuse dans un réacteur PBR.

L’équation de design à utiliser est (mêmes étapes qu’étape 1 question
1.a) mais en différenciant sur la masse de catalyseur W) :

𝒅𝑿 𝒓′𝑨
= −
𝒅𝑾 𝑭𝑨𝟎

Étape 2 : Exprimer la vitesse de réaction (-rA)

D’après les informations de l’énoncé :

−𝒓′𝑨 = 𝒌𝑪𝟎𝑨 𝑪𝑩 = 𝒌𝑪𝑩

Étape 3 : Détailler la stœchiométrie : la concentration CB en fonction

de la conversion X

A est le réactif limitant, la concentration CB s’exprime comme

suit (phase gazeuse, système en écoulement) :

𝒃
(𝜽𝑩 − 𝒂 𝑿) 𝑷 𝑻𝟎
𝑪𝑩 = 𝑪𝑨𝟎 ∗ ∗ ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻
𝒚𝑩𝟎 𝒃
Avec 𝜽𝑩 = = 𝟏 (alimentation équimolaire), 𝒂 = 𝟐, et en sachant que
𝒚𝑨𝟎
la réaction est isotherme, l’expression devient :

Page 12 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

(𝟏 − 𝟐𝑿) 𝑷
𝑪𝑩 = 𝑪𝑨𝟎 ∗ ∗( )
(𝟏 + 𝜺𝑿) 𝑷𝟎

Si 𝜺 = 𝜹 ∗ 𝒚𝑨𝟎 = (𝟐 + 𝟏 − 𝟐 − 𝟏) ∗ 𝟎. 𝟓 = 𝟎, CB s’écrit :

𝑷
𝑪𝑩 = 𝑪𝑨𝟎 ∗ (𝟏 − 𝟐𝑿) ∗ ( )
𝑷𝟎

Étape 4 : Trouver l’expression de 𝑷/𝑷𝟎

On sait que (avec 𝒑 = 𝑷/𝑷𝟎 ) :

𝒅𝒑 𝜶 𝑻 𝑭𝑻
= − ( )( )
𝒅𝑾 𝟐𝒑 𝑻𝟎 𝑭𝑻𝟎
𝑭
Avec 𝑭 𝑻 = 𝟏 + 𝜺𝑿 = 𝟏 (𝒄𝒂𝒓 𝜺 = 𝟎) et réaction isotherme :
𝑻𝟎

𝒅𝒑 𝜶
= −
𝒅𝑾 𝟐𝒑

𝟐𝒑𝒅𝒑 = −𝜶𝒅𝑾

Sachant qu’à W=0, p=1 car P=P0 :

𝒑 𝑾
∫ 𝟐𝒑𝒅𝒑 = ∫ −𝜶𝒅𝑾
𝟏 𝟎

[𝒑𝟐 ]𝒑𝟏 = [−𝜶𝑾]𝑾

𝟎

On obtient alors :

𝒑𝟐 − 𝟏 = −𝜶𝑾
𝒑𝟐 = 𝟏 − 𝜶𝑾
𝟏
𝒑 = (𝟏 − 𝜶𝑾)𝟐

Page 13 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 5 : Combiner l’expression de la vitesse réactionnelle, la

stœchiométrie et l’expression de la perte de charge dans l’équation
de design

𝒅𝑿 𝒓′ 𝑨 𝒌𝑪𝑩 𝑪𝑨𝟎 𝑷
=− = = ∗ 𝒌 ∗ (𝟏 − 𝟐𝑿) ∗ ( )
𝒅𝑾 𝑭𝑨𝟎 𝑭𝑨𝟎 𝑭𝑨𝟎 𝑷𝟎

𝑪 𝟏
Avec 𝑭𝑨𝟎 = 𝒗 et 𝒑 = 𝑷/𝑷𝟎 :
𝑨𝟎 𝟎

𝒅𝑿 𝟏
= ∗ 𝒌 ∗ (𝟏 − 𝟐𝑿) ∗ 𝒑
𝒅𝑾 𝒗𝟎
𝟏
D’après l’étape 3, 𝒑 = (𝟏 − 𝜶𝑾)𝟐 :

𝟏
𝒅𝑿 𝟏
= 𝒗 ∗ 𝒌 ∗ (𝟏 − 𝟐𝑿) ∗ (𝟏 − 𝜶𝑾)𝟐
𝒅𝑾 𝟎

Étape 6 : Résoudre l’équation pour déterminer la valeur de la

conversion (X)

𝒅𝑿 𝒌
= ∗ (𝟏 − 𝜶𝑾)𝟏/𝟐 𝒅𝑾
(𝟏 − 𝟐𝑿) 𝒗𝟎

𝑿 𝑾
𝒅𝑿 𝒌
∫ = ∗ ∫ (𝟏 − 𝜶𝑾)𝟏/𝟐 𝒅𝑾
𝟎 (𝟏 − 𝟐𝑿) 𝒗𝟎 𝟎

𝑿 𝟑 𝑾
𝟏 𝒌 𝟐
[− ∗ 𝐥𝐧(𝟏 − 𝟐𝑿)] = ∗ [ (𝟏 − (𝟏 − 𝜶𝑾)𝟐 )]
𝟐 𝟎 𝒗𝟎 𝟑𝜶 𝟎

𝟒∗𝒌 𝟑
𝐥𝐧(𝟏 − 𝟐𝑿) = − (𝟏 − (𝟏 − 𝜶𝑾)𝟐 )
𝟑 ∗ 𝒗𝟎 ∗ 𝜶

𝟒∗𝒌 𝟑
𝟏 − 𝟐𝑿 = 𝐞𝐱𝐩 (− (𝟏 − (𝟏 − 𝜶𝑾)𝟐 ))
𝟑 ∗ 𝒗𝟎 ∗ 𝜶

Page 14 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝟒∗𝒌 𝟑
𝟐𝑿 = 𝟏 − 𝐞𝐱𝐩 (− (𝟏 − (𝟏 − 𝜶𝑾)𝟐 ))
𝟑 ∗ 𝒗𝟎 ∗ 𝜶

𝟒∗𝒌 𝟑
𝟏 − 𝐞𝐱𝐩 (− 𝟑 ∗ 𝒗 ∗ 𝜶 (𝟏 − (𝟏 − 𝜶𝑾)𝟐 ))
𝟎
𝑿=
𝟐

𝟒 ∗ 𝟎. 𝟎𝟑 𝟑
𝟏 − 𝐞𝐱𝐩 (− (𝟏 − (𝟏 − 𝟎. 𝟎𝟏 ∗ 𝟖𝟎)𝟐 ))
𝟑 ∗ 𝟎. 𝟓 ∗ 𝟎. 𝟎𝟏
𝑿=
𝟐

𝑿 ∈ [𝟎. 𝟒𝟓 ; 𝟎. 𝟓𝟓]

b) Quelle est alors la pression à la sortie du réacteur? (0.5 point)

On sait que (voir étape 4, question a)) :

𝟏
𝒑 = (𝟏 − 𝜶𝑾)𝟐
𝟏
𝒑 = (𝟏 − 𝟎. 𝟎𝟏 ∗ 𝟖𝟎)𝟐
𝒑 = 𝟎. 𝟒𝟓
𝑷
Puisque 𝒑 = 𝑷 :
𝟎

𝑷 = 𝒑 ∗ 𝑷𝟎
𝑷 = 𝟎. 𝟒𝟓 ∗ 𝟏
𝑷 ∈ [𝟒. 𝟓 ; 𝟓] 𝒃𝒂𝒓𝒔

c) Calculez maintenant le ratio de masse de catalyseur (𝑊1 /𝑊2) entre les 10

premiers % de conversion (de 0 à 10 % où 𝑊1 = 𝑊10% − 𝑊0% ) et les 10
derniers % de conversion (en considérant la conversion trouvée en a)
comme étant la borne supérieure des 10 derniers % de conversion où 𝑊2 =
𝑊𝑠𝑜𝑢𝑠 𝑞𝑢𝑒𝑠𝑡𝑖𝑜𝑛 𝑎) − 𝑊𝑋𝑎) −10%). (1 point)

D’après l’étape 6 de la question a) :

𝑿 𝟑 𝑾
𝟏 𝒌 𝟐
[− ∗ 𝐥𝐧(𝟏 − 𝟐𝑿)] = ∗ [ (𝟏 − (𝟏 − 𝜶𝑾)𝟐 )]
𝟐 𝟎 𝒗𝟎 𝟑𝜶 𝟎

−𝒌 𝟒 𝟑
𝐥𝐧(𝟏 − 𝟐𝑿) = ∗ (𝟏 − (𝟏 − 𝜶𝑾)𝟐 )
𝒗𝟎 𝟑𝜶

Page 15 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

−𝟑𝜶 ∗ 𝒗𝟎 ∗ 𝐥𝐧(𝟏 − 𝟐𝑿) 𝟑

= 𝟏 − (𝟏 − 𝜶𝑾)𝟐
𝟒∗𝒌
𝟑 −𝟑𝜶 ∗ 𝒗𝟎 ∗ 𝐥𝐧(𝟏 − 𝟐𝑿)
(𝟏 − 𝜶𝑾)𝟐 = 𝟏 −
𝟒∗𝒌
𝟐
−𝟑𝜶∗𝒗𝟎 ∗𝐥𝐧(𝟏−𝟐𝑿) 𝟑
𝟏 − 𝜶𝑾 = (𝟏 − )
𝟒∗𝒌

𝟐
−𝟑𝜶 ∗ 𝒗𝟎 ∗ 𝐥𝐧(𝟏 − 𝟐𝑿) 𝟑
𝜶𝑾 = 𝟏 − (𝟏 − )
𝟒∗𝒌

𝟐
−𝟑𝜶 ∗ 𝒗𝟎 ∗ 𝐥𝐧(𝟏 − 𝟐𝑿) 𝟑
𝟏 − (𝟏 − )
𝟒∗𝒌
𝑾=
𝜶

Ce qui donne :

X= 0 W=0
X = 0.1 W = 1.87 kg
X = 0.4 W = 13.9 kg
𝑾𝟏
𝑾𝟐
0.028

𝑾𝟏
Réponse : 𝑾𝟐 ∈ [𝟎. 𝟎𝟐 ; 𝟎. 𝟎𝟑]

d) Commentez le résultat en c) d’un point de vue conception et économique.

(2 points)

Ce faible ratio veut dire que la vitesse de réaction ralentie et que les
derniers % sont très difficiles à aller chercher.

Économique : Dans le cas où le catalyseur est très dispendieux, il

n’est peut-être pas rentable d’aller chercher les derniers %, la
purification du produit est peut-être la meilleure option.
Conception : Plus de catalyseur = plus de perte de charge. Moins la
perte de charge est élevée, moins le compresseur nécessaire sera
dispendieux, moins de matériau pour le réacteur, etc.

Page 16 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

e) Proposez une modification au procédé pour atteindre la même conversion,

mais avec moins de catalyseur. (0.5 point)

Par exemple :
- Diminuer le débit volumétrique
- Jouer avec la taille de catalyseur (diminuer 𝜶), mais
ATTENTION, d’autres problèmes peuvent survenir.
- Etc.

Page 17 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Question 3 (4 points)
Soit la réaction isotherme en phase gazeuse suivante :

A+B→C

La réaction est produite dans un réacteur de type PBR avec une alimentation en
continu à une température de 28°C. L’alimentation est équimolaire en A et B avec
une pression de 125 psia à la sortie du réacteur. La réaction est du deuxième ordre
en A et d’ordre zéro en B. Une conversion de 50% est atteinte à l’aide de 110 kg
de catalyseur et un débit volumétrique de 1.67 L/s. Sachant que le paramètre de
perte de charge, α, est de 0.0095 kg −1 et que l’énergie d’activation est de 43
kJ/mol, calculez la constante de réaction à 450 K.

Étape 1 : Établir le bilan molaire (équation de design)

La réaction est réalisée en phase gazeuse dans un réacteur PBR.

L’équation de design à utiliser est (mêmes étapes qu’étape 1 question a)
mais en différenciant par rapport à la masse de catalyseur W) :

𝒅𝑿 𝒓′𝑨
= −
𝒅𝑾 𝑭𝑨𝟎

Étape 2 : Établir la vitesse de réaction (−𝒓𝑨)

D’après l’énoncé, on sait que −𝒓′𝑨 = 𝒌 ∗ 𝑪𝟐𝑨 ∗ 𝑪𝟎𝑩 = 𝒌 ∗ 𝑪𝟐𝑨

Étape 3 : Détailler la stœchiométrie : la concentration CA en fonction de

la conversion X

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷 𝑻𝟎
𝑪𝑨 = ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻

En sachant que la réaction est isotherme, l’expression devient :

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷
𝑪𝑨 = ( )
(𝟏 + 𝜺𝑿) 𝑷𝟎

Si 𝜺 = 𝜹 ∗ 𝒚𝑨𝟎 = (𝟏 + 𝟏 − 𝟏) ∗ 𝟎. 𝟓 = −𝟎. 𝟓, (réaction équimolaire)

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷
𝑪𝑨 = ( )
(𝟏 − 𝟎. 𝟓𝑿) 𝑷𝟎

Page 18 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 4 : Trouver l’expression de 𝑷/𝑷𝟎

On sait que (avec 𝒑 = 𝑷/𝑷𝟎 ) :

𝒅𝒑 𝜶 𝑻 𝑭𝑻
= − ( )( )
𝒅𝑾 𝟐𝒑 𝑻𝟎 𝑭𝑻𝟎
𝑭𝑻
Avec = 𝟏 + 𝜺𝑿 = 𝟏 − 𝟎. 𝟓𝑿 (𝒄𝒂𝒓 𝜺 = −𝟎. 𝟓) et réaction isotherme :
𝑭𝑻𝟎

𝒅𝒑 𝜶
= − (𝟏 − 𝟎. 𝟓𝑿)
𝒅𝑾 𝟐𝒑

𝟐𝒑𝒅𝒑 = −𝜶(𝟏 − 𝟎. 𝟓𝑿)𝒅𝑾

Sachant qu’à W=0, p=1 car P=P0 :

𝒑 𝑾
∫ 𝟐𝒑𝒅𝒑 = (𝟏 − 𝟎. 𝟓𝑿) ∗ ∫ −𝜶𝒅𝑾
𝟏 𝟎

[𝒑𝟐 ]𝒑𝟏 = (𝟏 − 𝟎. 𝟓𝑿) ∗ [−𝜶𝑾]𝑾

𝟎

On obtient alors :

𝒑𝟐 − 𝟏 = − (𝟏 − 𝟎. 𝟓𝑿)𝜶𝑾
𝒑𝟐 = 𝟏 − (𝟏 − 𝟎. 𝟓𝑿)𝜶𝑾
𝟏
𝒑 = (𝟏 − (𝟏 − 𝟎. 𝟓𝑿)𝜶𝑾)𝟐

D’après l’énoncé, pour atteindre X=0.5, une masse W = 110 kg de

catalyseur est nécessaire. On en déduit donc la valeur de p :
𝟏
𝒑 = (𝟏 − (𝟏 − 𝟎. 𝟓 ∗ 𝟎. 𝟓) ∗ 𝟎. 𝟎𝟎𝟗𝟓 ∗ 𝟏𝟏𝟎)𝟐

𝒑 =∈ [𝟎. 𝟒𝟓 ; 𝟎. 𝟓]

Page 19 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 5 : Combiner l’expression de la vitesse réactionnelle, la

stœchiométrie et l’expression de la perte de charge dans l’équation
de design

𝟐
𝒅𝑿 𝒓′ 𝑨 𝒌𝑪𝟐𝑨 𝒌 𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷
=− = = ∗( ( ))
𝒅𝑾 𝑭𝑨𝟎 𝑭𝑨𝟎 𝑭𝑨𝟎 (𝟏 − 𝟎. 𝟓𝑿) 𝑷𝟎

𝑪 𝟏
Avec 𝑭𝑨𝟎 = 𝒗 et 𝒑 = 𝑷/𝑷𝟎 :
𝑨𝟎 𝟎

𝒅𝑿 𝒌 ∗ 𝑪𝑨𝟎 (𝟏 − 𝑿)𝟐
= ∗ 𝒑𝟐
𝒅𝑾 𝒗𝟎 (𝟏 − 𝟎. 𝟓𝑿)𝟐

Étape 6 : Résoudre l’équation pour déterminer la valeur de la

température

(𝟏 − 𝟎. 𝟓𝑿)𝟐 𝒌 ∗ 𝑪𝑨𝟎 𝟐
𝟐
𝒅𝑿 = 𝒑 𝒅𝑾
(𝟏 − 𝑿) 𝒗𝟎

𝑿 (𝟏
− 𝟎. 𝟓𝑿)𝟐 𝒌 ∗ 𝑪𝑨𝟎 𝟐 𝑾
∫ 𝒅𝑿 = 𝒑 ∫ 𝒅𝑾
𝟎 (𝟏 − 𝑿)𝟐 𝒗𝟎 𝟎

𝑿
(𝟏 − 𝟎. 𝟓)𝟐 ∗ 𝑿 𝒌 ∗ 𝑪𝑨𝟎 𝟐
[𝟐 ∗ (−𝟎. 𝟓) ∗ (𝟏 − 𝟎. 𝟓)𝒍𝒏(𝟏 − 𝑿) + (−𝟎. 𝟓)𝟐 ∗ 𝑿 + ] = 𝒑 ∗𝑾
𝟏−𝑿 𝟎
𝒗𝟎

𝒗𝟎 𝟐
(𝟏 − 𝟎. 𝟓)𝟐 ∗ 𝑿
𝒌= [𝟐 ∗ ( −𝟎. 𝟓) ∗ ( 𝟏 − 𝟎. 𝟓) 𝒍𝒏 ( 𝟏 − 𝑿 ) + ( −𝟎. 𝟓) ∗ 𝑿 + ]
𝑪𝑨𝟎 ∗ 𝒑𝟐 ∗ 𝑾 𝟏−𝑿

Calculer la concentration CA0

𝑷𝟎 𝑷
𝑪𝑨𝟎 = 𝒚𝑨𝟎 ∗ = 𝟎. 𝟓 ∗
𝑹𝑻𝟎 𝒑 ∗ 𝑹 ∗ 𝑻𝟎

𝑷𝟎 𝟏𝟐𝟓 ∗ 𝟔 𝟖𝟗𝟓
𝑪𝑨𝟎 = 𝒚𝑨𝟎 ∗ = 𝟎. 𝟓 ∗
𝑹𝑻𝟎 𝟎. 𝟒𝟕 ∗ 𝟖. 𝟑𝟏𝟒 ∗ (𝟐𝟖 + 𝟐𝟕𝟑. 𝟏𝟓)

Page 20 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝒎𝒐𝒍
𝑪𝑨𝟎 ∈ [𝟑𝟒𝟎 ; 𝟑𝟖𝟓]
𝒎𝟑

Calculer la valeur de 𝒌(𝑻𝟏 = 𝟐𝟖°𝑪 = 𝟑𝟎𝟏°𝑲)

𝟎. 𝟎𝟎𝟏𝟔𝟕 𝟐
(𝟏 − 𝟎. 𝟓)𝟐 ∗ 𝟎. 𝟓
𝒌= [ 𝟐 ∗ ( −𝟎. 𝟓 ) ∗ (𝟏 − 𝟎. 𝟓 ) 𝒍𝒏 (𝟏 − 𝟎. 𝟓 ) + ( −𝟎. 𝟓 ) ∗ 𝟎. 𝟓 + ]
𝟑𝟕𝟎 ∗ 𝟎. 𝟒𝟕𝟐 ∗ 𝟏𝟏𝟎 𝟏 − 𝟎. 𝟓

𝒌(𝑻𝟏 = 𝟐𝟖°𝑪 = 𝟑𝟎𝟏°𝑲) ∈ [ 𝟏 ∗ 𝟏𝟎−𝟕 ; 𝟐 ∗ 𝟏𝟎−𝟕 ] 𝒎𝒊𝒏−𝟏

Calculer la valeur de 𝒌(𝑻𝟐 = 𝟒𝟓𝟎°𝑲)

En utilisant l’équation d’Arrhénius pour T1 et T2, on obtient le système

d’équation suivant :

𝑬
𝒌𝟏 (𝑻𝟏 ) = 𝑨 ∗ 𝐞𝐱𝐩 (− )
𝑹𝑻𝟏
𝑬
𝒌𝟐 (𝑻𝟐 ) = 𝑨 ∗ 𝐞𝐱𝐩 (− )
{ 𝑹𝑻𝟐

𝒌𝟏 (𝑻𝟏 ) 𝑬 𝑬
= 𝐞𝐱 𝐩 (− + )
𝒌𝟐 (𝑻𝟐 ) 𝑹𝑻𝟏 𝑹𝑻𝟐

𝒌𝟏 (𝑻𝟏 ) 𝑬 𝟏 𝟏
= 𝐞𝐱 𝐩 ( ( − ))
𝒌𝟐 (𝑻𝟐 ) 𝑹 𝑻𝟐 𝑻𝟏

𝒌𝟏 (𝑻𝟏 )
𝒌𝟐 (𝑻𝟐 ) =
𝑬 𝟏 𝟏
𝐞𝐱 𝐩 (𝑹 (𝑻 − 𝑻 ))
𝟐 𝟏

𝟏. 𝟒 ∗ 𝟏𝟎−𝟕
𝒌𝟐 (𝟒𝟓𝟎𝑲) =
𝟒𝟑 𝟎𝟎𝟎 𝟏 𝟏
𝐞𝐱 𝐩 ( 𝟖. 𝟑𝟏𝟒 ( − ))
𝟒𝟓𝟎 (𝟐𝟖 + 𝟐𝟕𝟑. 𝟏𝟓

𝒌𝟐 (𝟒𝟓𝟎°𝑲) ∈ [ 𝟐 ∗ 𝟏𝟎−𝟓 ; 𝟔 ∗ 𝟏𝟎−𝟓 ] 𝒎𝒊𝒏−𝟏

Page 21 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Section 2 : Vrai ou Faux et réponses courtes

Question 5 (5 points)
Pour les énoncés suivants, répondez vrai ou faux. Si vous répondez faux, corrigez
l’énoncé pour qu’il devienne vrai.

i. Dans un réacteur CSTR, plus le temps de réaction est grand, plus la

conversion est grande.
□ Vrai □ Faux

Pour un CSTR, c’est le volume qui augmente quand la conversion augmente.

ii. Plus l’énergie d’activation est grande, plus un changement de la

température de la réaction a d’effet sur la constante de réaction.
□ Vrai □ Faux

________________________________________________________________

iii. La vitesse de réaction est différente selon le type de réacteur.

□ Vrai □ Faux

La vitesse de réaction est indépendante du type de réacteur.

iv. La vitesse de réaction est une équation algébrique.

□ Vrai □ Faux

________________________________________________________________

v. Un graphique de Levenspiel permet d’estimer très rapidement l’aire

effective d’un réacteur.
□ Vrai □ Faux

Un graphique de Levenspiel permet d’estimer très rapidement le volume

d’un réacteur.

Page 22 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Question 6 (2 points)
Expliquez en quelques phrases quel serait l’effet d’omettre la perte de charge lors
de la conception d’un PBR.

Si la perte de charge est omise dans la conception d’un PBR, la

concentration 𝑪𝒋 , la conversion et la vitesse de réaction seront légèrement
surestimées. De plus, le débit volumétrique (𝒗𝟎 ) sera, quant à lui, sous-
estimé.

Question 7 (2 points)
Nommez une méthode permettant l’estimation du volume d’un PFR (ou PBR) à
l’aide de la courbe FA0/-rA en fonction de X sans avoir recours à des méthodes
d’intégrations numériques (ex. Simpson 1/3). Pourquoi cette méthode est-elle
valable?

Utiliser plusieurs petits CSTRs en série. Cette méthode n’est valable

uniquement si on utilise de petits volumes de CSTRs et en grand nombre.

D’autres réponses ont aussi été acceptées lorsque justifiées.

Question 8 (2 points)
Lors de votre stage chez A → B Inc., votre superviseur et vous réalisez que la
concentration en A à la sortie de réacteur (type PBR) augmente continuellement
par rapport à la valeur habituelle. Une vérification avec les opérateurs vous permet
de valider qu’il n’y a eu aucune variation du débit d’alimentation, de la température
et de la pression d’opération. De plus, il n’y a eu aucune fuite ou purge. Une
vérification avec les techniciens en contrôle qualité vous permet de valider qu’il n’y
a eu aucune variation de la qualité de la matière première reçue. Selon vous, quel
pourrait être le problème?

Il pourrait y avoir empoisonnement des catalyseurs.

D’autres réponses ont aussi été acceptées lorsque justifiées.

Page 23 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Section BONUS :

Question 4 (3 points)
Soit la réaction isotherme en phase gazeuse suivante :

A+B→C+D

La réaction est produite dans un réacteur de type PBR de 40 mm de longueur et

1 mm de diamètre. Le catalyseur est sous forme de particules sphériques de 40
𝜇𝑚 de diamètre. La pression partielle en A à l’entrée est de 240 kPa. L’alimentation
en réactifs est équimolaire. Le débit molaire en A est égal à 0.0003 mol/s. Le débit
volumétrique est de 0.00029 L/s. Considérez l’équation de vitesse suivante :

−𝑟𝐴′ = 𝑘𝐴 𝐶𝐴 𝐶𝐵

Déterminez la conversion à la sortie du réacteur. Le tableau suivant présente de

l’information additionnelle :

Paramètres Unité Valeur

W kg 3.5 × 10−6
T °C 120
α 𝑘𝑔−1 3.55 × 105
𝛷 % 40
k (à 120°C) m6/[Link] 0.004
𝜌 kg/m3 7
𝜇 Pa.s 1.94 × 10−5
AC m2 1.96 × 10−7
G kg/m2.s 10

Étape 1 : Établir le bilan molaire (équation de design)

La réaction est réalisée en phase gazeuse dans un réacteur PBR.

L’équation de design à utiliser est (mêmes étapes qu’étape 1 question a)
mais en différenciant par rapport à la masse de catalyseur W) :

𝒅𝑿 𝒓′𝑨
= −
𝒅𝑾 𝑭𝑨𝟎

Étape 2 : Établir la vitesse de réaction (−𝒓𝑨)

D’après l’énoncé, on sait que −𝒓′ 𝑨 = 𝒌 ∗ 𝑪𝟏𝑨 ∗ 𝑪𝟏𝑩 = 𝒌 ∗ 𝑪𝑨 ∗ 𝑪𝑩

Page 24 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 3 : Détailler la stœchiométrie : les concentrations CA et CB en

fonction de la conversion X

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷 𝑻𝟎
𝑪𝑨 = ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻

𝒃
𝑪𝑨𝟎 (𝜽𝑩 − 𝒂 𝑿) 𝑷 𝑻𝟎
𝑪𝑩 = ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻

𝒚𝑩𝟎 𝒃
Avec 𝜽𝑩 = = 𝟏 (alimentation équimolaire) et 𝒂 = 𝟏 :
𝒚𝑨𝟎
𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷 𝑻𝟎
𝑪𝑩 = ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻

En sachant que la réaction est isotherme, l’expression devient :

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷
𝑪𝑨 = ( )
(𝟏 + 𝜺𝑿) 𝑷𝟎

𝑪𝑨𝟎 (𝟏 − 𝑿) 𝑷 𝑻𝟎
𝑪𝑩 = ( )( )
(𝟏 + 𝜺𝑿) 𝑷𝟎 𝑻

Si 𝜺 = 𝜹 ∗ 𝒚𝑨𝟎 = (𝟏 + 𝟏 − 𝟏 − 𝟏) ∗ 𝟎. 𝟓 = 𝟎, (réaction équimolaire)

𝑷
𝑪𝑨 = 𝑪𝑨𝟎 (𝟏 − 𝑿) ∗ ( )
𝑷𝟎

𝑷
𝑪𝑩 = 𝑪𝑨𝟎 (𝟏 − 𝑿) ∗ ( )
𝑷𝟎

Page 25 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Étape 4 : Trouver l’expression pour la perte de charge

L’équation Ergun permet d’estimer la perte de charge dans un PBR :

𝒅𝑷 𝑮 𝟏 − 𝝓 𝟏𝟓𝟎(𝟏 − 𝝓)𝝁
= − ( 𝟑 )[ + 𝟏. 𝟕𝟓𝑮]
𝒅𝒛 𝝆𝒈𝒄 𝑫𝒑 𝝓 𝑫𝒑

Le réacteur PBR est opéré en régime permanent, de ce fait il y a

conservation du débit massique :

𝒎̇ 𝟎 = 𝒎̇

𝒎̇
Sachant que 𝝆 = :
𝒗

𝝆𝟎 𝒗𝟎 = 𝝆𝒗

On sait également, qu’en phase gazeuse :

𝑷 𝟎 𝑻 𝑭𝑻
𝒗 = 𝒗𝟎 ( )
𝑷 𝑻𝟎 𝑭𝑻𝟎

D’où :

𝝆𝟎 𝒗𝟎 𝑷 𝟎 𝑻 𝑭𝑻
𝝆= =𝝆 ( )
𝒗 𝑷 𝑻𝟎 𝑭𝑻𝟎

En introduisant l’expression de 𝝆 dans l’équation d’Ergun :

𝒅𝑷 𝑮 𝟏 − 𝝓 𝟏𝟓𝟎(𝟏 − 𝝓)𝝁 𝑷 𝟎 𝑻 𝑭𝑻
= − ( 𝟑 )[ + 𝟏. 𝟕𝟓𝑮] ( )
𝒅𝒛 𝝆𝟎 𝒈𝒄 𝑫𝒑 𝝓 𝑫𝒑 𝒑 𝑻𝟎 𝑭𝑻𝟎

Avec :

𝑮 𝟏 − 𝝓 𝟏𝟓𝟎(𝟏 − 𝝓)𝝁
𝜷𝟎 = ( 𝟑 )[ + 𝟏. 𝟕𝟓𝑮]
𝝆𝟎 𝒈𝒄 𝑫𝒑 𝝓 𝑫𝒑

𝒅𝑷 𝑷 𝟎 𝑻 𝑭𝑻
= −𝜷𝟎 ( )
𝒅𝒛 𝑷 𝑻𝟎 𝑭𝑻𝟎

Page 26 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Pour un réacteur PBR, c’est la masse de catalyseur qui est le paramètre

d’intérêt. La masse de catalyseur en fonction de la longueur du lit
s’exprime :

𝑾 = (𝟏 − 𝝓)𝑨𝒄 𝒛 ∗ 𝝆𝒄

Avec 𝝆𝒃 = 𝝆𝒄 (𝟏 − 𝝓)

On peut donc exprimer z en fonction de W selon :

𝑾
𝒛=
(𝟏 − 𝝓)𝑨𝒄 ∗ 𝝆𝒄

𝒅𝑷
En introduisant l’expression de z dans l’équation :
𝒅𝒛

𝒅𝑷 𝜷𝟎 𝑷 𝟎 𝑻 𝑭𝑻
= − ( )
𝒅𝑾 (𝟏 − 𝝓)𝑨𝒄 𝝆𝒄 𝑷 𝑻𝟎 𝑭𝑻𝟎

𝒅𝑷 𝜷𝟎 𝑷 𝟎 𝑻 𝑭𝑻
= − ( )
𝒅𝑾 (𝟏 − 𝝓)𝑨𝒄 𝝆𝒄 𝑷 𝑻𝟎 𝑭𝑻𝟎

𝒅𝑷 𝜷𝟎 𝑷 𝟎 𝑻 𝑭𝑻
= − ( )
𝒅𝑾 (𝟏 − 𝝓)𝑨𝒄 𝝆𝒄 𝑷 𝑻𝟎 𝑭𝑻𝟎

Cette expression se simplifie en :

𝒅𝑷 𝜶 𝑷𝟎 𝑻 𝑭𝑻
= − ( )
𝒅𝑾 𝟐 𝑷/𝑷𝟎 𝑻𝟎 𝑭𝑻𝟎
𝑷
En posant 𝒑 = 𝑷 ,
𝟎
𝒅𝒑 𝜶 𝑻 𝑭𝑻
= − ( )
𝒅𝑾 𝟐𝒑 𝑻𝟎 𝑭𝑻𝟎

Puisque la valeur du coefficient de perte de charge (𝜶) est fournie dans

𝒅𝑷
l’énoncé, l’utilisation de l’équation d’Ergun 𝒅𝒛 n’est pas utile.

𝑭
Avec 𝑭 𝑻 = 𝟏 + 𝜺𝑿 = 𝟏 (𝒄𝒂𝒓 𝜺 = 𝟎) et réaction isotherme :
𝑻𝟎

𝒅𝒑 𝜶
= −
𝒅𝑾 𝟐𝒑

𝟐𝒑𝒅𝒑 = −𝜶𝒅𝑾
Page 27 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

Sachant qu’à W=0, p=1 car P=P0 :

𝒑 𝑾
∫ 𝟐𝒑𝒅𝒑 = ∫ −𝜶𝒅𝑾
𝟏 𝟎

[𝒑𝟐 ]𝒑𝟏 = [−𝜶𝑾]𝑾

𝟎

On obtient alors :

𝒑𝟐 − 𝟏 = − 𝜶𝑾
𝒑𝟐 = 𝟏 − 𝜶𝑾
𝟏
𝒑 = (𝟏 − 𝜶𝑾)𝟐

Étape 5 : Combiner l’expression de la vitesse réactionnelle, la

stœchiométrie et l’expression de la perte de charge dans l’équation de
design

𝒅𝑿 𝒓′ 𝑨 𝒌𝑪𝑨 𝑪𝑩 𝑪𝟐𝑨𝟎 𝑷 𝟐
=− = =𝒌∗ (𝟏 − 𝑿)𝟐 ∗ ( )
𝒅𝑾 𝑭𝑨𝟎 𝑭𝑨𝟎 𝑭𝑨𝟎 𝑷𝟎
𝑪 𝟏
Avec 𝑭𝑨𝟎 = 𝒗 et 𝒑 = 𝑷/𝑷𝟎 :
𝑨𝟎 𝟎

𝒅𝑿 𝒌 ∗ 𝑭𝑨𝟎
= ∗ (𝟏 − 𝑿)𝟐 ∗ 𝒑𝟐
𝒅𝑾 𝒗𝟐𝟎

𝒅𝑿 𝒌 ∗ 𝑭𝑨𝟎
= 𝟐
∗ (𝟏 − 𝑿)𝟐 ∗ (𝟏 − 𝜶𝑾)
𝒅𝑾 𝒗𝟎

Étape 6 : Résoudre l’équation pour déterminer la valeur de la conversion

𝒅𝑿 𝒌 ∗ 𝑭𝑨𝟎
= (𝟏 − 𝜶𝑾)𝒅𝑾
(𝟏 − 𝑿)𝟐 𝒗𝟐𝟎

𝑿
𝒅𝑿 𝒌 ∗ 𝑭𝑨𝟎 𝑾
∫ = ∫ (𝟏 − 𝜶𝑾)𝒅𝑾
𝟎 (𝟏 − 𝑿)
𝟐
𝒗𝟐𝟎 𝟎

Page 28 de 29
 ÉCOLE POLYTECHNIQUE DE MONTRÉAL
DÉPARTEMENT DE GÉNIE CHIMIQUE
GCH3110 - CALCUL DES RÉACTEURS CHIMIQUES
CONTRÔLE PÉRIODIQUE
ÉTÉ 2020

𝑾
𝑿 𝑿 𝒌 ∗ 𝑭𝑨𝟎 𝜶 ∗ 𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
[ ] = 𝑾 − [ ]
𝟏−𝑿 𝟎 𝒗𝟐𝟎 𝒗𝟐𝟎 𝟐 𝟎

𝑿 𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
= [ 𝑾 − 𝜶 ∗ ]
𝟏−𝑿 𝒗𝟐𝟎 𝟐

𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐 𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
𝑿= [𝑾−𝜶∗ ]− [𝑾−𝜶∗ ]∗𝑿
𝒗𝟐𝟎 𝟐 𝒗𝟐𝟎 𝟐

𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐 𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
𝑿 (𝟏 + [𝑾−𝜶∗ ]) = [𝑾−𝜶∗ ]
𝒗𝟐𝟎 𝟐 𝒗𝟐𝟎 𝟐

𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
[ 𝑾 − 𝜶 ∗ ]
𝒗𝟐𝟎 𝟐
𝑿=
𝒌 ∗ 𝑭𝑨𝟎 𝑾𝟐
(𝟏 + [ 𝑾 − 𝜶 ∗ ])
𝒗𝟐𝟎 𝟐

−𝟏𝟐
𝟎. 𝟎𝟎𝟒 ∗ 𝟎. 𝟎𝟎𝟎𝟑 −𝟔 −𝟓 𝟑. 𝟓 ∗ 𝟏𝟎
[ 𝟑. 𝟓 ∗ 𝟏𝟎 − 𝟑. 𝟓𝟓 ∗ 𝟏𝟎 ∗ ]
𝟐. 𝟗 ∗ 𝟏𝟎−𝟏𝟒 𝟐
𝑿= −𝟏𝟐
𝟎. 𝟎𝟎𝟒 ∗ 𝟎. 𝟎𝟎𝟎𝟑 −𝟔 −𝟓 𝟑. 𝟓 ∗ 𝟏𝟎
(𝟏 + [ 𝟑. 𝟓 ∗ 𝟏𝟎 − 𝟑. 𝟓𝟓 ∗ 𝟏𝟎 ∗ ])
𝟐. 𝟗 ∗ 𝟏𝟎−𝟏𝟒 𝟐

𝑿 ∈ [ 𝟎. 𝟗 ; 𝟎. 𝟗𝟗]

Bon succès!

Chargé de cours
Alexandre Al-haiek

Page 29 de 29