La conduction thermique transitoire dans un élément chauffant électrique cylindrique est modélisée à l'aide de l'équation de la chaleur radiale, d'une génération interne de chaleur par effet Joule et d'une condition aux limites convective à la surface. Pour une longue tige pleine de rayon (R), présentant des propriétés constantes et une génération électrique uniforme, la température adimensionnelle est représentée par un profil radial stationnaire plus une série infinie de modes de fonctions de Bessel décroissants. Les taux de décroissance sont déterminés par des valeurs propres qui dépendent du nombre de Biot de surface.
L'idée essentielle consiste à séparer la température en un état stationnaire dû à la génération de chaleur et une correction transitoire. La symétrie cylindrique donne des fonctions propres spatiales (J_0(\lambda_n r)), tandis que la convection à la surface de la tige détermine les valeurs propres par la relation (\lambda_n J_1(\lambda_n)=\mathrm{Bi},J_0(\lambda_n)).
Définition du problème de chauffage cylindrique
Hypothèses physiques
Un modèle couramment utilisé pour une installation pilote suppose que :
- L'élément chauffant est un long cylindre, de sorte que les effets axiaux et d'extrémité sont négligeables.
- La génération de chaleur électrique est uniforme dans le solide.
- La conductivité thermique (k), la masse volumique (\rho) et la capacité calorifique (C_p) sont constantes.
- La température initiale est uniforme.
- La chaleur quitte la surface extérieure par convection vers le fluide de procédé environnant.
Ces hypothèses sont appropriées lorsque la longueur de la tige est très supérieure à son rayon et lorsque les propriétés électriques et matérielles ne varient pas de manière significative avec la température.
Génération électrique volumique
Le taux de génération interne de chaleur est noté (q'''), avec des unités de (\mathrm{W/m^3}). Pour un élément résistif, il peut être relié aux conditions électriques par
[ q'''=\sigma E^2 ]
où (\sigma) est la résistivité électrique et (E) le champ électrique. Des expressions équivalentes peuvent être formulées à l'aide du courant, de la tension, de la résistance et du volume de l'élément.
Équation énergétique gouvernante
Avec une symétrie cylindrique, l'équation thermique transitoire s'écrit
[ \rho C_p\frac{\partial T}{\partial t}
k\left[ \frac{1}{r}\frac{\partial}{\partial r} \left(r\frac{\partial T}{\partial r}\right) \right]+q''' ]
ou, en utilisant la diffusivité thermique (\alpha=k/(\rho C_p)),
[ \frac{\partial T}{\partial t}
\alpha\left[ \frac{1}{r}\frac{\partial}{\partial r} \left(r\frac{\partial T}{\partial r}\right) \right] +\frac{q'''}{\rho C_p}. ]
La coordonnée radiale vérifie (0\le r\le R).
Application des conditions aux limites et initiales
Condition sur l'axe central
La température doit rester finie sur l'axe central. La symétrie impose également un flux thermique radial nul :
[ \left.\frac{\partial T}{\partial r}\right|_{r=0}=0. ]
Cette condition élimine la fonction de Bessel de seconde espèce (Y_0), qui est singulière en (r=0).
Condition de surface convective
À la surface extérieure, la conduction dans la tige s'équilibre avec la convection vers le fluide environnant :
[ -k\left.\frac{\partial T}{\partial r}\right|_{r=R}
h,[T(R,t)-T_\infty], ]
où (h) est le coefficient de transfert thermique convectif et (T_\infty) la température du fluide de procédé en volume.
Il s'agit d'une condition aux limites de Robin, ou mixte, car elle relie à la fois la température et son gradient.
Condition initiale
Pour une tige préchauffée uniformément ou initialement froide,
[ T(r,0)=T_i. ]
La condition initiale détermine la contribution de chaque mode propre transitoire à la réponse ultérieure.
Formulation adimensionnelle
Variables adimensionnelles
Définissons le rayon, le temps et la température adimensionnels par
[ r^*=\frac{r}{R}, \qquad Fo=\frac{\alpha t}{R^2}, \qquad \Theta=\frac{k(T-T_\infty)}{q'''R^2}. ]
Le nombre de Fourier (Fo) mesure le temps écoulé par rapport au temps caractéristique de conduction interne.
Le nombre de Biot de surface est
[ \mathrm{Bi}=\frac{hR}{k}. ]
Il compare la résistance convective externe à la résistance conductrice interne.
Équation adimensionnelle
L'équation gouvernante devient
[ \frac{\partial \Theta}{\partial Fo}
\frac{1}{r^} \frac{\partial}{\partial r^} \left( r^\frac{\partial \Theta}{\partial r^} \right)+1. ]
Les conditions aux limites adimensionnelles sont
[ \left.\frac{\partial \Theta}{\partial r^}\right|_{r^=0}=0 ]
et
[ -\left.\frac{\partial \Theta}{\partial r^}\right|_{r^=1}
\mathrm{Bi},\Theta(1,Fo). ]
Construction de la solution thermique
Profil de température à l'état stationnaire
À l'état stationnaire, la dérivée temporelle est nulle :
[ \frac{1}{r^} \frac{d}{dr^} \left( r^\frac{d\Theta_{ss}}{dr^} \right)+1=0. ]
L'application des conditions sur l'axe central et à la surface convective donne
[ \Theta_{ss}(r^*)
\frac{1}{2\mathrm{Bi}} + \frac{1}{4}(1-r^{*2}). ]
Le premier terme représente l'élévation de température associée à la convection entre la surface et le fluide de procédé. Le second terme représente l'élévation radiale de température causée par la conduction de la chaleur du centre vers la surface.
Correction transitoire
Écrivons la solution complète sous la forme
[ \Theta(r^*,Fo)
\Theta_{ss}(r^)+\theta(r^,Fo). ]
La correction transitoire satisfait l'équation de diffusion homogène
[ \frac{\partial \theta}{\partial Fo}
\frac{1}{r^} \frac{\partial}{\partial r^} \left( r^\frac{\partial \theta}{\partial r^} \right). ]
La séparation des variables donne
[ \theta(r^*,Fo)
\sum_{n=1}^{\infty} A_n J_0(\lambda_n r^*) e^{-\lambda_n^2 Fo}. ]
Le champ de température complet est donc
[ \boxed{ \Theta(r^*,Fo)
\frac{1}{2\mathrm{Bi}} + \frac{1}{4}(1-r^{2}) + \sum_{n=1}^{\infty} A_n J_0(\lambda_n r^) e^{-\lambda_n^2 Fo} } ]
où (J_0) est la fonction de Bessel de première espèce et d'ordre zéro.
Équation aux valeurs propres
L'application de la condition convective à la limite à un mode (J_0(\lambda r^*)) donne
[ \lambda_n J_1(\lambda_n)
\mathrm{Bi},J_0(\lambda_n), ]
ou, de manière équivalente,
[ \lambda_n
\mathrm{Bi}, \frac{J_0(\lambda_n)}{J_1(\lambda_n)}. ]
Les racines positives (\lambda_n) sont déterminées numériquement. Chaque racine définit un mode radial de température.
Coefficients des modes
Les coefficients sont obtenus en développant l'écart initial par rapport à l'état stationnaire dans la base de fonctions de Bessel :
[ A_n
\frac{ \displaystyle \int_0^1 r^* \left[\Theta(r^,0)-\Theta_{ss}(r^)\right] J_0(\lambda_n r^),dr^ }{ \displaystyle \int_0^1 r^J_0^2(\lambda_n r^),dr^* }. ]
Le facteur de pondération (r^*) est requis par la géométrie cylindrique et la relation d'orthogonalité des fonctions propres.
Interprétation du modèle dans une installation pilote
Prédiction du temps de chauffage
Le mode le plus lent, associé à la plus petite valeur propre (\lambda_1), contrôle généralement la réponse aux temps longs :
[ \theta\sim A_1J_0(\lambda_1r^*)e^{-\lambda_1^2Fo}. ]
Cela permet aux ingénieurs d'estimer le temps nécessaire à l'élément pour atteindre sa température stationnaire avant d'effectuer un calcul complet en série.
Estimation des températures de surface et du centre
Le modèle fournit la température à n'importe quelle position radiale. En particulier :
- La température sur l'axe central s'obtient en posant (r^*=0), pour lequel (J_0(0)=1).
- La température de surface s'obtient en posant (r^*=1).
- La température maximale se trouve généralement au centre, car la chaleur doit se conduire vers l'extérieur à travers la tige.
Ces prédictions sont importantes lorsque le chauffage est immergé dans un fluide de procédé, intégré dans une paroi ou utilisé pour fournir de la chaleur à une opération unitaire en enceinte ou en conduite.
Relation entre le modèle et l'écoulement du fluide
Le nombre de Biot dépend du coefficient de convection (h). Dans une installation pilote, (h) peut varier avec la vitesse du fluide, l'agitation, la viscosité et le régime d'écoulement.
Ces effets sont couramment représentés par des corrélations faisant intervenir les nombres de Reynolds et de Nusselt. Une variation de la vitesse du fluide de procédé modifie donc (\mathrm{Bi}), les valeurs propres, la température de surface et le temps global de chauffage.
Évaluation de l'approximation à capacité concentrée
Lorsque
[ \mathrm{Bi}<0.1, ]
les gradients de température internes sont souvent suffisamment faibles pour que la tige puisse être approximée par une température uniforme. La solution complète en série de fonctions de Bessel devient nécessaire lorsque la résistance à la conduction interne est importante et que les températures au centre et à la surface diffèrent sensiblement.
Extension du modèle aux équipements réels
Éléments chauffants de longueur finie
Le modèle du long cylindre néglige les pertes de chaleur par les extrémités. Pour les éléments courts, la température dépend à la fois de (r) et de la coordonnée axiale (z), ce qui nécessite une formulation bidimensionnelle ou tridimensionnelle.
La solution implique alors des fonctions de Bessel radiales couplées à des fonctions propres axiales.
Propriétés dépendant de la température
La résistivité électrique, la conductivité thermique et la capacité calorifique peuvent varier avec la température. Dans ce cas, (q'''), (k) et (\alpha) ne sont plus constants, et la solution analytique en série peut devoir être remplacée par une intégration temporelle numérique.
Systèmes composites d'isolation et de paroi
Si l'élément chauffant est entouré d'une isolation, d'une gaine métallique ou de la paroi d'une cuve de procédé, chaque couche nécessite sa propre équation de conduction. À chaque interface, le modèle doit imposer
[ T_1=T_2 ]
et
[ -k_1\frac{\partial T_1}{\partial r}
-k_2\frac{\partial T_2}{\partial r}. ]
Ces conditions assurent la continuité de la température et du flux thermique à travers la structure composite.
Couplage avec le fluide de procédé
L'hypothèse d'une température limite constante (T_\infty) est appropriée lorsque le fluide environnant est bien mélangé ou se comporte comme un grand réservoir thermique. Si la température du fluide varie de manière significative, le modèle du chauffage doit être couplé à un bilan énergétique du fluide faisant intervenir la masse, la capacité calorifique, le débit, la surface d'échange thermique et le coefficient global de transfert thermique.
Comprendre les compromis
Solution analytique ou modèle numérique
La solution en série de fonctions de Bessel est rapide, transparente et utile pour les calculs de conception et les expériences pédagogiques. Elle devient moins pratique lorsque la géométrie, les propriétés, les taux de génération ou les conditions aux limites varient avec le temps ou la position.
Les méthodes des différences finies, des volumes finis ou des éléments finis prennent plus naturellement en charge ces complications, mais nécessitent une discrétisation numérique et une vérification.
Génération uniforme ou non uniforme
Le profil analytique ci-dessus suppose une génération de chaleur volumique uniforme. Une charge électrique non uniforme, une résistivité dépendant de la température ou un chauffage localisé modifie le terme source et, par conséquent, le profil stationnaire ainsi que la réponse transitoire.
Coefficient de convection constant ou conditions d'écoulement variables
L'utilisation d'un coefficient (h) constant simplifie le calcul, mais les conditions d'écoulement d'une installation pilote peuvent changer pendant le fonctionnement. La vitesse d'agitation, la viscosité du fluide, les changements de phase et l'encrassement peuvent modifier (h) ; les valeurs mesurées ou issues de corrélations doivent donc être vérifiées sur la plage de fonctionnement.
Température au centre ou température moyenne
La température au centre constitue souvent la variable critique pour la sécurité ou la qualité du matériau, mais elle n'est pas identique à la température moyenne volumique. Les bilans énergétiques du procédé peuvent nécessiter la température moyenne du solide, tandis que les limites imposées au matériau peuvent dépendre de la température maximale sur l'axe central.
Faire le choix approprié selon votre objectif
Le modèle approprié dépend de la précision requise, des mesures disponibles et des conditions de fonctionnement.
- Si votre objectif principal est d'estimer le temps de réponse du chauffage : utilisez les nombres de Fourier et de Biot adimensionnels avec les premiers modes de Bessel, en vérifiant si l'approximation à capacité concentrée est valide.
- Si votre objectif principal est de prédire le point le plus chaud : évaluez la solution radiale complète sur l'axe central, car la température maximale se trouve généralement en (r=0).
- Si votre objectif principal est de déterminer les performances de convection : utilisez des valeurs de (h) mesurées ou issues de corrélations pour calculer (\mathrm{Bi}), puis comparez les températures de surface et les vitesses de chauffage prédites aux données de l'installation pilote.
- Si votre objectif principal est de modéliser un élément court, isolé ou composite : utilisez un modèle numérique multidimensionnel ou multicouche avec continuité de la température et du flux thermique aux interfaces.
- Si votre objectif principal est la régulation du procédé : couplez le modèle de conduction cylindrique au bilan énergétique du fluide ou de la cuve au lieu de considérer la température du fluide environnant comme constante.
Un chauffage électrique cylindrique est donc modélisé en combinant la génération interne par effet Joule, la conduction radiale, les pertes convectives à la surface et les modes transitoires de fonctions de Bessel afin de prédire son évolution thermique pendant toute l'exploitation de l'installation pilote.
Tableau récapitulatif :
| Élément du modèle | Objectif | Relation clé |
|---|---|---|
| Génération interne de chaleur | Représente le chauffage électrique résistif | (q'''=\sigma E^2) |
| Équation thermique radiale | Décrit la température transitoire à l'intérieur du cylindre | (\rho C_p\frac{\partial T}{\partial t}=k\nabla_r^2T+q''') |
| Convection de surface | Modélise le transfert thermique vers le fluide de procédé | (-k\frac{\partial T}{\partial r}=h(T-T_\infty)) |
| Nombre de Fourier | Mesure la progression de la diffusion thermique interne | (Fo=\alpha t/R^2) |
| Nombre de Biot | Indique l'importance des gradients de température internes | (Bi=hR/k) |
| Fonctions propres de Bessel | Représentent les modes de température transitoire cylindriques | (J_0(\lambda_n r^*)) |
| Équation aux valeurs propres | Détermine les taux de décroissance modale | (\lambda_nJ_1(\lambda_n)=BiJ_0(\lambda_n)) |
| Vérification de la capacité concentrée | Identifie les situations où un modèle à température uniforme peut être suffisant | Généralement valide lorsque (Bi<0.1) |
| Modèle numérique étendu | Traite les longueurs finies, les couches composites, les propriétés variables ou les conditions variables du fluide | Utiliser les méthodes des différences finies, des volumes finis ou des éléments finis |
Concevez le système de chauffage adapté à votre installation pilote
Besoin d'un contrôle fiable de la température pour une application de génie chimique, de bioprocédés, de biotechnologie ou de traitement de l'environnement et de l'eau ? Contactez LABPARK pour discuter de vos besoins.
LABPARK fournit des installations pilotes d'opérations unitaires à vocation pédagogique et professionnelle pour les universités, les instituts de recherche et les entreprises. Notre équipe peut vous aider à sélectionner et à configurer des systèmes tenant compte du temps de réponse du chauffage, de la température sur l'axe central, des performances de transfert thermique, des conditions d'écoulement du fluide, de l'isolation, des parois composites et des exigences de régulation du procédé. Nous proposons des solutions pratiques pour les installations pilotes destinées à l'enseignement, à la recherche, au développement des procédés et à la formation industrielle, afin de vous aider à obtenir un fonctionnement reproductible et des données expérimentales pertinentes.
Produits associés
Les gens demandent aussi
- Pourquoi appliquer la correction de la DMLT dans les pilotes à tubes et calandre et comment la déterminer
- Pourquoi régler l'espacement des chicanes dans les pilotes d'échangeurs de chaleur à calandre et tubes ? Maîtrisez le transfert thermique et la perte de charge.
- Comment estimer la perte de charge côté coque dans une installation pilote d'échangeur de chaleur à calandre et tubes : explication des paramètres clés
- Comment le facteur d'encrassement (Rd) est-il évalué ? Principaux aperçus de l'usine pilote pour les étudiants
- Comment le facteur d'encrassement est-il démontré à l'aide de pilotes à calandre et tubes ? Guide de laboratoire pratique
