Écoulement de Stokes

Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.Si ce bandeau n'est plus pertinent, retirez-le. Cliquez ici pour en savoir plus.

Cet article ne cite pas suffisamment ses sources (mai 2023).

Si vous disposez d'ouvrages ou d'articles de référence ou si vous connaissez des sites web de qualité traitant du thème abordé ici, merci de compléter l'article en donnant les références utiles à sa vérifiabilité et en les liant à la section « Notes et références ».

En pratique : Quelles sources sont attendues ? Comment ajouter mes sources ?

Un écoulement de Stokes (ou écoulement rampant) caractérise un fluide visqueux qui s'écoule lentement en un lieu étroit ou autour d'un petit objet, dont les effets visqueux dominent alors sur les effets inertiels. On parle parfois de fluide de Stokes par opposition à fluide parfait. Il est en effet régi par une version simplifiée de l'équation de Navier-Stokes, l'équation de Stokes, dans laquelle les termes inertiels sont absents. Le nombre de Reynolds mesure le poids relatif des termes visqueux et inertiel dans l'équation de Navier-Stokes et est faible (beaucoup plus petit que 1) dans un écoulement de Stokes.

L'équation de Stokes permet en particulier de décrire la décantation de très petites particules dans les liquides voire les gaz (cas des cristaux de glace décantant dans la haute atmosphère ou du brouillard commun), ainsi que les écoulements de liquide dans les dispositifs microfluidiques. Les écoulements de Couette et de Poiseuille sont aussi décrits par cette équation.

Équation de Stokes

L'équation de Stokes, qui décrit l'écoulement d'un fluide newtonien incompressible en régime permanent et à faible nombre de Reynolds, s'écrit :

μ Δ v → = g r a d → p − ρ f → {\displaystyle \mu \Delta {\vec {v}}={\overrightarrow {\mathrm {grad} }}\,p-\rho {\vec {f}}}

où :

Conditions d'application

L'équation de Stokes est une forme simplifiée de l'équation de quantité de mouvement contenue dans les équations de Navier-Stokes. Pour un fluide newtonien, celle-ci s'écrit :

