Les fonctions de base mathématiques sont les harmoniques sphériques selon l’angle et les fonctions de Bessel sphériques selon le rayon, choisies pour satisfaire les conditions aux limites de l’équipement. Pour les modules sphériques axisymétriques, la dépendance angulaire est développée en polynômes de Legendre (P_\ell(\cos\theta)), tandis que la dépendance radiale utilise les fonctions de Bessel sphériques ou leur forme équivalente sinus-sur-rayon. Les valeurs propres admissibles sont déterminées selon que la sphère possède une valeur de surface imposée, un flux nul ou un transfert thermique ou massique convectif.

Le principe directeur consiste à choisir un ensemble de fonctions propres orthogonales qui satisfasse à la fois l’opérateur différentiel sphérique et les conditions aux limites physiques. Les conditions de Dirichlet, de Neumann et de Robin produisent des équations caractéristiques, des valeurs propres et des coefficients de développement différents.

Comment le problème de transport sphérique est formulé

La diffusion thermique et la diffusion de masse ont la même structure mathématique

Pour des propriétés constantes et en l’absence de génération interne, la température ou la concentration peut être représentée par un champ générique (\Theta(r,\theta,\phi,t)). L’équation transitoire prend la forme diffusive

[ \frac{\partial \Theta}{\partial t}

D\nabla^2\Theta, ]

où (D) est la diffusivité pertinente.

Pour la conduction thermique,

[ D=\alpha=\frac{k_T}{\rho \hat C_p}, ]

tandis que pour la diffusion de masse, (D) est la diffusivité de l’espèce.

Les coordonnées sphériques se séparent naturellement

Le laplacien en coordonnées sphériques se sépare en parties radiale et angulaire. Une solution séparée peut s’écrire schématiquement sous la forme

[ \Theta(r,\theta,\phi,t)

R(r)Y(\theta,\phi)G(t). ]

Le facteur dépendant du temps prend généralement la forme

[ G(t)=\exp(-D\lambda^2t), ]

de sorte que les modes spatiaux d’ordre supérieur décroissent plus rapidement.

L’axisymétrie réduit la base angulaire

Si l’équipement et le chargement sont indépendants de l’angle azimutal (\phi), le problème angulaire ne dépend que de

[ \eta=\cos\theta. ]

L’équation angulaire devient alors l’équation de Legendre, dont les solutions régulières sont les polynômes de Legendre

[ P_\ell(\eta)=P_\ell(\cos\theta). ]

Ces polynômes sont orthogonaux sur (-1\leq\eta\leq1), ce qui permet de développer le champ sous la forme

[ \Theta(r,\theta,t)

\sum_{\ell=0}^{\infty} \sum_{n} A_{\ell n} R_{\ell n}(r) P_\ell(\cos\theta) e^{-D\lambda_{\ell n}^2t}. ]

Pour un chargement entièrement tridimensionnel, les polynômes de Legendre sont remplacés par les harmoniques sphériques (Y_\ell^m(\theta,\phi)).

Quelles fonctions de base radiales sont utilisées ?

Les fonctions de Bessel sphériques constituent la base radiale générale

L’équation aux valeurs propres radiales contient l’indice de séparation angulaire (\ell). Sa solution régulière au centre de la sphère est la fonction de Bessel sphérique

[ R_\ell(r)=j_\ell(\lambda r). ]

La seconde solution indépendante, la fonction de Neumann sphérique (y_\ell(\lambda r)), est singulière en (r=0) et est donc exclue pour une sphère solide contenant l’origine.

Pour le mode à symétrie sphérique (\ell=0),

[ j_0(\lambda r)=\frac{\sin(\lambda r)}{\lambda r}. ]

Ainsi, la forme couramment utilisée

[ R_0(r)\propto\frac{\sin(\lambda r)}{r} ]

est spécifiquement la solution radiale (\ell=0), ou le résultat de la transformation

[ \Theta(r,t)=\frac{\Psi(r,t)}{r}. ]

