L’équation de Stefan est la relation issue de la théorie de la lubrification entre la force de compression, la viscosité, la taille des disques, l’entrefer et la vitesse de rapprochement. Pour deux disques circulaires parallèles de rayon (R) se rapprochant l’un de l’autre, l’équation, selon la convention du demi-entrefer (H(t)), est

[ F(t)=\frac{3\pi \mu U(t)R^4}{8H(t)^3}, ]

où (F) est la force normale, (\mu) la viscosité dynamique, (U=-dH/dt) la vitesse de rapprochement de chaque disque, et (H) le demi-écartement entre les disques. Dans une unité pilote pédagogique, la mesure de (F), (H) et (U) permet aux étudiants de tester la théorie de la lubrification et d’estimer la viscosité à partir des données force-déplacement.

Le résultat central est la forte dépendance inverse au cube de l’entrefer : la réduction du demi-entrefer augmente la force de compression requise selon (H^{-3}), tandis que l’augmentation du rayon des disques accroît la force selon (R^4). L’équation est fiable lorsque l’entrefer est étroit, l’écoulement lent et les effets visqueux dominants par rapport à l’inertie.

Ce que décrit l’équation de Stefan

L’expérience physique

L’expérience utilise deux disques circulaires parallèles séparés par une fine couche de fluide. Les disques se rapprochent, forçant le liquide à s’écouler radialement vers l’extérieur à travers l’entrefer qui se rétrécit.

Le dispositif représente une opération simplifiée d’écoulement par compression, pertinente pour la lubrification, la transformation des polymères, le revêtement, le formage et la mesure rhéologique.

Les variables clés

Les variables sont :

(R) : rayon du disque
(H(t)) : demi-écartement entre les disques
(U(t)) : vitesse de rapprochement de chaque disque
(F(t)) : force normale appliquée
(\mu) : viscosité dynamique
(\rho) : masse volumique du fluide
(r) : distance radiale par rapport au centre du disque

L’épaisseur totale de la couche de fluide est (2H). Cette convention est importante, car d’autres références définissent souvent l’entrefer total comme (h=2H), ce qui produit un coefficient numérique d’apparence différente.

Les hypothèses fondamentales

La dérivation repose sur les hypothèses suivantes :

Écoulement rampant, avec une inertie négligeable
Un entrefer étroit, (H_0/R \ll 1)
Un comportement newtonien et incompressible du fluide
L’absence de glissement aux surfaces des disques
Des disques parallèles et un écoulement axisymétrique
Une vitesse de rapprochement uniforme sur chaque disque
Des effets de bord et des fuites de fluide négligeables

Le nombre de Reynolds pertinent est

[ Re=\frac{H_0U_0}{\nu}\ll 1, ]

où (\nu=\mu/\rho) est la viscosité cinématique.

Comment l’équation est dérivée

Commencer par l’écoulement radial de compression

Lorsque les disques se rapprochent, le fluide ne peut pas traverser les disques solides. Il s’écoule donc vers l’extérieur dans la direction radiale.

Comme l’entrefer est faible par rapport au rayon des disques, les gradients de vitesse à travers l’entrefer dominent les gradients de vitesse dans la direction radiale. Il s’agit de l’approximation centrale de la théorie de la lubrification.

Réduire l’équation de quantité de mouvement

Dans les conditions d’écoulement rampant, les termes inertiels des équations de Navier-Stokes sont négligeables. Le bilan radial de quantité de mouvement se réduit à

[ 0=-\frac{\partial P}{\partial r} +\mu\frac{\partial^2 u_r}{\partial z^2}, ]

où (u_r) est la vitesse radiale et (z) la coordonnée à travers la couche de fluide.

La pression varie principalement avec le rayon, tandis que la vitesse radiale varie fortement à travers l’entrefer étroit.

Appliquer les conditions aux limites des disques

