Cours Mcaniquedes Fluides

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 92

See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/312291584

Introduction à la mécanique des fluides -Cours (Introduction to fluid


mechanics)

Book · December 2016

CITATIONS READS

0 2,473

1 author:

Abdelhakim Benslimane
Université de Béjaïa
17 PUBLICATIONS   23 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Rheologie of complex fluids View project

Multiaxial Thermomechanical Behavior of Charged Elastomers: Experimentation, Modeling & Validation. View project

All content following this page was uploaded by Abdelhakim Benslimane on 18 January 2017.

The user has requested enhancement of the downloaded file.


UNIVERSITE A. MIRA – BEJAIA
FACULTE DE TECHNOLOGIE
DEPARTEMENT GENIE MECANIQUE

INTRODUCTION A LA
MECANIQUE DES FLUIDES
Cours

Abdelhakim BENSLIMANE
Docteur en Mécanique des Fluides
Laboratoire Mécanique, Matériaux et Energétique
Département de Génie Mécanique- Faculté de Technologie
Université A. MIRA de Bejaia
[email protected] ; [email protected]
Ce document polycopié correspond aux notes de cours de mécanique des fluides dispensé aux

étudiants du premier cycle à l‟université de Bejaia. Il constitue une introduction à la mécanique des

fluides et l‟essentiel de ce qu‟un étudiant doit savoir. Il fournit les bases nécessaires à la

compréhension physique et au calcul des phénomènes en mécanique des fluides.

Dans la préparation de ce cours beaucoup d‟ouvrages classiques de la Mécanique des Fluides ont été

consultés librement. Toutefois, dans la rédaction proprement dite, des cours dispensés à l‟ENGEES et

à l‟UdS m‟ont été d‟un grand apport et dont mon passage se lit entre les lignes de ce document. Je

remercie au passage, K. Bekkour, P. François, R. Mosé, J. Vasquez, M. Dufresne, A. Wanko, F.

Lawniczak, V. Schmitt, Q. Araud, I. Bahlouli, ainsi que toute l‟équipe du Département Mécanique du

laboratoire Icube (ex IMFS). Je reste toujours reconnaissant envers tous mes anciens enseignants de

l‟université de Bejaia, qui m‟ont donné l‟élan intellectuel nécessaire pour réussir dans mon parcours.

Afin d‟améliorer la qualité de ce cours nous encourageons toutes les critiques et remarques et nous

vous remercions par avance.

A. BENSLIMANE, Décembre 2016


J‟essaye d‟écrire un texte, un poème
…J‟essaye d‟écrire, écrire pour te revoir
… Mais je ne vois que le noir dans le noir

Karim Bekkour – Les Fantômes de Laura

Et ainsi de suite, et ainsi de suite. Avant d'être oubliés, nous serons changés en kitsch. Le
kitsch, c'est la station de correspondance entre l'être et l'oubli.

Milan Kundera - L'insoutenable légèreté de l'être


Table des matières
I. GENERALITES ............................................................................................................................. 9
I.1. INTRODUCTION.......................................................................................................................... 9

I.2. DEFINITIONS.............................................................................................................................. 9

I.3. PRESSION : .............................................................................................................................. 12

I.4. STATIQUE ET DYNAMIQUE DES FLUIDES : ............................................................................... 13

I.5. FORCES DE VOLUME ET FORCES DE SURFACE : ....................................................................... 13

II. STATIQUE DES FLUIDES .................................................................................................... 15


II.1. INTRODUCTION : ..................................................................................................................... 15

II.2. PRINCIPE FONDAMENTAL DE LA STATIQUE DES FLUIDES : ..................................................... 15

II.3. VARIATION DE LA PRESSION DANS UN FLUIDE INCOMPRESSIBLE : ......................................... 18

II.4. APPLICATION DU PRINCIPE FONDAMENTAL DE LA STATIQUE DES FLUIDES :.......................... 20

II.5. FORCES DE PRESSION .............................................................................................................. 21

III. DYNAMIQUE DES FLUIDES PARFAITS INCOMPRESSIBLES ................................... 31


III.1. NOTIONS GENERALES SUR L‟ECOULEMENT ......................................................................... 31

III.2. EQUATION DE CONTINUITE .................................................................................................. 33

III.3. EQUATION DE BERNOULLI................................................................................................... 37

III.4. APPLICATIONS AUX MESURES DES DEBITS ET DES VITESSES:.............................................. 41

IV. DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLES .......................................... 45


IV.1. INTRODUCTION .................................................................................................................... 45

IV.2. TENSEUR DES CONTRAINTE DANS UN FLUIDE VISQUEUX .................................................... 45

IV.3. PRINCIPE FONDAMENTAL DE LA DYNAMIQUE .................................................................... 48

IV.4. ANALYSE DU MOUVEMENT D‟UN ELEMENT DE VOLUME FLUIDE – DEFORMATIONS .......... 50

IV.5. EQUATION DE NAVIER-STOKES POUR UN FLUIDE NEWTONIEN .......................................... 52

IV.6. APPLICATIONS ..................................................................................................................... 53

V. THEOREMES DE TRANSPORT ET EQUATIONS DE BILAN ....................................... 61


V.1. INTRODUCTION........................................................................................................................ 61
V.2. FORMULATION DES THEOREMES DE TRANSPORT .................................................................... 61

VI. COUCHE LIMITE ................................................................................................................... 67


VI.1. INTRODUCTION .................................................................................................................... 67

VI.2. LES EQUATIONS DE LA COUCHE LIMITE LAMINAIRE............................................................ 68

VI.3. RESOLUTION DES EQUATIONS DE LA COUCHE LIMITE ......................................................... 73

VII. ANALYSE DIMENSIONNELLE ET SIMILITUDE ........................................................ 81


VII.1. INTRODUCTION .................................................................................................................... 81

VII.2. ANALYSE DIMENSIONNELLE................................................................................................ 81

VII.3. SIMILITUDE .......................................................................................................................... 87


CHAPITRE I : GENERALITES 9

I. GENERALITES

I.1. Introduction

D‟une manière générale, la mécanique des fluides étudie les lois physiques régissant l'écoulement des

liquides et des gaz et aide à reconnaître les causes et les effets de ces écoulements afin de déterminer

leurs paramètres caractéristiques comme le champ de pression ou le champ de vitesse en tenant

compte des différentes propriétés du fluide telles que la densité et la viscosité et principalement les

relations existant entre elles dans différentes situations.

La mécanique des fluides, comme science, a de nombreuses applications dans différents domaines

entre autre : l'aéronavale, l'aéronautique, la météorologie, la climatologie ou encore l'océanographie,

etc…

I.2. Définitions

I.2.1. Fluide

Le physicien distingue classiquement trois états de la matière, solide, liquide et gazeux, en regroupant

sous le vocable fluide les gaz et la plupart des liquides. À l‟échelle microscopique, ce qui caractérise

les fluides, c‟est que les molécules ne sont pas bloquées dans leurs orientations relatives; elles ont ce

degré de liberté (de désordre) que n‟ont pas les molécules dans les solides.

Leurs propriétés communes sont qu‟ils n‟ont pas de forme propre, c‟est-à-dire qu‟ils sont dépourvus

de rigidité; les forces nécessaires pour engendrer des déformations par glissement et assez lentes sont

extrêmement petites. Cette distinction entre solides et fluides n‟est pas parfaitement nette, puisqu‟on

trouve des corps comme les gelées, les peintures, les pâtes, certaines solutions concentrées de

polymères, qui manifestent à la fois des comportements de solides (pendant des temps courts) et des

comportements de liquides (pendant des temps longs).

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE I : GENERALITES 10

En mécanique on appelle fluide un milieu considéré comme continu, à l‟échelle macroscopique, sans

rigidité, pouvant subir de grandes déformations sous l‟action de forces qui peuvent être très faibles.

- Les liquides, incompressibles, assez denses qui épousent la forme des récipients dans lesquels

ils sont placés et présentent une surface supérieure plane quand elle est libre entre autre : l‟eau,

l‟alcool, l‟huile, le sang, en sont des exemples (à température ambiante).

- Les gaz et les vapeurs, compressibles et indéfiniment expansibles, qui occupent tout le volume

qui leur est offert. Citons l‟air que nous respirons, la vapeur d‟eau ou le dioxyde de carbone

qui produit les bulles dans les boissons gazeuses.

Tous les fluides ne sont pas purs liquides ou gaz. On rencontre des fluides où deux ou plusieurs phares

peuvent coexister :

- Les dispersions : sont des mélanges de particules très fines. Ce sont souvent des particules

colloïdales telles que les argiles. Les dispersions ne sédimentent pas spontanément.

- Les suspensions : sont des mélanges de particules fins ou grossières. Contrairement aux

dispersions, les suspensions sédimentent plus ou moins rapidement et peuvent être filtrées

mécaniquement.

- Les émulsions : sont des mélanges de fines gouttelettes d‟un liquide dans un autre. Les

mousses sont des cas particuliers d‟émulsion où les gouttelettes sont les bulles de gaz.

Toutes ces substances ont généralement des comportements particuliers dont l‟étude relève d‟un

domaine de la mécanique appelée rhéologie.

I.2.2. Masse volumique :

La masse volumique   est le rapport :

M Masse  kg 
   m 3  I.1
V Volume

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE I : GENERALITES 11

Pour les liquides, le volume est pratiquement insensible aux variations de pression.

 eau  1000 kg
m3
 mercure  13546 kg
m3
 air  1.205 kg
m3

Remarque :

Contrairement aux liquides, les gaz sont fortement compressibles. La variation de masse volumique

dépend de la température et de la pression :   f  p, T  .

I.2.3. Densité :

La densité est le rapport de la masse volumique d‟un corps à la masse volumique d‟un autre corps pris

comme référence. Le corps de référence pour les liquides est de l‟eau pure à 4°C.

 liquide
  I.2
 eau

I.2.4. Compressibilité :

La compressibilité des fluides est la propriété qui traduit la diminution du volume en réponse à un

accroissement de la pression.

En général, pour un fluide compressible la masse volumique dépend de la température et de la

pression :

  f  p, T 

- Les gaz sont des fluides fortement compressibles.

Les fluides incompressibles ont la propriété d‟avoir leur masse volumique indépendante de toute

variation de pression et/ou de température.

ρ  constante

- Les liquides contrairement au gaz sont peu voir non compressible (incompressible).

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE I : GENERALITES 12

I.2.5. Viscosité :

La notion de viscosité est associée à la résistance qu‟oppose tout fluide à sa mise en mouvement

(écoulement).

- Fluide parfait : un fluide dont la viscosité n‟est pas prise en compte lors de l‟étude de son

écoulement (non visqueux).

- Fluide réel : un fluide dont la viscosité, cette fois-ci, est prise en compte (fluide

visqueux).

Le coefficient de viscosité  appelé viscosité dynamique s‟exprime en SI en [Pa.s] = [N.m-2.s]. Cette

unité porte le nom de Poiseuille [Pl]. On trouve aussi le Poise [Po] avec 1 [Po] = 0.1 [Pl].

On définit la viscosité cinématique  par :




m 2
 s 1  I.3

I.3. Pression :

La pression est une notion physique fondamentale, C‟est l‟intensité de la composante normale de la

force qu‟exerce le fluide sur l‟unité de surface. On peut la voir comme une force rapportée à la surface

sur laquelle elle s‟applique.

p
F
Pa  I.4
S

L‟unité internationale de pression est le Pascal : 1 Pa = 1 N/m². En mécanique des fluides on utilise

encore très souvent le bar. Le bar est égal à peu près à la pression atmosphérique moyenne : 1 bar =

105 Pa.

I.3.1. Pression absolue – pression relative :

- La pression absolue est définie par rapport à la pression dans le vide qui correspond à la

pression nulle. On en déduit donc que la pression minimale possible est zéro.

- La pression relative se définie par rapport à une pression de référence que l‟on choisit le

plus souvent égale à la pression atmosphérique.


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE I : GENERALITES 13

Avec :

pabsolue  prelative  patm I.5

Echelle pression Echelle pression


absolue relative

Patm
0

I.4. Statique et dynamique des fluides :

La mécanique des fluides a pour objet, d‟une part, l‟étude des comportements statiques et dynamiques

des fluides, et d‟autre part l‟étude des interactions entre fluides et solides.

- En statique des fluides, le fluide étudié est supposé au repos (sans mouvement). La vitesse en

chaque point est nulle. Les grandeurs intensives : pression P, température T et masse

volumique ρ restent constantes dans le temps et en chaque point mais bien évidemment

