Adloun

Procédés industriels continus : réacteurs ouverts

Cours complet · chimie : thermodynamique et électrochimie (PC), chapitre 5 · CPGE PC (2e année)

Travailler ce chapitre 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 thermiques 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é, ainsi que les compromis entre optimisation cinétique et optimisation 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, chauffage, distillation, filtration, séparation membranaire, …). 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. 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 .
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).

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. On suppose les débits volumiques d'entrée et de sortie égaux, (mélange incompressible sans variation de quantité de matière globale).

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

Il représente l'ordre de grandeur du temps de séjour moyen d'un élément de fluide dans le réacteur.

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.

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.

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

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

Figure : Taux de conversion en fonction du temps de passage (ordre 1) : le piston convertit plus

vite que le RPAC à donné.</div>

5.4 Sélectivité, rendement et taux de conversion

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 et sa stabilité.

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

En régime permanent, sur un RPAC de température uniforme , alimenté à , siège d'une réaction unique de vitesse , le premier principe en système ouvert donne (débits volumiques égaux, masse volumique et capacité thermique massique supposées constantes) :

où pour une réaction exothermique, le coefficient global d'échange et la température du fluide caloporteur.

Définition 5.12Chaleur produite et chaleur évacuée

On regroupe le bilan sous la forme avec, en utilisant et l'expression du RPAC pour l'ordre 1 :

où regroupe et . Avec , la courbe de production a une forme en S (sigmoïde), tandis que l'évacuation est affine en .

ImportantPoints de fonctionnement et stabilité

Les points de fonctionnement sont les températures pour lesquelles (intersections des deux courbes). Selon 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 / emballement).

Figure : Points de fonctionnement d'un RPAC exothermique : trois intersections, celle du milieu

est instable.</div>

Proposition 5.13Fonctionnement adiabatique

Si le réacteur est adiabatique (), le bilan thermique se réduit à une relation affine entre et la conversion, la droite adiabatique :

La température de fonctionnement se lit à l'intersection de cette droite avec la courbe imposée par la cinétique du RPAC.

5.6 Optimisation : cinétique contre thermodynamique

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 étages de catalyseur avec refroidissement intermédiaire (« cold shot ») du procédé Haber-Bosch.

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


import numpy as np
from scipy.optimize import brentq

# --- parametres physiques (unites SI) ---
R      = 8.314          # J/mol/K
Ea     = 7.0e4          # J/mol   (energie d'activation)
A      = 5.0e10         # 1/s     (facteur pre-exponentiel)
dHr    = -8.0e4         # J/mol   (exothermique : dHr < 0)
CAe    = 2.0e3          # mol/m3  (concentration d'entree)
Q      = 1.0e-3         # m3/s    (debit volumique)
V      = 1.0            # m3      (volume du reacteur)
tau    = V / Q          # s       (temps de passage)
rho_cp = 4.0e6          # J/m3/K  (rho * c_p du melange)
hA     = 3.0e3          # W/K     (echange thermique global)
Te     = 300.0          # K       (entree)
Tc     = 300.0          # 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):                   # chaleur liberee (W)
    return (-dHr) * Q * CAe * X(T)

def Qevac(T):                   # chaleur evacuee (W), affine en T
    Tstar = (Q * rho_cp * Te + hA * Tc) / (Q * rho_cp + hA)
    return (Q * rho_cp + hA) * (T - Tstar)

def bilan(T):
    return Qprod(T) - Qevac(T)

