54 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

De là, quand aD0.0025 m et D1ð10 3Ns/m2 ,

kD6,24 Nsn /m2

PROBLÈME 3.39
Deux liquides de densités égales, l'un newtonien et l'autre non newtonien.
loi" fluide, s'écoulant à des débits volumiques égaux le long de deux grandes surfaces verticales de la même
largeurs. Le fluide non newtonien a un indice de loi de puissance de 0,5 et a la même apparente
la viscosité, en unités SI, en tant que fluide newtonien lorsque son taux de cisaillement est de 0,01 s 1 . Montrez que,
pour des vitesses de surface égales des deux fluides, l'épaisseur du film pour le fluide newtonien
est 1,125 fois celle du fluide non newtonien.

Solution
Pour un fluide de loi de puissance :

n
de
RDk (équation 3.121)
dy
n 1
du dex
Dk D dux /dy (équation 3.122)
dy dy a

oùunest la vitesse apparente Dk dux /dy n 1

Pour un fluide newtonien :
dex
RD (équation 3.3)
dy
WhennD0.5 et dux /dy D0.01, D unet :

Dk dux /dy n 1Dk 0,01 0.5D10kD etkD0.1 .

unD
L'équation d'état du fluide à loi de puissance est donc :
RD0.1 dux /dy 0.5

Pour un fluide s'écoulant le long d'une surface verticale, longueur l, largeur w et épaisseur de film S,
à une distance de la surface solide, un équilibre de forces donne :
S y wl gDRwlDk dux /dy n wl

dex g 1/n
ou : S y
dy D k 1/n

g nC1 n
et : uxD S y n
Const.
k 1/n nC1
g nC1 n
Quand y a-t-il ? xD0 et la constante D S n
K 1/n nC1
g n nC1 nC1
et : uxD S n S y n
k 1/n nC1
 ÉCOULEMENT DANS LES TUYAUX ET CANAUX 55

À la surface libre où yDS :

g n nC1
usD S n (i)
k 1/n nC1
Le débit volumique Q est donné par :
s
g n nC1 nC1
QD wdy S n
s y n
0 k 1/n nC1
g n 2nC1
Dw S n (ii)
k 1/n 2nC1
Pour le fluide non newtonien, kD0,1, nD0,5 et l'équation (ii), lorsqu'elle est exprimée en S.I.
unités, devient :
g 2 gs2
Dw ð0,25 s4D25w (iii)
0.1 2

Pour le fluide newtonien, nD1 et kD et en substituant dans l'équation (ii) :

g
QDw ð0.33S3N (iv)

oùSNc'est l'épaisseur du film newtonien.

Pour des débits égaux, à partir des équations (iii) et (iv) :

gs2 g
25w 2 D0.33w S3N

g
ou : S3ND75 S4 (v)

Pour des vitesses de surface égales, le terme g/K dans l'équation (i) peut être substitué à partir de
équation (v) et :

g 2
Pour le fluide non newtonien : usD 0.1 ð0.33S3

S3N
D100 75S4 2 ð0.33S3

D0.00592S6N/S5
g
Pour le fluide newtonien : usD ð0.5S2N
S3N
D 75S4 ð0.5S2ND0.0067S5N/S4

et : SN /SD 0.00667/0.00592 D1.126

56 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

PROBLÈME 3.40
Un fluide qui présente un comportement non newtonien s'écoule dans un tuyau de diamètre 70 mm
et la chute de pression sur une longueur de 2 m de tuyau est de 4ð104N/m2 Lorsque le débit est
doublée, la baisse de pression augmente d'un facteur de 1,5. Un tube de Pitot est utilisé pour mesurer
le profil de vitesse à travers la section transversale. Confirmez que les informations données ci-dessous sont
conforme à l'écoulement laminaire d'un fluide à loi de puissance.
Toute équation utilisée doit être dérivée de la relation de base entre le stress de cisaillement R et
taux de cisaillement

RDk ,P n

distance radiale du centre vélocité

de tuyau (smm) (m/s)
0 0,80
10 0,77
20 0,62
30 0,27

Solution
Pour un fluide à loi de puissance :

Au débit initial : 4ð104Dkun

Avec un flux de : 6ð104Dk 2u n

Diviser : 1.5D2n
et donc : nD0,585
Pour le fluide de loi de puissance :
n
dex
RDk
dy
Un équilibre de forces sur un noyau fluide de rayon dans un tuyau de rayon donne :

Rs 2 slD P s2
n
dux Ps
ou : RsDk D
ds 2l
1/n
dex P
s1/n
ds D 2kl
1/n
P n nC1
Intégration : uxD s n
Cconstant
2ks nC1
 ÉCOULEMENT DANS LES TUBES ET CANAUX 57

Quand es-tu? xD0 (la condition sans glissement) :

1
P n n nC1
et donc la constante D r n
2kl nC1
1/n
P n nC1 nC1
Substitution : uxD r n
s n
2kl nC1
1/n
P n nC1
Sur la ligne médiane : uCL D r n
2kl nC1
nC1
ux s n
et donc :
uCL D1 r
s 2,71
à quand D0.585 : D1 r
Les données suivantes sont obtenues pour un rayon de tuyau de 35 mm :