La transformation (\Psi/r) simplifie l’équation radiale

Pour un problème à symétrie sphérique, l’équation radiale de diffusion est

[ \frac{\partial \Theta}{\partial t}

D\left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial\Theta}{\partial r} \right) \right]. ]

En utilisant

[ \Theta(r,t)=\frac{\Psi(r,t)}{r} ]

on transforme l’opérateur radial sphérique en équation de diffusion unidimensionnelle

[ \frac{\partial\Psi}{\partial t}

D\frac{\partial^2\Psi}{\partial r^2}. ]

Les fonctions propres spatiales résultantes sont des combinaisons de (\sin(\lambda r)) et de (\cos(\lambda r)). La régularité au centre élimine la contribution en cosinus dans la formulation habituelle d’une sphère solide.

L’orthogonalité permet les solutions en séries

Les fonctions propres radiales et angulaires forment des ensembles orthogonaux avec le poids sphérique approprié. Pour les modes radiaux, l’élément de volume introduit un poids en (r^2) :

[ \int_0^l r^2 R_{\ell n}(r)R_{\ell m}(r),dr=0, \qquad n\ne m. ]

Cette orthogonalité permet de projeter la distribution initiale de température ou de concentration sur les fonctions propres, ce qui produit les coefficients (A_{\ell n}).

Comment les conditions aux limites déterminent les valeurs propres

Le centre exige une condition de régularité ou de symétrie

Au centre d’une sphère solide, le champ doit rester fini. Le flux radial physique doit également s’annuler :

[ \left.\frac{\partial\Theta}{\partial r}\right|_{r=0}=0. ]

Il s’agit de la condition de symétrie sphérique. Ce n’est pas une simple commodité numérique facultative : elle impose la continuité au centre et élimine les solutions radiales singulières.

Les conditions de Dirichlet imposent une valeur de surface fixée

Une température ou une concentration de surface imposée est représentée par

[ \Theta(l,t)=\Theta_s. ]

Après soustraction de la valeur limite stationnaire, la condition homogène sur la fonction propre devient

[ R_{\ell n}(l)=0. ]

Pour le mode à symétrie sphérique, on obtient

[ \sin(\lambda_n l)=0, ]

et donc

[ \lambda_n l=n\pi, ]

avec un indice choisi selon que le mode zéro est inclus ou non.

Cette condition modélise une sphère dont la surface est maintenue par un bain, une jaquette ou un réservoir externe de concentration soigneusement régulé.

Les conditions de Neumann imposent un flux spécifié

Une surface isolée ou imperméable possède un flux normal nul :

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}=0 ]

pour le transfert thermique, ou

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}=0 ]

pour la diffusion de masse.

Pour la solution transformée à symétrie sphérique, cela produit l’équation caractéristique

[ \lambda_n l=\tan(\lambda_n l). ]

Plus généralement, la dérivée de la fonction de Bessel sphérique appropriée doit s’annuler :

[ \left.\frac{d}{dr}j_\ell(\lambda_{\ell n}r)\right|_{r=l}=0. ]

Une limite isolée ou imperméable conserve l’énergie thermique totale ou la masse totale de l’espèce dans la sphère, sous réserve de toute source ou tout puits interne.

Les conditions de Robin représentent un transfert convectif

Une surface convective est décrite par une condition mixte. Pour le transfert thermique,

[ -k_T\left.\frac{\partial\Theta}{\partial r}\right|_{r=l}

h_T\left[\Theta(l,t)-\Theta_\infty(t)\right]. ]

Pour le transfert de masse, la condition analogue est

[ -D\left.\frac{\partial C}{\partial r}\right|_{r=l}

k_m\left[C(l,t)-C_\infty(t)\right]. ]

Après adimensionnement, le coefficient de transfert apparaît par l’intermédiaire d’un paramètre de type Biot. Pour la solution transformée à symétrie sphérique, l’équation caractéristique peut s’écrire sous une forme telle que