La condition d’absence de glissement exige que la vitesse radiale soit nulle aux surfaces des deux disques :

[ u_r=0 \quad \text{at} \quad z=-H \text{ and } z=H. ]

L’intégration deux fois de l’équation simplifiée de quantité de mouvement à travers l’entrefer donne un profil de vitesse radiale parabolique. Le débit volumique associé à ce profil doit correspondre au volume de fluide déplacé par les disques qui se rapprochent.

Utiliser la continuité

Pour deux disques se rapprochant symétriquement, le demi-entrefer diminue à la vitesse (U) :

[ \frac{dH}{dt}=-U. ]

La variation de l’entrefer déplace le fluide à un débit proportionnel à l’aire locale du disque. L’équation de continuité axisymétrique relie alors l’écoulement radial au taux de réduction de l’entrefer.

La combinaison de l’équation de continuité et du profil de vitesse parabolique donne le gradient de pression radial

[ \frac{\partial P}{\partial r}

-\frac{3\mu U r}{2H^3}. ]

Le signe négatif indique que la pression diminue du centre vers le bord extérieur.

Intégrer le champ de pression

En prenant la pression au bord du disque comme pression de référence, (P(R)=0), l’intégration de (r) à (R) donne

[ P(r)=\frac{3\mu U}{4H^3}\left(R^2-r^2\right). ]

La pression est donc maximale au centre et diminue quadratiquement vers le bord.

Intégrer la pression pour obtenir la force

La force normale totale est la pression intégrée sur la surface circulaire du disque :

[ F=2\pi\int_0^R P(r),r,dr. ]

En substituant la distribution de pression, on obtient

[ F

2\pi\int_0^R \frac{3\mu U}{4H^3} \left(R^2-r^2\right)r,dr, ]

et donc

[ \boxed{ F=\frac{3\pi\mu U R^4}{8H^3} }. ]

Il s’agit de l’équation de Stefan selon la convention du demi-entrefer et du rapprochement symétrique indiquée.

Pourquoi la mise à l’échelle est importante

La force augmente rapidement lorsque l’entrefer se ferme

L’équation contient le facteur

[ F\propto \frac{1}{H^3}. ]

Réduire de moitié le demi-entrefer multiplie par huit la force prédite si (U), (R) et (\mu) restent constants.

Cette forte dépendance explique pourquoi de petites erreurs de mesure de l’entrefer peuvent entraîner de grandes erreurs dans la viscosité calculée ou la force prédite.

Le rayon du disque a un effet encore plus marqué

La force varie selon la quatrième puissance du rayon :

[ F\propto R^4. ]

Un disque plus grand produit une résistance nettement supérieure, car il déplace davantage de fluide et augmente la distance radiale sur laquelle la pression agit.

La viscosité et la vitesse de rapprochement ont des effets linéaires

La force est directement proportionnelle à la viscosité et à la vitesse de rapprochement :

[ F\propto \mu U. ]

Un liquide plus visqueux ou une compression plus rapide exige une force normale plus importante, à condition que le fluide reste newtonien et que les hypothèses d’écoulement rampant demeurent valides.

Application de l’équation dans une unité pilote

Préparer et inspecter l’appareil

Installez les disques circulaires, vérifiez qu’ils sont parallèles et mesurez le rayon du disque. Remplissez soigneusement l’entrefer afin d’éviter la présence d’air piégé, car la compressibilité de l’air peut fausser les mesures de force.

La température doit être surveillée, car la viscosité peut varier considérablement avec la température.

Contrôler le mouvement

Définissez un profil connu de vitesse de rapprochement (U(t)) ou imposez un historique de déplacement contrôlé (H(t)). Le contrôleur doit enregistrer la position réelle du disque plutôt que de se fier uniquement à la position commandée.

Le demi-entrefer est calculé à partir de l’écartement mesuré entre les disques :

[ H(t)=\frac{\text{écartement total des disques}}{2}. ]

