La diffusion transitoire solide-liquide est régie par les lois de Fick, la conservation de la matière et l’évolution de la différence de concentration entre les phases solide et liquide. À l’intérieur du solide, le soluté se déplace par diffusion vers la surface ; l’agitation transporte ce soluté dans le liquide et maintient approximativement uniforme la concentration du liquide en vrac. Comme le volume de liquide est fini, sa concentration augmente au fur et à mesure de l’extraction, réduisant la force motrice jusqu’à ce que l’équilibre soit atteint ou qu’un taux de récupération spécifié soit obtenu.

Point essentiel : Un modèle d’extraction à l’échelle pilote doit coupler la diffusion transitoire à l’intérieur du solide avec un bilan de matière dynamique de la phase liquide. La diffusion interne, la taille des particules, la diffusivité effective, l’équilibre interfacial, le volume de liquide et l’agitation déterminent ensemble la vitesse d’extraction et la durée de traitement requise.

Le modèle physique régissant le procédé

La première loi de Fick détermine le flux local

La première loi de Fick relie le flux de soluté au gradient de concentration :

[ \mathbf{N}A=-D{\mathrm{eff}}\nabla C_A ]

Ici, (D_{\mathrm{eff}}) est la diffusivité effective du soluté (A) à travers la matrice solide, et (C_A) est sa concentration dans le solide.

Le signe négatif indique que le soluté se déplace des régions à forte concentration vers les régions à faible concentration. Dans une particule poreuse, (D_{\mathrm{eff}}) tient compte des effets de la structure des pores, de la porosité, de la tortuosité, de la saturation, de la température et de la viscosité du liquide.

La deuxième loi de Fick décrit l’épuisement transitoire

Pour une diffusivité constante et en l’absence de réaction dans le solide, le champ de concentration transitoire satisfait :

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}}\nabla^2 C_A ]

Pour une particule sphérique, géométrie courante dans les modèles de lixiviation, cela devient :

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}} \left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial C_A}{\partial r} \right) \right] ]

La géométrie est importante, car les plaques planes, les cylindres et les sphères présentent des longueurs de diffusion et des profils de concentration différents.

La condition initiale représente le solide chargé

Le modèle commence par une distribution initiale du soluté, souvent supposée uniforme :

[ C_A(r,0)=C_{A,0} ]

Les solides plus complexes peuvent présenter un profil de concentration radial dû à un séchage, une imprégnation, une réaction ou une précipitation antérieure. Cette condition initiale peut influencer considérablement la vitesse d’extraction initiale.

La symétrie impose une contrainte au centre de la particule

Pour une particule sphérique, aucun flux net ne traverse le centre :

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

Cette condition empêche l’apparition d’une singularité artificielle au centre et reflète la symétrie physique de la particule.

Comment le volume fini de liquide intervient dans le modèle

Le liquide en vrac est souvent considéré comme parfaitement agité

Avec une agitation énergique, la concentration du liquide en vrac peut être supposée uniforme dans l’espace :

[ C_L=C_L(t) ]

Cette hypothèse évite de devoir résoudre un champ complet de concentration dans la phase liquide. Elle ne supprime pas le transfert de matière ; elle signifie que les gradients de concentration du côté liquide sont représentés par une résistance interfaciale ou supposés suffisamment faibles.

Le bilan du liquide couple les deux phases

Le soluté libéré par le solide s’accumule dans le volume fini de liquide :

[ V_L\frac{dC_L}{dt}

\dot{m}_{A,\mathrm{solid\ to\ liquid}} ]

Pour (N_p) particules identiques de surface (A_p), cela peut s’écrire :

[ V_L\frac{dC_L}{dt}

N_p A_p N_A ]

où (N_A) est le flux de soluté sortant à la surface de la particule.

Le bilan global de soluté fournit une importante vérification de cohérence :

[ M_{A,0}

M_A^{\mathrm{solid}}(t) + V_L C_L(t) ]

en supposant qu’il n’y ait ni perte de soluté, ni dégradation, ni prélèvement d’échantillons.

Le volume fini crée une force motrice décroissante

Une condition interfaciale couramment utilisée est :

[ -D_{\mathrm{eff}} \left. \frac{\partial C_A}{\partial r} \right|_{r=R}

k_L\left(C_{A,i}-C_L\right) ]

où (C_{A,i}) est la concentration interfaciale équivalente dans le liquide et (k_L) le coefficient de transfert de matière du côté liquide.

