Mécanique

Chapitre 6

Écoulements de fluides
PC⋆ , Fabert (Metz)

Écoulements de fluides
Après avoir vu quelques mouvements particuliers de fluide dans le chapitre précédent, nous allons
désormais nous intéresser à l’aspect local des écoulements. Après avoir déterminé les équations du
mouvement en terme de champ de vitesse, nous déterminerons explicitement ce dernier dans quelques
cas particuliers : pour des écoulements dits visqueux dans la première partie et pour des écoulements
parfaits dans la deuxième.
Ces calculs ne devront pas faire oublier qu’ils ne s’agit là que de cas très particuliers et que, dans
le cas général, il est impossible de trouver analytiquement l’expression du champ de vitesse d’un
écoulement.

© Matthieu Rigaut 2 / 70 Version du 11 mars 2014

PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

Table des matières

Biographies succinctes 6

I Écoulements visqueux 9
I·1 Actions au sein d’un fluide . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
I·1·i rappel du modèle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
I·1·ii modélisation des forces tangentielles : fluide newtonien . . . . . . . . . . . . . 9
situation canonique : écoulement laminaire 1D . . . . . . . . . . . . . . . . . 9
de l’interprétation à l’expression . . . . . . . . . . . . . . . . . . . . . . . . . 10
I·1·iii viscosité dynamique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
unité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
valeurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
I·1·iv équivalent volumique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
résultante des forces de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . 11
situation canonique : écoulement laminaire 1D . . . . . . . . . . . . . . . . . 12
I·2 Équation de Navier – Stokes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
I·2·i le PFD renommé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
« démonstration » . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
I·2·ii les conditions aux limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
constatation expérimentale . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
traduction technique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
I·2·iii 1re interprétation du nombre de Reynolds : accélération convective et effet
de viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
termes de l’équation de Navier – Stokes . . . . . . . . . . . . . . . . . . . 14
viscosité cinématique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
I·3 Écoulement ce Couette plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
I·3·i étude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
dispositif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
première simplification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
mise en équation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
le régime stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
I·3·ii 2e interprétation du nombre de Reynolds : épaisseur de la couche limite . . 18
deux situations équivalentes . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
où l’équation de diffusion refait une apparition . . . . . . . . . . . . . . . . . 19
exemple numérique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
I·4 Écoulement de Poiseuille plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
I·4·i qu’est-ce qu’un écoulement de Poiseuille ? . . . . . . . . . . . . . . . . . . 21
I·4·ii écoulement de Poiseuille plan . . . . . . . . . . . . . . . . . . . . . . . . . 21
dispositif, analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
mise en équation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
résolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
profil des vitesses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
débit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
interprétation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
I·4·iii 3e interprétation du nombre de Reynolds : transport de la quantité de mou-
vement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
deux types de transport . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
transport convectif sur l’exemple de Poiseuille plan . . . . . . . . . . . . . 25
© Matthieu Rigaut 3 / 70 Version du 11 mars 2014
PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

transport diffusif sur l’exemple de Poiseuille plan . . . . . . . . . . . . . . 25

comparaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
un terme diffusif étrange . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
I·5 écoulement de Poiseuille cylindrique . . . . . . . . . . . . . . . . . . . . . . . . . . 27
I·5·i champ des vitesses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
dispositif, analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
liminaires . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
mise en équation par le PFD – version locale . . . . . . . . . . . . . . . . . . 28
mise en équation par le TCI – version globale . . . . . . . . . . . . . . . . . . 30
I·5·ii quelques aspects de l’écoulement de Poiseuille cylindrique . . . . . . . . . 31
débit volumique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
pourquoi une dépendance en R4 ? . . . . . . . . . . . . . . . . . . . . . . . . 32
exemple qualitatif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
une nouvelle analogie diffusive . . . . . . . . . . . . . . . . . . . . . . . . . . 32
puissance perdue par viscosité . . . . . . . . . . . . . . . . . . . . . . . . . . 33

II Écoulements parfaits 35
II·1 Modèle de l’écoulement parfait . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
II·1·i définition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
II·1·ii en pratique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
II·2 Équation d’Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
II·2·i à partir de l’équation de Navier – Stokes . . . . . . . . . . . . . . . . . . 35
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
« démonstration » . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
II·2·ii lecture . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
équation non linéaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
retrouver la RFSF . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
tout est désormais possible . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
II·3 Premières applications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
II·3·i surface libre d’un tourbillon . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
dispositif, modèle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
idée, plan de bataille . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
pour r < a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
pour r > a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
détermination complète du champ de pression . . . . . . . . . . . . . . . . . 41
surface libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
retour sur l’accélération particulaire . . . . . . . . . . . . . . . . . . . . . . . 42
II·3·ii pression et forme des lignes de courant . . . . . . . . . . . . . . . . . . . . . 43
répartition de pression . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
effet Coanda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
II·3·iii pression dans un jet libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
jet libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
résultat utile . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
jet libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
II·4 Ondes sonores, l’écho . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
II·4·i rappel à la loi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
II·4·ii linéarisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
approximation acoustique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
l’accélération convective est négligeable . . . . . . . . . . . . . . . . . . . . . 47
l’équation d’Euler . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
© Matthieu Rigaut 4 / 70 Version du 11 mars 2014
PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

conservation de la masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
loi de comportement phénoménologique . . . . . . . . . . . . . . . . . . . . . 48
II·4·iii équation de propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
équations de couplage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
équation de propagation en surpression . . . . . . . . . . . . . . . . . . . . . 48
équation de propagation en vitesse particulaire . . . . . . . . . . . . . . . . . 49
II·5 Relations de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
II·5·i écoulement irrotationnel et homogène . . . . . . . . . . . . . . . . . . . . . . 49
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
II·5·ii écoulement non irrotationnel et homogène . . . . . . . . . . . . . . . . . . . . 50
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
démonstration sans l’homogénéité . . . . . . . . . . . . . . . . . . . . . . . . 52
remarques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
généralisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
II·5·iii idoinotons vérificateurs d’hypothèses . . . . . . . . . . . . . . . . . . . . . . . 53
importance du homogène . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
importance du irrotationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
II·5·iv limite de vitesse pour un écoulement incompressible . . . . . . . . . . . . . . 55
II·6 Applications des relations de Bernoulli . . . . . . . . . . . . . . . . . . . . . . . . . 56
II·6·i effet Venturi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
qualitativement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
débimètre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
II·6·ii vidange d’un réservoir : formule de Torricelli . . . . . . . . . . . . . . . . 60
expérience . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
cas quasi-stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
vérification de la quasi-stationnarité . . . . . . . . . . . . . . . . . . . . . . . 62
II·6·iii tube de Pitot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
dispositif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
modélisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
pressions en A et B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
mesure de vitesse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
II·7 Régime non stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
II·7·i idée . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
II·7·ii oscillations dans un tube en U . . . . . . . . . . . . . . . . . . . . . . . . . . 65
dispositif . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
équation d’évolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68

Fiche de révision 69

© Matthieu Rigaut 5 / 70 Version du 11 mars 2014

PC⋆ , Fabert (Metz) Biographies succinctes

Biographies succintes

Evangelista Torricelli
(1608 Faenza – 1647 Florence)

Fils d’un modeste artisan, Evangelista étudie pendant deux ans à l’école des
Jésuites avant d’être envoyé à Rome sur recommandation de ses professeurs. Il suit
alors les cours d’un ancien élève de Galilée, Castelli. Ce dernier remarque le
talent d’Evangelista et le prend comme secrétaire avant de le proposer à Galilée
comme assistant personnel. En 1641 le savant accepte mais meurt en 1642. C’est
alors Evangelista qui prend sa place comme philosophe et mathématicien privé
du grand duc de Toscane, ce qu’il restera jusqu’à la fin de sa vie.

Henri Pitot
(1695 Aramon, Languedoc – 1771 Aramon)

Alors qu’il n’aime pas les études, Henri rentre en tant qu’assistant dans le la-
boratoire de Réaumur à Paris. Il s’intéresse alors aux mathématiques et à la
mécanique et en vient à diriger la construction d’un canal. Par la suite il conti-
nue sa carrière d’ingénieur en se préoccupant à la fois des aspects théoriques et
pratiques. Il invente pour le cela le tube qui porte son nom.

Leonhard Euler
(1707 Bâle – 1783 Saint-Pétersbourg)

Fils de pasteur, Leonhard ne rencontre que peu de difficulté scolaire puisqu’il

rentre à l’université de Bâle à 14 ans et publie son premier article à 19 ans.
Aidé par son professeur, Johann Bernoulli dont le fils n’est autre que Daniel, il
rejoindra ce dernier à l’université de Saint-Pétersbourg pour y effectuer l’essentiel
de sa carrière (il passera quelques années à Berlin aussi). Leonhard est à la fois un
grand physicien et un des plus grands mathématiciens de l’histoire. Il travaillera
entre autre sur l’optique, l’hydrodynamique et la mécanique.

Giovanni Battista Venturi

(1746 Reggio Emilia – 1822 Reggio Emilia)

Giovanni est ordonné prêtre en 1769 et commence aussitôt à enseigner la logique.

Vers 1773 il devient professeur de géométrie et philosophie à l’université de Mo-
dène puis, en 1776, il devient professeur de physique. Il s’installe à Paris à partir
de 1796 et c’est là qu’il écrit un ouvrage où il décrit l’effet qui porte son nom.

© Matthieu Rigaut 6 / 70 Version du 11 mars 2014

PC⋆ , Fabert (Metz) Biographies succinctes
Claude Louis Marie Henri Navier
(1785 Dijon – 1835 Paris)

Claude, jeune, suit son père qui travaille à l’assemblé constituante mais ce dernier
meurt en 1793. Claude est alors pris en charge par son oncle ingénieur. Il rentre
à l’école Polytechnique en 1802 et à l’école des Ponts et Chaussées en 1804.
Ses résultats sur l’application des mathématiques aux sciences de l’ingénieur,
notamment en ce qui concerne la déformation élastique des solides, lui permettent
d’obtenir, en 1819, une charge de cours à l’école des Ponts et Chaussées. Il se
tourne alors vers l’écoulement des liquides et découvre en 1822 l’équation qui
porte son nom et qui sera redécouverte indépendamment par Stokes en 1845.

Jean-Louis Marie Poiseuille

(1797 Paris – 1869 Paris)

Fils d’un charpentier, Jean-Louis est élève de l’école Polytechnique de 1815 à

1816. Sa thèse Recherches sur la force aortique du coeur soutenue en 1828 est très
remarquée et même couronnée par la médaille d’or de l’Académie des Sciences car
c’est dans celle-ci que, pour la fois, est mesurée la pression sanguine. Dans la suite,
Jean-Louis s’intéresse principalement à la physiologie de la circulation sanguine.
Cela l’amène non seulement à s’intéresser à la physique des écoulements. En 1842
il est élu à l’Académie de médecine.

Heinrich Gustav Magnus

(1802 Berlin – 1870 Berlin)

Fils d’un riche marchant, l’un des cinq frères d’Heinrich (Eduard 1799 - 1872)
deviendra un célèbre peintre. Heinrich étudie à Berlin, Stokholm puis à Paris avec
Gay-Lussac. Il obtient son doctorat en 1831 et devient aussitôt présentateur. Il
est ensuite nommé assistant professeur puis professeur en 1834 et 1845. Heinrich
commence par s’intéresser à la chimie (sels de Magnus) avant d’étudier les gaz,
la vapeur d’eau, l’induction, l’aérodynamique, . . . Il prend sa retraite en 1869.

George Stokes
(1819 Sligo, Islande – 1903 Cambridge)

George fait ses études à Cambridge où il obtient son diplôme en 1841 pour y de-
venir professeur de mathématique en 1849. Son activité de recherche se concentre
d’abord sur l’hydrodynamique des fluides visqueux. Il s’intéresse après à la pro-
pagation du son et de la lumière. Il explique le phénomène de fluorescence vers
1852. Il est à noter que sa productivité de chercheur a notablement diminué à
partir de 1857, année de son mariage.
© Matthieu Rigaut 7 / 70 Version du 11 mars 2014
PC⋆ , Fabert (Metz) Biographies succinctes
Osborne Reynolds
(1842 Belfast – 1912 Watchet, Grande-Bretagne)

Le père d’Osborne fut, d’après ce dernier, son premier professeur. Directeur

d’école et prêtre, il a déposé de nombreux brevets pour améliorer des machines
agricoles. Osborne fait ses études à Cambridge et devient quelques temps après
l’un des premiers « professeurs en ingénieurie » dans une nouvelle école qui de-
viendra l’université de Manchester. En parallèle de ses cours, il fait des recherches
en mécaniques des fluides et c’est en 1883 qu’il publie un article dans lequel il
introduit le nombre qui porte son nom.

Maurice Couette
(1858 Tours – 1943 Dijon)

Maurice né en 1858 et est le seul fils de la famille dont le père est marchand de
vêtements. Il fait ses études chez les frères des écoles chrétiennes et obtient deux
baccalauréats (en science et en humanités) la même année 1874. Il suit ensuite
deux licences de mathématiques et de physique qu’il obtient en 1877 et 1879. En
1881 il s’inscrit à la Sorbonne pour préparer l’agrégation. Dans la suite il enseigne
un peu partout et notamment dans l’université catholique de l’ouest (Anger). Il
sera membre de la société de physique.

Henri Coanda
(1886 Bucarest – 1972 Bucarest)

Le père d’Henri est général, professeur de mathématiques et attaché d’embassade

en France. C’est pourquoi Henri commence ses études en France mais est envoyé
à 13 ans en Roumanie dans un lycée militaire. Sa passion pour tout ce qui touche
au vent et à l’aviation le conduit à contruire un avion-fusée pour l’armée en 1905.
Toutefois il quitte l’armée et voyage à travers le monde. À son retour, il reprend
des études d’ingénieur à Supaéro et en sort major en 1910. Il travaillera ensuite
dans l’aviation jusqu’en 1969 où il retourne en Roumanie pour la fin de sa vie.

© Matthieu Rigaut 8 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I – Écoulements visqueux

I – Écoulements visqueux
I·1 – Actions au sein d’un fluide
I·1·i – rappel du modèle
✧ Considérons une particule de fluide.

✧ Cette particule de fluide est soumise à deux types de forces :

