Le tenseur des contraintes visqueuses fait le lien entre ce que mesurent les instruments des unités pilotes et la manière dont les fluides transportent la quantité de mouvement. Dans un fluide newtonien, les contraintes visqueuses sont calculées à partir des gradients de vitesse locaux et de la viscosité dynamique. Dans les unités pilotes de génie chimique, cette formulation est utilisée pour interpréter le cisaillement à la paroi, les profils de vitesse, les pertes de charge, le comportement du mélange, l'écoulement dans les échangeurs de chaleur et les décisions de changement d'échelle.

Le tenseur des contraintes convertit la déformation du fluide en contrainte mécanique. Pour les liquides newtoniens incompressibles, la pression représente la partie isotrope de la contrainte totale, tandis que le tenseur visqueux décrit le cisaillement et la déformation causés par les gradients de vitesse.

Ce que représente le tenseur des contraintes visqueuses

Contrainte en un point du fluide

Un élément de fluide subit des forces de surface provenant du fluide environnant. Ces forces varient selon l'orientation de la surface ; elles sont donc représentées par un tenseur du second ordre plutôt que par un scalaire unique.

Le tenseur des contraintes de Cauchy total est généralement décomposé comme suit :

[ \boldsymbol{\sigma}=-p\mathbf{I}+\boldsymbol{\tau}, ]

où (p) est la pression thermodynamique, (\mathbf{I}) est le tenseur identité et (\boldsymbol{\tau}) est le tenseur des contraintes visqueuses.

La convention de signe ci-dessus traite la pression positive comme une compression. D'autres textes utilisent des signes différents, la convention doit donc rester cohérente tout au long d'une expérience ou d'un modèle.

Symétrie du tenseur

Pour les fluides conventionnels sans couples volumiques distribués, le tenseur des contraintes est symétrique :

[ \tau_{xy}=\tau_{yx},\qquad \tau_{xz}=\tau_{zx},\qquad \tau_{yz}=\tau_{zy}. ]

Bien qu'un tenseur (3 \times 3) contienne neuf entrées, cette symétrie laisse six composantes de contrainte indépendantes : trois contraintes normales et trois contraintes de cisaillement.

Cette propriété est importante dans les calculs des unités pilotes car la contrainte agissant sur une surface et l'équilibre des moments correspondant doivent être physiquement compatibles. Les fluides polaires spécialisés ou les fluides soumis à des couples volumiques importants nécessitent une théorie plus générale, mais cela sort du cadre de la plupart des expériences d'opérations unitaires conventionnelles.

Comment les gradients de vitesse produisent des contraintes visqueuses

Équation constitutive newtonienne

Pour un fluide newtonien isotrope, la contrainte visqueuse est linéairement liée au tenseur du taux de déformation :

[ \boldsymbol{\tau}

2\mu\mathbf{D} + \lambda(\nabla\cdot\mathbf{v})\mathbf{I}, ]

avec

[ \mathbf{D}

\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right]. ]

Ici, (\mu) est la viscosité dynamique, (\lambda) est le second coefficient de viscosité et (\mathbf{v}) est le champ de vitesse.

Le tenseur (\mathbf{D}) capture la partie du gradient de vitesse associée à la déformation du fluide. La rotation de corps rigide change l'orientation d'un élément de fluide mais ne crée pas, en soi, de contrainte visqueuse newtonienne.

Opérations sur les liquides incompressibles

Pour un écoulement incompressible,

[ \nabla\cdot\mathbf{v}=0. ]

L'équation constitutive se réduit donc à

[ \boldsymbol{\tau}=2\mu\mathbf{D}. ]

Pour un écoulement de cisaillement simple dans la direction (x) avec une vitesse variant selon (y),

[ v_x=v_x(y), ]

la contrainte de cisaillement pertinente devient

[ \tau_{xy}

\mu\frac{dv_x}{dy}. ]

C'est la relation la plus directement appliquée dans les démonstrations d'écoulement laminaire en conduite, les bancs de frottement à la paroi, les expériences de viscosimètre et les analyses de mélange simplifiées.

Composantes normales et de cisaillement

Les termes hors diagonale, tels que (\tau_{xy}), décrivent les contraintes de cisaillement. Ils apparaissent lorsque des couches de fluide adjacentes se déplacent à des vitesses différentes.

Les termes diagonaux décrivent les contraintes normales visqueuses. Dans un écoulement de cisaillement simple, stationnaire et incompressible avec une viscosité constante, ces termes s'annulent souvent ou deviennent négligeables, mais ils peuvent être importants dans les écoulements extensionnels, les passages convergents, les jets et les géométries de mélangeurs complexes.

