Procédés industriels continus : réacteurs ouverts
Cours · chimie : thermodynamique et électrochimie (PC), chapitre 5 · CPGE PC (2e année) · le début du chapitre, la suite sur Adloun
Les transformations réalisées au laboratoire mettent en jeu de faibles quantités de matière et sont conduites en réacteur fermé (bécher, ballon) : on introduit les réactifs, on laisse la réaction évoluer, puis on récupère le mélange. À l'échelle industrielle, les quantités manipulées sont considérables et l'on privilégie un fonctionnement continu, en réacteur ouvert : les réactifs sont injectés en permanence, les produits soutirés sans interruption. Ce chapitre réinvestit les acquis de cinétique et de thermodynamique pour établir les bilans de matière et les bilans énergétiques des deux modèles idéaux de réacteurs continus : le réacteur parfaitement agité continu (RPAC) et le réacteur en écoulement piston. On en déduit le dimensionnement, les points de fonctionnement et leur stabilité, et les règles de sécurité thermique ; deux compléments abordent la sélectivité et le compromis entre cinétique et thermodynamique.
5.1 Du protocole de laboratoire au procédé industriel
Un procédé est un enchaînement d'opérations unitaires : chaque opération réalise une action physique ou chimique élémentaire (mélange, réaction, échange thermique, compression, séparation gaz-liquide, distillation, évaporation, filtration, …). Le schéma de procédé les représente par des blocs reliés par des courants (flèches) de matière et d'énergie.
Un procédé est dit discontinu (ou en batch) lorsque la matière est traitée par lots : on charge, on transforme, on décharge ; la composition du réacteur dépend du temps. Il est dit continu lorsque les courants d'entrée et de sortie circulent en permanence. En régime permanent (ou stationnaire), les grandeurs (débits, compositions, température) sont indépendantes du temps en tout point de l'installation.
Pour un courant traversant une section, on définit :
- le débit volumique (en ou ), volume traversant la section par unité de temps ;
- le débit massique (en ) ;
- le débit molaire (en ) de l'espèce , relié à sa concentration et au débit volumique par , ou à sa fraction massique et au débit massique par .
Pour un système ouvert en régime permanent, l'accumulation de matière est nulle : le débit massique total entrant est égal au débit massique total sortant,
La masse totale se conserve ; en revanche le débit molaire d'une espèce donnée peut varier par l'effet des réactions chimiques.
Sur un système ouvert quelconque, en régime permanent, le bilan de l'espèce s'écrit
où est la quantité de produite par la réaction et par unité de temps ( : nombre stœchiométrique algébrique, : vitesse volumique de réaction). Pour un réactif, et (consommation).
Méthode : Exploiter un schéma de procédé et boucler les bilans
- Nommer chaque opération unitaire et son rôle ; repérer les courants d'entrée, de sortie, de recyclage et de purge.
- Choisir le système avant d'écrire un bilan : la boucle entière (son contour ne coupe que les courants externes) donne les débits externes sans rien savoir du recyclage ; le réacteur fait intervenir la conversion par passe ; un nœud de mélange relie courants frais et recyclés.
- Commencer par l'espèce dont le bilan a le moins de termes : un inerte n'est ni produit ni consommé et ne sort que par la purge, qui est donc indispensable.
- Distinguer la conversion par passe (dans le réacteur) et la conversion globale (sur la boucle entière, fixée par les pertes à la purge).
5.2 Le réacteur parfaitement agité continu (RPAC)
Le réacteur parfaitement agité continu (RPAC, en anglais CSTR) est un réacteur ouvert muni d'une agitation supposée idéale : le mélange y est instantanément homogène. En conséquence, la composition (et la température) sont uniformes dans tout le volume et égales à celles du courant de sortie. Le milieu réactionnel est un liquide de masse volumique constante : la conservation de la masse impose alors l'égalité des débits volumiques d'entrée et de sortie, , même si la réaction change la quantité de matière totale.
5.2.1 Bilan de matière et temps de passage
Le temps de passage du réacteur de volume parcouru par le débit volumique est
À masse volumique constante, c'est exactement le temps de séjour moyen d'un élément de fluide dans le réacteur, quel que soit le modèle d'écoulement ; dans un RPAC, les temps de séjour individuels sont très dispersés autour de cette moyenne.
Le mélange étant homogène, la vitesse est uniforme et le terme de production vaut simplement . Pour le réactif (), le bilan en régime permanent donne
où est évaluée aux conditions de sortie . C'est une équation algébrique (et non différentielle) : là réside toute la simplicité du RPAC. Pour une équation de nombre stœchiométrique , le terme de disparition est .
Le taux de conversion du réactif est la fraction de consommée :
Ainsi , avec .
En reportant dans le bilan et avec :
Cas particuliers () :
Pour l'ordre 1, la conversion tend vers lorsque , mais ne l'atteint jamais.
Le bilan du RPAC donne la vitesse sans dériver : , à la composition de sortie. En faisant varier le débit, un RPAC de laboratoire mesure une loi de vitesse. La concentration d'entrée est celle après mélange de tous les courants qui alimentent le réacteur.
Méthode : Dimensionner un RPAC
Pour atteindre un taux de conversion visé avec un débit donné :
- exprimer la vitesse à la sortie, en fonction de ;
- écrire le bilan et en tirer le volume
- vérifier l'ordre de grandeur du temps de passage .
5.3 Le réacteur en écoulement piston
Dans le réacteur en écoulement piston (en anglais PFR), le fluide progresse par tranches successives sans mélange axial : chaque tranche se comporte comme un petit réacteur fermé qui avance à vitesse constante. La composition varie continûment le long du réacteur, de l'entrée vers la sortie. On suppose l'écoulement isotherme et les débits volumiques égaux ( constant).
Sur une tranche de volume située à l'abscisse repérée par le volume cumulé , le bilan sur en régime permanent s'écrit
En introduisant le taux de conversion (, donc ), on obtient l'équation de dimensionnement
En PC, seule la loi de vitesse d'ordre 1 est exigible pour le réacteur piston ; l'équation de dimensionnement reste valable pour toute loi de vitesse.
Pour une réaction d'ordre 1 () à constant, et donnent
On retrouve la cinétique du réacteur fermé, le temps de passage y jouant le rôle du temps de réaction.
Pour une cinétique d'ordre positif, atteindre une même conversion demande moins de volume en écoulement piston qu'en RPAC. En effet, le RPAC opère entièrement à la faible concentration de sortie (donc à faible vitesse), tandis que le piston bénéficie des fortes concentrations le long du réacteur.
L'écart croît avec la conversion visée : pour , pour , pour . Une cascade de RPAC identiques en série ( à l'ordre
- ou un réacteur à recyclage sont intermédiaires entre les deux modèles idéaux.
5.4 Complément (hors programme) : sélectivité, rendement et taux de conversion
Les notions de cette section dépassent la lettre des programmes de PSI et de PC ; elles sont d'usage courant dans l'industrie et éclairent les choix de procédé.
Pour une réaction produisant un produit désiré à partir du réactif :
- Taux de conversion : (fraction de transformée) ;
- Sélectivité : (discrimine réaction utile et réactions parasites) ;
- Rendement : (fraction de effectivement valorisée en ).
Un fort taux de conversion n'est pas synonyme de bon procédé : si les réactions secondaires s'intensifient à forte conversion, la sélectivité chute et le rendement peut se dégrader. L'optimisation industrielle arbitre entre conversion (recyclage du réactif non transformé) et sélectivité (pureté, coût de séparation).
5.5 Bilan thermique d'un RPAC et points de fonctionnement
Une réaction exothermique dégage de la chaleur ; sans évacuation suffisante, la température s'élève, ce qui accélère la réaction (Arrhenius) et dégage encore plus de chaleur : c'est le risque d'emballement thermique. Le bilan énergétique permet de prévoir la température de fonctionnement, le flux thermique à échanger, et la stabilité du fonctionnement.
On considère un RPAC de température uniforme , alimenté à , siège d'une réaction unique d'enthalpie standard , de taux de conversion ; débits volumiques égaux, masse volumique et capacité thermique massique constantes ; le réacteur cède le flux thermique au fluide caloporteur. Alors
soit
Avec une paroi d'échange de coefficient global et un caloporteur à , .
Démonstration (Démonstration)
Système : le contenu du réacteur, ouvert, en régime stationnaire. Le premier principe pour un système ouvert en écoulement stationnaire s'écrit, par unité de temps, , où est la puissance thermique reçue (on néglige le travail de l'agitateur et les variations d'énergie cinétique et potentielle). L'enthalpie étant une fonction d'état, on calcule sur un chemin fictif :
- échauffer le flux entrant, sans réaction, de à : ;
- faire réagir à avec l'avancement par unité de temps : .
L'enthalpie de réaction en solution est confondue avec (solutions diluées, effets de mélange négligés), indépendante de sur l'intervalle étudié. On obtient le bilan annoncé.
Le flux à évacuer n'est pas la puissance libérée par la réaction : une alimentation plus froide que le réacteur en absorbe une partie. Si est assez basse, : il faut chauffer. Pour , le réacteur est autotherme.
On regroupe le bilan sous la forme avec, en utilisant et l'expression du RPAC pour l'ordre 1 :
Avec , la courbe de production a une forme en S (sigmoïde), qui plafonne à ; l'évacuation est affine en , de pente .
Les points de fonctionnement sont les températures pour lesquelles (intersections des deux courbes). Selon la position et la pente de la droite d'évacuation, il peut y en avoir un ou trois. Un point est stable si, autour de lui,
(l'évacuation « rattrape » toute élévation de température). Lorsqu'il existe trois points, celui du milieu est instable : une infime perturbation fait basculer le réacteur vers le point bas (extinction) ou vers le point haut (allumage, et risque d'emballement).
Pourquoi ce critère. Si la température du contenu s'écarte de du point de fonctionnement, le contenu reçoit la puissance nette , qui ramène la température si . Pour un RPAC adiabatique, où la composition suit la température, l'argument est rigoureux ; avec échange de chaleur, la condition est nécessaire, mais pas toujours suffisante (des oscillations sont possibles, hors programme). Le point du milieu : est positive en et négative à haute température ; avec trois racines, le signe de alterne (, , ), et la racine du milieu est instable.
paramètres de la capacité numérique, C) : les trois intersections sont
calculées ; celle du milieu est instable.Si le réacteur est adiabatique (), le bilan énergétique se réduit à une relation affine entre et la conversion, la droite adiabatique :
, élévation adiabatique de température, est l'échauffement d'un fluide qui convertirait tout son réactif sans rien échanger ; elle ne dépend que de l'alimentation.
Méthode : Température de fonctionnement d'un RPAC adiabatique
- Bilan de matière : (ordre 1 : ), courbe en S.
- Bilan énergétique : droite adiabatique .
- La température de fonctionnement est l'intersection des deux, cherchée par itérations, graphiquement ou numériquement ; s'il y a trois intersections, discuter la stabilité.
- Ne pas lire la température sur la seule droite adiabatique avec une conversion « mesurée à une autre température » : la conversion dépend elle-même de .
Quand varie, la droite d'évacuation se translate. Les deux positions où elle est tangente à la sigmoïde délimitent une fenêtre de températures d'entrée où trois points coexistent : en réchauffant l'alimentation, le réacteur reste éteint jusqu'à la tangence basse, puis saute sur la branche haute (allumage) ; en la refroidissant, il reste allumé jusqu'à la tangence haute, puis s'effondre (extinction). L'état du réacteur dépend de son histoire : c'est une hystérésis.
5.5.1 Sécurité des réacteurs : flux thermique et régulation
- Puissances : en chaque situation (marche normale, panne d'alimentation, panne de refroidissement), comparer la puissance libérée par la réaction à la puissance évacuable ; l'alimentation froide est elle-même un moyen de refroidissement, et sa perte peut suffire à emballer un réacteur.
- Accumulation : le réactif non converti contenu dans le réacteur réagira même en cas de panne ; la température maximale atteignable se calcule sur ce qui est accumulé, comme dans un réacteur fermé adiabatique.
- Régulation : la température d'entrée et le débit de caloporteur sont régulés ; un refroidissement plus fort raidit la droite d'évacuation et supprime la multiplicité des points ; une sécurité coupe l'alimentation en réactif sur température haute.
- Changement d'échelle : à géométrie semblable, la surface d'échange par unité de volume décroît comme ; la puissance évacuable par la paroi ne suit pas la puissance produite.
Capacité numérique : Points de fonctionnement d'un RPAC
On cherche les températures stationnaires d'un RPAC exothermique du premier ordre en résolvant , puis on discute leur stabilité. Les paramètres sont ceux de l'hydratation de l'oxyde de propylène en propane-1,2-diol, catalysée par l'acide sulfurique, dans une unité modèle (réacteur adiabatique de , alimentation à ).
import numpy as np
from scipy.optimize import brentq
# --- parametres (unites SI) : hydratation de l'oxyde de propylene ---
R = 8.314 # J/mol/K
Ea = 75.3e3 # J/mol (energie d'activation)
A = 4.71e9 # 1/s (facteur pre-exponentiel)
dHr = -92.6e3 # J/mol (exothermique : dHr < 0)
CAe = 2.5e3 # mol/m3 (concentration d'entree)
Q = 9.6 / 3600 # m3/s (debit volumique)
V = 1.2 # m3 (volume du reacteur)
tau = V / Q # s (temps de passage)
rho_cp = 3.8e6 # J/m3/K (rho * c_p du melange)
hA = 0.0 # W/K (echange thermique : adiabatique)
Te = 290.15 # K (entree, 17 C)
Tc = 288.15 # K (caloporteur)
def X(T): # conversion RPAC, ordre 1
k = A * np.exp(-Ea / (R * T))
return k * tau / (1.0 + k * tau)
def Qprod(T): # puissance liberee (W)
return (-dHr) * Q * CAe * X(T)
def Qevac(T): # puissance evacuee (W), affine en T
return Q * rho_cp * (T - Te) + hA * (T - Tc)
def bilan(T):
return Qprod(T) - Qevac(T)
# --- recherche des racines par changement de signe ---
Tgrid = np.linspace(270.0, 400.0, 4000)
f = bilan(Tgrid)
racines = []
for i in range(len(Tgrid) - 1):
if f[i] * f[i + 1] < 0.0:
racines.append(brentq(bilan, Tgrid[i], Tgrid[i + 1]))
print("Points de fonctionnement (T, X) et stabilite :")
for T in racines:
# stable si pente evacuation > pente production (derivees numeriques)
dT = 1.0e-2
dprod = (Qprod(T + dT) - Qprod(T - dT)) / (2 * dT)
devac = (Qevac(T + dT) - Qevac(T - dT)) / (2 * dT)
etat = "stable" if devac > dprod else "INSTABLE"
print(f" T = {T:6.1f} K, X = {X(T):.3f} -> {etat}")
Le programme localise les intersections par balayage puis les affine (brentq). La stabilité est jugée en comparant les pentes locales. Son exécution donne :
Points de fonctionnement (T, X) et stabilite :
T = 296.3 K, X = 0.101 -> stable
T = 314.9 K, X = 0.406 -> INSTABLE
T = 344.3 K, X = 0.888 -> stable
Trois points : éteint et stable (C, ), instable (C), allumé et stable (C, ). Avec un échange et un caloporteur à , il ne reste qu'un point, éteint ( K).
5.6 Complément (hors programme) : optimisation, cinétique contre thermodynamique
Cette section dépasse la lettre des programmes de PSI et de PC ; elle relie ce chapitre à l'optimisation des équilibres étudiée dans le cours sur les équilibres.
Pour une réaction exothermique et renversable (par exemple la synthèse de l'ammoniac ), deux effets s'opposent avec la température :
- cinétique : élever accélère la réaction (constante de vitesse croissante) ;
- thermodynamique : élever déplace l'équilibre dans le sens endothermique (loi de van 't Hoff, décroissant), donc abaisse la conversion maximale accessible.
Le taux de conversion à l'équilibre décroît quand augmente, alors que la vitesse d'approche de l'équilibre croît. Il existe donc une température (ou un profil de température) optimale maximisant la production.
Méthode : Trajectoire optimale de température
En pratique, on fait fonctionner le procédé à température élevée au début (là où l'on est loin de l'équilibre : la vitesse commande), puis on refroidit progressivement à mesure que l'on s'approche de l'équilibre (là où la thermodynamique commande). C'est le principe des convertisseurs à plusieurs lits de catalyseur du procédé Haber-Bosch, refroidis entre deux lits soit par des échangeurs, soit par trempe, c'est-à-dire injection de gaz de synthèse froid (quench ou cold shot).
5.7 Synthèse du chapitre
- Débits : ; en régime permanent, bilan sur : , sur un système choisi (boucle, réacteur, nœud) ; un inerte impose une purge.
- RPAC : composition uniforme = sortie ; bilan algébrique ; temps de passage ; ordre 1 : .
- Piston : bilan sur une tranche ; dimensionnement ; ordre 1 : .
- À conversion égale, (ordre positif), d'autant plus nettement que la conversion est forte.
- Bilan énergétique du RPAC : ; droite adiabatique ; le flux à évacuer n'est pas la chaleur de réaction.
- Points de fonctionnement : ; 1 ou 3 ; stable si ; le point du milieu est instable ; allumage, extinction, hystérésis.
- Sécurité : puissances comparées dans chaque situation, température maximale calculée sur ce qui est accumulé, surface d'échange par unité de volume en .
- Compléments : rendement ; exothermique renversable : compromis vitesse (haute ) / équilibre (basse ).
Exercices du chapitre
- Hydratation de l'éthylène : lire le schéma, boucler les bilans corrigé
- Le RPAC, instrument de cinétique : saponification de l'acétate d'éthyle corrigé
- Bassin de chloration : pourquoi on le cloisonne corrigé
- Réacteur à boucle : du piston au RPAC par le recyclage corrigé
- Établir le bilan énergétique d'un RPAC : l'hydrolyse de l'anhydride acétique corrigé
- Sécurité d'un RPAC de polymérisation du styrène : accumulation et changement d'échelle corrigé