Utilisez un bilan thermique à frontière mobile : l'approximation de l'état pseudo-stationnaire traite la conduction thermique à travers la couche solide existante comme instantanée par rapport au mouvement plus lent de l'interface solide-liquide. Le taux de croissance du solide est alors obtenu en égalant la chaleur latente libérée à l'interface à la chaleur transférée à travers le solide, la paroi du tube et le film de refroidissement interne. À mesure que la couche s'épaissit, sa résistance thermique augmente, et le taux de croissance diminue.

Le calcul pratique est un bilan de Stefan basé sur la résistance : déterminez le flux thermique de la température de changement de phase à la température globale du liquide de refroidissement, puis divisez ce flux par la chaleur latente nécessaire pour créer du solide supplémentaire.

Ce que suppose l'approximation de l'état pseudo-stationnaire

Les profils thermiques s'ajustent plus rapidement que le mouvement de l'interface

Le champ de température dans le solide est approximé comme stationnaire à chaque instant, bien que la frontière solide-liquide soit en mouvement. Cela est raisonnable lorsque la diffusion thermique à travers le solide est beaucoup plus rapide que le temps nécessaire pour que l'épaisseur de la couche change de manière significative.

L'approximation ne signifie pas que le processus est réellement stationnaire. L'épaisseur du solide et la résistance au transfert thermique varient toujours avec le temps.

La chaleur latente domine le bilan à frontière mobile

Pour la congélation ou la cristallisation, le principal besoin énergétique associé au mouvement de l'interface est la chaleur latente de fusion ou de cristallisation. L'accumulation de chaleur sensible dans le solide et le liquide est négligée ou traitée comme une correction secondaire.

Le bilan interfacial est donc basé sur

[ \text{chaleur évacuée}=\text{chaleur latente libérée par le nouveau solide}. ]

L'interface reste proche de la température de changement de phase

L'interface solide-liquide est généralement associée à la température de changement de phase à l'équilibre, (T_m), ou à une température de cristallisation effective lorsque les effets de concentration ou de pression sont importants.

Cela fournit la force motrice thermique pour évacuer la chaleur à travers l'ensemble solide et tube.

Construction du modèle de résistance thermique

Définir la géométrie du tube et du solide

Soit

(R_i) le rayon intérieur du tube,
(R_t) le rayon extérieur du tube,
(R_s(t)) le rayon de l'interface solide-liquide en mouvement,
(k_w) la conductivité thermique de la paroi du tube,
(k_s) la conductivité thermique de la couche solide,
(h_i) le coefficient de transfert thermique côté liquide de refroidissement interne,
(T_c) la température globale du liquide de refroidissement,
(\rho_s) la densité du solide,
(L) la chaleur latente de solidification ou de cristallisation.

Pour un dépôt croissant vers l'extérieur depuis la surface du tube, le solide occupe la région entre (R_t) et (R_s), et son épaisseur est

[ \delta(t)=R_s(t)-R_t. ]

Exprimer la résistance par unité de longueur de tube

Pour un tube cylindrique, les résistances du film de refroidissement, de la paroi du tube et de la couche solide par unité de longueur sont