Comment le tenseur entre dans le bilan de quantité de mouvement

Équation locale de la quantité de mouvement

L'application d'un bilan de quantité de mouvement à un élément de fluide infinitésimal donne

[ \rho\frac{D\mathbf{v}}{Dt}

\rho\mathbf{g} + \nabla\cdot\boldsymbol{\sigma}, ]

où (\rho) est la masse volumique, (\mathbf{g}) est l'accélération gravitationnelle et (D\mathbf{v}/Dt) est l'accélération matérielle.

L'utilisation de la décomposition des contraintes donne

[ \rho\frac{D\mathbf{v}}{Dt}

\rho\mathbf{g}

\nabla p + \nabla\cdot\boldsymbol{\tau}. ]

Le gradient de pression représente la force isotrope qui entraîne ou résiste au mouvement global du fluide. La divergence du tenseur visqueux représente la force visqueuse nette par unité de volume résultant des variations spatiales de la contrainte visqueuse.

Lien avec l'équation de Navier–Stokes

Pour un fluide newtonien à viscosité constante, le terme visqueux se simplifie en

[ \nabla\cdot\boldsymbol{\tau}

\mu\nabla^2\mathbf{v} ]

lorsque l'écoulement est incompressible. L'équation de la quantité de mouvement devient alors

[ \rho\frac{D\mathbf{v}}{Dt}

\rho\mathbf{g}

\nabla p + \mu\nabla^2\mathbf{v}. ]

Cette équation est la base mathématique pour interpréter les mesures de vitesse et de pression dans les équipements à l'échelle pilote.

Traction de surface et opérations tensorielles

La traction agissant sur une surface avec une normale unitaire (\mathbf{n}) est obtenue à partir du tenseur des contraintes :

[ \mathbf{t}=\boldsymbol{\sigma}\cdot\mathbf{n}. ]

Le résultat est un vecteur contenant à la fois des composantes de force normale et tangentielle sur cette surface. Cette opération est essentielle lors du calcul des forces sur les parois des conduites, les surfaces de garnissage, les pales d'agitateur, les parois des échangeurs de chaleur et d'autres limites d'équipement.

Comme la multiplication tenseur-vecteur dépend de l'ordre de l'opération, les conventions d'indices et de matrices adoptées doivent être utilisées de manière cohérente. Pour les tenseurs des contraintes symétriques des fluides newtoniens conventionnels, les relations physiques pertinentes sont cohérentes avec l'équilibre des couples.

Applications dans les opérations unitaires des unités pilotes

Expériences d'écoulement en conduite et de frottement fluide

Dans une conduite circulaire pleinement développée, le gradient de pression axial est équilibré par la variation radiale de la contrainte de cisaillement. La contrainte de cisaillement à la paroi est directement liée à la perte de charge sur une longueur de conduite mesurée :

[ \tau_w

-\frac{R}{2}\frac{dp}{dx}, ]

pour une conduite de rayon (R), sous réserve des hypothèses habituelles d'écoulement stationnaire et pleinement développé.

En régime laminaire, la combinaison de cet équilibre avec la relation constitutive newtonienne produit le profil de vitesse parabolique. Les mesures de débit et de perte de charge en unité pilote peuvent donc être comparées à la contrainte à la paroi, au facteur de frottement et à la demande de la pompe prédits.

Pompes et réseaux de tuyauterie

La formulation tensorielle identifie comment la quantité de mouvement est dissipée par les gradients de vitesse près des parois solides et dans les raccords. Ces pertes locales contribuent à la hauteur manométrique que les pompes doivent fournir.

À l'échelle pilote, les données de perte de charge résultantes aident à déterminer si une pompe, un diamètre de conduite, une disposition de vanne ou un débit de fonctionnement sont appropriés pour un procédé plus vaste. Le modèle de contrainte est particulièrement utile pour comparer des fluides ayant des viscosités différentes.

Échangeurs de chaleur

Dans les canaux des échangeurs de chaleur, les contraintes visqueuses déterminent la distribution de la vitesse et la résistance hydraulique. Le champ de vitesse influence ensuite le temps de séjour, le développement de la couche limite et le coefficient de transfert thermique convectif.

Un modèle de contrainte seul ne prédit pas le transfert thermique. Il fournit le champ de quantité de mouvement qui doit être couplé à l'équation de l'énergie et, le cas échéant, à la viscosité et à la masse volumique dépendantes de la température.

Cuves de mélange et réacteurs