➜ les forces à distances ;
➜ les forces de contact.
✧ En ce qui concerne les forces à distances qui sont « naturellement » volumique, nous trouvons :
➜ la pesanteur ;
➜ la force de Laplace ;
➜ par extension, les forces d’inertie.
✧ Au niveau des forces de contact, nous avons :
−−→
➜ les actions normales, pressantes, dont la résultante s’écrit df~press = −grad P dτ ;
➜ les action tangentielles associées au frottement.

I·1·ii – modélisation des forces tangentielles : fluide newtonien

✬ situation canonique : écoulement laminaire 1D
✧ Considérons un écoulement laminaire unidimensionnel, i.e. un écoulement où toutes les lignes de
courant sont parallèles et ne dépendent que d’une seule coordonnée cartésienne de l’espace du type
~v (M,t) = v (y,t) ~ux .

~uy
~ux

✧ Dans ces conditions, nous pouvons isoler, par la pensée, une particule de fluide (la ② sur le schéma
ci-dessous) et regarder ce qui va se passer.
© Matthieu Rigaut 9 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·1 – Actions au sein d’un fluide

③
~uy
~ux

✧ Imaginons que, comme le représente la figure ci-dessus, la vitesse soit croissante avec y. Alors :
➜ la particule de fluide ①, qui a une vitesse plus grande que la ②, va entraîner cette dernière vers
l’avant, i.e. elle va exercer une force suivant +~ux ;
➜ la particule de fluide ③, qui a une vitesse plus faible que la ②, va ralentir cette dernière et,
donc, va exercer une force suivant −~ux .
✧ Au delà de cette considération qualitative, nous pouvons dire, de manière naturelle, que la force de
frottement entre particules de fluide dépend :
➜ de l’écart de vitesse entre les deux particules, i.e. de la manière dont varie la vitesse dans
l’espace ;
➜ de la surface de contact entre les deux particules de fluide (plus la surface est grande, plus
l’interaction l’est).
✧ Ensuite, nous pouvons naturellement faire l’approximation que cette dépendance, au moins au pre-
mier ordre, est linéaire. Autrement dit, nous pouvons dire que la force est proportionnelle à la variation
spatiale de la vitesse et à la surface de contact entre les deux particules de fluide.

✬ de l’interprétation à l’expression

Dans la situation représentée ci-dessous, la force que la particule du haut exerce sur celle
du bas s’écrit
∂vx
f~h→b = +η S ~ux où :
∂y
➜ η est la viscosité dynamique (en Pl) ;
➜ S est la surface de contact entre les deux particules de fluide.

~uy
~ux

! Attention ! La forme exacte de la force n’est pas à connaître par cœur mais à savoir retrouver car
elle dépend :
➜ du repérage ;
➜ de qui exerce une force sur qui.

Un fluide tel que les forces tangentielles qui s’exercent entre particules de fluide
∂vx
s’écrivent « η S » est dit newtonien.
∂y

© Matthieu Rigaut 10 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·1 – Actions au sein d’un fluide

✧ En fait la plupart des fluides usuels sont newtoniens : l’eau, l’air, l’huile, le glycérol. . .

I·1·iii – viscosité dynamique

✬ unité

η est appelé la viscosité dynamique.

✧ η est aussi appelée plus simplement « viscosité », l’ajout du qualificatif « dynamique » n’étant là que
pour la discrimination avec « viscosité cinématique » que nous verrons plus loin.

L’unité de la viscosité dynamique η est le poiseuille, noté Pl ou, parfois, le Pa.s.

✧ Pour le montrer, partons de l’expression de la force, sachant qu’une force est une pression multipliée
par une surface.

∂vx ∂vx η ℓ
f ≡η× S≡PS η× ≡P ≡P× ≡Pτ
∂y ∂y v
✧ Rappelons aussi, pour la culture, que la force de traînée exercée par un fluide sur une sphère s’écrit
(formule de Stokes)

f~ = −6 π η r ~v avec :

➜ r le rayon de la sphère ;
➜ ~v la vitesse de la sphère par rapport au fluide ;
➜ η la viscosité dynamique.

✬ valeurs
✧ Le poiseuille est une « grande » unité car, un peu comme le tesla, l’orde de grandeur usuel est plutôt
en sous-multiple du poiseuille.

fluide air eau huile glycérol

η (Pl) 2.10 −5
10 −3
0,1 1

I·1·iv – équivalent volumique

✬ résultante des forces de viscosité

La résultante des forces de viscosité qui s’exerce sur une particule de fluide de volume dτ
dans un écoulement incompressible s’écrit

~ ~v dτ
df~visc = η △

© Matthieu Rigaut 11 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·1 – Actions au sein d’un fluide

✬ situation canonique : écoulement laminaire 1D

✧ Nous allons démontrer le résultat précédent dans un cas particulier et admettre la généralisation.
✧ Reprenons l’écoulement laminaire unidimensionnel précédent et cherchons la résultante des forces de
viscosité s’exerçant sur une particule de fluide.

y + dy

~uy
~ux

✧ Intéressons-nous à la particule de fluide en jaune et cherchons la résultante des forces de viscosité

s’exerçant sur elle.
✧ Il y a, a priori, 6 particules de fluide qui touchent PF, donc qui exercent des forces sur elle :
➜ les particules de fluide devant (contact en x + dx) et derrière (contact en x) vont à la même
vitesse (puisque celle-ci ne dépend que de y), donc elles ne peuvent pas exercer de force de
frottement sur PF ;
➜ il en est de même avec les particules de fluide à droite (contact en z + dz) et à gauche (contact
en z) ;
➜ seules les deux particules de fluide coloriées en bleu clair au-dessus (contact en y + dy) et
au-dessous (contact en y) exercent des forces.
✧ En notant ① et ② les particules de fluide au-dessus et en-dessous, nous avons

df~visc = df~1 →PF + df~2 →PF

✧ Or, en faisant très attention aux signes des forces, nous avons :
∂vx ∂vx
df~1 →PF = +η (y + dy) dS ~
ux et df~1 →PF = −η (y) dS ~
ux
∂y ∂y
✧ Étant donné que la surface de contact sur le dessus et le dessous est la même et s’écrit dS = dx dz
nous avons

df~visc = df~1 →PF + df~2 →PF

∂vx ∂vx
= +η (y + dy) dx dz ~
ux − η (y) dx dz ~
ux
∂y ∂y
!
∂vx ∂vx
= η dx dz (y + dy) − (y) ~ux
∂y ∂y
!
∂ ∂vx
= η dx dz × (y) dy ~
ux
∂y ∂y
∂ 2 vx
= η dx dy dz × (y) ~
ux
∂y 2
✧ Ce qui donne bien
© Matthieu Rigaut 12 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·2 – Équation de Navier – Stokes

∂ 2 vx
df~visc = η dτ × (y) ~
ux
∂y 2

I·2 – Équation de Navier – Stokes

I·2·i – le PFD renommé
✬ énoncé

équation de Navier – Stokes

Dans un fluide newtonien, pour un écoulement incompressible, le champ de vitesse obéit
aux équations (équivalentes) suivantes :

d~v −−→ ~ ~v
µ = −grad P + f~v,tot + η △
dt
  !
∂~v −−→ −−→ ~ ~v
µ + ~v · grad ~v = −grad P + f~v,tot + η △
∂t
!
∂~v −−→ v 2 − →  −−→ ~ ~v
µ + grad + rot ~v ∧ ~v = −grad P + f~v,tot + η △
∂t 2

➜ µ est la masse volumique ;

➜ f~v,tot est la densité volumique de force à distance (poids, électromagnétisme, inertie) ;
➜ η est la viscosité dynamique.

✬ « démonstration »
✧ La démonstration est quasiment immédiate.
✧ Écrivons le PFD sur une particule de fluide de volume dτ .
✧ Cette particule de fluide subit les forces :
➜ force à distance df~dist = f~v,tot dτ ;
−−→
➜ force de contact normales (les forces pressantes) df~press = −grad P dτ ;
~ ~v dτ .
➜ force de contact tangentielles (les forces de viscosité) df~visc = η △
✧ Le PFD donne donc

a = df~press + df~dist + df~visc

dm ~e
✧ Ce qui s’écrit, en remplaçant
−−→ ~ ~v dτ
a = −grad P dτ + f~v,tot dτ + η △
µ dτ ~e
✧ Ce qui donne bien le résultat en simplifiant par dτ car l’accélération lagrangienne n’est autre que la
dérivée particulaire de la vitesse
 
e d~v ∂~v −−→ ∂~v −−→ v 2 − → 
~a = = + ~v · grad ~v = + grad + rot ~v ∧ ~v
dt ∂t ∂t 2

© Matthieu Rigaut 13 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·2 – Équation de Navier – Stokes

I·2·ii – les conditions aux limites

✧ Nous avons déjà rencontré les conditions aux limites naturelles qui ne faisaient que traduire l’imper-
méabilité des parois.
✧ Pour les fluides newtoniens, il existe d’autres conditions aux limites.

✬ constatation expérimentale

Les particules de fluide d’un fluide newtonien adhèrent aux parois.

✧ Il s’agit là d’une constatation expérimentale, nous ne pouvons rien y faire, si ce n’est l’admettre et
le traduire analytiquement.

✬ traduction technique
✧ La traduction est simple : puisque la particule de fluide adhère à la paroi, cela signifie qu’elle (la
particule de fluide) a tout le temps la même vitesse que la paroi.

Pour un fluide newtonien, en un point I au niveau d’une paroi

~v (I ∈ fluide)|R = ~v (I ∈ paroi)|R

✧ Très fréquemment (mais pas systématiquement) la paroi sera immobile, ce qui donnera

Pour un fluide newtonien, en un point I au niveau d’une paroi immobile

~v(I ∈ fluide)|R = ~0

I·2·iii – 1re interprétation du nombre de Reynolds : accélération convec-

tive et effet de viscosité
✬ termes de l’équation de Navier – Stokes
✧ Reprenons l’équation de Navier – Stokes
  !
∂~v −−→ −−→ ~ ~v
µ + ~v · grad ~v = −grad P + f~v + η △
∂t
✧ Il y a deux
 termestout particuliers :
−−→
➜ µ ~v · grad ~v est l’accélération convective ;
~ ~v est l’action de viscosité.
➜ η△

Le nombre de Reynolds peut s’écrire

accélération convective
Re =
action de viscosité

✧ La relation précédente est à prendre, bien sûr, en terme d’ordre de grandeur.

✧ Nous avons ainsi
© Matthieu Rigaut 14 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·2 – Équation de Navier – Stokes


  

 −−→ V
accélération convective = µ ~v · grad ~v ≡ µ × V ×
L

 ~ V
 action de viscosité = η △ ~v ≡ η L2

✧ En regroupant, cela donne

V
accélération convective µ×V × accélération convective µV L
= L = = Re
action de viscosité V action de viscosité η
η 2
L

✧ Nous pouvons donc dire que si Re est petit, alors l’action de viscosité prédomine sur l’accélération
convective, ce qui se traduit par
 
−−→ ~ ~v
µ ~v · grad ~v ≪ η △

✬ viscosité cinématique

La viscosité cinématique ν est définie par

η
µ,
µ

✧ Dans ces conditions nous avons immédiatement

Le nombre de Reynolds s’écrit

VL
Re =
ν

✧ L’intérêt de la viscosité cinématique provient du fait que, quand le fluide change, la viscosité dyna-
mique η change en même temps que la masse volumique µ. En d’autre termes, il n’est pas facile de
changer notablement η sans changer µ et réciproquement et c’est pourquoi il vaut mieux parler du
η
rapport ν = .
µ
✧ Nous avons ainsi

fluide η (Pl) ν (m2 .s−1 )

eau 10−3 10−6
air 1,8.10−5 1,4.10−5

✧ Nous voyons donc que, bien que l’eau soit considérablement plus visqueuse que l’air, c’est quand
même elle qui va avoir tendance à créer des écoulement à grand nombre de Reynolds, i.e. des
écoulement où la viscosité est négligeable. Cela provient, évidemment, de son inertie, qui est, elle
aussi, énorme.
© Matthieu Rigaut 15 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·3 – Écoulement ce Couette plan

I·3 – Écoulement ce Couette plan

I·3·i – étude
✬ dispositif
✧ Considérons deux plans infinis, l’un immobile et l’autre mobile, entre lesquels s’écoule un fluide.

~v0

e ~g

✧ Cette situation peut fort bien correspondre à un zoom sur l’espace entre deux pièces mécaniques
lubrifiées.
✧ Nous allons chercher, ici, l’équation régissant le champ de vitesse et, si possible, nous allons trouver
explicitement celui-ci.
✧ Au point de vue des grandeurs pertinentes nous avons :
➜ l’inertie µ ;
➜ la géométrie e ;
➜ les actions g pour la pesanteur, η pour la viscosité ;
➜ v0 pour les contraintes.
✧ N’oublions pas de repérer le tout.

~v0

e ~g
~uy
~ux
O

✬ première simplification
✧ Pour simplifier un peu l’ensemble, nous allons considérer :
➜ que les plans sont infinis en ~ux et ~uz ;
➜ que l’écoulement est laminaire, i.e. que les lignes de courant sont portées uniquement par ~ux .
✧ Dans ces conditions, le champ de vitesse s’écrit

~v (M,t) = vx (M,t) ~ux ~v (M,t) = vx (✚ ✁ ~

x,y,z,t) ux ~v (M,t) = vx (y,t) ~ux
✧ De même l’invariance par translation suivant ~ux et ~uz nous permet de dire que

P (M,t) = P (y,t)

© Matthieu Rigaut 16 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·3 – Écoulement ce Couette plan

✬ mise en équation
✧ Commençons par reprendre l’équation de Navier – Stokes
  !
∂~v −−→ −−→ ~ ~v
µ + ~v · grad ~v = −grad P + f~v + η △
∂t

✧ Ici, la seule force à distance qui s’exerce est la pesanteur, ce qui donne
  !
∂~v −−→ −−→ ~ ~v
f~v,tot = µ ~g µ + ~v · grad ~v = −grad P + µ ~g + η △
∂t

✧ Projetons sur ~ux l’équation précédente. Comme nous sommes en coordonnées cartésiennes, l’opération
se fait facilement et nous avons, puisque le gradient de pression n’est porté que par ~uy ,
  !
∂vx −−→
µ (y,t) + ~
v · grad vx (y,t) = 0 + 0 + η △ vx (y,t)
∂t

✧ De plus
   
   
 vx   0 
   
   
   
−−→     ∂ 
~v · grad =  0  ·  
  ∂y  = 0
   
   
   
 0   0 
   

✧ Mais aussi