Si l’équilibre local est représenté par une relation de partage,

[ C_{A,i}=K_{\mathrm{eq}}C_{A,s} ]

la concentration à la surface du solide et la concentration du liquide sont alors liées par (K_{\mathrm{eq}}). Lorsque (C_L) augmente, la différence entre la concentration d’équilibre du liquide et (C_L) diminue, ce qui ralentit l’extraction.

Les principaux principes adimensionnels

Le nombre de Fourier mesure l’avancement de la diffusion

Le temps de diffusion adimensionnel est le nombre de Fourier :

[ \mathrm{Fo}=\frac{D_{\mathrm{eff}}t}{R^2} ]

Ici, (R) est le rayon caractéristique de la particule. De petites valeurs indiquent que la diffusion n’a pénétré qu’une partie du solide, tandis que des valeurs plus élevées indiquent un épuisement plus étendu.

Le temps de diffusion caractéristique varie comme suit :

[ t_D\sim \frac{R^2}{D_{\mathrm{eff}}} ]

Cette dépendance quadratique rend la taille des particules particulièrement importante dans la conception d’une installation pilote. Diviser par deux le rayon des particules peut réduire d’environ un facteur quatre le temps caractéristique de diffusion interne, si les autres limitations restent inchangées.

Le nombre de Biot compare les résistances externe et interne

Un nombre de Biot de transfert de matière compare la résistance du côté liquide à la diffusion interne dans le solide :

[ \mathrm{Bi}m=\frac{k_L R}{D{\mathrm{eff}}} ]

Une grande valeur de (\mathrm{Bi}_m) indique généralement que la résistance externe est relativement faible par rapport à la résistance de diffusion interne. Une faible valeur de (\mathrm{Bi}_m) indique que la couche limite liquide peut contrôler ou influencer significativement la vitesse d’extraction.

La définition exacte peut varier selon la géométrie et les conventions de partage, mais la question sous-jacente reste la même : où se produit la majeure partie de la chute de concentration ?

Les solutions de similitude expliquent le comportement aux temps courts

Pour une diffusion semi-infinie dans un solide ou hors de celui-ci, la variable de similitude est :

[ \eta=\frac{x}{2\sqrt{D_{\mathrm{eff}}t}} ]

L’équation de diffusion peut alors produire des profils de concentration en fonction d’erreur, tels que :

[ \Theta(\eta)=\operatorname{erfc}(\eta) ]

Ces solutions montrent que la profondeur de pénétration augmente avec la racine carrée du temps :

[ \delta(t)\propto \sqrt{D_{\mathrm{eff}}t} ]

Elles sont utiles pour comprendre l’extraction initiale à partir de solides volumineux ou de couches superficielles minces, bien que des solutions pour particules finies soient nécessaires lorsque le front de diffusion atteint le centre de la particule.

Le flux initial décroît souvent comme (t^{-1/2})

Pour un solide semi-infini initialement uniforme exposé à une concentration de surface constante, le flux à la surface s’écrit :

[ N_A\propto \left(\frac{D_{\mathrm{eff}}}{\pi t}\right)^{1/2} \Delta C ]

La dépendance en (t^{-1/2}) décrit la diminution caractéristique du flux de diffusion transitoire. Dans un extracteur à volume fini, la diminution est généralement plus prononcée ou modifiée d’une autre manière, car la concentration du liquide varie et la différence de concentration disponible diminue.

Application du modèle aux installations pilotes

La durée d’extraction est une variable de conception couplée

Une installation pilote doit fournir un temps de contact suffisant pour que la fraction ciblée de soluté quitte le solide. Cette durée dépend de (D_{\mathrm{eff}}), du rayon des particules, de la température, des propriétés du solvant, du rapport liquide-solide, du comportement à l’équilibre et de l’agitation.

Une estimation simple fondée uniquement sur (R^2/D_{\mathrm{eff}}) est utile pour un dimensionnement initial, mais elle ne peut pas prédire la récupération finale lorsque l’accumulation dans le liquide et les limitations d’équilibre sont importantes.

L’agitation contrôle principalement le transport externe

L’agitation réduit la couche limite liquide entourant les particules et augmente le coefficient de transfert de matière du côté liquide (k_L). Elle peut également maintenir les particules en suspension, empêcher leur sédimentation, améliorer l’uniformité de la température et réduire les gradients de concentration dans la cuve.