[ \mathcal{R}_i'=\frac{1}{h_i,2\pi R_i}, ]

[ \mathcal{R}_w'=\frac{\ln(R_t/R_i)}{2\pi k_w}, ]

[ \mathcal{R}_s'=\frac{\ln(R_s/R_t)}{2\pi k_s}. ]

Des résistances supplémentaires peuvent être incluses lorsque cela est justifié, comme une résistance de contact à l'interface tube-solide ou une résistance côté liquide à l'extérieur de la couche en croissance.

La résistance totale par unité de longueur est donc

[ \mathcal{R}_{\text{tot}}'

\frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s}. ]

Calculer le taux d'évacuation thermique instantané

Sous l'approximation pseudo-stationnaire, la chaleur évacuée par unité de longueur de tube est

[ q'(R_s)

\frac{T_m-T_c}{\mathcal{R}_{\text{tot}}'}. ]

Ce flux thermique est évalué en utilisant le rayon d'interface actuel. Par conséquent, le calcul est quasi-stationnaire mais capture tout de même la dépendance temporelle causée par la croissance de la couche.

Dérivation du taux de croissance du solide

Écrire le bilan de chaleur latente

Lorsque l'interface se déplace vers l'extérieur de (dR_s), le volume de solide nouvellement formé par unité de longueur de tube est

[ dV'=2\pi R_s,dR_s. ]

La chaleur latente correspondante libérée est

[ dQ_{\text{latent}}'

\rho_s L,2\pi R_s,dR_s. ]

En divisant par le temps et en égalant ce taux à la chaleur évacuée, on obtient

[ q'(R_s)

2\pi R_s\rho_s L\frac{dR_s}{dt}. ]

La substitution de l'expression de la résistance produit l'équation de croissance régissant le système :

[ \boxed{ \frac{dR_s}{dt}

\frac{T_m-T_c} { 2\pi R_s\rho_s L \left[ \frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s} \right] } } ]

Cette équation donne le taux de croissance radial instantané de la couche solide.

Interpréter la diminution du taux de croissance

À mesure que (R_s) augmente, la résistance du solide

[ \frac{\ln(R_s/R_t)}{2\pi k_s} ]

augmente également. Le flux thermique atteignant la frontière de phase diminue donc, et l'interface avance plus lentement.

Le facteur (2\pi R_s) dans le terme de chaleur latente compte également : un incrément radial donné crée plus de volume solide par unité de longueur de tube à mesure que la circonférence de l'interface grandit.

Simplifications utiles pour les estimations techniques

Utiliser une approximation planaire pour une couche mince

Si l'épaisseur du solide est faible par rapport au rayon du tube, (\delta \ll R_t), la résistance cylindrique peut être approximée par une résistance planaire. L'équation de croissance devient

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L \left[ \frac{1}{h_i} + \frac{t_w}{k_w} + \frac{\delta}{k_s} \right]}, ]

où (t_w) est l'épaisseur de la paroi du tube.

Si les résistances du film de refroidissement et de la paroi sont regroupées en une résistance constante (R_0), alors

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L\left(R_0+\delta/k_s\right)}. ]

Intégrer pour estimer l'épaisseur en fonction du temps

Pour des propriétés et une différence de température constantes, l'intégration donne

[ t

\frac{\rho_s L}{T_m-T_c} \left[ R_0\delta+\frac{\delta^2}{2k_s} \right]. ]

Ce résultat montre deux régimes limites.

À faible épaisseur, les résistances du film de refroidissement et de la paroi dominent, donc la croissance est approximativement linéaire dans le temps. À plus grande épaisseur, la conduction dans le solide domine, et le comportement familier de type diffusion apparaît :

[ \delta \propto \sqrt{t}. ]

Utiliser la forme cylindrique pour des dépôts substantiels

Lorsque l'épaisseur de la couche n'est pas faible par rapport au rayon du tube, l'approximation planaire peut sous-estimer la géométrie et la résistance changeantes. L'équation cylindrique doit alors être intégrée numériquement, ou analytiquement après avoir développé les termes de résistance.

Pour une unité pilote basée sur des tubes, cette distinction est importante car les dépôts peuvent devenir assez épais pour que la surface de transfert thermique locale et la résistance de conduction logarithmique influencent toutes deux le résultat.

Application du modèle à une unité de traitement basée sur des tubes

Estimer le coefficient de transfert thermique effectif

Le modèle de résistance identifie quelle partie du système limite la croissance. Le film de refroidissement interne est représenté par (1/h_i), tandis que le solide en croissance contribue par (\delta/k_s) dans la limite planaire.

Un coefficient global effectif peut être défini à partir de la résistance totale, sous réserve d'une zone de référence cohérente :