∂ 2 vx
△ vx = 0 + (y,t) + 0
∂y 2
✧ Il reste donc

∂vx ∂ 2 vx
µ (y,t) = η (y,t)
∂t ∂y 2
✧ Il s’agit ni plus ni moins que d’une équation de diffusion, autrement dit :
➜ le régime stationnaire est « lent » à venir ;
➜ il n’y a pas forcément de solution analytique.

✬ le régime stationnaire
Y simplification
✧ Supposons, ici et seulement ici, que le régime soit stationnaire.
✧ Dans ces conditions, le champ des vitesses obéit à l’équation

d2 vx d2 vx
0=η (y) (y) = 0
dy 2 dy 2

© Matthieu Rigaut 17 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·3 – Écoulement ce Couette plan

Y solution
✧ La solution est immédiate

vx (y) = α y + β
✧ Reste à déterminer α et β. Pour cela, regardons les conditions aux limites.

~v0

~v0

e ~g
~uy
immobiles ~ux
O

✧ Parce que le fluide adhère aux parois, nous devons avoir

v0
vx (0) = 0 et vx (e) = v0 vx (y) =×y
e
✧ Nous pouvons constater que quelques grandeurs pertinentes n’interviennent pas :
➜ µ, ce qui est normal puisqu’il s’agit d’une grandeur inertielle alors même qu’ici, toutes les
trajectoires des particules de fluide sont rectilignes uniformes ;
➜ g qui a une action verticale, i.e. dans le même sens que le gradient de vitesse, ce qui signifie
que c’est ce dernier qui prédomine ;
➜ et surtout η, ce qui prouve que, quelle que soit la viscosité du fluide, à la fin, la situation
cinématique est la même. En revanche, ce n’est pas forcément le cas au niveau des efforts.

Y profil des vitesses

✧ Représentons schématiquement le profil des vitesses dans une section droite du fluide.

~v0

e ~g
~uy
~ux
O

✧ Nous voyons qu’il s’agit là d’un profil affine.

I·3·ii – 2e interprétation du nombre de Reynolds : épaisseur de la couche

limite
✬ deux situations équivalentes
✧ Imaginons la situation où une plaque en dessous d’un fluide infini serait initialement immobile, puis
serait mise en mouvement à la vitesse v0 constante vers la gauche à partir de t > 0.
© Matthieu Rigaut 18 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·3 – Écoulement ce Couette plan

t < 0, immobile t > 0, en mouvement

~v0

✧ Cette situation est, en fait, à peu près équivalente à la situation ci-dessous où du fluide arrive à la
vitesse v0 sur la droite sur une paroi immobile.

L L h
~v0 L t0 =
t0 = v0
v0
~v0

✧ En effet, pour une particule de fluide située à la cote h (cf. schémas), parcourir la distance L au
L
delà du front de la plaque revient à « attendre » la durée t0 = au dessus d’une plaque mise en
v0
mouvement.

✬ où l’équation de diffusion refait une apparition

✧ Notons δ l’épaisseur de la couche limite après la longueur L d’obstacle.

L
~v0

✧ Dans ces conditions, δ est aussi la distance caractéristique sur laquelle s’est diffusée la vitesse v0 dans
la situation équivalente.

fluide encore immobile

L δ fluide en
t0 = mouvement
v0

~v0
© Matthieu Rigaut 19 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·3 – Écoulement ce Couette plan

✧ Il s’agit là d’une équation de diffusion comme nous l’avons montré dans le cas de l’écoulement de
Couette. Nous avons donc, en ordre de grandeur,

∂vx ∂ 2 vx v0 v0
µ =η µ =η 2
∂t ∂y 2 τ δ
L
✧ Et ainsi, puisque d’après la première situation τ =
v0
2 η L 2 η L2
δ ≡ × δ ≡ × L2 2
δ ≡
µ v0 µ L v0 Re
✧ Finalement nous pouvons dire que

L
L’épaisseur de la couche limite est en √ où L est la longueur parcourue sur l’obstacle.
Re

✬ exemple numérique
✧ Considérons une voiture qui roule à environ 100 km/h ou 30 m.s−1 .
✧ Le nombre de Reynolds associé vaut, pour une voiture de 3 mètres de long

µLV 1 × 30 × 3
Re = = −5
∼ 107
η 10
✧ L’écoulement est clairement turbulent.
✧ En revanche, en ce qui concerne l’épaisseur de la couche limite, tant qu’elle ne s’est pas décollée (i.e.
avant l’arrière de la voiture), son épaisseur est au maximum de

L 3
δ=√ =√ ∼ 1 mm
Re 107
✧ Autant dire que cette épaisseur est totalement négligeable géométriquement. Ceci étant, bien que
cette épaisseur soit négligeable, une grande partie des phénomènes physiques se déroulent à l’intérieur,
ce qui interdit de négliger l’effet de cette couche limite. Ainsi, si son épaisseur est négligeable, son
influence ne l’est pas.
✧ Ceci étant, la couche limite n’est pas la seule responsable de toute la physique. Il faut aussi prendre
en compte les turbulences à l’arrière de la voiture.

© Matthieu Rigaut 20 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

I·4 – Écoulement de Poiseuille plan

I·4·i – qu’est-ce qu’un écoulement de Poiseuille ?
✧ Dans l’écoulement de Couette nous avions un fluide mis en mouvement (ou dont le mouvement
était entretenu) par du cisaillement, i.e. par une force tangentielle.
✧ Dans l’écoulement de Poiseuille, le « moteur » du mouvement va être une force normale et dans
le sens du mouvement : un gradient de pression.
✧ Notons que ces forces (les forces pressantes), existent aussi pour les fluides dans lesquels il n’y a pas
de viscosité (ou, du moins, pour lesquels l’effet de viscosité est négligeable).
✧ Dans ces conditions, nous pourrons avoir des écoulements de type « Poiseuille » pour les fluides
non visqueux alors qu’il ne sera pas possible d’avoir, pour ces derniers, des écoulement de type
« Couette ».

I·4·ii – écoulement de Poiseuille plan

✬ dispositif, analyse
✧ Le dispositif ressemble (un peu) à celui de Couette mais possède une différence notable : il est de
taille finie L suivant ~ux de manière à pouvoir pousser le fluide.

P0 + ∆P e P0

O
~uy
~ux

✧ En revanche, nous allons toujours supposer l’invariance par translation suivant ~uz .
✧ Pour simplifier, nous allons aussi supposer :
➜ que la pesanteur est négligeable ;
➜ que l’écoulement est laminaire ;
➜ que le régime est stationnaire.
✧ Dans ces conditions les champs de vitesse et de pression s’écrivent

~v (M,t) = v (x,y) ~ux et P (M,t) = P (x,y)

✧ Les grandeurs pertinentes de ce problème sont :

➜ µ pour l’inertie ;
➜ e et L pour la géométrie ;
➜ η pour les efforts ;
➜ P0 et ∆P pour les contraintes.

© Matthieu Rigaut 21 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

✬ mise en équation
✧ L’équation de Navier – Stokes s’écrit, compte-tenu du fait que nous négligeons la pesanteur et
que l’écoulement est stationnaire
  !
−−→ −−→ ~ ~v
µ ~0 + ~v · grad ~v = −grad P + ~0 + η △

✧ Cherchons d’abord à simplifier un peu cette équation.

✧ Traduisons l’incompressibilité de l’écoulement

∂vx ∂vy ∂vz ∂vx

div ~v = 0 + + =0 =0
∂x ∂y ∂z ∂x
✧ Autrement dit, le champ de vitesse ne dépend pas de x ce qui donne ~v = vx (y) ~ux .
✧ Commençons par projeter l’équation de Navier – Stokes sur ~uy . Cela donne, puisqu’il n’y a de
vitesse que sur ~ux

∂P
0=− +0
∂y
✧ Nous pouvons donc dire que la pression ne dépend pas de y, ce qui fait qu’il reste P (x).
✧ Projetons l’équation de Navier – Stokes sur ~ux

−−→ dP
µ ~v · grad vx = − (x) + △ vx (y)
dx
✧ En simplifiant, cela donne, comme pour l’écoulement de Couette
   
   
 vx   0 
   
   
   
−−→     ∂  d2 vx
 
~v · grad =  0  ·  =0 et △ vx (y) = (y)
∂y  dy 2
   
   
   
 0   0 
   

✧ Il reste donc

dP d2 vx
0=− (x) + η (y)
dx dy 2

✬ résolution
✧ L’équation précédente peut s’écrire

dP d2 vx
(x) = η (y)
dx dy 2
✧ Comme les deux membres de cette équations sont des fonctions de variables indépendantes, nous en
déduisons

dP d2 vx te not
(x) = η (y) = C = K
dx dy 2

© Matthieu Rigaut 22 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

✧ Commençons par résoudre l’équation en pression

dP
(x) = K P (x) = a x + b
dx
✧ Nous pouvons alors trouver très vite la solution en tenant compte des conditions aux limites
 
x
P (0) = P0 + ∆P et P (L) = P0 P (x) = P0 + ∆P 1 −
L
✧ Nous avons alors
dP ∆P ∆P
(x) = − <0 K=−
dx L L
✧ Attaquons-nous, à présent, à l’équation en vitesse

d2 vx ∆P ∆P 2
η (y) = K = − vx (y) = − y + αy + β
dy 2 L 2Lη
✧ Les conditions aux limites donnent
 
