Corrections

November 10, 2025

[1]: import numpy as np

import [Link] as plt
%matplotlib widget

1 Exercice 4.1
Débit de vinaigre.
Le débit volumique 𝑄 d’un fluide coulant avec une vitesse uniforme 𝑈 dans une conduite de section
𝑆 vaut 𝑄 = 𝑈 𝑆.
2
Attention, la section d’un tube circulaire de diamètre 𝜙 vaut 𝑆 = 𝜋 𝜙4 . Il faut aussi penser à utiliser
les unités S.I. (m, s, kg)

[2]: phia=150e-3 #en m

sa=[Link]*phia**2/4 #en m2
phib=75e-3 #en m
sb=[Link]*phib**2/4 #en m2
q=100/1000 #en m3/s
va=q/sa #en m/s
vb=q/sb #en m/s
print("Vitesse dans la section de diamètre 150mm V={:.2f} m/s".format(va))
print("Vitesse dans la section de diamètre 75mm V={:.2f} m/s".format(vb))

Vitesse dans la section de diamètre 150mm V=5.66 m/s

Vitesse dans la section de diamètre 75mm V=22.64 m/s
Quand on divise la section par 2, on multiplie la vitesse par 4. La deuxième valeur (22.6 m.s−1 ) est
irréaliste!

2 Exercice 4.2
On change un peu l’énnoncé: le profil de vitesse est tel que:

4𝑉0
𝑣= 𝑧 (ℎ − 𝑧)
ℎ2

[3]: V0=1
h=1

1
 z=[Link](0,h,100)
v=4*V0*h**-2*z*(h-z)
fi,ax=[Link](1,1)
[Link](v/V0,z/h,label="Profil")
ax.set_xlim(left=0,right=1.2)
ax.set_ylim(bottom=0,top=1)
[Link]()
ax.set_xlabel(r"$\frac{v}{V_0}$")
ax.set_ylabel(r"$\frac{z}{h}$")
[Link]([2/3,2/3],[0,1],label="Vmoy")
[Link]()

[3]: <[Link] at 0x783e59744560>

Le profil est parabolique.

Les vitesses sont nulles sur les parois. C’est compatible avec un écoulement de fluide réel. En fait,
ce profil est une solution de l’équation de Navier-Stokes en régime permanent, laminaire.
La vitesse maximale est atteinte en 𝑧 = ℎ/2 par symétrie.
Pour calculer le débit, il faut sommer sur tous les éléments de surface élémentaire 𝑑𝑆 = 𝐿𝑑𝑧:

2
 ℎ
𝑄 = ∫ 𝑣𝑑𝑆 = ∫ 𝑣(𝑧)𝐿𝑑𝑧
𝑆 0

ℎ
𝑄 = ∫ 𝑣(𝑧)𝐿𝑑𝑧
0

𝑧
4𝑉0
𝑄=∫ 𝑧 (ℎ − 𝑧) 𝐿𝑑𝑧
0 ℎ2

𝑧
4𝑉 𝐿
𝑄 = 02 ∫ 𝑧 (ℎ − 𝑧) 𝑑𝑧
ℎ 0

𝑧
4𝑉 𝐿 ℎ𝑧2 𝑧3
𝑄 = 02 [ − ]
ℎ 2 3 0

4𝑉0 𝐿 3 1 1
𝑄= ℎ ( − )
ℎ2 2 3

2
𝑄 = 𝑉0 ℎ𝐿
3

La vitesse moyenne débitante vaut 2/3 de 𝑉0 .

3 Exercice 4.3
C’est pareil mais en cylindrique: l’élément d’intégration vaut 2𝜋𝑟𝑑𝑟. Le résultat est que la vitesse
moyenne est cette fois-ci la moitié de la vitesse maximale au centre de la conduite.
Cet exercice sous une version étendue est corrigé dans le document “Demo…pdf”, en section 5: à
lire ! Il est très important pour la compréhension de nombreuses notions.

4 Exercice 5.1
On utilise le théorème de Bernoulli entre 𝐴 et 𝐵 (sans pertes):

𝑝𝑡𝑜𝑡𝑎𝑙𝑒,𝐴 = 𝑝𝑡𝑜𝑡𝑎𝑙𝑒,𝐵

