La conduction thermique microscopique autour de particules sphériques est modélisée en résolvant l'équation de Laplace stationnaire autour d'une sphère isolée, puis en moyennant sa réponse thermique sur un volume représentatif de suspension. Pour une particule sphérique de conductivité (k_B) dispersée dans un liquide continu de conductivité (k_A), l'application d'un gradient de température uniforme produit un champ de température perturbé autour de la particule. La perturbation décroît comme (r^{-2}), et le moyennage volumique du flux thermique résultant donne une conductivité effective globale qui peut être utilisée dans les calculs de transfert thermique des installations pilotes.

Conclusion essentielle : Pour une suspension diluée de particules sphériques, le problème local de conduction se réduit à la relation de Maxwell pour la conductivité effective : [ \frac{k_{\mathrm{eff}}}{k_A}

1+ 3\left(\frac{\gamma-1}{\gamma+2}\right)\phi, \qquad \gamma=\frac{k_B}{k_A} ] où (\phi) est la fraction volumique de particules. Le résultat est fiable lorsque les particules sont approximativement sphériques, bien dispersées et suffisamment diluées pour que les interactions entre particules restent faibles.

Définition du problème microscopique

Les deux phases

La phase A est la matrice liquide continue de conductivité thermique (k_A). La phase B est la particule solide sphérique de conductivité thermique (k_B) et de rayon (a).

La suspension est initialement représentée par un liquide contenant une particule dans un milieu par ailleurs illimité. Ce problème de particule isolée constitue la base de la modélisation d'une suspension diluée.

Le gradient macroscopique imposé

Un gradient de température uniforme est appliqué loin de la particule. En coordonnées sphériques, la température non perturbée du liquide peut s'écrire

[ T_\infty = -G r\cos\theta ]

où (G) est l'amplitude du gradient de température imposé, (r) est la distance par rapport au centre de la particule et (\theta) est mesuré par rapport à la direction du gradient.

La particule modifie ce champ parce que sa conductivité diffère de celle du liquide environnant.

Résolution de la conduction autour de la sphère

L'équation directrice

En régime stationnaire et en l'absence de production interne de chaleur, la température satisfait l'équation de Laplace dans les deux régions :

[ \nabla^2 T_A=0 ]

dans le liquide, et

[ \nabla^2 T_B=0 ]

à l'intérieur de la particule.

Cela suppose que la conduction domine le problème local à l'échelle de la particule et que l'interface particule-liquide est à l'équilibre thermique.

Développement en polynômes de Legendre

Comme le champ imposé est axisymétrique autour de la direction du gradient de température, les champs de température peuvent être développés à l'aide de polynômes de Legendre. Pour le problème du gradient uniforme du premier ordre, seul le premier polynôme est nécessaire :

[ P_1(\cos\theta)=\cos\theta ]

La température du côté liquide prend donc la forme

[ T_A=-G r\cos\theta + C\frac{\cos\theta}{r^2} ]

et la température du côté particule prend la forme

[ T_B=-D r\cos\theta ]

Le premier terme à l'extérieur de la sphère est le champ de température non perturbé. Le second terme est la perturbation induite par la particule.

Pourquoi la perturbation décroît comme (r^{-2})

Le terme de perturbation,

[ C\frac{\cos\theta}{r^2}, ]

est la contribution dipolaire générée par le contraste de conductivité entre la particule et le liquide. Son amplitude diminue rapidement avec la distance, de sorte que la particule influence fortement le champ de température à proximité, mais exerce une influence progressivement plus faible à mesure que l'on s'en éloigne.

Cette décroissance permet d'intégrer la solution locale dans une procédure de moyennage macroscopique.

Conditions à l'interface

À la surface de la particule, (r=a), deux conditions sont appliquées.

La température est continue :

[ T_A=T_B ]

Le flux thermique normal est également continu :

[ -k_A\frac{\partial T_A}{\partial r}

-k_B\frac{\partial T_B}{\partial r} ]

La résolution de ces conditions détermine les constantes (C) et (D) en fonction de (G), (a), (k_A) et (k_B).

L'intensité de la perturbation dépend du rapport de conductivité

[ \gamma=\frac{k_B}{k_A}. ]

Conversion de la solution locale en conductivité globale

Moyennage du flux thermique sur un volume représentatif

Une suspension contient à la fois du liquide et des particules, de sorte que le flux thermique microscopique varie spatialement. À proximité d'une particule, les lignes de flux se courbent vers ou à l'écart de la particule selon que (k_B) est supérieur ou inférieur à (k_A).

Pour obtenir une propriété globale, le flux thermique local est moyenné sur un macro-volume représentatif contenant les proportions appropriées de liquide et de particules. Le flux moyenné est ensuite relié au gradient de température moyenné par la loi de Fourier :