(  
 β=0  β=0
v (0) = 0
∆P 2 ∆P
v (e) = 0  −2 η L e + α e = 0
  α = −2 η L e


✧ Finalement
∆P  
v (y) = e y − y2
2ηL

✬ profil des vitesses

✧ Le profil des vitesses à travers une section droite de l’écoulement est parabolique.

P0 + ∆P e P0

~uy
~ux

✬ débit
✧ Comme, ici, la vitesse n’est pas uniforme sur une section de l’écoulement, nous ne pouvons pas écrire
Dv = S v, ce qui n’aurait, d’ailleurs, aucun sens puisque nous ne saurions pas quel v prendre.
✧ Nous allons donc revenir à la définition du débit volumique en se rappelant que la densité surfacique
de courant de volume en volume 1 n’est autre que ~v = ~v . Cela donne
1. De tels noms ne cessent de faire penser à l’auteur que, comme pour la nomenclature en chimie, la rigueur –
nécessaire – peut créer des assemblages très bizarres.

© Matthieu Rigaut 23 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

ZZ ZZ
Dv = ~P
~v (P ) · dS Dv = ~P
~v (P ) · dS
P ∈S P ∈S

e ~uy ~P
dS

~ux ~uz
ℓ

✧ Etant donné que nous prenons une section droite, nous avons
ZZ ZZ
~P = dy dz ~ux

∆P  
dS Dv = ~v (y) ~ux · dy dz ~ux Dv = e y − y 2 dy dz
P ∈S P ∈S 2ηL
✧ Nous avons donc

Z e  Z ℓ
∆P 2 ∆P e3
Dv = × e y − y dy × dz (· · · ) Dv = × ×ℓ >0
2ηL 0 0 2ηL 6

✬ interprétation
✧ Nous voyons que
1 ∆P
Dv ∝
× × e3
η L
✧ C’est ainsi que, contrairement à l’écoulement de Couette la viscosité intervient : plus elle est grande,
plus l’écoulement est lent car « freiné » par les actions internes.
✧ De même nous pouvons constater sans trop de surprise que plus le gradient de pression est élevé,
plus le débit est grand.
✧ Ce n’est qu’en terme de dépendance spatiale qu’intervient la surprise. Le débit est proportionnel à
e3 ce qui est d’interprétation assez complexe surtout qu’il s’agit d’un cas idéal avec une invariance
par translation suivant ~uz . Retenons essentiellement que lorsque l’écoulement est visqueux, diminuer
la section par un facteur 2 diminue le débit de bien plus qu’un facteur 2.

I·4·iii – 3e interprétation du nombre de Reynolds : transport de la quan-

tité de mouvement
✬ deux types de transport
✧ Il y a deux moyens principaux de transporter de la quantité de mouvement.

Le transport convectif de quantité de mouvement correspond à la quantité de

mouvement emportée avec une particule de fluide lors de son mouvement.

Le transport diffusif de quantité de mouvement correspond à la quantité de mouvement

transférée entre particules de fluide par viscosité.

© Matthieu Rigaut 24 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

✬ transport convectif sur l’exemple de Poiseuille plan

✧ Reprenons la situation de l’écoulement de Poiseuille et cherchons la quantité de mouvement trans-
férée par convection.
✧ Comme les particules de fluide avancent globalement suivant +~ux , nous allons chercher la quantité
de mouvement qui passe à travers une surface élémentaire dS~ = dS ~ux pendant la durée dt.

~
dS

~uy
~ux

✧ Entre t et t + dt, des particules de fluide traversent la surface S, chacune transportant la quantité de
mouvement dmPF v.
✧ La quantité de mouvement total transportée vaut donc dmtot v puisque, sur une ligne de courant,
toutes les particules de fluide ont la même vitesse.
✧ Durant dt, la masse qui s’écoule vaut, par définition du débit massique

dmtot = µ v dS dt
✧ Donc, finalement, la quantité de mouvement transportée s’écrit

dpconv = dmtot v = µ v 2 dS dt
✧ Ce que nous pouvons écrire sous la forme

dpconv = jconv dS dt avec jconv = µ v 2

✬ transport diffusif sur l’exemple de Poiseuille plan

✧ Toujours sur l’exemple de l’écoulement de Poiseuille, cherchons la quantité de mouvement trans-
férée par diffusion.
✧ Comme les interactions de viscosité se font entre particules de fluide l’une au dessus de l’autre (i.e.
sur une même verticale), nous pouvons dire que le transfert diffusif est globalement suivant +~uy .
✧ Dans ces conditions, nous allons chercher la quantité de mouvement qui passe à travers une surface
élémentaire dS ~ = dS ~uy pendant la durée dt.

~
dS

~uy
~ux

© Matthieu Rigaut 25 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·4 – Écoulement de Poiseuille plan

✧ Entre t et t + dt, la quantité de mouvement qui « passe » de bas en haut s’écrit, grâce au PFD

d~p = f~bas→haut dt

✧ Ici la force est la force de viscosité, donc nous avons

∂vx
d~pdiff = − η dS dt
∂y
✧ Résultat que nous pouvons écrire sous la forme

∂vx
dpdiff = jdiff dS dt avec jdiff = −η
∂y

✬ comparaison
✧ Calculons le rapport suivant, en ordre de grandeur

jconv µ v2 jconv µV2 µV L

= ≡ ≡ = Re
jdiff ∂vx jdiff V η
−η η
∂y L

Le nombre de Reynolds représente le rapport des termes de transport convectif sur le

transport diffusif de quantité de mouvement.

✬ un terme diffusif étrange

✧ Reprenons le terme de transport diffusif de quantité de mouvement.

∂vx
jdiff = −η
∂y
✧ Nous pouvons remarquer que si nous changeons v −→ −v alors jdiff −→ −jdiff .
✧ Dans le même temps, même si v −→ −v alors jconv = µ v 2 reste identique.
✧ Cela signifie que, dans une situation où le transport diffusif prédomine, si nous changeons de signe
toutes les contraintes, alors le transport diffusif va, lui aussi, changer de signe, i.e. se faire dans l’autre
sens.
✧ Il est donc possible de revenir au point de départ avec un écoulement à très faible nombre de Rey-
nolds comme le montre l’expérience des goutte d’encre étallées puis réassemblées.

Montrer vidéo expérience.

✧ Dans l’expérience, il faut bien comprendre que, dans tous les cas, la transformation est irréversible
car elle obéit à l’équation de diffusion

∂~v ~ ~v
µ = η△
∂t
✧ En revanche, nous venons de prouver que la transformation pouvait être renversable.
© Matthieu Rigaut 26 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

I·5 – écoulement de Poiseuille cylindrique

✧ Dans le paragraphe précédent, nous avons vu l’écoulement dans une canalisation de largeur infinie, ou
du moins très grande devant son épaisseur, ce qui est possible, mais moins fréquent que l’écoulement
dans un tuyau de section circulaire.

I·5·i – champ des vitesses

✬ dispositif, analyse
✧ Prenons un dispositif presqu’analogue au précédent, sauf que cette fois la section est circulaire.

P0 + ∆P R P0 R
O

~uy
~ux

✧ Nous avons toujours :

➜ un écoulement incompressible ;
➜ une longueur finie en L ;
➜ un régime stationnaire ;
➜ un écoulement laminaire ;
➜ une gravité négligée.
✧ En terme de grandeur pertinentes, nous avons les mêmes que précédemment :
➜ µ pour l’inertie dont nous pouvons nous douter qu’elle n’interviendra pas puisque le régime est
stationnaire et que les lignes de courant sont rectilignes ;
➜ R et L pour la géométrie ;
➜ η pour les efforts ;
➜ P0 et ∆P pour les contraintes.
✧ Le soucis, ici, en reprenant l’équation de Navier – Stokes, c’est qu’elle fait intervenir un laplacien
vectoriel en cylindrique. . . que nous ne savons pas exprimer.
✧ Nous avons donc affaire avec une équation juste mais non utilisable. Il va falloir faire autrement.

✬ liminaires
✧ Commençons par simplifier les champs de vitesse et de pression qui s’écrivent, a priori, compte-tenu
du régime permanent stationnaire et de l’invariance par rotation

~v (M) = v (r,x) ~ux et P (M) = P (r,x)

✧ Comme l’écoulement est incompressible, nous avons
∂vx
div ~v = 0 =0 v (M) = v (r) ~ux
∂x
✧ De plus en projetant l’équation de Navier – Stokes sur ~ur nous avons
© Matthieu Rigaut 27 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

∂P
=0 P (x)
∂r

✬ mise en équation par le PFD – version locale

Y le système
✧ Considérons une particule de fluide en version cylindrique, i.e. une particule de fluide comprise entre
r et r + dr, x et x + dx et θ et θ + dθ.

r dθ
derrière devant

dx dr

✧ Les forces qui s’exercent sont :

➜ force à distance, le poids, mais est négligé ;
➜ force de contact normale, la pression ;
➜ force de contact tangentielle, la viscosité.
✧ Comme nous nous intéressons uniquement au mouvement suivant +~ux :
➜ les particules de fluide exerçant des forces pressantes sont celles situées devant (contact en
x + dx) et derrière (contact en x) ;
➜ les particules de fluide exerçant des forces de viscosité sont celles situées sur l’extérieur (contact
en r + dr) et à l’intérieur (contact en r).

Y les forces pressantes

−−→
✧ Nous savons que la résultante des forces pressantes est en −grad P dτ , donc, ici,

dP
df~press · ~ux = − (x) dx dr r dθ
dx

Y les forces de viscosité

✧ Pour les forces de viscosité, il faut faire attention au fait que la surface de contact en r + dr n’est
pas la même que la surface de contact en r.
✧ Dans ces conditions, en faisant (aussi) attention aux signes, nous avons

dfvisc = fext→PF + fint→PF

dvx dvx
= +η (r + dr) S (r + dr) − η (r) S (r)
dr ! dr
d dvx
= +η S (r) dr
dr dr

✧ Et comme la surface de contact S (r) entre deux particules de fluide s’écrit S (r) = r dθ dx, nous
arrivons à
© Matthieu Rigaut 28 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

!
d dvx
dfvisc = +η r (r) dθ dx dr
dr dr

Y regroupement
✧ Le PFD sur une particule de fluide s’écrit, sur ~ux ,

µ dτ e
ax = dfpress,x + dfvisc
✧ Comme les particules de fluide ont toutes une trajectoire rectiligne (portée par ~ux ) et uniforme
(indépendante de x et de t), nous pouvons dire

!
dP d dvx
ax = 0
e dfpress,x + dfvisc − (x) dx dr r dθ + η r (r) dθ dx dr
dx dr dr

✧ Ce qui donne, en simplifiant

!
dP η d dvx
(x) = × r (r)
dx r dr dr

Y résolution
✧ Comme nous sommes face à une équation dont les deux membres dépendent de deux variables
indépendantes, nous pouvons écrire
!
dP η d dvx te not
(x) = × r (r) = C = K
dx r dr dr
✧ Commençons par l’équation en pression.
dP
(x) = K
dx
✧ Comme le gradient de pression est uniforme, nous pouvons tout de suite écrire
dP P (L) − P (0) dP ∆P
(x) = (x) = −
dx L dx L
✧ Reste à résoudre l’équation en vx
! !
η d dvx ∆P d dvx ∆P
× r (r) = − r (r) =− ×r
r dr dr L dr dr ηL
✧ En primitivant
dvx ∆P dvx ∆P α
r (r) = − × r2 + α (r) = − ×r+
dr 2ηL dr 2ηL r
✧ Et encore
∆P
vx (r) = − × r 2 + α ln r + β
4ηL
✧ Les deux constantes d’intégration sont α et r0 .
© Matthieu Rigaut 29 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

✧ Comme en r = 0 la vitesse est finie, cela impose α = 0.

✧ De plus en r = R, le fluide adhère à la paroi, ce qui donne

∆P  
vx (R) = 0 vx (r) = × R2 − r 2
4ηL

Y profil des vitesses

✧ Représentons le profil du champ de vitesse sur une section droite.

P0 + ∆P R P0

~uy
~ux

✧ Il s’agit, comme pour le cas précédent, d’un profil parabolique.

✬ mise en équation par le TCI – version globale

Y le système
✧ Considérons, comme système, non pas une particule de fluide, mais tout le fluide contenu dans le
cylindre de longueur L et de rayon r.

P0 + ∆P R P0
r

~uy
~ux

✧ Remarquons que, comme chaque particule de fluide a une trajectoire rectiligne uniforme, elle subit
une résultante de force nulle.
✧ La somme totale des forces subies par toutes les particules de fluide contenues dans le système ne
peut donc qu’être nulle aussi.
© Matthieu Rigaut 30 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

Y les forces pressantes

✧ Au niveau d’une section la pression est uniforme puisqu’elle ne dépend que de x.
✧ Nous avons donc

fpress = +P (0) π r 2 − P (L) π r 2 fpress = +∆P π r 2

Y les forces de viscosité

✧ Le système ne touche que des particules de fluide situées sur son « extérieur » donc nous pouvons
écrire (en faisant attention au signe)

dvx dvx
fvisc = +η S avec S = 2πrL fvisc = +η 2 π r L
dr dr

Y regroupement
✧ Le TCI donne tout de suite

dvx dvx ∆P
0 = fpress + fvisc +∆P π r 2 + η 2 π r L =0 =− r
dr dr 2ηL
✧ Ce qui se primitive facilement (en tenant compte de la condition au limite vx (L) = 0)

∆P  
vx (r) = × R2 − r 2
4ηL

Y morale
✧ Bien évidemment les deux méthodes arrivent au même résultat, le contraire eût été inquiétant.
✧ Avec la deuxième méthode (la globale) nous passons outre le fait que la vitesse ne doit pas diverger
en r = 0. Cela s’explique par le fait que le système contient le fluide en r = 0 et, donc, exclut de fait
toute divergence.

I·5·ii – quelques aspects de l’écoulement de Poiseuille cylindrique

✬ débit volumique
✧ Comme dans le cas du dispositif de Poiseuille plan, étant donné que la vitesse n’est pas uniforme
sur une section, nous devons reprendre la définition du débit pour pouvoir la calculer.
ZZ
Dv = ~P
~v (P ) · dS
P ∈S

~P = r dr dθ ~ux
✧ Ici dS

R
~P
dS

© Matthieu Rigaut 31 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

✧ Nous arrivons donc à

ZZ   Z R   Z 2π
∆P ∆P
Dv = × R2 − r 2 r dr dθ Dv = × 2 2
R − r r dr × dθ
P ∈S 4ηL 4ηL 0 0

✧ Cela donne, après calculs

∆P R4 π ∆P 1
Dv = × × 2π Dv = × × × R4
4ηL 4 8 L η

✬ pourquoi une dépendance en R4 ?

✧ En fait il faut voir R4 comme R2 × R2 .
✧ Le premier R2 vient « naturellement » de la proportionnalité du débit avec la section.
✧ Le deuxième R2 vient de la viscosité dont nous savons que la résultante est en η △~ ~v soit, en ordre
η
de grandeur en 2 puisqu’ici la viscosité diffuse suivant R.
R
✧ Comme la viscosité est au dénominateur, cela signifie que le facteur R2 monte au numérateur.
✧ Pour la situation Poiseille plan, nous avions une dépendance en ℓ × e3 qu’il fallait voir en :
➜ e × ℓ pour la section ;
➜ e2 pour les effets de viscosité.

✬ exemple qualitatif
✧ Considérons un certain nombre N de tubes dont la section totale est constante et vaut S0 .
✧ Comment varie le débit avec N ?
✧ Pour répondre, remplissons un tableau

grandeur section totale section unitaire rayon débit unitaire débit total
1 √
dépendance fonctionnelle S0 Sunit R4 N × R4
N
N =1 S0 S0 R0 R0 4 Dv
4
S0 R R0 Dv
N 6= 1 S0 √0
N N N2 N

✧ Nous voyons donc qu’avec un fluide newtonien, du moins avec un écoulement à faible nombre de
Reynolds, plus le nombre de conduite est multiplié (à section constante), plus le débit est faible.
✧ Cela est dû à la multiplication du nombre de surface immobiles et, donc, à l’augmentation de forces
de frottement.

✬ une nouvelle analogie diffusive

✧ Le débit volumique n’est qu’un flux de volume à travers une section.
✧ C’est un peu comme le flux thermique
T1 − T2
Φ1→2 =
Rth
✧ Ici, le « moteur » de l’écoulement de Poiseuille est la différence de pression. Nous pouvons donc
écrire
∆P 8Lη
Dv = avec Rhyd =
Rhyd π R4
© Matthieu Rigaut 32 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

✧ De là à transformer un réseau hydraulique en circuit électrocinétique, il n’y a qu’un pas. . . .

✧ Remarquons aussi que cette analogie explique pourquoi les diminutions de pression entre le début et
la fin d’une canalisation s’appellent des « pertes de charge » alors que dans le même temps quand
un dipôle est relié directement à un générateur il est souvent dénommé « charge ».

✬ puissance perdue par viscosité

Y par un bilan énergétique global
✧ Faisons un bilan énergétique entre t et t + dt sur l’ensemble du fluide contenu dans le tuyau de
longueur L et de rayon R.

dEc + dEp = δWnc,ext + δWnc,int

✧ Comme le régime est stationnaire, nous avons dEm = 0.

✧ La pesanteur étant négligée, dEp = 0.
✧ En terme de force extérieures non conservatives, nous avons :
➜ les forces de contact exercées par la paroi sur le fluide ;
➜ les forces pressantes.
✧ Remarquons que les forces de contact exercées par la paroi ne travaillent pas car elles s’exercent
sur des particules de fluide immobiles ! C’est la même chose que lorsque nous cherchons à déterminer
l’énergie transférée par la route sur la roue de voiture, comme la force s’exerce en un point immobile,
le transfert énergétique est nul.
✧ En revanche, pour les forces pressantes, nous avons 2

δWpress = +f (0) v dt − f (L) v dt δWpress = +P (0) π R2 v dt − P (L) π R2 v dt

✧ Ce qui donne

δWpress = +∆P π R2 v dt

✧ Nous reconnaissons l’expression du débit volumique Dv , ce qui donne

δWpress = Dv × ∆P dt

✧ Enfin, dans les interactions intérieures non conservatives, nous ne trouvons que les interactions de
viscosité, ce qui donne

δWnc,int = Pvisc dt

✧ En rassemblant, nous obtenons

0 + 0 = 0 + Dv × ∆P dt + Pvisc dt Pvisc = −Dv × ∆P

✧ La puissance fournie par les interactions de viscosité est bien négative, ouf.
✧ Le lecteur curieux pourra réfléchir à l’interprétation de cette puissance dissipée au regard de l’analogie
électrocinétique précédente.

2. Nous simplifions délibéremment l’expression en faisant apparaître une vitesse uniforme v sachant que la démons-
tration complète reviendrait au même résultat mais avec des intermédiaires un peu plus lourds.

© Matthieu Rigaut 33 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) I·5 – écoulement de Poiseuille cylindrique

Y par un bilan énergétique local

✧ L’idée est toute simple : il suffit de sommer la puissance perdue par chaque particule de fluide à cause
de la viscosité.
ZZZ
Ptot = dPPF

✧ Or, nous savons que la résultante des forces de viscosité agissant sur une particule de fluide s’écrit
~ ~v dτ . Nous avons donc
η△
  ZZZ  
~
dPPF = η △ ~v dτ · ~v Ptot = ~
η △ ~v · ~v dτ