Enregistrer la force et le déplacement

Pendant l’étape de compression, enregistrez :

La force normale (F(t))
La position du disque ou l’entrefer (H(t))
La vitesse de rapprochement (U(t))
La température du fluide
Le temps

Une courbe force-déplacement est particulièrement utile, car elle montre comment la résistance augmente lorsque l’entrefer se rétrécit.

Estimer la viscosité

La réorganisation de l’équation de Stefan donne

[ \boxed{ \mu

\frac{8F H^3}{3\pi U R^4} }. ]

Les étudiants peuvent calculer une viscosité apparente instantanée à partir de chaque point de mesure ou ajuster le modèle à l’ensemble des données force-temps ou force-entrefer.

Pour une viscosité constante, les valeurs calculées devraient rester approximativement cohérentes dans la plage de fonctionnement valide.

Valider la mise à l’échelle prédite

L’expérience peut tester le modèle en faisant varier un paramètre à la fois :

Tracer (F) en fonction de (U) pour tester la dépendance linéaire à la vitesse.
Tracer (F) en fonction de (\mu) en utilisant des fluides de viscosités différentes.
Tracer (F) en fonction de (H^{-3}) pour tester la dépendance à l’entrefer.
Comparer des disques de rayons différents pour tester la dépendance en (R^4).

Un graphique log-log peut faciliter l’identification expérimentale des exposants de loi de puissance prédits.

Comprendre les compromis

L’approximation d’entrefer étroit peut devenir inapplicable

Le modèle exige

[ \frac{H_0}{R}\ll 1. ]

Si l’entrefer n’est pas suffisamment faible par rapport au rayon du disque, les effets radiaux et de bord deviennent importants, et l’approximation unidimensionnelle de la lubrification perd en précision.

L’inertie peut perturber les essais à grande vitesse

La condition d’écoulement rampant exige

[ Re=\frac{H_0U_0}{\nu}\ll 1. ]

À des vitesses de rapprochement élevées ou avec des fluides de faible viscosité, les forces d’inertie peuvent devenir mesurables. La force de compression appliquée doit satisfaire la condition pratique indiquée

[ F\ll \frac{\mu^2}{\rho} \left(\frac{R}{H_0}\right)^4, ]

ce qui contribue à maintenir faibles les interférences dues à l’inertie dans l’expérience en unité pilote.

Les fuites au bord réduisent la force mesurée

La dérivation idéale suppose une limite circulaire bien définie et une pression égale à la pression de référence au bord du disque. Si le fluide s’échappe de manière irrégulière, le champ de pression et la force mesurée seront inférieurs aux valeurs prédites.

Un mauvais parallélisme peut produire le même type d’erreur en créant un entrefer non uniforme.

Le comportement newtonien peut ne pas s’appliquer

L’équation de Stefan utilise une viscosité dynamique constante. Les polymères fondus, les suspensions concentrées, les gels et de nombreux fluides structurés peuvent être rhéofluidifiants, rhéoépaississants, viscoélastiques ou présenter une contrainte seuil.

Pour ces matériaux, la force mesurée peut représenter une viscosité apparente ou dépendant du taux de cisaillement plutôt qu’une constante unique du matériau.

La mesure de l’entrefer est particulièrement sensible

Comme la force varie selon (H^{-3}), une petite erreur de mise à zéro ou une incertitude sur la séparation des disques peut entraîner une importante erreur de viscosité.

L’appareil doit être étalonné pour le déplacement, la force, le parallélisme des disques et l’entrefer initial avant la collecte des données.

Les effets transitoires peuvent être importants

L’équation décrit la réponse visqueuse à la compression dans le cadre d’hypothèses quasi stationnaires. Au début et à la fin d’un cycle de mouvement, l’accélération de l’actionneur, le temps de réponse de l’instrument, l’air piégé et l’inertie du fluide peuvent créer des écarts transitoires.

