La méthode analytique consiste à résoudre le champ d’écoulement de Stokes autour de l’objet en rotation, à calculer la contrainte de cisaillement à la surface, puis à intégrer cette contrainte afin d’obtenir le couple. Pour une sphère immergée de rayon (R) tournant avec une vitesse angulaire (\omega) dans un fluide newtonien de viscosité (\mu), la solution à faible nombre de Reynolds donne le couple résistant

[ G=-8\pi\mu\omega R^3 ]

où le signe négatif indique que le couple hydrodynamique s’oppose à la rotation.

En régime rampant, le couple est entièrement déterminé par les contraintes visqueuses : il faut résoudre les équations linéaires de Stokes avec la condition aux limites imposée par la surface en rotation, évaluer la contrainte tangentielle à la surface et intégrer son moment sur l’objet.

Pourquoi la méthode de l’écoulement de Stokes s’applique

Inertie du fluide négligeable

La méthode s’applique lorsque le nombre de Reynolds est suffisamment faible pour que les termes inertiels puissent être négligés devant les termes visqueux. L’équation de quantité de mouvement devient

[ -\nabla p+\mu\nabla^2\mathbf{v}=0, \qquad \nabla\cdot\mathbf{v}=0 ]

Pour une sphère en rotation dans un fluide par ailleurs immobile, l’écoulement est axisymétrique et purement azimutal.

Relation linéaire entre le couple et la vitesse

Comme les équations de Stokes sont linéaires, le couple est directement proportionnel à la vitesse angulaire. Il est également proportionnel à la viscosité du fluide et au cube du rayon de la sphère :

[ |G|=8\pi\mu\omega R^3 ]

Cette proportionnalité est l’analogue rotationnel de la loi linéaire de traînée pour une sphère en translation.

Résolution du champ d’écoulement rotationnel

Appliquer les conditions aux limites

Utilisez les coordonnées sphériques avec la sphère centrée à l’origine. Pour une sphère tournant autour de son axe, les conditions de vitesse sont les suivantes :

[ v_r=0, \qquad v_\theta=0, \qquad v_\phi(R,\theta)=\omega R\sin\theta ]

La dernière condition représente la vitesse tangentielle de la surface de la sphère. Loin de la sphère, le fluide est immobile :

[ v_\phi(r,\theta)\rightarrow 0 \quad \text{as} \quad r\rightarrow\infty ]

Déterminer la vitesse azimutale

La résolution des équations de Stokes soumises à ces conditions donne

[ v_\phi(r,\theta)

\frac{\omega R^3\sin\theta}{r^2} ]

La vitesse décroît selon (1/r^2), ce qui reflète le caractère localisé de la perturbation rotationnelle en régime rampant.

Interpréter la distribution de l’écoulement

La vitesse du fluide est maximale à l’équateur de la sphère, où (\sin\theta=1), et s’annule sur l’axe de rotation. L’ensemble du champ d’écoulement est généré par le couplage visqueux entre la surface en rotation et le fluide environnant.

Calcul de la contrainte de cisaillement à la surface

Évaluer la contrainte tangentielle

Pour un écoulement azimutal axisymétrique, la composante pertinente du tenseur des contraintes en coordonnées sphériques est

[ \tau_{r\phi}

\mu \left( \frac{\partial v_\phi}{\partial r}

\frac{v_\phi}{r} \right) ]

En substituant le champ de vitesse et en évaluant à (r=R), on obtient

[ \tau_{r\phi}(R,\theta)

-3\mu\omega\sin\theta ]

Le signe négatif montre que la contrainte du fluide s’oppose à la rotation de la sphère.

Prendre en compte le bras de levier du couple

Un élément de surface situé à l’angle polaire (\theta) se trouve à une distance perpendiculaire de l’axe de rotation égale à

[ R\sin\theta ]

Le couple différentiel est donc égal à la force tangentielle exercée sur l’élément de surface, multipliée par ce bras de levier.

Intégration du couple

Intégrer sur la sphère

L’élément de surface est

[ dA=R^2\sin\theta,d\theta,d\phi ]

Ainsi, le couple total est

[ G

\int_0^{2\pi} \int_0^\pi \tau_{r\phi}(R,\theta) (R\sin\theta) R^2\sin\theta ,d\theta,d\phi ]

En substituant la contrainte de cisaillement, on obtient

[ G

-3\mu\omega R^3 \int_0^{2\pi}d\phi \int_0^\pi\sin^3\theta,d\theta ]

En utilisant

[ \int_0^{2\pi}d\phi=2\pi, \qquad \int_0^\pi\sin^3\theta,d\theta=\frac{4}{3} ]

on obtient

[ \boxed{G=-8\pi\mu\omega R^3} ]

Le couple appliqué requis a la même magnitude et un signe opposé :

[ \boxed{G_{\text{applied}}=8\pi\mu\omega R^3} ]

Relier le couple à la puissance

La puissance visqueuse nécessaire pour maintenir la rotation est

