Procédé chimiqueProcédés et eau

Réacteur chimique (CSTR)

Un réacteur agité exothermique et son bestiaire de dynamique non linéaire — multiplicité des états stationnaires, hystérésis d’amorçage/extinction, oscillations entretenues en cycle limite, emballement thermique, et la régulation par rétroaction qui apprivoise un point de fonctionnement instable en boucle ouverte.

DjiniousLabModélisation et simulation

Une courbe d’états stationnaires en forme de S, température du réacteur en fonction de la température du réfrigérant, avec des branches stable inférieure et supérieure, une branche instable médiane, et deux points de pli marquant l’amorçage et l’extinction

états stationnaires (une cuve): 3états stationnaires (une cuve)
hystérésis d’amorçage/extinction: 5,1 Khystérésis d’amorçage/extinction
oscillation en cycle limite: ~134 Koscillation en cycle limite
à-coup d’alimentation qui s’emballe: +20 Kà-coup d’alimentation qui s’emballe
taux de conversion de dimensionnement: 90.1%taux de conversion de dimensionnement
notebooks de conception: 10notebooks de conception

Une cuve, trois états stationnaires, et un emballement.

Placez une réaction exothermique dans une cuve agitée refroidie, et le réacteur le plus simple de la chimie se transforme en vitrine de dynamique non linéaire. La chaleur produite croît fortement avec la température tandis que la chaleur évacuée ne croît que linéairement — et l’endroit où ces courbes se croisent peut donner un point de fonctionnement ou trois, un amorçage et une extinction brusques, des oscillations auto-entretenues, ou un emballement thermique. Ce programme cartographie tout cela à partir des bilans couplés de matière et d’énergie, puis montre comment la régulation par rétroaction maintient le réacteur là où il ne pourrait sinon pas tenir.

Chaleur produite contre chaleur évacuée.

La vitesse de réaction croît exponentiellement avec la température, donc la chaleur qu’elle génère suit une courbe en S ; la double enveloppe de refroidissement évacue la chaleur selon une droite. Là où le S croise la droite se trouve un état stationnaire — et un S peut croiser une droite en trois points. Le croisement du milieu est instable, les deux extrêmes sont stables, et à mesure que l’on fait varier la température du réfrigérant, les états stables apparaissent et disparaissent en des points de pli. C’est le diagramme de van Heerden, et c’est pourquoi le même réacteur peut se trouver froid et à peine réactif, ou chaud et presque complet, selon son seul historique.

Le diagramme de van Heerden — une courbe de génération de chaleur en forme de S croisant une droite d’évacuation de chaleur en un ou trois points — Le diagramme de van Heerden : la courbe en S exponentielle de génération de chaleur face à la droite linéaire d’évacuation de chaleur. Une seule intersection signifie un point de fonctionnement unique ; trois signifient une multiplicité. Incliner ou décaler la droite d’évacuation — en changeant la température ou le débit du réfrigérant — est ce qui fait naître et disparaître les états stationnaires.

La courbe en S des états stationnaires, température du réacteur en fonction de la température du réfrigérant, avec des branches stables, une branche de selle instable, et des points de pli d’amorçage/extinction — La courbe en S de multiplicité, tracée par continuation et classée par stabilité : une branche stable froide, une branche stable chaude, et une selle instable entre les deux. Les points de pli sont l’amorçage (303,2 K) et l’extinction (298,1 K) — une fenêtre de 5,1 K où l’état du réacteur dépend de son point de départ.

Il peut osciller, et il peut s’emballer.

La multiplicité n’est pas la seule surprise. Autour de l’état médian instable, le réacteur peut se caler sur un cycle limite auto-entretenu — la température oscillant de plus de cent degrés, indéfiniment, sans aucune variation d’entrée — une bifurcation de Hopf. Et si l’on pousse trop l’alimentation, la génération de chaleur dépasse complètement le refroidissement : une perturbation de +10 K est absorbée, mais une perturbation de +20 K déclenche un emballement thermique à plus de 460 K. Le programme trace l’orbite du cycle limite dans le plan de phase et cartographie la frontière d’emballement — l’enveloppe de sécurité à l’intérieur de laquelle un réacteur réel doit rester.

Une orbite fermée de cycle limite dans le plan de phase concentration-température, avec des trajectoires transitoires s’enroulant vers elle depuis l’intérieur et l’extérieur — Le cycle limite de Hopf : une orbite fermée dans le plan concentration-température que le réacteur parcourt indéfiniment, avec des transitoires qui s’y enroulent des deux côtés autour de l’état stationnaire instable situé en son centre. Une oscillation auto-entretenue à partir d’une alimentation constante — le réacteur comme horloge chimique.

Une carte de sécurité de la température maximale du réacteur en fonction de l’amplitude de la perturbation et des conditions de fonctionnement, avec une zone d’emballement — La carte de sécurité d’emballement : température maximale en fonction de l’amplitude d’une perturbation d’alimentation et des conditions de fonctionnement. La perturbation tolérable diminue à mesure que la double enveloppe se réchauffe — la frontière de sensibilité paramétrique qui sépare un réacteur qui se rétablit d’un réacteur qui s’emballe.

Un régulateur PI maintenant le réacteur à un point de consigne instable en boucle ouverte qui oscillerait autrement, se stabilisant proprement — La régulation ferme la boucle : un régulateur PI agissant sur le réfrigérant maintient le réacteur au point de fonctionnement de 350 K, instable en boucle ouverte — où, livré à lui-même, il oscillerait. Le point de fonctionnement le plus productif est souvent celui que l’on ne peut atteindre que par rétroaction.

Chaque valeur est une valeur que vous pouvez recalculer.

Le notebook de validation finale re-dérive chaque exigence à partir des mêmes bilans couplés que construit le programme.

Résultat

États stationnaires (Tc = 301 K) : 3
Hystérésis d’amorçage/extinction : 5,1 K
Amplitude du cycle limite : ~134 K
Seuil d’emballement : perturbation de +20 K
Exigences vérifiées : 7 / 7

Exigence

États stationnaires (Tc = 301 K) : multiplicité
Hystérésis d’amorçage/extinction : 303,2 / 298,1 K
Amplitude du cycle limite : entretenue (Hopf)
Seuil d’emballement : +10 K absorbé
Exigences vérifiées : RÉUSSI

Une seule réaction, tout le bestiaire.

Le réacteur fait tourner une seule réaction exothermique du premier ordre A→B dans une cuve parfaitement mélangée à paramètres localisés et propriétés physiques constantes — le modèle canonique de Seborg, et non une thermochimie détaillée ni un réseau multi-réactions. Ce modèle minimal est délibérément ce qui montre proprement tout le bestiaire de la dynamique non linéaire : multiplicité, hystérésis, cycle limite de Hopf, emballement, et stabilisation par rétroaction, le tout traçable à la main. C’est exactement la fidélité dont ont d’abord besoin la conception d’un réacteur et l’intuition de sécurité — trouver les points de fonctionnement, cartographier l’enveloppe d’emballement, et concevoir la régulation — avant un modèle cinétique détaillé ou un modèle CFD.