✧ Après ? Il suffit de calculer le laplacien vectoriel de la vitesse, de multiplier scalairement par la vitesse
locale et de sommer le tout dans l’espace. . .
✧ Le lecteur avide de tester ses performances calculatoires peut le faire. L’auteur, quant à lui, préfère
utiliser la première méthode.

© Matthieu Rigaut 34 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II – Écoulements parfaits

II – Écoulements parfaits
II·1 – Modèle de l’écoulement parfait
II·1·i – définition

Un écoulement est dit parfait quand tous les phénomènes de diffusion peuvent être
négligés.

✧ Cela implique, entre autre que :

➜ il n’y a pas de diffusion de quantité de mouvement, i.e. il n’y a pas de viscosité ;
➜ il n’y a pas de diffusion thermique, i.e. chaque particule de fluide subit une transformation
adiabatique.

II·1·ii – en pratique
✧ Techniquement, si considérer un écoulement parfait revient à faire η = 0, nous allons surtout étudier
ce qui se passe en dehors de la couche limite.
L
✧ Rappelons que la couche limite est d’épaisseur √ . Négliger la viscosité revient à faire
Re

η −→ 0 Re −→ ∞ δ −→ 0

✧ Nous allons donc négliger l’épaisseur géométrique de la couche limite.

✧ Qualitativement un écoulement de Poiseuille qui deviendrait de plus en plus « parfait » ressem-
blerait à

Re petit Re plus grand

écoulement parfait

II·2 – Équation d’Euler

II·2·i – à partir de l’équation de Navier – Stokes

© Matthieu Rigaut 35 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·2 – Équation d’Euler

✬ énoncé

équation d’Euler
Dans un écoulement parfait, le champ de vitesse obéit aux équations (équivalentes)
suivantes :

d~v −−→
µ = −grad P + f~v,tot
dt
  !
∂~v −−→ −−→
µ + ~v · grad ~v = −grad P + f~v,tot
∂t
!
∂~v −−→ v 2 − →  −−→
µ + grad + rot ~v ∧ ~v = −grad P + f~v,tot
∂t 2

➜ µ est la masse volumique ;

➜ f~v,tot est la densité volumique de force à distance (poids, électromagnétisme, inertie) ;

✧ La plupart du temps les forces à distance extérieures se réduiront au poids, ce qui fait que nous
écrirons directement

  !
∂~v −−→ −−→
µ + ~v · grad ~v = −grad P + µ ~g
∂t
!
∂~v −−→ v 2 − →  −−→
µ + grad + rot ~v ∧ ~v = −grad P + µ ~g
∂t 2

✬ « démonstration »
✧ Il suffit d’écrire l’équation de Navier – Stokes et de faire η = 0.

II·2·ii – lecture
✬ équation non linéaire
 
−−→
✧ Le terme d’accélération convective ~v · grad ~v implique que l’équation d’Euler est non linéaire.
✧ Pour pouvoir le négliger (et rendre l’équation linéaire), il faut
 
∂~v −−→
≫ ~v · grad ~v
∂t
✧ Cela donne, en ordre de grandeur

V V2 L
≫ T ≪
T L V
✧ Il s’agit là d’une situation rare car cela implique un déplacement lent et plus le déplacement est lent,
plus le nombre de Reynolds diminue, ce qui contrarie le caractère parfait de l’écoulement.
© Matthieu Rigaut 36 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

✬ retrouver la RFSF
✧ En statique des fluides, nous avons

~e
a=
d~v ~
=0 ~0 = −−−→
grad P + f~v,tot
−−→
grad P = f~v,tot
dt

✧ Ce qui n’est autre que la relation de la statique des fluides.

✬ tout est désormais possible

✧ Théoriquement, lorsque l’écoulement est parfait, tout est déterminable.
✧ En effet, un fluide est caractérisé par :
➜ un champ de vitesse, qui est un champ vectoriel ;
➜ un champ de pression, champ scalaire ;
➜ un champ de température, champ scalaire ;
➜ un champ de masse volumique, champ scalaire.
✧ Nous avons donc 6 champs scalaires à déterminer. Pour cela il nous faut 6 lois scalaires. Les voici :
➜ l’équation d’Euler qui, en tant qu’équation vectorielle, « apporte » 3 équations scalaires ;
➜ la conservation de la masse, qui est une équation scalaire ;
➜ l’équation d’état du fluide (µ = Cte pour un liquide, P V = n R T pour un gaz parfait) ;
➜ l’équation traduisant le comportement adiabatique du fluide (T = Cte pour un liquide, P V γ =
Cte pour un gaz parfait).
✧ 6 champs scalaires à trouver, 6 lois scalaires à utiliser, le compte est bon, théoriquement ça passe.
✧ Le problème c’est que dire « théoriquement » ça passe est un peu abusé car, dans l’état actuel des
sciences, il n’est pas possible d’affirmer mathématiquement que la donnée de ces 6 lois scalaires :
➜ permettent d’avoir au moins une solution ;
➜ ne permettent d’avoir qu’une seule solution.
✧ Ne parlons pas des conditions aux limites minimales à donner pour pouvoir assurer l’unicité de la
solution, si l’unicité existe.
✧ Il s’agit de problèmes mathématiques très difficiles et encore largement ouverts.

II·3 – Premières applications

✧ L’équation d’Euler étant une équation fondamentalement non linéaire, pour la résoudre, il n’y a
pas véritablement de méthode, il n’y a presque que des cas particuliers.

II·3·i – surface libre d’un tourbillon

✬ dispositif, modèle
Y comme la tornade
✧ Imaginons un tourbillon dans de l’eau, tourbillon provoqué, pourquoi pas, par un agitateur magné-
tique comme sur la photo ci-dessous 3 .

3. Source : [Link]

© Matthieu Rigaut 37 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·3 – Premières applications

✧ La question que nous nous poserons est celle de la forme de la surface du fluide dont tout le monde
a déjà pu constater qu’elle se « creusait ».
✧ L’écoulement sera modélisé par :
➜ un écoulement incompressible et parfait ;
➜ un vecteur tourbillon uniforme Ω ~ =Ω ~ 0 dans le cylindre de rayon a ;
➜ un écoulement irrotationnel en dehors de cette zone.

✧ L’écoulement est considéré stationnaire.

Y rappel des résultats

✧ Nous retrouvons là le modèle de la tornade vu dans le chapitre 4.
✧ Ce modèle peut tout à fait s’appliquer ici puisque la seule hypothèse que nous avions faite alors était
l’incompressibilité de l’écoulement. À aucun moment nous avions utilisé le fait que le fluide était
liquide ou gazeux.
✧ Rappelons les résultats obtenus

zone ~
Ω ~v
r<a ~0
Ω Ω0 r ~uθ

~0 a2
r>a Ω0 ~uθ
r

© Matthieu Rigaut 38 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·3 – Premières applications

✬ idée, plan de bataille

✧ Pour déterminer la surface libre, nous allons faire comme pour le vase tournant vu au chapitre 5, à
savoir :
➜ chercher le champ de pression à l’aide d’une loi (ici l’équation d’Euler) ;
➜ imposer P (surface) = P0 et chercher la relation entre r et zsurface .
✧ Bien sûr, il faudra faire les deux cas r < a et r > a.

✬ pour r < a
Y accélération particulaire
✧ L’équation d’Euler s’écrit, au choix

  !
∂~v −−→ −−→
µ + ~v · grad ~v = −grad P + µ ~g
∂t
!
∂~v −−→ v 2 − →
 −−→
µ + grad + rot ~v ∧ ~v = −grad P + µ ~g
∂t 2

✧ Nous avons ainsi, pour commencer, parce que le régime est stationnaire

∂~v ~
=0
∂t
 
−−→
✧ Ensuite, comme le calcul de la dérivée convective via ~v · grad ~v n’est pas, a priori, aisé en coor-
données cylindro-polaire, nous allons utiliser la deuxième relation.

 
2 2 2
−−→ v −−→ Ω0 r 
grad = grad 
2 2
   
2 2 2 2
∂  Ω0 r  ∂  Ω0 r 
= ~ur + ~uz
∂r 2 ∂z 2

= Ω0 2 r ~ur

✧ Pour le deuxième terme, le calcul du rotationnel de la vitesse ne pose aucun problème puisque nous
connaissons le vecteur tourbillon. Ainsi

−→   
~ ∧ Ω0 r ~uθ
rot ~v ∧ ~v = 2 Ω
= 2 Ω0 ~uz ∧ Ω0 r ~uθ
= −2 Ω0 2 r ~ur

✧ Et en rassemblant

∂~v −−→ v 2 −
→

+ grad + rot ~v ∧ ~v = ~0 + Ω0 2 r ~ur − 2 Ω0 2 r ~ur = − Ω0 2 r ~ur
∂t 2
© Matthieu Rigaut 39 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

Y champ de pression
✧ L’équation d’Euler conduit donc à
−−→
−µ Ω0 2 r ~ur = −grad P − µ g ~uz
✧ Cela donne, en projetant sur ~uz

∂P
− −µg = 0 P (r,z) = −µ g z + f (r)
∂z
✧ La projection sur ~ur donne, quant à elle

∂P
−µ Ω0 2 r = −
∂r
✧ Soit, en introduisant le début de solution trouvé,

µ Ω0 2 2
−µ Ω0 2 r = 0 − f ′ (r) f (r) = + r +κ avec κ = Cte
2
✧ Et en rassemblant

µ Ω0 2 2
P (r,z) = −µ g z + r +κ
2

✬ pour r > a
Y accélération particulaire
✧ Utilisons la même technique.
✧ Nous avons toujours

∂~v ~
=0
∂t
✧ Ensuite

!
−−→ v 2 −−→ a4
grad = grad Ω0 2
2 2 r2
!
∂ a4
= Ω0 2 ~ur + 0 ~uz
∂r 2 r2
Ω0 2 a4
=− ~ur
r3
✧ Pour le deuxième terme, c’est très rapide
    Ω a2  
−→ ~ ∧ 0 ~uθ ~ = ~0 −→
rot ~v ∧ ~v = 2 Ω et Ω rot ~v ∧ ~v = ~0
r
✧ Et en rassemblant

∂~v −−→ v 2 −
→
 Ω0 2 a4
+ grad + rot ~v ∧ ~v = ~0 − ~ur + ~0
∂t 2 r3
© Matthieu Rigaut 40 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

Y champ de pression
✧ L’équation d’Euler donne

Ω0 2 a4 −−→
−µ ~ur = −grad P − µ g ~uz
r 3

✧ Cela donne, en projetant sur ~uz

∂P
− −µg = 0 P (r,z) = −µ g z + g (r)
∂z
✧ La projection sur ~ur donne, quant à elle

a4 ∂P
−µ Ω0 2
=−
r 3 ∂r
✧ Soit, en introduisant le début de solution trouvé,

a4 µ Ω0 2 a4
µ Ω0 2 = 0 + g ′(r) f (r) = − + κ′ avec κ′ = Cte
r3 2 r2
✧ Et en rassemblant

µ Ω0 2 a4
P (r,z) = −µ g z − + κ′
2 r2

✬ détermination complète du champ de pression

Y les conditions aux limites
✧ Il y a tout d’abord continuité de la pression en r = a, ce qui se traduit par

P (a− ,z) = P (a+ ,z)

✧ De plus nous pouvons dire que « loin » du tourbillon, l’écoulement n’est pas pertubé. La surface libre
est alors plane et nous en profiterons pour la prendre comme origine.
✧ Tout cela se traduit par

P (∞,0) = P0

Y résolution
✧ La deuxième condition nous permet de trouver κ′ et cela donne

0 + 0 + κ′ = P0 κ′ = P0
✧ La première condition se traduit par

µ Ω0 2 2 µ Ω0 2 a4 2
−µ g z + a + κ = −µ g z − r + P0 κ = P0 − µ Ω0 2 a2
2 2 a2
✧ Finalement

µ Ω0 2  2  µ Ω0 2 a4
P (r < a,z) = −µ g z + r − 2 a2 + P0 et P (r > a,z) = −µ g z − + P0
2 2 r2
© Matthieu Rigaut 41 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

✬ surface libre
Y expression
✧ Traduisons le fait que

P (r,zsurf ) = P0
✧ Cela donne, après calculs




 Ω0 2  2 
 pour r < a zsurf = r − 2 a2
2g

 Ω0 2 a4 1
 pour r > a zsurf = −

 2 g r2

Y graphiquement
✧ Schématiquement la surface libre ressemble, en coupe, à la forme ci-dessous.

1/r 2 parabole 1/r 2

✬ retour sur l’accélération particulaire

✧ En fait, nous aurions pu trouver « facilement » les accélérations particulaire des particules de fluide
−−→ v 2  
−→
sans passer par les calculs avec grad et rot ~v ∧ ~v .
2
✧ En effet, vu de dessus, chaque particule de fluide a une trajectoire circulaire.

~v

~a

✧ De plus le mouvement est uniforme donc nous pouvons écrire, conformément à la cinématique du
point matériel sur une trajectoire circulaire
© Matthieu Rigaut 42 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

v2
~e
a = « − ~ur »
r
✧ Nous trouvons ainsi, pour r < a
2 2
Ω0 r
v = Ω0 r ~uθ ~e
a=− ~e
a = −Ω0 2 r ~ur
r
✧ Et pour r > a

a2 Ω0 2 a4 1 a4
v = Ω0 ~uθ ~e
a=− 2 × ~e
a = −Ω0 2 3 ~ur
r r r r
✧ Dans les deux cas, c’est bien ce que nous avions trouvé.

II·3·ii – pression et forme des lignes de courant

✬ répartition de pression

La pression est plus faible au niveau de la concavité des lignes de courant.

✧ Pour le voir, imaginons une ligne de courant en régime stationnaire. Il s’agit alors, forcément, d’une
trajectoire particulaire.

f~

✧ Nous pouvons considérer une particule de fluide comme un point matériel qui aurait alors la même
trajectoire.

f~

✧ Dans ces conditions, nous savons que l’accélération, et donc la résultante des forces, est dirigée vers
la concavité de la trajectoire.

~a

✧ Nous en déduisons que la résultante des forces qui s’exercent sur la particule de fluide est aussi dirigée
vers la concavité.
✧ Or les forces qui s’exercent sont les forces pressantes et les force volumiques
© Matthieu Rigaut 43 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·3 – Premières applications

X −−→
f~ = −grad P + f~v,tot

✧ Écrivons la pression sous la forme

not
P = pression hydrostatique + surpression = Phyd + p

✧ La relation de l’hydrostatique nous permet de simplifier l’expression de la résultante des forces

