MecaFlu EMSE
MecaFlu EMSE
MecaFlu EMSE
2
Introduction
Le présent document est destiné aux étudiants de l'Ecole Nationale Supérieure des Mines de Saint-Etienne.
Il présente les bases de la mécanique des uides.
Bonne lecture.
Olivier Bonnefoy
Nota Bene : ce document est en cours d'élaboration. Il peut évidemment comporter des inexactitudes ou des
erreurs. Merci de bien vouloir en avertir l'auteur (bonnefoy@emse.fr). Il vous en sera reconnaissant et intégrera vos
i
Introduction
ii
Table des matières
1 Analyse dimensionnelle 1
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Méthodes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2.1 Former des grandeurs adimensionnées . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2.2 Adimensionner une loi physique inconnue . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2.3 Adimensionner une loi physique connue . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.4 Transposer des résultats d'une échelle à l'autre . . . . . . . . . . . . . . . . . . . . . . 7
1.2.5 Identier l'importance d'un phénomène . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3 Application à la mécanique des uides et aux transferts . . . . . . . . . . . . . . . . . . . . . 10
1.3.1 Echelles caractéristiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.2 Phénomènes les plus fréquents et temps caractéristiques associés . . . . . . . . . . . . 10
1.3.3 Nombres adimensionnés les plus fréquents . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.4 Adimensionnement de l'équation de Navier-Stokes . . . . . . . . . . . . . . . . . . . . 11
1.4 Solutions des exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
iii
Introduction
4 Bilans locaux 51
4.1 Forme générale d'un bilan local . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.2 Bilan local de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.2.1 Pour chaque constituant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.2.2 Pour le uide dans son ensemble . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
4.2.3 Cas particulier d'un uide incompressible . . . . . . . . . . . . . . . . . . . . . . . . . 54
4.3 Bilan local de quantité de matière . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
4.3.1 Liens entre les bilans de masse et de quantité de matière . . . . . . . . . . . . . . . . . 55
4.4 Bilan local de quantité de mouvement (linéaire) . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.4.1 Cas général . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
4.4.2 Cas particulier des uides newtoniens (équation de Navier-Stokes) . . . . . . . . . . . 57
4.5 Bilan local d'énergie mécanique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.5.1 Cas général . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
4.5.2 Commentaire sur le sens physique du terme σ T : (∇~v ) . . . . . . . . . . . . . . . . . . 60
4.5.3 Cas particulier pour un écoulement incompressible et stationnaire . . . . . . . . . . . 61
4.6 Bilan local d'énergie interne et d'enthalpie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
4.6.1 Bilan local d'énergie interne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
4.6.2 Bilan local d'enthalpie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.7 Conditions aux limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.7.1 Mécaniques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
4.7.2 Thermiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.7.3 Chimiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.8 Synthèse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.9 Solutions des exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
5 Bilans macroscopiques 69
5.1 Les diérents volumes de contrôle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
5.2 Expression générale du bilan macroscopique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
5.2.1 Volume de contrôle quelconque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
5.2.2 Volume de contrôle matériel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.2.3 Volume de contrôle imaginaire xe . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.3 Bilan macroscopique de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.4 Bilan macroscopique de quantité de mouvement . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.5 Bilan macroscopique d'énergie (totale) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
iv
Introduction
6 Compléments 75
6.1 Corrélations empiriques pour écoulements typiques . . . . . . . . . . . . . . . . . . . . . . . . 75
6.1.1 Ecoulement en conduite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.1.2 Ecoulement le long d'un plan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
6.1.3 Ecoulement autour d'une sphère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
7 Annexes mathématiques 79
7.1 Symboles de Kronecker et Levi-Civita, convention d'Einstein . . . . . . . . . . . . . . . . . . 79
7.2 Tenseurs et produits tensoriels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
7.2.1 Tenseurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
7.2.2 Produits tensoriels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
7.3 Opérateurs diérentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
7.3.1 Dénitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
7.3.2 Formulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
7.4 Conventions sur la construction du tenseur des contraintes . . . . . . . . . . . . . . . . . . . . 83
7.4.1 Normale entrante/sortante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
7.4.2 Assemblage vertical/horizontal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
7.4.3 Synthèse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
7.5 Démonstrations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
7.5.1 Bilan des forces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
7.5.2 Bilan des moments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
7.5.3 Systèmes électro-magnétiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.5.4 Systèmes en rotation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.5.5 Théorème Π de Vaschy-Buckingham . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.5.5.1 Notations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.5.5.2 Enoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.5.5.3 Démonstration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
7.6 Théorèmes importants . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
7.6.1 Théorème de Green-Ostrogradsky . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
7.6.2 Théorème du gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.6.3 Théorème de transport de Reynolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
7.7 Bilans en coordonnées cartésiennes (~ex , ~ey , ~ez ) . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
7.7.1 Produits tensoriels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
7.7.2 Opérateurs diérentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
7.7.3 Bilan local de masse et de quantité de matière . . . . . . . . . . . . . . . . . . . . . . 94
7.7.3.1 Bilan local de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
7.7.3.2 Loi de Fick pour la diusion binaire . . . . . . . . . . . . . . . . . . . . . . . 94
7.7.3.3 Bilan local de quantité de matière . . . . . . . . . . . . . . . . . . . . . . . . 94
7.7.4 Bilan local de quantité de mouvement (cas général) . . . . . . . . . . . . . . . . . . . . 94
7.7.5 Tenseur des taux de cisaillement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
7.7.6 Tenseur des contraintes de cisaillement pour un uide newtonien . . . . . . . . . . . . 95
7.7.7 Bilan local de quantité de mouvement (équation de Navier-Stokes) . . . . . . . . . . . 95
7.7.8 Bilan local d'énergie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
7.7.8.1 Equation de la chaleur (cas général) . . . . . . . . . . . . . . . . . . . . . . . 96
7.7.8.2 Loi de Fourier pour la diusion de chaleur . . . . . . . . . . . . . . . . . . . 96
7.7.8.3 Equation de la chaleur (uide newtonien incompressible) . . . . . . . . . . . 96
7.8 Bilans en coordonnées cylindriques (~er , ~eθ , ~ez ) . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
7.8.1 Produits tensoriels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
7.8.2 Opérateurs diérentiels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
7.8.3 Bilan local de masse et de quantité de matière . . . . . . . . . . . . . . . . . . . . . . 101
7.8.3.1 Bilan local de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
7.8.3.2 Loi de Fick pour la diusion binaire . . . . . . . . . . . . . . . . . . . . . . . 102
v
Introduction
vi
Chapitre 1
Analyse dimensionnelle
I have often been impressed by the scanty 1 attention paid even by original workers in physics
to the great principle of similitude. It happens not infrequently that results in the form of "laws"
are put forward as novelties on the basis of elaborate experiments, which might have been predicted
a priori after a few minutes' consideration.
1.1 Introduction
L'analyse dimensionnelle est un ensembe de techniques permettant de :
diminuer le nombre de grandeurs nécessaires à la description d'un phénomène physique
minimiser le nombre de mesures nécessaires à la détermination d'une loi physique inconnue
étudier les solutions asymptotiques d'une loi physique en identiant les termes prépondérants
établir un lien entre des mesures obtenues sur le dispositif réel et sur un modèle réduit
Dénitions
Dans un problème de physique, plusieurs grandeurs physiques xi interviennent : diamètre, vitesse,
durée, masse, . . . Elles sont liées par une loi physique qui, bien souvent, permet de connaître la réponse
du système lorsqu'une sollicitation lui est appliquée :
f (x1 , . . . , xn ) = 0
Dans certains cas particuliers, la théorie permet d'établir une expression analytique de cette loi. Par
exemple,
√ la vitesse de vidange d'un réservoir sous l'eet de la gravité est donnée par la formule de Torricelli
v = 2gh et l'intensité du courant électrique I induit par la diérence de potentiel U est donnée par la loi
d'Ohm I = U/R. Dans la grande majorité des situations intéressantes pour un ingénieur, toutefois, on sait
que les grandeurs sont liées mais il n'est pas possible d'écrire explicitement une équation ou de la résoudre.
Par exemple, la puissance P requise pour qu'un avion vole à la vitesse ~v dans une zone étendue de turbulences
n'est pas calculable par une équation simple. On peut simplement écrire f (~v , P, . . . ) = 0.
Chaque grandeur peut s'exprimer avec un nombre inni d'unités. Par exemple, le diamètre d'une sphère
peut s'exprimer en µm, cm, pouces, miles, km, année-lumière, . . . Il est important d'être très vigilant sur les
unités car les conséquences d'une erreur peuvent être funestes. A titre d'illustration, la sonde spatiale "Mars
Climate Orbiter" (193 millions de dollars) n'a pas atteint sa cible car les unités étaient Foot-Pound-Second
chez le constructeur (Lockheed Martin) et Mètre-Kilogramme-Seconde chez le gestionnaire des communica-
tions entre le satellite et le centre de contrôle de Pasadena (Jet Propulsion Laboratory).
1
Chapitre 1. Analyse dimensionnelle
En revanche, il existe un nombre ni de dimensions. Pour toute la physique, on en recense seulement
sept. Le choix des sept dimensions de base comporte une part d'arbitraire. Celles retenues par le Système
International sont : masse (M ), longueur (L), temps (T ), température (θ), intensité électrique (I ), intensité
lumineuse (J ) et quantité de matière (N ).
L'équation aux dimensions d'une grandeur quelconque x relie sa dimension, notée [x] à celles des
dimensions de base. Dans le Système International, la forme générale est :
De manière triviale, une grandeur sans dimension est telle que tous les exposants ai sont nuls 2 .
Table 1.1 Equations aux dimensions des grandeurs usuelles en mécanique des uides.
Exercice 1. en utilisant le Système International, donner l'équation aux dimensions du produit ρ.V 2
et du produit ρ.ν 2 . Commenter.
2. On remarquera qu'il existe des grandeurs sans dimension mais avec une unité. Par exemple, l'angle est le rapport de deux
longueurs (l'arc divisé par le rayon) et s'exprime en radian.
2
1.2. Méthodes
1.2 Méthodes
1.2.1 Former des grandeurs adimensionnées
Méthode
On décrit ici la procédure pour construire un ensemble de p grandeurs adimensionnées à partir de n grandeurs
physiques. Cette procédure est inspirée du théorème Π de Vaschy-Buckingham (démonstration en page 89).
1. identier les n grandeurs physiques xi du problème. Il peut s'agir de paramètres constants, de fonctions
ou de variables. Il est important de toutes les recenser.
2. écrire l'équation aux dimensions de chaque grandeur xi :
M
Y
ami
[xi ] = Dm
m=1
aM 1 ... aM n
Dans un problème donné, il est possible d'avoir plusieurs échelles caractéristiques de même dimension
(deux temps caractéristiques, deux longueurs caractéristiques, . . .).
7. former n grandeurs sans dimension en prenant le ratio xi /Xi :
8. éliminer les r grandeurs adimensionnées sans intérêt : constantes, multiples ou puissances d'un autre
nombre adimensionné, . . .
3
Chapitre 1. Analyse dimensionnelle
Figure 1.1 Nécessité de choisir des grandeurs liées par une relation de cause à eet. Source.
Exemple
Considérons un objet sphérique chutant dans le vide sous l'eet de la gravité et essayons de prévoir la
durée de chute. Les grandeurs physiques a priori pertinentes sont la durée de la chute τ , la hauteur initiale
h0 , la vitesse initiale v0 , le diamètre d, la masse m et l'accélération de la pesanteur g . Le nombre n de
grandeurs est égal à 6. La loi physique s'énonce ainsi :
f (τ, h0 , d, m, g, v0 ) = 0
τ h0 d m g v0
M 0 0 0 1 0 0
L 0 1 1 0 1 1
T 1 0 0 0 -2 -1
θ 0 0 0 0 0 0
I 0 0 0 0 0 0
J 0 0 0 0 0 0
N 0 0 0 0 0 0
Le rang r de cette matrice est égal à 3. Nous choisissons autant de grandeurs de base dimensionnellement
indépendantes ; par exemple m, g et v0 . En utilisant ces grandeurs de base, formons maintenant autant de
nombres adimensionnés qu'il y a de variables physiques 3 :
N ≡ v0 .gτ −1
1
h0
N2 ≡ v02 .g−1
N3 ≡ v2 .gd −1
0
N4 ≡ m
m
N5 ≡ gg
N6 ≡ vv00
4
1.2. Méthodes
Les trois dernières grandeurs adimensionnées sont des constantes, égales à l'unité, et n'ont aucun intérêt.
Le nombre p de grandeurs adimensionnées pertinentes est donc égal à 3, ce qui était prévisible (p ≡ n − r).
Exercice 3. trouver la vitesse caractéristique d'un uide s'écoulant dans un milieu poreux sous l'eet
d'un gradient de pression. Application à de l'eau s'égouttant par gravité dans un pot de eur.
f (x1 , . . . , xn ) = 0 (1.1)
Si cette loi est inconnue, on cherchera à établir une corrélation empirique entre ces grandeurs. Pour minimiser
le nombre d'expériences ou de simulations numériques, on peut avantageusement former p grandeurs adi-
mensionnées selon la procédure décrite dans la section précédente. En eet, la loi physique sera équivalente 4
à l'équation suivante, dite "loi physique adimensionnée" :
dont le nombre de variables est strictement inférieur à n. Dans le cas particulier où il n'y a qu'un nombre
adimensionné (p = 1), alors cette équation s'écrit simplement N1 = constante. De manière générale, lorsque
les diérents phénomènes sont découplés, les corrélations ont souvent la forme d'un produit de fonctions,
elles-mêmes souvent lois de puissance :
Lorsque les phénomènes sont couplés, l'expression analytique des corrélations est plus dicile à établir.
Exemple
Dans le prolongement de l'exemple de la section précédente, on peut écrire les équivalences suivantes :
τ.g
N1 ≡ v0
f (τ, h0 , d, m, g, v0 ) = 0 ⇔ N1 = h (N2 , N3 ) avec N2 ≡ hv02.g
0
N ≡ d.g
3 v02
4. Pour établir cette équivalence, on peut essayer d'isoler la grandeur adimensionnée faisant apparaître la sollicitation du
système. Cette dernière est en principe unique. La réponse (physique) du système à cette sollicitation est souvent unique mais
peut présenter diérentes solutions si le système est dans un état métastable ou instable.
5
Chapitre 1. Analyse dimensionnelle
Dans le cas où la loi physique fait intervenir une dérivée ou une intégrale, on utilise les règles suivantes :
dx1 X1 dx∗1
= .
dx2 X2 dx∗2
Z Z
x1 .dx2 = X1 .X2 . x∗1 .dx∗2
L'adimensionnement d'une loi physique connue permet réduire le nombre de grandeurs et donc de faciliter son
calcul. Quand la loi se présente sous la forme d'une somme de plusieurs termes, l'adimensionnement permet
également d'évaluer l'ordre de grandeur des diérents termes et ainsi, plus facilement, de considérer les cas
asymptotiques où l'un ou l'autre des termes, et donc des phénomènes physiques associés, est négligeable ou
prépondérant devant les autres.
Exemple
Toujours dans le prolongement de l'exemple de la section précédente, la théorie permet de relier le dépla-
cement de l'objet à la durée de chute. Il sut d'appliquer le principe fondamental de la dynamique et de
réaliser une intégration :
dx
x(t) tel que = g.t + v0
dt
Parmi les 4 grandeurs physiques (x1 , x2 , x3 , x4 ) = (x, t, g, v0 ), les deux premières sont des fonctions ou des
variables tandis que les deux dernières sont des paramètres constants. Ces quatre grandeurs font intervenir
deux dimensions (L et T ) et le rang de la matrice dimensionnelle est r = 2. Pour les deux grandeurs de base,
il est naturel de choisir les paramètres constants : (b1 , b2 ) = (g, v0 ). Cela permet de construire les échelles
caractéristiques Xi et les grandeurs sans dimension suivantes :
v2
∗
X1 ≡ g0 x ≡ x
t∗ ≡ Xt 1
v0
X2 ≡ g
et X2
X 3 ≡ g
N3 ≡ Xg3
v0
N4 ≡ X
X4 ≡ v0
4
En remplaçant les grandeurs physiques xi par les produits x∗i .Xi ou Ni .Xi , on peut réécrire la loi physique :
dx d (x∗ .X1 )
= g.t + v0 ⇔ = (N3 .X3 ) . (t∗ .X2 ) + (N4 .X4 )
dt d (t∗ .X2 )
Ce qui conduit, en remplaçant les échelles caractéristiques Xi par leur valeur puis en divisant membre à
membre par v0 à la loi physique adimensionnée :
dx∗
x∗ (t∗ ) tel que = t∗ + 1
dt∗
Cette équation permet de voir que pour les temps courts (t∗ 1), la vitesse reste quasiment identique à sa
valeur initiale et, que pour des temps longs (t∗ 1) la vitesse croît proportionnellement au temps. De plus,
cette démarche permet de dénir le temps caractéristique vg0 .
Exercice 4. évaluer le temps caractéristique de la décharge d'un condensateur dans un circuit RC.
6
1.2. Méthodes
Exercice 5. évaluer le temps caractéristique des oscillations d'un pendule oscillant obéissant à la loi
θ̈ + g
L . sin θ = 0.
Considérons un problème de physique dont la loi physique adimensionnée, connue ou inconnue, s'énonce
de la manière suivante :
N1 = h (N2 , . . . , Np ) (1.3)
Puisque la loi physique est la même quelle que soit la taille du système considéré, l'égalité des nombres
adimensionnés N2 , . . . , Np aux deux échelles entraîne l'égalité de N1 aux deux échelles. Dans la littérature, le
système à pleine échelle est souvent appelé prototype et les grandeurs associées sont désignées par l'exposant
P ou bien 1 (écoulement principal, premier écoulement) tandis que le système modèle est appelé maquette
(exposant M ou 2 pour "second écoulement").
P
N2 = N2M
.. ..
. . ⇒ N1P = N1M
NpP = NpM
Similitude complète : on dit que deux systèmes sont similaires, ou en complète similitude, lorsque les
conditions suivantes sont remplies :
tous les nombres adimensionnés sont identiques
les conditions aux limites adimensionnées sont identiques
les conditions initiales adimensionnées sont identiques (pour les états instationnaires)
De la première condition, on peut déduire que la forme géométrique est invariante : dans une similitude
complète, les deux systèmes sont homothétiques l'un de l'autre. Parmi toutes les invariances, certains au-
teurs distinguent la similitude géométrique (invariance des rapports de longueur), la similitude cinématique
(invariance des rapports de vitesse) et la similitude dynamique (invariance des rapports de force).
7
Chapitre 1. Analyse dimensionnelle
Similitude incomplète : on parle de similitude incomplète lorsqu'au moins une condition de similitude
n'est pas remplie. Cette situation se rencontre notamment lorsque :
les matériaux imposés par la similitude sont trop malcommodes d'emploi ou bien tout simplement
n'existent pas
les sollicitations ne peuvent pas être contrôlées. Par exemple, l'intensité de la force de gravitation
est la même quelle que soit la taille du système (sauf chute libre, centrifugation et changement de
planète).
un phénomène est négligeable à une échelle et prépondérant, ou au moins inuant, à une autre échelle.
Par exemple, l'eet de la tension supercielle jouera probablement un rôle important pour un système
millimétrique mais sera négligeable pour un système de la taille de l'océan.
L'extrapolation des résultats n'est alors pas justiable rigoureusement d'un point de vue mathématique. En
pratique, il arrive toutefois que cette diculté puisse être négligée ou contournée :
si la fonction N1 = h (N2 , . . . , Np ) reste quasi-constante quelle que soit la valeur de Np (sur R ou un
intervalle de R), l'extrapolation sera correcte même si Np n'est pas identiques pour les deux systèmes,
en renonçant à une pure homothétie et en introduisant une distorsion géométrique (par exemple une
échelle de largeur et une de longueur), on peut tout de même obtenir des résultats.
Exemple
Cherchons à évaluer la puissance consommée par un sous-marin opérant à VP =2,6 m.s−1 dans une eau
de mer de masse volumique ρP =1010 kg.m−3 et de viscosité cinématique νP =1,3.10−6 m2 .s−1 . Pour cela,
on construit un modèle réduit dans un rapport LM /LP = 1/20 que l'on place dans un canal d'eau douce
(ρM =988 kg.m−3 et νM =0,65.10−6 m2 .s−1 ). En se plaçant dans des conditions de similitude, on mesure la
force de poussée FM qui s'exerce sur la maquette ; elle est de 500 kN.
Le problème est décrit par n = 5 grandeurs physiques : la taille L et la vitesse V du sous-marin, la masse
volumique ρ et la viscosité cinématique ν du uide, la force de traînée F . Ces grandeurs font intervenir
r = 3 dimensions (L,M,T). La loi physique du problème relie donc p = n − r = 2 nombres adimensionnés.
Choisissons r = 3 grandeurs de base dimensionnellement indépendantes ; par exemple la taille L [L], la vitesse
V [L.T−1 ] et la force F [M.L.T−2 ]. Nous pouvons construire n = 5 nombres sans dimension :
N1 ≡ L
L
N2 ≡ VV
N3 ≡ F.L−2ρ.V −2
ν
N ≡ L.V
4
F
N5 ≡ F
On voit que seuls deux nombres adimensionnés sont pertinents : N3 et N4 . La similitude entre le sous-marin
à pleine échelle et la maquette impose l'égalité des nombres adimensionnés pertinents :
ρP
= F .Lρ−2
(
P M
N3 = N3M FP .L−2 −2 −2
P M soit ν
P .VP M
ν
M .VM
N4 = N4 P
L .V = L .V
M
P P M M
De la deuxième équation, on tire la vitesse (très rapide) que doit avoir la maquette :
νM LP
VM = VP . . d'où VM = 26 m.s−1
νP LM
De la première équation, on tire la force de traînée du sous-marin à pleine échelle :
2
ρP νP
FP = FM . . d'où FP = 2045 kN
ρM νM
En utilisant la formule P = F.V , on trouve que la puissance consommée par le sous-marin réel est de 5,32
MW. On verra dans la section 1.3 que, à un coecient près, les nombres adimensionnés N3 et N4 sont les
inverses du coecient de traînée Cx et du nombre de Reynolds Re respectivement.
8
1.2. Méthodes
Exercice 7. on considère un modèle réduit au tiers de la pompe étudiée dans un exercice de la section
1.2.1. Calculer la puissance PP nécessaire pour faire tourner le prototype à 300 rpm ainsi que son débit
sachant que la puissance PM consommée par la maquette est de 2 hp lorsque cette dernière fonctionne
edans des conditions de similitude caractérisées par : ΩM = 900 rpm, DM = 5 in, ∆HM = 10 ft,
QM = 3 ft3 .s−1 .
Dans cette situation de découplage, le comportement du système est piloté par un seul des deux phéno-
mènes, c'est-à-dire qu'une variation d'un paramètre intervenant dans le temps caractéristique de ce phéno-
mène a une inuence directe (le plus souvent proportionnelle ou inversement porportionnelle) sur le compor-
tement de ce système. Deux cas principaux sont à distinguer :
Deux phénomènes en parallèle (1/τ = 1/τrapide + 1/τlent ≈ 1/τrapide ) : le système est piloté par le
phénomène le plus rapide. On peut négliger la contribution du phénomène lent.
Deux phénomènes en série (τ = τrapide + τlent ≈ τlent ) : le système est piloté par le phénomène le plus
lent. On peut négliger la contribution du phénomène rapide.
Ces réexions sont à la base des stratégies d'intensication des procédés : adapter la taille des systèmes, les
matériaux et les conditions opératoires de manière à maximiser l'ecacité 5 .
Exemple
Considérons le transfert de chaleur d'un solide vers un uide. On distingue trois phénomènes :
1. la conduction dans le solide avec un temps caractéristique : τ1 ≡ Ls
αs
Lf
2. la conduction dans le uide avec un temps caractéristique : τ2 ≡ αf
Lf .ρf .cpf
3. la convection dans le uide avec un temps caractéristique : τ3 ≡ h avec λ = α.ρ.cp
L'importance relative de ces phénomènes peut être évaluée par les trois nombres adimensionnés suivants :
h.L
N u ≡ λf f = ττ23 : nombre de Nusselt, utile pour comparer la conduction dans le uide à la convection
dans le uide
h.L ρ .c
Bi ≡ λsf = ττ13 . ρfs .cps
pf
: nombre de Biot, utile pour comparer la conduction dans le solide à la
convection dans le uide
α
N ≡ αfs = ττ12 : nombre sans nom particulier, utile pour comparer la conduction dans le uide à la
conduction dans le solide
5. Lire par exemple page 84 de "Green Process Engineering : From Concepts to Industrial Applications" By Martine Poux,
Patrick Cognet, Christophe Gourdon 2015, ISBN : 1482208172.
9
Chapitre 1. Analyse dimensionnelle
(N1 .N2 , N2 /N1 ). Pour certains nombres adimensionnés N , le tableau 1.4 indique la terminologie spécique
utilisée pour désigner l'état du système lorsque les valeurs de N sont très petites ou très grandes devant
l'unité (situation de découplage).
Les grandeurs utilisées dans ces tableaux sont les suivantes : vitesse V , distance L, diamètre d, coecient
de diusion DAB , viscosité cinématique ν , masse volumique du uide ρ, pression p, gravité g , vitesse du son
10
1.3. Application à la mécanique des uides et aux transferts
Peclet Pe V.L
DAB τdi mat /τconv
Schmidt Sc ν
DAB τdi mat /τvisc
Euler Eu 1
∆p
2 n/a
2 ρV
Froude Fr √V
g.L
τgrav /τconv
Mach Ma = M V
vs τson /τconv
∆ρ.g.L2 2
Bond Bo σ (τsurf /τgrav )
2
Capillaire Ca µ.V
σ (τsurf ) / (τconv .τvisc )
Ohnesorge Oh √ µ
ρ.σ.L
τsurf /τvisc
Fluide+Solide Stokes St
ρS .d2S .V
µ.L τfrein /τconv
F/A
Coe de traînée Cx = CD 1 2 n/a
2 ρ.V
Fluide+Thermique Prandtl Pr ν
α τdi therm /τvisc
Lewis Le α
DAB τdi mat /τdi therm
h.Lf
Nusselt Nu λf τdi therm /τconv therm
11
Chapitre 1. Analyse dimensionnelle
Stokes St Particules pilotées Deux phases décou- Couplage entre les deux phases d'une
par uide plées suspension
Froude Fr Subcritique (uvial) Supercritique (tor- Régime d'écoulement à surface libre,
rentiel) hydraulique uviale et maritime
Mach Ma Subsonique Supersonique Eet de la compressibilité du uide, aé-
rodynamisme
Peclet Peθ Particules soumises à Particules entraî- Rhéologie
agitation brownienne nées par uide
Les grandeurs physiques sont au nombre de n = 7 : ~v , p, x, t pour les fonctions ou variables et ρ, µ, ~g pour
les paramètres constants. Comme le rang r de la matrice dimensionnelle est égal à 3, il faut trouver trois
grandeurs de base bi . Comme on va le constater, il s'avère pratique de choisir L, V et ρV 2 . Les échelles
caractéristiques Xi et les grandeurs adimensionnées sont alors :
v ∗ ≡ Xv1 avec X1 ≡ V
p∗ ≡ Xp avec X2 ≡ ρV 2
2
x∗ ≡ Xx3 avec X3 ≡ L
t∗ ≡ Xt4 avec X4 ≡ L
V
2
N5 ≡ Xρ5 avec X5 ≡ ρV
V2
2
N6 ≡ Xµ6 avec X6 ≡ ρVV .L
V2
N ≡ g avec X7 ≡
7 X7 L
On s'aperçoit que les nombres adimensionnés N6 et N7 sont liés à des nombres classiques :
1 ρ.V.L
N6 = avec Re ≡
Re µ
1 V
N7 = avec Fr ≡ √
F r2 gL
En remplaçant les grandeurs physiques xi par les produits x∗i .Xi ou Ni .Xi , l'équation de Navier-Stokes
devient :
∂~v ∗ 1 1 ~g
+ (~v ∗ · ∇∗ ) ~v ∗ = −∇∗ p∗ + .∆∗~v ∗ + .
∂t∗ Re F r2 k~g k
avec les opérateurs ∇∗ ≡ L.∇ et ∆∗ ≡ L2 .∆. L'étude des solutions asymptotiques est facilitée. Par exemple,
en régime stationnaire et en négligeant la gravité, l'équation de Navier-Stokes devient :
1
∇∗ p∗ = − (~v ∗ · ∇∗ ) ~v ∗ + .∆∗~v ∗
Re
Pour un nombre de Reynolds très grand devant l'unité, l'équation de Navier-Stokes devient l'équation d'Euler,
valable pour un uide parfait (µ = 0) mais aussi pour un uide réel dans un écoulement turbulent (Re → ∞) :
∇∗ p∗ = −~v ∗ ∇∗~v ∗
12
1.3. Application à la mécanique des uides et aux transferts
Pour un nombre de Reynolds très petit devant l'unité, l'écoulement est dominé par les forces visqueuses et
l'on parle d'écoulement rampant (creeping ow) ou écoulement de Stokes. Dans ce cas, l'équation se réduit
à:
1 1 ~g
∇∗ p∗ = .∆∗~v ∗ + .
Re F r2 k~g k
13
1.4 Solutions des exercices
Exercice 1 l'équation aux dimensions du produit ρ.V 2 est [ρ.V 2 ] = M.L−1 .T −2 . On constate que ce dernier
produit est homogène à une pression. Celle du produit ρ.ν 2 est [ρ.ν 2 ] = M.L.T −2 . On constate que ce dernier
produit est homogène à une force.
P ρ µ Ω D Q ∆H
M 1 1 1 0 0 0 0
L 2 -3 -1 0 1 3 1
T -3 0 -1 -1 0 -1 0
Le rang r de cette matrice est égal à 3. Nous choisissons autant de grandeurs de base dimensionnellement
indépendantes ; par exemple ρ, Ω et D. En utilisant ces grandeurs de base, formons maintenant autant de
nombres adimensionnés qu'il y a de variables physiques :
N1 ≡ ρ.ΩP3 .D5
N2 ≡ ρρ
N ≡ µ
3 ρ.Ω.D 2
Ω
N4 ≡ Ω
N5 ≡ D
D
Q
N6 ≡ Ω.D 3
N ≡ ∆H
7 D
Les grandeurs adimensionnées N2 , N4 et N5 sont des constantes, égales à l'unité, et n'ont aucun intérêt. Le
nombre p de grandeurs adimensionnées pertinentes est donc égal à 4, ce qui était prévisible (p ≡ n − r). Pour
caractériser le fonctionnement de ce type de pompe, on cherchera la corrélation entre les quatre grandeurs
suivantes :
P Q ∆H µ
= f , ,
ρ.Ω3 .D5 Ω.D3 D ρ.Ω.D2
Exercice 3 les grandeurs physiques sont :
matériau liquide : viscosité dynamique µ
matériau solide : taille des pores d
sollicitation : gradient de pression ∆p/L ou bien gravité g
réponse : vitesse v
∆p/L
L'analyse montre qu'il y a un nombre adimensionné : µ.v/d 2 . Il est donc constant. Si l'on considère qu'il
∆p d2
est proche de l'unité, alors une première estimation de la vitesse est donnée par v ≈ L . µ . Application
2
2 1000.9,81.(10−4 )
numérique : v ≈ ρgd
µ ≈ 10 −3 ≈ 10 cm.s−1 .
Exercice 5
q
τ ≡ Lg .
1.4. Solutions des exercices
Exercice 6 Réponse courte : les transferts de matière par diusion sont souvent très lents par rapport aux
∗
transferts convectifs. Il sut de résoudre ∂c ∗ ∗ ∗
∂t∗ + V .∇ c = −N6 .c
∗
avec N6 = V0 /L k
0
.
Réponse complète : les grandeurs physiques sont au nombre de n = 6 : (x1 , x2 , x3 , x4 , x5 , x6 ) = (x, t, c, V, D, k).
Ce sont respectivement l'espace x (caché dans le gradient et le laplacien), le temps t, un scalaire c qui peut
être la concentration (en mol.m−3 , dimension [N.L−3 ]), la vitesse V du uide (en m.s−1 , dimension [L.T−1 ]),
le coecient de diusion D (en m2 .s−1 , dimension [L2 .T−1 ]) et la constante cinétique de réaction k (en
s−1 , dimension [T−1 ]). Il y a r = 2 dimensions. On s'attend à voir apparaître p = n − r = 4 nombres
adimensionnés. Pour adimensionner l'équation, il convient de dénir r = 2 grandeurs de base. Arbitraire-
ment, choisissons une longueur caractéristique L0 et une vitesse caractéristique V0 . On s'assure qu'elles sont
dimensionnellement indépendantes. Pour chacune des grandeurs xi , on forme une échelle caractéristique Xi
de même dimension en utilisant ces deux grandeurs de base :
x → X1 ≡ L0
t → X2 ≡ L 0
V0
N
c → X3 ≡ L30
0
V → X4 ≡ V0
D → X5 ≡ L0 .V0
V0
k
→ X6 ≡ L0
Partant de là, on peut dénir autant de grandeurs adimensionnées (les variables sont désignées par une
astérisque tandis que les constantes sont désignées par la lettre N ) :
∗
x
≡ Xx1
t∗ ≡ Xt2
∗
≡ Xc3
c
V∗ ≡ V
X4
D
N5 ≡
X5
k
N6
≡ X6
Dans l'EDP, on remplace chacune des 6 grandeurs xi par le produit x∗i .Xi ou Ni .Xi . Elle devient :
∂c∗ .(1/L30 ) ∗ ∂c∗ .(1/L30 ) ∂ 2 c∗ .(1/L30 )
+ V .V 0 . = (N 5 .L0 .V 0 ). − N6 .(V0 /L0 ).c∗ .(1/L30 )
∂t∗ .(L0 /V0 ) ∂x∗ .L0 ∂(x∗ .L0 )2
En multipliant l'équation par L40 /V0 , on obtient l'équation adimensionnée que l'on cherchait :
D
∂c∗ N =
5
+ V ∗
.∇ ∗ ∗
c = N 5 .∇ ∗2 ∗
c − N 6 .c∗
avec L0 .V0
∂t∗ N = k
6 V0 /L0
Il y a trois phénomènes physiques : la convection (V.∇c), la diusion (D.∇2 c) et la réaction (k.c). L'EDP
montre que chacun de ces trois termes est homogène à ∂c/∂t. Par conséquent, on peut écrire :
[c] [c]
Convection : ≡ L
conv
[t] = [V ] . [x]
τ V0
0
L2
Diusion : [c]
[t] = [D] . [x][c]
2 d'où τdi mat ≡ D0
Réaction : [c]
≡ 1 chim
τ
= [k] . [c]
[t] k
15
Chapitre 1. Analyse dimensionnelle
Avec ces dénitions, les deux nombres adimensionnés N5 et N6 qui apparaissent dans l'EDP adimensionnée
peuvent s'écrire :
N = D τconv
5 L0 .V0 = τdi mat
N =
6
k
= τconv
V0 /L0 τchim
Avec les valeurs numériques données en entrée, le nombre N5 est de l'ordre de 3.10−7 , ce qui est très inférieur
à l'unité. Cela signie que le temps caractéristique de la convection est très court devant celui de la diusion.
Par conséquent, le comportement du système est imposé par la convection : une variation de vitesse du uide
aura une inuence directe tandis qu'une variation du coecient de diusion n'aura aucun impact notable.
Les transferts de matière par diusion étant extrêmement lents par rapport aux transferts convectifs, on
pourra négliger ce terme et s'attacher à résoudre l'équation suivante, plus simple :
∂c∗
+ V ∗ .∇∗ c∗ = −N6 .c∗
∂t∗
Exercice 7 Réponse courte : la puissance nécessaire est de 18 hp et le débit est de 27 ft3 .s−1 .
Réponse complète : on a vu plus haut que les trois nombres adimensionnés sont :
N ≡ ρ.ΩP3 .D5
1
Q
N5 ≡ Ω.D 3
N ≡ ∆H
6 D
16
Chapitre 2
17
Chapitre 2. Cinématique, forces et contraintes
Si l'on note ux la vitesse de déplacement du point de mire M , ce dernier passe de la position x à l'instant t
à la position x + ux .dt à l'instant t + dt. Par conséquent, la vitesse de variation est :
df f (x + ux .dt, t + dt) − f (x, t)
=
dt dt
Pour dt susamment petit, cette vitesse est égale :
df f (x + ux .dt, t + dt) − f (x, t + dt) f (x, t + dt) − f (x, t)
= +
dt dt dt
Le deuxième terme du membre de droite est égal à la dérivée partielle ∂t .
∂f
Concernant le premier, on peut
écrire le développement limité :
∂f
f (x + ux .dt, t + dt) = f (x, t + dt) + .ux .dt
∂x
En généralisant à trois dimensions, cela permet d'écrire :
df ∂f
M quelconque ⇒ = + (~u · ∇) f (2.1)
dt ∂t
Le dernier membre de cette expression est (~u · ∇) f qui traduit l'application de l'opérateur ~u · ∇ (produit
scalaire entre ~u et l'opérateur ∇) à la grandeur f . L'utilisation de parenthèses permet d'éviter la confusion
avec la grandeur ~u · ∇f qui peut être interprétée à tort comme le produit entre le vecteur ~u et la grandeur
∇f . En eet, cette expression n'a pas de sens mathématique (sauf si f est scalaire). La forme développée de
cet opérateur, dit opérateur d'advection est :
∂ ∂ ∂
~u · ∇ = ux . + uy . + uz . (2.2)
∂x ∂y ∂z
Cette dénition est valable quel que soit l'ordre tensoriel de f (scalaire, vecteur, tenseur d'ordre 2, . . .).
L'opérateur (~u · ∇) s'applique sur chacun des termes du tenseur f (par exemple ρ, vx ou bien σyz ). Les
18
2.2. Dérivées partielle, totale et particulaire
sections 7.7, 7.8 et 7.9 en page 92 et suivantes donnent l'expression développée de l'opérateur d'advection
dans les systèmes de coordonnées cartésiens, cylindriques et sphériques respectivement.
Cas particulier n°1 : le point de mire M est xe (~u = ~0). On est alors dans une description eulérienne
et la vitesse de variation de f est égale à la dérivée partielle, notée ∂f /∂t :
df ∂f
M xe ⇒ = (2.3)
dt ∂t
Cas particulier n°2 : le point de mire est attaché à une molécule de uide (~u = ~v ). On est alors dans une
description lagrangienne et la vitesse de variation de f est égale à ce qu'on appelle la dérivée matérielle ou
dérivée particulaire ou dérivée en suivant le mouvement, notée Df /Dt. Nous utiliserons préférentiellement
l'appellation dérivée matérielle.
df Df Df ∂f
M uide ⇒ = avec ≡ + (~v · ∇) f (2.4)
dt Dt Dt ∂t
Interprétation : l'équation 2.4 fait apparaître deux contributions dans la variation de la quantité f attachée
à une parcelle de uide qui se déplace :
une variation temporelle à position donnée
∂t . Lorsque le uide est au repos, le taux de variation de
∂f
19
Chapitre 2. Cinématique, forces et contraintes
Figure 2.4 Illustration de l'opérateur d'advection. A gauche, en noir, la trajectoire d'une parcelle de uide.
A droite, la température de cette même parcelle au cours du temps.
Exercice 8.
2 2
le champ de température dans un uide au repos est donné par T(x,y) = T . x R+y 2 . On
utilise un capteur infra-rouge pour mesurer à distance la température en un point M particulier à
la surface du uide. Par l'intuition, prédire les variations puis trouver l'expression de la fonction
TM (t) mesurée par le capteur lorsque la position (xM , yM ) de la mire est donnée par xM (t) = A.t2 et
yM (t) = 1 − xM (t). En déduire l'expression de TM (~rM ).
20
2.3. Représentations graphiques de l'écoulement
vitesse ~v = (vx , vy , vz ), il sut d'écrire l'équation suivante puis de réaliser une intégration :
dxP
dt = vx (xP , yP , zP , t)
d~rP
= ~v(~rP ) soit dyP
dt = vy (xP , yP , zP , t) (2.5)
dt
dzP = v (x , y , z , t)
dt z P P P
La photo de gauche sur la gure 2.1 a été obtenue avec un long temps d'exposition. Elle montre diérentes
trajectoires rendues visibles par la présence de poussières incandescentes qui sont emportées par les gaz
chauds dégagés par un feu de camp.
La résolution se fait par intégration, après avoir rendu explicite la dépendance de vx , vy et vz vis-à-vis de
x, y, z . Pour dénir les trois constantes d'intégration, il sut d'identier un point particulier (x0 , y0 , z0 ) par
lequel passe la ligne de courant. En deux dimensions, puisque dz = 0 et vz = 0, seule la première équation
scalaire est intéressante. Lorsque la ligne de courant peut être décrite par une fonction de la forme y(x),
alors on écrit simplement :
vy (x, y)
y 0 (x) = (2.7)
vx (x, y)
Un tube de courant est un ensemble de lignes de courant s'appuyant sur un contour fermé (voir schéma
de droite sur la gure 2.2). Par construction, en chaque point de la surface du tube, la vitesse est tangente
à la surface.
21
Chapitre 2. Cinématique, forces et contraintes
On notera que, dans le cas général, les lignes de courant peuvent se calculer mais ne peuvent pas s'observer
expérimentalement. Dans le cas particulier d'un écoulement stationnaire (et seulement dans ce cas), elles se
confondent avec les lignes d'émission et les trajectoires qui, elles, sont observables expérimentalement.
Table 2.2 Représentations eulériennes : vecteurs (gauche), lignes de courant (centre) et tube de courant
(droite).
Exercice 11. un écoulement est caractérisé par le champ de vitesse suivant : ~v = Axt.~ex − Byt2 .~ey
avec A = 0, 2 s et B = 0, 1 s−3 . Calculer la vitesse du uide en (2, 8) à l'instant t = 2 puis trouver
−2
l'équation de la ligne de courant passant par le point (2, 8). Établir l'expression donnant la trajectoire
de la parcelle de uide située en (2, 8) à l'instant t = 2.
Exercice 12. en utilisant le Système International, donner l'équation aux dimensions du produit
ρ.V 2 et du produit ρ.ν 2 . Commenter.
Exercice 13. un écoulement est caractérisé par le champ de vitesse ~v(x,y,t) = Ax.~ex − Ay.~ey et le
2 2
champ de température T(x,y,t) = T . x R+y2 . Calculer l'évolution temporelle de la température de la
parcelle de uide présente en (0, R) à t = 0 puis calculer l'accélération qu'elle subit.
22
2.4. Décomposition du champ de vitesse
Exercice 14. calculer la trajectoire d'une parcelle de uide située initialement en (x0 , y0 ) dans un
écoulement dont le champ de vitesse est : ~v = −ω.y.~ex + ω.x.~ey avec ω un scalaire positif.
Exercice 15. on considère un champ de vitesse instationnaire ~v(x,y,t) = V. sin (ω.t) .~ex + V.~ey avec
V et ω deux scalaires positifs. Dans cet écoulement, une cheminée d'usine positionnée en (xE , yE )
relargue de la fumée entre les instants ti et tf . Calculer l'équation de la ligne d'émission à un instant
t > tf .
Ecrivons un développement limité de la vitesse au 1er ordre. L'opérateur gradient sera noté ∇ (on trouve
aussi la notation "grad") :
∂vx ∂vx ∂vx
∂x ∂y ∂z
~v (~x + d~x) = ~v (~x) + (∇~v ) .d~x + Résidu avec ∇~v ≡ ∂vy ∂vy ∂vy
(2.8)
∂x ∂y ∂z
∂vz ∂vz ∂vz
∂x ∂y ∂z
Note : le gradient de vitesse s'exprime parfois en faisant intervenir le produit dyadique. On a alors la relation
T
∇~v = grad ~v = (∇ ⊗ ~v ) . Le gradient de vitesse ∇~v est un tenseur d'ordre 2 qu'il est possible d'exprimer
sous la forme d'une somme d'un tenseur symétrique d et d'un tenseur anti-symétrique r :
h i h i
T ∂vj
d ≡ 1
2. ∇~v + (∇~v ) soit dij ≡ 12 . ∂vi
∂xj + ∂xi
∇~v = d + r avec (2.9)
h i h i
T ∂vj
≡ 1
2 . ∇~
v − (∇~v ) soit rij ≡ 21 . ∂vi
−
r
∂xj ∂xi
Pour des raisons que nous allons préciser ci-après, le tenseur d est appelé tenseur des taux de défor-
mations et le tenseur r le tenseur des taux de rotation. La raison du choix de cette décomposition va
apparaître en analysant la gure 2.6.
23
Chapitre 2. Cinématique, forces et contraintes
2.4.1 Rotation
La gure 2.6 fait apparaître les angles α, β, γxy qui sont tous les trois négatifs (attention au sens trigonomé-
trique adopté). L'angle entre l'axe des abscisses et la bissectrice de l'angle BAC
\ ou bac c sera nommé θ. Il est
positif par construction. Ces conventions se traduisent notamment par les égalités :
α+γ +β −π
xy = 2
(2.10)
γxy
θ = −α − 2
On dénit maintenant le vecteur tourbillon (ou simplement tourbillon) par la relation suivante :
∂vz ∂vy
∂y − ∂z
~ ≡ 1 .rot ~v avec
Ω rot ~v ≡ ∂vx
− ∂vz (2.11)
2 ∂z ∂x
∂vy ∂vx
∂x − ∂y
Le sens physique de cette grandeur est simple : la norme du vecteur tourbillon est la vitesse angulaire θ̇.
En d'autres termes, le vecteur tourbillon permet de quantier la vitesse à laquelle la parcelle uide
h tournei
autour de son centre de gravité. En eet, en utilisant la relation 2.10, on peut écrire que θ̇ = 2 . β̇ − α̇ .
1
∂v ∂v
Par ailleurs, on a sin(dα) ≈ dα et sin(dα) ≈ − ∂xy .dt d'où α̇ ≈ − ∂xy et, par analogie β̇ ≈ − ∂v
∂y . Avec ces
x
éléments, on trouve facilement que Ω~ = θ̇.~ez . Cela signie qu'une valeur positive de Ωz est associée à une
rotation dans le sens trigonométrique.
On note que le vecteur tourbillon ainsi déni est lié au tenseur r des taux de déformation en rotation par la
24
2.4. Décomposition du champ de vitesse
relation :
rzy
~ =
Ω (2.12)
rxz
ryx
On parle d'écoulement tourbillonaire lorsqu'il existe au moins un endroit où le vecteur tourbillon n'est
pas nul. A contrario, on parle d'écoulement non tourbillonaire lorsque le vecteur tourbillon est nul
partout dans le volume étudié.
Compte-tenu de la dénition du tenseur r (équation 2.9) et du vecteur Ω
~ (équation 2.11), on obtient :
~ ∧ d~x
r.d~x = Ω (2.13)
Le calcul de la vorticité revêt une grande importance pratique car c'est souvent dans les tourbillons que
l'énergie se dissipe le plus. Le design des appareils peut alors être optimisé.
Table 2.3 Par imitation avec un faucon en vol à faible vitesse, ajout d'un aileron pour contrôler la
séparation des tourbillons et augmenter la portance d'un facteur 2. La couleur représente ωz (la vorticité
selon z ) en rouge pour une vorticité positive (sens trigonométrique) et en bleu pour une vorticité négative.
Source : web.
25
Chapitre 2. Cinématique, forces et contraintes
2.4.2 Déformations
En mécanique des uides, les taux de déformations sont usuellement caractérisés par le tenseur γ̇ plutôt que
le tenseur des taux de déformation d déni par l'équation 2.9 en page 23 :
∂vi ∂vj
γ̇ ≡ 2.d et donc γ̇ij = + (2.15)
∂xj ∂xi
Rigoureusement parlant, le tenseur γ̇ est un tenseur décrivant les taux de déformation. Toutefois, on utilise
souvent le terme tenseur des taux de cisaillement. Comme on va le voir par la suite, cette dénomination
est abusive car ces deux grandeurs ne sont identiques que lorsque les uides sont incompressibles.
Quoiqu'il en soit, il est naturel de décrire la déformation totale (d) comme la somme d'une expansion isotrope
(dexp ) et d'un cisaillement à volume constant (dcis ) :
Expansion. La gure 2.6 montre qu'au cours du laps de temps dt, le volume de matière V ≡ dx.dy.dz
évolue de la quantité suivante :
∂vx ∂vy ∂vz
dV = .dx.dt .dy.dz + .dy.dt .dx.dz + .dz.dt .dx.dy
∂x ∂y ∂z
La variation relative de volume s'écrit :
dV ∂vx ∂vy ∂vz dV 1 dV
= + + .dt = div ~v .dt soit = div ~v .dt ou encore div ~v = .
V ∂x ∂y ∂z V V dt
On retiendra la formulation suivante, qui montre qu'une divergence de la vitesse positive correspond à une
diminution locale de masse volumique :
−1 dρ
div ~v = . (2.17)
ρ dt
Naturellement, un uide incompressible est caractérisé par div ~v = 0.
La composante du tenseur des taux de déformations d associée à l'expansion isotrope sera notée dexp et,
logiquement, proportionnelle à la divergence de la vitesse. Puisqu'on travaille en trois dimensions, on aura :
Trace(d) div ~v
dexp ≡ .I soit dexp,ij ≡ .δij (2.18)
3 3
div ~v
1 ∂vi ∂vj
dcis ≡ d − dexp soit dcis,ij ≡ . + − .δij (2.19)
2 ∂xj ∂xi 3
Puisque la trace de d et la trace de dexp sont toutes deux égales à la divergence de la vitesse, le tenseur dcis
est de trace nulle. Comme d, le tenseur dcis est symétrique. Son sens physique est associé aux déformations
angulaires (cisaillement à volume constant). Le choix du symbole grec γ̇ pour désigner le tenseur des taux de
cisaillement (équation 2.15) n'est pas le fruit du hasard : sa composante γ̇xy se trouve être égale à la vitesse
de variation de l'angle γxy déni sur la gure 2.6. En eet, la deuxième équation 2.10 permet d'écrire :
d (α + β)
γ̇xy = −
dt
∂v
Comme on l'a déjà montré dans la section 2.4.1, on a dα dt = − ∂x et dt = − ∂y . Cela nous permet bien de
y dβ ∂vx
retrouver que le taux de cisaillement γ̇xy n'est rien d'autre que la vitesse à laquelle l'angle γxy évolue, c'est-
à-dire la vitesse avec laquelle la parcelle carrée de uide se transforme en losange sous l'eet des gradients
26
2.4. Décomposition du champ de vitesse
de vitesse.
Une valeur positive du taux de cisaillement se traduit par un angle γ dont la valeur absolue devient de plus
en plus faible. Une valeur élevée de γ̇ est associée à des déformations élevées. C'est pour cette raison que
le tenseur dcis est associé aux déformations angulaires. Le tableau 2.5 donne l'ordre de grandeur de γ̇ pour
quelques procédés ou instruments.
Figure 2.7 Schéma d'une extrudeuse bi-vis pour granulation de polymère. Le taux de cisaillement γ̇ est
très élevé (entre 10 et 105 s−1 ). Source : web.
2.4.3 Synthèse
Dans ce qui précède, on a décomposé la déformation d'une parcelle de uide en diérentes composantes.
On peut maintenant réécrire l'équation 2.8 sous la forme :
div ~v
~v (~x + d~x) = ~v (~x) + ~ d~x
| ∧
Ω + .d~x + dcis .d~x + Résidu (2.20)
| 3{z }
{z }
Rotation
| {z }
Cisaillement
| {z }
Translation Expansion
Exercice 16. On considère un écoulement incompressible dans le plan rθ. La composante radiale de
la vitesse est de la forme vr = V. cos θ. Identier les expressions possibles de la composante tangentielle
vθ .
27
Chapitre 2. Cinématique, forces et contraintes
Figure 2.8 Schéma d'une atomisation par spray-drying. Une suspension (particules dans liquide) est
pulvérisée par une buse et pénètre dans une tour à haute température où le solvant s'évapore. Le taux de
cisaillement γ̇ est très élevé (entre 103 et 106 s−1 ). Source : web.
Exercice 18. On considère l'écoulement stationnaire d'un uide newtonien de masse volumique ρ et
de viscosité µ entre deux plans parallèles xes situés en y = ±H , entre la section x = 0 et x = L.
La pression en (x, y, z) = (0, 0, 0) est notée p1 et celle en (L, 0, 0) est notée p2 . La gravité est notée
~g = −g.~ey . On peut montrer que le prol de vitesse est parabolique :
1 p1 − p2
~v = vx .~ex avec . H 2 − y2
vx (x, y, z, t) = .
2µ L
Décomposer le tenseur des gradients de vitesse en ses composantes élémentaires.
28
2.4. Décomposition du champ de vitesse
Figure 2.9 Illustration des lignes de courant, champs de vitesse, pression, taux de cisaillement et vorticité
pour un écoulement 2D convergent.
29
Chapitre 2. Cinématique, forces et contraintes
dF~ext/int
T~ ≡ lim (2.21)
dS→0 dS
30
2.5. Forces et contraintes
De la sorte, on peut calculer le vecteur-contrainte T~ pour une orientation ~nout quelconque par :
T~ = σ · ~nout (2.24)
et en déduire la force F~ext/int qu'exerce l'environnement extérieur sur l'intérieur du volume de contrôle à
travers la surface d'aire S (et de normale unitaire sortante ~nout ) :
x σ · ~n x Xσ 3
F~ext/int = out .dS d'où Fext/int,i = ij .nout,j .dS (2.25)
j=1
Avec les choix qui sont les nôtres 4 , nous retiendrons la convention suivante, illustrée sur la gure 2.10 :
Le terme σij .dS est la composante selon la direction ~ei de la force dF~ qu'exerce l'environnement
extérieur sur un volume de matière à travers l'élément de surface d'aire dS et de normale sortante de
direction ~ej .
Considérons maintenant un volume de matière délimité par une surface dont la normale sortante est de
direction ~ex . Lorsque σxx > 0, alors l'environnement exerce sur le volume de matière une force dont la
composante selon ~ex est positive. Cela signie que l'environnement "tire" la matière qui se trouve "en
traction". Inversement, lorsque σxx < 0, l'environnement "pousse" la matière qui se trouve "en compression".
Dans ce qui suit, on va voir que l'expression de la force F~ext/int contient deux contributions : une force
normale et une force tangentielle, dite de cisaillement.
Figure 2.10 Notations utilisées pour le tenseur des contraintes : σij est la composante selon l'axe i de la
force exercée sur la face de normale j .
4. D'autres choix sont possibles. Ils sont exposés dans la section 7.4 en page 83.
31
Chapitre 2. Cinématique, forces et contraintes
On voit que I1 = Trace (σ) et que I3 = Det (σ). Les dimensions de I1 , I2 et I3 sont respectivement Pa, Pa2
et Pa3 .
La pression hydrostatique p est une grandeur particulièrement utile en mécanique des uides. Elle
est construite à partir du premier invariant du tenseur des contraintes (sa trace), ce qui rend sa valeur
indépendante du choix du repère. Plus précisément, compte-tenu de la convention "normale sortante" adoptée
dans ce document, sa dénition 5 s'énonce :
−Trace (σ)
p≡ (2.27)
3
Un matériau en compression est caractérisé par une valeur positive de la pression, identique à la pression
dénie en thermodynamique. Notons également qu'un uide est toujours en compression : contrairement à
un solide qui peut résister à la traction, un uide ne le peut pas.
Le tenseur des contraintes a été décomposé en σ = τ − p.I . Il s'ensuit que l'équation 2.25 qui donnait
l'expression de la force F~ext/int s'exerçant sur un élément de surface dS peut être ré-écrite comme la somme
d'une force tangentielle (due au cisaillement) et d'une force normale (due à la pression) :
x τ · ~n x Xτ 3
F~ext/int = out .dS − p.S.~nout d'où Fext/int,i = ij .nout,j .dS − p.S.nout,i (2.29)
j=1
Le tenseur τ des contraintes de cisaillement possède lui aussi trois invariants, notés J1 , J2 , J3 . Ils sont dénis
de la même manière que pour le tenseur des contraintes, comme coecients du polynôme caractéristique.
J ≡ τxx + τyy + τzz
1
J2 ≡ τxx .τyy + τyy .τzz + τxx .τzz − τxy 2 2
− τyz 2
− τzx (2.30)
J ≡ τ .τ .τ + 2.τ .τ .τ − τ 2 .τ − τ 2 .τ − τ 2 .τ
3 xx yy zz xy yz zx xy zz yz xx zx yy
5. Si l'on avait adopté la convention "normale entrante" (voir discussion en page 83), il aurait fallu dénir la pression par
p≡
Tr(σ) pour garantir sa positivité.
3
6. Si l'on avait adopté la convention "normale entrante", il aurait fallu écrire σ ≡ τ + p.I .
32
2.5. Forces et contraintes
Par construction, le premier invariant J1 est nul. Cette condition particulière, et le fait que le tenseur des
contraintes de cisaillement τ est symétrique (démonstration en page 87), permet d'exprimer les invariants
d'une autre manière :
J1 = 0
h i
2 2 2
J2 = 13 .I12 − I2 = 16 . (σxx − σyy ) + (σyy − σzz ) + (σzz − σxx ) + σxy 2 2
+ σyz 2
+ σzx (2.31)
J = 2 .I 3 − 1 .I .I + I
3 27 1 3 1 2 3
Il est parfois utile de résumer le tenseur des contraintes par un scalaire unique, que l'on nomme contrainte
eective. Diérentes propositions ont été formulées : Tresca, Rankine et von Mises. Ci-dessous l'expression
de la contrainte eective de von Mises ; c'est probablement la plus utilisée en mécanique des uides.
vM ≡ 3.J
p
σe 2
Exercice 19. On considère trois situations où le tenseur des contraintes σ est connu. Que peut-on
dire dans chacun des cas ?
−2 0 0 1 0 1 1 0 4
σa ≡ 0 −2 0 σb ≡ 0 2 0 σc ≡
0 0 0
0 0 −2 −1 0 3 4 0 −4
Exercice 20. On considère l'écoulement stationnaire d'un uide newtonien de masse volumique ρ et
de viscosité µ entre deux plans parallèles xes situés en y = ±H , entre la section x = 0 et x = L.
La pression en (x, y, z) = (0, 0, 0) est notée p1 et celle en (L, 0, 0) est notée p2 . La gravité est notée
~g = −g.~ey . On peut montrer que le prol de vitesse est parabolique et que le tenseur des contraintes
se met sous la forme :
−p τxy 0
τ
xy = −y. p1 −p 2
Calculer la force exercée par le uide sur la paroi supérieure. Commenter les diérents termes.
33
2.6 Solutions des exercices
Exercice 8 La mire se déplace le long de la droite y = 1 − x depuis le point (0, 1) vers le bas. Le champ de
température est axisymétrique autour de l'axe [Oz) perpendiculaire au plan (Oxy), présente un minimum
en (0, 0) et croît quand on s'éloigne de l'origine. Qualitativement, on peut s'attendre à ce que la fonction
TM (t) soit décroissante puis croissante. Pour les calculs, on va utiliser la dénition de la dérivée totale
(équation 2.1) qui s'énonce dT ∂T
uM · ∇) T . La vitesse du point de mire s'obtient en dérivant
dt M = ∂t + (~
sa position par rapport au temps. On trouve une vitesse de mire croissante ~uM = 2A.t. (~ex − ~ey ). Comme
le système étudié est stationnaire (vitesse du uide et champ de température indépendants du temps),
on écrit que la dérivée partielle de la température par rapport au temps hest nulle. En utilisant
i la formule
développée de l'opérateur d'advection (équation 2.2), on a : dT ∂T ∂T
soit dT
dt M = 2A.t. ∂x M − ∂y dt M =
M
R2 .t.(x(t) − y(t) ) ou encore dt M = T . R2 .t.(2A.t − 1). Par intégration par rapport au temps, on obtient
T . 4A dT 4A
2
TM (t) = T0 + T . 2A R2 .(At − t ). La constante d'intégration T0 se calcule en disant que pour t = 0, la mire est
4 2
en (0, 1) où la température est connue. On en déduit que T0 = T /R2 . La fonction ainsi trouvée se comporte
qualitativement comme intuité.
Exercice 9 La mire se déplace le long d'un cercle centré à l'origine. Le champ de température est axi-
symétrique autour de l'axe [Oz) perpendiculaire au plan (Oxy), présente un minimum en (0, 0) et croît
quand on s'éloigne de l'origine. Qualitativement, on peut s'attendre à ce que la fonction TM (t) soit une
constante. Pour les calculs, on va utiliser la dénition de la dérivée totale (équation 2.1) qui s'énonce
dT ∂T
uM · ∇) T . La vitesse du point de mire s'obtient en dérivant sa position par rapport au
dt M = ∂t + (~
temps. On trouve une vitesse de mire croissante ~uM = − sin(t).~ex + cos(t).~ey . Comme le système étudié est
stationnaire (vitesse du uide et champ de température indépendants du temps), on écrit que la dérivée
partielle de la température par rapportau temps est nulle. En utilisant la formule développée de l'opérateur
d'advection (équation 2.2), on a : dT dt M = ... = 0. Par intégration par rapport au temps, on obtient que
TM (t) est une constante, conformément à ce qu'on a pu intuiter.
Exercice 10 Dans l'approche lagrangienne, l'accélération ~aP d'une parcelle de uide P se déduit de sa
2
position ~rP par dérivation : ~aP = ddt~r2P . Ici, la vitesse ~uP de la parcelle est égale à la vitesse (connue) ~v
du uide. Par dénition, on a dx dt = uP x soit dt = V. 1 + L . Cette équation diérentielle s'intègre
dxP xP
P
facilement et donne xP (t) = L. eV t/L − 1 . Par double dérivation, on obtient la composante selon x de
2 2
l'accélération de la parcelle de uide : aP x = ddtx2P = VL .eV t/L .
Dans l'approche eulérienne, l'accélération ~aP d'une parcelle de uide est égale à la dérivée matérielle du champ
de vitesse ~v . On a donc ~aP = D~
Dt . D'après la dénition donnée dans l'équation 2.4 et en ne s'intéressant qu'à
v
V 2 V t/L
la composante selon x, on peut écrire aP x = ∂v ∂t + vx . ∂x + vy . ∂y + vz . ∂z soit aP x = vx . ∂x = L .e
x ∂vx ∂vx ∂vx ∂vx
.
Exercice 11 à l'instant t = 2, la vitesse du uide en (2, 8) est ~v = 0, 8.~ex − 3, 2.~ey . Sa norme est d'environ
3,3 m.s−1 . Pour établir l'expression des lignes de courant, on utilise l'équation 2.7 pour trouver dy −Bt dx
y = A . x
−Bt
de quoi on tire y = C.x A (la grandeur t est considérée comme constante). La constante d'intégration C
est donnée par la condition initiale ; on obtient C = 16. A chaque instant t, la ligne de courant passant
−t
par le point (2, 8) a pour équation y(x) = 16.x 2 . Pour calculer la trajectoire d'une parcelle de uide, on
utilise l'équation 2.5 qui se décline sous la forme de dx
dt = A.xP .t et dt = −B.yP .t . Après séparation des
P dyP 2
Exercice 12 l'équation aux dimensions du produit ρ.V 2 est [ρ.V 2 ] = M.L−1 .T −2 . On constate que ce
dernier produit est homogène à une pression. Celle du produit ρ.ν 2 est [ρ.ν 2 ] = M.L.T −2 . On constate que
ce dernier produit est homogène à une force.
2.6. Solutions des exercices
Exercice 13 dans un premier temps, on va calculer la trajectoire de cette particule, c'est-à-dire calculer sa
position ~rP = (xP , yP ) au cours du temps. On utilise pour cela l'équation 2.5 qui se décline sous la forme
dt = A.xP et dt = −A.yP . Après séparation des variables, on peut intégrer par rapport au temps et
de dx P dyP
Ce résultat aurait pu être obtenu en utilisant l'équation 2.4 pour calculer le champ d'accélération (partout
dans l'espace) à partir du champ de vitesse :
∂vx
∂~v vx . ∂v ∂vx
∂x + vy . ∂y
x
A2 .x
~a = + (~v · ∇) ~v d'où ~a = ∂t + =
∂t ∂vy
vx .
∂vy
+ vy .
∂vy
A2 .y
∂t ∂x ∂y
Exercice 15 La ligne d'émission à un instant t donné est une courbe nie : l'extrémité la plus en aval (resp.
amont) est occupée par la parcelle de uide qui était en E à l'instant ti (resp. tf ). Considérons un point P
de coordonnées (xP , yP ) appartenant à la ligne d'émission cherchée. Par dénition, la parcelle de uide qui
se trouve en P à l'instant t était en E à l'instant t0 .
35
Chapitre 2. Cinématique, forces et contraintes
Calculons sa trajectoire entre t0 et t. On utilise l'équation 2.5 et une intégration temporelle pour obtenir :
x (t) = −V . cos (ω.t) + C
P ω 1
y (t) = V.t + C
P 2
Comme on sait que la parcelle était en E à l'instant t0 , on peut déterminer les deux constantes d'intégration
C1 et C2 . On obtient :
x (t) = x + V . cos (ω.t0 ) − V . cos (ω.t)
P E ω ω
y (t) = y + V. (t − t0 )
P E
La ligne d'émission recherchée est une portion de la courbe xP (yP ) ainsi trouvée.
(vz = 0), on a ∂v ∂θ = −V. cos θ . On en déduit que la composante radiale de la vitesse doit être de la forme
θ
Exercice 17 La vitesse n'ayant qu'une composante selon x, il est facile de trouver la position des points abcd
à l'instant t0 . Pour i ∈ {O, a, b, c, d}, on écrit x0i = xi + U. yHi .t0 . On trouve respectivement x0O = 5, x0a = 4,
x0b = 6, 5, x0c = 6 et x0d = 3, 5 mm. En traçant les nouveaux points, on peut déduire que la déformation est
causée à la fois par une déformation angulaire (car cOb d 6= c\ 0 O 0 b0 ) et par une rotation (car l'angle entre la
notre cas particulier, la composante verticale vy est identiquement nulle et, puisque ∂v ∂y = H , le taux de
x U
La dénition du vecteur tourbillon (équation 2.11) est rappelée Ω ~ ≡ 1 .rot ~v . Puisque vx est la seule compo-
2
sante non nulle et ne dépend que de y , on a Ω ~ = Ωz .~ez avec Ωz = − 1 . ∂vx = − 1 . U = −0, 5 s−1 . Le signe
2 ∂y 2 H
négatif indique que la rotation se fait dans le sens anti-trigonométrique (=sens horaire), ce que conrme la
lecture du schéma.
36
2.6. Solutions des exercices
D'où :
0 y 0 0 y 0
−1 p1 − p2 −1 p1 − p2
∇~v = d + r avec d= . . y et r= . . −y
0 0 0 0
2µ L 2µ L
0 0 0 0 0 0
On voit que la composante selon z du vecteur tourbillon est positive pour 0 < y < H . Cela signie que
la rotation dans cette région se fait dans le sens trigonométrique (l'angle négatif θ déni sur la gure 2.6
augmente au cours du temps / sa valeur absolue diminue).
Concernant l'expansion, la relation 2.18 permet d'écrire :
0 0 0
div ~v
1 ∂vx ∂vy ∂vz
dexp ≡ .I = . + + .I d'où dexp = 0 0 0
3 3 ∂x ∂y ∂z
0 0 0
On voit que le uide ne subit pas d'expansion. Cela signie qu'il est incompressible.
Concernant le cisaillement, la dénition 2.19 permet d'écrire que, dans notre cas particulier, le tenseur dcis
est égal au tenseur d.
Exercice 19 Cas a : les contraintes de cisaillement (termes non diagonaux) sont tous nuls. Les gradients de
vitesse étant nuls, le uide est au repos (ou en translation uniforme - écoulement piston). La pression p est
égale à l'opposé du tiers de la trace (équation 2.27), soit 2.
Cas b : ce cas ne peut pas exister car le tenseur des contraintes est symétrique pour des uides communs
(équation 7.38 en page 87).
Cas c : la pression est égale à 1. L'écoulement se fait dans le plan xz car les gradients de vitesse selon y sont
nuls.
Exercice 20 La force F~f luide/paroi que le uide exerce sur la paroi peut se calculer grâce à l'équation 2.25
Le vecteur normal unitaire ~nout est dirigé vers l'extérieur de la zone uide et donc s'écrit ~nout = ~ey pour la
face supérieure. En considérant une aire S = 1 m2 , on obtient :
−p τxy 0 0 −τxy
F~f luide/paroi = −σ · ~ey .S = − τxy −p 0 · 1 .S = p .S
0 0 −p 0 0
37
Chapitre 2. Cinématique, forces et contraintes
38
Chapitre 3
La deuxième section présente les propriétés de transfert des matériaux. Dans la notion de transfert se
cache la notion de cinétique (et non plus de thermodynamiques). En eet, ces propriétés caractérisent la
vitesse à laquelle un matériau réagit à une sollicitation (exemple : intensité du ux de chaleur pour une
diérence de température imposée).
La troisième section propose une typologie sommaire des diérents types d'écoulement rencontrés dans
la nature ou en situation industrielle.
39
Chapitre 3. Fluides particuliers et lois constitutives, classication des écoulements
On parle de uide incompressible lorsque sa masse volumique ρ peut être considérée comme quasiment
constante (au maximum quelques pourcents de variation sur la plage de pression et température caractéri-
sant l'écoulement). En général, la nature compressible/incompressible du uide a moins d'importance que
la nature compressible/incompressible de l'écoulement (voir page 47).
ρ (T, p) = ρ0 . [1 − α0 . (T − T0 )] (3.3)
où α0 est le coecient de dilatation thermique déni dans l'équation 3.1 et calculé pour la température de
référence T0 . On utilise cette expression de la masse volumique dans le bilan de quantité de mouvement
4.24 qui prend alors le nom d'équation de Boussinesq. Pour un gaz parfait 2 , on peut montrer que le
coecient de dilatation thermique α est égal à 1/T (soit environ 3,4.10−3 K−1 à 20°C). Pour l'eau, il est
égal à 2,07.10−4 K−1 à 20°C.
ρ (T, p) = ρ0 . [1 + κT . (p − p0 )] (3.4)
La thermodynamique nous apprend qu'on peut faire apparaître l'enthalpie massique h̄ et l'énergie interne
massique ū :
∂ h̄
c̄p = ∂T
p
(3.6)
∂ ū
c̄V =
∂T V
40
3.2. Equations constitutives
Le coecient adiabatique γ vaut 5/3 pour un gaz parfait monoatomique et 7/5 pour un gaz parfait diato-
mique. Pour une substance absolument quelconque, on peut montrer la relation :
κT
γ=
κS
La vitesse du son est faiblement inuencée par la pression mais peut varier signicativement en fonction de
la température. Pour l'eau liquide à (20°C, 1 bar), la vitesse du son est de 1482 m.s−1 . Elle est égale à 1484
m.s−1 pour (20°C, 10 bars) et 1526 m.s−1 pour (40°C, 1 bar). Pour un gaz parfait de masse molaire M̃ ,
l'équation d'état pV = N RT permet d'établir que :
R
q
vs = γ.R̃.T avec R̃ ≡ (3.8)
M̃
Pour de l'air (proche d'un gaz parfait diatomique), la vitesse du son est égale à 343,4 m.s−1 pour (20°C, 1
bar), à 344,5 m.s−1 pour (20°C, 10 bar) et à 355 m.s−1 pour (40°C, 1 bar).
Lorsque les champs de température, de vitesse et de concentration sont uniformes dans un système (i.e.
même valeur en chaque point de l'espace), alors ce dernier est au repos et son état ne varie pas au cours
du temps. En revanche, s'il existe des inhomogénéités spatiales, alors on observe l'apparition de ux qui
tendent à les réduire. Par exemple, un gradient de température induira un ux de chaleur, un gradient de
concentration induira un ux de quantité de matière, . . . Ainsi, les lois constitutives établissent un lien entre
une force motrice et un ux pour un matériau donné. Pour chaque loi constitutive, il y a, plus ou moins
caché, un couple de grandeurs conjuguées.
Pour de "faibles" inhomogénéités, les ux sont diusifs et on se trouve dans le domaine linéaire, c'est-
à-dire que le ux diusif d'une grandeur extensive est proportionnel à la force motrice (le gradient de sa
grandeur conjuguées). Les relations force-ux ont été identiées par les grands noms de la Science dans des
congurations simples. Si les homogénéités dépassent un certain seuil, il peut y avoir des ux convectifs
également.
41
Chapitre 3. Fluides particuliers et lois constitutives, classication des écoulements
de proportionnalité est notée λ et s'appelle conductivité thermique. Elle s'exprime en W.m−1 .K−1 . Les
gures 3.1 et 3.2 donnent la valeur de la conductivité thermique pour diérents matériaux courants.
42
3.2. Equations constitutives
sique (ou la concentration) des constituants, plus facile à mesurer ou calculer que le potentiel chimique.
Plusieurs mécanismes diusifs ont été identiés : la diusion de concentration (également appelée diu-
sion de masse ou diusion ordinaire ), la diusion thermique qui résulte d'un gradient de température, la
diusion de pression qui résulte d'un gradient de pression et la diusion forcée qui est causée par des forces
extérieures s'appliquant inégalement sur chaque constituant (exemple : force électrostatique). Dans les cas
usuels, la diusion de concentration est prépondérante.
Dans le cas particulier d'un système binaire A-B, la loi de Fick dresse un lien de proportionnalité entre
le ux diusif de matière et le gradient de composition. La constante de proportionnalité est notée DAB et
s'appelle coecient de diusion moléculaire. Elle s'exprime en m2 .s−1 . Dans sa variante massique, la
loi de Fick s'écrit :
Loi de Fick (massique) : ϕ~ di,A = −ρ.DAB .∇ωA (3.10)
où ϕ
~ di,A est le ux diusif massique de A (en kg.m−2 .s−1 ) et ωA sa fraction massique (voir section 4.3 en
page 54). Dans sa variante molaire, la loi de Fick s'écrit :
où ϕ
~ ∗di,A représente le ux diusif molaire de A (en mol.m−2 .s−1 ) et xA sa fraction molaire (voir section 4.2
en page 53).
87)
par construction, le tenseur des contraintes de cisaillement τ possède les mêmes symétries que σ
43
Chapitre 3. Fluides particuliers et lois constitutives, classication des écoulements
Le tenseur des contraintes de cisaillement peut également être vu comme la somme de ces deux contributions :
τ = τcis + τexp
On appelle uide newtonien un uide caractérisé par une simple relation de proportionnalité entre τcis et
dcis d'une part et entre τexp et dexp d'autre part :
τ ≡ 2µ.dcis
cis
(3.13)
≡ 3µV .dexp
τ
exp
Pour un uide subissant un cisaillement pur dans une cellule de Couette plane (écoulement entre deux plans
parallèles), le champ de vitesse est de la forme ~v = vx .~ex avec vx (x, y, z) = V0 . H
y
. Pour cet écoulement, que le
uide soit compressible ou incompressible, la divergence de la vitesse est nulle et l'on peut écrire τxy = µ. ∂v ∂vy .
x
Cela montre que la grandeur µ n'est rien d'autre que la célèbre viscosité dynamique, et s'exprime en Pa.s.
Deux autres unités se rencontrent : le Poiseuille (symbole PI) et la Poise (symbole P ou Po). Les équivalences
sont les suivantes : 1 Pa.s = 1 PI = 10 P = 1000 cP.
Pour un uide compressible subissant une compression pure, le champ de vitesse peut se mettre sous la
forme ~v = vx .~ex avec vx (x, y, z) = V0 . L
x
. Dans ce cas, on trouve que le tenseur des taux de cisaillement γ̇
est nul et l'on peut écrire τ = µ .div ~v .I avec µ0 ≡ µV − 23 µ. Le paramètre µV apparaît donc étroitement
0
On vient d'établir la forme particulière de la relation τ = f (γ̇) pour un uide newtonien. On voit que, pour
un uide newtonien et incompressible (div ~v = 0), on a une relation simple de proportionnalité entre la force
motrice (le taux de cisaillement γ̇ ) et la densité de ux de quantité de quantité de mouvement. Il s'agit de
la loi de Newton :
Loi de Newton : τ = µ.γ̇ (3.16)
Cette expression fait écho à la loi de Fourier et à la loi de Fick. Contrairement à ces deux dernières, on
notera l'absence de signe "-" dans la loi de Newton. C'est à mettre sur le compte de la convention choisie
sur l'orientation de la normale : elle est sortante dans ce document (voir tableau 7.3).
44
3.3. Classication des écoulements
Un écoulement externe est autour d'un objet immergé dans un uide sans limite (sous-marin en plongée,
avion en vol, balle de golf, fusée,. . .).
45
Chapitre 3. Fluides particuliers et lois constitutives, classication des écoulements
A contrario, il sera qualié d'instationnaire s'il existe un ou plusieurs endroits où les propriétés de l'écou-
lement varient dans le temps.
3.3.2.2 Isovolume
Un écoulement est dit isovolume ou isochore ou solénoïdal lorsque la divergence de la vitesse est nulle. Les
anglos-saxons parlent de "divergence-free ow". Il est dit "non-solénoïdal" dans le cas contraire. Rappelons
qu'un uide incompressible (son volume reste constant) est caractérisé par une divergence nulle de la vitesse
(équation 2.17 en page 26).
A contrario, dans un écoulement irrotationnel, le rotationnel de la vitesse est nul. Si le domaine de l'écou-
lement est dénué de "trous" (simplement connexe), alors il existe une fonction scalaire Φ appelée potentiel
des vitesses, vériant :
~v = −∇Φ
Si, de plus, l'écoulement est à divergence nulle (div ~v = 0), alors le laplacien du potentiel est nul. Il s'agit de
l'équation de Laplace qui facilite la résolution de certains problèmes :
∆Φ = 0
46
3.3. Classication des écoulements
inférieure (resp. supérieure) en tout point à la vitesse du son. On utilise souvent le nombre de Mach, dont la
dénition est donnée dans les tableaux 1.2 et 1.3 en page 10 :
τson τ
son ≡ vLs
Ma ≡ avec
τconv τ
conv ≡ V
L
c'est-à-dire :
V √
Ma ≡ ou bien M a = V. ρ.κS (3.17)
vs
Le tableau 3.1 en page 48 donne la valeur de la vitesse du son et de la compressibilité pour un certain nombre
de matériaux.
Un écoulement est turbulent lorsque le uide s'écoule de manière désorganisée, avec des tourbillons qui
provoquent un mélange du uide et/ou des couches uides qui se décollent les unes des autres.
Le nombre de Reynolds permet de quantier l'importance relative des forces visqueuses (qui tendent à
rendre l'écoulement laminaire) et des forces inertielles (qui tendent à rendre l'écoulement turbulent). On
rappelle sa dénition, donnée dans les tableaux 1.2 et 1.3 en page 10 :
L2
visc ≡ ν
τvisc τ
Re ≡ avec
τconv τ
conv ≡ L V
où τvisc est le temps caractéristique du transport diusif de la quantité de mouvement et τconv celui du
transport convectif de la quantité de mouvement. Compte-tenu de la dénition de la viscosité cinématique
ν ≡ µ/ρ, on peut aussi écrire le nombre de Reynolds en faisant intervenir la viscosité dynamique µ :
ρ.V.L
Re ≡ (3.18)
µ
Le choix de la longueur caractéristique L dépend du problème étudié. Il y a une part d'adaptation et
donc d'arbitraire. Toutefois, bien souvent, il s'agit de la plus petite dimension physique. Ainsi, on retiendra
le diamètre d'une longue canalisation, le diamètre d'une sphère solide dans un grand volume, l'épaisseur d'un
prol pour une aile d'avion, . . .
47
Chapitre 3. Fluides particuliers et lois constitutives, classication des écoulements
Pour un uide compressible, en utilisant les dénitions de la compressibilité adiabatique κS (équation 3.2
en page 39) et de la vitesse du son vs (équation 3.7 en page 41), on peut faire apparaître l'écart de pression
δp observé dans le système en écoulement 4 :
1 δp
L'écoulement est incompressible si et seulement si . 1 (3.20)
vs2 ρ
Ecoulement inertiel. Dans le cas d'un écoulement rapide dominé par l'inertie (nombre de Reynolds
élevé), les écarts de pression δp sont de l'ordre de grandeur de ρ.V 2 . Il s'ensuit que δρ
ρ ∝ M a et donc :
2
Ecoulement visqueux. Dans le cas d'un écoulement lent dominé par la viscosité (nombre de Reynolds
M a2
ρ ∝ Re et donc :
faible), les écarts de pression sont de l'ordre de grandeur de µ.V /L. Il s'ensuit que δρ
Cette condition est plus restrictive que celle pour les écoulements inertiels.
ρ µ λ DAB c̄p α κT vs
kg.m −3
Pa.s W.m−1
.K −1
m .s2 −1
J.K −1
.kg −1
m .s 2 −1
Pa m.s−1
Air (20°C) 1,204 1,825.10−5 0,0251 10−5 1007 2,07.10−5 1
p ≈ 10−5 343
Eau (20°C) 998 1,002.10 −3
0,598 2,3.10 −9
4182 1,43.10 −7
4,6.10 −10
1482
Ethanol (20°C) 789 1,144.10−3 0,169 1,1.10−9 2399 8,93.10−8 1,1.10−9 1144
Glycérol (20°C) 1264 1,519 0,286 2386 9,48.10 −8
2,5.10 −10
1920
Huile moteur 888 0,837 0,145 1881 8,68.10 −8
7.10 −10
1461
Mercure (25°C) 13534 1,534.10 −3
8,51 1,6.10 −9
139 4,51.10 −6
3,7.10 −11
1450
Métaux liquides 2.10 -2.10
3 4
10−4
-10 −3
1-10 2
10 −9
-10
−8
10 3
10 −6
-10 −4
Liquides organiques 103 10−4 -10−3 0,15 10−10 -10−7 103 -3.103 10−8 -10−7
Sels fondus 2.103 10−3 10−1 -1 10−10 103 -4.103 10−7 1850
Huiles silicones 10 3
10 −2
-103
0,1 10−13
-10 −9
2.10 3
10 −7
Table 3.1 Valeurs typiques de la masse volumique ρ, de la viscosité dynamique µ, de la conductivité ther-
mique λ, du coecient de diusion moléculaire DAB , de la capacité calorique massique c̄p , de la diusivité
thermique α = λ/(ρ.c̄p ), de la compressibilité isotherme κT et de la vitesse du son (adiabatique) vs .
4. Cette approche basée sur l'analyse dimensionnelle est inspirée de l'ouvrage d'Etienne Guyon et al., voir page 113.
48
3.4. Données numériques
49
3.5 Solutions des exercices
\nobreak
Chapitre 4
Bilans locaux
Lorsque l'on cherche à résoudre un problème où interviennent des transferts de masse, énergie, quantité
de mouvement,. . . il est naturel de faire un bilan de ces grandeurs extensives. Ces bilans portent parfois le
nom d'équations de conservation même si certaines grandeurs ne se conservent pas (ex :dissipation d'éner-
gie mécanique par frottement visqueux). Dans ce chapitre, les bilans sont établis à l'échelle locale. Dans le
chapitre suivant, ils seront intégrés sur des volumes macroscopiques.
Un tel bilan peut être établi pour un volume de contrôle macroscopique (ni) ou un volume de contrôle
microscopique (innitésimal). Dans le premier cas, le terme de variation est dF/dt tandis que dans le deuxième
cas, nous allons voir qu'il est égal à ∂f /∂t où la densité volumique f est dénie par :
f ≡ lim
F
c'est-à-dire F(t) =
yf (4.2)
(t) .dV
V →0 V
Ω(t)
Convention sur la normale à une surface. En s'appuyant sur la gure 4.1, on considère un volume Ω
de l'espace, séparé de son environnement par la surface ∂Ω. Le vecteur ~n désigne le vecteur unitaire normal
à la surface ∂Ω. Traditionnellement, en physique, mécanique des uides et simulation numérique, la normale
est sortante tandis qu'en génie des procédés et en thermodynamique, elle est entrante. Dans ce document,
lorsque le vecteur ~n apparaît sans indice, il faut considérer que la normale est sortante. Autant que faire
se peut, nous recommandons l'usage d'indice. Nous désignons par ~nin la normale entrante et par ~nout la
normale sortante, l'une et l'autre étant liées par ~nin = −~nout .
51
Chapitre 4. Bilans locaux
Figure 4.1 Représentation schématique d'un volume de contrôle Ω, de sa frontière ∂Ω avec l'environnement
et les normales entrante (~nin ) et sortante (~nout ) à la frontière.
ux net entrant à travers la frontière ∂Ω et le terme Σ qui est la source nette au sein du volume Ω :
dF
=Φ+Σ (4.3)
dt
Par dénition de la densité volumique f , on a :
F ≡
y f.dV (4.4)
Ω
Par dénition de la densité de ux ϕ, le ux net entrant à travers la frontière ∂Ω est :
Φ≡
{ ϕ.~n (4.6)
in .dS
∂Ω
Φ=−
y (div ϕ) .dV (4.7)
Ω
L'équation 4.3 se présente nalement sous la forme d'une égalité faisant intervenir trois intégrales triples.
Cette égalité étant vraie pour tout volume Ω, cela signie que l'on peut écrire :
∂f
= −div ϕ + s (4.8)
∂t
On peut vérier que ce résultat est indépendant de la convention choisie pour l'orientation de la normale à
une surface.
L'expérience montre qu'il existe trois modes de transfert d'une grandeur extensive :
la convection (macroscopique) : il s'agit d'un simple transport de la grandeur extensive qui est "atta-
chée" à la matière en mouvement. Par exemple, les passagers d'un wagon sont convectés (advectés)
par le train. Ils ont la vitesse du train.
la diusion (microscopique) : elle peut s'ajouter à la convection et sera décrite plus en détail dans ce
qui suit.
le rayonnement ou l'induction : ils ne seront pas décrits dans cet ouvrage.
Considérons le cas où le ux total est la somme d'un ux diusif et d'un ux convectif.
52
4.2. Bilan local de masse
Le bilan local peut prendre les deux formes suivantes. La première est issue des équations 4.8 et 4.9. La
deuxième fait, en plus, usage de la dénition de la dérivée matérielle (équation 2.4).
(forme eulérienne) ∂f
∂t = −div (f~v + ϕdi ) + s
Bilan local : (formes lagrangiennes) Df
Dt = −f.div ~v − div ϕdi + s (4.10)
Df
ρ Dtρ = −div ϕdi + s pour f 6= ρ
On précise que la première formulation lagrangienne (Df /Dt = ...) est inconditionnellement correcte tandis
que la deuxième (ρ.D fρ /Dt = ...) n'est vraie que lorsque l'équation de continuité 4.15 est observée, ce qui
est tout de même observé dans de très nombreuses situations.
Lavoisier a postulé en 1789 que la masse totale est une grandeur conservative. Cela est vrai dans la plupart
des systèmes physiques, à l'exception des systèmes relativistes et des réactions nucléaires pour lesquels la
masse peut évoluer en vertu de l'équivalence masse-énergie. Par conséquent, dans des conditions standard,
la somme des termes sources est nulle : X
sk = 0 (4.11)
k
On note ~vk la vitesse du constituant k dans le référentiel xe. Il ne s'agit pas de la vitesse instantanée d'une
molécule particulière de Ak mais de la vitesse moyenne de toutes les molécules de Ak dans un petit volume.
Par construction, la vitesse moyenne (en masse) associée au mélange est égale à :
X
~v ≡ ωk .~vk (4.12)
k
La vitesse ~v ainsi dénie est celle, mesurable par un tube de Pitot ou par anémométrie Laser, qui intervient
dans les bilans de quantité de mouvement.
53
Chapitre 4. Bilans locaux
Le ux massique ϕk du constituant k en kg.m−2 .s−1 s'écrit sous la forme d'une somme de deux termes, l'un
convectif et l'autre diusif :
ϕk
≡ ρk .~vk
ϕk ≡ ϕconv,k + ϕdi,k avec (4.13)
ϕconv,k ≡ ρk .~v
En dénissant le ux total massique ϕtot comme la somme des ux massiques, on arrive à trois relations très
pratiques :
ϕtot = ρ.~v
X
ϕtot ≡ ϕk d'où ϕconv,k = ωk .ϕtot
ϕk = ϕdi,k + ωk .ϕtot
∂ρ.ωk
X j
= −div (ρ.ωk .~v + ϕdi,k ) + νk .M̃k .rj
∂t j
B.L. de masse du constituant k : (4.14)
Dωk X j
div
ρ. = − ϕdi ,k + νk .M̃k .rj
Dt
j
Dans le cas d'un mélange binaire AB , on montre facilement que ϕdi,A = −ϕdi,B . La section 3.2.2 en page
42 présente une courte discussion sur les diérents mécanismes de diusion, à savoir l'équation de Fick liant
ϕdi,k aux concentrations.
div ~v = 0 (4.16)
54
4.3. Bilan local de quantité de matière
Si la grandeur F est la quantité de matière Nk du constituant k , alors la densité volumique f associée est la
concentration ck en mol.m−3 . On dénit la concentration c du mélange par la relation suivante :
X
c≡ ck
On note ~vk la vitesse du constituant k dans le référentiel xe. Il ne s'agit pas de la vitesse instantanée d'une
molécule particulière de Ak mais de la vitesse moyenne de toutes les molécules de Ak dans un petit volume.
Par construction, la vitesse moyenne molaire ~v ∗ du mélange est égale à :
X
~v ∗ ≡ xk .~vk (4.17)
k
La vitesse ~v ∗ ainsi dénie n'est pas facilement mesurable ; elle est diérente de la vitesse moyenne massique
~v qui intervient dans les bilans de quantité de mouvement. On notera en particulier qu'on peut avoir une
vitesse moyenne molaire ~v ∗ non nulle pour un mélange au repos (~v nul).
Le ux total molaire ϕ∗k en mol.m−2 .s−1 s'écrit sous la forme d'une somme de deux termes, l'un convectif et
l'autre diusif :
∗
ϕk
≡ ck .~vk
ϕk ≡ ϕconv,k + ϕdi,k avec
∗ ∗ ∗
ϕ∗conv,k ≡ ck .~v ∗
En dénissant le ux total molaire ϕ∗tot comme la somme des ux molaires, on arrive à trois relations très
pratiques :
ϕ∗tot = c.~v ∗
X
ϕ∗tot ≡ ϕ∗k d'où ϕ∗conv,k = xk .ϕ∗tot (4.18)
∗
ϕk = ϕ∗di,k + xk .ϕ∗tot
délicate et sujette à erreur. Dans la mesure du possible, il est préférable d'utiliser un bilan de masse. Si une
écriture molaire du bilan est indispensable, on pourra se référer à la page 851 de l'ouvrage de Bird, Lightfoot
et Stewart (2nde édition).
= ~v ∗ + ωk . (~vk − ~v ∗ )
P
~v
55
Chapitre 4. Bilans locaux
Les compositions molaires et massiques, ainsi que les ux, sont liés par :
ρk (xk ) = ck .M̃k
= M̃k .ϕ∗k
ϕk
ωk (xk ) = Pxkx.M̃ k
et (4.19)
j .M̃j
ϕconv,k = ωk . M̃k .ϕ∗k
P
xk (ωk ) = Pωkω/M̃
k
/M̃ j j
Exercice 21. la composante selon x de la vitesse d'un écoulement bidimensionnel est de la forme
vx (x, y, t) = Ax. Identier les diérentes composantes vy selon y qui décrivent un uide incompressible.
Exercice 22. on considère un écoulement radial unidimensionnel dans le plan rθ. Le champ de vitesse
est tel que vr = f (r) et vθ = 0. Déterminer les conditions que f (r) doit respecter pour garantir que
l'écoulement est incompressible.
Dans le chapitre précédent, on avait déni la pression p par la relation 2.27 et le tenseur des contraintes de
cisaillement τ par la relation 2.28. On rappelle leurs dénitions ci-dessous :
−Trace (σ)
p≡ et σ ≡ τ − p.I (4.23)
3
1. La deuxième équation s'obtient facilement en utilisant l'équation ωA = f (xA .M̃A , xB .M̃B ), la relation de fermeture
xA + xB = 1 et son corollaire ∇xA = −∇xB .
2. Cela peut être réalisé à titre d'exercice. Il sut d'utiliser la dénition de la dérivée matérielle (équation 2.4), l'équation
de continuité sous sa forme eulérienne et les relations 7.26 et 7.27 du formulaire en page 83.
3. Le produit ~u ⊗ ~v est parfois noté ~u~v sans symbole entre les vecteurs. Nous n'adoptons pas cette dernière notation car elle
porte à confusion.
56
4.4. Bilan local de quantité de mouvement (linéaire)
Avec ces notations, le bilan local de quantité de mouvement s'écrit de l'une ou l'autre manière :
∂t + div (ρ~v ⊗ ~v ) = div τ − ∇p + ρ~g
∂ρ~
v
B.L. de qté de mouvement (éq. de Cauchy) : ρ. ∂~ v
= div τ − ∇p + ρ~g (4.24)
v · ∇) ~v
∂t + (~
ρ D~v = div τ − ∇p + ρ~g
Dt
Cette expression est valable dans un référentiel inertiel (galiléen). Si l'on travaille dans un référentiel non
galiléen, il faut ajouter quelques forces dont la force de Coriolis. Ce cas est présenté en annexe (équation
7.39 en page 88).
L'équation 4.26 peut être simpliée dans le cas d'un écoulement incompressible (voir dénition en page
47). En eet, on a alors de très faibles variations relatives de masse volumique δρ/ρ 1 et l'on peut faire
l'approximation que la masse volumique est constante. D'après le bilan local de masse, équation 4.15 en
page 54, cela implique la nullité de la divergence du champ de vitesse. On en déduit, que pour un uide
incompressible et newtonien de viscosité constante, on a :
div τ = µ.∆~v (4.28)
Le bilan de quantité de mouvement 4.24, lorsque le uide est incompressible et newtonien de viscosité
constante, porte le nom d'équation de Navier-Stokes :
D~v
Equation de Navier-Stokes : ρ. = ρ~g − ∇p + µ.∆~v (4.29)
Dt
Le terme µ.∆~v est le produit de la viscosité et du laplacien de la vitesse qui est une mesure de la courbure
locale du champ de vitesse. Il est nul pour un prol de vitesse linéaire. La gure 4.2 illustre le fait que
les forces visqueuses ont pour eet de linéariser le prol de vitesse en accélérant/entraînant les parcelles
uides en retard et en freinant/ralentissant les parcelles uides en avance par rapport à un prol linéaire.
En d'autres termes, on retiendra que les forces visqueuses agissent comme une force de rappel et stabilisent
l'écoulement en linéarisant les prols de vitesse.
57
Chapitre 4. Bilans locaux
Figure 4.2 Les forces visqueuses tendent à régulariser le champ de vitesse en le rapprochant d'un prol
linéaire.
Cas particulier des écoulements à fort nombre de Reynolds. Ce cas se rencontre pour un écoulement
très turbulent (vitesse très élevée) ou bien un uide parfait (viscosité nulle). Dans l'équation de Navier-Stokes,
le terme visqueux est négligeable et l'on obtient l'équation d'Euler :
D~v
Equation d'Euler : ρ. = ρ~g − ∇p (4.30)
Dt
Notons qu'elle n'est pas valable à proximité d'une interface uide-solide. En eet, comme la vitesse relative du
uide est nulle au contact du solide, il existe nécessairement une zone de faible vitesse où les eets visqueux
ne peuvent être ignorés.
Cas particulier des écoulements stationnaires à faible nombre de Reynolds. Lorsque les eets
visqueux prédominent sur les eets inertiels et que l'écoulement est stationnaire, on parle d'écoulement
rampant (creeping ow). Dans ce cas, l'équation de Navier-Stokes se simplie et prend le nom d'équation
de Stokes :
Equation de Stokes : 0 = ρ~g − ∇p + µ.∆~v (4.31)
58
4.5. Bilan local d'énergie mécanique
Energie cinétique. Si l'on multiplie par ~v le bilan local de quantité de mouvement exprimé par la relation
4.21, on arrive à :
Dēc ec
ρ = (div σ + ρ~g ) .~v avec ēc ≡ (4.32)
Dt ρ
En utilisant l'expression générale d'un bilan local (4.10), on peut donner une version eulérienne du précédent
bilan :
∂ec
= −div (ec .~v ) + (div σ + ρ~g ) .~v (4.33)
∂t
Compte-tenu des égalités 7.23 et 7.24 du formulaire en page 83, on peut écrire le bilan local d'énergie cinétique
sous sa forme canonique :
∂ec
= −div ec .~v − σ T .~v + ρ~g .~v − σ T : (∇~v )
∂t
B.L. d'énergie cinétique : (4.34)
ρ. Dēc = −div −σ T .~v + ρ~g .~v − σ T : (∇~v )
Dt
où, par identication avec la forme canonique 4.10, le terme −σ T .~v représente le ux diusif d'énergie
cinétique. Cette expression est valable dans un référentiel inertiel (galiléen). Si l'on travaille dans un référentiel
non galiléen, il faut ajouter un terme lié à la force d'inertie. Ce cas est présenté en annexe (équation 7.40 en
page 88).
Energie potentielle. L'énergie potentielle Ep = m.Ψ est le produit de la masse m par le potentiel gra-
vitationnel Ψ. On note ep ≡ ρ.Ψ sa densité volumique. Le potentiel gravitationnel étant constant dans le
temps, on a :
∂ep ∂ρ
= Ψ.
∂t ∂t
D'après le bilan local de masse (équation 4.15), on a :
∂ep
= −Ψ.div (ρ~v )
∂t
On peut faire apparaître la divergence d'un produit et ainsi écrire :
∂ep
= −div (ρ~v Ψ) + ρ~v .∇Ψ
∂t
Dans l'objectif d'établir une formulation lagrangienne, on utilise le fait que la force de pesanteur dérive du
potentiel Ψ pour écrire ~g = −∇Ψ et donc :
∂ep
= −div (ep .~v ) − ρ.~g .~v (4.35)
∂t
59
Chapitre 4. Bilans locaux
Dans le membre de droite, on voit apparaître un terme de ux convectif et un terme source. Le bilan local
d'énergie potentiel prend donc la forme suivante :
∂ep
= −div (ep .~v ) − ρ.~g .~v e ≡ ρ.Ψ
∂t p
B.L. d'énergie potentielle : avec ~g ≡ −∇Ψ (4.36)
ρ. Dēp
= −ρ.~g .~v
ē ≡ ep /ρ
Dt p
Energie mécanique. En additionnant membre à membre les équations 4.34 et 4.36, on obtient la forme
générale du bilan local d'énergie mécanique.
∂em
= −div em .~v − σ T .~v − σ T : (∇~v )
∂t
B.L. d'énergie mécanique : (4.37)
ρ. Dēm = −div −σ T .~v − σ T : (∇~v ) = (divσ) .~v
Dt
Le deuxième terme du membre de droite est associé au mouvement d'expansion/compression. Il est positif en
cas d'expansion (champ de vitesse divergent, div ~v > 0) et négatif en cas de compression (champ de vitesse
convergent, div ~v < 0). Il est nul lorsque le uide est incompressible.
Le premier terme est associé au mouvement de cisaillement. Lorsque le tenseur des contraintes est symétrique
(absence de force magnétique), on peut jouer sur les indices et utiliser la dénition du tenseur γ̇ des taux de
cisaillement (équation 2.15 en page 26) pour écrire :
XX ∂vj XX ∂vi XX 1 ∂vi ∂vj
τ T : (∇~v ) = τji . = τij . = τij .γ̇ij avec γ̇ij ≡ +
i j
∂xi i j
∂xj i j
2 ∂xj ∂xi
Pour un uide incompressible et newtonien, l'équation 4.25 en page 57 se simplie en τ = µ.γ̇ . Dans ces
conditions :
1 XX 2
τ T : (∇~v ) = µ. γ̇ij (4.38)
2 i j
Comme la divergence de la vitesse est nulle pour un uide incompressible, on démontre ainsi que le terme
source du bilan d'énergie mécanique, pour un uide incompressible et newtonien, se met sous la forme :
1 XX 2
−σ T : (∇~v ) = − µ. γ̇ij < 0
2 i j
Ce terme "source" est toujours négatif pour un uide newtonien. Cela traduit le fait que le cisaillement d'un
uide visqueux dissipe de l'énergie mécanique en chaleur, de manière irréversible. On retiendra que dans la
pratique, il est négligeable sauf en cas de viscosité très élevée ou de gradients de vitesse très élevés comme
par exemple dans la lubrication, l'extrusion rapide, le malaxage de polymères et les vols d'aéronefs très
rapides.
60
4.6. Bilan local d'énergie interne et d'enthalpie
DP
= ~v .div τ
Dt
Le terme DP/Dt représente la dérivée matérielle de la pression généralisée, c'est-à-dire la variation temporelle
de cette grandeur associée à une parcelle de uide, au cours de son déplacement le long d'une ligne de courant
(tangente en tous points au vecteur vitesse).
Pour un uide incompressible newtonien de viscosité constante (équation 4.25), on a div τ = µ.∆~v . Si
l'on considère deux points 1 et 2 le long d'une ligne de courant, alors la variation de pression généralisée
∆P ≡ P2 − P1 est la suivante : Z 2
∆P = µ.∆~v .d~l (4.40)
1
Dans le cas général, on voit ainsi que la pression généralisée varie le long d'une ligne de courant. L'équation
4.40 montre que la pression généralisée reste constante dans deux situations. La première se rencontre lorsque
le champ de vitesse est à Laplacien nul (que le uide soit visqueux ou parfait). C'est par exemple le cas dans
les écoulements laminaires cisaillés plans. La deuxième situation se rencontre lorsque la viscosité est nulle,
c'est-à-dire lorsque le uide est parfait. L'équation prend alors le nom de théorème de Bernoulli 4 . Pour
un uide parfait, incompressible en écoulement permanent, il stipule que :
1 2
P≡ ρv + p + ρ.Ψ est constante le long d'une ligne de courant (4.41)
2
La gure 4.3 indique graphiquement dans quelles conditions le théorème de Bernoulli est applicable.
Exercice 23. On considère l'écoulement stationnaire d'un uide newtonien incompressible dans un
tube horizontal de rayon R et de longueur L. En fonction du débit volumique QV , calculer la perte
de charge ∆p dans le cas où l'écoulement est laminaire puis en utilisant le théorème de Bernoulli.
Commenter.
61
Chapitre 4. Bilans locaux
Figure 4.3 Zones de validité (en blanc) et d'invalidité (en rose) de la relation de Bernoulli. Source : Majid
Bahrami, Simon Fraser University (lien).
62
4.7. Conditions aux limites
Dū
ρ = −div ~q + τ T : (∇~v ) − p.div ~v
Dt
En utilisant les règles classiques de dérivation puis le bilan local de masse (équation 4.15), on peut établir :
D p 1 Dp p Dρ 1 Dp p
= . − 2. = . + .div ~v
Dt ρ ρ Dt ρ Dt ρ Dt ρ
Dh̄ Dp
ρ = −div ~q + τ T : (∇~v ) + (4.49)
Dt Dt
Il est parfois pratique de connaître l'évolution de la température. Pour ce faire, il sut de connaître l'équation
d'état h̄ = h̄ (T, p) spécique au matériau étudié. En eet, on peut écrire :
Dh̄ ∂ h̄ DT ∂ h̄ Dp
= . + .
Dt ∂T p Dt ∂p T Dt
∂∂ph̄ = v̄. [1 − α.T ] 1
1 + ∂∂ ln
ln ρ
= ρ. T
T p
où c̄p représente la capacité calorique massique et v̄ = 1/ρ le volume massique et α ≡ −1
ρ .
∂ρ
∂T le
p
coecient de dilatation. Par conséquent, le bilan d'enthalpie peut prendre la forme suivante. Notons bien
que la température n'étant pas une grandeur extensive, il ne s'agit pas d'un bilan local de température. Cette
équation est connue sous le nom d'équation de la chaleur.
DT ∂ ln ρ Dp
Equation de la chaleur : ρc̄p . = −div ~q + τ T : (∇~v ) − . (4.50)
Dt ∂ ln T p Dt
2
De la thermodynamique, on rappelle la relation de Mayer c̄p − c̄v = T.αρ.κT , la dénition du coecient adiaba-
c̄p
tique γ ≡ c̄v et son lien avec la vitesse du son adiabatique γ = vs .ρ.κT .
2
63
Chapitre 4. Bilans locaux
Dans certains cas de gure, la valeur d'une grandeur f telle que la vitesse, la température ou la concentration
d'une espèce est imposée. On parle alors de condition de Dirichlet :
Dans d'autres cas, c'est une valeur du ux de f qui est imposé (ux thermique, ux de matière ou ux de
force, c'est-à-dire une contrainte). On parle alors de condition de Neumann :
Pour un point ~r sur la frontière ∂Ω : (∇f )(~r) .~nout = une valeur imposée
On s'aidera de la gure 4.4 pour repérer les deux milieux, le vecteur unitaire ~n normal à l'interface et le
vecteur unitaire ~t tangent à l'interface. Le tableau 4.1 résume les diérentes situations qui seront explicitées
ci-dessous.
Table 4.1 Conditions aux limites de Dirichlet (valeur imposée) ou de Neumann (ux imposé).
4.7.1 Mécaniques
Conditions de Dirichlet. Si les milieux 1 et 2 ne traversent pas la frontière (absence d'interpénétration),
alors on peut écrire :
~v1 .~n = ~v2 .~n
Si les milieux 1 et 2 adhèrent l'un à l'autre (ce qui sera le cas pour les uides visqueux mais pas pour les
uides parfaits), alors on peut écrire :
~v1 .~t = ~v2 .~t
64
4.7. Conditions aux limites
Conditions de Neumann. Rappelons que la densité de ux de quantité de mouvement est une contrainte.
L'égalité des contraintes normales (la pression) s'exprime simplement par p1 = p2 si l'interface entre les mi-
lieux 1 et 2 est plane. Dans le cas général d'une interface courbe, il faut faire intervenir la tension supercielle
par la loi de Laplace :
1 1
p2 = p1 + γ12 . + 0
R R
où γ12 est la tension de surface entre les milieux 1 et 2 et où R et R0 sont les rayons de courbure principaux de
l'ellipsoïde qui tangente l'interface. La pression est toujours plus élevée à l'intérieur de la concavité. L'égalité
des contraintes tangentielles se traduit par :
T~1 .~t = T~2 .~t
où T~i est le vecteur-contrainte déni par l'équation 2.21 en page 30. Pour mémoire, ce vecteur contrainte est
lié au tenseur des contraintes par la relation 2.24 qui s'énonce T~ = σ · ~nout . Par analyse dimensionnelle, on
peut montrer que le vecteur contrainte est proportionnelle à la viscosité dynamique. Par conséquent, pour
un couple (X/uide parfait) ou bien (milieu dense/gaz), alors l'un des vecteurs contraintes est nul ou très
faible devant l'autre et l'on écrire T~ = ~0.
4.7.2 Thermiques
Conditions de Dirichlet. Sauf s'il y a une réaction chimique exothermique ou endothermique à l'interface
(exemple : oxydation d'un métal par un gaz oxydant), il y a continuité de la température :
T1 = T2
Conditions de Neumann. Sauf s'il y a une réaction chimique exothermique ou endothermique à l'inter-
face (exemple : oxydation d'un métal par un gaz oxydant), il y a continuité du ux de chaleur :
~q1 = ~q2
Si un des deux milieux est isolant thermique (conductivité thermique inniment faible), on parle de conditions
adiabatiques et l'on a ~q = ~0.
4.7.3 Chimiques
Conditions de Dirichlet. Sauf s'il y a une réaction chimique à l'interface, on écrit l'égalité des potentiels
chimiques. Dans certains congurations, cette égalité implique l'égalité des concentrations (voir cours de
thermodynamiques). On retiendra :
µ1 = µ2
65
Chapitre 4. Bilans locaux
Conditions de Neumann. On écrit l'égalité des ux de matières (en mol/m2 /s) sauf s'il y a une réaction
chimique à l'interface qui consomme ou produit le constituant étudié.
66
4.8. Synthèse
4.8 Synthèse
Les tableaux 4.2 et 4.3 donnent l'expression du bilan local pour plusieurs grandeurs extensives F . Ces
expressions sont très générales. Elles s'appliquent à tous les matériaux (solide, liquide, gaz) et ne font pas
d'approximations autres que l'absence de force électro-magnétique.
∂f D fρ
= −div (f~v + ϕdi ) + s ou bien ρ = −div ϕdi + s pour f 6= ρ
∂t Dt
QDM ∂ρ~
v
∂t = −div (ρ~v ⊗ ~v − σ) + ρ~g Dt = div σ + ρ~
ρ D~
v
g
E. cin. ∂ec
= −div ec .~v − σ T .~v + ρ~g .~v − σ T : (∇~v ) Dt = div σ .~
ρ Dē
T
∂t
c
v + ρ~g .~v − σ T : (∇~v )
∂ep Dēp
E. pot. ∂t = −div (ep .~v ) − ρ.~g .~v ρ Dt = −ρ~g .~v
E. int. ∂u
∂t = −div (u.~v + ~q) + σ T : (∇~v ) ρDDt = −div ~
ū
q + σ T : (∇~v )
E. tot. ∂e
= −div e.~v + ~q − σ T .~v Dt = −div ~
ρ Dē
∂t q − σ T .~v
Enth. ∂h
∂t = −div (h.~v + ~q) + τ T : (∇~v ) + Dp
Dt Dt = −div ~
ρD h̄
q + τ T : (∇~v ) + Dp
Dt
Dt = −div ~
∂ ln ρ
ρc̄p . DT q + τ T : (∇~v ) − ∂ ln T . Dp
Dt
p
Table 4.2 Bilans locaux, dans leur formulation eulérienne et lagrangienne. Les densités volumiques sont
ec ≡ 12 ρv 2 , ep ≡ ρΨ. Les grandeurs barrées f¯ désignent les densités massiques : f¯ ≡ f /ρ.
67
4.9 Solutions des exercices
Exercice 21 pour un uide incompressible, on sait par l'équation 4.16 que div ~v = 0. La dénition de la
∂vy
divergence (page 93) conduit à écrire ∂v
∂x + ∂y = 0. Compte-tenu de la forme imposée de vx , on déduit
x
∂v
que ∂yy = −A. Une intégration conduit à vy (x, y, t) = −A.y + f (x, t) où f est une fonction quelconque (par
exemple nulle).
Exercice 23 Pour un écoulement laminaire, on peut utiliser l'équation de Navier-Stokes 4.29 et la géométrie
du système (~v = v(r) .~ez et p = p(z) ) pour trouver que la perte de charge est p0 − pL = π8 . µ.Q
R4 .L.
V
Pour un écoulement turbulent, on ne peut pas utiliser l'équation de Navier-Stokes. En revanche, on peut
utiliser le théorème de Bernoulli 4.41 qui s'écrit p0 − pL = 12 ρ. vL − v02 . Puisque le prol de vitesse est
2
invariant par translation selon z , les vitesses v0 et vL sont identiques, ce qui conduit à p0 − pL = 0.
Les deux résultats sont diérents. Cela met en lumière le fait que le théorème de Bernoulli ne s'applique pas
aux écoulements où la viscosité joue un rôle important.
Chapitre 5
Bilans macroscopiques
Dans de nombreux contextes liés à la mécanique des uides, on cherche à calculer des grandeurs telles que
des débits, des températures, des forces, des puissances mécaniques ou thermiques. Dans ces situations, il
n'est pas particulièrement utile de connaître le champ de température ou de vitesse en tout point du volume
de contrôle. Au contraire, une grandeur macroscopique est souvent susante, au moins dans une première
étape de dimensionnement.
Le schéma 5.1 montre un système macroscopique délimité par des pointillés rouges. Le uide qu'il contient
reçoit une puissance mécanique externe Ẇ0 (par exemple par un ensemble arbre+hélice en rotation ou bien
par un barreau magnétique mis en mouvement à distance) ainsi qu'une puissance thermique externe Q̇0 (par
exemple par contact avec un échangeur de chaleur situé au coeur du système ou bien par excitation des
molécules par des micro-ondes). Le système est alimenté par plusieurs conduites. Les parois solides exercent
sur le uide interne une force F~s/f et chauent le uide (densité de ux de chaleur par conduction ~q).
Par dénition, l'approximation unidimensionnelle est valable lorsque l'écoulement est considéré
comme unidimensionnel, c'est-à-dire lorsque les grandeurs p, T, ~v , ρ sont considérées comme uniformes sur
la section S des conduites connectées au système étudié. L'indice i désigne le numéro de la conduite d'arri-
vée/sortie.
69
Chapitre 5. Bilans macroscopiques
dF y (s − div (ϕ { (5.1)
= di + f~v )) .dV + f~v∂Ω .~nout .dS + Fext
dt
Ω(t) ∂Ω(t)
dF
=
y s.dV + { ϕdi .~nin .dS +
{ f (~v − ~v∂Ω ) .~nin .dS + Fext (5.3)
dt
Ω(t) ∂Ω(t) ∂Ω(t)
70
5.3. Bilan macroscopique de masse
Malgré les apparences de simplicité, cette équation n'est pas évidente à mettre en oeuvre. En eet, dans de
nombreux cas pratiques, on ne connaît pas l'évolution des frontières du volume matériel (exemple typique :
écoulement turbulent).
Quoique plus longue, cette expression est très utilisée en pratique car les intégrales se calculent sur des
domaines connus et invariants dans le temps.
où ṁi,in désigne le débit massique du uide (de masse volumique ρi ) entrant par la section Si de la conduite
à la vitesse moyenne ~vi . Par construction, il s'agit d'une grandeur algébrique : le débit est positif si la vitesse
est dirigée vers l'intérieur du système et négatif sinon.
Lorsque la frontière imaginaire xe ∂Ω du système peut être considérée comme l'union de la frontière ∂Ωs
entre le uide interne et les parois solides et de la frontière ∂Ωf entre le uide interne et le uide externe,
on obtient :
dP~
= m~g +
{ σ.~nout .dS +
{ σ.~nout .dS +
{ (ρ~v ⊗ ~v ) .~nin .dS +
{ ~ ext
(ρ~v ⊗ ~v ) .~nin .dS +F (5.9)
dt
∂Ωs ∂Ωf ∂Ωs ∂Ωf
| {z } | } | {z } |
T1 T3
{z {z }
T2 T4
71
Chapitre 5. Bilans macroscopiques
Le terme T 1 représente la force F~s/f que les parois solides exercent sur le uide contenu dans le système
(voir équation 2.25 en page 31). C'est souvent la grandeur que l'on cherche à calculer. Compte-tenu du fait
que σ = τ − pI et que la contrainte de cisaillement τ est souvent négligeable en norme par rapport à la
pression p d'un uide 1 , le terme T 2 se simplie en une intégrale sur la surface de p.dS
~ . Lorsque les parois sont
imperméables, la vitesse du uide à leur contact est nulle. Par conséquent, le terme T 3 est identiquement
nul. On en déduit l'expression suivante du bilan macroscopique de quantité de mouvement :
Cas particulier n°1 : l'état stationnaire. Dans ce cas, la résultante F~tot/f des forces qui s'appliquent
sur le uide (gravité + parois solides + pression du uide + force extérieure) est égale et opposée au débit
de quantité de mouvement entrant dans le système. Il s'ensuit que l'équation 5.10 se simplie et prend la
forme suivante, appelée théorème d'Euler :
Th. d'Euler :
{ (ρ~v ⊗ ~v) .~n avec
{ p.~n ~ ext
F~tot/f = − in .dS F~tot/f ≡ m~g + F~s/f + in .dS + F
∂Ωf ∂Ωf
(5.11)
On notera que le sens physique du terme ṁi .~vi se perçoit facilement dans le cas d'une fusée qui décolle (sa
quantité de mouvement augmente) lorsqu'un gaz de combustion en sort à grande vitesse.
dE
= ΣE +
{ T
~q − τ − pI .~v .~nin .dS +
1 2 {
ρv + ρgz + ρū ~v .~nin .dS (5.13)
dt 2
∂Ω ∂Ω
où la source (macroscopique) d'énergie ΣE est égale à la somme de la puissance mécanique Ẇ0 (autre que
les forces de pression et de pesanteur) et de la puissance thermique Q̇0 . Comme la contrainte de cisaillement
τ est souvent négligeable, en norme, à la pression p d'un uide, on peut réorganiser les termes et donner
l'expression générale du bilan macroscopique d'énergie :
dE
= Ẇ0 + Q̇tot +
{
1 2
ρv + ρgz + ρh̄ ~v .~nin .dS (5.14)
dt 2
∂Ω
où l'intégrale sur la surface inclut le travail des forces de pression et de pesanteur. Dans l'équation 5.14,
l'enthalpie massique h̄ est liée à l'énergie interne massique ū par la relation h̄ = ū + p/ρ 2 et la puissance
1. Pour un uide newtonien 1000 fois plus visqueux que l'eau (µ = 1 Pa.s) et une vitesse de cisaillement typique d'un
procédé "gentil" (γ̇ = 102 s−1 ), la contrainte de cisaillement τ est de l'ordre de 102 Pa, ce qui est négligeable devant la pression
atmosphérique (105 Pa).
2. Cette relation s'obtient en divisant par m la relation connue H = U + pV .
72
5.6. Bilan macroscopique d'énergie mécanique
thermique totale Q̇tot contient le terme Q̇0 lié à la chaleur apportée directement au coeur du système et un
autre terme lié au ux conductif ~q intégré sur la surface de contrôle :
Q̇tot ≡ Q̇0 +
{ ~q.~n (5.15)
in .dS
∂Ω
Bilan d'énergie avec les variables T et p Lorsque les grandeurs T et p sont imposées, on peut considérer
l'enthalpie comme une fonction de T et p et écrire :
∂ h̄
∂T = c̄p
∂ h̄ ∂ h̄
dh̄ = .dT + .dp avec p (5.17)
∂T p ∂p T
∂p∂ h̄ 1 ∂ ln ρ
= ρ . 1 + ∂ ln T
T p
Bilan d'énergie avec les variables s et p Bien souvent, la température n'est pas connue à l'avance et sa
valeur dépend des irréversibilités dans l'écoulement (le cisaillement visqueux transforme l'énergie mécanique
en énergie thermique). Dans ce cas, il est plus judicieux de considérer l'enthalpie comme une fonction de la
pression p et de l'entropie s pour écrire :
∂ h̄
∂ h̄ ∂∂sh̄ = T
dh̄ = .ds + .dp avec p (5.18)
∂s p ∂p s ∂ h̄
= ρ1
∂p
s
Graphiquement, la température T étant strictement positive, les isobares h̄ (s, p0 ) sont des fonctions crois-
santes de s. Par ailleurs, la masse volumique étant aussi strictement positive, l'isobare (s, pHigh ) se situera
au dessus de l'isotherme (s, pLow ) dans le diagramme h̄ − s.
dEm
=
{
1 2
ρv .~v + ρ.Ψ.~v − σ T .~v .~nin .dS −
y
σ T : (∇~v ) .dV + Σext
Ec
dt 2
∂Ω Ω
où Σext
Ec désigne la source (macroscopique) d'énergie cinétique extérieure (par exemple amenée par une agi-
tation mécanique ou bien un champ électro-magnétique). Compte-tenu du fait que le tenseur des contraintes
peut être décomposé selon σ = τ − p.I et de l'égalité σ T : (∇~v ) = τ T : (∇~v ) − p.div ~v , on a :
dEm { P.~v.~n {τ T
y p.div ~v.dV − y τ T
: (∇~v ) .dV + Σext
= in .dS − .~v .~nin .dS + Ec
dt
∂Ω ∂Ω Ω Ω
où P désigne la pression généralisée dont on rappelle la dénition donnée dans l'équation 4.39 en page 61 :
1
P ≡ p + ρv 2 + ρ.Ψ (5.19)
2
73
Chapitre 5. Bilans macroscopiques
A l'aide du théorème de Green-Ostrogradsky (équation 7.48 en page 90) et de l'égalité 7.22 du formulaire en
page 83, ce bilan d'énergie mécanique peut aussi s'écrire :
dEm { P.~v.~n y p.div ~v.dV + y [~v.div τ ] .dV + Σ ext
(5.20)
= in .dS + Ec
dt
∂Ω Ω Ω
La première intégrale volumique fait apparaître la divergence de la vitesse (qui est nulle lorsque le uide est
incompressible, voir équation 4.16 en page 54). La deuxième intégrale volumique est toujours négative. Elle
représente la vitesse à laquelle l'énergie mécanique est dégradée en énergie thermique par les frottements
visqueux. Dans la section 5.6 en page 74, on verra que cette puissance dissipée est proportionnelle à la viscosité
dynamique lorsque le uide est simple (newtonien, incompressible, . . .). Pour simplier les équations, on pose :
Pvisc. ≡
y [~v.div τ ] .dV (5.21)
Ω
Cas particulier. Le bilan d'énergie mécanique prend une forme particulièrement simple dans le cas où les
trois conditions suivantes sont observées :
régime stationnaire (∂/∂t = 0 et d/dt = 0)
uide incompressible (div ~v = 0)
uide newtonien de viscosité constante (équation 4.25)
Dans le cas d'un écoulement stationnaire, pour un uide incompressible et en l'absence de source externe
d'énergie cinétique, le bilan macroscopique d'énergie mécanique (équation 5.22) se simplie en :
{ P.~v.~n
in .dS = −Pvisc.
∂Ω
Pour un uide incompressible et newtonien de viscosité constante, on peut utiliser la relation 4.28 en page
57 pour écrire la puissance dissipée (en [W]) :
Pvisc. = µ.
y ~v.∆~v.dV
Ω
Si, en plus de toutes les hypothèses précédentes, le uide est parfait (µ = 0) et s'il y a une entrée et une
sortie, alors le bilan d'énergie mécanique prend la forme suivante, que l'on peut qualier de théorème de
Bernoulli macroscopique :
1 2
Pentrée = Psortie avec P≡ ρv + p + ρ.Ψ (5.24)
2
74
Chapitre 6
Compléments
Figure 6.1 Longueur d'entrée hydrodynamique (cas laminaire et turbulent). Source : Larousse.
Notons ∆P la perte de charge dans une conduite cylindrique de longueur L et de diamètre d parcourue par
un uide de masse volumique ρ ayant une vitesse moyenne V . L'analyse dimensionnelle montre qu'elle peut
75
Chapitre 6. Compléments
où fF est le coecient de perte de charge de Fanning. Les deux coecients fD et fF sont sans
dimension et sont liés par la relation :
fD = 4.fF (6.3)
La gure 6.2 synthétise une grande masse de données expérimentales. Il s'agit du diagramme de Moody.
Figure 6.2 Diagramme de Moody permettant de calculer une perte de charge en conduite cylindrique.
fD = 64
Re pour un écoulement laminaire (théorie)
ε/d
= −2.Log10 . 2,51
pour un écoulement turbulent (Colebrook, 1939)
1
√
fD Re.
√ + 3,7
fD 1,11
√1 = −1, 8.Log10 . 6,9
+ ε/d pour un écoulement turbulent (Haaland, 1983)
fD Re 3,7
Considérons les transferts de chaleur entre un uide circulant dans le tube à température constante. La
longueur d'entrée thermique (nécessaire à l'établissement d'un prol de température invariant par translation)
est : (
Ltherm Re.P r
20 pour un écoulement laminaire
=
d 10 pour un écoulement turbulent
76
6.1. Corrélations empiriques pour écoulements typiques
Pour un écoulement laminaire, le coecient de transfert de chaleur moyen peut être calculé par :
1/3 0,14
Re.P r µbulk
N u = 1, 86. . pour z < Ltherm (Sieder et Tate, 1936)
L/D µparoi
0, 065.Re.P r. D
2/3 pour tout z (Mills)
L
N u = 3, 66 +
1 + 0, 04. Re.P r. D
L
Notons que la valeur N u∞ = 3, 66 est celle mesurée pour un écoulement laminaire complètement développé.
Pour un écoulement turbulent dans un tube lisse et long (L > 60.D), avec un uide tel que 0, 7 < P r < 160,
on utilise la formule :
(
0, 3 si le uide est refroidi
N u = 0, 023.Re .P r pour Re > 10000 (Dittus et Boelter) avec n =
0,8 n
0, 4 si le uide est réchaué
Toutes les propriétés physiques doivent être évaluées à la température bulk moyenne du uide (moyenne
arithmétique entre la température bulk à l'entrée et à la sortie).
77
Chapitre 6. Compléments
78
Chapitre 7
Annexes mathématiques
La convention d'Einstein permet d'alléger les écritures. Lorsqu'il y a des sommes de grandeurs avec des
indices, on ajoute les termes contenant des indices répétés (muets). Les autres indices sont dits explicites ou
réels. Quelques exemples :
X X X X
ai bi = ai bi = aj bj = aj bj et aij bj = aij bj et ~u = ui .~ei = ui .~ei
i j j i
79
Chapitre 7. Annexes mathématiques
Un tenseur d'ordre n est caractérisé par n indices. Si l'on désigne par (~e1 , ~e2 , ~e3 ) la base orthonormée directe
dans laquelle on travaille, alors un tenseur ~u d'ordre 1, un tenseur σ d'ordre 2 et un tenseur T d'ordre 3
s'écrivent :
~u = ui .~ei
σ = σij .~ei ⊗ ~ej
T = Tijk .~ei ⊗ ~ej ⊗ ~ek
où les grandeurs ~ei ⊗ ~ej et ~ei ⊗ ~ej ⊗~ek sont des produits dyadiques. Pour l'instant, on peut interpréter ~ei ⊗~ej
comme un élément constituant une base de l'espace vectoriel des tenseurs d'ordre 2 (idem pour les tenseurs
d'ordre 3).
80
7.2. Tenseurs et produits tensoriels
Produit dyadique
P
~u ⊗ ~v σ σ= ui .vj .~ei ⊗ ~ej σij = ui .vj
i,j
Produit contracté
P
~u · ~v s s= ui .vi s = ui .vi
i
P
σ · ~u ~v ~v =σij .uj .~ei vi = σij .uj
i,j
P
σ·τ Π Π= σik .τkj .~ei ⊗ ~ej Πij = σik .τkj
i,j,k
Table 7.2 Expression générale des diérents produits tensoriels. La colonne "Représentation" donne les
composantes des vecteurs ou matrices dans la base (~e1 , ~e2 , ~e3 ) (et uniquement dans cette base).Le symbole
εijk désigne le symbole de permutation de Levi-Civita (équation 7.2). Les indices (i, j, k) prennent chacun
leur valeur parmi (x, y, z), (r, θ, z) ou (r, θ, φ).
Propriétés
T
(~u ⊗ ~v ) = ~v ⊗ ~u (7.3)
σT : τ = σ : τ T (7.4)
81
Chapitre 7. Annexes mathématiques
Le gradient d'un tenseur T d'ordre n est un tenseur d'ordre n + 1 noté ∇T (ou grad T ) déni par :
dT ≡ ∇T · d~x (7.5)
où ε est le tenseur d'ordre 3 de Levi-Civita (voir équation 7.2) et où ∇~u désigne le gradient de vitesse (un
tenseur d'ordre 2). Le symbole " :" désigne le produit doublement contracté déni dans le tableau 7.2.
La divergence d'un tenseur T d'ordre n est un tenseur d'ordre n − 1 noté div T et déni par :
Le laplacien d'un tenseur T d'ordre n est un tenseur d'ordre n noté ∆T et déni par la divergence du
gradient :
∆T ≡ div (∇T ) (7.8)
7.3.2 Formulaire
Les symboles s et t désignent des scalaires, ~u et ~v des vecteurs et σ et τ des tenseurs d'ordre 2. Pour plus
de facilité, les équations suivantes peuvent être établies en coordonnées cartésiennes mais elles sont valables
quel que soit le système de coordonnées.
Gradient
∇ (s.t) = s.∇t + t.∇s (7.9)
∇ (s.~u) = s.∇~u + ~u ⊗ ∇s (7.10)
∇ (~u · ~v ) = (~u · ∇) ~v + (~v · ∇) ~u + ~u ∧ (rot ~v ) + ~v ∧ (rot ~u) (7.11)
Rotationnel
rot (∇s) = ~0 (7.12)
rot (s.~u) = s.rot ~u + (∇s) ∧ ~u (7.13)
rot (~u ∧ ~v ) = (div ~v ) .~u − (div ~u) .~v . + (~v · ∇) ~u − (~u · ∇) ~v (7.14)
rot (rot ~u) = ∇ (div ~u) − ∆~u (7.15)
rot (∆~u) = ∆ (rot ~u) (7.16)
Divergence
div (rot ~u) = 0 (7.17)
div (∇~u) = ∆~u (7.18)
div (s.~u) = s.div ~u + (∇s) · ~u (7.19)
div (~u ∧ ~v ) = (rot ~u) .~v − (rot ~v ) .~u (7.20)
82
7.4. Conventions sur la construction du tenseur des contraintes
σ ≡ out
ver σ
T~ = in
ver σ · ~nin (7.30)
Puisque ~nin = −~nout , on calcule bien la même grandeur physique T~ . Notons pour nir qu'avec la convention
"normale entrante", inver σxx > 0 implique que la matière est "en compression".
83
Chapitre 7. Annexes mathématiques
Avec le tenseur des contraintes ainsi déni, le vecteur-contrainte n'est pas calculé par l'équation 2.24 mais
par :
T~ = ~nout · out
hor σ (7.32)
Cette dernière expression n'est pas très correcte syntaxiquement. En Europe, on utilise l'écriture ( out
T
hor σ) ·~nout
pour désigner la même grandeur (ce qui est bien égal à σ · ~nout ).
Pour information, on trouve ci-dessous les conventions adoptées dans certains ouvrages de référence :
ce document : normale sortante et assemblage vertical
Bird : normale entrante et assemblage horizontal
Guyon : normale sortante et assemblage vertical + tenseur des contraintes de cisaillement noté σ 0 au
lieu de τ dans ce document (Guyon, p. 163)
Fox : normale sortante et assemblage horizontal
7.4.3 Synthèse
Quelle que soit la convention adoptée, on a :
les grandeurs dF~ et T~ sont identiques (car elles ont un sens physique immédiat)
un matériau en compression est caractérisée par une pression p positive
Par ailleurs, on rappelle les formules de passage :
~nin = −~nout
hor σ (ver
T
out = out σ)
(7.33)
hor σ = − (ver
out σ)
T
in
ver
ver
in σ = −out σ
84
7.4. Conventions sur la construction du tenseur des contraintes
Table 7.3 Diérentes conventions. Dans ce document, on adopte la convention en haut à droite : assemblage
vertical et normale sortante.
85
Chapitre 7. Annexes mathématiques
7.5 Démonstrations
7.5.1 Bilan des forces
Le bilan local de quantité de mouvement se résume à un bilan des forces. Il sut pour s'en convaincre de
reconnaître que la dimension d'une force est celle d'un débit de quantité de mouvement par unité de surface :
kg. m
1 N = 1 s.ms2 . Considérons donc le volume élémentaire décrit sur la gure 2.10 en page 31, on peut écrire :
D~v
dm. = dF~contact + dF~volumique
Dt
où dm est la masse et D~ Dt l'accélération du volume élémentaire. En page 31, on a déni le tenseur des
v
contraintes σ par un assemblage vertical, ce qui a permis d'écrire l'équation 2.25 établissant le lien entre la
force F~contact et le tenseur des contraintes σ . En projection sur l'axe Ox, on obtient :
dm ≡ ρ.dx.dy.dz
dFcontact,x ≡ σxx|x+dx − σxx|x .dy.dz + σxy|y+dy − σxy|y .dx.dz + σxz|z+dz − σxz|z .dx.dy
dFvolumique,x ≡ dm.gx
où l'on rappelle que σij est la composante selon l'axe i de la force exercée sur la face de normale j . Le bilan
des forces s'écrit alors :
Dvx ∂σxx ∂σxy ∂σxz
ρ. Dt = ∂x + ∂y + ∂z + ρ.gx
Dv ∂σyx ∂σyy ∂σyz
ρ. Dty = ∂x + ∂y + ∂z + ρ.gy
ρ. Dvz
∂σzx ∂σzy ∂σzz
Dt = ∂x + ∂y + ∂z + ρ.gz
On obtient nalement :
D~v
= div σ + ρ.~g
ρ. (7.34)
Dt
où div σ désigne la divergence du tenseur des contraintes.
86
7.5. Démonstrations
Par construction, le couple de volume est le produit d'une constante kz nie et du volume élémentaire :
Kz = kz .dx.dy.dz
Après division par dx.dy.dz , la vitesse angulaire ω̇z du cube autour de l'axe Oz est :
h i
2 2
ρ. (dx) + (dy) Dθ̇
z
. = (σxy − σyx ) + kz
12 Dt
Pour que la vitesse angulaire reste nie lorsque dx et dy tendent vers zéro, il faut que le membre de droite
soit identiquement nul. On a donc :
σyx = σxy + kz (7.37)
La force gravitationnelle est une force volumique qui ne crée par de couple. On a donc ~k = ~0. En revanche,
dans les uides magnétiques (gaz stellaires, noyau liquide de la terre, autres systèmes technologiques,. . .), le
terme kz peut être non nul. Nous retiendrons qu'en dehors des systèmes magnéto-hydrodynamiques, cette
constante est nulle et donc que le tenseur des contraintes est symétrique :
Le tenseur des contraintes σ est symétrique (sauf systèmes magnétiques) : σij = σji (7.38)
Comme la densité de courant électrique J~ (en A.m−2 ) est dénie par J~ ≡ ρe~v , on peut écrire :
f~ = ρe .E
~ + J~ ∧ B
~
L'évolution temporelle du champ magnétique est donnée par l'équation d'induction, dont la forme complète
est complexe. Lorsque le nombre de Reynolds magnétique est très grand devant 1 (ce qui est très souvent le
cas en astrophysique et dans les machines de fusion), alors on est dans le cadre de la MHD idéale et la loi
d'induction prend la forme simpliée suivante :
~
∂B
≈ rot ~v ∧ B
~
∂t
D'après les équations de Maxwell, on a rotB ~ = µ0 J~ + ε0 µ0 ∂ E~ . Dans le cas particulier où la vitesse caracté-
∂t
ristique de l'écoulement est très faible devant la vitesse de la lumière, les eets relativistes sont négligeables
et l'approximation suivante est valable :
rot B
~ ≈ µ0 J~
Lorsque l'on étudie l'écoulement d'un objet en rotation, il est souvent pratique de se placer dans le réfé-
rentiel lié à l'objet. Ce référentiel n'est pas inertiel et le bilan des forces doit être ré-écrit. Cette situation se
87
Chapitre 7. Annexes mathématiques
rencontre par exemple lorsque l'on étudie la circulation des océans ou bien le noyau uide de la terre ou bien
encore l'écoulement dans une éprouvette placée dans une centrifugeuse ou bien l'atomisation d'une nappe de
liquide résultant de l'impact d'un jet sur une surface plane tournante.
On considère tout d'abord un repère xe R d'origine O et un repère mobile R0 d'origine O0 . Dans le
repère xe, le centre O0 du repère R0 est à la position ~r(O0 ) , se déplace à la vitesse ~v(O0 ) et possède une
accélération ~a(O0 ) . On appelle Ω ~ le vecteur rotation instantanée de R0 par rapport à R, avec Ω~ = θ̇.~k et θ̇ la
vitesse angulaire de rotation. Précisons que la rotation de R0 par rapport à R est le changement de direction
des seuls axes (O0 x0 , O0 y 0 , O0 z 0 ) par rapport aux axes (Ox, Oy, Oz), sans rapport avec la trajectoire du point
O0 . On considère ensuite un point M particulier de l'espace dont la position, la vitesse et l'accélération sont
notées ~r, ~v et ~a dans le repère xe R et ~r 0 , ~v 0 et ~a 0 dans le repère mobile R0 . On peut établir trois lois de
composition de mouvement :
Position : ~r = ~r 0 + ~rO0
Vitesse : ~v ~ ∧ ~r 0
= ~v 0 + ~vO0 + Ω
~
Accélération : ~ ∧ Ω
~a = ~a 0 + ~aO0 + Ω ~ ∧ ~r 0 + dΩ ~ ∧ ~v 0
∧ ~r 0 + 2Ω
dt
m~a = F~vraie
où F~vraie est la résultante des vraies forces. En vertu de la loi de composition des accélérations, on a :
m~a 0 = F~vraie + F~Centrifuge + F~Euler + F~Coriolis − m.~aO0
Application 1 : dans le référentiel inertiel, le bilan local de quantité de mouvement (équation 4.24 en page
57) s'énonce :
D~v
ρ = div σ + ρ~g
Dt
Si l'on souhaite travailler dans le référentiel mobile, il prendra donc la forme suivante :
D~v 0
ρ. = div σ + ρ~g + F~Centrifuge + F~Euler + F~Coriolis − m.~aO0 (7.39)
Dt
Application 2 : dans le référentiel inertiel, le bilan local d'énergie cinétique (équation 4.34 en page 59) peut
s'écrire :
Dec
= (∇ · σ + ρ~g ) .~v
Dt
En multipliant l'équation 7.39 par ~v 0 , on trouve :
Dē0c
ρ. = ~v 0 . div σ + ρ~g + F~Centrifuge + F~Euler − m.~aO0 (7.40)
Dt
Notons ici que le terme ~v 0 .F~Coriolis n'apparaît pas parce qu'il est identiquement nul par construction : la
force de Coriolis ne travaille pas.
88
7.5. Démonstrations
Par exemple, si x est une vitesse, on aura [x] = L.T −1 et (a1 , . . . , aM ) = (0, 1, −1, 0, 0, 0, 0).
La dimension du produit de deux grandeurs x et y ainsi que d'une puissance de x s'établit facilement :
k
(7.43)
k
x = [x]
Par ailleurs, on peut établir les équations aux dimensions suivantes :
dx ∂x [x]
= = (7.44)
dy ∂y [y]
Z
y(x) .dx = [y] . [x] (7.45)
[x]
[∇x] = (7.46)
L
[x]
(7.47)
2
∇ x = [∆x] = 2
L
Par exemple, la dimension du produit d'un diamètre par une vitesse est celle d'une longueur au carré divisée
par un temps soit [d.V ] = L2 .T −1 . De même, la dimension d'une vitesse élevée au cube est celle d'une longueur
au cube divisée par un temps au cube soit [V 3 ] = L3 .T −3 et l'ensemble des exposants est (0, 3, −3, 0, 0, 0, 0).
7.5.5.2 Enoncé
Considérons un ensemble
QM deami
n grandeurs physiques x1 , x2 , . . . , xn . Pour chacune d'elle, l'équation aux di-
mensions est [xi ] = m=1 Dm . Notons r le rang de la matrice A dénie par :
a11 ... a1n
. .. ..
A≡ . .
. .
aM 1 ... aM n
7.5.5.3 Démonstration
Considérons un ensemble de n coecients k1 , . . . , kn et formons la grandeur N en multipliant des puis-
sances de xi :
Yn
N≡ xki i
i=1
89
Chapitre 7. Annexes mathématiques
Cherchons les conditions à réunir pour que N soit sans dimension. Pour cela, écrivons son équation aux
dimensions :
" n # n
Y Y k
ki
[N ] = xi = [xi ] i
i=1 i=1
ou encore :
!ki n
n
Y M
Y n Y
Y M M Y
Y n M
Y
P
ami .ki
ami ami .ki ami .ki i=1
[N ] = Dm = Dm = Dm = Dm
i=1 m=1 i=1 m=1 m=1 i=1 m=1
Les jeux de coecients ki rendant la grandeur N sans dimension sont tels que l'exposant de chacun des Dm
est nul, c'est-à-dire lorsque les M équations suivantes sont observées :
Pn
a1i .ki = 0
i=1
... = 0
Pn
a .kMi i = 0
i=1
Il s'agit d'un système sous-déterminé (le nombre n d'inconnues est supérieur au nombre M d'équations) et
homogène (second membre nul). En algèbre linéaire, on a vu que, si le déterminant de la matrice A est non
nul, alors il existe n − r solutions linéairement indépendantes, où r est le rang de la matrice.
2. On trouve aussi les dénominations théorème de Green ou théorème de Gauss car sa découverte fut faite indépen-
damment par ces chercheurs
90
7.6. Théorèmes importants
91
Chapitre 7. Annexes mathématiques
σ : τ = (σxx .τxx + σxy .τyx + σxz .τzx ) + (σyx .τxy + σyy .τyy + σyz .τzy ) + (σzx .τxz + σzy .τyz + σzz .τzz )
Le gradient d'un vecteur ~u est le tenseur d'ordre 2 (voir note de bas de page 3 ) :
∂ux ∂ux ∂ux
∂x ∂y ∂z
X ∂ui
.~ei ⊗ ~ej d'où ∂uy ∂uy ∂uy
(7.55)
∇~u ≡ ∇~u = ∂x ∂y ∂z
∂xj
i,j∈{x,y,z} ∂uz ∂uz ∂uz
∂x ∂y ∂z
3. Le gradient d'un vecteur peut être noté ∇~u ou grad~u. On peut montrer que ∇~u = (∇ ⊗ ~u)T .
92
7.7. Bilans en coordonnées cartésiennes (~ex , ~ey , ~ez )
93
Chapitre 7. Annexes mathématiques
94
7.7. Bilans en coordonnées cartésiennes (~ex , ~ey , ~ez )
Note : si le vecteur normal était orienté vers l'intérieur, alors il y aurait un signe "-" devant le deuxième
terme du membre de droite. Par ailleurs, les indices de τ seraient permutés si le tenseur des contraintes était
construit par un assemblage horizontal.
Notation développée :
γ̇xx = 2. ∂v
∂x
x
∂vy
γ̇yy = 2 ∂y
γ̇zz = 2 ∂v
∂z
z
∂vy
γ̇xy = γ̇yx = ∂vx
∂y + ∂x
∂vy
γ̇yz = γ̇zy = ∂z + ∂vz
∂y
= ∂x + κ − 3 .div ~
µ. 2 ∂v 2µ
τxx x
v
h i
∂v ∂vy
= µ. 2 ∂yy + κ − 2µ 3 .div ~ avec div ~v ≡ ∂vx ∂vz
τyy v ∂x + ∂y + ∂z
= µ. 2 ∂z + κ − 2µ 3 .div ~
∂vz
τzz v
h i
∂vy
τxy = τyx = µ. ∂v∂y + ∂x
x
h i
∂v
τyz = τzy = µ. ∂zy + ∂v ∂y
z
= =
∂vz ∂vx
τzx τxz µ. ∂x + ∂z
95
Chapitre 7. Annexes mathématiques
Notation développée :
2
∂ 2 vx ∂ 2 vx
∂vx ∂vx ∂vx ∂vx ∂p ∂ vx
ρ + vx + vy + vz =− +µ + + + ρgx
∂t ∂x ∂y ∂z ∂x ∂x2 ∂y 2 ∂z 2
2
∂ 2 vy ∂ 2 vy
∂vy ∂vy ∂vy ∂vy ∂p ∂ vy
ρ + vx + vy + vz =− +µ + + + ρgy
∂t ∂x ∂y ∂z ∂y ∂x2 ∂y 2 ∂z 2
2
∂ 2 vz ∂ 2 vz
∂vz ∂vz ∂vz ∂vz ∂p ∂ vz
ρ + vx + vy + vz =− +µ + + + ρgz
∂t ∂x ∂y ∂z ∂z ∂x2 ∂y 2 ∂z 2
Notation développée :
∂T ∂T ∂T ∂T ∂qx ∂qy ∂qz T ∂ ln ρ Dp
ρc̄p . + vx + vy + vz =− + + + τ : (∇~v ) − .
∂t ∂x ∂y ∂z ∂x ∂y ∂z ∂ ln T p Dt
avec :
∂vx ∂v
+ τyx . ∂xy + τzx . ∂v
τ T : (∇~v ) ≡ τxx . ∂x ∂x
z
∂vy
+τxy . ∂v
∂y
x
+ τyy . ∂y + τzy . ∂vz
∂y
∂vx
∂vy ∂vz
+τxz . ∂z + τyz . ∂z + τzz . ∂z
et :
Dp ∂p ∂p ∂p ∂p
≡ + vx + vy + vz
Dt ∂t ∂x ∂y ∂z
~q = −λ.∇T
Notation développée :
∂T
∂x
~q = −λ. ∂T
∂y
∂T
∂z
96
7.7. Bilans en coordonnées cartésiennes (~ex , ~ey , ~ez )
Notation développée :
et
∂vx ∂vy ∂vz
div ~v = + +
∂x ∂y ∂z
97
Chapitre 7. Annexes mathématiques
On note (~ex , ~ey , ~ez ) la base orthonormée directe en coordonnées cartésiennes et (~er , ~eθ , ~ez ) la base orthonormée
directe en coordonnées cylindriques. On note que la première est indépendante de la position du point M
considéré tandis que la deuxième l'est. On a les relations suivantes :
~er = cos θ.~ex + sin θ.~ey + 0.~ez
~ex = cos θ.~er + (− sin θ).~eθ + 0.~ez
~eθ = − sin θ.~ex + cos θ.~ey + 0.~ez et ~ey = sin θ.~er + cos θ.~eθ + 0.~ez (7.68)
~e = 0.~ex + 0.~ey + 1.~ez ~e = 0.~er + 0.~eθ + 1.~ez
z z
On peut déduire les 9 égalités suivantes, qui seront utiles pour trouver l'expression des opérateurs diérentiels
en coordonnées cylindriques :
98
7.8. Bilans en coordonnées cylindriques (~er , ~eθ , ~ez )
σrr .ur + σrθ .uθ + σrz .uz
σ · ~u =
σθr .ur + σθθ .uθ + σθz .uz
σzr .ur + σzθ .uθ + σzz .uz
σrr .τrr + σrθ .τθr + σrz .τzr σrr .τrθ + σrθ .τθθ + σrz .τzθ σrr .τrz + σrθ .τθz + σrz .τzz
σ·τ =
σθr .τrr + σθθ .τθr + σθz .τzr σθr .τrθ + σθθ .τθθ + σθz .τzθ σθr .τrz + σθθ .τθz + σθz .τzz
σzr .τrr + σzθ .τθr + σzz .τzr σzr .τrθ + σzθ .τθθ + σzz .τzθ σzr .τrz + σzθ .τθz + σzz .τzz
σ : τ = (σrr .τrr + σrθ .τθr + σrz .τzr ) + (σθr .τrθ + σθθ .τθθ + σθz .τzθ ) + (σzr .τrz + σzθ .τθz + σzz .τzz )
Compte-tenu des expressions de r, θ, z en fonction de x, y, z (voir équations 7.67), cela permet d'obtenir :
∂•
= cos θ. ∂• sin θ ∂•
∂r θz − r . ∂θ rz
∂x yz
∂• ∂• cos θ ∂• (7.70)
∂y = sin θ. ∂r θz
+ r . ∂θ rz
xz
∂•
∂•
∂z xy = ∂z rθ
99
Chapitre 7. Annexes mathématiques
7.68). On utilise ensuite la règle de dérivation des fonctions composées (équations 7.70) et l'on simplie pour
obtenir :
∂• ~eθ ∂• ∂•
∇• = ~er . + . + ~ez . (7.71)
∂r θz r ∂θ rz ∂z rθ
Le gradient d'un scalaire s est un vecteur. Pour obtenir son expression, la première méthode, lente, consiste
à partir de son expression en coordonnées cartésiennes (équation 7.54 de gauche), à exprimer les vecteurs
~ex , ~ey et ~ez en fonction de ~er , ~eθ et ~ez (équations 7.68) puis à utiliser l'équation 7.70. La deuxième méthode,
rapide, consiste à appliquer l'opérateur ∇ à la grandeur scalaire s :
∂s
∂r
∂s 1 ∂s ∂s
.~ez soit
∇s = .~er + . .~eθ + ∇s = 1 ∂s
r . ∂θ
∂r θz r ∂θ rz ∂z rθ
∂s
∂z
Le gradient d'un vecteur ~u est un tenseur d'ordre 2. Pour obtenir son expression, la première méthode,
lente, consiste à partir de son expression en coordonnées cartésiennes (équation 7.55 de gauche), à remplacer
les vecteurs ~ex , ~ey et ~ez en fonction de ~er , ~eθ et ~ez (équations 7.68), à exprimer le vecteur ~u sous la forme
ur .~er + uθ .~eθ + uz .~ez puis à utiliser les relations 7.70. La deuxième méthode, plus rapide, consiste à partir de
T
l'égalité générale ∇~u = (∇ ⊗ ~u) . On utilise l'expression 7.71 pour l'opérateur ∇, la dénition du produit
dyadique ⊗ donnée dans le tableau 7.2 et les relations 7.69. On obtient :
∂ur 1 ∂ur uθ ∂ur
∂r r ∂θ − r ∂z
∇~u = ∂uθ 1 ∂uθ
+ ur ∂uθ
∂r r ∂θ r ∂z
∂uz 1 ∂uz ∂uz
∂r r ∂θ ∂z
La divergence d'un vecteur ~u est un scalaire. Elle s'obtient comme la trace du gradient de ce vecteur :
1 ∂(r.ur ) 1 ∂uθ ∂uz
div ~u = tr (∇~u) = . + . + (7.72)
r ∂r r ∂θ ∂z
100
7.8. Bilans en coordonnées cylindriques (~er , ~eθ , ~ez )
∂uθ
τ : (∇~u) = τrr . ∂u ∂uz
∂r + τrθ . ∂r + τrz . ∂r
r
∂uθ
+ τθr . 1r . ∂u 1 1 ∂uz
∂θ − uθ + τθθ . r . ∂θ + ur + τθz . r . ∂θ
r
1 ∂uθ
+ τzr . ∂u ∂uz
∂z + τzθ . r . ∂z + τzz . ∂z
r
∂ρ Dρ
+ div (ρ~v ) = 0 ou bien = −ρ.div ~v
∂t Dt
Notation développée :
∂ρ 1 ∂ (rρvr ) 1 ∂ (ρvθ ) ∂ (ρvz ) Dρ 1 ∂ (rvr ) 1 ∂vθ ∂vz
+ + + = 0 ou bien = −ρ. + +
∂t r ∂r r ∂θ ∂z Dt r ∂r r ∂θ ∂z
101
Chapitre 7. Annexes mathématiques
Notation développée :
∂ωA
∂r
~ di,A
ϕ = −ρ.DAB . 1 ∂ωA
r . ∂θ
∂ωA
∂z
1 ∂ 2 ωA ∂ 2 ωA
∂ωA ∂ωA vθ ∂ωA ∂ωA 1 ∂ ∂ωA
ρ. + vr + + vz = ρDAB . r + 2 + + rA
∂t ∂r r ∂θ ∂z r ∂r ∂r r ∂θ2 ∂z 2
On peut préférer utiliser les fractions molaires xA au lieu des fractions massiques ωA . Pour cela, il sut
N
d'utiliser l'équation précédente en remplaçant (ρ, ωα , ~v , rα ) par respectivement (c, xα , ~v ∗ , Ralpha − xα Rβ )
P
β=1
N
où ~vα est la vitesse de α dans le référentiel xe, ~v ∗ = ~vα la vitesse molaire moyenne, cα la concentration
P
α=1
N
molaire de α, c = cα et Rα le taux volumique de production molaire de α.
P
α=1
vθ2
h i
1 ∂(rτrr )
ρ ∂vr
∂t + vr ∂v
∂r +
r vθ ∂vr
r ∂θ − r + vz ∂v
∂z
r
= − ∂p
∂r + r ∂r + 1 ∂τrθ
r ∂θ − τθθ
r + ∂τrz
∂z + ρgr
2
1 ∂ (r τθr )
∂vθ
+ vr ∂v vθ ∂vθ vr vθ
+ vz ∂v = − 1r ∂p 1 ∂τθθ ∂τθz τrθ −τθr
ρ ∂t ∂r +
θ
r ∂θ + r ∂z
θ
∂θ + r2 ∂r + r ∂θ + ∂z + r + ρgθ
h i
1 ∂(rτzr )
∂vz
+ vr ∂v vθ ∂vz
+ vz ∂v = − ∂p 1 ∂τzθ ∂τzz
ρ ∂t ∂r +
z
r ∂θ ∂z
z
∂z + r ∂r + r ∂θ + ∂z + ρgz
Note : si vecteur normal était orienté vers l'intérieur, alors il y aurait un signe "-" devant le deuxième
terme du membre de droite. Par ailleurs, les indices de τ seraient permutés si le tenseur des contraintes était
construit par un assemblage horizontal.
102
7.8. Bilans en coordonnées cylindriques (~er , ~eθ , ~ez )
Notation développée :
γ̇rr = 2. ∂v
∂r
r
= 1 ∂vθ vr
γ̇θθ 2. r ∂θ + r
γ̇zz = 2. ∂v
∂z
z
Notation développée :
= ∂r + κ − 3 µ .div ~
µ. 2 ∂v 2
r
τrr v
1 ∂(rvr )
= µ. 2 r ∂θθ + vrr + κ − 32 µ .div ~v avec div ~v =
1 ∂v 1 ∂vθ ∂vz
τθθ r ∂r + r ∂θ + ∂z
= ∂z + κ − 3 µ .div ~
µ. 2 ∂v 2
z
τzz v
= = µ. 1r ∂v vθ ∂vθ
τrθ τθr ∂θ − r + ∂r
r
= = µ. ∂v 1 ∂vz
τθz τzθ ∂z + r ∂θ
θ
= = µ. ∂v ∂vr
z
τzr τrz ∂r + ∂z
v2 1 ∂ 2 vr ∂ 2 vr
∂vr ∂vr vθ ∂vr ∂vr ∂p ∂ 1 ∂ 2 ∂vθ
ρ + vr + − θ + vz =− +µ (rvr ) + 2 − + + ρgr
∂t ∂r r ∂θ r ∂z ∂r ∂r r ∂r r ∂θ2 r2 ∂θ ∂z 2
103
Chapitre 7. Annexes mathématiques
1 ∂ 2 vθ ∂ 2 vθ
∂vθ ∂vθ vθ ∂vθ vr vθ ∂vθ 1 ∂p ∂ 1 ∂ 2 ∂vr
ρ + vr + + + vz =− +µ (rvθ ) + 2 + + +ρgθ
∂t ∂r r ∂θ r ∂z r ∂θ ∂r r ∂r r ∂θ2 r2 ∂θ ∂z 2
1 ∂ 2 vz ∂ 2 vz
∂vz ∂vz vθ ∂vz ∂vz ∂p 1 ∂ ∂vz
ρ + vr + + vz =− +µ r + 2 + + ρgz
∂t ∂r r ∂θ ∂z ∂z r ∂r ∂r r ∂θ2 ∂z 2
Notation développée :
∂T ∂T vθ ∂T ∂T 1 ∂(r.qr ) 1 ∂qθ ∂qz ∂ ln ρ Dp
ρc̄p . + vr + + vz =− + + + τ T : (∇~v ) − .
∂t ∂r r ∂θ ∂z r ∂r r ∂θ ∂z ∂ ln T p Dt
avec :
∂uθ
τ T : (∇~u) = τrr . ∂u ∂uz
∂r + τθr . ∂r + τzr . ∂r
r
∂uθ
+ τrθ . 1r . ∂u 1 1 ∂uz
∂θ − uθ + τθθ . r . ∂θ + ur + τzθ . r . ∂θ
r
1 ∂uθ
+ τrz . ∂u ∂uz
∂z + τθz . r . ∂z + τzz . ∂z
r
et :
Dp ∂p ∂p vθ ∂p ∂p
≡ + vr + + vz
Dt ∂t ∂r r ∂θ ∂z
~q = −λ.∇T
Notation développée :
∂T
∂r
~q = −λ. 1 ∂T
r . ∂θ
∂T
∂z
Notation développée :
1 ∂2T ∂2T
∂T ∂T vθ ∂T ∂T 1 ∂ ∂T
ρc̄p . + vr + + vz = λ. r + 2 2 + + µ.Φ
∂t ∂r r ∂θ ∂z r ∂r ∂r r ∂θ ∂z 2
104
7.8. Bilans en coordonnées cylindriques (~er , ~eθ , ~ez )
1 ∂vr 2
2 ∂vr ∂vz 2 2
vθ
+ r ∂θ + r ∂θ + ∂v − 23 . [div ~v ]
∂ 1 ∂vz
+ r ∂r r ∂z
θ
+ ∂z + ∂r
et
1 ∂(r.vr ) 1 ∂vθ ∂vz
div ~v = . + . +
r ∂r r ∂θ ∂z
105
Chapitre 7. Annexes mathématiques
En utilisant la règle de dérivation des fonctions composées, on peut dériver la grandeur • avec :
∂• cos θ. cos φ sin φ
sin θ. cos φ. ∂• ∂• ∂
= + . − .
∂x ∂r θφ r ∂θ rφ r. sin θ ∂φ
yz
rθ
∂• ∂• cos θ. sin φ ∂• cos φ ∂•
∂y = sin θ. sin φ. ∂r θφ + r . ∂θ rφ + r. sin θ . ∂φ
xz rθ
∂• ∂•
sin θ ∂•
∂z xy = cos θ. ∂r θφ − r . ∂θ rφ
On note (~ex , ~ey , ~ez ) la base orthonormée directe en coordonnées cartésiennes et (~er , ~eθ , ~eφ ) la base orthonor-
mée directe en coordonnées sphériques. On note que la première est indépendante de la position du point M
considéré tandis que la deuxième l'est. On a les relations suivantes :
~er = sin θ. cos φ.~
e x + sin θ. sin φ.~e y + cos θ.~
e z
~ex = sin θ. cos φ.~er + cos θ. cos φ.~eθ − sin φ.~eφ
~eθ = cos θ. cos φ.~ex + cos θ. sin φ.~ey − sin θ.~ez et ~ey = sin θ. sin φ.~er + cos θ. sin φ.~eθ + cos φ.~eφ
~e = − sin φ.~e + cosφ.~e + 0.~e ~e = cos θ.~e − sin θ.~e + 0.~e
φ x y z z r θ φ
106
7.9. Bilans en coordonnées sphériques (~er , ~eθ , ~eφ )
107
Chapitre 7. Annexes mathématiques
∂uθ ∂u
τ : (∇~u) = τrr ∂u
∂r + τrθ ∂r + τrφ ∂r
r φ
uθ ∂u
+ τθr 1r ∂u + τθθ 1r ∂u ur
+ τθφ 1r ∂θφ
∂θ − r ∂θ + r
r θ
1 ∂ur uφ 1 ∂uθ uφ cos θ 1 ∂uφ ur uθ cos θ
+ τφr r. sin .
θ ∂φ − r + τφθ .
r. sin θ ∂φ − .
r sin θ + τφφ r. sin θ . ∂φ + r + r . sin θ
∂ρ Dρ
+ div (ρ~v ) = 0 ou bien = −ρ.div ~v
∂t Dt
Notation développée :
∂ρ 1 ∂ r2 ρvr 1 ∂ (ρvθ sin θ) 1 ∂ρvφ
+ 2 + + =0
∂t r ∂r r sin θ ∂θ r sin θ ∂φ
ou bien
" #
Dρ 1 ∂ r 2 vr 1 ∂ (vθ sin θ) 1 ∂vφ
= −ρ. 2 + +
Dt r ∂r r sin θ ∂θ r sin θ ∂φ
108
7.9. Bilans en coordonnées sphériques (~er , ~eθ , ~eφ )
Notation développée :
∂ωA
∂r
~ di,A
ϕ = −ρ.DAB . 1 ∂ωA
r . ∂θ
1 ∂ωA
r. sin θ . ∂φ
∂ 2 ωA
∂ωA ∂ωA vθ ∂ωA vφ ∂ωA 1 ∂ ∂ωA 1 ∂ ∂ωA 1
ρ. + vr + + = ρDAB . 2 r2 + 2 sin θ + 2 2 +rA
∂t ∂r r ∂θ r sin θ ∂φ r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
On peut préférer utiliser les fractions molaires xA au lieu des fractions massiques ωA . Pour cela, il sut
N
d'utiliser l'équation précédente en remplaçant (ρ, ωα , ~v , rα ) par respectivement (c, xα , ~v ∗ , Ralpha − xα Rβ )
P
β=1
N
où ~vα est la vitesse de α dans le référentiel xe, ~v ∗ = ~vα la vitesse molaire moyenne, cα la concentration
P
α=1
N
molaire de α, c = cα et Rα le taux volumique de production molaire de α.
P
α=1
!
∂vθ ∂vθ vθ ∂vθ vφ ∂vθ vr vθ − vφ2 cot θ
ρ + vr + + +
∂t ∂r r ∂θ r sin θ ∂φ r
1 ∂p 1 ∂ 1 ∂ 1 ∂τθφ (τrθ − τθr ) − τφφ cot θ
r3 τθr +
=− + 3 (τθθ sin θ) + + + ρgθ
r ∂θ r ∂r r sin θ ∂θ r sin θ ∂φ r
109
Chapitre 7. Annexes mathématiques
∂vφ ∂vφ vθ ∂vφ vφ ∂vφ vr vφ + vφ vθ cot θ
ρ + vr + + +
∂t ∂r r ∂θ r sin θ ∂φ r
1 ∂p 1 ∂ 3
1 ∂ 1 ∂τφφ (τrφ − τφr ) + τθφ cot θ
=− + 3 r τφr + (τφθ sin θ) + + + ρgφ
r sin θ ∂φ r ∂r r sin θ ∂θ r sin θ ∂φ r
Note : si vecteur normal était orienté vers l'intérieur, alors il y aurait un signe "-" devant le deuxième
terme du membre de droite. Par ailleurs, les indices de τ seraient permutés si le tenseur des contraintes était
construit par un assemblage horizontal.
γ̇rr = 2. ∂v
∂r
r
= 1 ∂vθ vr
γ̇θθ 2. r ∂θ + r
∂vφ
γ̇φφ = 2. 1
r. sin θ . ∂φ + vr +vθ cot θ
r
= = ∂ vθ 1 ∂vr
γ̇rθ γ̇θr r ∂r r + r ∂θ
vφ
γ̇θφ = γ̇φθ = sin θ ∂
r ∂θ sin θ + 1 ∂vθ
r sin θ ∂φ
vφ
γ̇φr = γ̇rφ = ∂
r ∂r r + 1 ∂vr
r sin θ ∂φ
= ∂r + κ − 3 µ .div ~
µ. 2 ∂v 2
r
τrr v
= µ. 2 r ∂θ + vrr + κ − 23 µ .div ~v
1 ∂vθ
τθθ
h i
∂vφ
= 1 vr +vθ cot θ
+ κ − 32 µ .div ~v
τφφ µ. 2 r sin θ ∂φ + r
= =
∂ vθ 1 ∂vr
τrθ τθr µ. r ∂r r + r ∂θ
h i
vφ
τθφ = τφθ = µ. sinr θ ∂θ
∂ 1
sin θ + r sin θ ∂φ
∂vθ
h i
vφ
τφr = τrφ = 1
µ. r sin ∂vr
θ ∂φ + r ∂
∂r r
avec :
1 ∂(r2 .vr ) 1 ∂(vθ sin θ) 1 ∂vφ
div ~v = . + . + .
r2 ∂r r. sin θ ∂θ r. sin θ ∂φ
110
7.9. Bilans en coordonnées sphériques (~er , ~eθ , ~eφ )
1 ∂2 ∂ 2 vr
∂p 1 ∂ ∂vr 1
+ µ 2 2 r 2 vr + 2
=− sin θ + 2 2 + ρgr
∂r r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
!
∂vθ ∂vθ vθ ∂vθ vφ ∂vθ vr vθ − vφ2 cot θ
ρ + vr + + +
∂t ∂r r ∂θ r sin θ ∂φ r
∂ 2 vθ
1 ∂p 1 ∂ ∂vθ 1 ∂ 1 ∂ 1 2 ∂vr 2 cos θ ∂vφ
=− +µ 2 r2 + 2 (vθ sin θ) + 2 2 + − + ρgθ
r ∂θ r ∂r ∂r r ∂θ sin θ ∂θ r sin θ ∂φ2 r2 ∂θ r2 sin2 θ ∂φ
∂vφ ∂vφ vθ ∂vφ vφ ∂vφ vr vφ + vφ vθ cot θ
ρ + vr + + +
∂t ∂r r ∂θ r sin θ ∂φ r
∂ 2 vφ
1 ∂p 1 ∂ 2 ∂vφ 1 ∂ 1 ∂ 1 2 ∂vr 2 cos θ ∂vθ
=− +µ 2 r + 2 (vφ sin θ) + 2 2 + 2 + 2 2 +ρgφ
r sin θ ∂φ r ∂r ∂r r ∂θ sin θ ∂θ r sin θ ∂φ2 r sin θ ∂φ r sin θ ∂φ
Notation développée :
1 ∂(r2 .qr )
∂T ∂T vθ ∂T vφ ∂T 1 ∂ 1 ∂qφ T ∂ ln ρ D
ρc̄p . + vr + + =− 2
+ (q θ . sin θ) + +τ : (∇~
v )− .
∂t ∂r r ∂θ r sin θ ∂φ r ∂r r sin θ ∂θ r sin θ ∂φ ∂ ln T p
avec :
∂vθ ∂v
τ T : (∇~v ) = τrr ∂v
∂r + τθr ∂r + τφr ∂r
r φ
vθ 1 ∂vθ 1 ∂vφ
+ τrθ 1r ∂v vr
∂θ − r + τθθ r ∂θ + r + τφθ r ∂θ
r
1 ∂vr vφ 1 ∂vθ vφ cos θ 1 ∂vφ vr vθ cos θ
+ τrφ r. sin .
θ ∂φ − r + τ θφ .
r. sin θ ∂φ − .
r sin θ + τφφ r. sin θ . ∂φ + r + r . sin θ
et :
Dp ∂p ∂p vθ ∂p vφ ∂p
≡ + vr + +
Dt ∂t ∂r r ∂θ r sin θ ∂φ
111
Chapitre 7. Annexes mathématiques
~q = −λ.∇T
Notation développée :
∂T
∂r
~q = −λ. 1 ∂T
r . ∂θ
1 ∂T
r. sin θ . ∂φ
Notation développée :
∂2T
∂T ∂T vθ ∂T vφ ∂T 1 ∂ ∂T 1 ∂ ∂T 1
ρc̄p . + vr + + = λ. 2 r2 + 2 sin θ + 2 2 +µ.Φ
∂t ∂r r ∂θ r sin θ ∂φ r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
avec :
2
∂vr 2 1 ∂vθ vr 2 1 ∂vφ vr +vθ cot θ
Φ = 2 ∂r + r ∂θ + r + r sin θ ∂φ + r
h i2 h i2
vθ 1 ∂vr 2 vφ ∂vθ ∂vr vφ
∂ sin θ ∂ 1 1 ∂
+ r ∂r r + r ∂θ + r ∂θ sin θ + r sin θ ∂φ + r sin θ ∂φ + r ∂r r
− 23 div ~v
112
7.10. Ressources documentaires
Introduction to Fluid Mechanics. Philip J. Pritchard et John W. Mitchell, Wiley, 8th Edition, 12 mars
2015, 672 pages. https://books.google.fr/books?id=vXu5BwAAQBAJ
Hydrodynamique physique. Étienne Guyon, Jean-Pierre Hulin et Luc Petit, EDP Sciences, 3ème édi-
tion, septembre 2012, 724 pages. https://laboutique.edpsciences.fr/produit/595/9782759808939/
Hydrodynamique%20physique
Transport Phenomena. R. Byron Bird, Warren E. Stewart, Edwin N. Lightfoot, John Wiley & Sons,
Revised 2nd edition, 2007, 905 pages. https://books.google.fr/books?id=L5FnNlIaGfcC
Lien Google vers le livre de Bird, Lightfoot et Stewart (2nd edition)
Visualisation des écoulements. National Committee for uid mechanics lms, Flow visualization, S.J.
Kline, Stanford University
www.emse.fr/~bonnefoy/Public/MFTVideos/FlowVisualizationShapiro.pdf
7.10.3 Vidéos
Série de lms de Shapiro. Ascher Shapiro, National Committee for Fluid Mechanics Films, 1961, 547
pages. http://web.mit.edu/hml/ncfmf.html
7.10.5 Tutoriels
Beginner's guide to aeronautics. Site web de la NASA https://www.grc.nasa.gov/www/k-12/airplane/
bga.html
113
Chapitre 7. Annexes mathématiques
114
Chapitre 8
Méthodologie et formulaire
115
Chapitre 8. Méthodologie et formulaire
Pour répondre à une question de Mécanique des Fluides et Transferts, on peut s'inspirer de la méthodologie
suivante, en l'adaptant aux spécicités du problème étudié :
Méthodologie
1. dénir le système étudié et son état initial
• dénir sans ambiguïté un volume de contrôle Ω par sa frontière ∂Ω avec l'environnement
(gure 4.1 en page 52). En particulier, bien clarier si cette frontière est immobile, advectée
par le uide ou bouge de manière dé-corrélée du uide (gure 5.2 en page 70).
• dénir les conditions aux limites (page 63)
• dénir l'état initial du système
• identier la nature de l'information recherchée : locale (chapitre 4) ou macroscopique
(chapitre 5)
3. résoudre mathématiquement
• identier les invariances (pour chaque grandeur successivement) :
invariance temporelle : stationnaire/instationnaire
invariance par translation selon une direction de l'espace
invariance par rotation autour d'un axe ou d'un point
invariance par retournement autour d'un plan
• choisir le système de coordonnées pertinents : cartésien, cylindrique ou sphérique
• choisir un volume élémentaire qui simplie la résolution mathématique tout en capturant la
physique intéressante
• tenir compte des conditions aux limites
116
MASSE QUANTITE DE ENERGIE ENERGIE
MOUVEMENT MECANIQUE INTERNE
Eq. de continuité Eq. de Cauchy Bilan d'énergie mécanique Bilan d'énergie interne
(4.15 page 54) (4.24 page 57) (4.37 page 60) (4.46 page 62)
Dρ v m ū
= −ρ.div ~v ρ D~ = ρ~g + div σ ρ Dē = [div τ − ∇p] .~v ρD = −div ~q + σ T : (∇~v )
Dt Dt Dt Dt
ρ cst.
ρ cst µ cst. ep =cst h̄ ≡ ū+ ρp
Fluide newt.
Fluide parf.
Re = 0 Fluide no phase
parfait ρ cst. change
stat.
stat.
Eq. de la chaleur
Eq. de Stokes Eq. d'Euler Th. de Bernoulli
(4.50 page 63)
(4.31 page 58) (4.30 page 58) (4.41 page 61)
v 1 ρc̄p . DT
Dt
= −div ~q + τ T :
ρ~g − ∇p + µ.∆~v = ~0 ρ. D~ = ρ~g − ∇p ρv 2 + p + ρgz = cst. ln ρ
Dp
Dt 2 (∇~v ) − ∂∂ ln .
T p Dt
117
Chapitre 8. Méthodologie et formulaire
Formulaire
Cinématique
h i h i
T ∂vj
d ≡ 1
2. ∇~v + (∇~v ) soit dij ≡ 21 . ∂vi
∂xj + ∂xi
∇~v = d + r avec p. 23
h i h i
T ∂vj
1
soit rij ≡ 21 . ∂vi
r ≡ 2. ∇~v − (∇~v ) ∂xj − ∂xi
~ ≡ 1 .rot ~v = ω
Ω ~ et ~ ∧ d~x
r.d~x = Ω p. 24
2
Contraintes
Le terme σij .dS est la composante selon la direction ~ei de la force dF~
σ ≡ τ − p.I p. 32
Divers
α≡ −1 ∂ρ
ρ .
∂T
dρ
ρ = −α.dT + κT .dp avec
p p. 39
κ ≡ 1 . ∂ρ
T ρ ∂p
r T
q
vs ≡ ∂p
∂ρ = 1
ρ.κS et κT = γ.κS p. 41
S
118
Index
Approximation thermo-élastique, 39
de Boussinesq, 40 Colburn-Schlichting (corrélation de), 77
unidimensionnelle, 69 Colebrook (corrélation de), 76
Assemblage horizontal/vertical, 83 Compressibilité
adiabatique κS , 39, 47
Bernoulli (théorème de), 61, 74 isotherme κT , 39, 47
Bilan Conditions aux limites, 63
forme générale, 51 chimiques, 65
Bilan local de Dirichlet, 64
d'énergie cinétique, 59 de Neumann, 64
d'énergie interne, 62 de Robin, 64
d'énergie mécanique, 60, 61 mécaniques, 64
d'énergie potentielle, 60 thermiques, 65
d'enthalpie, 63 Conductivité thermique λ, 42
de masse, 53, 54 Constituant, 53, 54
de quantité de matière, 54 Continuité (équation de), 54
de quantité de mouvement, 56, 57, 86 Contraintes
forme générale, 53 eective de von Mises, 33
Bilan macroscopique de cisaillement τ , 32
d'énergie, 72 pression hydrostatique p, 32
d'énergie mécanique, 74 tenseur σ des, 30
de masse, 71 vecteur-contrainte T~ , 30
de quantité de mouvement, 71 Convention
forme générale, 70 de sommation d'Einstein, 79
Bond (nombre de), 11 sur la normale entrante/sortante, 31, 51, 83
Boussinesq (approximation de), 40 sur le tenseur des contraintes, 83
Corrélation
Capacité calorique cp et cv , 40 de Churchill et Rose, 77
Capillaire (nombre), 11 de Colburn-Schlichting, 77
Cauchy de Colebrook, 76
équation de, 56, 57 de Dittus et Boelter, 77
tenseur des contraintes σ ), 30 de Haaland, 76
tenseur des contraintes de, 30 de Morrisson, 77
Champ de vitesse, 20, 23 de Newton, 77
Churchill et Rose (corrélation de), 77 de Schiller et Nauman, 77
Cisaillement, 26 de Sieder et Tate, 77
Coecient Couple de forces, 86
adiabatique γ , 40 Courant
d'échange de chaleur h, 65 ligne de, 61
de compressibilité κ, 39 tube de, 21
de diusion moléculaire, 43 Courant (ligne de), 21
de dilatation α, 39 Courbe maîtresse, 5
de perte de charge de Darcy fD , 76
de perte de charge de Fanning fF , 76 Déformation, 26
de traînée, 11, 77 angulaire, 26
119
Index
Dérivée de continuité, 54
matérielle (D/Dt), 17 de la chaleur, 63
particulaire (D/Dt), 17 de Laplace, 46
partielle (∂/∂t), 17 de Navier-Stokes, 11, 57
totale (d/dt), 17 de Stokes, 58
Damköhler (nombres de), 11 Equation d'état, 39
Darcy (coecient de), 76 de Tait-Murnaghan, 40
Densité d'un gaz parfait, 40
de ux ϕ, 52 Eulérienne (formulation), 17
volumique f , 51 Euler
volumique de source s, 52 équation de, 58
Diagramme de Moody, 76 nombre de, 11
Diusion théorème de, 72
de chaleur, 42 Expansion, 26
de matière, 42
de quantité de mouvement, 44 Fanning (coecient de), 76
Dimension, 2 Fick (loi de), 43
Dirichlet (conditions aux limites), 64 Fluide
Dittus et Boelter (corrélation de), 77 complexe, 44
incompressible, 40, 47
Echelle newtonien, 44, 57
caractéristique, 3, 10 parfait, 45
changement de, 7 Fluvial (écoulement), 12
Ecoulement Flux
externe, 45 convectif ϕconv , 52
uvial, 12 diusif ϕdi , 52
incompressible, 47 total ϕ, 52
inertiel, 48 Force, 31, 71
interne, 45 Formulation
iso-volume, 46, 54 eulérienne, 17
isochore, 46 lagrangienne, 17
laminaire, 12, 47 Fourier (loi de), 42
rotationnel, 46 Froude (nombre de), 11
solénoïdal, 46
stationnaire, 46 Gaz parfait, 41
subsonique/supersonique, 12, 46 Gradient (théorème du), 91
torrentiel, 12 Grandeur
tourbillonaire, 24 de base, 3
turbulent, 12, 47 physique xi , 1
Einstein (convention de sommation), 79
Haaland (corrélation de), 76
Electro-magnétique (système), 87
Horizontal/vertical (assemblage), 83
Emission (ligne d'), 21
Energie Incompressible
cinétique, 59 écoulement, 47
interne, 61 uide, 40, 47
mécanique, 59 Invariants, 32
potentielle, 59 Invariants (d'un tenseur), 32
totale, 61 Iso-volume (écoulement), 46, 54
Enthalpie h, 63, 72 Isochore (écoulement), 46
Equation
aux dimensions, 2 Kronecker (symbole de), 79
constitutive, 41
d'Euler, 58 Lagrangienne (formulation), 17
de Cauchy, 56, 57 Laminaire (écoulement), 12, 47
120
Index
121
Index
Unité, 1
122
Ecole Nationale Supérieure des Mines de Saint-Etienne
Travaux Dirigés
124
Chapitre 9
Travaux dirigés
Sommaire
9.1 Exercice d'application sur la cinématique . . . . . . . . . . . . . . . . . . . . . . 127
9.2 Convection libre et forcée simultanées . . . . . . . . . . . . . . . . . . . . . . . . 128
9.3 Dimensionnement d'une canalisation cylindrique . . . . . . . . . . . . . . . . . 129
9.4 Résistance à l'avancement d'un navire . . . . . . . . . . . . . . . . . . . . . . . . 130
9.5 Maquette d'un camion en souerie . . . . . . . . . . . . . . . . . . . . . . . . . 131
9.6 Transfert de chaleur d'un lm ruisselant vers une surface solide froide . . . . 132
9.7 Echauement dans un viscosimètre de Couette . . . . . . . . . . . . . . . . . . 134
9.8 Ecoulement de deux uides immiscibles adjacents . . . . . . . . . . . . . . . . . 135
9.9 Prol de température dans une ailette de radiateur . . . . . . . . . . . . . . . 136
9.10 Refroidissement par transpiration . . . . . . . . . . . . . . . . . . . . . . . . . . 138
9.11 Prol de température dans une résistance électrique . . . . . . . . . . . . . . . 139
9.12 Ecoulement d'un uide de Bingham en conduite . . . . . . . . . . . . . . . . . 140
9.13 Tuyère et aube . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
9.14 Rendement d'une éolienne . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
9.15 Obstacle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
9.16 Absorbeur-neutralisateur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144
9.17 Solidication d'une goutte sphérique par refroidissement . . . . . . . . . . . . 146
9.18 Ressaut hydraulique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 147
9.19 Dimérisation catalytique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148
Chapitre 9. Travaux dirigés
126
9.1. Exercice d'application sur la cinématique
L'objectif de cet exercice, qui compte une dizaine de questions, est de manipuler les grandeurs de base de
cinématique et d'établir un lien fort entre les descriptions lagrangienne et eulérienne.
Question 3 : rappeler la condition dans laquelle les lignes de courant se confondent avec les trajectoires.
Question 4 : considérant cette condition réalisée, calculer la trajectoire d'une micro-particule injectée au
point (0,L/2).
127
Chapitre 9. Travaux dirigés
N u = f (Re, P r, Gr)
On envisage, pour déterminer cette fonction f , de réaliser une expérimentation sur une maquette réduite au
dixième (k ≡ L2 /L1 = 0, 1). On opérera avec les mêmes T , ∆T et g . On considérera le même uide (l'air)
mais, si nécessaire, avec une pression et une vitesse diérentes. On supposera que les grandeurs µ, β , cp et
λ restent constantes. L'air est assimilé à un gaz parfait.
Question : indiquer dans quelles conditions de vitesse et pression, on mènera l'expérimentation et comment
on en déduira les pertes thermiques (énergie calorique perdue par unité de temps) de l'installation réelle.
On aectera l'indice 1 à l'installation réelle et l'indice 2 à la maquette.
128
9.3. Dimensionnement d'une canalisation cylindrique
Question 2 : En utilisant l'équation ci-dessous, quelle diérence de pression entre l'amont et l'aval faut-il
imposer pour obtenir ce débit ?
L 1 2
∆p = fD . . ρV (9.1)
d 2
Question 3 : mêmes questions pour une conduite en acier devant transporter un sirop. Pour l'application
numérique, on prendra : Q̇V,2 = 10 l.s−1 , ρ2 = 103 kg.m−3 , µ2 = 10−2 Pa.s, d2 = 40 mm, L2 = 1 m, ε2 = 50
µm.
Question 4 : en admettant que la perte de charge ∆p est proportionnelle à la longueur de la conduite pour
un écoulement incompressible, établir l'équation 9.1 en utilisant l'analyse dimensionnelle. En particulier, on
mettra en oeuvre le contenu des sections "Former des grandeurs adimensionnées" (page 3) et "Adimensionner
une loi physique inconnue" (page 5).
Question 5 : comparer l'écoulement dans les deux conduites. Rééchir à la possibilité d'utiliser la deuxième
(sirop, L = 1 m) comme maquette de la première (eau, L = 100 m). En d'autres termes, comment prédire
la perte de charge ∆p1 de la grande conduite à partir de la mesure expérimentale de ∆p2 faite sur la petite
conduite, sans utiliser le diagramme de Moody ?
129
Chapitre 9. Travaux dirigés
Question I-1 : conduire l'analyse dimensionnelle du problème : identier les nombres sans dimension
pertinents et établir la forme de la fonction qui les lie.
Question I-2 : pour établir la loi qui prédit la résistance à l'avancement du navire, on veut se servir d'une
maquette plus petite d'un facteur λ = 10. Ecrire les conditions de similitude correspondantes. En déduire la
vitesse à laquelle la maquette doit se déplacer. La similitude peut-elle être complète ?
Question I-3 : dans l'expérience sur maquette réalisée au bassin des carènes, le uide utilisé est néces-
sairement de l'eau. Du fait du résultat obtenu à la question précédente, on utilise l'hypothèse simplicatrice
suivante. En première approximation, la résistance à l'avancement R peut se décomposer en deux forces
indépendantes R = Rf + Ro avec :
Rf : force liée à la viscosité du uide et traduisant l'eet du frottement de l'eau le long de la coque
Ro : force représentant la résistance d'onde, c'est-à-dire la résistance opposée par le système de vagues
Question I-4 : en tenant compte du fait qu'il existe de nombreuses formules permettant de calculer la
résistance de frottement, quelle sera alors la marche à suivre pour déterminer la résistance totale qui s'exerce
sur un navire à partir d'une expérience sur maquette ?
Question II-1 : On considère maintenant l'application au cas d'un transbordeur rapide sur l'eau (viscosité
dynamique µ = 10−3 ). Calculer la résistance à l'avancement R1 du transbordeur réel à la vitesse de service V1
de 24 noeuds (c'est-à-dire 12,35 m.s−1 ). Ce bateau a les caractéristiques suivantes : longueur à la ottaison
L1 = 155 m et surface mouillée S1 = 4214 m2 . Une maquette au 1/21 de ce transbordeur a été réalisée et a
fait l'objet d'essais au bassin des carènes. Les résultats expérimentaux obtenus pour la résistance totale R2
en fonction de la vitesse V2 sont indiqués dans le tableau ci-dessous. On donne l'expression suivante pour la
Question II-2 : quelle est la puissance développée pour propulser le bateau à cette vitesse ?
130
9.5. Maquette d'un camion en souerie
On utilise une souerie pour caractériser un attelage "camion porteur + remorque". Il apparaît que la force
de traînée FD dépend de la surface frontale A du véhicule, de sa vitesse V ainsi que des propriétés de l'air.
Question 2 : expliquer comment prédire la force de traînée à l'échelle 1 :1 en utilisant les résultats des
essais en souerie.
faut choisir pour les essais en souerie. Sachant que la force de traînée mesurée en souerie est de 2,46 kN,
calculer la puissance du moteur du camion (en kW).
131
Chapitre 9. Travaux dirigés
9.6 Transfert de chaleur d'un lm ruisselant vers une surface solide
froide
Dans de nombreux contextes industriels, il est nécessaire de refroidir un liquide avant de le ré-utiliser (cf.
centrale thermique à vapeur pour produire de l'électricité) ou de le relâcher dans la Nature (cf. préservation
de la faune aquatique / uviale). Ce processus sera d'autant plus rapide que l'aire de l'interface entre le
liquide chaud et la source de froid sera plus importante. C'est ce qui explique que de nombreuses solutions
consistent à disperser un liquide sous forme de nes gouttelettes dans un gaz plus froid (tour aéroréfrigérante)
ou bien de l'étaler en lm mince sur des parois solides refroidies par ailleurs. Dans cet exercice, on va étudier
le transfert de chaleur dans cette deuxième conguration : un lm de liquide chaud ruisselle en lm mince
sur une surface solide froide. De cette manière, le liquide est refroidi.
Partie I - Question 1 : établir le prol de vitesse à travers un lm liquide ruisselant le long d'une paroi
plane verticale. Expliciter les hypothèses faites et la méthodologie utilisée (s'inspirer de la page 116). Dans
cette question, on considère la température comme uniforme. On utilisera les notations données sur la gure.
Partie I - Question 2 : quelle est la force F~ qu'exerce le uide sur la paroi solide (prendre une aire de 1
m )?
2
Partie I - Question 3 : calculer la vitesse moyenne < vz > en fonction de l'épaisseur δ du lm liquide.
Partie I - Question 4 : proposer un critère permettant de prédire les conditions pour lesquelles l'hypothèse
selon laquelle ~v = vz (x).~ez ne devient plus valable. En déduire une épaisseur de lm maximale dans le cas
de l'eau.
Si la hauteur H du lm est susamment petite (tout en étant grande devant δ ), la chaleur provenant de
la paroi solide ne diuse que sur une petite distance dans le lm liquide. Dans cette sous-région du lm,
donner une approximation du prol de vitesse puis établir l'équation aux dérivées partielles dont la résolution
permettra de déterminer le prol de température. Justier les hypothèses permettant d'arriver à une EDP
sous la forme de l'égalité de deux produits "M.N=P.Q".
Partie II - Question 2 : proposer trois conditions aux limites raisonnables pour résoudre l'EDP établie.
Les traduire sous forme d'équations.
132
9.6. Transfert de chaleur d'un lm ruisselant vers une surface solide froide
Partie II - Question 3 : T −T
on pose θ ≡ T0 −Tff et η ≡ (9β.z) x
1/3 et β ≡ ρ2 .gµ.λ
z .c̄p .δ
. Montrer que l'EDP
peut se mettre sous la forme d'une EDO faisant intervenir la fonction θ(η) et transformer ses conditions aux
limites.
Partie II - Question 5 : donner l'expression du ux de chaleur ~q de la paroi vers le uide et calculer sa
valeur moyenne sur la hauteur H .
133
Chapitre 9. Travaux dirigés
On considère un uide newtonien incompressible placé entre deux cylindres verticaux et coaxiaux. Le cylindre
extérieur, de rayon Rext , est immobile alors que le cylindre intérieur, de rayon Rint ≡ κ.Rext avec κ < 1, est
animé d'un mouvement de rotation de vitesse angulaire ω autour de son axe (voir gure ci-dessous).
Question 1 : déterminer le prol des vitesses dans le uide. On supposera que l'écoulement de uide entre
les deux cylindres est laminaire.
134
9.8. Ecoulement de deux uides immiscibles adjacents
Deux liquides incompressibles et immiscibles s'écoulent selon la direction z dans une ne fente horizontale
de longueur L et de largeur 2b sous l'eet d'un gradient de pression ∆p ≡ (p0 − pL )/L > 0. Les débits sont
ajustés de telle manière que la fente est remplie à moitié par le uide I et à moitié par le uide II. Pour
éviter l'instabilité de Rayleigh-Taylor, on dispose le uide le plus dense en bas. Les vitesses d'écoulement sont
susamment faibles pour que l'écoulement soit stable (ni instabilité de cisaillement de Kelvin-Helmholtz, ni
turbulence,. . .) et donc que l'interface reste parfaitement plane.
Question : déterminer les champs de pression et de contrainte en l'absence de gravité. En déduire le prol
de vitesse lorsque les uides sont newtoniens.
135
Chapitre 9. Travaux dirigés
Question 1 : rappeler les diérents modes de transfert de la chaleur et les équations constitutives associées.
Question 2 : identier l'équation bilan pertinente pour calculer le prol de température. Proposer une
stratégie de résolution. Indiquer pourquoi elle n'est pas pratique à mettre en oeuvre analytiquement.
136
9.9. Prol de température dans une ailette de radiateur
Question 3 : en s'aidant d'une petite analyse dimensionnelle, faire une approximation 0D pour calculer
la puissance thermique dissipée.
137
Chapitre 9. Travaux dirigés
Dans cet exercice, on va regarder comment cette technique d'isolation peut s'appliquer au stockage
cryogénique et, plus précisément, comment elle peut aider à réduire la consommation énergétique nécessaire
pour maintenir de l'azote liquide (en bleu) à la température Tint = −200°C, alors que le réservoir, supposé
sphérique, est au contact de l'atmosphère à température Text = 20°C. La résolution du problème avec une
enveloppe poreuse épaisse allant de r = Rint à r = Rext fait intervenir les écoulements dans les milieux
poreux et dépasse le cadre de ce cours d'introduction à la mécanique des uides. Pour cette raison, nous
allons étudier un dispositif où l'enveloppe poreuse épaisse est remplacée par un espace rempli de gaz et bornée
par deux enveloppes (poreuses) très minces en r = Rint et en r = Rext . On note κ ≡ Rint /Rext ∈ [0; 1].
Que l'isolation soit passive ou active, la chaleur qui vient de l'atmosphère (chaude) vers l'azote (froid) est
évacuée par un serpentin de réfrigération.
Question 1 : les deux enveloppes sphériques sont poreuses. De l'air sec initialement à Tint est soué depuis
le coeur de l'enveloppe interne vers l'extérieur avec un débit massique Wr [kg.s−1 ]. En mettant en oeuvre
la méthodologie donnée dans la section 8 en page 116, évaluer l'inuence de l'écoulement sur le prol de
température. Si besoin, on pourra utiliser les valeurs numériques suivantes pour estimer l'ordre de grandeur
des diérents termes : Rext = 10 cm, ρ = 1, 3 kg.m−3 , κ = 12 et une vitesse typique inférieure à 10 mm.s−1 .
On pourra aussi prendre c̄p = 1000 J.kg−1 .K−1 et λ = 0, 026 W.m−1 .K−1 .
Question 2 : en fonction du débit d'air sec, calculer la puissance de réfrigération Q̇ qu'il faut établir pour
que la température de l'enveloppe interne reste Tint . Comparer à la puissance Q̇0 calculée dans le cas où le
uide est au repos.
138
9.11. Prol de température dans une résistance électrique
Question : identier l'équation bilan pertinente pour calculer le prol de température. Proposer une
stratégie de résolution et la mettre en oeuvre.
139
Chapitre 9. Travaux dirigés
Question 1 : dessiner les contraintes σθr , σzr et σzθ . Représenter par une èche la force qui s'exerce sur
les surfaces élémentaires associées. Se reporter à la section 2.5 en page 30 pour les notations.
Question 2 : identier les équations bilans permettant de calculer le champ de pression p(r,θ,z) et le champ
de contrainte de cisaillement τ(r,θ,z) .
Question 3 : identier les simplications possibles, résoudre et tracer le prol de contrainte de cisaillement.
Calculer sa valeur maximale τmax .
Question 4 : utiliser le résultat précédent pour calculer le prol de vitesse dans le tube lorsque le uide
est décrit par le modèle de Bingham :
pour τ ≤ τc : γ̇ = 0
pour τ ≥ τc
: γ̇ = τ −τc
ou bien τ = µ.γ̇ + τc
µ
Question 5 : vérier que l'on obtient le prol de vitesse d'un uide newtonien en prenant τc = 0.
140
9.13. Tuyère et aube
Question 2 : On suppose que l'aube s'éloigne de la tuyère de manière rectiligne (en restant sur l'axe
de la tuyère) à une vitesse ~vroue . Donner l'expression de la force F~/roue qui s'exerce sur l'aube de la roue.
Celle-ci comporte deux contributions : la force F~eau/roue exercée par le jet liquide et la force F~air/roue due
à la pression atmosphérique. Déterminer la direction et l'intensité de la force totale. Donnée numérique :
k~vroue k = 10 m.s−1
Question 3 : Calculer la puissance P reçue par la roue à aube, celle qui peut être récupérée sur son axe.
Comparer cette puissance à la puissance cinétique PEc du jet (son débit d'énergie cinétique). Optimiser
l'installation.
141
Chapitre 9. Travaux dirigés
Question 2 : On utilise les indices 1 et 2 pour désigner un plan à l'amont lointain et un plan à l'aval
lointain. On utilise les indices A et B pour désigner les plans juste avant et juste après le disque capteur
(l'éolienne). Montrer que la vitesse de l'air varie continument au passage de l'hélice. Montrer également
que le tube de courant qui s'appuie sur le cercle décrit par l'extrémité des pales de l'éolienne s'élargit. On
note v1 la vitesse du vent non perturbé par l'éolienne. Il s'avère avantageux de découper le système en trois
sous-systèmes : 1-A, A-B et B-2.
Question 3 : exprimer la force F~f /h exercée par le uide sur l'hélice en fonction de S1 , v1 et v2 .
Question 6 : Proposer une dénition pour le rendement cp d'une éolienne ? Trouver vh qui maximise la
puissance reçue. Quel est alors la vitesse v2 loin en aval de l'éolienne ?
142
9.15. Obstacle
9.15 Obstacle
Un jet d'eau vient impacter une plaque. La plaque est prolée en son centre pour aider l'eau à prendre
le virage. L'eau quitte la plaque perpendiculairement à sa direction initiale, en respectant une symétrie
cylindrique. Quelle(s) force(s) subit la plaque ? Application numérique : v = 5 m.s−1 , djet = 1 cm, µ = 10−3
Pa.s et R1 = 10 cm.
143
Chapitre 9. Travaux dirigés
9.16 Absorbeur-neutralisateur
Un absorbeur-neutralisateur (scrubber en anglais) est un dispositif destiné à piéger des particules en sus-
pension dans un gaz en les mettant en contact avec des gouttelettes de liquide. Ce type de dispositif permet
de nettoyer un euent gazeux avant de le rejeter dans l'atmosphère ou l'envoyer dans un autre équipement
industriel. L'eau polluée par les particules "nocives" est récupérée puis traitée. Parmi les diérentes techno-
logies, on trouve le système ci-dessous.
Un jet d'eau liquide propre à grande vitesse (vitesse vL et débit massique ṁL dans tube de rayon RL ) est
injecté en (1). Il se fragmente en minuscules gouttelettes qui se dispersent dans le gaz en se déplaçant dans
le tube. En (2), ces gouttelettes et le gaz entraîné se déplacent à la même vitesse v2 , vitesse par ailleurs
approximativement uniforme sur la section (2). La présence du jet liquide (1) a pour conséquence d'aspirer
le gaz en (0), entre les deux disques.
Pour simplier, on considère que le rayon Rd de ces disques est grand devant le rayon Rt du tube horizontal.
On suppose aussi que la pression est uniforme dans chaque section du tube. La pression p1 dans le liquide
injecté est en particulier aussi celle du gaz qui arrive de l'espace entre les disques en Rt . Cela signie que l'on
ne tient pas compte de l'eet du virage pris par ce dernier. On néglige l'eet de la pesanteur et on étudie le
régime stationnaire. L'objectif de l'exercice est d'établir la relation suivante, donnant le débit massique de
144
9.16. Absorbeur-neutralisateur
gaz entraîné :
ṁG = f (vL , ṁL , Rt , Rd , h, ρG , ρL )
Question 1 : établir une relation entre la pression p1 et le débit de gaz ṁG en s'intéressant à l'écoulement
du gaz entre les deux disques (pour r compris entre Rt et Rd ). Indiquer la raison pour laquelle on ne peut
pas utiliser un bilan macroscopique de quantité de mouvement pour établir cette relation. Sans les vérier,
on explicitera les hypothèses faites.
Question 3 : combiner les deux résultats pour exprimer ṁG en fonction de ṁL , h, RL , ρG , ρL .
Question 4 : vérier a posteriori si les hypothèses adoptées dans les questions 1 et 2 sont valables. On
prendra les valeurs numériques suivantes : Rt = 10 cm, Rd = 50 cm, h = 1 cm, ṁL = 4 kg.s−1 , vL = 40
m.s−1 , ρL = 1000 kg.m−3 , ρG = 1, 3 kg.m−3 , µG = 2.10−5 Pa.s, p0 = patm = 105 Pa et c = 330 m.s−1 .
145
Chapitre 9. Travaux dirigés
Question 2 : on se place dans le cas où le transfert par conduction est limitant (= le transfert de chaleur
à l'interface sphère-gaz est comparativement beaucoup plus rapide). Ceci permet d'adopter une stratégie
de résolution en deux étapes. Dans un premier temps, on considère donc que le régime est stationnaire.
Ecrire l'EDP permettant d'établir le prol de température dans la zone solide, entre la surface située en R et
l'interface liquide-solide en Req , puis l'intégrer (en considérant que la température de surface est constante).
Question 4 : montrer que le ux de chaleur QR peut se mettre sous la forme QR = H.S. (T0 − Tair ) où
H est le coecient d'échange global en [W.m−2 .K−1 ]. Expliciter
H en fonction de h , R et Req .
1 1
Question 5 : dans un deuxième temps, on considère connu le ux de chaleur QR quittant la sphère. Donner
−dR
une expression permettant de calculer la vitesse de propagation de l'interface solide-liquide dt eq en m.s−1 .
On notera L la chaleur latente de solidication en J.kg .
−1
146
9.18. Ressaut hydraulique
La gure de gauche ci-dessous représente un ressaut hydraulique obtenu en dirigeant un jet liquide vertical
sur un plan horizontal. En régime stationnaire, on observe un bourrelet de faible dimension qui sépare la
région centrale où l'épaisseur du lm liquide est faible de la région périphérique où la hauteur du lm est plus
importante. A droite, on observe une coupe de ce ressaut. La suite de cet exercice se fera dans un référentiel
cartésien.
Question 3 : les conditions pour écrire le théorème de Bernoulli entre les sections (1) et (2) sont-elles
réunies ? Justier.
Question 6 : montrer que, sur chaque section (1) ou (2), la pression généralisée est indépendante de la
profondeur. Comparer ces deux pressions généralisées.
Question 7 : pensez-vous qu'un écoulement de même géométrie pourrait exister avec la vitesse orientée
dans l'autre direction ? Justier.
147
Chapitre 9. Travaux dirigés
Dans cet exercice, on considère que le régime est stationnaire et les transferts de chaleur sont négligeables
(bien qu'en réalité les réactions solide-gaz soient souvent très rapides et exothermiques). Le lm est supposé
au repos et les constituants à l'état gazeux. La gure de gauche ci-dessous représente un lit catalytique et
celle de droite les notations.
Question 1 : exprimer le ux massique de A qui est converti en B en fonction de l'épaisseur δ de la couche
de diusion. On suppose ici que la réaction chimique se déroule instantanément et que la fraction molaire
de A à la limite de la zone de diusion est égale à xA0 . Cette valeur de xA0 décroît quand la position y dans
le lit catalytique s'éloigne de l'entrée.
Question 2 : reprendre la question 1 en considérant que la réaction chimique n'est pas inniment rapide
mais obéit à une loi cinétique de premier ordre. A l'interface solide-uide, le ux massique de consommation
de A s'exprime par :
ϕA = νA .M̃A .rA
350°C
où νA = −2 est le coecient stoechiométrique de Ag dans la réaction 2Ag −−−−→ Bg et M̃A sa masse molaire.
métal
La vitesse de réaction rA est toujours positive et donnée par :
rA = k1 .cA avec cA = xA .c
Puisque la réaction est surfacique, ϕA est en [kg.m−2 .s−1 ], rA est en [mol.m−2 .s−1 ] et k1 en [m.s−1 ] 1 .
1. Dans le cas usuel d'une réaction chimique , on écrit la même équation ϕA = M̃A .k1 .cA mais ϕA s'exprime en
volumique
[kg.m−3 .s−1 ] et k1 en [s−1 ].
148