Les fonctions de Bessel modifiées fournissent le profil radial exact de concentration pour une diffusion stationnaire avec une réaction du premier ordre dans un cylindre solide. Pour une pastille de catalyseur cylindrique ou un réacteur tubulaire modélisé radialement, la concentration adimensionnelle du réactif est généralement

[ \Theta(\rho)=\frac{C_A(r)}{C_{As}} =\frac{I_0(\phi \rho)}{I_0(\phi)}, \qquad \rho=\frac{r}{R}, ]

où (I_0) est la fonction de Bessel modifiée de première espèce, (R) est le rayon du cylindre et (\phi=R\sqrt{k/D_e}) est le module de Thiele cylindrique. La solution correspondante de diffusion-réaction met en évidence les gradients de concentration internes, les limitations de diffusion et l’efficacité du catalyseur.

Idée essentielle : Dans un cylindre solide, le bilan de diffusion-réaction du premier ordre produit une équation de Bessel modifiée. La solution physiquement acceptable utilise (I_0), car celle-ci reste finie au centre ; le profil obtenu et le facteur d’efficacité montrent dans quelle mesure la diffusion intraparticulaire réduit la vitesse de réaction observable.

Formation du modèle cylindrique de diffusion-réaction

Bilan stationnaire sur une coquille

Considérons une pastille cylindrique de coordonnée radiale (r), de rayon (R), de diffusivité effective (D_e) et siège d’une réaction du premier ordre

[ -r_A=kC_A. ]

À l’état stationnaire, la diffusion radiale vers l’intérieur de la pastille doit équilibrer la consommation du réactif. En supposant (D_e) constant, aucune variation axiale et aucune variation de la vitesse de réaction autre que celle due à (C_A), l’équation gouvernante est

[ D_e\left[\frac{1}{r}\frac{d}{dr} \left(r\frac{dC_A}{dr}\right)\right]-kC_A=0. ]

De manière équivalente,

[ \frac{d^2C_A}{dr^2} +\frac{1}{r}\frac{dC_A}{dr} -\frac{k}{D_e}C_A=0. ]

Pourquoi l’équation est une équation de Bessel modifiée

Introduisons le rayon adimensionnel

[ \rho=\frac{r}{R} ]

et définissons le module de Thiele cylindrique

[ \phi=R\sqrt{\frac{k}{D_e}}. ]

L’équation devient

[ \frac{d^2\Theta}{d\rho^2} +\frac{1}{\rho}\frac{d\Theta}{d\rho} -\phi^2\Theta=0, ]

où

[ \Theta=\frac{C_A}{C_{As}}. ]

Il s’agit de l’équation de Bessel modifiée d’ordre zéro.

Application des conditions aux limites

Symétrie sur l’axe central

Au centre du cylindre, aucune direction radiale n’est privilégiée. Par conséquent,

[ \left.\frac{dC_A}{dr}\right|_{r=0}=0. ]

La concentration doit également rester finie en (r=0).

Concentration à la surface

Pour une pastille dont la résistance externe au transfert de matière est négligée, la concentration à la surface est imposée :

[ C_A(R)=C_{As}, ]

ou, sous forme adimensionnelle,

[ \Theta(1)=1. ]

Si la résistance du film fluide externe est importante, la concentration à la surface n’est pas connue à l’avance et doit être déterminée à partir d’une condition d’égalité des flux.

Sélection de la fonction de Bessel physiquement valide

La solution générale est

[ \Theta(\rho)=A I_0(\phi\rho)+B K_0(\phi\rho), ]

où (I_0) et (K_0) sont les fonctions de Bessel modifiées de première et de deuxième espèces.

La fonction (K_0) diverge lorsque son argument tend vers zéro. Elle est donc inadmissible pour un cylindre solide contenant l’axe central, qui exige une concentration finie. En posant (B=0) et en appliquant la condition à la surface, on obtient

[ \boxed{ \Theta(\rho)= \frac{I_0(\phi\rho)}{I_0(\phi)} } ]

et donc

[ \boxed{ C_A(r)=C_{As} \frac{I_0\left(r\sqrt{k/D_e}\right)} {I_0\left(R\sqrt{k/D_e}\right)} }. ]

La concentration est minimale au centre et maximale à la surface, car la réaction consomme le réactif dans toute la pastille.

Ce que révèle le profil de concentration

Faible limitation par la diffusion

Lorsque (\phi\ll 1), la diffusion est rapide par rapport à la réaction. La concentration reste presque uniforme :

[ C_A(r)\approx C_{As}. ]

La pastille fonctionne à une vitesse proche de sa vitesse cinétique intrinsèque, de sorte qu’une augmentation de la diffusivité apporte peu de bénéfice supplémentaire.

Forte limitation par la diffusion

Lorsque (\phi\gg 1), la réaction est rapide par rapport à la diffusion radiale. Le réactif ne pénètre que sur une distance limitée depuis la surface, ce qui produit un fort gradient de concentration vers le centre.

