Le tenseur du taux de déformation convertit les gradients de vitesse en étirement, compression et cisaillement locaux du fluide. Dans une unité pilote de génie chimique, les ingénieurs obtiennent le champ de vitesse à partir de mesures ou de la CFD, calculent ses gradients spatiaux et forment le tenseur symétrique

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

où (\mathbf{v}=(u,v,w)). Ses termes extradiagonaux quantifient le comportement du taux de cisaillement, tandis que ses termes diagonaux décrivent l'extension ou la compression ; pour les fluides newtoniens, ces grandeurs déterminent les contraintes visqueuses, la dissipation d'énergie et d'importants paramètres de changement d'échelle.

À retenir : Le tenseur distingue la déformation réelle d'un élément fluide de la rotation d'un corps rigide. Dans les unités pilotes, cela permet de relier les gradients de vitesse mesurés ou simulés à la contrainte de cisaillement locale, à la perte de pression, aux performances de mélange, aux transferts thermique et massique, ainsi qu'aux dommages potentiels causés aux matériaux sensibles au cisaillement.

Ce que représente le tenseur du taux de déformation

Il décrit le mouvement relatif entre des éléments fluides voisins

Un élément fluide peut se translater, tourner, s'étirer, se comprimer ou subir un cisaillement lorsqu'il se déplace dans les conduites et les équipements.

Le tenseur du gradient de vitesse contient toutes les dérivées spatiales locales :

[ \nabla \mathbf{v}= \begin{bmatrix} \partial u/\partial x & \partial u/\partial y & \partial u/\partial z\ \partial v/\partial x & \partial v/\partial y & \partial v/\partial z\ \partial w/\partial x & \partial w/\partial y & \partial w/\partial z \end{bmatrix} ]

Le tenseur du taux de déformation en est la partie symétrique :

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

Ce tenseur rend compte des variations des distances entre les points matériels voisins.

Il doit être distingué de la rotation du fluide

La partie antisymétrique du gradient de vitesse représente la rotation locale d'un corps rigide :

[ \mathbf{W}=\frac{1}{2}\left(\nabla \mathbf{v}-(\nabla \mathbf{v})^T\right) ]

La rotation modifie l'orientation d'un élément fluide, mais ne modifie pas en elle-même sa forme et ne génère pas de contrainte visqueuse newtonienne.

Cette distinction est particulièrement importante à proximité des agitateurs, des dispositifs rotor-stator, des coudes et des écoulements tourbillonnaires, où des gradients de vitesse élevés pourraient sinon être entièrement assimilés à un cisaillement élevé.

Comment calculer le taux de déformation de cisaillement

Les termes extradiagonaux quantifient la déformation locale par cisaillement

Pour un écoulement bidimensionnel dans le plan (x)-(y), la composante tensorielle pertinente est

[ D_{xy}=\frac{1}{2}\left(\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right) ]

Il s'agit du taux de déformation de cisaillement tensoriel.

Le taux de cisaillement d'ingénierie correspondant s'écrit généralement

[ \dot{\gamma}{xy}=2D{xy} =\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x} ]

Le facteur deux est important. Certains rapports d'unités pilotes appellent (D_{xy}) le taux de cisaillement, tandis que d'autres utilisent (\dot{\gamma}_{xy}) ; la définition doit toujours être précisée.

Un exemple simple d'écoulement parallèle

Pour un écoulement pleinement développé entre des surfaces fixes et mobiles, la vitesse dominante est souvent (u(y)), avec un (v) négligeable. L'expression se réduit à

[ D_{xy}\approx \frac{1}{2}\frac{\partial u}{\partial y} ]

[ \dot{\gamma}_{xy}\approx \frac{\partial u}{\partial y} ]

Ainsi, le taux de cisaillement local est déterminé par la rapidité avec laquelle la vitesse varie à travers l'espace entre les surfaces ou la couche limite.

Utilisation du tenseur dans l'analyse d'une unité pilote

1. Obtenir le champ de vitesse

Les ingénieurs mesurent ou calculent d'abord les composantes de la vitesse dans l'ensemble de l'équipement.

Les sources courantes comprennent :

