Calcul de perte de charge fluide non newtonien
Estimez la perte de charge dans une conduite circulaire pour un fluide non newtonien de type loi de puissance. L’outil calcule le nombre de Reynolds généralisé, le facteur de friction, la chute de pression, la perte de charge en mètres de fluide et des grandeurs utiles au dimensionnement des pompes et réseaux process.
Résultats du calcul
Guide expert du calcul de perte de charge pour un fluide non newtonien
Le calcul de perte de charge d’un fluide non newtonien est une étape décisive dans le dimensionnement des réseaux de pompage en industrie agroalimentaire, cosmétique, pharmaceutique, minière et pétrolière. Contrairement à un fluide newtonien, dont la viscosité reste constante lorsque le cisaillement varie, un fluide non newtonien change de comportement dès que le débit, le diamètre, la température ou la formulation évoluent. Dans la pratique, cela signifie qu’une conduite qui paraît correctement dimensionnée avec une hypothèse simplifiée peut devenir sous-dimensionnée, énergivore ou incapable de garantir la qualité produit si l’on ignore la rhéologie réelle.
Un exemple classique est celui des sauces, des crèmes, des boues de forage, des peintures, des pâtes concentrées et des suspensions polymériques. Ces fluides ont souvent un comportement pseudoplastique: plus le cisaillement augmente, plus la viscosité apparente diminue. Le phénomène est favorable au pompage à haut débit, mais il complique le calcul parce que le facteur de friction ne dépend plus seulement du nombre de Reynolds classique. Il faut alors utiliser un cadre plus robuste, fondé sur la loi de puissance et le nombre de Reynolds généralisé de Metzner-Reed.
Pourquoi les fluides non newtoniens exigent un calcul spécifique
Dans un tuyau, la perte de charge représente l’énergie dissipée par frottement. Pour un fluide newtonien, cette dissipation se calcule facilement avec Darcy-Weisbach, la viscosité dynamique μ étant unique pour tout le champ de cisaillement. Pour un fluide non newtonien, la viscosité apparente dépend du taux de cisaillement. Or le taux de cisaillement lui-même dépend de la vitesse moyenne, du diamètre et du profil d’écoulement. Le calcul devient donc couplé: la vitesse modifie la viscosité apparente, la viscosité apparente modifie le nombre de Reynolds généralisé, et ce nombre de Reynolds influence le facteur de friction.
En ingénierie, les conséquences sont concrètes:
- mauvaise sélection de pompe si la chute de pression est sous-estimée;
- surconsommation électrique si le diamètre est trop petit;
- dégradation produit si la vitesse ou le cisaillement sont excessifs;
- instabilité de procédé lorsque la température fait évoluer K et n;
- erreurs de mise à l’échelle entre pilote et production.
Les principaux modèles rhéologiques utilisés
Le modèle le plus simple et le plus répandu pour les calculs rapides est la loi de puissance, aussi appelée modèle d’Ostwald-de-Waele. Elle relie la contrainte de cisaillement τ et le taux de cisaillement γ̇ par l’équation τ = Kγ̇n. Le coefficient K s’exprime en Pa.sn et l’indice n est sans dimension.
- n < 1 : fluide pseudoplastique ou amincissant au cisaillement, comme beaucoup de sauces, peintures ou solutions polymériques.
- n = 1 : comportement newtonien, comme l’eau ou certains solvants.
- n > 1 : fluide dilatant ou épaississant au cisaillement, comme certaines suspensions concentrées.
La loi de puissance fonctionne bien dans une plage de cisaillement définie. En revanche, pour les fluides à seuil d’écoulement, comme des boues, pâtes ou gels structurés, les modèles de Bingham ou Herschel-Bulkley sont souvent plus adaptés. Le calculateur ci-dessus cible le cas loi de puissance, qui reste le plus fréquent pour une estimation rapide ou un pré-dimensionnement.
Équations de base pour le calcul
Le point de départ est le débit volumique Q et le diamètre intérieur D. On en déduit la vitesse moyenne:
V = 4Q / (πD²)
Pour un fluide loi de puissance, on emploie ensuite le nombre de Reynolds généralisé:
Reg = ρ Dn V2-n / [K 8n-1 ((3n+1)/(4n))n]
Ce nombre de Reynolds étend la logique de Reynolds au cas non newtonien. En régime laminaire, on peut utiliser une forme analogue à la loi classique:
f = 64 / Reg
avec f le facteur de Darcy. En régime turbulent, les corrélations deviennent plus complexes. Une référence courante est la relation de Dodge-Metzner pour les conduites lisses. Le calculateur l’emploie en mode automatique et bascule sur une approximation turbulente lisse si nécessaire pour conserver une bonne stabilité numérique.
La chute de pression linéaire vaut ensuite:
ΔPlin = f (L/D) (ρV²/2)
Les pertes singulières dues aux coudes, vannes et changements de direction s’ajoutent:
ΔPsing = KΣ (ρV²/2)
Finalement, la perte totale est:
ΔPtot = ΔPlin + ΔPsing
Valeurs typiques observées pour plusieurs fluides non newtoniens
Le tableau suivant donne des ordres de grandeur couramment rencontrés pour des fluides industriels ou alimentaires dans une plage de température proche de l’ambiante. Ces valeurs sont indicatives. Elles varient selon la concentration, la formulation, la teneur en solides et le niveau de cisaillement de référence.
| Fluide | Type de comportement | K typique (Pa.sⁿ) | n typique | Masse volumique typique (kg/m³) |
|---|---|---|---|---|
| Ketchup industriel | Pseudoplastique | 8 à 25 | 0,20 à 0,45 | 1100 à 1300 |
| Yaourt brassé | Pseudoplastique | 0,6 à 3,0 | 0,35 à 0,65 | 1020 à 1080 |
| Peinture latex | Pseudoplastique | 0,3 à 2,5 | 0,50 à 0,90 | 1100 à 1400 |
| Boue de forage à base d’eau | Pseudoplastique à seuil possible | 0,15 à 1,5 | 0,45 à 0,90 | 1050 à 1600 |
| Solution CMC concentrée | Pseudoplastique | 0,2 à 5,0 | 0,30 à 0,80 | 1000 à 1150 |
Exemple de calcul appliqué
Supposons une suspension pseudoplastique avec les paramètres suivants: ρ = 1050 kg/m³, K = 0,8 Pa.sn, n = 0,6, longueur de conduite de 50 m, diamètre intérieur de 50 mm et débit de 8 m³/h. Avec ces données, la vitesse moyenne dépasse 1 m/s, le taux de cisaillement au mur augmente, la viscosité apparente chute, et l’écoulement peut se rapprocher d’une zone de transition selon la formulation exacte. Le résultat utile n’est pas seulement la perte de charge totale, mais aussi le régime d’écoulement, la contrainte pariétale et la viscosité apparente au cisaillement de service.
Cet exemple montre pourquoi deux fluides de densité voisine peuvent nécessiter des pompes très différentes. Un fluide à n faible profite d’un amincissement rapide avec la vitesse. À débit égal, il peut parfois présenter une perte de charge plus basse qu’un fluide presque newtonien, malgré une viscosité apparente élevée à faible cisaillement. À l’inverse, les fluides dilatants deviennent plus pénalisants lorsque le débit augmente.
Comparaison d’impact du débit sur la perte de charge
Le tableau ci-dessous illustre un cas représentatif de conduite lisse de 50 mm et 50 m, avec ρ = 1050 kg/m³, K = 0,8 Pa.sn, n = 0,6 et pertes singulières modérées. Les chiffres sont cohérents avec un modèle loi de puissance et montrent surtout la sensibilité élevée de ΔP au débit.
| Débit (m³/h) | Vitesse moyenne (m/s) | Reynolds généralisé | Facteur de Darcy | Perte totale estimée (kPa) |
|---|---|---|---|---|
| 4 | 0,57 | 115 | 0,56 | 95 à 110 |
| 6 | 0,85 | 177 | 0,36 | 150 à 175 |
| 8 | 1,13 | 235 | 0,27 | 205 à 245 |
| 10 | 1,41 | 289 | 0,22 | 255 à 310 |
| 12 | 1,70 | 341 | 0,19 | 310 à 385 |
Méthode pratique de dimensionnement
- Mesurer ou récupérer les paramètres rhéologiques K et n à la température réelle de service.
- Vérifier le diamètre intérieur effectif de la conduite et non le diamètre nominal commercial.
- Convertir le débit dans une unité cohérente, idéalement m³/s pour le calcul.
- Calculer la vitesse moyenne et le nombre de Reynolds généralisé.
- Déterminer le facteur de friction avec une corrélation adaptée au régime.
- Ajouter les pertes singulières, parfois sous-estimées sur les réseaux compacts.
- Comparer la chute de pression au point de fonctionnement de la pompe, en tenant compte de la marge d’exploitation.
- Valider le cisaillement admissible si le produit est sensible, par exemple en bioprocédés ou sur produits texturés.
Erreurs fréquentes à éviter
- Utiliser une viscosité unique mesurée à un seul cisaillement pour tout le réseau.
- Oublier les singularités, notamment vannes de régulation, échangeurs, filtres et raccords sanitaires.
- Négliger l’effet de température, souvent majeur pour les polymères et produits alimentaires.
- Employer le diamètre nominal au lieu du diamètre intérieur, ce qui fausse V et ΔP.
- Ignorer la variabilité lot à lot, importante pour certaines formulations viscoélastiques.
Comment interpréter correctement le résultat
Une perte de charge élevée n’est pas forcément un problème si la pompe est correctement sélectionnée et si le produit supporte le cisaillement. L’enjeu est d’optimiser simultanément l’énergie, la qualité produit et la robustesse du procédé. Par exemple, augmenter le diamètre de 40 à 50 mm peut faire chuter fortement la perte de charge, mais au prix d’un coût d’installation supérieur. Inversement, conserver un petit diamètre peut être acceptable pour une ligne courte, un fluide peu sensible et un temps de pompage faible.
En conception avancée, il est utile de comparer plusieurs scénarios: débit nominal, débit maximal, phase de démarrage à froid, fin de lot plus concentrée et nettoyage en place. C’est précisément dans ces écarts de fonctionnement que les hypothèses simplificatrices montrent leurs limites. Un bon calcul de perte de charge sur fluide non newtonien ne doit donc pas être vu comme un simple chiffre, mais comme un outil d’aide à la décision.
Références utiles et ressources d’autorité
Pour compléter un pré-dimensionnement, il est judicieux de consulter des sources institutionnelles et académiques sur les unités, la mécanique des fluides et les systèmes de pompage:
- NIST.gov: système international d’unités et cohérence métrologique
- MIT.edu: cours de mécanique des fluides et transferts
- Energy.gov: principes de performance des systèmes de pompage
Conclusion
Le calcul de perte de charge fluide non newtonien est indispensable dès que la viscosité varie avec le cisaillement. En utilisant la loi de puissance, le nombre de Reynolds généralisé et une corrélation de friction adaptée, on obtient une estimation utile pour la conception de tuyauteries, le choix des pompes et l’analyse énergétique. Le calculateur proposé sur cette page permet d’aller vite tout en respectant les principes essentiels de la rhéologie appliquée. Pour un projet critique, la meilleure pratique reste de compléter ce pré-dimensionnement par des mesures rhéologiques sur produit réel, aux températures et cisaillements d’exploitation.