Le comportement dynamique non linéaire d'un CSTR non isotherme est linéarisé en développant ses équations de bilan non linéaires autour d'un point de fonctionnement à l'état stationnaire et en ne conservant que les termes du premier ordre. La principale non-linéarité provient de la vitesse de réaction d'Arrhenius, (k(T)c_A = k_0 e^{-E/(RT)}c_A), car la température affecte la vitesse de réaction de manière exponentielle tandis que la concentration l'affecte directement. Le modèle obtenu, exprimé en variables d'écart, est un ensemble couplé d'équations différentielles linéaires adapté à l'analyse de la stabilité, de la réponse transitoire et du contrôle.

Un CSTR non isotherme est linéarisé localement et non globalement : les bilans de matière et d'énergie non linéaires sont remplacés par leur approximation de Taylor du premier ordre autour d'un état stationnaire sélectionné. Cela produit un modèle linéaire d'état-espace qui décrit avec précision les faibles écarts de concentration, de température, de débit, de conditions d'alimentation et de puissance de refroidissement.

Pourquoi le modèle du CSTR est non linéaire

Le bilan de composant

Pour une réaction du premier ordre, le bilan du réactif peut s'écrire

[ \frac{dc_A}{dt}

\frac{F_i}{V}(c_{Ai}-c_A)

k_0e^{-E/(RT)}c_A ]

Le terme non linéaire est la vitesse de réaction,

[ r_A=k_0e^{-E/(RT)}c_A ]

La dépendance exponentielle à la température est la principale source de non-linéarité.

Le bilan d'énergie

Le bilan énergétique du réacteur est

[ \frac{dT}{dt}

\frac{F_i}{V}(T_i-T) + \frac{(-\Delta H_r)k_0e^{-E/(RT)}c_A}{\rho c_p}

\frac{Q}{\rho c_pV} ]

Pour une réaction exothermique, une augmentation de la température accroît la vitesse de réaction, ce qui génère une chaleur supplémentaire. Cette chaleur peut à son tour augmenter davantage la température, créant la rétroaction responsable de la sensibilité thermique, de l'existence de plusieurs états stationnaires et d'une instabilité potentielle.

Le bilan de volume

Le bilan massique total est

[ \frac{dV}{dt}=F_i-F ]

Si les débits d'entrée et de sortie sont maintenus à des valeurs égales, le volume du réacteur est approximativement constant et les principales variables d'état dynamiques sont la concentration et la température. Si les dynamiques de niveau ou de débit sont importantes, (V) est conservé comme variable d'état supplémentaire lors de la linéarisation.

Sélection du point de fonctionnement

Conditions d'état stationnaire

Choisissez un état stationnaire nominal défini par

[ \bar{c}_A,\quad \bar{T},\quad \bar{V},\quad \bar{F}i,\quad \bar{c}{Ai},\quad \bar{T}_i,\quad \bar{Q} ]

À l'état stationnaire, les dérivées des variables d'état sont nulles :

[ \left.\frac{dc_A}{dt}\right|{\mathrm{ss}}=0, \qquad \left.\frac{dT}{dt}\right|{\mathrm{ss}}=0, \qquad \left.\frac{dV}{dt}\right|_{\mathrm{ss}}=0 ]

Les bilans à l'état stationnaire établissent le point de référence autour duquel le modèle est approché.

Variables d'écart

Chaque variable de procédé est exprimée comme la somme de sa valeur à l'état stationnaire et d'un faible écart :

[ c_A=\bar{c}_A+\Delta c_A ]

[ T=\bar{T}+\Delta T ]

[ V=\bar{V}+\Delta V ]

Les entrées et les perturbations sont traitées de manière similaire :

[ Q=\bar{Q}+\Delta Q, \qquad c_{Ai}=\bar{c}{Ai}+\Delta c{Ai}, \qquad T_i=\bar{T}_i+\Delta T_i ]

Les variables d'écart décrivent l'éloignement du procédé par rapport à sa condition nominale de fonctionnement.

Linéarisation de la vitesse de réaction d'Arrhenius

Développement de Taylor du premier ordre