Les simulations CFD de conduites, vannes, mélangeurs et réacteurs
La vélocimétrie par image de particules ou les mesures par Doppler laser
Les mesures de débit et de pression combinées à des modèles d'écoulement
Les données relatives à la vitesse de l'agitateur, à la géométrie et au jeu
Les mesures du profil de vitesse à travers des canaux ou des cuves

Le champ de vitesse doit être résolu avec une précision suffisante, car la dérivation amplifie le bruit de mesure et les erreurs numériques.

2. Calculer les gradients spatiaux de vitesse

L'étape suivante consiste à évaluer des dérivées telles que

[ \frac{\partial u}{\partial y},\qquad \frac{\partial v}{\partial x},\qquad \frac{\partial u}{\partial x} ]

Ces gradients identifient les zones où l'écoulement varie rapidement dans l'espace.

Les gradients élevés apparaissent généralement près des parois solides, dans les jeux étroits, au niveau des pales d'agitateur, des restrictions de vannes, des élargissements brusques et des zones de recirculation.

3. Former le tenseur de déformation

Les gradients sont assemblés dans (\mathbf{D}). Pour un écoulement tridimensionnel,

[ \mathbf{D}= \begin{bmatrix} \partial u/\partial x & \frac{1}{2}(\partial u/\partial y+\partial v/\partial x) & \frac{1}{2}(\partial u/\partial z+\partial w/\partial x)\ \frac{1}{2}(\partial v/\partial x+\partial u/\partial y) & \partial v/\partial y & \frac{1}{2}(\partial v/\partial z+\partial w/\partial y)\ \frac{1}{2}(\partial w/\partial x+\partial u/\partial z) & \frac{1}{2}(\partial w/\partial y+\partial v/\partial z) & \partial w/\partial z \end{bmatrix} ]

Les composantes extradiagonales décrivent la déformation par cisaillement. Les composantes diagonales décrivent l'extension ou la compression normale.

4. Convertir la déformation en contrainte visqueuse

Pour un fluide newtonien isotrope, la contrainte visqueuse est liée au tenseur de déformation par

[ \boldsymbol{\tau} =2\mu\mathbf{D} +\lambda(\nabla\cdot\mathbf{v})\mathbf{I} ]

où (\mu) est la viscosité dynamique et (\lambda) un paramètre de viscosité volumique.

Pour un liquide incompressible,

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

la contrainte visqueuse déviatorique se simplifie donc en

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

Par exemple,

[ \tau_{xy}=2\mu D_{xy} =\mu\left(\frac{\partial u}{\partial y} +\frac{\partial v}{\partial x}\right) ]

Cette relation permet aux ingénieurs d'estimer la contrainte de cisaillement à la paroi, la contrainte de cisaillement interne et les sollicitations mécaniques subies par les particules en suspension, les gouttelettes ou les cellules biologiques.

Ce que révèlent les termes diagonaux

Opérations avec des liquides incompressibles

Les termes diagonaux sont

[ D_{xx}=\frac{\partial u}{\partial x},\qquad D_{yy}=\frac{\partial v}{\partial y},\qquad D_{zz}=\frac{\partial w}{\partial z} ]

Leur somme correspond au taux d'expansion volumique :

[ \nabla\cdot\mathbf{v}

\frac{\partial u}{\partial x} +\frac{\partial v}{\partial y} +\frac{\partial w}{\partial z} ]

Pour un liquide incompressible, cette somme est nulle. Les termes diagonaux individuels peuvent néanmoins être non nuls, mais l'extension locale dans une direction doit être compensée par une contraction dans une autre.

Opérations avec des gaz compressibles

Dans la compression des gaz, le transport pneumatique et les autres unités pilotes fonctionnant avec des écoulements compressibles, la divergence n'est généralement pas nulle.

La relation

[ \nabla\cdot\mathbf{v}

-\frac{1}{\rho}\frac{D\rho}{Dt} ]

relie la déformation volumique aux variations de la masse volumique du fluide. Elle est nécessaire pour modéliser l'accélération, la compression et la détente des gaz, ainsi que les variations de masse volumique induites par la pression.

Applications dans les principales unités pilotes

Conduites, raccords et vannes

Dans un écoulement en conduite, le tenseur identifie les distributions du taux de cisaillement à travers le rayon de la conduite et près de la paroi.

Ces résultats permettent de calculer :