peuvent varier d‟un point à un autre.

- En dynamique des fluides, le fluide est en mouvement. La vitesse en chaque point du fluide

est, par conséquent, fonction non seulement des coordonnées de l‟espace mais aussi du temps.

I.5. Forces de volume et forces de surface :

Les forces qui agissent sur un volume fini de fluide sont deux types :

- Les forces de volumes,

- Les forces de surfaces.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE I : GENERALITES 14

I.5.1. Les forces de volume :

Elles se composent des forces suivantes :

- Les forces de pesanteur provenant de la gravité.

- Les forces d‟accélération pure :

Elles proviennent de la variation de la vitesse (V) d‟une masse (m) dans le temps.

V
Faccélération pure  m I.6
t

- Les forces d‟accélération convective :

Elles proviennent de la variation de la vitesse (u, v, w) dans l‟espace.

 u v w 
Faccélération convective  m u  v  w I.7
 x y z 

I.5.2. Les forces de surface :

- Les forces de pression :

La pression p est le rapport entre une force F agissant perpendiculairement à la surface S d‟un

fluide :

F
p
S

- Les forces de frottement de viscosité :

Nous avons vu précédemment qu‟un fluide, dont les particules sont en mouvement relatif, génèrent

des forces de frottement dues à la viscosité. La force de frottement s‟écrit :

 U 
F  S   I.8
 y 

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 15

II. STATIQUE DES FLUIDES

II.1. Introduction :

La statique des fluides ou ce qu‟on appelle communément l‟hydrostatique étudie les conditions

d‟équilibre des liquides au repos. Ce chapitre aborde l‟étude de la répartition de la pression,

notamment en fonction de la distance verticale, ainsi que les forces qui en résultent.

II.2. Principe fondamental de la statique des fluides :


Fz  dz  px, y, z  dz dS
dS

dz

dy
dx


Fz  px, y, z dS
y

Figure II.1. Pression sur un parallélépipède fluide

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 16

Considérons dans le réservoir un fluide au repos. Isolons un volume parallélépipède d‟eau (Figure 1)


avec : dV  dx  dy  dz  . On suppose un champ de force de volume G g x , g y , g z . 
On fera la démonstration pour la direction z et on généralisera aux autres directions x et y.

Suivant la direction z :

 Les forces de volumes sont :

Concernant les forces de volume, il n‟en existe qu‟une seule : la force de pesanteur. Elle s‟écrit de la

manière suivante :
 
Fv  z   g z dV  ez II.1

Les forces d‟inertie n‟existent pas, le fluide est au repos (vitesse nulle).

 Les forces de surface sont :

La pression agit sur la face supérieure et inférieure de l‟élément. Ces forces s‟écrivent de la manière

suivante :
 
Fs  z   px, y, z dS  ez II.2

 
Fs  z  dz  px, y, z  dz dS   ez  II.3

D‟où :
  
Fs  z   Fs  z dz   px, y, z   px, y, z  dz dS  ez  II.4

  
Fs  z   Fs  z dz   px, y, z  dz   px, y, z dS  ez  II.5

  p 
Fs  z   Fs  z  dz   dz  dS  e z  II.6
dz
  p 
Fs  z   Fs  z  dz    dx  dy  dz  e z  II.7
dz

L‟équation de l‟hydrostatique est déterminée en écrivant l‟équilibre de l‟ensemble des forces. La

condition d‟équilibre des forces selon z s‟écrit :

p  
 dx  dy  dz  e z  g z  dx  dy  dz  e z  0 II.8
dz

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 17

D‟où :

p
  g z  0 II.9
dz

De manière analogue on peut écrire l‟équation d‟équilibre selon les autres directions :

 p
 dx  g x  0

 p
  g y  0 II.10
 dy
 p
  g z  0
 dz

Sous forme vectorielle :



ρG   grad p 0 II.11


Si l‟on se place dans le champ gravitationnel on a : G  0,0, g 

 p
 dx  0

 p
 0 II.12
 dy
 p
   g
 dz

C‟est trois relations montrent que :

- La pression est indépendante de x et de y, la pression ne varie pas dans les directions x et y. La

pression est constante dans les plans horizontaux.

- La troisième relation traduit le principe fondamental de l‟hydrostatique, c'est-à-dire que la

pression croît de haut en bas. La pression ne dépend que de z, on écrit :

dp   gdz II.13

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 18

II.3. Variation de la pression dans un fluide incompressible :

Considérons un réservoir rempli d‟un fluide, exposé à la pression atmosphérique. Le fluide est supposé

incompressible   cste  .
z
3 Surface libre

2 h

Plan de référence

Figure II.2 Variation verticale de la pression

 Calcul de la différence de pression entre les points (1) et (2) :

L‟intégration de l‟équation fondamentale de l‟hydrostatique entre les points (1) et deux (2) s‟écrit :

p  ρgz  cste II.14

p1  ρgz1  p2  ρgz2  cste II.15

p2  p1  ρg z1  z 2  II.16

Cette dernière relation signifie que la variation de pression entre deux niveaux est proportionnelle

(linéaire) à la différence de hauteur entre ces deux niveaux, appelée théorème de Pascal.

La variation de la pression est proportionnelle à la différence de hauteur entre les deux niveaux. On

nomme P l‟énergie potentielle par unité de volume avec :

P  p  ρgz  cste II.17

Une autre écriture consiste à diviser l‟équation par g  :

P p
  z  cste II.18
ρg ρg

La quantité P g est la charge piézométrique appelée aussi hauteur piézométrique. Elle a la

dimension d‟une hauteur [m] et elle est constante dans un liquide au repos.
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE II : STATIQUE DES FLUIDES 19

 Calcul de la différence de pression entre les points (1) et (3) :

Le point (3) est situé à la surface libre du réservoir. La pression en ce même point (3) est égale à la

pression atmosphérique, ce qui nous permet d‟écrire :

p3  patm

En écrivant, l‟équation fondamentale de l‟hydrostatique (EFH) par rapport au plan de référence :

p  ρgz  cste II.19

p1  gz1  p3  gz3  cste II.20

Soit :

p1  p3  g z3  z1  II.21

p1  patm  g z3  z1  I.22

Il faut noter que :

 z 3  z1  h

 z1  z 3  h

Au Final :

p1  patm  ρgh II.23

D‟où :

p1  patm  ρgh II.24

La pression p1 au point (1) est mesurée en pression absolue. Dans la pratique, on préfère utiliser des

pressions mesurées en pression relative par rapport à la pression atmosphérique patm.

p1*  p1  patm  ρgh II.25

La pression p1* au point (1) est mesurée en pression relative. La plupart des instruments de mesure

fournissent une pression relative appelée également pression manométrique.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 20

II.4. Application du principe fondamental de la statique des fluides :

II.4.1. Surface libre d'un liquide

La surface libre d'un liquide au repos est plane et horizontale.

p2  p1  gh
1 h
patm  patm  gh  0

h0 2

II.4.2. Surface de séparation de deux liquides non miscibles

La surface de séparation de deux liquides non miscibles au repos est horizontale.

- dans le fluide I : p B  p A  1 gh
ρ2
- dans le fluide II : p B  p A   2 gh A h

D‟où : 1 gh   2 gh , gh1   2   0
B
ρ1
Or : g  0 et 1   2   0 donc h0

II.4.3. Vases communicants contenant plusieurs liquides non miscibles :

Les dénivellations de deux liquides non miscibles dans des vases communicants sont en rapport

inverse de leurs masses volumiques.

p2  p1  1 gh1 et p2  p3   2 gh2 ρ1 h1
3
Or : p1  p3  patm h2
2

D‟où : 1 gh1   2 gh2


ρ2
Si ρ1 = ρ2 Un fluide est à la même hauteur dans un tube en U.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 21

II.4.4. mesure de la pression atmosphérique (Torricelli, ~ 1643)

La pression atmosphérique (au niveau de la mer) vaut :

3
h
P0 = 1.013105 Pa (i.e. 1.013 bars) soit 76 cm de mercure.

P2 – p3 = ρHg.g.h = p1 = p2 = patm 1 2
ρH

g
II.4.5. Transmission des pressions (principe de Pascal)

Une variation de pression en un point d'un fluide incompressible est transmise intégralement en tout

autre point.

Points 1 et 2 fixes, fluide incompressible : p2 – p1 = ρgh = Cte

Si p1 → p1 + dp de l‟invariance de (p2 – p1=Cte) ⇒ p2 → p2 + dp

Exemple de la presse hydraulique :

Surpression p = F1/S1 transmise en tout point

F2 = p.S2 = (F1/S1). S2

Or : S2 >> S1

F2 >> F1

II.4.6. "Paradoxe" hydrostatique

A surface de fond identique (et même hauteur de liquide), la force de pression exercée par un liquide

sur le fond du récipient est indépendante de la forme du récipient.

II.5. Forces de pression

II.5.1. Théorème d’Archimède

Un corps solide complètement immergé dans un fluide subit de la part de celui-ci une poussée

verticale dirigée de bas en haut et égale au poids du fluide déplacé. Cette résultante est appelée

poussée d‟Archimède. Elle s‟applique au centre de masse du fluide déplacé.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 22

 Le raisonnement est toujours valable si le corps n‟est que partiellement immergé. Dans ce cas

la poussée est égale au poids en fluide du volume effectivement immergé.

Soit une surface fermée formant un corps solide de masse volumique ρs, de volume total V et de

volume immergé Vimmergé, se trouve immergée entièrement ou partiellement (Vimmergé ≤ V) dans un

liquide au repos de poids volumique ρ. Les forces verticales qui agissent sur le corps sont :

- Le poids du solide : ρsgV

- La poussée d‟Archimède : ρgVimmergé

II.5.2. Equilibre des corps immergés

Un corps est en équilibre si le poids P et la force d‟Archimède sont égaux, opposés et situés sur la

même ligne verticale. Dans le cas contraire, il en résulte un mouvement.

Appelons G le centre de gravité d‟un corps flottant et A le centre de poussée de ce même corps

immergée complètement ou partiellement. Le point M, situé à l‟intersection de la verticale passant par

le point A et de l‟axe de symétrie du solides, est appelé métacentre et dm est la distance métacentrique.

 G est en dessous de A : C‟est l‟équilibre stable. En inclinant légèrement le corps,

celui-ci est soumis à un couple de redressement qui le fait tourner jusqu‟à ce qu‟il

revienne à sa position initiale.

 G est au-dessus de A : C‟est l‟équilibre instable. En effet, en inclinant dans ce cas

légèrement le corps celui-ci est soumis à un couple déstabilisant qui le fait tourner

jusqu‟à ce qu‟il retrouve son équilibre stable.

II.5.3. Forces de pression sur les parois

Les forces hydrostatiques sur une surface proviennent des forces de pressions du fluide agissant sur

cette surface. Le calcul des forces hydrostatiques sur une surface quelconque immergée dans l‟eau,

consiste à déterminer les trois caractéristiques suivantes :

 L‟intensité de la force qui s‟applique sur la surface dF  p  dS  gh  dS


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE II : STATIQUE DES FLUIDES 23

 La direction de la force

 Et son point d‟application (centre de poussée).

II.5.3.1. Force de pression sur une paroi plane inclinée

On s‟intéresse aux surfaces planes de forme quelconque entièrement immergée dans l‟eau. La

figure suivante représente à gauche la surface immergée et à droite une vue A-A de cette

surface.

(ρ)
A

z dF1
zG dS

dF2

z x
G

Figure II.4 Force de pression sur un élément de paroi

Soit la surface pleine AB inclinée d‟un angle α par rapport à l‟horizontale et immergée dans

un fluide de masse volumique ρ. Le point G (xG, yG) est le centre de gravité de la section.

Etablissons l‟expression de la force F résultante des forces exercées par le fluide sur la surface

AB. Considérons pour cela la force élémentaire dF s‟exerçant sur un élément de surface dS.

dF  dF1  dF2 II.26

dF  gz  patm dS  patmdS II.27

dF  ρgz  dS II.28

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 24

Avec gz la pression relative à l‟ordonnée z.

F  gz  dS   gx sin   dS  g sin   x  dS


S AB S AB S AB
II.29

Avec :  x  dS  x
S AB
G S AB représente le moment statique de la surface SAB par rapport à l‟axe ox.

Et comme : xG sin   z G

Au final l‟expression de F s‟écrit :

F  g  z G S II.30

zG : profondeur du centre de gravité de la paroi immergée,

S : surface de la paroi immergée.

La résultante des forces de pression exercée par un liquide au repos sur une paroi plane