[ \langle \mathbf q\rangle

-k_{\mathrm{eff}} \langle\nabla T\rangle ]

Cela définit la conductivité thermique effective de la suspension.

Relation de Maxwell pour les suspensions diluées

Pour des particules sphériques dispersées dans une matrice continue, le résultat dilué du premier ordre est

[ \frac{k_{\mathrm{eff}}}{k_A}

1+ 3\left(\frac{\gamma-1}{\gamma+2}\right)\phi ]

où

(k_{\mathrm{eff}}) est la conductivité effective de la suspension,
(k_A) est la conductivité du liquide,
(k_B) est la conductivité des particules,
(\gamma=k_B/k_A),
(\phi) est la fraction volumique de particules.

L'équation montre que la variation de conductivité dépend à la fois de la charge en particules et du contraste entre la conductivité des particules et celle du liquide.

Interprétation du rapport de conductivité

Si (k_B>k_A), alors (\gamma>1), et les particules augmentent la conductivité effective.

Si (k_B<k_A), alors (\gamma<1), et les particules réduisent la conductivité effective.

Lorsque (\gamma) devient très grand, la contribution des particules tend vers une limite supérieure pour ce modèle dilué. L'effet n'augmente pas indéfiniment simplement parce que la conductivité des particules devient extrêmement élevée.

Vérifications des cas limites

Lorsque (\phi=0), l'expression donne

[ k_{\mathrm{eff}}=k_A, ]

comme requis pour un liquide sans particules.

Lorsque (k_B=k_A), (\gamma=1), et le résultat donne à nouveau

[ k_{\mathrm{eff}}=k_A. ]

Ces limites confirment que des particules ayant la même conductivité que le liquide ne provoquent aucune perturbation thermique.

Application du modèle aux opérations en installation pilote

Échangeurs de chaleur pour boues

Dans un échangeur de chaleur pour boues, (k_{\mathrm{eff}}) fournit une estimation de la propriété conductrice globale de la suspension. Cette valeur peut être utilisée avec une corrélation convective appropriée pour estimer les coefficients de film et le coefficient global de transfert thermique.

Le modèle de conductivité ne détermine pas à lui seul les performances complètes de l'échangeur. Le régime d'écoulement, la viscosité, la sédimentation des particules, l'encrassement des parois et les turbulences induites par les particules doivent également être pris en compte.

Cuves agitées et fermenteurs

Pour les suspensions agitées, la conductivité effective aide à caractériser le transport thermique au sein du volume mélangé. Elle peut contribuer aux calculs de bilan énergétique et au dimensionnement préliminaire des jaquettes, des serpentins ou des échangeurs de chaleur externes en recirculation.

Le modèle de conduction locale est distinct du modèle d'agitation. Le mélange détermine l'efficacité avec laquelle la suspension transporte la chaleur par mouvement global, tandis que (k_{\mathrm{eff}}) représente la réponse conductrice homogénéisée.

Relier la conductivité au transfert thermique global

La conductivité effective peut influencer les corrélations convectives par l'intermédiaire de groupes sans dimension tels que le nombre de Prandtl,

[ Pr=\frac{C_p\mu}{k_{\mathrm{eff}}}, ]

lorsque la suspension est traitée comme un fluide homogène doté de propriétés effectives.

Les coefficients convectifs obtenus peuvent ensuite être intégrés à un calcul global de transfert thermique tel que

[ Q=UA\Delta T_m. ]

Les valeurs des propriétés doivent être évaluées à des températures représentatives, car la conductivité, la viscosité et la capacité thermique du liquide peuvent varier considérablement au cours du fonctionnement de l'installation pilote.

Comprendre les compromis

L'hypothèse de dilution

La relation de Maxwell est fondamentalement un résultat applicable aux suspensions diluées. Elle suppose que la perturbation de température produite par une particule ne se superpose pas fortement à celle des particules voisines.

À des fractions volumiques plus élevées, les interactions entre particules, l'encombrement et les modifications des chemins suivis par le liquide peuvent rendre l'expression du premier ordre inexacte.

Forme des particules et distribution granulométrique

La dérivation suppose des particules sphériques. Les particules irrégulières, allongées, poreuses ou en forme de paillettes créent des champs thermiques différents et nécessitent généralement des modèles dépendant de la forme ou une caractérisation expérimentale.

La taille des particules n'apparaît pas explicitement dans l'équation de Maxwell pour les suspensions diluées. Cependant, elle peut tout de même affecter le comportement réel en installation pilote par l'intermédiaire de la sédimentation, du mouvement brownien, de la résistance interfaciale, de l'agglomération et des changements du régime d'écoulement de la suspension.

Résistance thermique interfaciale