ρ ∂ v → ∂ t + ρ ( v → ⋅ g r a d → ) v → = ρ f → − g r a d → p + μ Δ v → + ( μ ′ + μ ) g r a d → ( d i v v → ) {\displaystyle \rho {\frac {\partial {\vec {v}}}{\partial t}}+\rho \left({\vec {v}}\cdot {\overrightarrow {\mathrm {grad} }}\right){\vec {v}}=\rho {\vec {f}}-{\overrightarrow {\mathrm {grad} }}\,p+\mu \Delta \,{\vec {v}}+\left(\mu '+\mu \right){\overrightarrow {\mathrm {grad} }}(\mathrm {div} \,{\vec {v}})}

où :

NB : pour établir cette formule, on doit supposer que les variations spatiales de μ {\displaystyle \mu } et μ ′ {\displaystyle \mu '} sont négligeables.

Si de plus le fluide est incompressible (bonne approximation pour les liquides), alors div v → = 0 {\displaystyle {\text{div}}\,{\vec {v}}=0} et l'équation se simplifie :

ρ ∂ v → ∂ t + ρ ( v → ⋅ g r a d → ) v → = ρ f → − g r a d → p + μ Δ v → {\displaystyle \rho {\frac {\partial {\vec {v}}}{\partial t}}+\rho \left({\vec {v}}\cdot {\overrightarrow {\mathrm {grad} }}\right){\vec {v}}=\rho {\vec {f}}-{\overrightarrow {\mathrm {grad} }}\,p+\mu \Delta \,{\vec {v}}}

Dans la plupart des cas, on écrit l'équation de Stokes pour un écoulement stationnaire (ou quasi-stationnaire), où le terme ρ ∂ v → ∂ t {\displaystyle \rho {\frac {\partial {\vec {v}}}{\partial t}}} est négligeable devant le terme visqueux μ Δ v → {\displaystyle \mu \Delta \,{\vec {v}}} . L'équation se simplifie encore : ρ ( v → ⋅ g r a d → ) v → = ρ f → − g r a d → p + μ Δ v → {\displaystyle \rho \left({\vec {v}}\cdot {\overrightarrow {\mathrm {grad} }}\right){\vec {v}}=\rho {\vec {f}}-{\overrightarrow {\mathrm {grad} }}\,p+\mu \Delta \,{\vec {v}}} .

On peut évaluer l'ordre de grandeur des termes inertiels et visqueux dans cette équation. Si la vitesse caractéristique du liquide est U {\displaystyle U} , et si l'échelle spatiale typique de variation de la vitesse est L {\displaystyle L} (celle-ci pourra être imposée par les dimensions du canal dans lequel le liquide s'écoule, les dimensions de l'objet autour duquel le liquide s'écoule, etc.), alors :

ρ ( v → ⋅ g r a d → ) v → ∼ ρ U 2 L {\displaystyle \rho \left({\vec {v}}\cdot {\overrightarrow {\mathrm {grad} }}\right){\vec {v}}\sim {\frac {\rho U^{2}}{L}}} :

et

μ Δ v → = ρ ν Δ v → ∼ ρ ν U L 2 {\displaystyle \mu \Delta \,{\vec {v}}=\rho \nu \Delta \,{\vec {v}}\sim {\frac {\rho \nu U}{L^{2}}}} :

où ν {\displaystyle \nu } est la viscosité cinématique du liquide.

Le terme inertiel ρ ( v → ⋅ g r a d → ) v → {\displaystyle \rho \left({\vec {v}}\cdot {\overrightarrow {\mathrm {grad} }}\right){\vec {v}}} sera donc négligeable devant le terme visqueux μ Δ v → {\displaystyle \mu \Delta \,{\vec {v}}} si :

ρ U 2 L ≪ ρ ν U L 2 {\displaystyle {\frac {\rho U^{2}}{L}}\ll {\frac {\rho \nu U}{L^{2}}}} , soit U L ν ≪ 1 {\displaystyle {\frac {UL}{\nu }}\ll 1} :

où l'on reconnaît l'expression du nombre de Reynolds.

Propriétés des solutions de l'équation de Stokes

Contrairement à l'équation de Navier-Stokes, l'équation de Stokes est linéaire (le terme inertiel, non-linéaire, est en effet négligeable). Les écoulements solutions de cette équation possèdent par conséquent des propriétés bien particulières :

Solution pour une sphère

L'équation de Stokes peut être résolue analytiquement pour calculer la force exercée sur une sphère de rayon r :

F → = − 6 π μ r v → {\displaystyle {\vec {F}}=-6\pi \mu r{\vec {v}}}

Cependant, selon l’habitude des ingénieurs quand il s’agit des sphères et des cylindres, il est plus commode d’utiliser le diamètre D de la sphère que son rayon (les ingénieurs, même en aérodynamique, utilisent très peu le rayon car ils ne savent pas le mesurer) ; la traînée de la sphère en régime de Stokes s’écrit alors :

F → = − 3 π μ D v → {\displaystyle {\vec {F}}=-3\pi \mu D{\vec {v}}}

Au vu de cette équation, il est clair que la traînée de la sphère est proportionnelle à la viscosité dynamique, à la vitesse v et au diamètre D de la sphère. Il est alors naturel de considérer que le scalaire 3π caractérise la forme de la sphère et est le coefficient adimensionnel de traînée de la sphère (en référence à son diamètre D). Horace Lamb en a eu l’intuition en 1911 et a caractérisé la traînée de ses cylindres infinis par un coefficient adimensionnel que Zdravkovich nomme coefficient de résistance de Lamb.

Cependant, afin de relier le résultat de G. G. Stokes pour la sphère avec tous ceux obtenus pour les plus hauts nombres de Reynolds (toujours pour la sphère), Wieselsberger, en dessinant sa courbe (voir ci-dessous la courbe unifiée de Wieselsberger après prise en compte des mesures les plus récentes) a dû adopter, pour la plage de Stokes, le coefficient adimensionnel de traînée classique (que nous qualifierons de quadratique) :

C x = | F | ( 1 / 2 ) ρ v 2 π r 2 {\displaystyle C_{x}={\frac {|F|}{(1/2)\rho v^{2}\pi r^{2}}}}

L’adoption de ce coefficient quadratique conduit, après simplification à :

C x = 24 R e D {\displaystyle C_{x}={\frac {24}{R_{eD}}}}

ReD étant le nombre de Reynolds basé sur le diamètre D = 2r de la sphère, soit :

R e D = ρ v D μ {\displaystyle R_{eD}={\frac {\rho vD}{\mu }}}

Ce choix de Wieselsberger tend cependant, au vu de la courbe unifiée qui en résulte, « à donner l’impression que la résistance aux bas Reynolds prend des valeurs énormes », selon Zdravkovich. Ceci s’explique par le fait que ce Cx quadratique, si utile et éloquent aux forts Reynolds, n’a plus aucune signification en régime de Stokes : en effet, contrairement à ce qui se passe aux hauts Reynolds, la traînée en régime de Stokes n’est nullement proportionnelle au carré de la vitesse, pas plus qu’à la surface frontale du corps et encore moins à la masse volumique du fluide (qui est d’ailleurs censée ne pas agir aux Reynolds très faibles).

Variations du Cx d'une sphère en fonction du nombre de Reynolds.

Au contraire, on peut songer, suivant la proposition d’Horace Lamb, à caractériser les mouvements de la sphère (et de tous les autres corps) en régime de Stokes par un coefficient linéaire adimensionnel de traînée qui sera constant en régime de Stokes. Ce coefficient linéaire adimensionnel de traînée peut être nommé Cx linéaire (mais on peut aussi l’appeler Coefficient de Lamb) et est défini ainsi :

C x Lin Réf.L = | F | μ L v {\displaystyle C_{x{\text{Lin Réf.L}}}={\frac {|F|}{\mu Lv}}}

définition où F est la Traînée (ou la projection de la résultante hydrodynamique sur l’axe des x), μ {\displaystyle \mu } la viscosité dynamique, v est la vitesse sur l’axe des x et L une longueur caractéristique du corps.

Pour la sphère, comme dit précédemment, cette longueur caractéristique gagne à être prise comme le diamètre D, ce qui conduit (en référence à ce diamètre) au Cx linéaire 3π, constant en régime de Stokes.

Pour tous les autres corps dans ce même régime, le Cx linéaire sera de même constant : c’est une des caractéristiques des écoulements de Stokes. Citons par exemple le Cx linéaire du cube qui vaut 4π en référence à son côté (dans toutes les directions) ou les Cx linéaires du disque qui valent 8 et 5,333 en mouvements frontaux et coplanaires (en référence à son diamètre). Quand on connaît le Cx linéaire d’un corps quelconque la Traînée en est alors facilement obtenue en multipliant leur Cx linéaire par la longueur caractéristique qui a présidé à son établissement, par µ et par la vitesse v de l’écoulement.

De la même façon qu’aux hauts Reynolds la surface de référence qui a servi à l’établissement d’un Cx quadratique doit toujours être précisée, la longueur qui a servi à l’établissement d’un Cx linéaire doit être toujours précisée.

Solution pour d'autres corps que la sphère

C x {\displaystyle C_{x}} linéaire de quelques particules en régime de Stokes (cinq tableaux).

L’équation de Stokes a été résolue pour d’autres corps que la sphère. Ainsi Oberbeck a obtenu dès 1876 la valeur exacte de la traînée de l’ellipsoïde de révolution (allongé ou aplati). Ce qui donne, en aplatissant au maximum cet ellipsoïde, le Cx linéaire du disque (8) en déplacements frontaux ou 5,333 en déplacements dans son propre plan (ces deux Cx linéaires en référence au diamètre du disque). De même, en allongeant au maximum l’ellipsoïde, on obtient le Cx linéaire des "aiguilles ellipsoïdales" qui fut longtemps utilisé pour les cylindres de grandes longueurs ou "bâtonnets".

La traînée du tore (même jointif) a été de même calculée, ainsi que celle de la lentille bi-sphérique et des sphères fusionnées, de même que celle des calottes sphériques creuses ou celle de sphères décalottées symétriquement. Ainsi, le Cx linéaire de l’hémisphère creux est-il la moyenne de celui du disque (8) et celui de la sphère (3π) soit 8,7124 (en référence au diamètre de la sphère intègre et du disque).

Cox a calculé la traînée des « bâtonnets cylindriques » (c.-à-d. des cylindres d’élancement Longueur / Diamètre assez grand) en déplacements transverses ou axiaux. Il a étendu sa démonstration au cône et au bicône (toujours d’élancements assez forts).

Bartuschat et coll. ont plus récemment calculé, par la méthode lattice Boltzmann 3D, les caractéristiques de traînée de corps hémisphéro-cylindriques d’élancement 2 à 7 en déplacements transverses ou axiaux. De telles méthodes mixant les principes mathématiques et les moyens informatiques ont donné accès aux caractéristiques de traînée (toujours en régime de Stokes) des courts cylindres (Roger et Ui) ou des plaques rectangulaires ou carrées (Sunada, Tokutake et Okada) ; le Cx linéaire de la plaque carrée mince en ressort comme valant, en référence à son côté, 9,15 pour des déplacements frontaux et 6,13 pour des déplacements dans son propre plan (déplacements coplanaires).

Les caractéristiques de traînée des prismes rectangulaires ont été approchées par l'expérience (Heiss & Coull, Malaika). Plus récemment, S. Sunada, R. Ishida, H. Tokutake et S. Okada ont déterminé ces caractéristiques par le calcul. Une mise en panorama de ces résultats montre que le problème de la traînée des prismes rectangulaires est sur le point d'être résolue.

E. O. Tuck, mathématicien australien, a inventé une méthode inverse permettant de dessiner au hasard des corps de révolution quelconques (de grand élancement L/D) dont on connaît préalablement la traînée. Certains de ces corps de révolution quelconques peuvent être très proches de corps géométriques, comme les corps à génératrice circulaire ou les corps hémi-ellipsoïdo-cylindriques (dont les corps hémisphéro-cylindriques sont des cas particuliers).

La traînée du cube a été déterminée de façon satisfaisante par l’expérience (Cx linéaire = 4π, en référence à son côté).

Pour les corps de révolution (ou les corps axisymétriques comme les prismes à base carrée) d'assez grand élancement L/D, il est souvent dit que leur traînée transverse est proche du double de leur traînée axiale (à fluide et à vitesse identiques). Cela n'est vrai que pour les très forts élancements (ce quotient de traînées n'atteint ~ 1,75 que pour un cylindre circulaire ou un ellipsoïde d'élancement 1000, par exemple).

En ce qui concerne les corps 2D comme les cylindres de longueur infinie de section circulaire ou elliptique ou les palettes (également de longueur infinie), les mathématiciens, butant sur le paradoxe de Stokes, n’ont pu en déterminer la traînée qu’en régime décrit par les équations de Stokes-Oseen (pour les Reynolds autour de l’unité où l’inertie commence à se faire sentir ).

Autres solutions de l'équation de Stokes

La simplicité relative de la loi de Stokes autorise le changement d'échelle (prise de moyenne volumique ou homogénéisation), conduisant à des résultats théoriques importants :

Notes

  1. Voir par exemple l'image : Cx linéaire des prismes rectangulaire en régime de Stokes d'après Sunada et coll.

Références

  1. (en) E. M. Purcell, « Life at low Reynolds number », American Journal of Physics 45, p. 3-11.
  2. (en) George K. Batchelor, An Introduction to Fluid Mechanics, Cambridge/New York, Cambridge University Press, 2000, 615 p. (ISBN 0-521-66396-2)
  3. (en) Momchilo M. Zdravkovich, « A critical remark on use of drag coefficient at low reynolds numbers », Recueil des travaux de l'Institut Mathématique, Nouvelle série, t. 3, no 11,‎ 1979 (lire en ligne ).
  4. (en) Carl Wieselsberger, « Further information on the laws of fluid resistance », NACA T. N. N°121,‎ 1922 (lire en ligne).
  5. (en) Stephen Whitaker, The Method of Volume Averaging, Kluwer Academic Publishers, 2010, 471 p. (ISBN 978-3-642-05194-4).
  6. Y. Achdou, P. Le Tallec, F. Valentin et O. Pironneau, « Constructing wall laws with Domain Decomposition or asymptotic expansion techniques », Computer Methods in Applied Mechanics and Engineering, vol. 151, nos 1-2,‎ 1998, p. 215-232.

Bibliographie

Articles connexes