quelconque, est égale au poids d'une colonne de liquide ayant pour base la surface de la paroi

et pour hauteur la profondeur du centre de gravité de la surface immergée.

II.5.3.2. Position du point d‟application de la force F

Le point d‟application de la force résultante des pressions C (xc, yc) est appelé : centre de

poussée. La position de ce point est définie par la position du barycentre des surfaces

élémentaires (dS) pondérées par la pression sur chaque surface, ce qui revient à calculer le

moment équivalent des forces de pression, c‟est-à-dire :

M / O F    xdF  xc F II.31
S

xc F   xdF   gx
2
sin  dS  g sin    x 2 dS II.32
S AB S AB S AB

x dS  I yy représente le moment d‟inertie de la surface AB par rapport à l‟axe oy.


2
Avec :
S AB

g sin xG S AB xC  g sin I yy II.33

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 25

I yy
xC  II.34
xG S AB

L‟utilisation du théorème de Huygens nous permet d‟écrire : I yy  I GG  S AB xG2

I GG
xC   xG II.35
xG S AB

D‟une autre manière :

I 'GG
zC   zG II.36
zG S xAB

Avec : I ' GG le moment d‟inertie de la surface S xAB . S xAB la projection de S AB sur l‟axe oz.

II.5.3.3. Force de pression sur paroi plane horizontale

Considérons une paroi de largeur unitaire et de surface S immergée horizontalement à une

profondeur z. Les pressions (relatives) aux points (1) et (2) sont égales : p = pl = p2 = pG
z=zc=zG

1 C,G 2

S F
Figure II.5. Forces de pression sur une paroi plane horizontale.

Si on calcul l‟intensité de la force résultante sur la paroi horizontale :

F  g  zG S  gV II.37

Remarque :

o La force F est évidemment égale au poids de la colonne d‟eau au-dessus de la paroi.

o Le point d‟application, C, est confondu avec le centre de gravité G de la surface.

o Quelle que soit la forme des réservoirs, s‟ils sont remplis du même liquide (), jusqu‟à la

même hauteur z, alors les fonds de même surface, S, sont soumis à la même force de pression.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 26

II.5.3.4. Force de pression sur une paroi à surface gauche

x
A
(ρ) Fx

Fz F

z
Figure II.6 Force de pression hydrostatique sur un élément de surface gauche.

Le calcul de la force résultante F sur une paroi gauche AB consiste à écrire F sous forme de deux

composante suivant les axes x et z, respectivement, sous la forme :


  
F  Fx ex  Fz ez II.38

La projection Fx sur la direction horizontale Ox est égale à la poussée hydrostatique s‟exerçant sur la

projection Sx de la surface SAB sur le plan vertical perpendiculaire à Ox.

Fx  gz G S x II.39

Sx : projection de la surface SAB sur un plan vertical. Où Fx passe par le centre de poussée calculée par

rapport à cette surface Sx.

La projection Fz sur la direction verticale Oz de F est égale en grandeur, au poids d‟une colonne de

liquide verticale limitée vers le bas par la surface S et vers le haut par le plan de la surface libre.

Fz  gV II.40

V : le volume de la colonne verticale de liquide limitée par la surface S et la surface libre. Où Fz passe

par le centre de gravité de cette colonne liquide.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 27

II.5.3.5. Applications

a- Surface plane verticale

Déterminer la résultante des forces des forces hydrostatique pression agissant sur une plaque plane

rectangulaire verticale. Puis calculer son point d‟application (centre de poussée).

Calcul de la force résultante :

F  g  z G S II.41

Avec :

 h
 zG 
 2 II.42

S  h  l

1
F g  l  h 2 II.43
2

Le centre de poussée :

lh 3
I h
z C  z G  GG   12 II.44
zG S 2 h
 hl
2

2
zC  h II.45
3

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 28

b- Vanne rectangulaire

Dans cet exemple, comme pour le précédent, il est à déterminer la résultante des forces des forces

hydrostatique pression agissant sur une vanne rectangulaire verticale. Puis de calculer son point

d‟application (centre de poussée).

Calcul de la force résultante :

F  g  z G S II.46

Avec :

 a
 zG  h 
 2 II.47

S  a  l

1  a
F g   h  a  l II.48
2  2

Le centre de poussée :

l  a3
I a 12
z C  z G  GG  h   II.49
zG S 2  a
h   a l
 2

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE II : STATIQUE DES FLUIDES 29

a a2
zC  h   II.50
2  a
12 h  
 2

c- Paroi cylindrique

Calcul de la force résultante F :

F  Fx2  Fz2 II.51

Pour le calcul de F, on passe par le calcul de :

Fx : résultante des forces de pression s‟appliquant sur la projection de la paroi sur l‟axe vertical, et

Fz : le poids du fluide au-dessus de la surface concernée.

Fx  gz G S x II.52

Sx : projection de la surface cylindrique sur l‟axe vertical

zG : profondeur du centre de gravité de Sx

R2
Fx  g b II.53
2

La composante verticale s‟écrit :

R 2
Fz  gV  g b II.54
4
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE II : STATIQUE DES FLUIDES 30

D‟où F est donnée par :

1 2
F  gR 2 b  II.55
4 16

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 31

III. DYNAMIQUE DES FLUIDES

PARFAITS INCOMPRESSIBLES

III.1. Notions générales sur l’écoulement

III.1.1. Régime permanent

L‟écoulement du fluide est dit permanent ou encore stationnaire si les composantes du vecteur vitesse

sont indépendantes de la variable temps t (les composantes sont constantes dans le temps), en chacun

des points de la conduite. D‟une autre manière, la dérivée des grandeurs physiques par rapport au

temps est nulle  t  0 . Il est dit non-permanent ou encore instationnaire si cette condition n‟est pas

réalisée.

III.1.2. Trajectoire et lignes de courant

Trajectoire : On appelle trajectoire la courbe orientée décrite par une particule au cours de son

mouvement, c‟est-à-dire l‟ensemble de ses positions occupée successivement entre deux instants.

z
y
x

Figure III.1 Trajectoire d‟une particule fluide

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 32

Ligne de courant: C‟est la courbe tangente en tout point de l‟espace au vecteur vitesse, à un instant

donné.

z
y
x
M2
M1 
 V2 t0 
V1 t0 

Figure III.2 Ligne de courant d‟une particule fluide

Un tube de courant est l‟ensemble des lignes de courants s‟appuyant sur un contour fermé.

En écoulement permanent, la trajectoire et la ligne de courant sont confondues.

z
y
Tube de courant
x

Lignes de courant

Figure III.3 Tube de courant

III.1.3. Représentation eulérienne et représentation lagrangienne

 Représentation lagrangienne

Le mouvement de la particule fluide peut être défini par ses coordonnées x, y, z, t, tout le long de la

trajectoire - représentation de Lagrange. En d‟autres termes, la méthode revient à suivre une particule

en écoulement.

La méthode de Lagrange s‟avère dans la plupart des cas délicate et peu employée, car il n‟est pas

facile de suivre les particules.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 33

 Représentation eulérienne

Dans ce cas, le mouvement de la particule fluide peut être défini par la vitesse des particules qui

passent en un point donné de l'espace - méthode d'Euler. La description eulérienne consiste donc à

définir les grandeurs physiques en des points fixes.



Nous nous plaçons dans ce qui suit en représentation eulérienne, la vitesse V est déterminée par ses

composantes u, v, w.

Pendant une durée dt, la distance parcourue par la particule s'exprime, en projection sur les axes, par :

dx  udt

dy  vdt III.1
dz  wdt

D‟où :

dx dy dz
dt    III.2
u v w

L‟équation de la ligne de courant est obtenue avec la résolution de l‟équation différentielle.

III.2. Equation de continuité

Soit un tube de courant dans lequel l'écoulement peut être considéré comme unidimensionnel (1D).

 Le fluide est incompressible   cste  .

 L‟écoulement est permanent  t  0 .

Visualisons par la pensée, dS1 et dS2 deux sections transversales droites, situées aux abscisses x1 et x2,

respectivement. z
y

dS1
dS2

V2 dx2 dx1 V1

Figure III.4 Tube de courant

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 34

On notera V1 et V2 les vitesses du fluide sur les sections respectives dS1 et dS2. Durant le temps dt le

fluide qui traverse dS1 et dS2 avance respectivement de dx1  V1dt et dx2  V2 dt .

Dans l‟hypothèse d‟un fluide incompressible, et lorsque le régime est permanent, d‟après le principe

de conservation de masse la masse du fluide qui traverse une section du tube par unité de temps

dm dt  cste  est constante. C‟est-à-dire que la masse dm1 qui traverse la section dS1 pendant un

laps de temps dt est égale à la masse dm2 qui traverse la section dS2 pendant le même laps de temps dt.

dm1 dm2 dm
   cste III.3
dt dt dt

Pendant un temps dt les masses qui traversent, respectivement, dS1 et dS2 :

dm1   dS1  dx1     dS1  V1dt



dm2   dS 2  dx2     dS 2  V2 dt
III.4

On peut avoir :

dm1 dm2 dm
     dS1V1    dS 2V2  dqm III.5
dt dt dt

Avec :

dm  kg 
dq m   s  III.6
dt

dq m est le débit massique qui traverse le tube de courant.

On peut aussi écrire, puisque dm  dV  :

1 dm dV
  dS1  V1  dS 2  V2  dQv III.7
 dt dt

C‟est l‟équation de continuité d‟un écoulement de fluide incompressible.

Avec :

dV  m3 
dQv    III.8
dt  s 

dQv est le débit volumique qui traverse le tube de courant.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 35

Dans la cas de l‟hypothèse du fluide parfait ou fluide non-visqueux   0 , la vitesse est uniforme

sur toute la section de la conduite, donc l‟équation de continuité peut s‟écrire après intégration :

Qv  S1  V1  S 2  V2  cste III.9

Cette relation montre que le produit de la section par la vitesse est constant tout au long d‟une

conduite.

Puisque   dm dV  :

dqm    dQv III.10

III.2.1. Generalisation

Imaginons un parallélépipède de volume élémentaire dV tel que : dV  dx  dy  dz  .

A l‟instant t le parallélépipède contient une masse de fluide égale à :

dmt  dV    dx  dy  dz III.11

Après un laps de temps dt cette masse devient :

  
dmt  dt     dt dx  dy  dz III.12
 t 

D‟une part, la variation de masse durant cet instant est de :


dmt  dt  dmt  dt  dx  dy  dz III.13
t

D‟une autre part, si nous nous intéressons à la variation de la masse fluide traversant les différentes

faces du parallélépipède pendant la durée dt, c‟est-à-dire la différence de masse entre ce qui rentre et

ce qui sort pendant le temps dt. Nous pouvons procéder de la sorte :

Suivant l‟axe x :

On définit le débit massique suivant x:

qmassiquex  ρ  Qvolumiquex  ρu  dy  dz III.14

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 36

La différence de masses fluide entrant par la face de gauche et sortant par la face de droite

pendant l‟intervalle de temps dt est donnée par :

  u    u 
 u  dy  dz   Vx  dx dy  dz  dt   dx  dy  dz  dt III.15
  x   x

Suivant l‟axe y :

  v    v 
 v  dx  dz   Vy  dy  dx  dz  dt   dx  dy  dz  dt III.16
 y  y
   

Suivant l‟axe z :

  w   w
 w  dx  dy   Vz  dz  dx  dy  dt   dx  dy  dz  dt III.17
  z   z

La variation totale de masse de fluide dans l‟élément de volume dV pendant un temps dt est égale à la

masse de fluide qui entre par les différentes faces dans ce volume à laquelle on soustrait la masse de

fluide sortante par les autres faces, on a :

 u  v  w


    dx  dy  dz  dt III.18
 x y z 

D‟où :

  u  v  w


dt  dx  dy  dz       dx  dy  dz  dt III.19
t  x y z 

Donc :

 u  v  w


   0 III.20
t x y z

Sous autre forme :


 

 div V  0 III.21
t

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 37

 Fluide incompressible :

- Si le fluide est incompressible,   cste :


 u v w
div V   
x y z
0 III.22

 Ecoulement permanent :

- Si l‟écoulement est permanent,  t  0 :

 u  v  w


 
div V 
x

y

z
0 III.23

III.3. Equation de Bernoulli

Soit le système fermé contenu à l'instant t entre x1 et x2. A t + Δt entre x1 + l1 et x2 + l2.

 Fluide incompressible, parfait en écoulement permanant

Si on applique le principe de conservation de l‟énergie mécanique :

ΔEc  ΔE p  W III.24