[ U=\frac{1}{R_{\text{total}}}. ]

Cela permet de comparer les taux de croissance mesurés avec les valeurs calculées et peut révéler si la limitation dominante est la convection côté liquide de refroidissement, la conduction de la paroi du tube ou la conduction de la couche solide.

Évaluer les conditions de fonctionnement du liquide de refroidissement

La température et le débit du liquide de refroidissement affectent la force motrice thermique et le coefficient interne (h_i). Abaisser (T_c) augmente généralement la force motrice (T_m-T_c), tandis qu'augmenter la vitesse du liquide de refroidissement peut réduire la résistance du film interne.

Le modèle peut donc être utilisé pour prédire comment les changements dans les conditions d'alimentation en liquide de refroidissement affectent le taux de croissance initial et le temps nécessaire pour atteindre une épaisseur de couche cible.

Prendre en compte les effets de température et de composition

Dans les systèmes de cristallisation, (T_m) peut varier avec la concentration, la pression ou le rejet de soluté à l'interface. Si le liquide devient progressivement concentré, la température de changement de phase et la force motrice disponible peuvent changer pendant l'opération.

Dans ce cas, le bilan doit être mis à jour en fonction de la position, du temps ou de la composition locale plutôt que d'utiliser une seule constante (T_m).

Coupler le modèle au comportement axial du tube

Une unité de transfert thermique réelle peut avoir une température de liquide de refroidissement qui change le long du tube. La force motrice locale est alors

[ \Delta T(z,t)=T_m(z,t)-T_c(z,t), ]

et l'équation de croissance doit être résolue le long de la longueur du tube ainsi que dans le temps.

Cela produit une épaisseur locale (R_s(z,t)), qui peut identifier les régions où un dépôt excessif, un blocage ou une cristallisation non uniforme est probable.

Comprendre les compromis

L'approximation omet le chauffage sensible transitoire

Le modèle pseudo-stationnaire néglige le temps nécessaire pour refroidir le solide nouvellement formé et le liquide proche. Cela est plus défendable lorsque le nombre de Stefan est faible, ce qui signifie que la chaleur latente est grande par rapport à la chaleur sensible associée au changement de température.

Lorsque les effets sensibles sont significatifs, un modèle complet de frontière mobile transitoire est plus approprié.

Des propriétés constantes peuvent être irréalistes

La conductivité thermique, la densité, la chaleur latente, la viscosité et les coefficients de transfert thermique peuvent varier avec la température et la composition. Les traiter comme constants simplifie l'analyse mais peut biaiser l'épaisseur et les taux de croissance prédits.

La variation des propriétés doit être incluse lorsque la plage de fonctionnement est large ou lorsque les données expérimentales montrent des écarts systématiques par rapport au modèle simple.

Les conditions de surface peuvent modifier le comportement de l'interface

Le modèle suppose généralement une interface bien définie et un flux thermique radial unidimensionnel. La rugosité, l'encrassement, le comportement de nucléation, la convection naturelle, le cisaillement et la morphologie des cristaux peuvent violer ces hypothèses.

Pour la cristallisation, le dépôt peut également être poreux ou dendritique plutôt qu'une couche solide dense. Sa conductivité effective et sa surface interfaciale peuvent alors différer considérablement des valeurs idéales.

La croissance peut réduire la capacité utile du processus

Une couche solide plus épaisse augmente la résistance thermique et peut réduire la capacité de transfert thermique effective de l'unité. Elle peut également restreindre les passages de flux, augmenter la perte de charge ou provoquer des changements de qualité du produit en raison d'un temps de séjour excessif ou d'une taille de cristal non uniforme.

Par conséquent, maximiser le taux de croissance instantané n'est pas toujours identique à l'optimisation de l'opération globale de l'unité.

Faire le bon choix pour votre objectif

