La base théorique est l'équation de diffusion combinée à la superposition des fonctions de Green. Une libération continue de soluté est représentée comme une séquence de libérations ponctuelles instantanées infinitésimales, chacune se propageant selon l'équation de diffusion. L'intégration de ces solutions instantanées sur l'historique de libération donne le champ de concentration transitoire, tandis que la limite à long terme produit la distribution en régime permanent lorsque la source reste constante et que les conditions environnantes permettent un flux stationnaire.
Une libération continue est l'intégrale temporelle des réponses de diffusion instantanées. Pour une libération à source ponctuelle constante dans un milieu tridimensionnel infini et stagnant, la solution transitoire approche (C(r)=\dot{m}/(4\pi Dr)) à mesure que le temps augmente.
Comment les libérations continues sont modélisées
L'équation de diffusion directrice
Pour la diffusion moléculaire avec une source ponctuelle, la concentration est régie par la deuxième loi de Fick :
[ \frac{\partial C}{\partial t}=D\nabla^2 C ]
Ici, (C) est la concentration, (t) est le temps, (D) est la diffusivité moléculaire, et (\nabla^2 C) décrit la courbure spatiale du champ de concentration.
Pour une source située à l'origine, la symétrie sphérique réduit le problème à la distance radiale (r). La concentration dépend alors de la manière dont le soluté se propage sur les surfaces sphériques entourant la source.
La solution de la source ponctuelle instantanée
Une libération instantanée de masse (m) produit une réponse de diffusion décrite par une fonction de Green gaussienne dans un milieu tridimensionnel non borné :
[ C_{\text{inst}}(r,t)= \frac{m}{(4\pi Dt)^{3/2}} \exp\left(-\frac{r^2}{4Dt}\right) ]
Cette solution représente une masse finie libérée à un instant donné. Elle se propage vers l'extérieur avec le temps, tandis que sa concentration maximale diminue car la masse occupe un volume de plus en plus grand.
Superposition sur le temps de libération
La diffusion est un processus de transport linéaire selon les hypothèses énoncées. Par conséquent, la réponse à de nombreuses libérations est égale à la somme des réponses aux libérations individuelles.
Pour un taux de libération (\dot{m}(\tau)), où (\tau) est le temps auquel chaque incrément de masse infinitésimal est libéré, la concentration est obtenue en intégrant la solution instantanée :
[ C(r,t)= \int_0^t \frac{\dot{m}(\tau)} {\left[4\pi D(t-\tau)\right]^{3/2}} \exp\left[ -\frac{r^2}{4D(t-\tau)} \right]d\tau ]
Cette intégrale est le fondement théorique pour calculer les distributions de concentration transitoires à partir de libérations continues.
La distribution de concentration transitoire
Libération à taux constant
Pour une source qui commence à (t=0) et libère de la masse en continu à un taux constant (\dot{m}), l'évaluation de l'intégrale de superposition donne :
[ C(r,t)= \frac{\dot{m}}{4\pi Dr} \operatorname{erfc} \left( \frac{r}{2\sqrt{Dt}} \right) ]
La fonction d'erreur complémentaire, (\operatorname{erfc}), décrit le développement du front de diffusion à une distance (r) après un temps (t).
Signification physique du terme de la fonction d'erreur
Le facteur
[ \operatorname{erfc} \left( \frac{r}{2\sqrt{Dt}} \right) ]
est une correction dépendante du temps au profil radial permanent. Il tient compte du fait que le soluté n'a eu qu'un temps fini pour diffuser de la source vers le lieu d'observation.
Aux premiers instants, (r/(2\sqrt{Dt})) est grand à une distance donnée, donc la fonction d'erreur complémentaire est petite. La concentration est donc bien inférieure à sa valeur finale en régime permanent.
Échelle de temps de diffusion
La comparaison pertinente se fait entre le temps écoulé (t) et le temps de diffusion caractéristique :
[ t_D\sim\frac{r^2}{D} ]
Lorsque (t) est petit par rapport à (r^2/D), la mesure est fortement transitoire. Lorsque (t) est grand par rapport à (r^2/D), la concentration à cette distance est proche de sa valeur en régime permanent.
La distribution de concentration en régime permanent
Prise de la limite à long terme
Lorsque (t) approche l'infini,
[ \frac{r}{2\sqrt{Dt}}\rightarrow 0 ]
et
[ \operatorname{erfc}(0)=1 ]
L'expression transitoire devient donc :
[ C(r)=\frac{\dot{m}}{4\pi Dr} ]
Il s'agit de la distribution de concentration en régime permanent pour une libération ponctuelle constante dans un milieu tridimensionnel infini.
Interprétation du flux
Le résultat en régime permanent découle directement de la conservation de la masse. Au rayon (r), le taux de libération total doit traverser une surface sphérique d'aire (4\pi r^2).
Le flux diffusif radial est donc :
[ J_r=\frac{\dot{m}}{4\pi r^2} ]
En utilisant la première loi de Fick,
[ J_r=-D\frac{dC}{dr} ]
et en intégrant avec la condition que la concentration approche zéro loin de la source, on obtient :
[ C(r)=\frac{\dot{m}}{4\pi Dr} ]
La dépendance en (1/r) reflète la dilution géométrique sur des surfaces sphériques en expansion.
Pourquoi le régime permanent n'est pas un équilibre de masse finie
Un régime permanent ne signifie pas que le soluté a cessé de se déplacer. Cela signifie que la concentration à chaque emplacement ne change plus avec le temps car le flux diffusif sortant est constant.
La source continue d'ajouter de la masse, et cette masse doit quitter la région modélisée ou s'étendre indéfiniment dans le milieu environnant. Sans exutoire, puits limite, réaction ou autre mécanisme d'élimination, un véritable régime permanent peut ne pas être physiquement atteignable dans un système pilote fermé.
Interprétation des mesures d'expériences pilotes
Comparaison du temps d'observation avec la distance
La question de savoir si une mesure doit être traitée comme transitoire ou en régime permanent dépend à la fois du temps écoulé et de la distance par rapport au point de libération.
Un emplacement proche de la source peut approcher son profil permanent relativement rapidement, tandis qu'un emplacement plus éloigné peut rester transitoire pendant la même expérience.
Estimation du régime de transport pertinent
La quantité adimensionnelle
[ \eta=\frac{r}{2\sqrt{Dt}} ]
fournit un diagnostic direct.
Un grand (\eta) indique que le front de diffusion ne s'est pas complètement développé à l'emplacement de mesure. Un petit (\eta) indique que la correction transitoire est faible et que l'approximation du régime permanent est de plus en plus appropriée.
Utilisation du modèle transitoire complet
L'équation transitoire doit être utilisée lorsque les données de concentration sont collectées pendant le démarrage de la source, lorsque la durée de libération est comparable au temps de diffusion, ou lorsque les emplacements spatiaux sont suffisamment éloignés pour que le front de diffusion n'ait pas atteint une condition proche du régime permanent.
C'est également le modèle approprié lorsque l'objectif est d'estimer le moment de percée, l'accumulation de concentration ou la vitesse à laquelle le système pilote approche son profil stationnaire.
Comprendre les compromis
Hypothèses derrière les équations
Les solutions énoncées supposent une source ponctuelle, une diffusivité constante, une symétrie sphérique, un milieu infini ou suffisamment grand, et aucune advection, réaction, sorption ou effet de bord.
L'équipement pilote réel peut violer plusieurs de ces hypothèses. L'écoulement, les parois des cuves, le mélange, la géométrie finie et les propriétés non uniformes des matériaux peuvent modifier matériellement la distribution de concentration.
La signification du « régime permanent » dans une unité pilote
L'équation (C(r)=\dot{m}/(4\pi Dr)) décrit un champ permanent idéalisé contrôlé par la diffusion. Elle ne doit pas être automatiquement interprétée comme l'état final d'un réacteur fini ou d'une cuve expérimentale.
Si les limites réfléchissent le soluté, si le système est fermé ou si la libération s'arrête, le comportement à long terme peut différer de la solution pour source ponctuelle en milieu infini.
Libération continue versus libération par impulsion
Une libération continue accumule les effets de toutes les libérations antérieures, ce qui produit la solution de la fonction d'erreur. Une libération par impulsion suit plutôt la solution gaussienne instantanée et présente généralement un maximum de concentration mobile plutôt qu'un profil radial stationnaire.
Appliquer l'équation de source continue à une impulsion finie peut donc fausser à la fois l'amplitude de la concentration et la dépendance temporelle.
Diffusion versus advection
Les formules décrivent uniquement la diffusion moléculaire. Si le mouvement global du fluide est significatif, le transport doit également inclure un terme d'advection, et le champ de concentration peut être régi par une équation d'advection-diffusion plutôt que par l'équation de diffusion pure.
Faire le bon choix pour votre objectif
Utilisez le modèle qui correspond à la question physique et à l'échelle de temps de l'expérience pilote.
- Si votre objectif principal est le comportement transitoire : Utilisez (C(r,t)=\frac{\dot{m}}{4\pi Dr}\operatorname{erfc}\left(\frac{r}{2\sqrt{Dt}}\right)) et comparez le temps écoulé avec (r^2/D).
- Si votre objectif principal est la distribution radiale à long terme : Utilisez (C(r)=\frac{\dot{m}}{4\pi Dr}) uniquement après avoir confirmé que les emplacements pertinents ont atteint une condition approximativement permanente.
- Si votre objectif principal est la validité du modèle : Vérifiez si l'advection, les limites, les réactions, les dimensions finies de la source ou les taux de libération non constants nécessitent un modèle de transport plus complet.
- Si votre objectif principal est l'interprétation des mesures : Ajustez d'abord la solution transitoire lorsque les effets de démarrage ou de percée sont visibles, puis évaluez si les données supportent une approximation en régime permanent.
La distinction est régie par le même principe tout au long : les champs de concentration à libération continue sont des réponses de diffusion intégrées dans le temps, et le régime permanent est leur limite à flux constant à long terme.
Tableau récapitulatif :
| Modèle ou concept | Équation clé | Meilleure utilisation pour |
|---|---|---|
| Équation de diffusion | ∂C/∂t = D∇²C | Décrire la diffusion moléculaire |
| Libération instantanée | Fonction de Green gaussienne | Modéliser des libérations par impulsion ou de masse finie |
| Libération continue | C(r,t) = [ṁ/(4πDr)] erfc[r/(2√Dt)] | Prédire l'accumulation transitoire et la percée |
| Libération en régime permanent | C(r) = ṁ/(4πDr) | Estimer les profils de concentration radiaux à long terme |
| Échelle de temps de diffusion | tD ≈ r²/D | Déterminer si les mesures sont transitoires ou proches du régime permanent |
| Vérifications de validité | Considérer l'advection, les limites, les réactions, la sorption et la géométrie de la source | Sélectionner un modèle de transport plus complet si nécessaire |
Modélisez votre processus de diffusion avec LABPARK
Besoin de données fiables pour l'accumulation de concentration transitoire, le timing de percée ou le transport en régime permanent ? Contactez LABPARK pour discuter d'une unité pilote d'opérations unitaires éducative ou professionnelle adaptée à votre application. LABPARK soutient les universités, les instituts de recherche et les entreprises avec des usines pilotes pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que le traitement de l'eau et de l'environnement. Nos systèmes vous aident à étudier la diffusion, le transfert de masse, le comportement d'advection-diffusion, les effets de bord et le passage à l'échelle des processus en utilisant des expériences pratiques et contrôlées.
Produits associés
- Pilote d'Unités Opératoires Pédagogique pour la Mesure du Facteur Efficace de Diffusion Intraparticulaire
- Usine pilote pédagogique de distillation spéciale multifonctionnelle
- Installation pilote éducative pour la distillation continue, discontinue et extractive
- Installation pilote de rectification bimodale pour la formation pratique aux opérations unitaires
- Installation pilote éducative pour la détermination de la distribution des temps de séjour et des performances de mélange de cuves agitées en série
Les gens demandent aussi
- Comment les installations pilotes d'opérations unitaires éducatives combinent-elles théorie et conception ? Combler le fossé de l'ingénierie
- En quoi les exigences en matière d'opérations unitaires diffèrent-elles entre les procédés chimiques traditionnels et les bioprocédés ? Construire de meilleurs programmes pour les usines pilotes
- Comment les installations pilotes d'opérations unitaires éducatives abordent-elles la sécurité et la gestion des déchets lors de la mise à l'échelle ?
- Comment les usines pilotes d'opérations unitaires éducatives font-elles le pont entre la théorie et le développement des procédés industriels ? Apprentissage pratique
- Comment les usines pilotes éducatives enseignent-elles la mise à l'échelle ? Combler le fossé vers la fabrication commerciale.