Avec : E c la variation de l‟énergie cinétique, E p la variation de l‟énergie potentielle de pesanteur

et W le travail de pression.
S2 l2= v2 ∆t

v2

l1= v1 ∆t
S1
z2
v1

z1

Figure III.5 Ecoulement dans un tube de courant

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 38

- La variation d‟énergie cinétique :

ΔEc 
1
2

Δm v22  v12  III.25

- La variation de l‟énergie potentielle de pesanteur:

E p  m  g  z 2  z1  III.26

- Le travail de pression :

W  p1 S1v1  t  p2 S 2 v2  t  p1Vol1  p2Vol 2 III.27

Le principe de l‟énergie mécanique s‟écrit :

1
2
 
m v22  v12  m  g  z 2  z1   p1Vol1  p2Vol 2 III.28

En divisant par m :


1 2
2
 p

p
v2  v12  g  z 2  z1   1  2

III.29

v22 p 2 v12 p
 gz 2    gz1  1  cste III.30
2  2 

1 2 1
v1  gz1  p1  v22  gz 2  p2  cste III.31
2 2

1 2
v + gz + p = pt
2
pression pression statique pression totale III.32
dynamique

Ce qui est connu sous l‟appellation de l‟équation de Bernoulli !!!

En divisant par g , il vient :

1 v2 p
 z   H
2 g g
hauteur due côte du hauteur due charge totale III.33
à la vitesse point à la pression
hauteur piézométrique

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 39

III.3.1. Enoncé du théorème de Bernoulli sans échange de travail

En écoulement stationnaire, dans un tube de courant parcouru par un liquide parfait, la charge totale

est constante en tout point.

1 2 1
v1  gz1  p1  v22  gz 2  p2  cste III.34
2 2

L‟équation de Bernoulli est valable en tout point du fluide incompressible en mouvement permanent.

III.3.2. Equation de Bernoulli avec échange de travail

Dans un écoulement, l‟énergie mécanique totale par unité de volume de fluide, peut être modifiée

d‟une section à l‟autre en introduisant dans le circuit une machine hydraulique.

Lorsque le fluide traverse une machine hydraulique, il échange de l‟énergie avec cette machine sous

forme de travail ∆W pendant une durée ∆t.

La puissance P échangée est :

W
P= [J s] ≡ [W ] III.35
t

Si P > 0 : la machine fournie de l‟énergie au fluide (pompe) ;

Si P< 0 : la machine extrait de l‟énergie au fluide (turbine).

Ainsi l‟expression de l‟équation de Bernoulli s‟écrit :

1 2 1 P
ρv1  ρgz1  p1  ρv22  ρgz 2  p 2  III.36
2 2 Qv

III.3.3. Equation de Bernoulli pour les fluides réels

Nous avons vu que pour le cas d‟un fluide réel et en régime permanent, d‟autres forces interviennent,

notamment les forces dues au frottement, qui font apparaître une dissipation de l‟énergie mécanique en

énergie thermique.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 40

Le long d‟une ligne de courant, l‟énergie mécanique diminue du fait de l‟augmentation de l‟énergie

thermique ; on appelle ce phénomène la perte de charge due aux frottements dans un liquide. Après

intégration entre deux points sur la même ligne de courant l‟équation de Bernoulli s‟écrit :

1 2 1
ρv1  ρgz1  p1  ρv22  ρgz 2  p 2  htot III.37
2 2

htot : la perte de charge totale entre le point (1) et (2)

htot  hl  hs III.38

La perte de charges totale n‟est que la somme de la perte de charge linière hl et la perte de charge

singulière hs .

- la perte de charge linéaire hl : représentant l‟énergie perdue entre les deux points,

- la perte de charge singulière hs : intervient lorsque l‟écoulement uniforme est

localement perturbé.

III.3.3.1. Expression de la perte de charge linéaire

Compte tenu des difficultés pour résoudre l‟équation de Navier-Stokes, la perte de charge

sera traduite par une équation empirique du type :

L v2
hl   III.39
D 2g

λ est un coefficient de perte de charge. Il est sans dimension et est fonction du nombre de Reynolds et

de la rugosité de la paroi. L est la longueur de la conduite, D le diamètre, v la vitesse moyenne dans la

conduite.

 En régime laminaire Re < 2000

En régime laminaire, seules les forces de viscosité interviennent. La vitesse est très faible et l‟état de

surface de la paroi n‟intervient pas dans le calcul du coefficient λ

64
 III.40
Re
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 41

Le nombre de Reynolds est donné par :

vD
Re  III.41

 Ecoulement turbulent (Re > 2000).

Le coefficient de perte de charge est donné par la relation de Blasius :

  0.316 Re 0.25 III.42

III.3.3.2. Expression de la perte de charge singulière

La perte de charge singulière, localisée dans une section de la conduite, est provoquée par un

changement de direction et d‟intensité de la vitesse. Une telle non-uniformité de la vitesse peut être

provoquée par :

- un branchement de section de la conduite,

- un changement de direction (coude),

- un branchement ou raccordement,

Comme pour les pertes de charge linéaire, les pertes de charges singulières se traduisent par la

relation :

v2
hs    III.43
2g

 est le coefficient de perte de charge singulière.

III.4. Applications aux mesures des débits et des vitesses:

III.4.1. Le tube de pitot

Le tube de Pitot, représenté ci-dessus, permet de mesurer la vitesse (ou le débit) d‟un fluide dans une

canalisation de section S.

Le fluide circulant dans la canalisation est animé d‟une vitesse v et possède une masse volumique ρ

alors que le tube en U est rempli d‟un liquide de masse volumique ρ0.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 42

- Le fluide est incompressible   cste  et parfait ( = 0) .

- L‟écoulement est permanent (∂∂


t = 0)

B A

B1
h
A1

ρ0

Figure III.6 Tube de Pitot

Nous pouvons donc appliquer la relation de Bernoulli sous la forme:

1 v2 p
+z+ = cste III.44
2 g g

- L‟ouverture A est dirigée face au jet et constitue un point d‟arrêt (la vitesse du fluide en A

est nulle, v A = 0 ).

- z A  zB

- v A1  v A  v B1  0

- La vitesse dans le tube est nulle (pas de mouvement), ce qui nous laisse écrire en

appliquant le principe fondamental de la statique des fluide ans cette partie du système :

 p A1  pB1  ρ0 gh


 pB1  pB  ρg zB  zB1  III.45


 p A1  p A  ρg z A  z A1 
1 v B2 p B 1 v A2 p
 zB    zA  A III.46
2 g g 2 g g

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 43

v B2
pA    pB III.47
2

 
p A1  g z A  z A1  
v B2
2

 p B1  g z B  z B1  III.48

 
v B  2 gh   0  1 III.49
  

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE III : DYNAMIQUES DES FLUIDES PARFAITS INCOMPRESSIBLES 44

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 45

IV. DYNAMIQUE DES FLUIDES REELS

INCOMPRESSIBLES

IV.1. Introduction

Lorsqu‟une particule fluide est soumise à un système de forces non équivalent à zéro, elle se met en

mouvement, et on applique sur elle les lois de la mécanique classique.

Sur un élément de volume, les forces appliquées sont :

- Forces de volume proportionnelles au volume

- Forces de pression proportionnelles aux surfaces et normales à ces surfaces

- Forces d‟inertie proportionnelles à l‟accélération

- Forces de frottement dues à la viscosité

IV.2. Tenseur des contrainte dans un fluide visqueux

Pour un fluide réel (visqueux) en mouvement, les forces de surface ne sont plus seulement normales à

la surface : il existe des contraintes tangentielles dues à la viscosité (frottements).

En un point M d‟une surface dS, la force de surface s‟exprime comme :


 
dF  Tn dS IV.1
 
Où Tn est la contrainte qui s‟exerce sur la surface de normale n .

Considérons maintenant une surface dont l‟orientation est quelconque. Dans le repère cartésien, sa

normale peut se décomposer en :


   
n  nx ex  n y e y  nz ez IV.2

 
      
Tn dS  Tx dS x  Ty dS y  Tz dS z  n xTx  n yTy  n z Tz dS IV.3

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 46

z
ζzz

ηxz ηyz

ηzy ζyy

ηzx ηxy y

ηyx
ζxx

D‟où :
   
Tn  nxTx  n yTy  nz Tz IV.4

 
Considérons une surface perpendiculaire à l‟axe x. La normale à cette surface est : n  e x


La contrainte exercée sur cette surface est alors notée : T x et peut se décomposée comme suite :

   
Tx   xx ex   xy e y   xz e z IV.5

On remarque que la contrainte exercée sur la surface se décompose en une composante normale et

deux composantes tangentielles. On peut de même considérer la surface perpendiculaire à l‟axe y.

On a ainsi la contrainte :
   
Ty   xy e x   yy e y   zy ez IV.6

Et pour la surface ⊥ à l‟axe z la contrainte s‟exprime :


   
Tz   xz ex   yz e y   zz ez IV.7

En développant, on obtient :

n x xx  n y xy  n z xz e x

 
Tn  n x yx  n y yy  n z yz e y

IV.8

 n x zx  n y zy  n z zz e z

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 47

Ce qui revient à effectuer le produit d‟une matrice par la normale :



 Tx     xy  xz   n x 
     xx  
Tn   T y     yx  yy  yz   n y  IV.9
  
 T    zx  zy  zz   n z 
 z
 
Tn  T n IV.10

Avec : T tenseur des contrainte au point M. Ce tenseur des contraintes peut alors se décomposer en

la somme d‟un tenseur sphérique et d‟un tenseur de trace nulle :

T  I  T ' IV.11

 1 0 0    ' xx  xy  xz 
   
T  a 0 1 0     yx  ' yy  ' yz  IV.12
0 0 1    zy  zz 
   zx

Ainsi :

 xx     ' xx

 yy     ' yy IV.13

 zz     ' zz

Avec :  ' xx  ' yy  ' zz  0 et d‟où :  


1
3
 xx   yy   zz   1 Tr T
3

Cette décomposition permet alors de reformuler la contrainte :
   
Tn  T n  I n  T ' n IV.14

Il parait alors deux termes dans l‟expression de la force :


   
dF  Tn dS  T ndS  I ndS  T ' ndS
pression contraintes IV.15
hyrostatique de viscosité

- Le premier terme correspond évidemment à une force purement normale à la surface

qu‟on peut facilement identifier à la force de pression hydrostatique ou (α = -p).

- Le second terme rend compte spécifiquement des forces dues à la viscosité, c‟est le

tenseur des contraintes visqueuse.


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 48

 
dF   pndS  T ' ndS
pression contraintes IV.16
hyrostatique de viscosité

IV.3. Principe Fondamental de la Dynamique

De la même façon que l‟on a établi l‟équation fondamentale de la statique des fluides, nous allons

établir l‟équation fondamentale de la Dynamique en appliquant le principe fondamental de la

dynamique à un élément de volume fluide en mouvement.

Soit un système d‟axes (O, x y, z) et un volume de fluide élémentaire (dV=dx.dy.dz). Il nous faut donc

faire le bilan des forces s‟exerçant en surface et en volume.


 
dv dv
dF  dm  dFs  dFv  dV IV.17
dt dt
 
La dérivée particulaire est donnée par :
dv v   

dt t
 v v  
Avec :

dFv  dVg
    IV.18
dFs  dFsx e x  dFsy e y  dFsz e z

Intéressons-nous à la composante suivant y, chacune des 6 faces est soumise à une contrainte dont une

des trois composantes contribue à dFsy dans la direction de e y . On aura :

     
dFsy   xy xdx   xy x  dydz   yy y dy   yy y  dxdz   yz z dz   yz z  dxdy IV.19

Pour un développement limité à l‟ordre 1 :

  xy
 xy x  dx   xy x  dx
 x
  yy
 yy y  dy   yy y  dy IV.20
 y
  yz
 yz z  dz   yz z  dz
 z

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 49

Ce qui conduit à :

 xy  yy  yz
dFsy  dxdydz  dxdydz  dxdydz IV.21
x y z

  yx  yy  yz 
dFsy     dV
z 
IV.22
 x y

Et par analogie, suivant les autres axes :

   xy  xz 
dFsx   xx   dV IV.23
 x y z 

   zy  zz 
dFsz   zx   dV
z 
IV.24
 x y

Ce qui peut s‟écrire comme :

σ τ xy τ xz 
   xx  
dFs   τ yx σ yy τ yz  dV  T dV IV.25
τ τ zy σ zz 
 zx

Par conséquent :

dv
dF  dFs  dFv  dV IV.26
dt
 
 dv