−−→ X −−→
grad Phyd = f~v,tot f~ = −grad p
−−→ −−→
✧ Nous pouvons donc en conclure que −grad p est dirigé vers la concavité de la trajectoire donc grad p
est dirigé vers l’extérieur.

pression élevée
−−→
grad p

pression faible

✧ Cela démontre bien que la surpression est plus faible à l’intérieur d’une ligne de courant.

✬ effet Coanda
✧ Imaginons un obstacle (comme une balle) qui contrarirait les lignes de courant.

pression « élevée »

pression « faible »

pression « élevée »

✧ Dans ces conditions la pression serait plus faible juste au dessus de la balle alors qu’en dessous, où
les particules de fluides sont quasi-immobiles, la pression est plus élevée.
✧ La balle subit donc une résultante de force vers le haut.
✧ Avec une balle légère il est possible de la maintenir en suspension comme le montre la photoa ci-
dessous issue d’une vidéo de l’expérience 4 .

4. Source :

[Link]

© Matthieu Rigaut 44 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·3 – Premières applications

✧ Les deux captures d’écran suivantes, toujours extraites de la même vidéo, montrent que le flux d’air
est intense au-dessus de la balle alors qu’il est quasiment nul en dessous.

✧ C’est aussi cet effet qui est responsable du fait que, contrairement à ce que laisserait le penser
l’intuition, un objet placé en dessous d’un jet d’eau (pas trop fort) n’est pas éjecté mais, au contraire,
est « aspiré » pour rester sous le jet.

II·3·iii – pression dans un jet libre

✬ jet libre

Un jet est un écoulement à lignes de courant parallèles où la vitesse ne dépend pas de la

ligne de courant.

✧ Autrement dit, en notant ~ux la direction du jet, nous avons ~v = v (x,t) ~ux .

✬ résultat utile

Dans un jet, la répartition de pression orthogonalement aux lignes de courant est

hydrostatique.

© Matthieu Rigaut 45 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·4 – Ondes sonores, l’écho

✧ En effet, en projetant l’équation d’Euler sur ~uy (ou ~uz ) qui sont les direction normales à l’écoulement
nous avons

∂P ∂P
0=− + fv,tot,y = fv,tot,y
∂y ∂y
✧ Ce qui est bien l’équation de statique des fluides.

✬ jet libre

Un jet libre est un jet à l’intérieur d’un autre fluide.

✧ Avec le résultat précédent, nous pouvons dire immédiatement

Dans tout jet libre à l’intérieur de l’atmosphère, la pression n’est autre que la pression
atmosphérique.

II·4 – Ondes sonores, l’écho

✧ Les ondes sonores étant un cas particulier d’écoulement de fluide, il paraît normal de pouvoir les
retrouver à l’aide des équations de la mécanique des fluides.

II·4·i – rappel à la loi

✧ Nous disposons de trois lois.
✧ Il y a tout d’abord l’équation d’Euler puisque, pour les ondes sonores, l’écoulement est considéré
parfait
  !
∂~v −−→ −−→
µ + ~v · grad ~v = −grad P + µ ~g
∂t
✧ Il y a ensuite la loi locale de conservation de la masse

dµ
+ µ div ~v = ~0
dt
✧ Il y a, enfin, la loi comportementale du fluide

µ = µ(P )
✧ Cette dernière loi représente la manière dont le fluide est compressible.
✧ Insistons : les ondes sonores ne sont qu’un écoulement particulier, parfait mais, pour une fois, non
incompressible.

II·4·ii – linéarisation
✬ approximation acoustique
✧ Rappelons que l’approximation acoustique consiste à dire que les ondes sonores ne sont qu’une faible
perturbation de l’état de repos, i.e.
© Matthieu Rigaut 46 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·4 – Ondes sonores, l’écho

P = P 0 + p1 avec p1 ≪ P 0
µ = µ0 + µ1 avec µ1 ≪ µ0
~v = ~0 + ~v1 avec k~v1 k ≪ c

✧ Dans notre cas, nous allons aussi supposer la pression uniforme au repos, ce qui revient à négliger
l’influence du poids.

✬ l’accélération convective est négligeable

✧ Dans la dérivée particulaire, il y a deux termes à prendre en compte : l’accélération locale et l’accé-
lération convective.
✧ Comparons-les en ordre de grandeur en sachant que la longueur caractéristique de propagation pour
les ondes est la longueur d’onde λ

∂
1
∂t variation locale λ/T c
≡ T ≡ ≡ ≫1
−−→ 1 variation convective V V
~v · grad V
λ
✧ Nous pouvons donc en conclure que, dans tous les cas, nous pouvons négliger le terme de variation
convective.

✬ l’équation d’Euler
✧ L’équation d’Euler s’écrit

∂~v1 −−→

µ0 + µ1 = −grad P0 + p1
∂t
∂~v1
✧ Comme le gradient de P0 est nul (pression uniforme) et que le terme µ1 est d’ordre 2, nous
∂t
arrivons à

∂~v1 −−→
µ0 = −grad p1
∂t

✬ conservation de la masse
✧ La loi de conservation de la masse, dans sa deuxième version, donne

dµ
+ µ div ~v = 0
dt
✧ Comme la variation convective est négligeable devant la variation locale, la dérivée particulaire de la
masse volumique se réduit à la dérivée locale et ainsi

∂µ ∂


+ µ div ~v = 0 µ0 + µ1 + µ0 + µ1 div ~v1 = 0
∂t ∂t
✧ µ0 est une constante et µ1 div ~v1 est d’ordre 2, donc nous arrivons à

∂µ1
+ µ0 div ~v1 = 0
∂t
© Matthieu Rigaut 47 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·4 – Ondes sonores, l’écho

✬ loi de comportement phénoménologique

✧ Celle-ci ne change pas.
✧ Le coefficient de compressibilité isentropique s’écrit

1 ∂µ
χS =
µ ∂P
✧ La loi consiste à dire χS = Cte .
✧ Nous avons ainsi, en assimilant les petites variation à l’écart à l’état de repos
1 µ1
χS = p1 (µ0 + µ1 ) χS = µ1
µ 0 + µ 1 p1
✧ Et en négligeant le terme d’ordre 2

p1 µ0 ,χS = µ1

II·4·iii – équation de propagation

✬ équations de couplage
✧ Commençons par supprimer la variable µ1 pour faire des équations de couplage entre p1 et ~v1 .
✧ Pour l’équation d’Euler il n’y a rien à faire

∂~v1 −−→
= −grad p1
µ0
∂t
✧ Pour l’équation de conservation de la masse, ça se passe tout seul

∂µ1 ∂

+ µ0 div ~v1 = 0 et µ1 = p1 µ0 ,χS p1 µ0 χS + µ0 div ~v1 = 0
∂t ∂t
✧ Et ainsi, en simplifiant par µ0 ,

∂p1
χS = −div v~1
∂t

✬ équation de propagation en surpression

✧ Il faut « éliminer » ~v1 .  
−−→
✧ Commençons par calculer △ p1 = div grad p1

  !
−−→ ∂~v1 ∂~v1
div grad p1 = div −µ0 = −µ0 div
∂t ∂t
∂

= −µO div ~v1
∂t !
∂ ∂p1
= −µ0 −χs
∂t ∂t

✧ Et finalement nous arrivons bien à

∂ 2 p1 1 ∂ 2 p1 1
△ p1 = µ0 χS △ p1 = avec c= √
∂t2 c2 ∂t2 µ0 χS

© Matthieu Rigaut 48 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

✬ équation de propagation en vitesse particulaire

✧ Comme nous faisons depuis toujours, pour trouver l’équation vérifié, en 3D, par un champ vectoriel,
calculons le rotationnel du rotationnel.
✧ Ici, pour ~v1 , cela donne d’abord, avec l’aide de la définition du laplacien vectoriel

−→ −→  −−→

~ ~v1
rot rot ~v1 = grad div ~v1 − △
!
−−→ ∂p1 ~ ~v1
= grad −χS −△
∂t
!
−−→ ∂p1 ~ ~v1
= −χS grad −△
∂t
 
∂ −−→ ~ ~v1
= −χS grad p1 − △
∂t
!
∂ ∂ v~1 ~ ~v1
= −χS −µ0 −△
∂t ∂t
∂ 2~v1 ~
= +χS µ0 − △ ~v1
∂t2

−→
✧ Reste à calculer rot ~v . Pour cela reprenons l’équation d’Euler linéarisée et calculons son rotationnel
!  
− → ∂~v1 → −−→
− ∂ − →

rot µ0 = −rot grad p1 rot ~v1 = ~0
∂t ∂t
✧ Nous pouvons donc dire que le rotationnel de la vitesse est constant. Sauf qu’à l’état initial, l’état
de repos, celui-ci est nul puisque tout est immobile.
−
→
✧ De là nous en déduisons rot ~v1 = ~0 puis
2
~0 = +χS µ0 ∂ ~v1 − △
~ ~v1 1 ∂ 2~v1 ~ ~v1 1
=△ avec c= √
∂t2 c2 ∂t2 µ0 χS

II·5 – Relations de Bernoulli

II·5·i – écoulement irrotationnel et homogène
✬ énoncé

Dans un écoulement Parfait Homogène Irrotationnel Stationnaire Inccompressible

(PHISIc), nous pouvons écrire
P v2
+gh+ = Cte dans tout le fluide
µ 2

✧ Commençons par insister : la constante est valable dans tout le fluide. Il existe en effet une autre
relation de Bernoulli où tel ne sera pas le cas.
✧ Pour retenir toutes les hypothèses (5), nous pouvons remarquer que, rassemblées, elles forment l’acro-
nyme PHISIc.
© Matthieu Rigaut 49 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

✧ Le terme g z représente l’énergie potentielle de pesanteur massique. Si d’autre énergies potentielles

sont présentes, alors la relation de Bernoulli s’écrirait

P v2
+ ep + = Cte
µ 2

✬ démonstration
✧ Commençons par reprendre l’équation d’Euler, puisque l’écoulement est parfait.
!
∂~v −−→ v 2 − →  −−→
µ + grad + rot ~v ∧ ~v = −grad P + µ ~g
∂t 2
✧ Sachant que l’écoulement est stationnaire, il se simplifie en
!
−−→ v 2   −−→
−
→
µ ~0 + grad + rot ~v ∧ ~v = −grad P + µ ~g
2
✧ Puis l’écoulement est irrotationnel
!
−−→ v 2 ~ −−→
µ grad + 0 = −grad P + µ ~g
2
✧ De plus nous pouvons, en notant ~uz le vecteur unitaire ascendant écrire

−−→ v 2 −−→ −−→

µ grad = −grad P − µ grad (g z)
2
✧ Divisons par la masse volumique

−−→ v 2 1 −−→ −−→

grad = − grad P − grad (g z)
2 µ
✧ Comme l’écoulement est homogène, nous pouvons passer la masse volumique « dans » le gradient

−−→ v 2 −−→ P −−→

grad = −grad − grad (g z)
2 µ
✧ Enfin, en regroupant les gradients
!
−−→ v2 P v2 P
grad + +gz =0 + + g z = Cte
2 µ 2 µ
✧ Ce qui est bien la loi recherchée en renommant z en h pour rendre l’écriture intrinsèque, i.e. indé-
pendant du système de coordonnées.

II·5·ii – écoulement non irrotationnel et homogène

✬ énoncé

Dans un écoulement Parfait Stationnaire Inccompressible et Homogène (PSIcH), nous

pouvons écrire
2
P v
+gh+ = Cte sur une ligne de courant
µ 2

© Matthieu Rigaut 50 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

✧ Bien que la formulation de la relation soit la même que précédemment, la loi est bien différente car,
cette fois, la constante n’est valable que le long d’une ligne de courant !
✧ La plupart du temps, comme nous ne saurons pas forcément si l’écoulement est rotationnel ou non
(ce qui n’est pas toujours évident à « voir » 5 ), nous utiliserons cette dernière forme de la relation de
Bernoulli.
✧ En revanche, il faudra bien faire attention que l’écoulement soit parfait tout le long de la ligne de
courant. En particulier il ne faut pas que la ligne de courant traverse une zone de turbulence.

✬ démonstration
✧ La démonstration est très similaire à la précédente, évidemment.
✧ Commençons par reprendre l’équation d’Euler, puisque l’écoulement est parfait.
!
∂~v −−→ v 2 −
→
 −−→
µ + grad + rot ~v ∧ ~v = −grad P + µ ~g
∂t 2

✧ Sachant que l’écoulement est stationnaire, cela nous donne, en divisant par la masse volumique

~0 + −−→ v 2 −
→  1 −−→
grad + rot ~v ∧ ~v = − grad P + ~g
2 µ

✧ Multiplions scalairement par d~ℓ, déplacement élémentaire le long d’une ligne de courant, tout en
−−→
transformant ~g en −grad (g z)
!    
−−→ v 2 −→  1 −−→ −−→
grad · d~ℓ + ~
rot ~v ∧ ~v · dℓ = − grad P · d~ℓ − grad (g z) · d~ℓ
2 µ

✧ Comme ~v est, par définition, tangent à la ligne de courant, nous avons

−  
→  −→ 
~v //d~ℓ rot ~v ∧ ~v ⊥ d~ℓ rot ~v ∧ ~v · d~ℓ = 0

✧ Comme l’écoulement est homogène, nous pouvons « rentrer » la masse volumique dans le gradient
et ainsi
! !
−−→ v 2 −−→ P −−→
grad · d~ℓ + 0 = − grad · d~ℓ − grad (g z) · d~ℓ
2 µ

✧ Utilisons la relation fondamentale du gradient

! !
v2 P
d = −d − d(g z)
2 µ

✧ Et regroupons
!
v2 P
d + +gz =0
2 µ

5. C’est ainsi que l’écoulements de Couette et de Poiseuille vus dans la partie précédente ne sont pas irrota-
tionnels. Ceci dit, ce n’est pas grave vu qu’ils ne sont pas parfaits non plus, nous n’aurions pas pu utiliser la relation
de Bernoulli.

© Matthieu Rigaut 51 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli
not v2 P
✧ En notant de manière intrinsèque z = h, cette relation signifie que la grandeur + + g h reste
2 µ
~
constante pour tout déplacement envisagé. Et comme nous avons imposé dℓ suivant une ligne de
courant, cela implique que la constante n’est vraie que sur une telle ligne.

✬ démonstration sans l’homogénéité