Dans ce régime, une grande partie de l’intérieur de la pastille contribue peu à la vitesse globale de réaction, car sa concentration locale en réactif est faible.

Interprétation du nombre de Damköhler

Le groupe adimensionnel pertinent s’écrit souvent

[ \mathrm{Da}=\frac{kR^2}{D_e}=\phi^2. ]

Il compare la vitesse caractéristique de réaction à la vitesse caractéristique de diffusion. Une valeur élevée de (\mathrm{Da}) indique une limitation interne par la diffusion plus forte, tandis qu’une faible valeur indique un champ de concentration presque uniforme.

La valeur numérique précise dépend de la longueur caractéristique choisie et de la définition utilisée par le modèle du réacteur ; la géométrie et la convention doivent donc toujours être précisées.

Calcul du facteur d’efficacité

Définition

Le facteur d’efficacité interne compare la vitesse de réaction réelle intégrée sur la pastille à la vitesse qui serait obtenue si toute la pastille était à la concentration de surface :

[ \eta= \frac{\text{vitesse totale réelle de réaction}} {\text{vitesse à } C_A=C_{As}\text{ partout}}. ]

Pour une réaction du premier ordre dans un cylindre solide, l’intégration du profil obtenu avec les fonctions de Bessel donne

[ \boxed{ \eta= \frac{2I_1(\phi)} {\phi I_0(\phi)} } ]

où (I_1) est la fonction de Bessel modifiée de première espèce et d’ordre un.

Comportement aux limites

Pour (\phi\ll 1),

[ \eta\approx 1, ]

de sorte que la quasi-totalité de la pastille est cinétiquement active.

Pour (\phi\gg 1),

[ \eta\approx \frac{2}{\phi}, ]

ce qui montre que la vitesse de réaction effective diminue lorsque le module de Thiele augmente.

Relation entre le facteur et la vitesse mesurée

La vitesse totale de réaction par volume de pastille peut s’écrire

[ \overline{r}_A

\eta k C_{As}. ]

Ainsi, un modèle d’installation pilote peut utiliser (I_0) pour calculer le profil interne et (\eta) pour convertir la constante cinétique intrinsèque en une vitesse observable à l’échelle de la pastille.

Utilisation du modèle dans l’analyse d’une installation pilote

Estimation des gradients de concentration internes

À partir de (k), (D_e), (R) et (C_{As}), les ingénieurs peuvent calculer (\phi), évaluer (I_0) et déterminer (C_A(r)) dans toute la pastille.

Cela permet de déterminer si le centre de la pastille reste actif ou si la réaction se concentre dans la région externe.

Comparaison des géométries de catalyseurs

Le rayon intervient dans le module de Thiele par l’intermédiaire de (R\sqrt{k/D_e}). Réduire la dimension radiale caractéristique diminue (\phi), améliore la pénétration du réactif et augmente généralement le facteur d’efficacité.

C’est pourquoi la géométrie du catalyseur, la taille des pastilles et les supports structurés sont des variables de conception importantes dans les réacteurs à l’échelle pilote.

Distinction entre limitations internes et externes

La solution faisant intervenir les fonctions de Bessel décrit la diffusion et la réaction internes. Elle ne tient pas compte, à elle seule, de la résistance du film fluide environnant.

En présence d’un film externe, le modèle combine le facteur d’efficacité interne avec une condition de flux telle que

[ N_A=k_g(C_{Ab}-C_{As})

-D_e\left.\frac{dC_A}{dr}\right|_{r=R}, ]

où (C_{Ab}) est la concentration dans le fluide en vrac et (k_g) le coefficient externe de transfert de matière.

Cette distinction évite d’interpréter à tort une résistance externe du film comme une faible diffusivité à l’intérieur de la pastille.

Géométries et fonctions de Bessel associées

Réacteurs cylindriques creux et pastilles annulaires

Si le domaine est un cylindre creux ou un anneau, l’axe central est exclu. Les deux solutions indépendantes peuvent alors être conservées :

[ C_A(r)=A I_0\left(r\sqrt{k/D_e}\right) +B K_0\left(r\sqrt{k/D_e}\right). ]

Les coefficients sont déterminés à partir des conditions imposées aux rayons intérieur et extérieur.

Comparaison avec les pastilles sphériques

Pour une sphère, l’équation gouvernante contient le terme radial sphérique

[ \frac{1}{r^2}\frac{d}{dr} \left(r^2\frac{dC_A}{dr}\right) -\frac{k}{D_e}C_A=0. ]

Sa solution peut être exprimée à l’aide de fonctions hyperboliques, ce qui donne

[ \Theta(\rho)= \frac{\sinh(\phi\rho)} {\rho\sinh(\phi)} ]

