- Fluide de Bingham
-
Le fluide de Bingham est un modèle théorique de milieu viscoplastique qui correspond à un comportement de solide parfait sous faibles contraintes, et à un comportement de fluide visqueux au delà d'une contrainte-seuil. Ce modèle porte le nom d’Eugène Bingham qui en a donné l'expression mathématique[1].
Ce modèle mathématique est d'emploi courant en génie pétrolier pour organiser le pompage des boues ou l'évacuation des rebuts[2]. Un autre bon exemple est le dentifrice[3], qu'on ne peut extruder que moyennant l'application d’une certaine pression sur le tube, et qui se dépose comme un gel quasi-solide.
Sommaire
Le modèle théorique de Bingham
Bingham, pour présenter ses observations relatives au comportement des peintures[4], a cherché à tracer le débit en fonction de la contrainte appliquée, comme sur la figure 1 ci-contre : il s'agit d'une comparaison du comportement d'un fluide de Bingham avec un fluide newtonien (en rouge). Le diagramme s’interprète de la façon suivante : avec un liquide newtonien, l’augmentation de la pression d'entrée d'une conduite crée un champ de contrainte dans le fluide qui le met en mouvement (cisaillement), et le débit augmente en proportion de cet accroissement. Mais si la conduite est remplie d'un fluide de Bingham (courbe caractéristique en bleu), il faut que le cisaillement atteigne une certaine intensité, dite « contrainte d'écoulement », avant que le débit augmente avec la pression. Ces propriétés confèrent aux fluides de Bingham une texture crémeuse avec en surface des rides et des fronces, contrairement à la surface unie que présente tout fluide newtonien.
La Figure 2 correspond à la représentation donnée dans les manuels actuels[5],[3] : la contrainte de cisaillement est en ordonnée, fonction de la vitesse de déformation ; cette représentation a éclipsé la précédente, car le débit n'est pas une donnée intrinsèque au matériau : il dépend du diamètre de la conduite ; tandis que la vitesse de déformation (plutôt que le débit, qui en est une résultante) permet de construire le profil de vitesse. Elle est proportionnelle au débit, mais ne dépend pas du diamètre de la conduite. Sur ce diagramme (comme sur le précédent), un liquide newtonien se déforme pour toute valeur non nulle de cisaillement (contrainte) imposée. Le fluide de Bingham, au contraire, ne se déforme pas (il n'y a donc pas de débit) en dessous d'un certain niveau de contrainte. La courbe caractéristique de réponse du fluide newtonien est une droite, et la pente de cette droite est la viscosité, qui suffit à caractériser ce liquide. Le fluide de Bingham requiert, lui, deux coefficients : la contrainte d’écoulement et la pente de la partie rectiligne de sa caractéristique, appelée « viscosité plastique » ou « viscosité de Bingham ».
En termes mathématiques, si l'on convient de noter τ la contrainte de cisaillement, τ0 la cohésion apparente, dite « contrainte d'écoulement » et ∂u/∂y la vitesse de déformation (cf. l'article viscosité), les relations ci-dessus s'écrivent :
Mécanisme moléculaire sous-jacent
Les fluides de Bingham sont des milieux composés de micro-particules (par ex. des feuillets d’argile) ou de molécules (par ex. les macromolécules d'un polymère) exerçant entre elles une loi d'interaction particulière, qui se traduit au plan macroscopique par une cohésion apparente. Cette cohésion détermine une structure lâche, naguère appelée semi-solide : il faut une certaine dépense d'énergie pour rompre sa cohésion ; une fois que les interactions de cohésion sont rompues, les particules n'exercent plus entre elles que des forces d'interaction « visqueuses », c'est-à-dire régie par la loi de Newton ; enfin si le champ de contrainte est annulé, les intéractions de cohésion redeviennent prépondérantes et les particules s'organisent en un édifice moléculaire stable (solide).
Application aux pertes de charge linéaires
L’un des problèmes de base de l’hydraulique consiste à évaluer la chute de pression entrée/sortie dans un réseau en charge[6]. Une fois que l’on connaît la rugosité des conduites, on dispose de plusieurs moyens d’en déduire les pertes de charges linéaires[7], et on peut alors calculer, par ex. le coût de pompage ou déterminer le débit pour une certaine pression d'entrée. Ainsi on obtient la perte de charge grâce à l’Équation de Darcy-Weisbach:
où :
- est la perte de charge linéaire (en m)
- est le coefficient de perte de charge (nombre sans dimension)
- est la longueur de conduite (en m)
- est l’accélération de la pesanteur (en m/s²)
- est le diamètre de la conduite (en m)
- est la vitesse moyenne d'écoulement (en m/s)
Seulement il est particulièrement difficile d'exprimer analytiquement le coefficient de perte de charge pour des fluides non-newtoniens : aussi s'en remet-on à des approximations.
Écoulement laminaire
La première formulation exacte des perte de charge linéaires dans un fluide de Bingham dans un écoulement laminaire en charge (celui d'une conduite sous pression) a été publiée par l'hydraulicien Buckingham[8]. Le résultat, dit équation de Buckingham-Reiner, s'exprime sous forme adimensionnelle ainsi :
où :
- est le coefficient de perte de charge laminaire (nombre sans dimension)
- est le nombre de Reynolds
- est le nombre de Hedström
Le nombre de Reynolds et le nombre d'Hedstrom sont définis par :
- , et
- ,
où :
- est le poids volumique du fluide ( kg/m3)
- est la viscosité cinématique du fluide (m²/s)
Écoulement turbulent
Darby et Melson ont proposé la loi empirique suivante[9] :
où:
- est le coefficient de perte de charge turbulent (sans dimension)
Approximations de l'équation de Buckingham-Reiner
Bien que l’équation de Buckingham-Reiner soit résoluble (c'est une équation algébrique du quatrième degré en λ), la complexité de la solution justifiait, jusqu’à l’avènement de l’informatique, la recherche de formules approchées, d’application plus simple.
Équation de Swamee-Aggarwal
L’équation de Swamee-Aggarwal, qui donne directement le coefficient de perte de charge de l'équation de Darcy–Weisbachpour l'écoulement laminaire des fluides de Bingham[10], n’est qu’une approximation des racines de l’équation Buckingham-Reiner, mais l'écart aux données expérimentales reste bien en deçà des incertitudes de mesure. Elle s’écrit :
La solution de Danish-Kumar
Danish et al. ont publié une méthode explicite pour calculer le coefficient λ en utilisant la méthode de décomposition d’Adomian[11]. Le coefficient de perte de charge fait intervenir deux facteurs :
où :
- , and
Formule donnant les pertes de charges dans un écoulement quelconque
En 1981, Darby et Melson, suivant une approche proposée précédemment par Churchill[12] ont obtenu une expression synthétique donnant le coefficient de perte de charge pour tous les régimes d'écoulement[9] :
où:
On peut combiner les équations de Swamee-Aggarwal et de Darby-Melson pour former une expression synthétique donnant le coefficient de perte de charge avec le fluide de Bingham pour tous les régimes. Il faut observer que la rugosité relative n'apparaît nulle part dans ces équations : le coefficient de perte de charge linéaire des fluides de Bingham n’en dépend pas.
Notes et références
- E.C. Bingham,(1916) U.S. Bureau of Standards Bulletin, 13, 309-353 "An Investigation of the Laws of Plastic Flow"
- M. A. Mian,Mohammed et A. Mian, Petroleum engineering handbook for the practicing engineer, vol. 2, PennWell Books, 1992, 688 p. (ISBN 0878143793), p. 368. Cf. par ex.
- ISBN 0-9632036-1-4 J. F. Steffe (1996) Rheological Methods in Food Process Engineering 2nd ed
- E. C. Bingham, Fluidity and Plasticity, New York, McGraw-Hill, 1922, p. 219 Cf.
- Paul Germain, Mécanique des milieux continus, vol. 1 : Théorie générale, Masson & Cie, 1972, 418 p. (ISBN 2-225-35937-4), p. 229, fig. 13-a D’après
- Ron Darby, Chemical Engineering Fluid Mechanics., Marcel Dekker, 1996 (ISBN 0824704444).
- Armando Lencastre, Hydraulique générale, Eyrolles, 1957 (réimpr. 1983,1995), 636 p. (ISBN 2-21201894-0), « 4 - Ecoulements en régime permanent »
- E. Buckingham, « on Plastic Flow through Capillary Tubes », dans ASTM Proceedings, no 21, 1921, p. 1154–1156. Référence :
- R. Darby, J. Melson, « How to predict the friction factor for flow of Bingham plastics », dans Chemical Engineering, no 28, 1981, p. 59–61. Référence :
- Swamee, P.K. and Aggarwal, N.(2011). "Explicit equations for laminar flow of Bingham plastic fluids". Journal of Petroleum Science and Engineering. doi:10.1016/j.petrol.2011.01.015.
- M. Danish, et al., « Approximate explicit analytical expressions of friction factor for flow of Bingham fluids in smooth pipes using Adomian decomposition method », dans Communications in Nonlinear Science and Numerical Simulation, no 16, 1981, p. 239–251.
- S.W. Churchill, « Friction factor equation spans all fluid-flow regimes », dans Chemical Engineering, 1977, p. 91–92, et S.W. Churchill et R.A. Usagi, « A general expression for the correlation of rates of transfer and other phenomena », dans AIChE Journal, vol. 18, no 6, 1972, p. 1121-1128. Cf.
Annexes
Articles connexes
Bibliographie
- Paul Germain (scientifique), Mécanique des milieux continus, vol. 1 : Théorie générale, Masson & Cie, 1972, 418 p. (ISBN 2-225-35937-4), p. 229, fig. 13-a
- J. Lemaitre J.-L. Chaboche, Mécanique des matériaux solides, Paris, Dunod, 1988, 544 p. (ISBN 2-04-18618-2)
- (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Bingham plastic » (voir la liste des auteurs)
Wikimedia Foundation. 2010.