expérimental
ux s 2,71
rayon(mm) ux(m/s) ux /uCL
uCL D1 r
0 1 0,80 1
10 0,966 0,77 0,96
20 0,781 0,62 0.77
30 0,341 0,27 0,34

Ainsi, les valeurs calculées et expérimentales de ux /uCLd'accord dans des limites raisonnables de
précision expérimentale.

PROBLÈME 3.41

Un fluide de type Bingham-plastique (tension d'écoulement 14,35 N/m) 2et la viscosité plastique 0,150 Ns/m2 ) est
s'écoulant à travers un tuyau de diamètre 40 mm et de longueur 200 m. En commençant par la rhéologique
équation, montrer que la relation entre le gradient de pression ∂P/l et le débit volumique
Qis
P r4 4 1
QD 1 XCX4
8l p 3 3

où est le rayon du tuyau, pest la viscosité plastique, et X est le rapport de la contrainte d'écoulement
au stress de cisaillement à la paroi du tuyau.
Calculez le débit pour ce pipeline lorsque la chute de pression est de 600 kN/m.2 Cela peut
on doit présumer que l'écoulement est laminaire.
58 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

Solution
Pour un matériau plastique de Bingham, le stress de cisaillementsau radiussis donné par :
dux
R s R YD C R s RY
p ds
dex
ds D0 Rs RY
Le noyau central non écorché a un rayon rcDr RY /R (où rD est le rayon du tuyau et R D
tension de cisaillement au mur) puisque la tension de cisaillement est proportionnelle au rayon.
Dans la région annulaire :
dux 1 1 s
D R RY D P RYd'un équilibre des forces
ds p p 2l
1 s2
u xD dexD P RY sCconstant
p 4l
Pour la condition sans glissement : uxD0, quand est-ce que Dr

1 r2
Ainsi : 0D P RY sCconstant
p 4l
1 P
et : usD r 2 s2 RY r s
p 4l
2R
Substituer : PD
l/r
1 R
usD r 2 s2 RY r s (i)
p2r

Le débit volumique à travers l'anneau élémentaire, dQUnVouss 2 sds

r
1 R
Ainsi : Q AD r 2 s2 RY r s 2 sds
rc p2r
r
2 1 r 2 s2 s4 RYrs2 s3
D R
p 2r 2 4 R 2 3 rc

RY RY
Écriture DXandr Dr, alors
R c R:
2 1 r4 r4 r3 r3 1 X2 r 4 X4 r 4
QUnD R X
p 2r 2 4 2 3 2r 2 4
r 3 X2 r 3 X3
CX 2 3
2 R3 1 1 1 1 1 1
r X 2
XCX 4CX3 4
D p 8 6 4 8 X 2 3
 FLUX DANS LES TUYAUX ET CANAUX 59

Rr3 5
D 4 1 4/3X 2X2C4X3 X4 (ii)
p 3
Dans la région centrale

Substitution : sDrcD RY /R rDXrin équation (i) pour la vitesse du cœur uCdonne :

1 R
ucD
p2r
r 2 X2 r 2 RYr

Rr 1 Rr
1 X2 X 1 X 2 4XC4X2 g
D p 2 4 p f2 1 X
Rr 2
D 4 p 2 4XC2X

Le débit à travers le noyau est : uc rc2Vousc X2 r 2DQc

Rr
Ainsi : QcD 4 X2 r 2 f2 4XC2X2 g
p

Rr3
D 4 p f2X2 4X3C2X4 g

Le débit total est : QUnCQc DQ

Rr3 4 4 1
et QD 1 XCX
3 3
4 p

Pr
Mettre RD alors :
2l
P r4 4 1
QD f1 XCX
3 4g 3
8l p

Quand : PD6ð105N/m2 ,lD200 m dD40 mm et rD0.02 m.

r 0,02
Alors : RD PD6ð 5 2
2l 400 ð10 D30 N/m

pD0.150 Ns/m2
RYD14,35 N/m2
RY 14,35
et : XD
R D 30 D0.478
6ð105 0.02 4 4 1 3
Ainsi : QD 1 ð 0.478C 0.478
8ð200ð0.150 3 3

D0.000503 m3 /s
 SECTION 4

Écoulement des fluides compressibles

PROBLÈME 4.1
Un gaz ayant un poids moléculaire de 13 kg/kmol et une viscosité cinématique de 0,25 cm2 /s
s'écoule à travers un tuyau de 0,25 m de diamètre interne et de 5 km de long à un débit de 0,4 m3 /s
et est livré à pression atmosphérique. Calculez la pression nécessaire pour maintenir cela
débit sous des conditions isothermes. Le volume occupé par 1 kmol à 273 K et
101,3 kN/m2est 22,4 m3 Quel serait l'effet sur la pression requise si le gaz
devaient être livrés à une hauteur de 150 m (i) au-dessus, et (ii) en dessous de son point d'entrée dans
la pipe ?

Solution
À partir de l'équation 4.57 et, comme première approximation, en omettant le terme d'énergie cinétique :

P2 P1 /vmC4 R/ u2 l/d G/A 2D0

À la pression atmosphérique et à 289 K, la densité D⟨13/22.4⟩⟨273/289⟩D 0.542 kg/m 3

Débit massique de gaz, GD 0,4(0,542)D0,217 kg/s.

Aire de section transversale, AD /4 0,25 2D0,0491 m2 .

Vitesse du gaz, uD 0.4/0.0491 D8.146 m/s

∴ G/AD 0.217/0.0491 D4.413 kg/m2 s

Nombre de Reynolds, ReD du/

D 0,25ð8,146/0,25ð10 4 D8.146ð104

De la Fig. 3.7, fore/dD0.002, R/ u2D0.0031

v2D 1/0.542 D1.845 m3 /kg

v1D 22.4/13 298/273 101.3/P1 D190.5/P1m3 /kg
et vmD 0.923P1C95.25 /P1m3 /kg
Substituer dans l'équation 4.57 :

P1 P1 101.3 103 / 0.923P1C95.25 D4 0.0031 5000/0.25 4.726 2

et : P1D111,1 kN/m2
60
 ÉCOULEMENT DES FLUIDES COMPRESSIBLES 61

Le terme d'énergie cinétique D G/A 2ln P1 /P2 D 4.413 2ln 111.1/101.3

D1,81 kg2 /m4 s2

Cela est négligeable par rapport aux autres termes qui équivalent à 5539 kg.2 /m4 s2donc
que l'approximation initiale est justifiée. Si le tuyau n'est pas horizontal, le terme gdzin
l'équation 4.49 doit être incluse dans le calcul. Si l'équation 4.49 est divisée parv2 , cela
2.
terme sur l'intégration devient g z/vm

∴ vmD 0,923ð111,1C95,25 /111,1D1,781 m3 /kg

vairD 24,0/29 D0,827 m3 /kg

Comme le gaz est moins dense que l'air,vmest remplacé parvair vm D 0,954 m3 /kg.
∴ 2 D 9.81ð150/0.9542 D1616 N/m2ou 0,16 kN/m2
g z/vm
(i) Si le point de livraison est à 150 m au-dessus du niveau d'entrée, alors puisque le gaz est moins dense,

P1D 111.1 0.16 D110.94 kN/m2

(ii) Si le point de livraison est situé 150 m en dessous du niveau d'entrée alors,

P1D 111,1C0,16 D111,26 kN/m2

PROBLÈME 4.2
Azote à 12 MN/m2la pression est transférée à travers un tuyau en acier doux de 25 mm de diamètre à un
usine d'ammoniac synthétique à un rythme de 1,25 kg/s. Quelle sera la chute de pression sur
un tuyau de 30 m de longueur pour un écoulement isotherme du gaz à 298 K ? Rugosité absolue
de la surface du tuyau D0.005 mm. Volume moléculaire en kilogrammes D22.4 m3 Viscosité de
azote D0,02 mN/s.m2 .

Solution
Poids moléculaire de l'azote D28 kg/kmol.
En supposant une pression moyenne dans le tuyau de 10 MN/m2 , le volume spécifique,vmà
10 MN/m2et 298 K est :

vmD 22.4/28 101.3/10ð103 298/273 D0.00885 m3 /kmol

Nombre de Reynolds, ud/ Dd G/A / .

AD /4 0.025 2D4.91ð10 3m2 .

∴ G/A D 1.25/4.91ð10 3 D2540 kg/m2 s

et : ReD 0.025ð2540/0.02ð10 3 D3.186

62 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

D'après la Fig. 3.7, pour ReD3.18(10)6ande/dD 0.005/25 D0.0002

R/ u2D0.0017
Dans l'équation 4.57 et en négligeant le premier terme :

P2 P1 /vmC4 R/ u2 l/d G/A 2D0

ou : P1 P2D4vm R/ u2 l/d G/A 2

D4ð0.00885 0.0017 30/0.025 2540 2

D466,000 N/m2ou 0,466 MN/m2

C'est petit en comparaison avec P1D12 MN/m2 , et la pression moyenne de

10 MN/m2est considéré comme trop bas. Une pression moyenne de 11,75 kN/m2est donc sélectionné
et le calcul répété pour donner une perte de pression de 0,39 MN/m2 La pression moyenne
est alors 12C11.61/2D11.8 MN/m2qui est suffisamment proche de la valeur supposée.
Il reste à vérifier si l'hypothèse selon laquelle le terme d'énergie cinétique est négligeable est
justifié.

Terme d'énergie cinétique D G/A 2ln P1 /P2 D 2540 2ln(12/11.61)D2.13×105 kg2/m4s2

Le terme P1 P2 /vm oùvmest le volume spécifique à la pression moyenne de

11,75 MN/m2D 0,39ð106 /0,00753D5,18ð107kg2 /m4 s.
D'où l'omission du terme d'énergie cinétique est justifiée

et la chute de pression D0,39 MN/m2

PROBLÈME 4.3
L'hydrogène est pompé d'un réservoir à 2 MN/m2pression à travers un mild horizontal propre
tuyau en acier de 50 mm de diamètre et 500 m de long. La pression en aval est également de 2 MN/m2
et la pression de ce gaz est élevée à 2,6 MN/m2par une pompe à l'extrémité amont de
le tuyau. Les conditions d'écoulement sont isothermes et la température du gaz est de 293 K.
What is the flowrate and what is the effective rate of working of the pump? Viscosity of
hydrogène D0,009 mN s/m2à 293 K.

Solution
En négligeant le terme d'énergie cinétique dans l'équation 4.55, alors :

P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0

oùP1D2.6 MN/m2etP2D2.0 MN/m2 .

Ainsi : v1D 22,4/2 293/273 0,1013/2,6 D0,468 m3 /kg

IfReD107ande/dD0.001, à partir de la Fig. 3.7,R/ u2D0,0023.

 ÉCOULEMENT DES FLUIDES COMPRIMABLES 63

Substitution :

2.02 2.62 1012/ 2ð2.6ð106ð0.468 C4 0.0023 500/0.05 G/A 2D0

à partir de quel G/AD111 kg/m2 s.

ReDd G/A / D 0.05ð111/ 0.009ð10 3 D6.2ð105

∴

Ainsi, la valeur choisie de Re était trop élevée. Si Re est pris comme 6,0 × 105et le problème
retraité2 s etReD6.03ð105qui est en bon accord.

AD /4 0.05 2D0.00197 m2
et: GD108ð0.00197D0.213 kg/s

La puissance requise est donnée par l'équation 8.71 comme 1/ÁGP1 v1ln P1 /P2
Si une efficacité de 60 % est supposée, alors la puissance requise est :

D 1/0.6 ð0.213ð2.6ð106ð0.468 ln 2.6/2

D 1.13ð105 W ou 113 kW

PROBLÈME 4.4
Dans une usine d'ammoniac synthétique, l'hydrogène est acheminé par un tuyau en acier de 50 mm vers le
convertisseurs. La chute de pression sur la longueur de 30 m de tuyau est de 500 kN/m2 la
pression
à l'extrémité aval étant de 7,5 MN/m2 Quelle puissance est nécessaire pour surmonter
pertes par friction dans le tuyau ? Supposer une expansion isotherme du gaz à 298 K. Quelle erreur
est introduit en supposant que le gaz soit un fluide incompressible de densité égale à cela
à la pression moyenne dans le tuyau ? D0,02 mNs/m2 .

Solution
Si la pression en aval est de 7,5 MN/m2et la chute de pression due à la friction D
500 kN/m2 , la pression amont D8.0 MN/m2et la pression moyenne D7,75 MN/m2 .
Le volume spécifique moyen est :vmD 22.4/2 298/273 0.1013/7.75 D0.16 m3 /kg
et : v1D 22,4/2 298/273 0,1013/8,0 D0,15 m3 /kg
Il est nécessaire d'assumer une valeur de R/ u2 , calculez G/A et le nombre de Reynolds et
vérifiez que la valeur de e/dis est raisonnable. Si le gaz est supposé incompressible
fluide de densité égale à la pression moyenne dans le tuyau et R/ u2D0.003, la pression
chute due à la friction D500 kN/m2est :

∴ 500ð103 /0.16 D4 0.003 30/0.05 G/A 2

et G/AD658 kg/m2 s.
ReDd G/A / D 0.05ð658/0.02ð10 3 D1.65ð106
64 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

D'après la Fig. 3.7, cela correspond à une valeur de e/dof d'environ 0,002, ce qui est
raisonnable pour un tube en acier.
Pour l'écoulement compressible :

G/A 2ln P1 /P2 C P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0 (équation 4.55)
Substitution :

G/A 2ln(8.0/7.5)C(7.5)2 8.02 1012/ 2ð8.0ð106ð0.15

C4 0,003 30/0,05 G/A 2D0

de quel poids : G/AD667 kg/m2 s et GD667ð /4 0.05 2D1,31 kg/s

Une petite erreur est faite par l'hypothèse de simplification dans ce cas particulier.
La puissance requise est donnée par l'équation 8.71 :

D 1/Á GP1 v1lnP1 /P2

Si l'efficacité du compresseur est de 60%,

exigence de puissance D 1/0.6 ð1.31ð8.0ð106ð0.15 ln 8/7.5

D 1.69ð105 W ou 169 kW

PROBLÈME 4.5
Une usine de distillation sous vide fonctionnant à 7 kN/m2la pression au sommet a un taux de montée de
0,125 kg/s de xylène. Calculez la perte de pression le long d'un tuyau de vapeur de 150 mm de diamètre utilisé
pour connecter la colonne au condenseur. La longueur du tuyau peut être considérée comme équivalente à
6 m,e/dD0.002 et D0.01 mN s/m2 .

Solution
À partir des données de pression de vapeur, la température de vapeur D338 K et le poids moléculaire de
xylène D106 kg/kmol.
Dans l'équation 4.55 :

G/A 2ln P1 /P2 C P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0

Aire de section transversale du tuyau, AD /4 0.15 2D1,76ð10 2m2

G/AD 0,125/1,76ð10 2 D7,07 kg/m2 s

Le nombre de Reynolds, est ud/ Dd G/A /

D 0,15ð7,07/ 0,01ð10 3 D1.06ð105

D'après la Fig. 3.7, avec/dD0.002 et ReD1.06ð105 , R/ u2 D0.003.

Volume spécifiquev1D 22.4/106 338/273 101.3/7.0 D3.79 m3 /kg.

 ÉCOULEMENT DES FLUIDES COMPRESSIBLES 65

Substituant dans l'équation 4.55 :

7.0 2ln 7/P2 C P22 72 ð106 /2ð7ð103ð3.79C4ð0.003 6/0.15 7.07 2D0

oùP2est-ce que la pression au condenseur est en kN/m 2 .

Résoudre par essais et erreurs :

P2D6,91 kN/m2
∴ P1 P2 D 7.0 6.91 D0.09 kN/m2ou 90 N/m2

PROBLÈME 4.6
Azote à 12 MN/m2la pression est introduite à travers un tuyau en acier doux de 25 mm de diamètre à un
usine d'ammoniac synthétique à un taux de 0,4 kg/s. Quelle sera la chute de pression sur
une longueur de tuyau de 30 m en supposant une expansion isotherme du gaz à 300 K ? Quelle est le
quantité moyenne de chaleur par unité de surface de tuyau qui doit passer à travers les parois dans
Dans le but de maintenir des conditions isothermes ? Quelle serait la chute de pression dans le tuyau si
il était parfaitement décalé? D0,02 mNs/m 2 .

Solution
À haute pression, le terme d'énergie cinétique dans l'équation 4.55 peut être négligé pour donner :

P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0

Volume spécifique à l'entrée du tuyauv1D 22,4/28 300/273 0,1013/12

D0,00742 m3 /kg
Aire de la section transversale du tuyau, AD /4 0,025 2D0,00049 m2
∴ G/AD 0.4/0.00049 D816 kg/m2 s.
Nombre de Reynolds, d G/A / D0.025ð816/ 0.02ð10 3 D1.02ð106
Ife/dD0.002 et ReD1.02ð106 R/ u2D0.0028 de la Fig. 3.7.
Substitution : 122 P22 1012 / 2ð12ð106ð0.00742 D4 0.0028 30/0.025 816 2

et:P2D11.93 MN/m2
et : chute de pression D 12.0 11.93 D0.07 MN/m2À70 kN/m2

La chaleur nécessaire pour maintenir un flux isotherme est donnée dans la Section 4.5.2 comme G u.2 /2.
La vitesse à l'extrémité haute pression du tuyau Ddébit volumétrique/surface

D G/A v1D 816ð0.0072 D6.06 m/s

et la vitesse dans l'usine est considérée comme nulle.

Ainsi : G u2 /2D0.4ð 6.06 2 /2D7.34 W

66 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

Zone extérieure du tuyau D 30ð ð0.025 D2.36 m2 .

Chaleur requise D 7.34/2.36 D3.12 W/m2
Cette faible valeur de la chaleur requise provient du fait que le changement d'énergie cinétique
est petit et les conditions sont presque adiabatiques. Si le tuyau était parfaitement isolé, le flux
serait adiabatique et la chute de pression serait alors calculée à partir des équations 4.77
et 4,72. Le volume spécifique à l'extrémité basse pressionv2à calculer à partir de :
2 2
1P1 A v1 C1 v
8 R/ u2 l/d D C v 1 dans 2
2 1 G v2 v1
(équation 4.77)
Pour l'azote, D1.4 et donc :
2 2
1.4 1 12 ð106 1 0,00742
8 0.0028 30/0.025 D C 1
2ð1.4 0,00742 816 v2
1.4C1 v2
ln
1.4 0,00742

Résoudre par essai et erreur,v2D0,00746 m3 /kg.

Ainsi :
2
1 G 1 G 2
v12C P1 v1D 2 2 P2 v2 (équation 4.72)
2 A 1 Un v2C 1
La substitution donne :

816 2 0,007422 /2C[1.4/ 1.4 1 ]12ð106ð0.00742

D 816 2 0,007462 /2C[1.4/ 1.4 1 ]P2ð106ð0.00746

et:P1D11,94 MN/m2
La chute de pression pour l'écoulement adiabatique D 12.0 11.94 D0.06 MN/m2ou 60 kN/m2

PROBLÈME 4.7
Air, à une pression de 10 MN/m2et une température de 290 K, s'écoule d'un réservoir
à travers un tuyau en acier doux de 10 mm de diamètre et de 30 m de long vers un deuxième réservoir à un
pressionP2 Tracez le taux massique de débit de l'air en fonction de la pression P2 Négligence
tous les effets attribuables aux différences de niveau et en supposant une expansion adiabatique de la
air. D0.018 mN s/m2 , D1.36.

Solution
G/A est requis en fonction de P2 . v2ne peut pas être trouvé directement puisque l'aval
températureT2est inconnu et varie en fonction du débit. Pour un flux adiabatique,
 ÉCOULEMENT DES FLUIDES COMPRIMABLES 67

v2peut être calculé à partir de l'équation 4.77 en utilisant des valeurs spécifiées de G/A et substitué dans
équation 4.72 pour obtenir la valeur de P2 De cette manière, les données requises peuvent être calculées.
2 2
1P1 A v1 C1 v2
8 R/ u2 l/d D C v 1 ln
2 1 G v2 v1
(équation 4.77)

0.5 G/A 2 v2C[ / 1 ]P1 v1D0.5 G/A 2 v2C[ / 1 ]P2 v2

1 2
(équation 4.72)

ou : 0.5 G/A 2 v21 v22 C[ / 1 ]P1 v1

[ / 1 ]v2 DP2

QuandP2DP1D10 MN/m2 ,G/AD0.

IfG/Ais 2000 kg/m2 s, alors :

ReD 0,01ð2000/0,018ð10 3 D1.11ð106

Quand/eD0.0002,R/ u2D0.0028 de la Fig. 3.7 et :

v1D 22.4/29 290/273 0.1013/10 D0.0083 m3 /kg

Substituer dans l'équation 4.77 :
2
0,36 10ð106 1
8 0.0028 30/0.01 D C 0,0083
2ð1.36 2000
2
0,0083 2,36 v2
ð1 ln
v2 1,36 0,0083

etv2D0,00942 m3 /kg.
Substituer pourv2dans l'équation 4.72 donne :

P2D[0.5 2000 2 0,008320,009422 C 1,36/0,36 10ð106ð0.0083] / 1.36/0.36

ð0.00942
et:P2D8,75 MN/m2 .
De la même manière, le tableau suivant peut être produit.

G/A kg/m2 s v2 m3 /kg P2 MN/m2

0 0,0083 10,0
2000 0,00942 8.75
3000 0,012 6,76
3500 0,0165 5,01
4000 0,025 3,37
4238 0,039 2.04
68 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

5000

4000

3000

Un

G
2000

1000

0
10 9 8 7 6 5 4 3 2 1 0
PressionP2(MN/m2)

Figure 4a.

Ces données sont tracées dans la Fig. 4a. Il est montré dans la Section 4.5.4, Volume 1, que le
la vitesse maximale qui peut se produire dans un tuyau dans des conditions d'écoulement adiabatique est le sonic
vitesse qui est égale à P2pv2 .
À partir du tableau ci-dessus pP2 v2le débit maximum est :

1.36ð2.04ð106ð0.039D329 m/s
La température dans cette condition est donnée par P2 v2DRT/M, et :

T2D 29ð0.039ð2.04ð106 /8314 D227 K

La vitesse du son dans l'air à 227 K est de 334 m/s, ce qui sert de vérification sur le
données calculées.

PROBLÈME 4.8
Sur une longueur de 30 m de tuyau à vide de 150 mm transportant de l'air à 293 K, la pression chute
à partir de 1 kN/m2à 0,1 kN/m2 . Si la rugosité relative / dis 0,002, quelle est l'approximative
débit?

Solution
Le volume spécifique de l'air à 293 K et 1 kN/m2est :

v1D 22,4/29 101,3/1,0 D83,98 m3 /kg

293/273
Il est nécessaire d'assumer un nombre de Reynolds pour déterminer R/ u2et ensuite calculer un
valeur de G/A qui devrait correspondre à la valeur supposée d'origine. Supposer un Reynolds
nombre de 1 à 105 .
 ÉCOULEMENT DES FLUIDES COMPRESSIBLES 69

Whene/dD0.002 etReD105 ,R/ u2D0.003 de la Fig. 3.7.

G/A 2ln P1 /P2 C P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0 (équation 4.55)
Substituant :

G/A 2ln(1.0/0.1) C(0.1)2 12 ð106 / 2ð1ð103ð83.98

C4 0,003 30/0,15 G/A 2D0

et: G/A D1,37 kg/m2 s.

La viscosité de l'air est de 0,018 mN s/m2 .

∴ ReD 0.15ð1.37 / 0.018ð10 3 D1.14ð104

Ainsi, la valeur choisie de Reis est trop élevée. Quand ReD1ð104 ,R/ u2D0.0041 et
G/AD1.26 kg/m2 s.
Renow équivaut à 1,04ð104ce qui correspond bien à la valeur supposée.

Ainsi : GD1.26ð /4 ð 0.15 2D0,022 kg/s

PROBLÈME 4.9
Un système sous vide est requis pour traiter 10 g/s de vapeur (masse moléculaire 56 kg/kmol)
afin de maintenir une pression de 1,5 kN/m2dans un vaisseau situé à 30 m du vide
pompe. Si la pompe est capable de maintenir une pression de 0,15 kN/m2à son point d'aspiration, quoi
Le diamètre du tuyau est-il nécessaire ? La température est de 290 K, et des conditions isothermes peuvent être.
soit supposé dans le tuyau, dont la surface peut être considérée comme lisse. La loi des gaz idéaux est
viscosité du gaz D0.01 mN s/m2 .

Solution
On utilise l'équation 4.55 pour résoudre ce problème. Il est nécessaire d'assumer une valeur de
le diamètre du tuyau afin de calculer les valeurs de G/A, le nombre de Reynolds et R/ u2 .
IfdD0.10 m,AD /4 0.10 2D0,00785 m2
∴ G/AD 10ð10 3 /0.00785 D1.274 kg/m2 s
et ReDd G/A / D0.10ð1.274/ 0.01ð10 3 D1.274ð104

Pour un tuyau lisse, R/ u2D0.0035, d'après la Fig. 3.7.

Volume spécifique à l'entrée,v1D 22.4/56 290/273 101.3/1.5 D28.7 m3 /kg

G/A 2ln P1 /P2 C P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0 (équation 4.55)
La substitution donne :

1.274 2ln(1.5/0.15)C(0.15)2 1.52 ð106 / 2ð1.5ð103ð28.7

C 0.0035 30/0.10 1.274 2D 16.3
et la valeur choisie de dis est trop grande.
70 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

Une valeur supposée supplémentaire de dD0,05 m donne une valeur du côté droit de
l'équation 4.55 de 25.9 et la procédure est répétée jusqu'à ce que cette valeur soit nulle.
Cela se produit lorsque dD0,08 m ou 80 mm.

PROBLÈME 4.10
Dans un système sous vide, l'air s'écoule isothermiquement à 290 K à travers un diamètre de 150 mm.
pipeline de 30 m de long. Si la rugosité relative de la paroi du tuyau est de 0,002 et le bas-
la pression de flux est de 130 N/m2 , quelle sera la pression amont si le débit d'air est
0,025 kg/s ? Supposons que la loi des gaz parfaits s'applique et que la viscosité de l'air soit constante
à 0,018 mN s/m2 .
Quelle erreur serait introduite si le changement d'énergie cinétique du gaz résultait
de l'expansion ont été négligées ?

Solution
En tant que volumes spécifiques en amont et moyensv1etvmsont nécessaires dans les équations 4.55 et
4.56 respectivement, on utilise l'équation 4.57 :

G/A 2ln P1 /P2 C P22 P12 / 2RT/M C4 R/ u2 l/d G/A 2D0

RD8.314 kJ/kmol K et donc :

2RT/MD 2ð8.314ð103ð290 /29D1.66ð105J/kg

Le deuxième terme a des unités de N/m.2 2 / J/kg Dkg2 /s2m4ce qui est cohérent avec le
autres termes.

AD /4 0.15 2D0.0176 m2
∴ G/AD 0,025/0,0176 D1,414
et ReDd G/A / D 0.15ð1.414 / 0.018ð10 3 D1.18ð104

Pour les tuyaux lisses et ReD1.18ð104 ,R/ u2D0.0040 de la Fig. 3.7. Remplaçant dans
l'équation 4.57 donne :

1.414 2ln₀P /130 C 1302 P12 /1.66ð105C4ð0.0040 30/0.15 1.414 2D0

1

Résoudre par essais et erreurs, la pression en amont, P1D1,36 kN/m2

Si le terme d'énergie cinétique est négligé, l'équation 4.57 devient :

P22 P12 / 2RT/M C4 R/ u2 l/d G/A 2D0 et P1D1.04 kN/m2

Ainsi, une erreur considérable serait introduite par cette hypothèse simplificatrice.
 ÉCOULEMENT DES FLUIDES COMPRESSIBLES 71

PROBLÈME 4.11
L'air s'écoule à un rythme de 30 kg/m2 à travers un tuyau lisse de 50 mm de diamètre et
Long de 300 m. Si la pression en amont est de 800 kN/m2 , quelle sera la pression en aval
soit le flux est isotherme à 273 K ? Prenez la viscosité de l'air comme 0,015 mN s/m2et
supposer que le volume occupe 22,4 m3 Quelle est la signification du changement de cinétique
énergie du fluide?

Solution

G/A 2ln P1 /P2 C P22 P12 /2P1 v1C4 R/ u2 l/d G/A 2D0 (éqation 4.55)
Le volume spécifique à la condition amont est :

v1D 22,4/29 273/273 101,3/800 D0,098 m3 /kg

G/AD30 kg/m2 s
ReD 0,05ð30 / 0,015ð10 3 D1.0ð105
∴

Pour un tuyau lisse, R/ u2D0.0032 de la Fig. 3.7.

La substitution donne :

30 2ln 800/P2 C P22 8002 ð106 / 2ð800ð103ð0.098

C4 0,0032 300/0,05 30 2D0

et la pression en aval, P2D793 kN/m2

Le terme d'énergie cinétique D G/A 2ln 800/793 D7,91 kg2 /m4 s2
C'est insignifiant par rapport à 69 120 kg2 /m4 s2quelle est la valeur de la
autres termes dans l'équation 4.55.

PROBLÈME 4.12
Si la température ne change pas avec l'altitude, estimez le point d'ébullition de l'eau à une hauteur.
à 3000 m au-dessus du niveau de la mer. La lecture du baromètre au niveau de la mer est de 98,4 kN/m 2et le
La température est de 288,7 K. La pression de vapeur de l'eau à 288,7 K est de 1,77 kN/m.2 . Le
le poids moléculaire effectif de l'air est de 29 kg/kmol.

Solution

La pression atmosphérique à 3000 m est P2et la pression au niveau de la mer, P1D98.4 kN/m2 .

vdPCgdzD0

vDv1 P/P1
72 INGÉNIERIE CHIMIQUE VOLUME 1 SOLUTIONS

dP
et : P1 v1 gdzD0
P C
et : P1 v1ln P2 /P1 Cg z2 z1 D0
v1D 22,4/29 288,7/273 101,3/98,4 D0,841 m3 /kg.
∴ 98,400ð0.841 ln P2 /98.4 C9.81 3000 0 D0
et : P2D68,95 kN/m2
La relation entre la pression de vapeur et la température peut être exprimée comme :
logPDaCbT
Quand, TD288,7, PD1,77 kN/m2
et quand, TD373, PD101.3 kN/m2
∴ logPD 5.773C0.0209T

QuandP2D68.95,TD364 K.

PROBLÈME 4.13
Une conduite de gaz de 150 mm est utilisée pour transférer du gaz (poids moléculaire 13 kg/kmol et kine-
viscosité matic 0,25 cm2 /s) à 295 K d'une usine à une station de stockage à 100 m de distance, à un
taux de 1 m3 /s. Calculez la chute de pression si le tuyau peut être considéré comme lisse.
Si la chute de pression maximale admissible est de 10 kN/m2 est-il possible d'augmenter le
débit de 25 % ?

Solution
Si le débit de 1 m3 /s est à l'état normal, le volume spécifique du gaz est :

22,4/13 D1,723 m3 /kg.

Le débit massique, GD 1.0/1.723 D0.58 kg/s.
Section transversale, AD /4 0.15 2D0.0176 m2
∴ G/AD32.82 kg/m2 s
/ D0,25 cm2 /sD0.25ð10 4m2 /s
et D 0.25ð10 4 1/1.723 D1.45ð10 5N s/m2
ReD 0.15ð32.82/1.45ð10 5 D3.4ð105
∴

Pour les tuyaux lisses, R/ u2D0.0017, de la Fig. 3.7.

La chute de pression due à la friction est :

4 R/ u2 l/d G/A 2D4 0.0017 100/0.15 32.82 2D4883 kg2 /m4s2

et: PD 4883/1.723 D2834 N/m2ou 2,83 kN/m2 .
 ÉCOULEMENT DE FLUIDES COMPRESSIBLES 73

Si le débit augmente de 25%, GD 1.25(0.58)D0.725 kg/s

G/AD41.19 kg/m2s
et : ReD0.15ð41.9/1.45ð105 D4.3ð105

et, à partir de la Fig. 3.7, R/ u2D0.00165

La chute de pression D4 0.00165 100/0.15 41.19 2 1.723D4.33 kN/m2(ce qui est moins
que 10 kN/m2 )
Il est donc possible d'augmenter le débit de 25 %.
 SECTION 5

Flux des mélanges multiphases

PROBLÈME 5.1
Il est nécessaire de transporter du sable d'une taille de particule de 1,25 mm et d'une densité de 2600 kg/m3à la
débit de 1 kg/s à travers un tuyau horizontal de 200 m de long. Estimez le débit d'air nécessaire.
le diamètre du tuyau et la chute de pression dans la canalisation.

Solution
Pour le transport pneumatique conventionnel dans les pipelines, un rapport masse de solides-gaz d'environ 5 est
employé. Débit massique d'air D 1/5 D0,20 kg/s
et, prenant la densité de l'air comme 1,0 kg/m3 , débit volumétrique de l'air D 1.0 (0.20)

D0.20 m3 /s
Afin d'éviter des chutes de pression excessives, une vitesse de l'air de 30 m/s semble raisonnable.
Ignorant le volume occupé par le sable (qui représente environ 0,2 % de celui occupé par le)
air), la section transversale de tuyau requise D 0.20/30 D0.0067 m2 ,
équivalent à un diamètre de tube pde 4ð0.0067/ D0.092 m ou 92 mm.
Ainsi, un diamètre de tuyau de 101,6 mm (100 mm) serait spécifié.
D'après le tableau 5.3 pour le sable de taille de particule 1,25 mm et de densité 2600 kg/m3 , le gratuit-
la vitesse de chute est :
u0D4,7 m/s
p
Dans l'équation 5.37, uG us D4.7/[0.468C7.25 4.7/2600 ]D6.05 m/s
La surface tranversale d'un tuyau de 101,6 mm de diamètre intérieur D ð0,10162 /4 D0,0081 m2 .

∴ vitesse de l'air, uGD 0,20/0,0081 D24,7 m/s

et : usD 24.7 6.05 D18.65 m/s
En prenant la viscosité et la densité de l'air comme 1,7 × 10 5N s/m2et 1,0 kg/m3respectivement,
le nombre de Reynolds pour l'écoulement de l'air seul est :

ReD 0,102ð24,7ð1,0 / 1,7ð10 5 D148,000

et à partir de la Fig. 3.7, le facteur de friction D0.002.

74