Une fois que la résistance externe devient faible, une agitation supplémentaire apporte des améliorations décroissantes, car la diffusion interne dans le solide est devenue la résistance dominante.

La taille des particules influence fortement la capacité de traitement

La réduction de la taille des particules raccourcit la longueur de diffusion caractéristique et augmente la surface externe disponible par unité de masse. Ces changements accélèrent généralement l’extraction, mais peuvent également accroître la consommation d’énergie, créer des fines, compliquer la filtration et provoquer des problèmes de manipulation ou de sécurité.

Des essais pilotes sont donc nécessaires pour déterminer la plus petite taille de particule pratique, plutôt que de supposer qu’un broyage maximal est optimal.

Le volume de liquide contrôle la force motrice

Un volume de solvant plus important peut maintenir la concentration du liquide en vrac plus faible pour une quantité donnée de soluté extrait. Cela préserve une force motrice plus élevée et peut augmenter la récupération finale lorsque le partage à l’équilibre limiterait autrement l’extraction.

Le compromis réside dans une consommation accrue de solvant, une charge de concentration plus importante en aval et, potentiellement, un volume d’équipement plus grand. Le rapport liquide-solide approprié est donc une variable d’optimisation plutôt qu’une constante universelle.

Les données pilotes relient la théorie au changement d’échelle

Une campagne pilote utile mesure la concentration du liquide en vrac en fonction du temps, tout en enregistrant la taille des particules, la température, la vitesse d’agitation, le volume de solvant et la charge en solide. Ces données peuvent être ajustées à un modèle de diffusion transitoire afin d’estimer (D_{\mathrm{eff}}), (k_L), les paramètres d’équilibre ou des combinaisons de ces paramètres.

La seule vitesse d’agitation ne constitue pas un critère de changement d’échelle transférable. Le changement d’échelle doit préserver le comportement de transport déterminant, en utilisant des grandeurs pertinentes telles que le temps de mélange, la puissance fournie, la qualité de la suspension, le nombre de Reynolds ou des coefficients de transfert de matière corrélés expérimentalement.

Comprendre les compromis

L’hypothèse d’un mélange parfait a ses limites

Un liquide peut sembler vigoureusement mélangé tout en présentant des gradients locaux de concentration près des particules, des parois de la cuve ou dans des régions mal mises en suspension. Si le volume de liquide est important ou si la suspension est dense, l’hypothèse d’un (C_L(t)) spatialement uniforme peut être inexacte.

Des mesures du temps de mélange, des essais avec traceur ou des modèles compartimentés et de mécanique des fluides numérique peuvent être utilisés lorsque l’uniformité globale est incertaine.

La diffusion interne peut ne pas être le seul mécanisme

La vitesse d’extraction apparente peut inclure un lavage rapide de la surface, une diffusion dans les pores, un gonflement, une rupture des particules, une dissolution, une réaction chimique, une adsorption ou une désorption. Considérer toute déviation par rapport au comportement fickien idéal comme une variation de la diffusivité peut conduire à des paramètres ajustés trompeurs.

Un modèle de réaction-diffusion ou un modèle de particule rétrécissante peut être nécessaire lorsque la structure du solide évolue pendant l’extraction.

La diffusivité effective est généralement une propriété apparente

(D_{\mathrm{eff}}) n’est pas nécessairement la diffusivité moléculaire mesurée dans un liquide pur. Elle intègre souvent la géométrie des pores, la tortuosité, la constriction, la saturation, le partage et les interactions avec la matrice solide.

Sa valeur peut varier avec la température, la teneur en humidité, la concentration du soluté ou le degré d’épuisement. Un coefficient ajusté à l’échelle du laboratoire doit donc être validé dans des conditions pilotes.

L’équilibre peut limiter la récupération finale

Une augmentation de l’agitation ne peut pas surmonter une relation d’équilibre défavorable. Si le liquide approche sa concentration d’équilibre, la vitesse d’extraction devient faible, même lorsque le mélange est excellent.

Le fonctionnement à contre-courant, le renouvellement du solvant, l’extraction par étapes, l’ajustement du pH ou une réaction chimique sélective peuvent améliorer la force motrice globale lorsque ces approches sont compatibles avec le produit et le procédé.

Les systèmes multicomposants nécessitent une attention supplémentaire