[ \lambda_n l

\left(1-A l^2\right)\tan(\lambda_n l), ]

selon la définition de (A), la convention de signe et la variable dépendante transformée.

Le point physique essentiel est que les conditions de Robin tiennent compte de la résistance externe au solide, notamment de la couche limite fluide, du film côté jaquette ou du film externe de transfert de masse.

Les conditions non homogènes nécessitent un ajustement à l’état stationnaire

Les développements en fonctions propres sont les plus simples lorsque les conditions aux limites sont homogènes. Si la température ou la concentration de surface varie avec le temps, le champ est généralement décomposé en une partie particulière ou quasi stationnaire et une partie transitoire satisfaisant des conditions homogènes.

Cette étape empêche que des termes de bord restent non résolus lors de la transformation de l’équation de diffusion.

Pourquoi les conditions aux limites doivent correspondre à la base

Les conditions aux limites définissent le problème aux valeurs propres

L’opérateur différentiel seul ne détermine pas la base. L’opérateur, associé à la condition au centre, à la condition de surface et aux éventuelles contraintes de symétrie, définit le problème aux valeurs propres de Sturm-Liouville.

Par conséquent, deux sphères ayant des propriétés matérielles identiques peuvent nécessiter des valeurs propres différentes si l’une possède une température de surface imposée et l’autre un chauffage convectif.

La correspondance des fonctions de base élimine les termes de bord indésirables

Lorsqu’on déduit des transformées de Fourier finies ou des transformées en fonctions propres par intégration par parties, des termes de bord apparaissent. Une base satisfaisant les mêmes conditions aux limites homogènes permet à ces termes de s’annuler ou de se convertir en constantes connues.

Utiliser, par exemple, une base sinus pour une limite convective ne permet généralement pas de représenter correctement la physique réelle de la surface. Le problème transformé qui en résulte peut présenter des valeurs propres, des coefficients, ou les deux, incorrects.

Les conditions aux limites doivent refléter le module réel

Pour les équipements d’une installation pilote, le choix mathématique doit suivre l’interface physique :

Un bain régulé ou une composition de surface imposée suggère un comportement de Dirichlet.
Une paroi isolée, un plan de symétrie ou une surface imperméable suggère un comportement de Neumann.
Une jaquette, un gaz en circulation, un film liquide ou une phase externe en vrac suggère généralement un comportement de Robin.
Une entrée, une sortie ou une interface réactive peut nécessiter une condition distincte de valeur ou de flux imposé, plutôt qu’une condition de surface sphérique.

Comprendre les compromis

Les valeurs propres sous forme fermée sont limitées

Les problèmes de Dirichlet et certains problèmes de Neumann produisent des équations caractéristiques simples et des suites explicites de valeurs propres. Les problèmes de Robin conduisent généralement à des équations transcendantes nécessitant la recherche numérique des racines.

Le travail numérique supplémentaire est justifié lorsque la résistance externe au transfert affecte sensiblement le procédé.

La forme sinus-sur-rayon a une portée limitée

Écrire la solution radiale sous la forme

[ \frac{\sin(\lambda r)}{r} ]

est approprié pour le mode à symétrie sphérique (\ell=0). Il ne s’agit pas de la base complète pour une variation angulaire arbitraire, pour laquelle il faut utiliser (j_\ell(\lambda r)).

Utiliser la forme simplifiée pour un chargement non axisymétrique ou fortement non uniforme peut omettre des modes angulaires physiquement importants.

Les hypothèses idéales peuvent ne pas être valables dans les équipements pilotes

La base analytique standard suppose notamment des propriétés constantes, une géométrie sphérique et une équation linéaire de diffusion ou de conduction. Les modules réels peuvent inclure une diffusivité variable, une réaction interne, un mouvement fluide, une résistance de contact, un transfert de surface non uniforme ou une forme sphérique imparfaite.

Ces effets peuvent nécessiter des méthodes numériques telles que les calculs en volumes finis, en éléments finis ou spectraux.