# --- recherche des racines par changement de signe ---
Tgrid = np.linspace(280.0, 900.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 : trois solutions apparaissent typiquement (point bas éteint stable, point milieu instable, point haut allumé stable).

5.7 Synthèse du chapitre

À retenir
  • Débits : ; en régime permanent, bilan sur : .
  • RPAC : composition uniforme = sortie ; bilan algébrique ; temps de passage ; ordre 1 : .
  • Piston : bilan sur une tranche ; dimensionnement ; ordre 1 : .
  • À conversion égale, (ordre positif).
  • Conversion / sélectivité / rendement : .
  • Bilan thermique : ; 1 ou 3 points de fonctionnement ; stable si .
  • Exothermique renversable : compromis vitesse (haute ) / équilibre (basse ) température optimale décroissante le long du procédé.

<i class="fa-solid fa-dumbbell mr-2" style="color:#2E7559"></i>5.8 Exercices résolus

Exemple 5.14Temps de passage d'un RPAC, ordre 1

Une réaction d'ordre 1 a pour constante . On veut atteindre dans un RPAC. Déterminer le temps de passage requis.

Démonstration

Pour l'ordre 1 en RPAC, , d'où Le temps de passage vaut .

Exemple 5.15Comparaison des volumes RPAC / piston

Même réaction (, ), . Calculer les volumes d'un RPAC et d'un réacteur piston, puis leur rapport.

Démonstration

RPAC : , soit .
Piston : , soit .
Rapport : le piston est près de quatre fois plus compact à forte conversion.

Exemple 5.16Dimensionnement pour un débit imposé

On traite d'une solution à . La réaction est d'ordre 1, . Quel volume de RPAC pour ?

Démonstration

. Avec , .

Exemple 5.17Réaction d'ordre 2 en RPAC

Pour d'ordre 2 (, ), , on souhaite . Déterminer .

Démonstration

Bilan RPAC : . Donc

Exemple 5.18Rendement et sélectivité

Un réactif donne le produit désiré () et un sous-produit (). En sortie, et l'on mesure . Calculer la sélectivité et le rendement en .

Démonstration

consommé : . transformé en : . Sélectivité . Rendement (cohérent avec ).

Exemple 5.19Température adiabatique de fonctionnement

Un RPAC adiabatique traite , entrée à , , . Quelle est la température de sortie si ?

Démonstration

Droite adiabatique : L'élévation adiabatique maximale () serait de .

Exemple 5.20Flux thermique à évacuer

Le RPAC précédent fonctionne maintenant en isotherme à avec le même dégagement. Le débit vaut et . Quel flux thermique le système de refroidissement doit-il évacuer ?

Démonstration

Chaleur libérée : En isotherme, tout ce flux doit être évacué par le caloporteur.

5.9 Exercices d'entraînement

  • Une réaction d'ordre 1 () est conduite en RPAC avec . Calculer le taux de conversion. Que devient-il pour ?
  • Établir, à partir du bilan sur une tranche, l'expression pour un réacteur piston d'ordre 1. Calculer pour si .
  • Un procédé traite , , réaction d'ordre 1, . Dimensionner un RPAC puis un piston pour ; comparer les volumes.
  • Pour une réaction d'ordre 0 (, ) en RPAC, exprimer en fonction de . À partir de quel la conversion serait-elle totale si ?
  • Démontrer que, pour une cinétique d'ordre 1, le rapport vaut et calculer sa valeur pour .
  • On associe en série deux RPAC identiques de temps de passage chacun (réaction d'ordre 1). Montrer que la conversion globale est et comparer à un RPAC unique de temps de passage .
  • Une réaction d'ordre 2 (, ) est conduite en piston. Établir la relation entre et , puis calculer pour .
  • Un RPAC exothermique adiabatique reçoit , , , . Établir la droite adiabatique et donner l'élévation maximale de température.
  • Écrire le bilan énergétique complet d'un RPAC non adiabatique et identifier les termes de chaleur produite et de chaleur évacuée. Préciser le critère de stabilité d'un point de fonctionnement.
  • Expliquer pourquoi un RPAC exothermique peut présenter trois points de fonctionnement. Lequel est instable ? Que se passe-t-il physiquement lors d'un « allumage » ou d'une « extinction » ?
  • Pour la synthèse de l'ammoniac (exothermique), justifier qualitativement l'existence d'une température optimale et décrire la stratégie industrielle de conduite en température (étages refroidis).
  • À l'aide de la capacité numérique du chapitre, étudier l'effet d'une augmentation du coefficient d'échange sur le nombre de points de fonctionnement. Interpréter en termes de sécurité du réacteur.

Continuer sur Adloun : animation, QCM, fiches, exercices