Définissons la fonction non linéaire de vitesse de réaction comme suit :

[ r(c_A,T)=k_0e^{-E/(RT)}c_A ]

Un développement de Taylor à deux variables autour de ((\bar{c}_A,\bar{T})) donne

[ r(c_A,T) \approx r(\bar{c}A,\bar{T}) + \left.\frac{\partial r}{\partial c_A}\right|{\mathrm{ss}}\Delta c_A + \left.\frac{\partial r}{\partial T}\right|_{\mathrm{ss}}\Delta T ]

Les termes faisant intervenir des produits tels que (\Delta c_A\Delta T), ainsi que les termes du deuxième ordre et des ordres supérieurs, sont négligés.

Évaluation des dérivées

La constante de vitesse nominale est

[ \bar{k}=k_0e^{-E/(R\bar{T})} ]

Par conséquent,

[ \left.\frac{\partial r}{\partial c_A}\right|_{\mathrm{ss}}

\bar{k} ]

[ \left.\frac{\partial r}{\partial T}\right|_{\mathrm{ss}}

\bar{k}\bar{c}_A \frac{E}{R\bar{T}^{2}} ]

La vitesse de réaction linéarisée devient

[ \Delta r \approx \bar{k}\Delta c_A + \bar{k}\bar{c}_A \frac{E}{R\bar{T}^{2}}\Delta T ]

ou, de manière équivalente,

[ r \approx \bar{r} + \bar{k}\Delta c_A + \bar{k}\bar{c}_A \frac{E}{R\bar{T}^{2}}\Delta T ]

Cette expression montre directement comment les perturbations de concentration et de température modifient la vitesse de réaction à proximité du point de fonctionnement choisi.

Construction du modèle dynamique linéarisé

Approximation à volume constant

Pour un réacteur à volume constant avec (F_i=F=\bar{F}), le bilan de concentration devient

[ \frac{dc_A}{dt}

\frac{\bar{F}}{\bar{V}}(c_{Ai}-c_A)

r(c_A,T) ]

Après substitution de l'approximation de Taylor et annulation des termes à l'état stationnaire,

[ \frac{d\Delta c_A}{dt}

\left( \frac{\bar{F}}{\bar{V}}+\bar{k} \right)\Delta c_A

\bar{k}\bar{c}A \frac{E}{R\bar{T}^{2}}\Delta T + \frac{\bar{F}}{\bar{V}}\Delta c{Ai} ]

Cette équation décrit la dynamique locale de la concentration comme une combinaison linéaire des écarts de concentration, de température et de concentration de l'alimentation.

Bilan d'énergie linéarisé

Définissons

[ \alpha=\frac{-\Delta H_r}{\rho c_p} ]

La contribution de la chaleur de réaction est (\alpha r). Le bilan d'énergie linéarisé est alors

[ \frac{d\Delta T}{dt}

\alpha\bar{k}\Delta c_A + \left[ -\frac{\bar{F}}{\bar{V}} + \alpha\bar{k}\bar{c}_A \frac{E}{R\bar{T}^{2}} \right]\Delta T + \frac{\bar{F}}{\bar{V}}\Delta T_i

\frac{1}{\rho c_p\bar{V}}\Delta Q ]

L'équation de température contient la rétroaction thermique locale. Le coefficient multipliant (\Delta T) reflète la compétition entre le refroidissement convectif et l'augmentation de la génération de chaleur provoquée par une hausse de température.

Prise en compte des dynamiques de volume ou de débit

Si le volume du réacteur n'est pas constant, le bilan de volume est également linéarisé :

[ \frac{d\Delta V}{dt}=\Delta F_i-\Delta F ]

Les facteurs (F_i/V) dans les bilans de concentration et d'énergie doivent alors être développés eux aussi. Par exemple, l'approximation du premier ordre du terme de dilution à l'entrée est

[ \frac{F_i}{V} \approx \frac{\bar{F}_i}{\bar{V}} + \frac{1}{\bar{V}}\Delta F_i

\frac{\bar{F}_i}{\bar{V}^{2}}\Delta V ]