dans le cas d’une concentration de surface imposée. Ce profil n’est pas identique à la solution cylindrique en (I_0), car la géométrie modifie la surface de diffusion disponible à chaque rayon.

Comprendre les compromis

Une pastille plus grande n’est pas toujours plus productive

L’augmentation du rayon de la pastille peut accroître la quantité de catalyseur par particule, mais elle augmente également (\phi) et rend la diffusion interne plus limitante.

La vitesse obtenue par volume de catalyseur peut alors être inférieure à celle attendue à partir de la seule cinétique intrinsèque.

La diffusivité effective dépend du modèle

La valeur de (D_e) intègre des effets tels que la structure des pores, la tortuosité, la porosité et, parfois, les hypothèses relatives à l’adsorption ou au transport multicomposant.

Utiliser une diffusivité moléculaire à la place d’une diffusivité effective dans un milieu poreux peut sous-estimer considérablement le gradient de concentration prévu.

Les hypothèses isothermes et de propriétés constantes peuvent échouer

Le profil standard suppose que (D_e) et (k) sont constants et que la température est uniforme. Des effets thermiques importants, une porosité variable, une désactivation, une cinétique non linéaire ou une diffusivité dépendant de la concentration nécessitent un modèle modifié.

Une concentration de surface fixe peut être irréaliste

Le profil normalisé en (I_0) suppose que (C_{As}) est connue. Dans un réacteur en fonctionnement, la concentration à la surface peut être réduite par un film externe ou par l’appauvrissement axial du fluide en vrac.

La solution interne doit donc être couplée aux bilans à l’échelle du réacteur lorsque ces effets sont significatifs.

Faire le bon choix selon votre objectif

Utilisez la solution cylindrique en fonctions de Bessel modifiées comme modèle de référence, puis étendez-la uniquement lorsque les conditions de fonctionnement exigent de prendre en compte des effets supplémentaires de transport ou de réaction.

Si votre objectif principal est d’estimer les profils de concentration internes : calculez (\phi=R\sqrt{k/D_e}) et utilisez (\Theta(\rho)=I_0(\phi\rho)/I_0(\phi)) avec des conditions aux limites de centre fini et de surface fixe.
Si votre objectif principal est de prévoir les performances du catalyseur : utilisez (\eta=2I_1(\phi)/[\phi I_0(\phi)]) pour convertir la cinétique intrinsèque du premier ordre en vitesse de réaction à l’échelle de la pastille.
Si votre objectif principal est de diagnostiquer les limitations de transfert de matière : comparez le module de Thiele interne à la résistance du film externe plutôt que d’attribuer toute perte de vitesse à la diffusion intraparticulaire.
Si votre objectif principal est d’analyser des équipements cylindriques creux : conservez à la fois (I_0) et (K_0), car l’exclusion de l’axe central n’élimine plus la solution (K_0).

La formulation à l’aide des fonctions de Bessel modifiées transforme la compétition entre la diffusion radiale et la réaction du premier ordre en une base quantitative pour interpréter et concevoir des systèmes cylindriques à l’échelle pilote.

Tableau récapitulatif :

Élément clé	Enseignement pratique
Modèle gouvernant	Diffusion radiale stationnaire avec une réaction du premier ordre
Profil de concentration	`Theta(rho) = I_0(phi rho) / I_0(phi)`
Module de Thiele	`phi = R sqrt(k / D_e)` ; des valeurs plus élevées indiquent une limitation par la diffusion plus forte
Solution physique	`I_0` est conservée, car elle reste finie au centre du cylindre ; `K_0` y diverge
Facteur d’efficacité	`eta = 2 I_1(phi) / [phi I_0(phi)]`
Conséquences pour la conception	Des pastilles plus petites et une diffusivité effective plus élevée améliorent la pénétration du réactif
Modèles étendus	La résistance du film externe, les géométries annulaires, les propriétés variables et les effets thermiques peuvent nécessiter des équations supplémentaires

Améliorer les performances des réacteurs et des catalyseurs à l’échelle pilote

Besoin d’évaluer la diffusion interne, les vitesses de réaction ou les risques liés au changement d’échelle dans des systèmes cylindriques ? Contactez LABPARK pour discuter de l’installation pilote d’opérations unitaires éducatives et professionnelles adaptée à votre application.

LABPARK accompagne les universités, les instituts de recherche et les entreprises avec des installations pilotes pour le génie chimique, les bioprocédés et les biotechnologies, ainsi que le traitement de l’environnement et de l’eau. Nos systèmes vous aident à étudier les limitations de transfert de matière, à valider les modèles cinétiques, à comparer les géométries de catalyseurs et à former des ingénieurs à l’aide de données pratiques d’exploitation. Indiquez-nous le type de réacteur, le catalyseur ou le procédé, les conditions de fonctionnement et vos objectifs d’essai : notre équipe pourra vous aider à identifier une solution appropriée à l’échelle pilote.