Utilisez le modèle pseudo-stationnaire comme premier outil d'ingénierie, puis augmentez la complexité du modèle uniquement lorsque les hypothèses affectent matériellement le résultat.

Si votre objectif principal est d'estimer l'épaisseur du solide au fil du temps : Utilisez le bilan de résistance cylindrique lorsque le dépôt est substantiel et l'équation intégrée planaire pour les estimations de couches minces.
Si votre objectif principal est d'identifier la limitation du transfert thermique : Séparez les résistances du film de refroidissement, de la paroi du tube et de la couche solide pour déterminer quelle variable opérationnelle aura le plus grand effet.
Si votre objectif principal est d'optimiser le fonctionnement du liquide de refroidissement : Évaluez la température et le débit du liquide de refroidissement ensemble, car la force motrice et la convection interne contrôlent toutes deux le taux d'évacuation de la chaleur latente.
Si votre objectif principal est de prédire une croissance non uniforme dans une unité pilote : Couplez l'équation de croissance locale aux bilans de température et de composition axiaux du liquide de refroidissement.
Si votre objectif principal est une prédiction haute fidélité : Incluez la chaleur sensible transitoire, les propriétés variables, la résistance externe du liquide et la morphologie du dépôt lorsque le comportement expérimental ne peut être expliqué par le modèle pseudo-stationnaire.

Un bilan de résistance à l'état pseudo-stationnaire convertit une frontière de solidification mobile en un calcul pratique de transfert thermique tout en clarifiant pourquoi la croissance ralentit à mesure que la couche solide devient plus épaisse.

Tableau récapitulatif :

Élément clé	Aperçu technique
Hypothèse centrale	Les profils thermiques s'ajustent plus rapidement que l'interface solide-liquide ne se déplace.
Force motrice	Température de changement de phase moins température globale du liquide de refroidissement, ΔT = Tm − Tc.
Chemin de transfert thermique	Film de refroidissement interne, paroi du tube et couche solide en croissance.
Bilan de croissance	La chaleur évacuée est égale à la chaleur latente libérée par le solide nouvellement formé.
Croissance cylindrique	Utilisez le rayon d'interface actuel et la résistance solide logarithmique pour les dépôts épais.
Estimation de couche mince	La croissance est initialement quasi-linéaire, puis approche δ ∝ √t à mesure que la résistance solide augmente.
Contrôles principaux	Température du liquide de refroidissement, vitesse du liquide de refroidissement, propriétés thermiques, chaleur latente et épaisseur du dépôt.
Limites du modèle	Le chauffage sensible transitoire, les propriétés variables, les effets de composition, l'encrassement et les dépôts poreux peuvent nécessiter une modélisation de plus haute fidélité.

Prédire et améliorer votre processus basé sur des tubes

Besoin de contrôler la congélation, la cristallisation ou le dépôt solide dans une unité pilote ? Contactez LABPARK pour discuter de votre application avec notre équipe.

LABPARK fournit des unités pilotes d'opérations unitaires éducatives et professionnelles pour les universités, les instituts de recherche et les entreprises, y compris des solutions pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que le traitement de l'eau et de l'environnement. Nous pouvons vous aider à évaluer les conditions du liquide de refroidissement, la résistance au transfert thermique, la croissance de l'interface, la variation axiale de la température et les exigences de mise à l'échelle afin que vous puissiez :

Estimer l'épaisseur de la couche solide et le temps de croissance de manière plus fiable.
Identifier si le débit du liquide de refroidissement, la conduction de la paroi du tube ou la résistance du dépôt limite les performances.
Évaluer les risques de blocage, de perte de charge et de croissance non uniforme avant l'opération pilote.
Sélectionner l'instrumentation et les conditions de fonctionnement pratiques pour l'enseignement, la recherche et le développement de la production.

Contactez LABPARK pour définir votre matériau cible, la géométrie du tube, les conditions du liquide de refroidissement et l'objectif du processus, et développer une approche d'unité pilote adaptée à vos besoins.