Cela introduit un couplage entre le débit, le volume de liquide, la concentration et la température.

Expression du modèle sous forme d'état-espace

Vecteurs d'état, d'entrée et de perturbation

Un vecteur d'état typique est

[ x= \begin{bmatrix} c_A\ T\ V \end{bmatrix} ]

Pour une analyse à volume constant, il peut être réduit à

[ x= \begin{bmatrix} c_A\ T \end{bmatrix} ]

Les variables manipulées peuvent inclure la puissance de refroidissement (Q), le débit de sortie (F) ou le débit de fluide caloporteur dans une jaquette. La concentration de l'alimentation, la température d'alimentation et le débit d'entrée sont généralement considérés comme des perturbations.

Représentation par matrice jacobienne

Le modèle non linéaire peut s'écrire de manière générale

[ \dot{x}=f(x,u,d) ]

où (x) contient les variables d'état, (u) les entrées manipulées et (d) les perturbations. La linéarisation autour de l'état stationnaire donne

[ \Delta\dot{x}

A\Delta x+B\Delta u+E\Delta d ]

avec

[ A= \left.\frac{\partial f}{\partial x}\right|{\mathrm{ss}}, \qquad B= \left.\frac{\partial f}{\partial u}\right|{\mathrm{ss}}, \qquad E= \left.\frac{\partial f}{\partial d}\right|_{\mathrm{ss}} ]

L'équation de sortie s'écrit de manière similaire

[ \Delta y=C\Delta x+D\Delta u ]

Cette forme d'état-espace constitue la base standard de l'analyse du contrôle des procédés.

Ce que révèle le modèle linéaire

Stabilité locale

Les valeurs propres de la matrice d'état (A) déterminent la stabilité locale. Si toutes les valeurs propres ont des parties réelles négatives, les petites perturbations s'atténuent et l'état stationnaire sélectionné est localement stable.

Si au moins une valeur propre possède une partie réelle positive, les petits écarts augmentent et l'état stationnaire est localement instable.

Réponse transitoire

Le modèle linéaire prédit la réponse de la concentration et de la température à de faibles variations de la puissance de refroidissement, de la concentration d'alimentation, de la température d'alimentation ou du débit. Il peut être utilisé pour estimer la vitesse de réponse, le dépassement, les interactions entre variables et les constantes de temps dominantes.

Contrôle du transfert thermique

La puissance de refroidissement (Q) apparaît comme une entrée du bilan d'énergie. Le modèle linéarisé est ainsi utile pour évaluer les stratégies de contrôle de la température, notamment le contrôle par rétroaction ou en cascade faisant intervenir la température du réacteur et le débit dans la jaquette ou du fluide caloporteur.

Multiples points de fonctionnement

Un CSTR non isotherme peut présenter plusieurs états stationnaires, car la génération de chaleur par la réaction augmente fortement avec la température tandis que l'évacuation de chaleur dépend des mécanismes d'écoulement et de refroidissement. Chaque état stationnaire possède sa propre matrice jacobienne et donc son propre comportement dynamique local.

Un modèle linéaire développé autour d'un point de fonctionnement ne doit pas être automatiquement utilisé pour représenter la dynamique à proximité d'un autre point de fonctionnement.

Comprendre les compromis

Le modèle n'est précis que localement

L'approximation de Taylor est fiable pour des écarts suffisamment faibles par rapport à l'état stationnaire sélectionné. De grandes variations de température ou de concentration peuvent rendre significatifs les termes d'ordre supérieur négligés.

Pour le démarrage, l'arrêt, les perturbations importantes ou les transitions entre états stationnaires, le modèle non linéaire d'origine doit être simulé.

Une forte sensibilité thermique peut limiter la précision

Comme le terme d'Arrhenius varie de manière exponentielle avec la température, même des écarts de température modérés peuvent entraîner des variations importantes de la vitesse de réaction. Un modèle linéaire peut donc perdre sa précision plus rapidement dans un système fortement exothermique ou à haute énergie d'activation que dans un système faiblement sensible à la température.