Dans les cuves agitées, les agitateurs créent de fortes variations spatiales de vitesse. Le tenseur des contraintes visqueuses décrit le taux de déformation local, tandis que le champ d'écoulement résultant affecte le mélange, la circulation, la suspension et la dissipation d'énergie.

Pour les fluides très visqueux, les prédictions de contrainte laminaire ou transitionnelle peuvent être étroitement liées au couple et à la puissance du mélangeur. Pour un fonctionnement turbulent, la contrainte d'écoulement moyen doit également inclure l'effet des fluctuations de vitesse.

Colonnes à garnissage et équipement de séparation

Dans les lits et colonnes garnis, des gradients de vitesse se produisent autour des surfaces de garnissage et dans les passages d'écoulement étroits. Les contraintes associées contribuent à la perte de charge par frottement et influencent la distribution du liquide.

Le cadre tensoriel est le plus utile lorsque la mauvaise distribution locale, le passage préférentiel (channeling) ou les effets de paroi sont importants. Les corrélations basées uniquement sur le débit moyen peuvent être adéquates pour la conception de routine, mais elles ne révèlent pas les concentrations de contraintes locales.

De l'écoulement laminaire à l'écoulement turbulent

Contraintes de Reynolds

Dans un écoulement turbulent, la vitesse instantanée est décomposée en parties moyenne et fluctuante :