dF  T dV  gdV  dV IV.27
dt
 
 dv
T  g   IV.28
dt

Au final :
  
dv 
  p  T '  g IV.29
dt

C’est l’équation fondamentale de la dynamique (équation locale).

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 50
 
En posant, g   ge z , la projection de l'équation de la dynamique sur les 3 axes du repère cartésien

donne :

  u u u u  p   ' xx  xy  xz 
    u v  w        
   t x  y z  x   x  y  z 

  v v v v  p   yx  ' yy  yz 
    u  v  w         IV.30
  t  x y z  y  x y z 
 
   w  u w  v w  w w    p    zx  zy   ' zz   g

  t x y z  
z  x y z 

IV.4. Analyse du mouvement d’un élément de volume fluide – Déformations

Au sein de l‟écoulement, chaque particule fluide subit des changements de position, d‟orientation et de

forme. Afin d‟analyser ces changements, considérons 2 points appartenant à la même particule fluide :

M(x,y,z) et M’(x+dx,y+dy,z+dz)

Soient vM (u,v,w) la vitesse au point M, et vM′ (u′,v′,w′) la vitesse au point M’. Exprimons vM′ en fonction

de vM et de dr=MM’

vM′ = v(r+dr)+ dv IV.31

Avec dv l‟accroissement de vitesse.

Alors :

 u u u 
 
     x
u ' u y z   dx 
     v v v   
 v'    v     dy  IV.32
 w'   w   x x x   
    w w w   dz 
 
 x x x 

Avec :

 u u u 
 
 x y z 
  v v v 
gradv   est le tenseur des taux de déformation
x x x 
 w w w 
 
 x x x 

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 51

Or, quel que soit le tenseur, il est toujours possible de le décomposer en deux, pour en faire la somme

d‟un tenseur symétrique et d‟un tenseur antisymétrique :

 u u u   u 1  u v 
  
1  u w  
  
   x 2  y x  2  z x  
 x y z 
  v v v   1  u v  v 1  v w  
gradv          
x x x   2  y x  x 2  z y  
 w w w   1 u w
    1  v w  w 
x   2  z  x     
 x x
   2  z y  x 
IV.33
 1  u v  1  u w  
 0      
 2  y x  2  z x  
 1  u v  1  v w  
      0   
 2  y x  2  z y  
 1  u w  1  v w  
         0 
 2  z x  2  z y  


gradv  D   IV.34

D tenseur des taux de déformations pures (élongation + déformation angulaire).

 tenseur des taux de rotations pures.

On peut en effet remarquer que :

 1  u v  1  u w  
 0      
 2  y x  2  z x  
 0  z y 
 1  u v  1  v w    
       0      z 0  x  IV.35
 2  y x  2  z y   
 1  u w    y x 0 
1  v w 
         0 
 2  z x  2  z y  

Par conséquent :

 0  z  y  dx    x   dx 
      
dr    z 0   x  dy     y    dy  IV.36
 
 y x 0  dz    z   dz 

Finalement :

   
v r  dr   v r   Ddr    dr IV.37

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 52

Ou bien :

 
vM '  vM  DMM '    MM ' IV.38

On dit que l‟écoulement est irrotationnel lorsque les particules fluides ne subissent pas de rotations

pures :   0

IV.5. Equation de Navier-Stokes pour un fluide Newtonien

Par définition, les fluides newtoniens sont ceux pour lesquels les composantes du tenseur des

contraintes de viscosité T ' dépendent linéairement des composantes du tenseur des taux de

déformation pure D .

1  vi v j 

Considérons les éléments tensoriels de D : D  
2  x j xi 

On rappelle que ce tenseur est symétrique.

On admettra alors que pour un fluide isotrope, les éléments tensoriels de T ' et D sont liés par la

relation suivante :

T '  2 D   ' Dxx  D yy  Dzz  ij IV.39

Avec :  la viscosité du fluide,  ' viscosité de dilatation souvent négligée et  ij le symbole de

Kronecker.

Dans le cas d‟un fluide incompressible :

T ' 2 D IV.40

Reprenons l‟équation fondamentale de la dynamique :


  
dv 
  p  T '  g IV.41
dt

Pour un fluide incompressible newtonien, cette équation devient donc :


  
dv 
  p  2 D  g IV.42
dt
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 53

Par conséquent :
 
dv  
  p  2v  g IV.43
dt

C‟est l‟équation fondamentale de la dynamique pour un fluide newtonien incompressible, ou équation

de Navier-Stokes (équation locale).


 
En posant, g   ge z , la projection de l'équation de Navier-Stokes sur les 3 axes du repère cartésien

donne :

  u u u u  p   2u  2u  2u 

    u  v  w 
      2  2  2 
  t  x  y z  x  x y z 
  v v v v  p   2v  2v  2v 

    u  v  w       2  2  2  IV.44
  t x y z  y  x y z 
  w w w w  p  2w 2w 2w 
  u v  w       2  2  2   g
  t x y z  z  x y z 

Les équations de Navier-Stokes et l‟équation de continuité sont les 4 équations simultanées nécessaires

pour résoudre les problèmes de la mécanique des fluides à 4 inconnues, à savoir les composantes de la

vitesse (u, v, w) et la pression p. Néanmoins la résolution de ce système requiert des conditions aux

limites. La condition d‟adhérence à la paroi doit être satisfaite, c‟est à dire que les composantes

normales et tangentielles de la vitesse à la paroi doivent être nulles. Nous allons voir ce que cela

implique dans notre système d‟équations.

IV.6. Applications

IV.6.1. Ecoulement laminaire entre deux plans parallèles

On étudie le cas d‟un fluide incompressible en écoulement permanent et bidimensionnel entre deux

plaques planes et parallèles distantes de 2H. Pour un écoulement unidirectionnel - on a pris ici la

direction selon x - V u, 0, 0 ,

- Ecoulement permanent u t  0 ,

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 54

- Ecoulement uniforme (établi), la vitesse u est indépendante de la direction x

u x  0
y
H
x

-H

Les équations de Navier-Stokes, pour cet écoulement unidirectionnel se résume à:

1 p   ²u 
   
 x  z ² 
1 p
 0 IV.45
 y
1 p
 g 0
 z

On note u = f(y) et p = f(x) car u et p ne dépendent, respectivement, que de y et x.

En intégrant deux fois la première équation on obtient:

u y  
1 dp y ²
 C1 y  C 2 V.46
 dx 2

Avec C1 et C2 les deux constantes d‟intégration, obtenues en utilisant les différentes conditions aux

limites.

Si on pose d‟une manière générale :

y  H  u0
 IV.47
 y  H  u  0

On obtient :

H 2  dp  y 2 
u y    
2  dx 

1  dp  2
y  H 2
    2  1
2  dx  H
IV.48

En tenant compte de la condition aux limites : y  0  u  u max

 y2 
u  y   u max  2  1 IV.49
H 
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 55

Avec :

H2  dp 
u max    IV.50
2  dx 

Le profil de pression est donné par :

 dp  2
   u max 2 IV.51
 dx  H

2u max
p ( x)   x  p0 IV.52
H2

En effet, la pression décroit linéairement avec x.

Le débit volumique est donné par :

H
Qv   u ( y )dS   u y dy
H
IV.53

2  p  3
Qv     H IV.54
3  x 

IV.6.2. Ecoulement la minaire dans une conduite cylindrique

Un écoulement laminaire, incompressible et permanent dans une conduite droite et cylindrique de

section circulaire est appelé écoulement de Hagen- Poiseuille. Cet écoulement peut être considéré

comme étant axisymétrique par rapport au centre de la conduite.

R
x

On va passer en coordonnées cylindriques (x, r, ); la vitesse V v x , vr , u  se réduit encore à

V v x , 0, 0 .
A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 56

Equations de continuité en coordonnées cylindriques est donné par:

  v x  1   rvr  1 v  IV.55


   0
t x r r r 

Equations de Navier-Stokes en coordonnées cylindriques sont données par:

vx v v v vx 1 P   ²vx  ²vx 1 vx 1  2vx 


 vx x  vr x          2 
t x r r   x  x ² r ² r r r  2 

IV.56
vr v v v vr v2 1 P   ²vr  ²vr 1 vr vr 1  2vr 2 v 
 vx r  vr r           2 2  2  
t x r r  r  r  x ² r ² r r r r  2
r  

v v v v v vr v 1 P   ²v  ²v 1 v v 1  2v 2 v 


 vx   vr            2 2  2 r 
t x r r  r   r    x ² r ² r r r r   2
r  

Les équations de Navier-Stokes, pour cet écoulement unidirectionnel se résume à:

 P   ²vx 1 vx 
 x    r ²  r r   0
  
 P IV.57
 0
 r
 1 P
 r   0

On note vx = f(r) et p = f(x) car v et p ne dépendent, respectivement, que de r et x.

La première équation peut être réécrite :

dp 1 d  dv x 
 r  IV.58
dx r dr  dr 

d  dv x  1  dp 
r r   IV.59
dr  dr    dx 

En intégrant, on obtient :

 dv x  1  dp  2
r   r  C1 IV.60
 dr  2  dx 

Ou :

dv x 1  dp  C
  r  1 IV.61
dr 2  dx  r

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 57

En intégrant une deuxième fois :

dv x 1  dp  2
  r  C1 ln r   C 2 IV.62
dr 4  dx 

Les constantes d‟intégration sont déterminées à partir des conditions aux limites :

r  o  r lim 0 v x r   v max v x r  est finie 


 IV.63
r  R  v x R   0

C1  0

 1  dp  2 IV.64
C 2   4  dx  R
  

v x r   
1  dp  2
 r R
4  dx 
2
 IV.65

R 2  dp  r 2 
v x r     2  1 IV.66
4  dx  R 

R 2  dp 
Or par symétrie, v x 0  v max     ceci conduit au profil de vitesse suivant :
4  dx 

 r2 
v x r   v max  2  1 IV.67
R 

Le profil de vitesse est donc parabolique.

Le débit volumique Q traversant la section est défini par :

R
Q   v x r 2rdr  IV.68
0

R 2  p  r 2 
R
Q   2  12rdr IV.69
0
4  x  R 


 p  4
Q   D IV.70
128  x 

Avec D le diamètre de la conduite D=2R.


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 58

La vitesse moyenne Vmoy est donnée par :


 p  4
  D
Q 128  x 
Vmoy   IV.71
S R 2

Q 1  p  2
Vmoy    R IV.72
S 8  x 

1
Vmoy  v max IV.73
2

Le profil de pression est donné par :

 dp  4
   v max 2 IV.74
 dx  R

4v max
p ( x)   x  p0 IV.75
R2

En effet, la pression décroit linéairement avec x.

Pour une longueur de conduite L, on a :

4v max
pL   L  p0 IV.76
R2

p L  p0 4v max
 IV.77
L R2

p p L  p 0 dp
Si    Cst on a
L L0 dx

p
p ( x)   x  p0 IV.78
L

La chute de pression (ou perte de charge) Δp entre la section d‟entrée et la sortie de la conduite

séparées d‟une distance L, est donnée par :

 p  8Vmoy 32Vmoy
    IV.79
 x  R2 D2

La chute de pression dans une conduite est provoquée par les forces de viscosité ce qui signifie que la

perte de charge est régulière ou linéaire, Δhl ; on écrit alors:

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 59

p 32LVmoy
hl   IV.80
g gD 2

2
Vmoy
32LVmoy 2  64L Vmoy
2
64
2
L Vmoy
hl   IV.81
gD 2 2
Vmoy gD 2Vmoy 2 DVmoy D 2 g
2 

 64  L  Vmoy   L  Vmoy 


2 2

hl        IV.82
 Re  D  2 g   D  2 g 

Avec : Re est le nombre de Reynolds pour une conduite circulaire. On appelle λ le coefficient

de frottement ou de perte de charge. Ces résultats sont valables seulement pour un écoulement

laminaire dans conduite circulaire.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE IV : DYNAMIQUE DES FLUIDES REELS INCOMPRESSIBLE 60

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 61

V. THEOREME DE TRANSPORT ET

EQUATIONS DE BILAN

V.1. Introduction

Jusqu‟à présent, les principes de conservation de masse et de quantité de mouvement pour des

systèmes fluides ont été exprimés pour un volume infinitésimal en utilisant une formulation des

équations dite locale. Pour beaucoup d'applications, une résolution complète des équations de bilan

locales n'est pas réalisable.

Dans ce qui va suivre une autre écriture des équations de bilan sera présentée une formulation dite

globale ou intégrale sur un volume de contrôle.

Cette multitude de formulation est au début perçue comme une complexité supplémentaire de la

