Écoulement de Poiseuille
Un article de Wikivisual, l'encyclopédie libre.
Un fluide visqueux, s'il est en écoulement lent dans un tuyau de petit diamètre ou entre deux plaques proches, est en écoulement de Stokes.
En première approximation, si le tuyau est cylindrique ou que les plaques sont parallèles, l'écoulement du fluide est partout parallèle aux parois (approximation de lubrification).
Le frottement aux parois implique qu'aux échelles macroscopiques, la vitesse du fluide y est nulle (condition de non-glissement).
Par ailleurs, la pression ne varie pas dans l'épaisseur de l'écoulement (approximation de lubrification).
Ces trois conditions impliquent que l'écoulement s'organise selon un champ de vitesse parabolique : vitesse nulle aux parois et maximale à mi-hauteur.
Ci-dessous, on considère deux problèmes différents qui donnent lieu à un écoulement de Poiseuille :
- l'écoulement dans un tube de section circulaire et de rayon constant <math> R </math>,
- l'écoulement entre deux plaques planes et parallèles, distantes de <math> h </math> ; ce calcul permet notamment d'évaluer la force entre deux objets (par exemple deux disques) immergés dans un fluide visqueux et s'approchant à une vitesse donnée.
[modifier] Champ de vitesse dans un tube
La vitesse est parallèle à l'axe du tube (noté <math> z </math>) : <math> \vec{v} = v \vec{u}_z </math>.
Équation du profil de vitesse :
<math> v(r,z,\theta) = v(r) = v_{\rm max}\;\left( 1-\frac{r^2}{R^2} \right) </math>
où la vitesse maximale (au centre du tube) est liée au gradient de pression, à la viscosité et au rayon : <math> v_{\rm max} = \frac{R^2}{4\;\eta} \; \frac{{\rm d} p}{{\rm d} z} </math>
La démonstration de ce résultat est donnée plus bas.
[modifier] Champ de vitesse entre deux plaques
On suppose que le gradient de pression est orienté selon l'axe <math> x </math> et que la normale aux plaques est orientée selon <math> z </math>, avec les plaques situées en <math> z = \pm h/2</math>. La vitesse est alors parallèle aux plaques, et plus précisément orientée selon l'axe <math> x </math> : <math> \vec{v} = v \vec{u}_x </math>.
Équation du profil de vitesse :
<math> v(x,y,z) = v(z) = v_{\rm max}\;\left( 1-\frac{4\,z^2}{h^2} \right) </math>
où la vitesse maximale (au centre du tube) est liée au gradient de pression, à la viscosité et à la distance entre les plaques : <math> v_{\rm max} = \frac{h^2}{8\;\eta} \; \frac{{\rm d} p}{{\rm d} x} </math>
La démonstration de ce résultat est similaire à celle donnée ci-dessous dans le cas du tube circulaire.
[modifier] Démonstration (dans le cas du tube)
1. Par symétrie, l'écoulement ne varie ni en <math> z </math>, ni en <math> \theta </math> :
<math> v(r,z,\theta) = v(r) </math>
2. Par conséquent, les seuls efforts de cisaillement sont des forces selon <math> z </math>
transmises radialement (selon <math> r </math>) :
<math> \sigma_{rz}(r,z,\theta) = \sigma_{rz}(r) = \eta\;\frac{{\rm d} v(r)}{{\rm d} r} </math>
3. Par symétrie également, la variation de la pression est constante le long de l'axe <math> z </math> :
<math> \frac{{\rm d} p}{{\rm d} z} = {\rm const} </math>
4. Considérons les efforts subis par une zone cylindrique
de rayon <math> r </math> et de longueur <math> \Delta z </math>.
Les efforts de pression sur les deux faces circulaires du cylindre ont une résultante égale à :
<math> F_{\rm faces} = \pi\,r^2 \; \Delta z \; \frac{{\rm d} p}{{\rm d} z} </math>
Les contraintes de cisaillement sur le bord du cylindre lui transmettent une force orientée selon son axe <math> z </math> :
<math> F_{\rm bord} = 2\pi\,r \; \Delta z \; \sigma_{rz}(r) </math>
La force totale exercée sur le cylindre de liquide est nulle puisque l'écoulement est permanent. Ainsi :
<math> \sigma_{rz}(r) = \frac{r}{2} \; \frac{{\rm d} p}{{\rm d} z} </math>
5. Il s'ensuit que le gradient de vitesse est linéaire en <math> r </math> :
<math> \frac{{\rm d} v(r)}{{\rm d} r} = \frac{\sigma_{rz}(r)}{\eta} = \frac{r}{2\;\eta} \; \frac{{\rm d} p}{{\rm d} z} </math>
6. Autrement dit, le champ de vitesse est parabolique :
<math> v(r) = {\rm const} + \frac{r^2}{4\;\eta} \; \frac{{\rm d} p}{{\rm d} z} </math>
7. Compte tenu de la condition de non-glissement (<math> v(R) = 0 </math>) :
<math> v(r) = -\frac{R^2}{4\;\eta} \; \frac{{\rm d} p}{{\rm d} z} \; \left( 1-\frac{r^2}{R^2} \right) </math>
La vitesse est plus importante au centre du conduit malgré le signe négatif, étant donné que la vitesse est orientée à l'encontre du gradient de pression. Écoulement dans le sens positif pour un gradient négatif... CQFD