[ u_i=\overline{u_i}+u_i'. ]

La moyenne des équations de la quantité de mouvement introduit les contraintes de Reynolds, telles que

[ \tau_{xy}^{R}

-\rho,\overline{u'v'}. ]

Ce ne sont pas des contraintes visqueuses moléculaires. Elles représentent le transport de quantité de mouvement causé par les fluctuations de vitesse turbulentes et doivent être modélisées ou mesurées lors de la résolution de l'écoulement moyen.

Transport apparent de quantité de mouvement

Les contraintes de Reynolds augmentent le transport effectif de quantité de mouvement par rapport à la viscosité moléculaire seule. C'est pourquoi les pertes de charge turbulentes, le frottement à la paroi et les performances de mélange ne peuvent généralement pas être prédits en utilisant l'équation constitutive laminaire sans modification.

Dans les modèles d'ingénierie, des fermetures de turbulence, des modèles de viscosité turbulente ou des corrélations de frottement ajustées expérimentalement sont utilisés pour représenter ce transport supplémentaire.

Comportement près de la paroi

La sous-couche visqueuse près d'une paroi imperméable est dominée par la viscosité moléculaire, tandis que plus loin de la paroi, le transport turbulent devient de plus en plus important. Les méthodes numériques de champ peuvent résoudre cette région directement ou utiliser des fonctions de paroi pour connecter les régions d'écoulement proche de la paroi et extérieur.

Les unités pilotes peuvent utiliser les pertes de charge mesurées, les profils de vitesse ou des sondes locales pour évaluer si le modèle de turbulence sélectionné représente l'équipement avec précision.

Comprendre les compromis

Les modèles tensoriels sont locaux, les instruments sont souvent moyennés

Le tenseur des contraintes est défini point par point, mais les instruments des unités pilotes rapportent généralement des quantités moyennées telles que le débit total, la pression d'entrée et de sortie, le couple ou la température globale.

Une perte de charge mesurée peut valider une conséquence intégrée du champ de contraintes, mais elle ne détermine pas de manière unique chaque valeur locale de contrainte ou de vitesse.

Les hypothèses newtoniennes peuvent échouer

La relation

[ \boldsymbol{\tau}=2\mu\mathbf{D} ]

suppose un fluide newtonien avec une viscosité indépendante de l'historique de déformation et du taux de cisaillement. Les solutions polymères, les boues, les émulsions, les pâtes et certains fluides biologiques peuvent être rhéofluidifiants, rhéoépaississants, viscoélastiques ou présenter un seuil d'écoulement.

Pour ces matériaux, l'équation constitutive doit être remplacée par un modèle rhéologique approprié. L'utilisation d'une viscosité constante peut produire des prédictions trompeuses de la perte de charge, de la puissance du mélangeur et du comportement de changement d'échelle.

L'incompressibilité a des limites

La condition (\nabla\cdot\mathbf{v}=0) est généralement adaptée aux opérations unitaires sur les liquides lorsque les changements de masse volumique sont faibles. Elle peut ne pas être adaptée aux gaz à haute vitesse, aux changements de pression importants, aux écoulements compressibles ou aux procédés avec une dilatation thermique substantielle.

La formulation correcte doit tenir compte de la variation de la masse volumique lorsque la compressibilité affecte l'expérience.

Le changement d'échelle ne préserve pas toutes les quantités

Le simple respect des dimensions géométriques ne garantit pas des champs de contraintes dynamiquement similaires. Le nombre de Reynolds, le régime d'écoulement, la viscosité, la rugosité, les conditions de l'agitateur, le temps de séjour et les groupes adimensionnels pertinents peuvent tous influencer le résultat.

Une expérience pilote est la plus précieuse lorsque ses conditions de fonctionnement et ses résultats mesurés sont sélectionnés pour soutenir un critère de changement d'échelle clairement défini.

Comment appliquer cela à votre projet

La formulation mathématique devient utile lorsqu'elle est reliée à des variables de fonctionnement mesurables et à une question de conception spécifique.

Si votre objectif principal est la prédiction de la perte de charge : Mesurez le débit, les propriétés du fluide, la géométrie de la conduite ou du canal, et la perte de pression, puis utilisez le bilan des contraintes visqueuses pour évaluer le cisaillement à la paroi et le comportement du frottement.
Si votre objectif principal est les profils de vitesse laminaire : Appliquez la relation constitutive newtonienne au gradient de vitesse mesuré et comparez le profil prédit avec les données de vitesse locales.
Si votre objectif principal est l'écoulement turbulent en conduite ou en canal : Incluez les contraintes de Reynolds via un modèle de turbulence approprié ou une corrélation validée plutôt que de traiter la viscosité moléculaire comme le seul mécanisme de transport de quantité de mouvement.
Si votre objectif principal est la performance de mélange : Reliez les taux de déformation locaux et la dissipation d'énergie à la vitesse de l'agitateur, au couple, à la viscosité et aux mesures de concentration ou de traceur.
Si votre objectif principal est la conception d'échangeurs de chaleur : Utilisez la formulation des contraintes pour obtenir le champ de vitesse, puis couplez-le à l'équation de l'énergie pour prédire le transfert thermique convectif.
Si votre objectif principal est le changement d'échelle : Préservez les groupes adimensionnels et le régime d'écoulement qui contrôlent le comportement pertinent en matière de contrainte, de perte de charge, de mélange ou de transfert thermique.

Une unité pilote transforme le tenseur des contraintes visqueuses d'une équation abstraite de mécanique des milieux continus en un outil pratique pour mesurer, valider et mettre à l'échelle le transport des fluides.

Tableau récapitulatif :

Application	Ce que le tenseur des contraintes aide à déterminer	Mesures ou entrées clés
Écoulement en conduite	Contrainte de cisaillement à la paroi, profils de vitesse, frottement et perte de charge	Débit, perte de pression, rayon, viscosité
Pompes et tuyauterie	Pertes hydrauliques et besoins en pompage	Disposition des conduites, raccords, débit, propriétés du fluide
Échangeurs de chaleur	Distribution de la vitesse, résistance hydraulique et conditions d'écoulement pour le transfert thermique	Géométrie du canal, débit, température, viscosité
Cuves de mélange et réacteurs	Taux de déformation, couple, puissance, circulation et comportement de mélange	Vitesse de l'agitateur, couple, viscosité, données de traceur
Colonnes à garnissage	Perte de charge par frottement, distribution du liquide et passage préférentiel	Géométrie du garnissage, débit, propriétés du fluide
Systèmes turbulents	Transport de quantité de mouvement supplémentaire dû aux contraintes de Reynolds	Nombre de Reynolds, données de vitesse, modèle de turbulence
Études de changement d'échelle	Similitude dynamique et préservation du comportement d'écoulement pertinent	Groupes adimensionnels, conditions de fonctionnement, géométrie

Reliez la théorie de la mécanique des fluides aux résultats des unités pilotes

Transformez les calculs du tenseur des contraintes visqueuses en données d'ingénierie fiables avec les unités pilotes d'opérations unitaires éducatives et professionnelles LABPARK. Nos systèmes aident les universités, les instituts de recherche et les entreprises à étudier la perte de charge, le cisaillement à la paroi, le mélange, le transfert thermique, l'écoulement turbulent, la rhéologie et le changement d'échelle dans des conditions de fonctionnement mesurables.

LABPARK fournit des solutions d'unités pilotes pour le génie chimique, les bioprocédés et la biotechnologie, ainsi que les applications de traitement de l'eau et de l'environnement. Utilisez des expériences pratiques pour valider des modèles, comparer le comportement des fluides, former les étudiants et les opérateurs, et prendre de meilleures décisions en matière d'équipement et de procédés.

Contactez LABPARK pour discuter de vos objectifs expérimentaux et identifier l'unité pilote adaptée à vos besoins d'enseignement, de recherche, de développement de procédés ou de formation industrielle.