Les fonctions de transfert sont dérivées en linéarisant le bilan énergétique transitoire du chauffage à réservoir agité autour d'un point de fonctionnement stable. Pour un réservoir parfaitement mélangé avec une masse de liquide constante, une température d'entrée (T_i), une température de vapeur (T_{st}), un débit (F), une surface d'échange thermique (A) et un coefficient global de transfert thermique (U), le modèle résultant est du premier ordre. Il produit des fonctions de transfert distinctes entre la température du fluide d'entrée et la température de la vapeur vers la température du réservoir.

La température du réservoir réagit aux deux perturbations par le biais des mêmes dynamiques de stockage thermique et d'évacuation de chaleur. Le chemin de la température d'entrée est régi par le remplacement du fluide, tandis que le chemin de la température de la vapeur est régi par la capacité de transfert thermique du serpentin (UA).

Commencer par le bilan énergétique instationnaire

Définir les hypothèses sur le réservoir et l'écoulement

Supposons que le réservoir soit parfaitement agité, de sorte que la température de sortie soit égale à la température du réservoir (T). Supposons également un volume de liquide, une densité, une capacité thermique, un débit et des propriétés de transfert thermique constants.

Le réservoir contient une masse de liquide (M), et le fluide de procédé a une capacité thermique (C_p). La chaleur entre par le serpentin de vapeur selon

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

où (U) est le coefficient global de transfert thermique et (A) est la surface effective de transfert thermique.

Écrire le terme d'accumulation

L'énergie thermique stockée dans le réservoir change au taux de

[ MC_p\frac{dT}{dt} ]

Ce terme représente le comportement dynamique : une masse de liquide ou une capacité thermique plus importante nécessite plus d'énergie pour modifier la température du réservoir.

Écrire les termes d'écoulement et de chauffage

L'énergie nette entrant par le flux de procédé est

[ FC_p(T_i-T) ]

La première partie, (FC_pT_i), représente l'énergie entrant avec l'alimentation. La seconde partie, (FC_pT), représente l'énergie sortant avec la sortie bien mélangée.

L'ajout de la puissance du serpentin de vapeur donne le bilan complet :

[ MC_p\frac{dT}{dt}

FC_p(T_i-T)+UA(T_{st}-T) ]

Réorganiser le modèle dynamique

Développer les termes de transfert thermique

Le développement du côté droit donne

[ MC_p\frac{dT}{dt}

FC_pT_i-FC_pT+UAT_{st}-UAT ]

Le regroupement des termes de température du réservoir produit

[ MC_p\frac{dT}{dt}+(FC_p+UA)T

FC_pT_i+UAT_{st} ]

Normaliser l'équation

La division par (MC_p) donne

[ \frac{dT}{dt} + \left( \frac{F}{M}+\frac{UA}{MC_p} \right)T

\frac{F}{M}T_i + \frac{UA}{MC_p}T_{st} ]

Définissons

[ a=\frac{F}{M}+\frac{UA}{MC_p} ]

[ b_i=\frac{F}{M} ]

[ b_{st}=\frac{UA}{MC_p} ]

Le modèle devient

[ \frac{dT}{dt}+aT=b_iT_i+b_{st}T_{st} ]

Le coefficient (a) détermine la vitesse à laquelle le réservoir atteint un nouvel équilibre. Son inverse,

[ \tau=\frac{1}{a}

\frac{MC_p}{FC_p+UA} ]

est la constante de temps thermique effective.

Convertir en variables d'écart

Définir le point de fonctionnement stable

Soient les valeurs initiales à l'état stationnaire (T_s), (T_{i,s}) et (T_{st,s}). À l'état stationnaire,

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

donc le bilan à l'état stationnaire est

[ aT_s=b_iT_{i,s}+b_{st}T_{st,s} ]

Ce point de fonctionnement est la référence utilisée pour étudier les perturbations.

Définir les écarts de température

Définissons les variables d'écart comme

[ T'=T-T_s ]

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

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

Comme les paramètres du système sont supposés constants, la soustraction de l'équation d'état stationnaire de l'équation dynamique donne

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

b_iT_i'+b_{st}T_{st}' ]

Cette étape élimine les termes constants du point de fonctionnement et ne laisse que la réponse aux changements des deux entrées.

Appliquer la transformée de Laplace

Transformer l'équation différentielle

Supposons un écart initial nul, ce qui signifie que le réservoir commence à sa température d'état stationnaire. La prise de la transformée de Laplace donne

