Adloun

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

Un réacteur ouvert parfaitement agité : les réactifs entrent, le mélange est homogène et les produits sortent en continu.

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.

Définition 5.1Procédés discontinu et continu

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.

Définition 5.2Débits

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 .
Proposition 5.3Bilan de matière global en régime permanent

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.

ImportantBilan de matière sur une espèce

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).

Figure : Schéma de procédé élémentaire : réaction, séparation, recyclage du réactif non converti et purge. Le contour en tirets ne coupe pas le recyclage.

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)

Définition 5.4Modèle du 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.

Figure : Réacteur parfaitement agité continu : composition uniforme, égale à celle de la sortie.

5.2.1 Bilan de matière et temps de passage

Définition 5.5Temps 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.

ImportantBilan sur le réactif dans un RPAC

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 .

Définition 5.6Taux de conversion

Le taux de conversion du réactif est la fraction de consommée :

Ainsi , avec .

Proposition 5.7Taux de conversion en fonction du temps de passage (RPAC)

En reportant dans le bilan et avec :

Cas particuliers () :

Pour l'ordre 1, la conversion tend vers lorsque , mais ne l'atteint jamais.

iRemarque

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

Définition 5.8Modèle du réacteur 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).

Figure : Réacteur piston : bilan sur une tranche élémentaire de volume .
ImportantBilan sur une tranche (réacteur piston)

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.

Proposition 5.9Conversion dans un réacteur piston, ordre 1

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.

Proposition 5.10Piston contre RPAC : comparaison

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

  1. ou un réacteur à recyclage sont intermédiaires entre les deux modèles idéaux.
Figure : Taux de conversion en fonction du temps de passage (ordre 1) : le piston convertit plus vite que le RPAC à donné.

5.4 Complément (hors programme) : sélectivité, rendement et taux de conversion

iRemarqueComplément

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é.

Définition 5.11Conversion, sélectivité, rendement

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 ).
iRemarque

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.

ImportantBilan énergétique d'un RPAC en régime permanent

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é.

Attention

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.

Définition 5.12Puissance produite et puissance évacuée

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 .

ImportantPoints de fonctionnement et stabilité

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).

iRemarque

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.

Figure : Points de fonctionnement d'un RPAC adiabatique (hydratation de l'oxyde de propylène,

paramètres de la capacité numérique, C) : les trois intersections sont

calculées ; celle du milieu est instable.
Proposition 5.13Fonctionnement adiabatique

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 .
iRemarqueAllumage, extinction, hystérésis

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

ImportantRaisonner en puissances et en accumulation
  • 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

iRemarqueComplément

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.

ImportantLe compromis exothermique

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

À retenir
  • 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

La suite du cours · Tous les exercices du chapitre