La contrainte de cisaillement à la paroi
Le comportement de la couche limite
La perte de pression par frottement
Les besoins de pompage
Les zones vulnérables à l'érosion ou aux dépôts

Le tenseur fournit la mécanique locale, tandis que la perte de charge globale nécessite toujours une intégration avec le bilan de quantité de mouvement, les propriétés du fluide, la géométrie et le régime d'écoulement.

Cuves de mélange et réacteurs agités

À proximité des agitateurs et des chicanes, les gradients de vitesse peuvent être très hétérogènes. Les ingénieurs utilisent (\mathbf{D}) pour cartographier les zones de fort cisaillement et les comparer aux régions mal mélangées ou faiblement déformées.

L'analyse permet d'évaluer :

Les performances de l'agitateur et des chicanes
L'uniformité du mélange
La qualité de la mise en suspension
La dispersion des gaz
Le comportement du temps de mélange
L'exposition au cisaillement des matériaux sensibles

Pour un fluide newtonien incompressible, le taux local de dissipation visqueuse est généralement évalué à partir du tenseur de déformation comme suit :

[ \Phi=2\mu,\mathbf{D}:\mathbf{D} ]

Cela indique où l'énergie mécanique est convertie en chaleur et où une déformation intense peut influencer le mélange ou la qualité du produit.

Broyeurs rotor-stator et dispositifs à faible jeu

Une première estimation utile pour un dispositif rotor-stator peut être obtenue à partir de la vitesse périphérique du rotor et de la taille du jeu.

Si (N) est la vitesse de rotation en tours par unité de temps, la vitesse périphérique est

[ v_{\text{tip}}=\pi D N ]

et un taux de cisaillement caractéristique dans le jeu est

[ \dot{\gamma}\approx\frac{v_{\text{tip}}}{h} ]

où (D) est le diamètre du rotor et (h) le jeu entre le rotor et le stator.

Si la vitesse angulaire (\omega) est plutôt donnée en radians par seconde, alors

[ v_{\text{tip}}=\omega\frac{D}{2} ]

Cette estimation du jeu est utile pour un changement d'échelle préliminaire, mais une CFD détaillée ou des mesures locales sont nécessaires lorsque l'écoulement dans le jeu est fortement tridimensionnel, turbulent ou influencé par une géométrie complexe des pales.

Échangeurs de chaleur et équipements de transfert

Les gradients de vitesse influencent le cisaillement à la paroi, l'épaisseur de la couche limite et les niveaux de turbulence dans les canaux des échangeurs de chaleur.

Le champ de déformation contribue donc à expliquer les variations de :

Transfert thermique convectif
Encrassement et dépôts
Perte de pression
Distribution des temps de séjour
Contraintes mécaniques sur les phases dispersées

Il doit être utilisé conjointement avec les bilans énergétiques et les corrélations de transport appropriées, et non en remplacement de ceux-ci.

Comprendre les compromis

Un cisaillement élevé peut améliorer le procédé mais endommager les produits

L'augmentation du cisaillement peut améliorer la dispersion, l'émulsification, la réduction de la taille des particules et le transfert de matière.

Cependant, une déformation locale excessive peut endommager les cellules biologiques, déstabiliser les gouttelettes fragiles, dégrader les polymères ou provoquer un échauffement indésirable. L'objectif approprié est généralement un cisaillement contrôlé, et non un cisaillement maximal.

Un seul taux de cisaillement moyen peut être trompeur

Les équipements pilotes présentent souvent d'importantes variations spatiales de (\mathbf{D}). Une valeur moyenne dans la cuve peut masquer de courtes régions soumises à un cisaillement très élevé près des agitateurs, des vannes ou des jeux.

Les décisions de conception doivent donc tenir compte des maxima locaux, des distributions pondérées par le volume, du temps de séjour dans les zones de fort cisaillement et du temps de réponse pertinent du matériau.

La relation newtonienne entre contrainte et déformation a ses limites

La relation (\boldsymbol{\tau}=2\mu\mathbf{D}) suppose un comportement newtonien avec une viscosité connue.

Pour les fluides non newtoniens, la viscosité peut dépendre du taux de cisaillement, du temps, de la température ou de l'historique de déformation. Dans ces systèmes, le tenseur de déformation décrit toujours la cinématique, mais le modèle constitutif de contrainte doit être choisi de manière appropriée.