[ P=|G|\omega ]

Pour la sphère en rotation,

[ \boxed{P=8\pi\mu\omega^2R^3} ]

Cela fournit une estimation directe de la puissance mécanique dissipée dans le fluide.

Étendre la méthode aux opérations unitaires

Utiliser des coefficients de résistance rotationnelle

Pour d’autres corps immergés, la même procédure s’exprime au moyen d’un coefficient de résistance rotationnelle :

[ G=-\zeta_r\omega ]

Pour une sphère,

[ \zeta_r=8\pi\mu R^3 ]

Le coefficient varie selon la forme de l’objet, le confinement, la proximité de parois et la présence d’autres surfaces en mouvement.

Appliquer le cadre aux agitateurs

Dans les opérations unitaires à faible vitesse et à forte viscosité, le couple requis par une turbine ou un agitateur rotatif peut être obtenu en résolvant le champ d’écoulement rampant ou en déterminant un coefficient de résistance approprié. Le couple obtenu permet de calculer la contrainte de cisaillement, la dissipation de puissance et les besoins du moteur.

Distinguer les grandeurs locales et totales

La distribution des contraintes décrit la charge hydrodynamique locale sur l’ensemble de la surface. Le couple intégré est la grandeur globale nécessaire au dimensionnement des équipements et aux calculs de puissance.

Comprendre les compromis

L’hypothèse de faible nombre de Reynolds est restrictive

Le résultat analytique suppose que l’inertie du fluide est négligeable. Si le nombre de Reynolds n’est plus faible, l’inertie convective modifie les champs de vitesse et de contrainte, et la relation (G\propto\mu\omega R^3) peut ne plus être suffisante.

Le résultat pour une sphère n’est pas universel

Le coefficient (8\pi) s’applique spécifiquement à une sphère isolée en rotation dans un fluide newtonien non borné. Il ne doit pas être appliqué directement aux turbines, aux cylindres, aux particules proches de parois ou aux cuves confinées.

Le comportement newtonien est supposé

La dérivation suppose une viscosité constante et un comportement constitutif newtonien. Les solutions de polymères, les suspensions, les fluides biologiques et autres matériaux non newtoniens peuvent nécessiter un modèle de contrainte modifié et une relation différente entre le couple et la vitesse.

Les effets de paroi peuvent être importants

Les parois et les surfaces proches modifient le champ de vitesse et augmentent la résistance hydrodynamique. Dans les équipements pilotes, le confinement peut donc nécessiter des facteurs de correction de paroi, des solutions numériques de l’écoulement de Stokes ou un étalonnage expérimental.

Faire le bon choix selon votre objectif

La solution analytique pour une sphère constitue une référence et un premier calcul pratique, mais la géométrie et le régime d’écoulement doivent être vérifiés avant de l’appliquer à un équipement.

Si votre objectif principal est d’estimer le couple d’une sphère isolée en rotation : Utilisez (G=-8\pi\mu\omega R^3), le signe représentant la résistance à la rotation imposée.
Si votre objectif principal est de dimensionner un agitateur à faible vitesse : Déterminez la résistance rotationnelle spécifique à la géométrie de l’agitateur plutôt que d’appliquer directement le coefficient d’une sphère.
Si votre objectif principal est de prédire le cisaillement local : Utilisez le champ de vitesse et la distribution des contraintes à la surface avant d’intégrer pour obtenir le couple total.
Si votre objectif principal est de fonctionner en dehors du régime rampant : Prenez en compte l’inertie du fluide et utilisez une corrélation, un modèle numérique ou une expérience approprié(e) à un nombre de Reynolds plus élevé.

Le principe directeur est simple : résoudre le champ d’écoulement visqueux, évaluer la contrainte à la surface et intégrer son moment afin d’obtenir le couple hydrodynamique.

Tableau récapitulatif :

Aspect clé	Résultat ou application
Modèle directeur	Équations linéaires de Stokes en régime rampant
Couple exercé sur une sphère	`G = -8πμωR³`
Résistance rotationnelle	`ζᵣ = 8πμR³` pour une sphère isolée
Puissance requise	`P = 8πμω²R³`
Méthode de calcul	Résoudre l’écoulement, évaluer la contrainte à la surface et intégrer son moment
Principales limites	La forme, les parois, l’inertie et les fluides non newtoniens nécessitent des corrections ou des modèles différents

Vous avez besoin d’équipements fiables de mélange à faible vitesse et de traitement des fluides pour votre application ? Contactez LABPARK pour discuter de vos besoins. LABPARK fournit des installations pilotes d’opérations unitaires à vocation éducative et professionnelle pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que le traitement de l’environnement et de l’eau. Nos solutions aident les universités, les instituts de recherche et les entreprises à évaluer le couple, le cisaillement, la consommation d’énergie et les performances des équipements grâce à des systèmes pilotes axés sur les applications, conçus pour la formation pratique, la recherche, le développement des procédés et la mise à l’échelle.