1 1
𝑝𝐴 + 𝜌𝑉𝐴2 + 𝜌𝑔𝑧𝐴 = 𝑝𝐵 + 𝜌𝑉𝐵2 + 𝜌𝑔𝑧𝐵
2 2

1
𝑝𝐵 = 𝑝𝐴 + 𝜌 (𝑉𝐴2 − 𝑉𝐵2 ) + 𝜌𝑔 (𝑧𝐴 − 𝑧𝐵 )
2

3
[4]: q=350/1000 #en m3/s
DA=35e-2 #en m
sa=[Link]*DA**2/4 #en m2
va=q/sa
print("Vitesse dans la section de diamètre 350mm V={:.2f} m/s".format(va))
DB=64e-2
sb=[Link]*DB**2/4 #en m2
vb=q/sb
print("Vitesse dans la section de diamètre 640mm V={:.2f} m/s".format(vb))
g=9.81
rho=1000
za=3
zb=7.5
delta_dyn=0.5*rho*(va**2-vb**2)
print("Delta pression dynamique={:.3f} bar".format(delta_dyn/1e5))
delta_stat=rho*g*(za-zb)
print("Delta pression hydrostatique={:.3f} bar".format(delta_stat/1e5))
print("P en b = P en a + {:.3f} bar".format((delta_stat+delta_dyn)/1e5))
# pdc avec kelargissement=0.1, kcoude=0.05 et L=10:
pdc=0.5*rho*vb**2 * (10/0.64 + 0.1 + 0.05)
print("Chute de pression totale irréversible = {:.3f} bar".format(-pdc/1e5))

Vitesse dans la section de diamètre 350mm V=3.64 m/s

Vitesse dans la section de diamètre 640mm V=1.09 m/s
Delta pression dynamique=0.060 bar
Delta pression hydrostatique=-0.441 bar
P en b = P en a + -0.381 bar
Chute de pression totale irréversible = -0.093 bar
Le ralentissement du à l’augmentation de section augmente la pression, mais l’effet d’élévation dans
le champ de pesanteur la fait baisser beaucoup plus. La perte de pression totale due au coude, à
l’élargissement et au frottement pourrait valoir en ordre de grandeur (pour un tube de 10 m) 1 à
1.5 fois l’effet dynamique et annuler celui-ci.

5 Exercice 5.2
Théorème de Bernoulli entre 0 et 1 (on place la référence d’altitude n’importe où):

1 1
𝑝0 + 𝜌𝑔𝑧0 + 𝜌𝑉02 = 𝑝1 + 𝜌𝑔𝑧1 + 𝜌𝑉12 (1)
2 2

Théorème de Bernoulli entre 0 et 2:

1 1
𝑝0 + 𝜌𝑔𝑧0 + 𝜌𝑉02 = 𝑝2 + 𝜌𝑔𝑧2 + 𝜌𝑉22 (2)
2 2

Conservation du débit :

4
 𝑆0 𝑉0 = 𝑆1 𝑉1 + 𝑆2 𝑉2

Mais 𝑆1 = 𝑆2 = 1/2𝑆0 , donc

(𝑉1 + 𝑉2 ) = 2𝑉0 (3)

Si on fait la différence des 2 équations de Bernoulli (2)-(1):

1
0 = (𝑝2 − 𝑝1 ) + 𝜌𝑔 (𝑧2 − 𝑧1 ) + 𝜌 (𝑉22 − 𝑉12 )
2
Mais 𝑝1 = 𝑝2 , et 𝑧1 − 𝑧2 = ℎ

1 2
(𝑉 − 𝑉12 ) = 𝑔ℎ
2 2

(𝑉2 − 𝑉1 ) (𝑉2 + 𝑉1 ) = 2𝑔ℎ

𝑔ℎ
(𝑉2 − 𝑉1 ) = (4)
𝑉0

Reste à résoudre le système (3) et (4). (3)+(4) donne:

𝑔ℎ
2𝑉2 = 2𝑉0 +
𝑉0

𝑔ℎ
𝑉2 = 𝑉0 (1 + )
𝑉02

(3)-(4) donne:

𝑔ℎ
𝑉1 = 𝑉0 (1 − )
𝑉02

𝑉1 est plus petite que 𝑉2 , et on vérifie que si ℎ = 0, il ne se passe rien (𝑉1 = 𝑉2 = 𝑉0 ).

6 Autres:
Ils sont corrigés dans le document “Demo…pdf” A LIRE !!!!