mécanique des fluides, mais à l‟usage, elle s‟avère fort pratique car cela permet une meilleure

compréhension physique et une résolution plus simple des problèmes.

V.2. Formulation des théorèmes de transport

En Mécanique des Fluides, l‟évolution des grandeurs matérielles est analysée à l‟aide d‟équations

intégrales de bilan sur des domaines fluides macroscopiques. Le transport de ces grandeurs dans

l‟écoulement est explicité en suivant le mouvement; il est par conséquent nécessaire d‟établir

l‟expression de la dérivée particulaire d‟une intégrale volumique.

V.2.1. Volumes et surfaces de contrôle

Un volume de contrôle est un volume imaginaire sur lequel au procède au bilan intégral d‟une

grandeur physique comme la masse, la quantité de mouvement ou encore l‟énergie. On appelle surface

de contrôle l‟enveloppe d‟un volume de contrôle.


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 62

V.2.2. Formulation du cas général



Soit un volume de contrôle V limité par une surface S. Soient f x , t  une fonction scalaire continue et

dérivable et I(t) son intégrale sur le volume V :



I t    f x , t dV V.1
V

Le taux de variation dans le temps de l‟intégrale I(t) est donné par :

 f  
I t    f x,t dV   dV   f v  n dS
d d
V.2
dt dt V V
t S

 
Où v désigne la vitesse locale de la surface de contrôle et n la normale extérieure.

Le second membre du théorème de transport peut être mis sous la forme d‟une intégrale de volume. Il

suffit d‟utiliser la formule de Green-Ostrogradski :

 f 
I t    f x,t dV   dV   div  f  v dV
d d
V.3
dt dt V V
t V

On obtient:

 f 
fdV     div  f  v  dV
d

dt V V 
t 
V.4

  
Le second membre peut être réécrit en utilisant l‟identité : div  f  v   f  div v   v  grad f

df f 
En introduisant aussi la dérivée particulaire :   v  grad f
dt t

On écrira alors :

 df 
fdV     fdiv v  dV
d

dt V V 
dt 
V.5

V.2.3. Cas particuliers


 
Examinons le cas particulier ou f x , t  est une fonction scalaire constante f x , t  =1, dans ce cas

précis l‟intégrale I(t) représente le volume du domaine considéré :

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 63

I t    dV  V t  V.6
V

dI t 
Le taux de variation du volume de contrôle n‟est que la dérivée particulaire soit :
dt

dV t  
  div v dV V.7
dt V

La divergence de la vitesse locale du fluide représente donc le taux de variation relative du volume de

la particule fluide.

On dit qu‟un écoulement de fluide est isovolume ou incompressible si le taux de variation volumique

est nul en tout point et on écrit :



div v   0 V.8
 
Examinons le cas maintenant le cas où f x , t  est la masse volumique  x, t  du fluide, dans ce cas

l‟intégrale I(t) représente la masse du domaine fluide considéré :



I t   M V     x , t dV V.9
V

La masse du domaine est conservée au cours du mouvement ce qui nous laisse écrire :

  
M V      div v  dV  0
d
V.10
dt V 
t 

Dans le cadre de l‟hypothèse de continuité cette expression reste vraie sur un volume infinitésimal. On

en déduit le principe de conservation de la masse sous sa forme locale ou bien ce que nous avons déjà

appelé plus haut équation de continuité :

 
 div v   0 V.11
t

V.2.4. Principe de bilan

Les lois générales de la mécanique des fluides s‟appuient sur le principe commun de l‟analyse de bilan

appelées souvent loi de conservation. Ce principe exprime qu‟une grandeur physique est soit

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 64

invariante (d‟où le vocable loi de conservation), soit varie d‟une manière définie en fonction des

influences extérieures (d‟où le vocable loi de bilan).

Remarquons d‟abord que les grandeurs physiques astreintes à suivre des lois de bilan sont des

grandeurs extensives (masse, quantité de mouvement, énergie, nombre de molécules d‟une espèce

particulière…). Il s‟agit donc de grandeurs matérielles dont on évalue la quantité intégrale sous la

forme d‟une intégrale matérielle:

I t    fdV V.12


V

dont nous savons calculer la dérivée particulaire. Les lois de bilans sont donc des relations entre la

variation dans le mouvement de I(t) et les mécanismes physiques qui génèrent ces variations.

Les mécanismes physiques qui produisent (sources) ou détruisent (puits) la grandeur matérielle

considérée relèvent de deux classes bien distinctes:

 des sources (ou puits), effectivement distribuées à l‟intérieur du domaine Vm mais

commandées par un système extérieur (c‟est le cas des forces de volume comme source de quantité

de mouvement),

 des sources distribuées sur la frontière S du volume de contrôle et commandées par le

matériau situé hors de V (c‟est le cas des forces de surfaces pour la quantité de mouvement).

On peut donc écrire un bilan sous la forme

d  
dt  fdV   dV     ndS
V V S
V.13


où  est la densité volumique des sources internes et  le flux (vecteur ou tenseur) surfacique des

sources externes. En appliquant le théorème de la divergence au dernier terme de la relation et le

théorème de transport, on peut écrire:


  dt dV     div dV
df
V.14
V V

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 65

V.2.5. Equation de bilan de masse



La masse d‟un système matériel étant conservative, les termes sources  et les flux  sont

identiquement nuls. L‟équation de conservation de la masse est obtenue en écrivant les formulations

pour f ≡1. Cette équation a déjà été établie plus haut :

  
  t  div v  dV  0 V.15
V 

V.2.6. Equation de bilan de quantité de mouvement

La quantité de mouvement totale contenue dans le volume de contrôle V est donnée par l‟intégrale

suivante :

I t    v dV V.16
V

Le principe de conservation de la quantité de mouvement n‟est autre que la loi de Newton qui s‟écrit

dans ce cas :

d 
dt  v dV  F
V
ext V.17

où Fext est la résultante des forces extérieures exercées sur le volume V. La force totale exercée sur le

volume V est la somme des contributions des forces à distance (ou volumiques) qui s‟exercent en tout

point du domaine et qui peuvent être représentées par une intégrale de volume et des forces de contact

(ou surfaciques) représentées par une intégrale de surface correspondant au flux du tenseur des

contraintes  à travers la surface S du volume de contrôle.

 
Fext   gdV     ndS V.18
V S

Ce qui permet „établir la forme intégrale de l‟équation de conservation de la quantité de mouvement

comme suit :

d   
dt  v dV   gdV     ndS
V V S
V.17

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE V : THEOREME DE TRANSPORT ET EQUATIONS DE BILAN 66

d   
dt  v dV   gdV     ndS
V V S
V.18

Ce qui peut être écrit :



v   
   div v v dV   gdV   div ( p)dV   div ' dV V.19
V
t V V V V
 

Ou encore :

v      
   v v ndS   gdV   pndS    '  ndS V.20
V
t S V S S

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 67

VI. COUCHE LIMITE

VI.1. Introduction

Dans les écoulements à grande vitesse autour d‟obstacle ou près d‟une paroi, le nombre de Reynolds

de l‟écoulement est le plus souvent très grand, ce qui fait que l‟écoulement peut être considéré à

l‟échelle macroscopique comme étant dans un régime turbulent et les effets de la viscosité sont

négligeables. Toutefois, près d‟une paroi solide, la condition d‟adhérence implique que la vitesse doit

tendre rapidement vers 0.

Si on définit un nombre de Reynolds local à l‟aide de la vitesse réelle (et non d‟une échelle de vitesse),

celui-ci tend également vers 0, ce qui veut dire que très localement, dans le voisinage de la paroi,

l‟écoulement est dans un régime laminaire et les effets de viscosité deviennent prédominants.

Cette zone de faible épaisseur accolée à la paroi s‟appelle une couche limite.

Cette notion a été proposée par Prandtl en 1905 :

- près d‟une paroi solide, il existe une couche de très faible épaisseur dans laquelle les

forces de viscosité sont prédominantes ;

- loin des parois, l‟écoulement peut être considéré comme turbulent ou non visqueux.

Cette décomposition permet de traiter un grand nombre de problème en découplant les effets à grande

échelle (liés à la turbulence) et ceux intervenant à petite échelle près d‟une paroi (et faisant jouer un

rôle crucial à la viscosité du fluide).

Considérons maintenant l‟écoulement d‟un fluide sur une paroi plane. Si ce fluide est parfait, il va

glisser sur la paroi sans que la présence de celle-ci perturbe l‟écoulement (voir Fig. VI.1 ci-dessous).

Cette situation ne correspond pas à l‟expérience physique ; en réalité, on observe la situation suivante :

le fluide réel adhère à la paroi (u = 0 en y = 0) puis sa vitesse augmente progressivement à mesure que

la distance à la paroi augmente pour atteindre finalement, à une distance y = δ de la paroi, une valeur

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 68

dite “vitesse extérieure”, semblable à la vitesse d‟un fluide parfait - de viscosité nulle donc -s‟écoulant

sur cette même paroi.

U
Fluide
parfait

U Couche Fluide
limite visqueux
u=0

Condition de glissement Condition d‟adhérence

Figure VI.1 Ecoulement d‟un fluide parfait et d‟un fluide visqueux sur une plaque plane

VI.2. Les équations de la couche limite laminaire

Dans le cas d‟un fluide newtonien en écoulement bidimensionnel, incompressible et permanent, les

équations de conservation de la masse et de la quantité de mouvement (Equations de Navier- Stokes)

s‟écrivent comme suit :

Equation de continuité :

u v
 0 VI.1
x y

Equations de quantité de mouvement

 u u 1 p   2u  2u 
u  v       
 x y  x  x 2 y 2 
 VI.2
u v  v v   1 p     v   v 
2 2

 x
 y  y  x 2 y 2 

Le fait que l‟épaisseur de la couche limite soit petite devant les autres dimensions caractéristiques de

l‟écoulement permet de faire des approximations dans l‟équation de Navier-Stokes.

On admet aussi que:

- la vitesse de l‟écoulement extérieur est d‟ordre U∞

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 69

- la longueur caractéristique est d‟ordre L

- l‟épaisseur caractéristique est d‟ordre δ (δ : épaisseur de la couche limite), avec δ ˂˂ L

- le gradient de la pression est d‟ordre ρU²∞

 
- v ˂˂ u ainsi que ˂˂
x y

VI.2.1. Normalisation des termes qui figurent dans les équations de Navier- Stokes:

 u x
u*  U x* 
L
 

 v y
v*  y*  VI.3
 V 
 p
 p*  U 2
 

V : une vitesse de normalisation suivant la direction transversale que l‟on déterminera plus tard.

Les dérivées partielles de l‟équation de continuité s‟écrivent :

 u U  u *
 
 x L x *
 v V v * VI.4
 

 y  y *

D‟où l‟équation de continuité :

U  u * V v *
 0 VI.5
L x *  y *

Pour que cette équation soit satisfaite, ses deux termes doivent être du même ordre de grandeur:

U V
 VI.6
L 

On peut donc déduire l‟ordre de grandeur de la vitesse transversale :

U
V  VI.7
L

Les équations de conservation de la quantité de mouvement sont :

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 70

U 2 u * U 2 u * U 2 p * U   2u * U   2u *
 u *  v *      
 L x * L y * L x * L2 x *2  2 y *2
 2 VI.8
U  u * v *  U  v * v *   U  p *   U   v *   U   v *
2 2 2 2


 L x * L y *  y * L3 x *2 L y *2

U 2 U2
En divisant la première équation par et la deuxième par  on aura :
L 

 u * u * p *   2u * L  2u *
u *  v *    
 x * y * x * U  L x *2 U  2 y *2
 2 VI.9
  u * v *  v * v *    p *    v *    v *
2 2 2

 L  x * y *  y * U  L3 x *2 U  L y *2

Dans un tel cas d‟écoulement, l‟accroissement de la vitesse axiale est plus important suivant y que

suivant x.

u u
Par conséquent: ˂˂
x y


Avec : ˂˂ 1
L

Il en résulte que le premier terme visqueux de la première équation est négligeable devant le

deuxième:

  2u * L  2 u *
˂˂ VI.10
U  L x *2 U  2 y *2

2
L‟existence du terme ˂˂ 1 fait que tout son premier membre tend vers zéro.
L2

Les coefficients de l‟inégalité ci-dessous (VI.10) sont du même ordre de grandeur

 L2 L
 1
U L  2
U  2

 2
avec  est très faible devant l‟unité.
U L L2

On en déduit que les deux termes visqueux de l‟équation (VI.8) sont nuls.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 71

 2