Le modèle de base suppose un contact thermique parfait entre le solide et le liquide. Si une résistance thermique de contact existe à l'interface particule-liquide, le flux thermique est discontinu sous la forme simple utilisée ci-dessus, et la conductivité effective peut être sensiblement différente.

Cela est particulièrement pertinent lorsque les particules sont revêtues, mal mouillées, poreuses ou séparées du liquide par une couche interfaciale.

Dispersion et agglomération des particules

Le modèle considère que les particules sont uniformément dispersées. La sédimentation, la floculation ou l'agglomération produisent des concentrations locales qui ne sont pas représentées par une seule valeur globale de (\phi).

Dans les installations pilotes, des échantillons prélevés à différents endroits peuvent donc présenter des conductivités apparentes différentes, même lorsque la formulation nominale reste inchangée.

Conduction et convection

La solution pour une particule sphérique décrit la conduction microscopique dans un référentiel local stationnaire. Elle ne tient pas compte de la circulation des fluides à grande échelle, de la migration des particules ni des turbulences.

Dans une unité en fonctionnement, les performances de transfert thermique mesurées ou déduites combinent conduction, convection, mélange, effets de paroi et éventuellement encrassement. Considérer (k_{\mathrm{eff}}) comme l'unique propriété déterminante peut conduire à des prévisions incorrectes des équipements.

Faire le bon choix selon votre objectif

Utilisez le modèle microscopique comme estimation de propriété, puis combinez-le avec les modèles hydrodynamiques et à l'échelle de l'équipement requis par l'opération en installation pilote.

Si votre objectif principal est d'estimer la conductivité globale de la suspension : Utilisez la relation de Maxwell avec des valeurs de (k_A), (k_B) et de la fraction volumique de particules (\phi) cohérentes avec la température, à condition que la suspension soit diluée et bien dispersée.
Si votre objectif principal est de dimensionner un échangeur de chaleur pour boues : Utilisez (k_{\mathrm{eff}}) comme l'une des données d'entrée des calculs de transfert thermique convectif et global, tout en tenant compte séparément de la viscosité, du nombre de Reynolds, de la sédimentation, de l'encrassement et de la configuration d'écoulement.
Si votre objectif principal est de modéliser des suspensions concentrées ou non idéales : Remplacez l'approximation de particules sphériques diluées par un modèle dépendant de la concentration, de la forme ou validé expérimentalement.
Si votre objectif principal est d'interpréter des mesures en installation pilote : Comparez la conductivité mesurée à la prédiction de Maxwell afin d'identifier les effets tels que l'agglomération, la résistance interfaciale, les variations de température ou la distribution non uniforme des solides.

Un modèle fiable pour une installation pilote commence par la solution de conduction autour d'une particule sphérique, mais sa fiabilité dépend de la validation de la propriété effective obtenue par rapport à la suspension réelle et aux conditions opératoires.

Tableau récapitulatif :

Facteur clé	Enseignement du modèle	Pertinence pour l'installation pilote
Physique directrice	Conduction stationnaire décrite par l'équation de Laplace autour d'une sphère isolée	Établit les champs microscopiques de température et de flux thermique
Modèle de propriété effective	Relation de Maxwell : `k_eff/k_A = 1 + 3[(gamma - 1)/(gamma + 2)]phi`	Estime la conductivité globale de la suspension
Rapport de conductivité	`gamma = k_B/k_A`	Prévoit si les particules augmentent ou diminuent la conductivité
Charge en particules	Fraction volumique `phi`	Des concentrations plus élevées peuvent nécessiter des modèles non dilués
Hypothèses clés	Particules sphériques et bien dispersées avec un contact thermique parfait	Les résultats peuvent varier en cas d'agglomération, de sédimentation ou de résistance interfaciale
Utilisation pour les équipements	Combiner `k_eff` avec des modèles de viscosité, d'écoulement, d'encrassement et de convection	Contribue à la conception d'échangeurs pour boues, de cuves, de fermenteurs et de systèmes de recirculation

Besoin d'un modèle fiable des propriétés thermiques pour votre installation pilote ? Contactez LABPARK pour discuter de votre application. LABPARK fournit des installations pilotes éducatives et professionnelles pour les opérations unitaires en génie chimique, bioprocédés et biotechnologies, ainsi que pour 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 des configurations d'installations pilotes adaptées, à intégrer des équipements de chauffage et de refroidissement et à évaluer le comportement du transfert thermique des suspensions dans des conditions opératoires pertinentes. Partagez les propriétés de vos particules, votre système liquide, la concentration en solides, la plage de températures et l'objectif du procédé afin que nous puissions vous aider à développer une solution mieux validée pour l'enseignement, la recherche, le changement d'échelle ou le développement de procédés industriels.