Les développements multipolaires fournissent une approximation en champ lointain de la concentration de soluté autour d’une particule en suspension. Les ingénieurs représentent le flux surfacique détaillé de la particule au moyen des moments d’une fonction de Green : le terme monopolaire, proportionnel à la libération ou à la consommation nette de soluté par la particule et décroissant comme (r^{-1}), ainsi que des termes d’ordre supérieur tels que le dipôle, qui décroît comme (r^{-2}) et traduit les variations directionnelles du flux. Lorsque la particule ne présente aucune génération nette de matière, le monopôle s’annule et le dipôle devient généralement la principale contribution à longue portée.

À retenir : Les développements multipolaires remplacent un calcul complexe du transfert de matière à la surface par un petit ensemble de moments sources ayant une signification physique. Le monopôle décrit l’échange total de soluté de la particule, tandis que le dipôle décrit comment un gradient de concentration imposé ou un flux surfacique non uniforme déforme le champ environnant.

Pourquoi les concentrations en champ lointain sont importantes dans les équipements pilotes

Le problème microscopique

Dans une cuve de lixiviation, un réservoir de suspension ou un bioréacteur, la concentration de soluté est régie par la diffusion et éventuellement par la convection autour des particules ou des cellules. La condition aux limites sur chaque surface particulaire irrégulière peut intégrer une cinétique de réaction locale, un partage, une perméabilité ou un flux de matière variable spatialement.

Résoudre cette condition aux limites en chaque point de surface pour chaque particule est généralement irréalisable dans une suspension à l’échelle pilote. L’analyse multipolaire condense la perturbation de concentration qui en résulte en termes pouvant être évalués à des distances grandes devant la taille de la particule.

Les échelles de longueur pertinentes

Soit (a) une taille caractéristique de la particule et (r) la distance entre le centre de la particule et le point d’observation. Le développement est particulièrement utile lorsque

[ r \gg a. ]

Cette condition sépare le champ proche, où la géométrie de la particule et l’écoulement local dominent, du champ lointain, où la particule peut être représentée approximativement par une source ponctuelle, un dipôle ou une combinaison de moments d’ordre supérieur.

Comment le développement représente le champ de concentration

À partir de la fonction de Green

Pour une phase continue dominée par la diffusion, la perturbation de concentration peut s’écrire schématiquement sous la forme d’une intégrale de surface faisant intervenir la fonction de Green (G) :