Puisque  ce qui nous permet d‟écrire :
U L L2

L L
  VI.11
U Re L


Ainsi pour que l‟approximation ˂˂ 1 soit vraie, on doit avoir le nombre de Reynolds basé sur la
L

distance L suffisamment grand, ce qui explique pourquoi on classe les écoulements de couche limite

comme des écoulements à grands nombres de Reynolds.

Cependant, si l‟on se place à une abscisse x quelconque, on peut toujours considérer que la plaque

s‟arrête là, et que sa longueur est x. L‟épaisseur de couche limite δ(x) n‟en est évidemment pas

affectée. Aussi, avec un abus d‟écriture, on définit le « nombre de Reynolds local Rex » par :

 x  
L
VI.12
Re x

En régime laminaire le nombre de Reynolds est Re x  5  10


5

VI.2.2. Equations de la couche limite - équations de Prandtl (1904)

Les équations réduites de Navier-Stokes s‟écrivent en non dimensionnel comme suit :

 u * u * p * L  2u *
 u *  v *   
 x * y * x * U  2 y *2
 VI.13
 p *  0
 y *

Les équations de mouvement pour la couche limite laminaire ont alors pour expression:

Equations de quantité de mouvement

 u u 1 p   2 u 
u  v     
 x y  x  y 2 
 VI.14
 p  0
 y

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 72

On en déduit que dans une couche limite, la pression ne varie que suivant la direction de l‟écoulement

p  p x  .

A l‟extérieur de la couche limite le fluide se comporte comme étant un fluide parfait, en raison de

l‟absence de frottement visqueux nous pouvons appliquer la relation de Bernoulli pour écrire l‟égalité

suivante :

U 2
p  ρgh  ρ  cst VI.15
2

U 2
p  ρ  cst VI.16
2

dp dU 
 ρU  0 VI.17
dx dx

1 dp dU 
  U VI.18
 dx dx

En remplaçant dans l‟équation de quantité de mouvement, on aura les équations de la couche limites

suivantes :

L‟équation de continuité :

u v
 0 VI.19
x y

L‟équation de quantité de mouvement :

u u dU    2u 
u v  U    2  VI.20
x y dx  y 

Et les conditions aux limites :

y  0  uv0
 VI.21
y    u  U

On va s‟intéresser dans la suite de ce cours à l‟analyse des solutions de cette équation et à l‟obtention

de telles solutions, exactes (dans certains cas particuliers) ou approchées (dans le cas général).

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 73

VI.3. Résolution des équations de la couche limite

L‟obtention d‟une solution des équations de couche limite suppose la résolution analytique du système

d‟équations aux dérivées partielles non-linéaires muni des conditions aux limites appropriées. On ne

sait pas résoudre ce problème dans le cas général et on se contente donc de solutions approchées. En

outre, leur résolution numérique complète est extrêmement ardue. Cependant, sous certaines

conditions de nature expérimentale, des simplifications sont possibles qui permettent de faciliter

grandement la résolution du problème où il est possible de déterminer des solutions "quasi-exactes"

des équations de la couche limite.

VI.3.1. Solution de Blasius

L‟écoulement sur une plaque plane, c‟est-à-dire à vitesse extérieure constante, représente l‟exemple de

la couche limite le plus simple. C‟est le cas de Blasius.

- En dehors de la couche limite, l‟écoulement est uniforme, sans gradient de vitesse :

U∞= Cst

dp x 
- Ce qui implique à partir de l‟équation d‟Euler que 0
dx

Les équations de couche limite sont donc:

L‟équation de continuité :

u v
 0 VI.22
x y

L‟équation de quantité de mouvement :

u u   2u 
u v    2  VI.23
x y  y 

C‟est le problème dont la solution a été obtenue par Blasius. Comme l‟écoulement est incompressible,

 
le champ de vitesse est déterminé par une fonction de courant, ψ, telle que u  et v   .
y x

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 74

L‟équation de continuité est alors satisfaite. La similitude de la solution requiert que:

u
 g ( ) VI.24
U

η est un paramètre sans dimension construit à partir de x et de y :

y x
où  x  
x
 
 x 
VI.25a,b
Re x U

Pour exprimer les différents termes de l‟équation de quantité de mouvement, il nous faut d‟abord

exprimer la fonction de courant ψ :

  y
  u x   U    x   g   VI.26
 y 

La fonction de courant est nécessairement de la forme :

  U   x   G  VI.27

Avec :

g    G'   
u
VI.28
U

D‟où :

u  U G'   VI.29

Exprimons ensuite la composante v de la vitesse :

   x  G    x  G   
v  U   G    x   U   G    x 
  x 
VI.30
x  x x   x

Avec :

  ' x 
  2  ' x   
y
(Théorème de dérivation des fonctions composées) VI.31
x  x   x 

D‟où :


v  U   ' x G    ' x G   VI.32
x

v  U  ' x G   G VI.33

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 75

Ainsi nous avons :

u  U  G 

v  U   ' G   G 
VI.34

Nous devons encore calculer les autres termes :

 u u  '
   U  G 
 x  x 
 u u  U 
   G  VI.35
 y  y 
 2u U
 2  2 G 
 y 

On injecte ensuite ces expressions dans l‟équation de quantité de mouvement suivant x et on en

déduit :

U ' '
 G   U 2 G G   U 2 G   G G  VI.36
 2
 

U '
 G   U 2 G   G G   G G  VI.37
 2

U '
 G   U 2 GG   0 VI.38
 2

 '
G   U  GG   0 VI.39

 '
Ceci impose que le coefficient du second terme soit constant U     cst

Cependant, le choix de la constante est arbitraire, et ne modifie pas le résultat final. La convention

usuelle, due à Blasius, consiste à prendre :

1

2

L‟équation de Blasius s‟écrira donc en pratique :

2G   GG   0 VI.40

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 76

Nous avons maintenant en main une équation différentielle à une inconnue F() et une seule variable

, qui remplace le système de deux équations aux dérivées partielles à deux inconnues u(x,y) et v(x,y).

Les conditions aux limites associées à l‟équation de Blasius sont :

- A la paroi :

 u ( 0)
u  0  G 0  0
y  0   0  U VI.41
v  0  G (0)  0

- Loin de la paroi :

y     u U  G   1 VI.42

Au final l‟équation de Blasius s‟écrit :

2G   GG   0

G 0  0 , G 0  0 VI.43
G    1

Cette équation n‟a pas de solution analytique. Elle doit donc être résolue par intégration numérique.

Pour ce faire, on la réécrit sous forme d‟un système de trois équations différentielles ordinaires du

premier ordre:


G    F  

 F    H    G   VI.44

 H     G  H    G 
1
 2

Le système est donc sous la forme :

 
 
 G     F   
d  
 F      H    VI.45
d    
 H     1 G  H  
 
 2 

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 77

Le système d‟équation peut être intégrer numériquement par la „‟méthode de tir‟‟. On débute

l‟intégration en   0 avec G0  0 et F 0  G 0  0 . En suite prendre une valeur pour

H 0  G 0 , intégrer numériquement jusqu‟à   1 en suite itérer la valeur H 0 jusqu‟à obtenir

F    G    1 quand   1 .

  
 G  , de vitesse,  G   et de contrainte de
u
Les profils de fonction de courant,
U  x  U

 x 
cisaillement,  G   , ainsi obtenus sont présentés à la figure VI.2.
U 

Figure VI.2 Solution de Blasius pour la couche limite avec U constant.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 78

η G G' G''
0 0 0 0,3
0,2 0 0,06 0,3
0,4 0,012 0,12 0,3
0,6 0,036 0,18 0,29964
0,8 0,072 0,239928 0,298561
1 0,119986 0,29964 0,296412
1,2 0,179914 0,358923 0,292855
1,4 0,251698 0,417494 0,287586
1,6 0,335197 0,475011 0,280348
1,8 0,430199 0,53108 0,270951
2 0,536415 0,585271 0,259294
2,2 0,653469 0,637129 0,245385
2,4 0,780895 0,686207 0,22935
2,6 0,918136 0,732077 0,21144
2,8 1,064552 0,774365 0,192027
3 1,219425 0,81277 0,171585
3,2 1,381979 0,847087 0,150661
3,4 1,551396 0,877219 0,12984
3,6 1,72684 0,903187 0,109697
3,8 1,907477 0,925127 0,090754
4 2,092503 0,943278 0,073443
4,2 2,281158 0,957966 0,058075
4,4 2,472752 0,969581 0,044827
4,6 2,666668 0,978547 0,033743
4,8 2,862377 0,985295 0,024745
5 3,059436 0,990244 0,017662
5,2 3,257485 0,993776 0,012258
5,4 3,45624 0,996228 0,008265
5,6 3,655486 0,997881 0,005408
5,8 3,855062 0,998963 0,003431
6 4,054855 0,999649 0,002109
6,2 4,254785 1,000071 0,001254
6,4 4,454799 1,000321 0,00072
6,6 4,654863 1,000466 0,000399
6,8 4,854956 1,000545 0,000213
7 5,055065 1,000588 0,00011

Tableau VI.1 Solutions numériques de l‟équation de Blasius

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 79

VI.3.2. Contrainte de cisaillement et coefficient de frottement à la paroi

La contrainte visqueuse locale à la paroi s‟écrit :

 u   u     U  U 
 p            G 0  G 0 VI.46
 y  y 0    0  y  y 0  x  x
U

1
   2 2
3
p   U  G 0 VI.47
 x 

La valeur numérique de G‟‟(0)=0.3 d‟où :

1
   2 2
3
 p  0.3    U VI.48
 x 

Le coefficient adimensionnel de frottement local s‟écrit alors :

1
p   2 
1
Cf   0.6     0.6 Re 2 VI.49
1
U 2  U  x 
2

VI.3.3. Epaisseur de la couche limite

Par convention l‟épaisseur de la couche limite  x  est définie par :

u  
 0.99 VI.50
U

u  
D‟autre part :  G  
U

 0.99
Posons :   
 x 

u  
 G     0.99 VI.51
U

La résolution numérique de l‟équation de Blasius donne :   5 , d‟où :

 0.99 1
 0.99  5 x 

 5  Re x 2 VI.52a,b
x

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VI : COUCHE LIMITE 80

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 81

VII. ANALYSE DIMENSIONNELLE ET

SIMILITUDE

VII.1. Introduction

Dans cette partie on regardera deux sujets complémentaires: l‟analyse dimensionnelle et la similitude.

L‟analyse dimensionnelle est une méthode qui s‟intéresse à l'obtention d‟équations décrivant des

phénomènes à partir de la simple supposition que ceux-ci peuvent être caractérisés par certaines

variables. Le résultat d'une analyse dimensionnelle est l‟obtention de groupes sans dimensions

équivalents au système dimensionnel. Cette opération réduit le nombre de variables décrivant le

problème.

La similitude, est une méthode de travail très utilisée en ingénierie permettant l‟étude de problèmes sur

la base d‟expérimentation et de comparaison. Elle est souvent appliquée lorsqu‟on ne peut pas

effectuer un essai en grandeur nature. Dans ce cas, on fait une étude sur maquette et la similitude

permet d‟extrapoler les résultats vers le cas où l‟objet réel, appelé le prototype.

VII.2. Analyse dimensionnelle

Les objectifs spécifiques de l‟analyse dimensionnelle sont:

1. Réduire des variables dimensionnelles sous la forme de groupes sans dimension pourvus d‟un

sens physique

2. Faciliter l‟échange d‟information entre des modèles à échelles différentes

3. Diminuer le nombre d‟expériences à réaliser

Plusieurs méthodes ont été développées pour l‟obtention de paramètres sans dimension.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 82

Ces méthodes reposent sur un système cohérent d‟unités. En mécanique des fluides, on utilise le

système Masse, Longueur et Temps (MLT) ou le système Force, Longueur et Temps (FLT).

On dit qu'un système d'unités est cohérent avec un groupe d‟équations quand celles-ci sont satisfaites

dans chaque cas particulier.

VII.2.1. Dimension fondamentales

Le système international d'unités définit sept unités fondamentales, associées justement aux sept

grandeurs de base que nous avons utilisées pour écrire l'équation aux dimensions d'une grandeur. Ces

unités sont définies par un phénomène particulier. Elles sont rappelées dans le tableau 1.

Grandeur Dimension SI Symbole


Longueur Mètre m L

Masse Kilogramme kg M
Temps Seconde s T

Température Kelvin K θ
Intensité électrique ampère A I
Intensité lumineuse candela cd J

Quantité de matière mole mol N

Tableau VII.1 Les dimensions fondamentales