Les données doivent être interprétées dans une région où le mouvement commandé et le mouvement mesuré concordent.

Comment appliquer cela à votre projet

L’expérience pédagogique la plus fiable combine un mouvement contrôlé, une mesure précise de l’entrefer, un contrôle de la température et des vérifications dimensionnelles avant l’interprétation des données de force.

Si votre objectif principal est de valider la théorie de la lubrification : Utilisez un fluide newtonien de faible viscosité ou bien caractérisé, un faible rapport entrefer-rayon et plusieurs vitesses de rapprochement afin de tester (F\propto U) et (F\propto H^{-3}).
Si votre objectif principal est de mesurer la viscosité : Réorganisez l’équation pour calculer (\mu=8FH^3/(3\pi U R^4)), puis comparez le résultat avec une mesure réalisée à l’aide d’un viscosimètre indépendant.
Si votre objectif principal est d’étudier les effets d’échelle : Faites varier systématiquement le rayon du disque ou l’entrefer initial et comparez la force mesurée aux dépendances prédites (R^4) et (H^{-3}).
Si votre objectif principal concerne la transformation des polymères ou les fluides non newtoniens : Considérez le résultat comme une mesure de viscosité apparente et examinez la variation de la viscosité déduite avec la vitesse de rapprochement et l’entrefer.
Si votre objectif principal est d’obtenir des données fiables et défendables : Confirmez (Re\ll1), (H_0/R\ll1), le parallélisme des disques, le confinement du fluide et la stabilité de la température avant d’ajuster le modèle de Stefan.

Lorsque la géométrie et les conditions de fonctionnement respectent ses hypothèses, l’équation de Stefan établit un lien direct entre la force de compression mesurée et la physique des écoulements en entrefer étroit dominée par la viscosité.

Tableau récapitulatif :

Thème	Informations clés
Équation	F = 3πμUR⁴/(8H³)
Variables principales	F : force ; μ : viscosité ; U : vitesse de rapprochement ; R : rayon du disque ; H : demi-entrefer
Mise à l’échelle fondamentale	La force augmente avec μ, U et R⁴, et croît fortement selon H⁻³
Hypothèses principales	Écoulement rampant, entrefer étroit, fluide newtonien, absence de glissement, disques parallèles, fuites négligeables
Mesures expérimentales	Force normale, écartement des disques, vitesse de rapprochement, température et temps
Calcul de la viscosité	μ = 8FH³/(3πUR⁴)
Tests de validation	Faire varier la vitesse, la viscosité, l’entrefer ou le rayon du disque et comparer les tendances de loi de puissance mesurées
Sources d’erreur courantes	Air piégé, mauvais parallélisme, erreurs de mesure de l’entrefer, fuites au bord, inertie et variations de température
Applications pédagogiques	Théorie de la lubrification, rhéologie, transformation des polymères, revêtement, formage et formation en mécanique des fluides

Appliquez la théorie de l’écoulement par compression avec LABPARK

Faites passer l’équation de Stefan du manuel à une expérience contrôlée et mesurable. LABPARK fournit des unités pilotes pédagogiques et professionnelles pour les opérations unitaires destinées à la formation en génie chimique, bioprocédés et biotechnologies, ainsi qu’en traitement de l’environnement et de l’eau.

Nos unités pilotes aident les universités, instituts de recherche et entreprises à enseigner et à évaluer la mécanique des fluides grâce au contrôle du mouvement, à la mesure de la force et du déplacement, au suivi de la température et à des conditions de fonctionnement reproductibles. Les étudiants et les chercheurs peuvent valider la théorie de la lubrification, estimer la viscosité, étudier les comportements newtonien et non newtonien, et comprendre comment l’entrefer, le rayon du disque et la vitesse de rapprochement influencent les performances du procédé.

Contactez LABPARK pour discuter d’une solution d’unité pilote adaptée à vos objectifs de formation, de recherche ou de formation industrielle.