[ sT'(s)+aT'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

La factorisation de la transformée de la température du réservoir donne

[ (s+a)T'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

Par conséquent,

[ T'(s)

\frac{b_i}{s+a}T_i'(s) + \frac{b_{st}}{s+a}T_{st}'(s) ]

Identifier les deux fonctions de transfert

La fonction de transfert de la température du fluide d'entrée vers la température du réservoir est

[ G_i(s)

\frac{T'(s)}{T_i'(s)}

\frac{b_i}{s+a}

\frac{F/M}{s+F/M+UA/(MC_p)} ]

avec l'écart de température de la vapeur maintenu à zéro.

La fonction de transfert de la température de la vapeur vers la température du réservoir est

[ G_{st}(s)

\frac{T'(s)}{T_{st}'(s)}

\frac{b_{st}}{s+a}

\frac{UA/(MC_p)}{s+F/M+UA/(MC_p)} ]

avec l'écart de température du fluide d'entrée maintenu à zéro.

La relation entrée-sortie complète est donc

[ T'(s)=G_i(s)T_i'(s)+G_{st}(s)T_{st}'(s) ]

Interpréter physiquement le modèle

Réponse à la température d'entrée

Le chemin de la température d'entrée a pour numérateur (F/M). Un débit plus élevé fait en sorte que les perturbations d'entrée influencent plus fortement le réservoir et réduit le temps de séjour du liquide.

À l'état stationnaire, le gain de température d'entrée est

[ K_i=G_i(0)=\frac{b_i}{a} ]

[ K_i

\frac{F/M}{F/M+UA/(MC_p)} ]

Réponse à la température de la vapeur

Le chemin de la température de la vapeur a pour numérateur (UA/(MC_p)). L'augmentation de la surface de transfert thermique ou du coefficient global de transfert thermique augmente l'influence des changements de température de la vapeur sur le réservoir.

Son gain à l'état stationnaire est

[ K_{st}=G_{st}(0)=\frac{b_{st}}{a} ]

[ K_{st}

\frac{UA/(MC_p)}{F/M+UA/(MC_p)} ]

Pour ce modèle simplifié,

[ K_i+K_{st}=1 ]

Cela reflète le fait que la température stable du réservoir est un équilibre pondéré entre la température du fluide d'entrée et la température de la vapeur.

Constante de temps et pôle

Les deux fonctions de transfert ont le même dénominateur et donc le même pôle :

[ s=-a ]

La constante de temps commune est

[ \tau=\frac{1}{a} ]

En utilisant (\tau), les fonctions de transfert peuvent également être écrites comme

[ G_i(s)=\frac{K_i}{\tau s+1} ]

[ G_{st}(s)=\frac{K_{st}}{\tau s+1} ]

Le pôle commun signifie que les deux entrées sont filtrées par le même stockage thermique du réservoir. Leurs effets à l'état stationnaire diffèrent, mais leur vitesse de réponse de base du premier ordre est la même sous les hypothèses énoncées.

Prédire les réponses échelon

Changement échelon de la température d'entrée

Pour un changement échelon d'amplitude (\Delta T_i) de la température d'entrée, avec une température de vapeur inchangée,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) ]

La température du réservoir atteint sa nouvelle valeur d'état stationnaire de manière exponentielle.

Changement échelon de la température de la vapeur

Pour un échelon de température de vapeur d'amplitude (\Delta T_{st}), avec une température d'entrée inchangée,

[ T'(t)

K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

La forme est identique, mais le gain est déterminé par la force relative du transfert thermique du serpentin.

Changements simultanés

Si les deux températures changent en même temps,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) + K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

Pour des entrées variables dans le temps plus générales, le modèle dans le domaine de Laplace est utilisé directement, ou l'équation différentielle originale peut être résolue numériquement et comparée aux données des capteurs de l'installation pilote.

Connecter les paramètres aux mesures de l'installation pilote

Estimer la capacité de transfert thermique

La quantité (UA) contrôle l'efficacité avec laquelle les changements de température de la vapeur sont transmis au réservoir. Elle peut être estimée à partir de

[ Q=UA\Delta T_m ]

où (\Delta T_m) est la différence de température moyenne appropriée, souvent la différence de température logarithmique moyenne pour les calculs de type échangeur.

La valeur mesurée de (U) intègre la convection, la conduction de paroi et les résistances d'encrassement.

Évaluer l'effet de l'agitation

L'agitation améliore le mélange et peut augmenter le coefficient de transfert thermique côté fluide. Elle soutient également l'hypothèse de mélange parfait utilisée dans le modèle.

Un mauvais mélange peut entraîner une déviation des réponses mesurées par rapport au modèle à température unique, car différentes régions du réservoir peuvent réagir à des vitesses différentes.

Utiliser les paramètres adimensionnels avec précaution

Le rapport

[ B=\frac{UA}{FC_p} ]

compare la capacité de transfert thermique du serpentin avec le débit de capacité thermique du flux de procédé.

Un grand (B) indique que le chauffage par vapeur est fort par rapport aux effets du débit de fluide. Un petit (B) indique que le flux d'entrée a une plus grande influence sur la température de sortie.

Comprendre les compromis

Un modèle du premier ordre est une approximation

Le modèle suppose (U), (A), (F), (M) et (C_p) constants. Dans une installation pilote réelle, ces paramètres peuvent varier avec le débit, la température, la viscosité, l'encrassement et la vitesse de l'agitateur.

Les fonctions de transfert sont donc plus fiables près du point de fonctionnement stable utilisé pour la linéarisation.

La température de la vapeur peut ne pas être égale à la puissance thermique

Le modèle traite la température de la vapeur comme une entrée et utilise

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

Ceci est approprié lorsque la vapeur se condense à une température approximativement constante et que la dynamique du serpentin est rapide par rapport à la dynamique du réservoir.

Si la pression de la vapeur, la position de la vanne, l'élimination des condensats ou la température de la paroi du serpentin ont une dynamique significative, la température de la vapeur seule peut ne pas représenter adéquatement l'entrée de chauffage.

La dynamique du transport et des capteurs peut être importante

Le bilan idéal n'a pas de retard explicite entre un changement de température d'entrée et son effet sur le réservoir. La tuyauterie réelle, les emplacements d'échantillonnage, les puits thermométriques et les capteurs de température peuvent introduire des retards de transport et des constantes de temps supplémentaires.

Ces effets doivent être inclus lorsque les fonctions de transfert ajustées ne correspondent pas aux données expérimentales des tests échelon.

La linéarisation limite la plage de validité

Le modèle à variables d'écart suppose des changements suffisamment petits pour que les paramètres du système restent approximativement constants. De grands changements de température ou de débit peuvent nécessiter un bilan énergétique non linéaire ou une intégration numérique.

Une validation fiable vérifie la condition initiale, l'état stationnaire à long terme et si la courbe transitoire prédite concorde avec les données mesurées.

Faire le bon choix pour votre objectif

Utilisez le modèle dérivé du bilan comme point de départ pour l'analyse et la conception du contrôle.

Si votre objectif principal est les perturbations de la température d'entrée : Utilisez (G_i(s)), en mettant l'accent sur le rapport débit/retenue (F/M) et le gain de réponse échelon d'entrée (K_i).
Si votre objectif principal est le contrôle de la température de la vapeur : Utilisez (G_{st}(s)), en mettant l'accent sur (UA/(MC_p)), la surface de transfert thermique, l'agitation et le coefficient global (U).
Si votre objectif principal est la validation de l'installation pilote : Effectuez des tests échelon, estimez les gains et la constante de temps, et comparez la réponse mesurée avec les prédictions du premier ordre.
Si votre objectif principal est la simulation haute fidélité : Conservez le bilan énergétique non linéaire et ajoutez la dynamique du transport, du serpentin, de la vanne et du capteur là où les données montrent que le modèle du premier ordre est insuffisant.

Une fonction de transfert est la conséquence dynamique compacte d'un bilan énergétique physique, reliant les changements de procédé mesurables à la réponse de la température du réservoir.

Tableau récapitulatif :

Élément clé	Résultat	Signification pratique
Modèle dynamique	`MCp(dT/dt) = FCp(Ti - T) + UA(Tst - T)`	Décrit l'accumulation de chaleur, le débit d'entrée et le chauffage par vapeur
Fonction de transfert de la température d'entrée	`Gi(s) = (F/M) / [s + F/M + UA/(MCp)]`	Prédit la réponse du réservoir aux changements de température du fluide d'entrée
Fonction de transfert de la température de la vapeur	`Gst(s) = [UA/(MCp)] / [s + F/M + UA/(MCp)]`	Prédit la réponse du réservoir aux changements de température de la vapeur
Gain d'entrée à l'état stationnaire	`Ki = (F/M) / [F/M + UA/(MCp)]`	Mesure l'influence de la température d'entrée sur celle du réservoir
Gain de vapeur à l'état stationnaire	`Kst = [UA/(MCp)] / [F/M + UA/(MCp)]`	Mesure l'influence chauffante de la température de la vapeur
Constante de temps commune	`tau = MCp / (FCp + UA)`	Indique la rapidité avec laquelle le réservoir atteint une nouvelle température
Facteurs de validation clés	`U`, `A`, `F`, `M`, `Cp`, mélange, retard et dynamique des capteurs	Détermine la précision du modèle par rapport aux données de l'installation pilote

Modélisez et validez votre installation pilote avec LABPARK

Besoin de données dynamiques fiables pour la recherche en génie chimique ou en contrôle de procédé ? Contactez LABPARK pour discuter d'une installation pilote d'opérations unitaires éducative ou professionnelle conçue pour votre application.

LABPARK soutient les universités, les instituts de recherche et les entreprises avec des systèmes à l'échelle pilote pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que le traitement de l'eau et de l'environnement. Notre chauffage à réservoir agité et nos équipements d'opérations unitaires associés vous aident à :

Démontrer les bilans énergétiques, les fonctions de transfert, la dynamique des procédés et les stratégies de contrôle.
Exécuter des tests échelon reproductibles sur la température d'entrée et la température de la vapeur.
Évaluer les coefficients de transfert thermique, le temps de séjour, l'agitation et la réponse thermique.
Générer des données pratiques pour l'enseignement, la recherche, le développement de procédés et la formation des opérateurs.
Construire une plateforme d'installation pilote configurable alignée sur vos besoins de laboratoire, de programme ou de production.

Indiquez-nous votre capacité requise, la température de fonctionnement, l'instrumentation et vos objectifs d'enseignement ou de recherche. LABPARK peut vous aider à identifier la configuration d'installation pilote appropriée et fournir un soutien technique pratique, de la spécification à l'application.