[ c'(\mathbf{x}) \sim \int_{S} G(\mathbf{x},\mathbf{y}),q(\mathbf{y}),dS_{\mathbf{y}}, ]

où (q(\mathbf{y})) est le flux local de soluté sur la surface de la particule (S), et (c') la perturbation de concentration par rapport au champ non perturbé.

Pour un fluide tridimensionnel non borné, la fonction de Green principale s’écrit

[ G(\mathbf{x},\mathbf{y}) \propto \frac{1}{|\mathbf{x}-\mathbf{y}|}. ]

Le développement de ce noyau autour d’un point de référence, généralement le centre de la particule, produit la série multipolaire.

Le terme monopolaire

Le coefficient monopolaire est déterminé par le flux surfacique intégré :

[ Q=\int_S q(\mathbf{y}),dS. ]

La concentration en champ lointain correspondante se comporte comme

[ c'_{\text{mono}}(\mathbf{x}) \propto \frac{Q}{D r}, ]

où (D) est la diffusivité de la phase continue. La dépendance en (r^{-1}) signifie que le monopôle domine suffisamment loin lorsque la particule présente un taux net non nul de libération ou d’absorption.

Physiquement, (Q) représente l’échange total de soluté entre la particule et le liquide environnant. Dans la lixiviation solide-liquide, il correspond au soluté net libéré par la particule ; dans le cas d’une particule absorbante ou réactive, il représente la consommation nette.

Le terme dipolaire

Le coefficient dipolaire dépend du premier moment spatial du flux surfacique, ainsi que du gradient de concentration imposé et de la condition aux limites de la particule. Sa contribution en champ lointain s’écrit

[ c'_{\text{dip}}(\mathbf{x}) \propto \frac{\mathbf{p}\cdot \mathbf{r}}{r^3}, ]

avec une amplitude globale proportionnelle à (r^{-2}).

Le dipôle décrit le comportement directionnel. Un côté de la particule peut libérer ou absorber davantage de soluté que l’autre en raison d’un gradient de concentration externe, de l’hydrodynamique locale ou de propriétés de surface asymétriques.

Le cas du monopôle nul

Si la particule ne présente aucune génération nette de matière,

[ Q=0, ]

la contribution monopolaire disparaît. Le dipôle devient alors le terme dominant en champ lointain, même si la particule peut continuer à redistribuer ou à perturber fortement le soluté localement.

Ce cas est particulièrement important pour les particules imperméables et les inclusions non réactives. Elles ne créent ni ne détruisent de soluté, mais modifient le transport en redirigeant les chemins de diffusion et en perturbant le gradient de concentration environnant.

Application de la méthode à des particules complexes

Représentation d’une géométrie de surface quelconque

La particule n’a pas besoin d’être sphérique pour que la méthode fonctionne. Sa géométrie irrégulière intervient par l’intermédiaire de la distribution du flux surfacique et des coefficients multipolaires qui en résultent.

Un calcul pratique se déroule généralement comme suit :

Résoudre ou approximer le problème local de transport à la surface.
Intégrer le flux surfacique afin d’obtenir le monopôle (Q).
Calculer les premiers moments spatiaux et les moments d’ordre supérieur afin d’obtenir les coefficients dipolaire et d’ordre supérieur.
Évaluer le développement tronqué aux points d’observation situés dans le champ lointain.

La géométrie détaillée est ainsi conservée dans les coefficients, plutôt que d’être introduite explicitement dans chaque calcul en champ lointain.

Séparation des comportements local et global

Le développement ne remplace pas une solution en champ proche. À proximité de la particule, les caractéristiques individuelles de la surface, les couches limites, les sillages et les chemins de diffusion internes peuvent être importants.

Les termes multipolaires offrent plutôt un raccord contrôlé entre le problème local à l’échelle de la particule et le modèle de suspension à plus grande échelle. Ils sont ainsi utiles pour coupler le transfert de matière à l’échelle particulaire aux bilans à l’échelle pilote.

Prise en compte des gradients de concentration imposés

Dans une suspension, les particules sont généralement soumises à une concentration de fond non uniforme plutôt qu’à une concentration parfaitement uniforme. La réponse de la particule à ce gradient de fond contribue au moment dipolaire.

Le développement permet ainsi de décrire non seulement la libération nette de soluté, mais aussi la manière dont les particules modifient un champ de concentration existant. Cette distinction est essentielle pour estimer les propriétés de transport effectives.

Relier les multipôles à la diffusivité à l’échelle de la suspension

D’une particule à plusieurs particules

Pour une suspension diluée, la perturbation en champ lointain de chaque particule peut être combinée afin d’estimer le champ de concentration moyen de l’ensemble. Le monopôle contrôle l’intensité moyenne de la source ou du puits, tandis que les termes dipolaires décrivent l’influence des particules sur les gradients.

À des concentrations plus élevées, les particules interagissent par l’intermédiaire de champs de perturbation qui se chevauchent. Un développement fondé sur une particule isolée nécessite alors des corrections prenant en compte la disposition des particules, les interactions hydrodynamiques et la polarisation collective de la concentration.

Modèles de transport effectif de type Maxwell

Pour des suspensions diluées contenant des particules capables de partager ou de diffuser le soluté, une relation de type Maxwell est

[ \frac{K_{\mathrm{eff}}D_{\mathrm{eff}}}{D}

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

où :

(D) est la diffusivité du fluide continu,
(D_{\mathrm{eff}}) est la diffusivité effective de la suspension,
(K_{\mathrm{eff}}) est le coefficient de partage effectif,
(\phi) est la fraction volumique de particules,
(K) est le coefficient de partage du soluté entre la particule et le fluide,
(D_s) est la diffusivité de la particule, et
(\gamma=KD_s/D).

La réponse dipolaire d’une inclusion individuelle constitue la principale grandeur microscopique à la base de ce type de description par milieu effectif. Elle mesure la réponse de l’inclusion à un gradient de concentration imposé à grande échelle.

Particules sphériques imperméables

Pour des sphères totalement imperméables, (\gamma=0). Le résultat prédit dans la limite diluée dépend de la convention utilisée pour la concentration :

[ \frac{D_{\mathrm{eff}}}{D}

1-\frac{3}{2}\phi ]

lorsque la référence est la concentration dans la phase fluide, ou

[ \frac{D_{\mathrm{eff}}}{D}

1-\frac{1}{2}\phi ]

lorsque la référence est la concentration moyennée sur le volume.

Ces résultats ne sont pas contradictoires. Ils utilisent des définitions différentes de la concentration moyenne et décrivent donc des coefficients effectifs différents.

Comprendre les compromis

Précision en champ lointain et détail géométrique

Un développement multipolaire tronqué est efficace, car il ne conserve que les moments dominants. Sa précision diminue lorsque le point d’observation se rapproche de la particule ou lorsque celle-ci devient très allongée, rugueuse ou autrement mal représentée par un développement d’ordre faible.

Les multipôles d’ordre supérieur améliorent la précision locale, mais augmentent la complexité du calcul et de l’interprétation.

Hypothèses de prédominance de la diffusion

Les formes simples (r^{-1}) et (r^{-2}) s’appliquent directement à une diffusion stationnaire tridimensionnelle dans un milieu effectivement non borné. Une advection importante, les parois de la cuve, un fonctionnement transitoire ou une agitation turbulente peuvent modifier la fonction de Green et la décroissance spatiale.

Les équipements pilotes nécessitent donc d’examiner le nombre de Péclet pertinent, la géométrie de la cuve et le régime de mélange avant d’appliquer la forme correspondant à un milieu non borné.

Suspensions diluées et concentrées

Les expressions de diffusivité effective fondées sur les multipôles de particules isolées sont particulièrement fiables pour de faibles fractions volumiques de particules. Dans les boues concentrées, le chevauchement des champs de diffusion et les contacts entre particules peuvent rendre importants les effets à deux particules et les effets collectifs.

Utiliser une formule pour suspension diluée sans vérifier (\phi), l’espacement entre les particules et les conditions opératoires peut conduire à attribuer à tort à la physique à l’échelle particulaire un coefficient apparent de transfert de matière.

Le flux net n’est pas identique au flux local

Un monopôle nul ne signifie pas que la particule n’a aucun effet sur le transfert de matière. Cela signifie seulement que le flux surfacique intégré est nul.

Un dipôle non nul peut néanmoins modifier considérablement les gradients de concentration et la diffusivité effective, en particulier autour de particules imperméables ou de particules exposées à un gradient externe.

Faire le bon choix selon votre objectif

Utilisez le développement en fonction de l’échelle et de la grandeur que vous devez prédire :

Si votre objectif principal est la libération ou l’absorption nette de soluté : Calculez le monopôle à partir du flux intégré à la surface de la particule et utilisez son champ en (r^{-1}) pour estimer la perturbation de concentration à longue portée.
Si votre objectif principal est la déformation du gradient autour de particules non réactives : Annulez le monopôle associé au flux net et évaluez la réponse dipolaire au gradient de concentration imposé.
Si votre objectif principal est le comportement de la suspension à l’échelle pilote : Combinez les coefficients multipolaires d’une particule avec un modèle de moyenne volumique ou de milieu effectif, tout en vérifiant l’hypothèse de suspension diluée.
Si votre objectif principal est une géométrie particulaire irrégulière : Utilisez un calcul détaillé en champ proche pour déterminer les moments, puis appliquez le développement multipolaire uniquement au-delà de la région à l’échelle de la particule.
Si votre objectif principal est le changement d’échelle expérimental : Comparez les profils de concentration en champ lointain et les diffusivités effectives prédites aux mesures réalisées avec la géométrie réelle de la cuve et le régime d’agitation utilisé.

Les développements multipolaires transforment la physique complexe du transfert de matière à l’échelle particulaire en moments sources et gradients interprétables, qui peuvent être reliés aux performances de la suspension à l’échelle pilote.

Tableau récapitulatif :

| Aspect | Idée clé |\n|---|---|\n| Monopôle | Représente la libération ou l’absorption nette de soluté ; décroît comme (r^{-1}). |\n| Dipôle | Traduit le flux directionnel et la déformation du gradient ; décroît comme (r^{-2}). |\n| Flux net nul | Le monopôle s’annule, faisant du dipôle le terme dominant en champ lointain. |\n| Applicabilité | Particulièrement utile lorsque la distance d’observation est très supérieure à la taille de la particule, (r \gg a). |\n| Particules complexes | La géométrie irrégulière est représentée par les moments du flux surfacique. |\n| Modélisation des suspensions | Les coefficients multipolaires contribuent aux modèles moyennés sur le volume et aux modèles de diffusivité effective. |\n| Principales limites | La précision dépend des conditions dominées par la diffusion, de la concentration en particules, de la géométrie de la cuve et de l’agitation. |

Besoin de relier le transfert de matière à l’échelle particulaire aux performances à l’échelle pilote ? Contactez LABPARK pour discuter des unités pédagogiques et professionnelles de pilotes d’opérations unitaires destinés au génie chimique, aux bioprocédés et aux biotechnologies, ainsi qu’au traitement de l’environnement et de l’eau. LABPARK aide les universités, les instituts de recherche et les entreprises à évaluer la libération, l’absorption et la diffusion des solutés, le comportement des suspensions et le changement d’échelle à l’aide d’équipements pilotes pratiques adaptés aux conditions réelles de fonctionnement.