VII.2.2. Homogéniété dimensionnelle

Tester l'homogénéité d'une expression est un critère permettant d'éliminer des résultats dont on sait

qu'ils sont nécessairement faux.

Une équation est homogène lorsque ses deux membres ont la même dimension.

Le critère de pertinence s'énonce ainsi :

Une expression non homogène est nécessairement fausse.

On peut énoncer les conséquences suivantes :

- On ne peut additionner que les termes ayant la même dimension.


A. BENSLIMANE 2016/2017 DGM/FT/UAMB
CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 83

- L'argument d'une fonction transcendante (sin, cos, tan, exp, ln, ch, sh, th) doit être

sans dimension.

On manipule les dimensions à l'aide des règles suivantes :

- La dimension du produit de deux grandeurs est le produit des dimensions de chacune

des grandeurs.

- La dimension de Ar est égale à [A]r, ou r est un nombre sans dimension.

Exemple :

L‟équation de Bernoulli et dimensionnellement homogène, en effet :

p v2
H  z  VII.1
g 2 g

H  L
ML1T 2

LT 1 
 LLL

2

VII.2
ML2T 2 LT 2

VII.2.3. Méthode de Rayleigh

On peut exprimer une variable (n) en fonction de n paramètres (yi) en utilisant l‟analyse

adimensionnelle :

n  y1 a y2 b ......yn s VII.3

Avec : a, b, s sont des coefficients à déterminer.

Exemple d‟application :

Variables Symboles Dimensions


Chute de pression par unité de longueur Δp/L M.L-2.T-2
Diamètre de la conduite D L
Vitesse moyenne du fluide Umoy L.T-1
Viscosité dynamique m M.L-1.T-1
Masse volumique ρ M.L-3
Rugosité ε L
Tableau VII.2 Variables du problème

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 84

Dans cet exemple nous allons nous intéresser à l‟étude de l‟écoulement d‟un fluide incompressible

dans une conduite rectiligne.

Les variables retenues et leurs dimensions sont listées dans le tableau ci-dessus.

En observant que le diamètre D, la vitesse de débit Umoy, la viscosité µ, la masse volumique ρ du fluide

et la rugosité ne peuvent dépendre de la longueur L de la conduite, on peut écrire a priori que le

gradient linéique de pression doit s‟exprimer en fonction des autres variables:

p
 f  ,  ,  , D,U moy  VII.4
L

Développons la fonction f en série de puissance du type:

p
  a  b cU moy
d
De VII.5
L

Cette équation doit être dimensionnellement correcte, c‟est-à-dire que dans l‟équation aux dimensions:

ML 2
 
T 2  k ML3  ML T  L LT  L
a 1 1 b c 1 d e
VII.6

Les exposants de chaque dimension fondamentale doivent être égaux. Par conséquent :

M   1  a  b

L   2  3a  b  c  d  e VII.7
T    2  b  d

On a trois équations à cinq inconnues !!!

Grandeurs privilégiées : le diamètre, la masse volumique et la vitesse.

Les grandeurs secondaires : la viscosité et la rugosité

De ce fait, on exprime : a, d, e en fonction de b, c

a  1  b

e  1  c  b VII.8
d  2  b

p
 k 1c  b cU moy
2 b
D 1cb VII.9
L

k : est un paramètre sans dimension

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 85

b
p U moy
2
     c
k     VII.10
L D  U moy D   D 

Au final, on voit apparaître le nombre de Reynolds et la rugosité relative. On pourra donc écrire:

p U moy
2

k f Re,  * VII.11
L D


Avec :  *  la rugosité relative (adimensionnelle)
D

VII.2.4. Théorème π de Vaschy-Buckingham

Le travail de Aimé Vaschy sur les lois de similitude en physique a été présenté en 1892 dans Annales

Télégraphiques 19 (1892), 25-28 Deux années auparavant, dans les préliminaires de son ouvrage:

Traité d‟électricité et de magnétisme, il explique la réduction du nombre de paramètres lorsqu'on

regroupe des quantités physiques pour générer des nombres adimensionnels

Selon le théorème de Vaschy-Buckingham (ou théorème π), dans un problème comprenant n

grandeurs physiques où il y a m dimensions fondamentales, on peut réécrire ces grandeurs physiques

en (n-m) paramètres adimensionnels indépendants.

VII.2.4.1. Enoncé

Soient Α1 , Α2 , Α3 ,..., Αn les différentes grandeurs physiques: comme la vitesse, la pression, la

viscosité, etc. Entre toutes ces quantités, il y a une relation de la forme:

  A1  A2  ......  An   0 VII.12

Si π1, π2, π3, …. représentent les quantités adimensionnelles parmi les quantités physiques Α1 , Α2 , Α3

,..., Αn , on peut alors écrire une équation de la forme:

f  1   2  ......   nm   0 VII.13

Ou :

 1  f  2  ......   nm   0 VII.14

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 86

VII.2.4.2. Méthodologie

Pour effectuer une analyse dimensionnelle, on doit considérer les étapes suivantes:

1. Dresser la liste de toutes les quantités physiques Ai et leur dimension correspondante.

2. Écrire la fonction :   A1  A2  ......  An   0

3. Choisir les variables répétitives. Ces variables doivent contenir toutes les m dimensions du

problème. Souvent, on retient une variable parce qu‟elle détermine l‟échelle, une autre, parce

qu‟elle détermine les conditions cinématiques; il faut une variable liée avec la masse ou les

forces du système. Exemple: on peut retenir D, v et ρ comme variables fondamentales.

4. Écrire les paramètres π en fonction des exposants inconnus:

 1  v x D y  z   LT 1  L y ML3  ML1T 1   M 0 L0T 0


x1 1
z1
1 1 1

S‟assurer que toutes les quantités Ai sont incluses dans les groupes πi.

5. Écrire les équations des paramètres π pour les exposants; on doit obtenir une somme

algébrique nulle pour chaque dimension.

6. Résoudre les équations simultanément.

7. Remplacer les exposants trouvés (x1 , y1 , z1 ,...) dans les expressions de π (étape 4) pour

obtenir les paramètres π sans dimension.

8. Déterminer la fonction: f  1   2  ......   nm   0

S‟assurer que tous les paramètres πi sont indépendants les uns des autres.

9. Mettre les résultats sous la forme de nombres sans dimension connus (Re, Fr, Eu,...).

VII.2.4.3. Application

Dans un écoulement externe sur une plaque plane, l‟épaisseur de la couche limite δ dépend de la

distance x, la vitesse U, la viscosité µ et la masse volumique ρ du fluide.

  , x,U  ,  ,    0 VII.15

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 87

Le nombre de variables est 5. Le nombre de produits Pi doit être 2 et ils ont la forme suivante :

 1   x    M L T
 a b c 0 0 0

 VII.16
 2   x  U   M L T
 a b c 0 0 0

Nous obtenons :

 
 1  x
 VII.17
 2  U  x
 

Le théorème de Vaschy-Buckingham implique :

 
f  1 ,  2   f  , Re x   0 VII.18
x 

D‟où :


 f Re x  VII.19
x

Ce qui correspond au résultat de la théorie présentée dans le chapitre précédent.

VII.3. Similitude

Alors que l‟étude expérimentale est le plus souvent nécessaire, les essais en vraie grandeur (par

exemple d‟un barrage, d‟une construction portuaire, de navire en bassin de carène, d‟une turbine ou

d‟un avion, etc ) sont rares et très coûteux, et ne sont possible que dans un nombre limité de cas. C‟est

pourquoi on fait souvent recours aux modèles aux échelles réduites (appelés maquettes) du système en

vraie grandeur (appelé prototype) à étudier.

Pour que les résultats obtenus expérimentalement sur maquette puissent être extrapolés au prototype,

un nombre des principes de similitude doit être respectés. Pour mettre en évidence ces principes on

écrit d‟abord les équations de Navier–Stokes, ainsi que les conditions aux limites et initiales, sous une

forme adimensionnelle. Ainsi, nous employons les symboles avec astérisque pour désigner les

grandeurs sans dimensions.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 88

On pose :

u  Uu * x  Dx* p  U 2 p *

v  Uv y  Dy *
*
U * VII.20
w  Uw* t t
 z  Dz * D

Les équations de Navier-Stockes s‟écriront :

  u * * u
*
u * * u 
*
p *   2u *  2u *  2u * 
 
  *  u  v  w * 
      *2  *2  *2 
  t x *
 y *
 z  x *
 x y z 
  v *
 * v * v * v  p *   2v*  2v*  2v* 
* * *
 
  *  u  v  w       *2  *2  *2  VII.21
  t x * y * z *  y *  x y z 
  w*
* w * w * w  p 1   2 w*  2 w*  2 w*  1
* * *
  u v w     *2  *2  
  t * x * y * z *  z Re  x *2 y z  Fr

Les équations résultantes sont sans dimensions et ne font intervenir que deux nombres sans

dimensions qui dépendent de l‟´ecoulement. Cela traduit alors le fait que tout écoulement n‟est défini

que par les valeurs de nombres sans dimensions qui lui caractérisent.

 UD
 Re 
 
 2
VII.22
 Fr  U
 gD

Fr est appelé nombre de Froude, et l‟on reconnaît le nombre de Reynolds Re.

Notons que les conditions aux limites s‟appliquent aux frontières du système, donc liées à la géométrie

du système, et les conditions initiales expriment l‟état du système à l‟instant initial.

Pour que deux écoulements relatifs respectivement à la maquette et au prototype soient semblables, il

faut qu‟ils soient décrits par la même équation dynamique ci-dessus. Les variables réduites étant sans

dimension, il suffit d‟avoir l‟égalité des nombres de Reynolds et l‟égalité des nombres de Froude entre

la maquette et le prototype

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 89

VII.3.1. Nombre adimensionnel classique en mécanique des fluides

ρV L ρV 2
Re   Nombre de Reynolds
μ μV/ L

F L2
Cx  Coefficien t de traînée
ρV 2 2

V2  V 
Fr   ou Fr  Nombre de Froude
gL  g L 

ΔP Nombre d' Euler, ou nombre de pression


Eu 
ρV 2 2

ρV 2 2  ρL
We   ou We  V  Nombre de Weber
σ L  σ 

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


CHAPITRE VII : ANALYSE DIMENSSIONNELLE ET SIMILITUDE 90

A. BENSLIMANE 2016/2017 DGM/FT/UAMB


REFERENCES 91

REFERENCES

ANCEY C., MECANIQUE DES FLUIDES UNE INTRODUCTION A L‟HYDRAULIQUE POUR LES

INGENIEURS CIVILS , 2016, ÉCOLE POLYTECHNIQUE FEDERALE DE L AUSANNE.

BENHAMOUDA R., NOTIONS DE MECANIQUE DES FLUIDES , COURS ET EXERCICES CORRIGES ,

2008, CENTRE DE PUBLICATION UNIVERSITAIRE , TUNIS.

CHASSAING P., MECANIQUE DES FLUIDES - É LEMENTS D 'UN PREMIER PARCOURS , 2000,

CEPADUES E DITIONS .

COMOLET R., MECANIQUE EXPERIMENTALE DES FLUIDES , 1982, MASSON.

GHAZI B., MECANIQUE DES FLUIDES UNE INTRODUCTION A LA DYNAMIQUE DES FLUIDES REELS

INCOMPRESSIBLES , 2015, ECOLE NATIONALE D‟INGENIEURS DE T UNIS.

RIDHA A., DYNAMIQUE DES FLUIDES REELS, 2008, UNIVERSITE DE CAEN

ROBINET J.C., AERODYNAMIQUE FONDAMENTALE, 2013, ARTS & METIERS P ARIS T ECH.

VAZQUEZ J., HYDRAULIQUE GENERAL , ECOLE NATIONALE DU GENIE DE L‟EAU ET DE

L‟ENVIRONNEMENT DE S TRASBOURG

VAZQUEZ J., DUFRESNE M., HYDRAULIQUE A SURFACE LIBRE, 2016, ECOLE NATIONALE DU

GENIE DE L‟E AU ET DE L‟ENVIRONNEMENT DE S TRASBOURG

WHITE F.M., FLUID MECHANICS , UNIVERSITY OF RHODE ISLAND, WCB MCGROW-HILL

WINCKELMANS G., ET LEGAT V., MECANIQUE DES FLUIDES ET T RANSFERTS I, 2009, ECOLE

POLYTECHNIQUE DE LOUVAIN.

A. BENSLIMANE 2016/2017 DGM/FT/UAMB

View publication stats

Vous aimerez peut-être aussi