✧ Ce petit paragraphe est totalement hors-programme, c’est uniquement culturel.
✧ En fait le caractère « homogène » n’est pas nécessaire à la démonstration. Que l’écoulement soit
incompressible suffit.
✧ En effet, l’intérêt de l’homogénéité de l’écoulement est de pouvoir écrire

1 −−→ −−→ P
grad P = grad
µ µ
✧ En fait, nous avons de manière plus générale et en supposant uniquement le caractère incompressible
de l’écoulement
  !
1 −−→ −−→ P
grad P · d~ℓ = grad · d~ℓ
µ µ
!
−−→ P
✧ Pour le prouver, calculons grad · d~ℓ. Les relations sur les opérateurs vectoriels donnent
µ
! ! !
−−→ P 1 −−→ −−→ 1
grad · d~ℓ = grad P · d~ℓ + P grad · d~ℓ
µ µ µ
!  
1 −−→ P −−→
= grad P · dℓ − 2 grad µ · d~ℓ
~
µ µ

✧ Calculons, maintenant, la dérivée particulaire de la masse volume.

dµ ∂µ −−→
= + ~v · grad µ
dt ∂t
✧ Comme l’écoulement est incompressible, les particules de fluide ont une masse volumique constante
et comme l’écoulement est stationnaire, nous obtenons
−−→
0 = 0 + ~v · grad µ
d~ℓ
✧ Comme ~v = avec d~ℓ le long d’une ligne de courant, nous avons donc
dt
−−→
d~ℓ · grad µ = 0
!
−−→ P
✧ En reportant ce résultat dans le développement de grad · d~ℓ, nous obtenons
µ
! !
−−→ P 1 − −→
grad · d~ℓ = grad P · d~ℓ + 0 C.Q.F.D.
µ µ

© Matthieu Rigaut 52 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

✬ remarques
Y où la mécanique du point refait surface
✧ Sans véritablement s’en rendre compte, le fait de multiplier scalairement par d~ℓ le long d’une ligne
de courant tout en étant en régime stationnaire transforme l’équation d’Euler en PFD. Autrement
dit, nous revenons à la vision de la mécanique du point.
✧ En effet, en régime stationnaire les trajectoires et les lignes de courant sont identiques. Donc si nous
nous focalisons sur une ligne de courant, c’est que nous regardons les trajectoires une à une donc, en
fait, chaque particule de fluide de manière isolée.
✧ En remarquant cela, nous comprenons mieux pourquoi le caractère « homogène » qui implique une
relation entre les particules de fluide n’est pas nécessaire pour la démonstration.

Y et la relation de Bernoulli vue avec les bilans énergétiques ?

✧ En fait, lorsque nous avons « démontré » une relation de Bernoulli dans le chapitre précédent,
c’est la version « PSIcH » que nous avons vue.
✧ Dans ce chapitre la démonstration a été faite via l’équation d’Euler, i.e. via, d’une certaine manière,
le PFD, alors que dans le chapitre précédent nous sommes partis de la conservation de l’énergie.

✬ généralisation

Suivant les hypothèses, les relations de Bernoulli peuvent s’écrire sous la forme
P v2
+ ep + = Cte
µ 2
La relation n’est valable que là où les hypothèses l’autorisent avec ep l’énergie potentielle
massique.

✧ L’énergie potentielle massique est le plus souvent l’énergie potentielle de pesanteur ep,pes = g h.
✧ Il n’est pas interdit d’imaginer des problèmes où l’énergie potentielle est celle associée à la force
d’inertie d’entraînement
1
ep,ie = − Ω2 HM 2
2

II·5·iii – idoinotons vérificateurs d’hypothèses

✬ importance du homogène
✧ Considérons un écoulement irrotationnel, incompressible et non homogène comme, par exemple,
l’atmosphère statique isotherme.
✧ Dans une telle atmosphère, la pression et la masse volumique suivent les lois respectives 6

RT
P (z) = P0 e −z/H et µ(z) = µ0 e −z/H avec H=
Mg
✧ Question la relation de Bernouilli est-elle vérifiée dans tout le fluide ?
✧ Nous avons

P v2 P0 e −z/H P v2 P0
+ +gz = +0+gz + +gz = + g z 6= Cte
µ 2 µ0 e −z/H µ 2 µ0
6. Voir chapitre précédent.

© Matthieu Rigaut 53 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

✧ Nous pouvons constater que la relation n’est pas constante dans tout le fluide.
✧ En revanche, comme nous l’avons vu, nous pourrions utiliser l’autre relation de Bernoulli, celle
qui n’impose pas l’homogénéité. Et là nous aurions, sur une ligne de courant

P v2
+ + g z = Cte
µ 2
✧ C’est bien le cas ici sauf que les lignes de courant se réduisent à un point. Autant dire que la relation
est inutile.

✬ importance du irrotationnel
Y modèle de la tornade
✧ Rappelons les résultats obtenus précédemment.
✧ Pour r < a :
➜ vecteur tourbillon : Ω~ = Ω0 ~uz ;
➜ champ de vitesse : ~v = r Ω0 ~uθ ;
1
➜ champ de pression : P = −µ g z + µ Ω0 (r 2 − 2 a2 ) + P0 .
2
✧ Pour r > a :
➜ vecteur tourbillon : Ω~ = ~0 ;
a2
➜ champ de vitesse : ~v = Ω0 ~uθ ;
r
1 a4
➜ champ de pression : P = −µ g z − µ Ω0 2 + P0 .
2 r

Y pour r < a
✧ Nous avons

!
P v2 1 P0
+ +gz = −g z + Ω0 (r 2 − 2 a2 ) + + r 2 Ω0 2 + g z
µ 2 2 µ
P0
= + Ω0 2 (r 2 − a2 )
µ

✧ Nous trouvons donc

P v2
+ + g z = f (r)
µ 2
✧ Autrement dit, nous avons bien une constante pourvu que nous restions à r fixé.
✧ Or, fixer r, revient à rester sur un cercle, i.e. sur une trajectoire particulaire ou encore une ligne de
courant.
✧ La constante n’est bien constante que sur une ligne de courant : c’est cohérent avec le fait qu’ici,
pour r < a, l’écoulement est rotationnel.

Y pour r > a
✧ Nous avons
© Matthieu Rigaut 54 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·5 – Relations de Bernoulli

!
P v2 1 a4 P0 a4
+ +gz = −g z − Ω0 2 + + Ω0 2 +gz
µ 2 2 r µ r2
P0
=
µ

✧ Cette fois, nous voyons bien que la constante est valable dans toute la zone où le rotationnel est nul,
conformément aux hypothèses de la relation de Bernoulli.

II·5·iv – limite de vitesse pour un écoulement incompressible

✧ Rappelons une loi vue prédemment.

Un écoulement ne peut être considéré comme incompressible que si les vitesses

particulaires sont très inférieures à celle du son.

✧ L’exemple de la propagation sonore nous montre tout de suite qu’il s’agit là d’une condition nécessaire
mais non suffisante. En effet :
➜ l’écoulement est compressible (cf. utilisation de χS ) ;
➜ les vitesses particulaires sont faibles (approximation acoustique).
✧ Pour le justifier, supposons un écoulement incompressible, stationnaire et parfait.
✧ Sur une ligne de courant, nous pouvons écrire, en négligeant la pesanteur

P v2
+ = Cte
µ 2
✧ Ainsi, en supposant que 0 < v < vmax , nous pouvons dire que la pression vérifie la double inégalité

vmax 2
Pmin < P < Pmin + µ
2
✧ Connaissant la variation de pression que subit une particule de fluide, cherchons sa variation de masse
volumique

1 ∂µ 1 δµ δµ
χS = χS = × = χS δP
µ ∂P µ0 δP µ0
✧ Or δP = Pmax − Pmin ce qui conduit à

δµ vmax 2
= χS × µ
µ0 2
✧ Nous reconnaissons l’expression de la célérité du son

1 δµ vmax 2
χS × µ = 2 =
c µ0 2 c2
✧ Dans ces conditions, nous pouvons dire que tant que les particules de fluides ont des vitesses faibles
devant la célérité du son, alors l’écoulement peut être considéré comme incompressible.
✧ Nous ne retrouvons là qu’un des aspects de l’approximation acoustique qui implique que la masse
volumique des particules de fluide ne varie que de manière infinitésimale.
✧ Sauf que dans le cas du son, cette variation infinitésimale correspond à un phénomène intéressant.
© Matthieu Rigaut 55 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

II·6 – Applications des relations de Bernoulli

II·6·i – effet Venturi
✬ qualitativement
✧ Qualitativement, l’effet Venturi correspond à la baisse de pression engendrée par une augmentation
locale de vitesse conformément à la relation de Bernoulli où l’effet de pesanteur est négligé

P v2
+ = Cte
µ 2

Y trompe à eau
✧ Ce phénomène est utilisé, par exemple, dans les trompe à eau en chimie (voir photos ci-dessous 7 ) :
➜ le resserement du tube dans lequel s’écoule l’eau provoque une diminution de pression dans
l’eau ;
➜ la pression étant une fonction continue de l’espace, la pression autour du jet diminue ;
➜ cela crée une dépression par rapport à la pression atmosphérique et, de là, une aspiration.

Y soulever une balle

✧ Autre manifestation (plus spectaculaire) de l’effet Venturi : quand un courant d’air « assez fort » ar-
rive par l’embouchure d’un entonnoir retourné au dessus d’une balle de ping-pong, celle-ci à tendance
à se soulever. Voir la capture d’écran ci-dessous d’une vidéo de l’expérience 8 .

✧ La raison est simple :

➜ la présence de la balle obture le passage pour l’air ;
7. Sources :
➜ [Link]
➜ [Link]

8. Source : [Link]

© Matthieu Rigaut 56 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

➜ comme l’écoulement est incompressible 9 , l’air voit sa vitesse augmenter, spécialement au dessus
de la balle, là où il y a le moins de place ;
➜ puisque la vitesse augmente au-dessus de la balle, la pression diminue ;
➜ et comme la pression est plus faible au dessus de la balle qu’en dessous, cette dernière subit
une poussée de bas en haut.

Y les toits des maisons

✧ L’effet Venturi existe de manière domestique, au niveau des toits.

✧ Ceux-ci, par leurs présences diminuent la place disponible pour le vent, donc celui-ci accélère et la
pression baisse :
➜ si le vent est faible, c’est plutôt une bonne chose pour le tirage des cheminées ;
➜ si le vent est fort, l’effet est tel que les tuiles peuvent s’envoler.

Y voitures de course

✧ Dernier exemple : les voitures de course.

✧ En s’arrangeant, aérodynamiquement pour créer un « goulot » sous la voiture, l’air qui s’y engoufre
doit accélérer ce qui fait chuter la pression. L’intérêt est de créer une force verticale sur la voiture de
haut en bas de manière à « plaquer » la voiture au sol.
✧ En effet, plus la voiture appuie sur le sol, plus la réaction normale est grande et plus la réaction
normale est grande, plus la réaction tangentielle peut l’être.
✧ Et comme seule la réaction fait accélerer, nous voyons bien qu’il s’agit là d’un effet recherché.
✧ L’effet est si important que, pour limiter la vitesse, il est contrôlé par les réglements.
K Remarque. Cet effet, dans le cas des voitures, est appelé « effet de sol »

✬ débimètre

✧ Nous allons voir, dans ce paragraphe, un exemple de dispositif basé sur l’effet Venturi et qui permet
de mesurer des débits.

Y expérience

✧ Voici ci-dessous la photo d’une expérience 10 .

9. L’air n’est pas envoyé à la vitesse du son !

10. Source : [Link]

© Matthieu Rigaut 57 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

✧ Comme cela ne se voit pas, la soufflerie crée un courant d’air dans le tube.
✧ Au niveau du resserement, comme l’écoulement est incompressible, la vitesse doit augmenter donc la
pression diminuer.
✧ Cette diminution de pression crée une « aspiration » nettement visible dans les pailles numéroté.
✧ En connaissant les caractéristiques géométriques de l’appareil, il est possible d’en déduire le débit.

Y dispositif étudié
✧ Ici nous n’allons pas étudier un écoulement d’air, mais un écoulement de liquide dans une canalisation
horizontale.

ℓA
~g
ℓB
SA
SB

✧ Comme nous le voyons sur le schéma précédent, il y a deux petit tuyaux greffés sur l’écoulement
principal.
✧ Nous allons montrer que la différence de hauteur permet de mesurer le débit.

Y relation de Bernoulli
✧ Choisisons la ligne de courant qui passe par l’axe et notons A et B les points sur cette ligne à la
verticale des prises de pression.
© Matthieu Rigaut 58 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

ℓA
~g
ℓB
SA
SB
A B

✧ La relation de Bernoulli s’écrit donc, pour cet écoulement PSIcH (Parfait, Stationnaire, Incom-
pressible et Homogène)

PA vA 2 PB vB 2
+ g hA + = + g hB +
µ 2 µ 2
✧ Comme les points A et B sont à la même cote, il reste

PA vA 2 PB vB 2
+ = +
µ 2 µ 2

Y pression en A
✧ Focalisons-nous sur ce qui se passe autour de A.

A2
ℓA ℓ2
A′1 ~g

A1
ℓ1

✧ Entre A et A1 , nous avons affaire à un jet donc la répartition de pression est hydrostatique, ce qui
donne

PA = PA1 + µ g ℓ1
✧ Entre A1 et A′1 , la pression est continue

′
PA1 = PA1
✧ Enfin entre A′1 et A2 , l’eau étant statique, la répartition de pression est hydrostatique donc

′
PA1 = PA 2 + µ g ℓ2
✧ Par continuité de la pression, nous avons
© Matthieu Rigaut 59 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

PA2 = Patm
✧ En rassemblant le tout

PA = Patm + µ g ℓ1 + µ g ℓ2 PA = Patm + µ g ℓA

Y pression en B
✧ Avec un raisonnement identique, nous trouvons

PB = Patm + µ g ℓB

Y expression du débit
✧ En remplaçant les expressions de la pression dans la relation de Bernoulli nous obtenons

vA 2 vB 2

Patm + µ g ℓA + µ = Patm + µ g ℓB + µ vB 2 − vA 2 = 2 g ℓA − ℓB
2 2
✧ Or, l’incompressibilité de l’écoulement nous permet d’écrire

Dv = SA vA et Dv = SB vB
✧ Ce qui donne

D v 2 Dv 2

SA 2 SB 2
2 − 2 = 2 g ℓA − ℓB Dv 2 = 2 g ℓ A − ℓ B × 2
SB SA SA − SB 2
✧ Et ainsi
s