Les conditions aux limites expérimentales sont rarement parfaites

La température d’un bain supposé constant peut varier, et un coefficient de film estimé peut changer avec le débit, la viscosité ou la concentration. La validation du modèle doit donc comparer les données mesurées à la surface et à l’intérieur, plutôt que de considérer la condition aux limites supposée comme exacte.

La condition aux limites constitue souvent l’une des plus grandes sources d’incertitude dans l’estimation de la diffusivité ou des coefficients de transfert.

Faire le bon choix pour votre objectif

La base doit être choisie à partir de la géométrie, de la symétrie et de la physique des interfaces du module pilote.

Si votre objectif principal est le comportement transitoire à symétrie sphérique : utilisez la base de Bessel sphérique (\ell=0), ou la représentation équivalente (\sin(\lambda r)/r), avec la condition de symétrie au centre et l’équation caractéristique de surface appropriée.
Si votre objectif principal est un chauffage ou une concentration non uniforme selon l’angle : développez la solution en polynômes de Legendre pour les systèmes axisymétriques ou en harmoniques sphériques pour les systèmes entièrement tridimensionnels, en utilisant (j_\ell(\lambda r)) selon le rayon.
Si votre objectif principal est une température ou une concentration de surface contrôlée : appliquez une condition de Dirichlet et utilisez des fonctions propres qui s’annulent à la surface après soustraction de la valeur imposée.
Si votre objectif principal est l’isolation, l’imperméabilité ou la symétrie : appliquez une condition de Neumann avec un flux normal nul et utilisez des fonctions propres compatibles avec la dérivée.
Si votre objectif principal est un transfert limité par le côté jaquette ou par un film : appliquez une condition de Robin en utilisant le coefficient de transfert thermique ou massique mesuré ou estimé, puis déterminez les valeurs propres numériquement.
Si votre objectif principal est l’accord avec les mesures de l’installation pilote : vérifiez que les conditions aux limites supposées reflètent les conditions réelles d’entrée, de sortie, de paroi, de symétrie et de fluide externe avant d’ajuster les paramètres de transport.

Associer correctement les fonctions de base de Bessel sphériques ou les harmoniques aux conditions aux limites physiques transforme le modèle transitoire en une représentation défendable de l’installation pilote, plutôt qu’en une simple solution mathématique formelle.

Tableau récapitulatif :

Aspect	Choix clé	Application
Base angulaire	Polynômes de Legendre ou harmoniques sphériques	Chargement axisymétrique ou 3D
Base radiale	Fonctions de Bessel sphériques (j_\ell)	Solutions régulières dans les sphères solides
Symétrie sphérique	Forme (\sin(\lambda r)/r)	Comportement transitoire (\ell=0)
Condition de Dirichlet	Température ou concentration de surface fixée	Bain ou réservoir contrôlé
Condition de Neumann	Flux nul ou spécifié	Surfaces isolées ou imperméables
Condition de Robin	Transfert thermique ou massique convectif	Jackettes, films et phases en circulation
Valeurs propres	Racines des équations caractéristiques	Déterminées par les conditions aux limites
Validation du modèle	Comparer les données prédites et mesurées	Estimation des paramètres de l’installation pilote

Besoin d’un modèle fiable pour le transfert thermique ou massique transitoire dans votre équipement pilote ? Contactez LABPARK pour discuter de votre application. LABPARK fournit des installations pilotes d’opérations unitaires éducatives et professionnelles pour le génie chimique, les bioprocédés et les biotechnologies, ainsi que le traitement environnemental et de l’eau. Nos systèmes aident les universités, les instituts de recherche et les entreprises à étudier le transport sphérique, à évaluer les performances des transferts thermique et massique et à relier les modèles mathématiques aux expériences pratiques. Nous pouvons vous aider à sélectionner les modules appropriés, à définir des conditions aux limites réalistes et à concevoir des dispositifs de formation ou de recherche adaptés aux objectifs de votre procédé.