Les hypothèses de volume constant peuvent masquer des dynamiques importantes

L'hypothèse (F_i=F) simplifie le modèle en supprimant la dynamique de volume. Cette hypothèse est inappropriée lorsque les variations de niveau, les différences entre les débits d'entrée et de sortie ou les dynamiques de contrôle du débit influencent de manière significative la réponse du réacteur.

La linéarisation ne détermine pas la stabilité globale

L'analyse des valeurs propres décrit le comportement à proximité d'un point de fonctionnement. Elle ne peut pas, à elle seule, établir le comportement du réacteur non linéaire dans toutes les conditions possibles de température, de concentration ou de débit.

Comment appliquer cette méthode à votre projet

La procédure de linéarisation appropriée dépend de l'objectif de l'analyse et des variables qui sont contrôlées activement.

Si votre objectif principal est la stabilité locale : Résolvez les bilans à l'état stationnaire, calculez la matrice jacobienne au point de fonctionnement et examinez les valeurs propres de (A).
Si votre objectif principal est le contrôle de la température : Conservez la dérivée de la vitesse d'Arrhenius par rapport à la température et utilisez la puissance de refroidissement ou le débit du fluide caloporteur comme entrée manipulée.
Si votre objectif principal est le contrôle de la composition : Incluez la concentration comme sortie contrôlée et analysez les effets de la concentration d'alimentation, du débit d'alimentation et des perturbations de température.
Si votre objectif principal est le contrôle du niveau ou du débit : Conservez le bilan de volume et linéarisez les termes (F_i/V) au lieu de supposer un volume constant.
Si votre objectif principal concerne les grands transitoires ou l'emballement thermique : Utilisez les bilans de matière et d'énergie entièrement non linéaires, car un seul modèle linéaire local ne peut pas représenter le comportement thermique global.

Le modèle de contrôle des procédés essentiel est obtenu en évaluant les premières dérivées des bilans non linéaires du CSTR à l'état stationnaire, ce qui transforme le réacteur en un système linéaire local pour une analyse rigoureuse en petits signaux.

Tableau récapitulatif :

Aspect clé	Résumé
Principale non-linéarité	Cinétique de réaction d'Arrhenius, dans laquelle la température affecte la vitesse de réaction de manière exponentielle.
Point de fonctionnement	État stationnaire sélectionné, défini par les conditions nominales de concentration, de température, de volume, de débit, d'alimentation et de refroidissement.
Méthode de linéarisation	Développement de Taylor du premier ordre des bilans non linéaires de matière et d'énergie.
Principales variables d'état	Concentration du réactif et température ; le volume du réacteur est inclus lorsque les dynamiques de niveau ou de débit sont importantes.
Modèle linéaire	Forme en variables d'écart : Δẋ = AΔx + BΔu + EΔd.
Entrées de contrôle	Puissance de refroidissement, débit du fluide caloporteur, débit de sortie ou autres variables manipulées.
Applications de l'analyse	Stabilité locale, réponse transitoire, constantes de temps, interactions entre procédés et conception du contrôle de la température.
Limitation principale	Le modèle est précis localement et ne doit pas remplacer le modèle non linéaire pour l'analyse des grands transitoires ou de l'emballement thermique.

Vous planifiez une installation pilote de CSTR à des fins d'enseignement, de recherche, de développement de procédés ou de validation industrielle ? Contactez LABPARK pour discuter de vos besoins. LABPARK fournit des installations pilotes d'opérations unitaires éducatives et professionnelles pour le génie chimique, les bioprocédés et les biotechnologies, ainsi que les applications environnementales et de traitement de l'eau. Notre équipe peut aider les universités, les instituts de recherche et les entreprises à sélectionner un système CSTR non isotherme approprié, à configurer les fonctions de mesure et de contrôle et à adapter l'installation pilote aux études portant sur la cinétique des réactions, le transfert thermique, le contrôle des procédés, la stabilité et le changement d'échelle. Contactez-nous pour développer une solution d'installation pilote pratique et fiable, qui favorise la formation pratique, les expériences reproductibles et la recherche avancée en contrôle des procédés.