Les bilans énergétiques dynamiques pour un réservoir agité chauffant sont établis en soustrayant le bilan en régime permanent du bilan transitoire. Pour un réservoir à volume constant, parfaitement mélangé, avec un débit et des propriétés constants, le modèle de température devient :

[ \frac{dT'}{dt}+aT'=\frac{1}{\tau}T_i'+KT_{st}' ]

où (T'), (T_i') et (T_{st}') sont les écarts de la température du réservoir, de la température d'entrée et de l'apport de chauffage par rapport à leurs valeurs nominales en régime permanent. Cette forme isole uniquement les changements dynamiques qui importent pour l'analyse de la réponse et la conception du contrôle.

Point clé : Définissez chaque variable comme son écart par rapport à un point de fonctionnement en régime permanent, substituez ces définitions dans le bilan énergétique dynamique et utilisez l'équation en régime permanent pour annuler les termes constants. Le modèle d'écart résultant décrit le comportement thermique transitoire sans conserver les décalages du point de fonctionnement nominal.

Commencer par le bilan énergétique transitoire

Définir le volume de contrôle et les hypothèses

Pour un réservoir agité chauffant parfaitement mélangé, la température du liquide est supposée uniforme dans tout le réservoir. Un modèle d'unité pilote courant suppose un volume de liquide constant (V), un débit volumétrique constant (F), une masse volumique constante (\rho) et une capacité thermique constante (c_p).

Sous ces hypothèses, l'énergie cinétique, l'énergie potentielle, le travail sur l'arbre et les effets de changement de volume sont négligés. Les termes dominants sont l'accumulation d'énergie, le flux d'enthalpie d'entrée et de sortie, et la chaleur fournie de l'extérieur.

Écrire le bilan énergétique dimensionnel

Pour un réservoir de liquide sans réaction avec apport de chaleur (Q),

[ \rho c_p V\frac{dT}{dt}

\rho c_p F(T_i-T)+Q ]

où :

(T) est la température du liquide dans le réservoir,
(T_i) est la température d'entrée,
(F) est le débit volumétrique,
(V) est le volume de liquide,
(Q) est la puissance de chauffage.

La température de sortie est (T) car le réservoir est parfaitement mélangé.

Mettre le bilan sous forme dynamique standard

En divisant par (\rho c_p V), on obtient

[ \frac{dT}{dt} + \frac{F}{V}T

\frac{F}{V}T_i + \frac{Q}{\rho c_p V} ]

Le coefficient (F/V) est l'inverse du temps de séjour. En définissant

[ \tau=\frac{V}{F} ]

on obtient

[ \frac{dT}{dt} + \frac{1}{\tau}T

\frac{1}{\tau}T_i + \frac{Q}{\rho c_p V} ]

Si l'apport du chauffage est représenté par une variable (T_{st}), le terme d'apport de chaleur peut être écrit comme (KT_{st}). Le modèle résultant est couramment exprimé sous la forme

[ \frac{dT}{dt}+aT=\frac{1}{\tau}T_i+KT_{st} ]

où (a), (1/\tau) et (K) regroupent les paramètres pertinents de débit, volume, capacité thermique et transfert thermique.

Convertir le modèle en variables d'écart

Définir le point de fonctionnement nominal

Soient les conditions de fonctionnement en régime permanent

[ T_s,\qquad T_{i,s},\qquad T_{st,s} ]

En régime permanent, la température du réservoir ne change pas, donc

[ \frac{dT_s}{dt}=0 ]

L'équation en régime permanent est donc

[ aT_s=\frac{1}{\tau}T_{i,s}+KT_{st,s} ]

Cette équation décrit le point de fonctionnement autour duquel le modèle dynamique est développé.

Définir les écarts par rapport au régime permanent

Les variables d'écart sont

[ T'=T-T_s ]

[ T_i'=T_i-T_{i,s} ]

[ T_{st}'=T_{st}-T_{st,s} ]

De manière équivalente,

[ T=T_s+T' ]

[ T_i=T_{i,s}+T_i' ]

[ T_{st}=T_{st,s}+T_{st}' ]

Comme (T_s) est constant,

[ \frac{dT}{dt}=\frac{dT'}{dt} ]

Substituer dans le bilan dynamique

En commençant par

[ \frac{dT}{dt}+aT=\frac{1}{\tau}T_i+KT_{st} ]

la substitution donne

[ \frac{dT'}{dt} +a(T_s+T')

\frac{1}{\tau}(T_{i,s}+T_i') + K(T_{st,s}+T_{st}') ]

En réarrangeant,

[ \frac{dT'}{dt}+aT' + aT_s

\frac{1}{\tau}T_{i,s} + \frac{1}{\tau}T_i' + KT_{st,s} + KT_{st}' ]

Les termes constants satisfont la relation de régime permanent

[ aT_s=\frac{1}{\tau}T_{i,s}+KT_{st,s} ]

et s'annulent donc. L'équation en variables d'écart est

[ \boxed{ \frac{dT'}{dt}+aT'

\frac{1}{\tau}T_i' + KT_{st}' } ]

C'est le bilan énergétique dynamique souhaité sous forme d'écart.

Interpréter le modèle d'écart

Séparer les entrées manipulées et les perturbations

L'apport de chauffage (T_{st}'), ou de manière équivalente (Q'), est généralement la variable manipulée d'écart. L'écart de température d'entrée (T_i') est couramment traité comme une perturbation.

L'écart de température du réservoir (T') est la sortie du processus. L'équation décrit donc comment les changements du chauffage et les perturbations à l'entrée affectent la température de sortie ou du réservoir.

Utiliser la forme de puissance thermique physique si approprié

Si l'entrée manipulée est la puissance de chauffage plutôt qu'une température de fluide caloporteur, définissez

[ Q'=Q-Q_s ]

Le modèle d'écart à propriétés constantes devient

[ \rho c_p V\frac{dT'}{dt}

\rho c_p F(T_i'-T')+Q' ]

[ \boxed{ \frac{dT'}{dt} + \frac{F}{V}T'

\frac{F}{V}T_i' + \frac{Q'}{\rho c_p V} } ]

Cette forme est souvent plus transparente pour les calculs d'unités pilotes car (Q') a la signification physique directe d'un changement de puissance de chauffage.

Reconnaître la variable d'état

Pour un réchauffeur à volume constant, (T') est la principale variable d'état thermique. Si le niveau de liquide ou le volume change de manière significative, le bilan massique doit également être inclus :

[ A\frac{dh}{dt}=F_i-F ]

et le bilan énergétique peut être écrit comme

[ Ah\frac{dT}{dt}

F_i(T_i-T)+\frac{Q}{\rho c_p} ]

Dans ce cas, la dynamique de la température est couplée à la dynamique du niveau, et une seule équation de température du premier ordre peut ne plus être suffisante.

Obtenir les fonctions de transfert pour l'analyse dynamique

Appliquer la transformée de Laplace

En prenant la transformée de Laplace de

[ \frac{dT'}{dt}+aT'

\frac{1}{\tau}T_i'+KT_{st}' ]

on obtient

[ s\bar T'(s)-T'(0)+a\bar T'(s)

\frac{1}{\tau}\bar T_i'(s) + K\bar T_{st}'(s) ]

Avec un écart de température initial nul,

[ T'(0)=0 ]

donc

[ (s+a)\bar T'(s)

\frac{1}{\tau}\bar T_i'(s) + K\bar T_{st}'(s) ]

Dériver la fonction de transfert de la température d'entrée

En fixant l'écart du chauffage à zéro, on obtient

[ \boxed{ \frac{\bar T'(s)}{\bar T_i'(s)}

\frac{1/\tau}{s+a} } ]

Cette fonction de transfert décrit la réponse de la température du réservoir aux perturbations de la température d'entrée.

Dériver la fonction de transfert de l'apport de chauffage

En fixant l'écart de la température d'entrée à zéro, on obtient

[ \boxed{ \frac{\bar T'(s)}{\bar T_{st}'(s)}

\frac{K}{s+a} } ]

Si la puissance de chauffage (Q') est utilisée directement, la fonction de transfert correspondante est

[ \boxed{ \frac{\bar T'(s)}{\bar Q'(s)}

\frac{1}{\rho c_p V} \frac{1}{s+F/V} } ]

Ces relations supportent l'analyse de réponse par échelon, impulsion et rampe, et fournissent une base pour la conception de contrôleurs à rétroaction.

Étendre le bilan lorsque le réchauffeur est une double enveloppe ou un serpentin

Inclure la force motrice du transfert thermique

Si la chaleur est fournie par une double enveloppe ou un serpentin, la puissance peut être modélisée comme

[ Q=UA(T_{st}-T) ]

où (U) est le coefficient global de transfert thermique et (A) est la surface d'échange.

Le bilan énergétique transitoire devient

[ \rho c_p V\frac{dT}{dt}

\rho c_p F(T_i-T) + UA(T_{st}-T) ]

En réarrangeant,

[ \frac{dT}{dt} + \left( \frac{F}{V}+\frac{UA}{\rho c_p V} \right)T

\frac{F}{V}T_i + \frac{UA}{\rho c_p V}T_{st} ]

Ainsi, dans cette représentation,

[ a=\frac{F}{V}+\frac{UA}{\rho c_p V} ]

[ K=\frac{UA}{\rho c_p V} ]

L'équation d'écart conserve la même structure :

[ \frac{dT'}{dt}+aT'

\frac{F}{V}T_i' + \frac{UA}{\rho c_p V}T_{st}' ]

Prendre en compte les effets de réaction si nécessaire

Pour un réservoir agité avec réaction, le bilan énergétique doit également inclure la chaleur générée ou consommée par les réactions. Un bilan d'enthalpie général peut être représenté comme

[ V\frac{d}{dt} \left(\sum_i C_i\bar H_i\right)

F\left[ \sum_i(C_i\bar H_i)_{\text{in}}

\sum_i C_i\bar H_i \right] +\dot Q ]

La chaleur de réaction peut rendre le bilan non linéaire car les vitesses de réaction et les chaleurs peuvent dépendre de la température et de la composition.

Un modèle d'écart linéaire est alors obtenu en linéarisant le bilan non linéaire autour du régime permanent choisi. Les coefficients dans l'équation d'écart représentent les sensibilités locales à ce point de fonctionnement.

Comprendre les compromis

Les modèles d'écart sont des modèles locaux

L'équation d'écart est généralement précise près du régime permanent sélectionné. De grands changements de température, des propriétés physiques changeantes, un débit variable, des changements de phase ou des vitesses de réaction fortement dépendantes de la température peuvent invalider l'approximation à coefficients constants.

Un nouveau point de fonctionnement peut nécessiter une nouvelle linéarisation et de nouveaux paramètres de modèle.

Le régime permanent doit être cohérent

L'annulation des termes constants ne fonctionne que lorsque (T_s), (T_{i,s}) et (T_{st,s}) satisfont le bilan énergétique réel en régime permanent. Si les valeurs de base sont incohérentes, des termes constants résiduels subsistent dans l'équation d'écart.

C'est une source courante d'erreurs de modélisation apparentes.

Une hypothèse de volume constant peut masquer les effets de niveau

Si les débits d'entrée et de sortie diffèrent, le niveau du réservoir change. Comme le volume de liquide détermine la capacité thermique, la dynamique du niveau altère la réponse en température.

Pour les unités pilotes à niveau contrôlé, le modèle à volume constant est souvent approprié. Cependant, lors des transitoires de niveau, les bilans de masse et d'énergie doivent être conservés.

La température du réchauffeur et la puissance de chauffage ne sont pas des entrées identiques

L'utilisation de (T_{st}') suppose que la température du fluide caloporteur est l'entrée sélectionnée et que la relation de transfert thermique a été intégrée dans (K). L'utilisation de (Q') traite la puissance thermique comme la variable manipulée directe.

Le choix correct dépend de l'actionneur physique et des mesures disponibles sur l'unité pilote.

Comment appliquer cela à votre projet

La formulation appropriée dépend de si l'unité pilote est modélisée pour une analyse de réponse de base, une conception de contrôle ou une simulation détaillée de l'équipement.

Si votre objectif principal est la réponse en température aux perturbations d'entrée : Utilisez (\frac{dT'}{dt}+aT'=\frac{1}{\tau}T_i'+KT_{st}'), en traitant (T_i') comme l'entrée de perturbation.
Si votre objectif principal est le contrôle du réchauffeur : Utilisez l'écart de puissance de chauffage (Q') et écrivez (\frac{dT'}{dt}+\frac{F}{V}T'=\frac{F}{V}T_i'+\frac{Q'}{\rho c_pV}).
Si votre objectif principal est la conception de la double enveloppe ou du serpentin : Conservez (UA) explicitement afin que le modèle montre comment la surface et le coefficient de transfert thermique affectent à la fois le gain et le temps de réponse.
Si votre objectif principal est le contrôle combiné du niveau et de la température : Formulez des bilans de masse et d'énergie d'écart couplés avec (h') et (T') comme variables d'état.
Si votre objectif principal est une unité pilote avec réaction : Commencez par le bilan d'enthalpie complet et linéarisez les termes de réaction, de transfert thermique et de propriétés autour du point de fonctionnement.

Un modèle d'écart correctement formulé transforme le problème du point de fonctionnement du réchauffeur à réservoir agité en une description claire et traitable de la manière dont les perturbations et les actions de contrôle modifient la température.

Tableau récapitulatif :

Aspect clé	Résumé
Hypothèse principale	Volume, débit, masse volumique et capacité thermique constants dans un réservoir parfaitement mélangé
Modèle dynamique	`dT/dt + aT = (1/tau)Ti + KTst`
Variables d'écart	`T' = T - Ts`, `Ti' = Ti - Ti,s`, et `Tst' = Tst - Tst,s`
Équation d'écart	`dT'/dt + aT' = (1/tau)Ti' + KTst'`
Perturbation courante	Écart de température d'entrée, `Ti'`
Variable manipulée courante	Température du fluide caloporteur, `Tst'`, ou écart de puissance de chauffage, `Q'`
Modèle de double enveloppe ou serpentin	Inclure `UA` pour représenter les effets de transfert thermique
Limites du modèle	Les grands changements de fonctionnement, les propriétés variables, les changements de niveau et les effets de réaction peuvent nécessiter un modèle non linéaire ou couplé

Prêt à améliorer la formation au contrôle de la température, la recherche ou le développement de procédés ? Contactez LABPARK pour discuter d'une unité pilote d'opérations unitaires éducative ou professionnelle adaptée à vos besoins. LABPARK fournit des systèmes d'unités pilotes pratiques pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que les applications de traitement de l'eau et de l'environnement. Les universités, les instituts de recherche et les entreprises peuvent utiliser notre équipement pour une formation pratique, l'analyse de procédés, le développement de systèmes de contrôle, des études de mise à l'échelle et la formation technique, avec des configurations adaptées à vos objectifs de procédé et aux exigences de votre laboratoire.