2 g ℓA − ℓB
Dv = S A S B ×
SA 2 − SB 2

II·6·ii – vidange d’un réservoir : formule de Torricelli

✬ expérience
✧ Prenons deux réservoirs de diamètres différents et remplissons-les à la même hauteur.
✧ Installons des petits tubes de vidange à différentes hauteurs sur ces réservoirs et regardons ce qui se
passe 11 .

11. Source : [Link]

© Matthieu Rigaut 60 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

✧ Nous constatons que plus la vidange se fait près de la surface, plus la vitesse est faible et ce, quelle
que soit la taille du réservoir.
✧ Nous allons montrer, dans ce paragraphe, que la vitesse d’éjection de l’eau ne dépend, en effet que
de la hauteur de l’eau au-dessus de l’éjection.

✬ cas quasi-stationnaire
✧ Supposons l’écoulement quasi-stationnaire

Y modélisation
✧ Modélisons le réservoir de la manière suivante.

✧ La section s de l’écoulement est très petite devant la section S du réservoir.

✧ Cherchons la vitesse au niveau de la vidange.

Y vitesse d’éjection
✧ Comme l’écoulement est PSIcH (Parfait, (quasi-)stationnaire, Incompressible et Homogène), nous
pouvons utiliser la relation de Bernoulli sur la ligne de courant dessinée ci-dessous.

S
A

s
B

✧ Ainsi, entre A et B, nous avons

PA vA 2 PB vB 2
+ + g hA = + + g hB
µ 2 µ 2
✧ Nous avons :
➜ PA = Patm par continuité de la pression ;
➜ PB = Patm parce qu’en B il y a un jet libre.
✧ Il reste

vA 2 vB 2
+ g hA = + g hB
2 2
© Matthieu Rigaut 61 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

✧ De plus, l’incompressibilité de l’écoulement impose

vA S = vB s et S≫s vA ≪ vB
✧ Nous pouvons donc simplifier encore la relation
q
p
vB 2
g hA = + g hB vB = 2 g hA − hB vB = 2gH
2

Y interprétation
✧ Nous reconnaissons là la vitesse acquise par un point matériel après une chute libre de hauteur H
lâché sans vitesse initiale.
✧ En fait, cela n’a rien d’étonnant car, dans le réservoir, les particules de fluides ne perdent pas d’énergie
(pas de viscosité).
✧ Il est donc normal que si elles perdent l’énergie potentielle dm g H, cette dernière ne peut se trans-
former qu’en énergie cinétique.

✬ vérification de la quasi-stationnarité
∂~v ~
✧ Lorsque le régime est stationnaire, cela signifie que = 0.
∂t
∂~v
✧ Donc, lorsque le régime est quasi-stationnaire, cela signifie que est négligeable.
∂t
✧ Négligeable
  devant quoi ? Devant l’autre terme d’accélération, à savoir l’accélération convective
−−→
~v · grad ~v.
✧ Comparons donc les deux accélérations.
✧ En ordre de grandeur, nous avons, en notant τ la durée caractéristique de vidange
 
∂~v vB −−→ vB 2
≡ et ~v · grad ~v ≡
∂t τ h
✧ Ce qui donne

∂~v
∂t h
  ≡
−−→ vB τ
~v · grad ~v

✧ Or, précisément, nous pouvons dire, en ordre de grandeur

h
= vA
τ
✧ Ce qui conduit à

∂~v
∂t vA
  ≡
−−→ vB
~v · grad ~v

✧ Nous pouvons donc en conclure que l’approximation de quasi-stationnarité est d’autant mieux vérifiée
que vA ≪ vB ou, ce qui est équivalent, que s ≪ S.
© Matthieu Rigaut 62 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

II·6·iii – tube de Pitot

✬ dispositif
✧ Un tube de Pitot ou « tube Pitot » est un dispositif très simple permettant de déterminer la
vitesse d’un écoulement.
✧ Ce dispositif est très utilisé dans l’aviation pour connaitre la vitesse de l’avion par rapport à l’air,
paramètre extrêmement important.
✧ Ci-dessous 12 , le tube fin au bout du nez de ce mirage 2000 est un tube Pitot.

✬ modélisation
✧ Vu en coupe, un tube Pitot peut être représenté de la manière suivante.

∆h

✧ Nous allons considérer que l’écoulement est PSIcH à savoir Parfait, Stationnaire, Incompressible et
Homogène.
✧ De plus, nous négligerons l’effet de pesanteur.
✧ Introduisons les points remarquables et traçons les lignes de courant utiles.

12. Source : [Link]

/nez/Mirage.2000C_17g.jpg

© Matthieu Rigaut 63 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) II·6 – Applications des relations de Bernoulli

B∞ B
PB ′
A
A∞

PA ′ ∆h

✬ pressions en A et B
✧ Étant donné les hypothèses, nous pouvons écrire la relation de Bernoulli entre A∞ et A, ce qui
donne

P0 v0 2 PA vA 2
+ = +
µair 2 µair 2
✧ Et comme A est un point d’arrêt, vA = 0 ce qui donne

v0 2
PA = P0 + µair
2
✧ Nous pouvons, de plus, écrire la relation de Bernoulli entre B∞ et B, ce qui donne

P0 v0 2 PB vB 2
+ = +
µair 2 µair 2
✧ Or tout est fait pour que le tube soit fin et qu’au niveau de B l’écoulement ne soit pas perturbé.
✧ En d’autre termes, nous avons vB ∼ v0 . Cela nous conduit à

PB = P0

✬ mesure de vitesse
✧ Le gaz contenu dans le tube étant immobile, nous pouvons dire que

PA = PA ′ et PB = PB ′
✧ Ensuite, par continuité de la pression, nous pouvons écrire, dans le fluide

PA = PB + µeau g H
✧ En rassemblant, cela donne

v0 2 v0 2
PA = PB + µair µeau g H = µair
2 2
✧ Nous pouvons ainsi en déduire l’expression de la vitesse
© Matthieu Rigaut 64 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·7 – Régime non stationnaire

s
2 µeau g H
v0 =
µair

II·7 – Régime non stationnaire

II·7·i – idée
✧ Il est possible d’étudier des cas non stationnaires lorsque ceux-ci sont parfaits, incompressibles et
homogènes.
✧ Pour cela, une fois les lignes de courant identifiées, nous emploierons la technique utilisée lors de la
démonstration de la seconde relation de Bernoulli, à savoir la multiplication par d~ℓ porté par une
ligne de courant.
✧ Ensuite,
Z nous sommerons le résultat sur toute une ligne et aviserons, au cas par cas, du devenir de
∂~v ~
· dℓ.
∂t
✧ En revanche les autres termes seront identiques.

II·7·ii – oscillations dans un tube en U

✬ dispositif
✧ Imaginons un tube en U rempli de liquide.

~g

✧ Imaginons qu’à un moment donné l’un des bras soit plus haut que l’autre. Des oscillations vont
naître. Le but va être de déterminer leur période.

✬ équation d’évolution
Y schématisaton
✧ Commençons par représenter ce qui se passe en introduisant des points remarquables.
© Matthieu Rigaut 65 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·7 – Régime non stationnaire

B
zB

zA = z < 0
A
~uz
~g

~u
M

✧ Notons que l’incompressibilité associée à l’uniformité de la section impose que la vitesse est la même
partout.
✧ De plus, vu qu’il n’y a pas de viscosité, il n’y a pas d’effet de bord et, donc, la vitesse est uniforme
sur une section.
✧ Les lignes de courants ressemblent donc aux lignes représentées sur le schéma.

Y travailler l’équation d’Euler

✧ Commençons par écrire l’équation d’Euler
!
∂~v −−→ v 2 −
→  −−→
µ + grad + rot ~v ∧ ~v = −grad P + µ ~g
∂t 2

✧ Multiplions scalairement par d~ℓ où d~ℓ est choisi sur une ligne de courant. Nous obtenons d’abord
 ! 
2
 
∂~v ~ −−→ v −−→
µ  · dℓ + grad ~
· dℓ + 0 = − grad P · d~ℓ + µ ~g · d~ℓ

∂t 2

✧ Et en utilisant la relation du gradient (sans oublier que l’homogénéité de l’écoulement permet de dire
que µ est une constante)
!
∂~v ~ v2

µ · dℓ + d µ = −dP − d µ g z
∂t 2
✧ Réarrangeons les termes
!
∂~v ~ v2
µ · dℓ = −d µ + P + µ g z
∂t 2
✧ Sommons de A à B la partie droite de l’équation
Z B !
v2 vA 2 vB 2
− d µ +P +µgz = µ + PA + µ g zA − µ − PB + µ g zB
A 2 2 2

✧ Et comme PA = PB = Patm et que vA = vB (la vitesse est la même le long d’une ligne de courant), il
reste
© Matthieu Rigaut 66 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·7 – Régime non stationnaire

Z !
B
v2
− d µ +P +µgz = µ g (zA − zB )
A 2

Y manipuler l’accélération locale

✧ Regardons de plus près le terme de gauche.
✧ Comme, en chaque point fixe de l’espace, ~v = v ~u et d~ℓ = dℓ ~u où ~u est un vecteur constant tangent
à la ligne de courant, le produit scalaire s’écrit

∂~v ~ ∂v
· dℓ = dℓ
∂t ∂t
✧ Et comme la norme de la vitesse ne dépend ni de la ligne de courant, ni de la position sur cette ligne,
nous pouvons écrire

∂v dv
=
∂t dt
✧ La somme de A à B de ce terme donne donc, puisque v est constant le long de la ligne
Z B Z B Z B Z B
∂~v ~ dv ∂~v ~ dv
· dℓ = dℓ · dℓ = dℓ
A ∂t A dt A ∂t dt A

✧ En notant L la longueur totale du tube, nous avons donc

Z B
∂~v ~ dv
· dℓ = ×L
A ∂t dt

Y rassemblement
✧ En reprenant les expressions trouvées, cela nous donne d’abord

dv
µ × L = µ g (zA − zB )
dt
✧ De plus, avec les notations choisies, nous avons (attention aux signes)

zA = −zB = z
✧ En ce qui concerne la vitesse, quand la ligne de courant va de A vers B (cf. schéma), nous avons en
même temps ~v dans le sens de ~u (donc v > 0) et un point A qui va vers le bas, donc avec vA < 0.
Ainsi

dzA
vA = = −v
dt
✧ Et ainsi, finalement

d2 z d2 z 2g
−µ × L = µg2z (t) + z (t)
dt2 dt 2 L
r
2g
✧ Il s’agit de l’équation d’un oscillateur harmonique de pulsation ω0 = .
L
© Matthieu Rigaut 67 / 70 Version du 11 mars 2014
 PC⋆ , Fabert (Metz) II·7 – Régime non stationnaire

✬ discussion
✧ Notons bien que les hypothèses sont toutes importantes ici.
✧ Sans le caractère parfait et incompressible, il n’aurait pas été possible d’écrire que la norme de la
vitesse était partout la même.
✧ D’ailleurs, pour que cette dernière approximation soit vraie, il est important que la ligne de courant
la plus courte (celle à l’intérieure) soit quasiment de la même longueur que la plus grande (celle à
l’extérieur).
✧ Cela rajoute donc une dernière hypothèse : le diamètre du tube doit être très inférieur au rayon du
tube.

© Matthieu Rigaut 68 / 70 Version du 11 mars 2014

 PC⋆ , Fabert (Metz) Mécanique n°6 2013 – 2014

Écoulements de fluides

Au niveau du cours
✬ Programme concerné
✧ Programme de 2e année :
➜ I.A.4. Équations dynamiques locales

✬ Les définitions
✧ Sont à savoir :
➜ viscosité dynamique, viscosité cinématique, fluide newtonien, écoulement laminaire ;
➜ équation de Navier – Stokes, d’Euler ;
➜ écoulement de Couette, de Poiseuille ;
➜ effet Venturi, Coanda ;
➜ tube de Pitot, formulle de Torricelli.

✬ Les grandeurs
✧ Connaître les petites relations suivantes ainsi que leurs interprétations :
∂vx (y)
➜ df~t = ±η ~ux avec [η] en Pl ;
∂y
η
➜ ν = avec ν en m2 .s−1 ;
µ
1 ∂µ
➜ χS = en Pa−1 .
µ ∂P S
✧ Connaître les valeurs de :
➜ viscosité dynamique de l’eau ηeau = 10−3 Pl, de l’air 1,8.10−5 Pl, de l’huile 0,1 Pl ;
➜ viscosité cinématique de l’eau νeau = 10−6 m2.s−1 et de l’air νair = 1,4.10−5 m2.s−1 .

✬ Les lois
✧ Sont à connaître :
➜ équation d’Euler ;
➜ expression surfacique des forces de viscosité pour un fluide newtonien ;
➜ connaître les deux formes de la dérivée particulaire ;
➜ connaître les relations de Bernoulli dans le cas d’un écoulement stationnaire homogène et
incompressible.

✬ la phénoménologie
✧ Connaître / savoir :
➜ interpréter le nombre de Reynolds compte tenu du type d’écoulement (en terme de rapport
d’effet convectif et de viscosité, d’épaisseur de couche limite ou de transport de quantité de
mouvement) ;
➜ savoir retrouver l’équation de Navier – Stokes à partir de l’expression surfacique des forces
de viscosité ;
➜ interpréter l’effet renversable d’un écoulement visqueux ;
➜ la définition d’un écoulement parfait ;
© Matthieu Rigaut Fiche de révision
 PC⋆ , Fabert (Metz) Mécanique n°6 2013 – 2014

➜ reconnaitre les gradients de pression à partir de la forme des lignes de courant ;

➜ quelques applications de l’effet Venturi.

Au niveau des savoir-faire

✬ exercices classiques
✧ Savoir refaire / retrouver :
➜ retrouver le champ des vitesses dans le cas d’un écoulement de Couette plan en régime
stationnaire ;
➜ retrouver le champ des vitesses dans le cas d’un écoulement de Poiseuille plan en régime
stationnaire ;
➜ mettre en équation l’écoulement de Poiseuille cylindrique et en déduire le champ des vitesses
en régime stationnaire ;
➜ retrouver l’équation de propagation des ondes sonores à partir des équations locales de la
dynamique des fluides ;
➜ retrouver l’expression du débit pour un débit-mètre à effet Venturi ;
➜ retrouver la formulle de Torricelli ;
➜ retrouver la vitesse d’un écoulement mesuré par un tube de Pitot.

© Matthieu Rigaut Fiche de révision