Les erreurs numériques et de mesure affectent les gradients

Les mesures de vitesse sont souvent moins bruitées que leurs dérivées spatiales. Une résolution spatiale insuffisante, des maillages CFD trop grossiers ou des conditions aux limites incohérentes peuvent produire des taux de déformation peu fiables.

Les calculs de gradients doivent être vérifiés au moyen des lois de conservation, d'un raffinement du maillage, d'estimations d'incertitude et de limites analytiques connues, telles que l'écoulement laminaire pleinement développé dans une conduite.

Faire le bon choix selon votre objectif

Le tenseur doit être considéré comme un lien entre les données du champ de vitesse et les décisions d'ingénierie à l'échelle de l'équipement.

Si votre objectif principal est la contrainte de cisaillement : Calculez le tenseur symétrique du gradient de vitesse, puis utilisez la relation constitutive appropriée pour obtenir les contraintes de cisaillement locales et pariétales.
Si votre objectif principal est le mélange : Cartographiez la distribution spatiale de la déformation et de la dissipation visqueuse plutôt que de vous appuyer uniquement sur la vitesse de l'agitateur ou sur le cisaillement moyen dans la cuve.
Si votre objectif principal est la perte de pression : Utilisez les informations sur le cisaillement pariétal dérivées de la déformation avec les bilans de quantité de mouvement, les corrélations de frottement, la géométrie et l'analyse du régime d'écoulement.
Si votre objectif principal est le changement d'échelle : Faites correspondre les caractéristiques pertinentes du taux de cisaillement local ou distribué, les temps de séjour et le comportement de dissipation d'énergie, et pas seulement la vitesse de rotation ou la puissance totale.
Si votre objectif principal concerne un procédé biologique ou sensible au cisaillement : Identifiez les zones de forte déformation et limitez l'exposition selon la sensibilité des cellules, des gouttelettes, des particules ou du produit.

Lorsqu'il est appliqué avec le modèle constitutif approprié et une résolution spatiale suffisante, le tenseur du taux de déformation transforme les gradients de vitesse d'une unité pilote en informations exploitables sur les performances de l'écoulement, la qualité du produit et les risques liés au changement d'échelle.

Tableau récapitulatif :

Domaine d'analyse	Ce que révèle le tenseur	Valeur pour l'unité pilote
Déformation par cisaillement	Composantes extradiagonales et taux de cisaillement local	Prévoit la contrainte de cisaillement et l'exposition des matériaux
Extension et compression	Composantes diagonales et divergence de l'écoulement	Évalue la déformation des liquides et l'expansion des gaz
Contrainte visqueuse	Contrainte résultant de la déformation et de la viscosité du fluide	Facilite l'analyse des contraintes pariétales et des équipements
Mélange et réacteurs	Zones de fort cisaillement et dissipation d'énergie	Améliore l'homogénéisation, la dispersion et le changement d'échelle
Conduites et vannes	Gradients pariétaux et points chauds de déformation	Oriente les études de perte de pression, d'érosion et de dépôts
Équipements de transfert thermique	Comportement de la couche limite et du cisaillement	Aide à évaluer le transfert thermique et l'encrassement
Procédés sensibles au cisaillement	Déformation locale et temps de séjour	Réduit les dommages causés aux cellules, aux gouttelettes et aux polymères

Appliquez l'analyse de la dynamique des fluides à votre unité pilote

Besoin d'informations fiables sur le cisaillement, le mélange, le transfert de matière ou les performances liées au changement d'échelle ? Contactez LABPARK pour discuter d'une unité pilote d'opérations unitaires à vocation éducative ou professionnelle, adaptée à votre application. LABPARK accompagne les universités, les instituts de recherche et les entreprises avec des unités pilotes destinées au génie chimique, aux bioprocédés et aux biotechnologies, ainsi qu'au traitement de l'environnement et de l'eau. Nos systèmes vous aident à évaluer les gradients de vitesse, les taux de déformation, les contraintes visqueuses, la dissipation d'énergie et les performances des équipements dans des conditions opératoires réalistes, afin de fournir à votre équipe des données pratiques pour l'enseignement, la recherche, le développement des procédés, la protection des produits et un changement d'échelle maîtrisé.