Pour une extraction diluée à soluté unique, un modèle fickien ou pseudo-binaire est souvent suffisant. Dans les systèmes concentrés ou multicomposants, les interactions entre espèces, l’adsorption compétitive, la thermodynamique non idéale et la diffusion couplée peuvent nécessiter des descriptions de transport de Stefan-Maxwell ou d’autres modèles multicomposants.

Le modèle doit correspondre à la plage de composition et au niveau de précision requis pour la décision concernant l’installation pilote.

Comment appliquer cela à votre projet

Un modèle d’extraction transitoire doit être développé à partir de propriétés mesurables et validé par rapport aux données de concentration en fonction du temps. Les priorités suivantes contribuent à aligner le modèle sur l’objectif d’ingénierie :

Si votre objectif principal est la durée d’extraction : Estimez (R^2/D_{\mathrm{eff}}), mesurez les distributions granulométriques et ajustez le modèle complet de diffusion transitoire dans le solide afin de déterminer la durée de traitement nécessaire pour atteindre la récupération requise.
Si votre objectif principal est la performance de l’agitation : Faites varier la vitesse de l’agitateur ou les conditions de fluidisation tout en mesurant la vitesse d’extraction, la qualité de la suspension et le temps de mélange afin de distinguer la résistance du côté liquide de la diffusion interne.
Si votre objectif principal est la consommation de solvant : Utilisez le bilan du liquide à volume fini et la relation d’équilibre pour choisir le rapport liquide-solide ou le nombre d’étapes d’extraction permettant d’atteindre la récupération souhaitée avec une consommation de solvant acceptable.
Si votre objectif principal est le changement d’échelle du pilote à la production : Déterminez si la diffusion interne, le transfert de matière externe, le mélange ou l’équilibre contrôle le procédé, puis préservez le comportement adimensionnel déterminant et validez-le à plus grande échelle.
Si votre objectif principal est la fiabilité du modèle : Vérifiez la conservation globale de la matière, mesurez la température et la viscosité, caractérisez la porosité et la taille des particules et vérifiez que (D_{\mathrm{eff}}) et (k_L) restent valables dans la plage opératoire.

La tâche centrale de l’ingénierie consiste à identifier la résistance déterminante et à concevoir l’installation pilote de façon que ses mesures permettent de distinguer les effets de la diffusion, du mélange et de l’équilibre.

Tableau récapitulatif :

Principe	Rôle dans la conception de l’extraction
Première loi de Fick	Calcule le flux de soluté à partir des gradients de concentration.
Deuxième loi de Fick	Prédit l’épuisement transitoire du soluté à l’intérieur des particules.
Bilan du liquide fini	Suit l’augmentation de la concentration globale et la diminution de la force motrice.
Nombre de Fourier	Relie l’avancement de la diffusion au temps et à la taille des particules.
Nombre de Biot de transfert de matière	Compare la résistance externe du liquide à la diffusion interne.
Rayon des particules	Influence fortement la durée d’extraction, approximativement proportionnelle à R^2/D_eff.
Agitation	Augmente le transfert de matière du côté liquide, la mise en suspension et l’uniformité de la température.
Relation d’équilibre	Détermine la récupération maximale pouvant être atteinte.
Données pilotes	Valident D_eff, k_L, les paramètres d’équilibre et les hypothèses de changement d’échelle.

Améliorez votre installation pilote d’extraction

Vous devez déterminer la durée d’extraction, le volume de solvant, les besoins en agitation ou les conditions de changement d’échelle ? Contactez LABPARK pour discuter d’une installation pilote d’opérations unitaires à vocation éducative ou professionnelle, adaptée à votre projet.

LABPARK aide les universités, les instituts de recherche et les entreprises à évaluer l’extraction solide-liquide à l’aide d’équipements pilotes pratiques destinés au génie chimique, aux bioprocédés et aux biotechnologies, ainsi qu’aux applications environnementales et de traitement de l’eau. Nos systèmes permettent d’étudier de manière contrôlée la taille des particules, la température, le rapport liquide-solide, l’agitation, le mélange, le transfert de matière et l’équilibre, afin que votre équipe puisse identifier la résistance déterminante et générer des données fiables de concentration en fonction du temps.

Travaillez avec LABPARK pour relier la théorie de la diffusion transitoire aux performances mesurables du procédé, améliorer l’efficacité de l’utilisation des solvants et de l’énergie, valider les paramètres de changement d’échelle et sélectionner des conditions opératoires permettant d’atteindre votre récupération cible.