Flux dans les tuyaux : laminaire et turbulent
Flux dans les tuyaux : laminaire et turbulent
226
i. Connaître la résistance au flux dans les tuyaux (pertes
par : frottement, élargissement, contractions,
vannes, etc.).
ii. Currents d'écoulement dans des tuyaux laminaires et turbulents.
iii. Effets du fluide à l'intérieur des tuyaux « coup de
Ariete
iv. Analyse de tuyaux simples et multiples.
v. Analyse des tuyaux dans un système de réseaux et méthodes
de résolution,"Méthode de Cross, Méthode de Newton-
Raphson
227
Il peut arriver que les particules aient des mouvements indéfinis, dans ce
caso le liquide ne circule plus en formant des filets parallèles, mais plutôt des trajectoires plus ou
Imaginons un tube par lequel va circuler de l'eau avec des vitesses qui
nous pouvons faire varier. Au début, nous allons faire en sorte que le liquide contenu dans le
tuyau, et qui est au repos, circule en passant par une section donnée, avec une
vitesse progressivement plus élevée. quand le liquide commence à bouger, comme
la vitesse est faible, le régime est laminaire ; en augmentant progressivement
228
la vitesse atteint une certaine valeur à partir de laquelle le régime devient turbulent.
Cette valeur que doit atteindre la vitesse pour passer du régime laminaire au
turbulent, ça s'appelle Vitesse Critique Élevée.
Supposons maintenant que l'eau circule à grande vitesse dans le tuyau,
de telle manière que le régime est turbulent et qui le fait diminuer
Paulatinement ; le régime reste turbulent jusqu'à ce qu'un moment arrive
que si la vitesse diminue encore, elle passe à laminaire. À cette vitesse pour la
La vitesse à laquelle le flux passe de turbulent à laminaire est appelée Vitesse Critique Basse.
Les lois de l'écoulement laminaire dans les tuyaux droits ont été déterminées
expérimentalement, de manière indépendante par Hagen et Poiseuille. La loi de
Hagen-Poiseville peut être déduit des principes fondamentaux,
considérez les forces qui agissent sur un cylindre de
Liquide de longueur "l" et de rayon "y", comme illustré dans la figure 5-3. Si le
le mouvement est uniforme, la force engendrée par la différence de pression sur
les extrémités du cylindre doivent s'équilibrer exactement avec la résultante de
les efforts de cisaillement en sa surface latérale, comme l'exprime l'équation
(5.2.1):
(p1 p 2) y 2 2 yl
5.2.1
En simplifiant et en introduisant la valeur de τ :
du
(p1 p 2) y 2
l
dy
5.2.2
Où :
u; c'est l'adresse selon la figure 5-3
μ; viscosidad (Kg-s/m2)
229
Le deuxième membre est négatif à cause du fait que du et dy sont de signes opposés.
(Fig5-3), en isolant l'équation 5.2.2 nous aurons :
(p p ) y 2
u du [Kg-s/m2]
5.2.3
(p1 p 2)d2
C
16µl
5.2.4
( p 1 p 2) 2 d2
u y
4ul 4
5.2.5
5.2.7
230
En substituant dans l'équation 5.2.5 :
( p 1 p 2) d/ 2 d 2 ydy
Q y 2 dy
2 l 0
4
(p1 p2)d 4
Q
[m3/s]128 l
5.2.8
(p p 2)d4
[m/s] 32 l
5.2.9
On peut voir par les équations 5.2.7 et 5.2.8 que dans un écoulement laminaire, la
la vitesse moyenne est la moitié de la vitesse maximale. L'expression pour la
la variation de la charge piézométrique sur une longueur « l » se l'obtiendra de la
de la manière suivante :
( p 1 p 2) 32 lv
h
w d2 g
5.2.10
231
64l v2
hf
Red2g
5.2.11
l v2
Ec. Darcy-
Weisbach
5.2.12
Où :
f; c'est le facteur de frottement
v; c'est la vitesse
d; diamètre de la tuyauterie
Re; nombre de Reynolds
hfperte par friction
Connue sous le nom de formule de Darcy-Weisbach pour l'écoulement dans les tuyaux,
étant h la perte d'énergie en l mètres de tuyau en kilogrammes par
kilogramme. Ce qui précède démontre que dans un courant laminaire, « f » est complètement
indépendant de la rugosité, car il varie seulement avec le nombre de Reynolds, est
dire, avec la résistance relative des forces visqueuse et d'inertie.
232
Dans les tuyaux très lisses, les valeurs de "f" varient avec R (ANNEXE).
On peut observer qu'elle ne devient jamais une droite horizontale, ce qui
démontre que les propriétés des liquides influencent le courant dans tout
l'intervalle des nombres de Reynolds. Les tubes en verre et ceux en métal étiré
avec une surface très lisse se trouvent dans cette catégorie.
La figure 5-4 illustre le débit dans des tuyaux rugueux à des valeurs élevées de
R par le biais de la zone située au-dessus et à droite de la ligne de traits, appelée
zone de turbulence complètement établie ou pleine. Dans cette zone, les courbes
«f» deviennent horizontales, prouvant ainsi que le courant est complètement
indépendant des propriétés des liquides. La troisième catégorie de
la turbulence apparaît lorsque les valeurs de "f" se trouvent dans le
zone comprise entre la courbe pour tuyaux lisses et la ligne de pointillés (Fig5-4).
Le courant dans les tubes commerciaux se produit généralement à l'intérieur de celui-ci.
categoría. En esta zona, las curvas de “f” para los diversos valores de la
la rugosité relative se sépare en points successifs de celles pour les tuyaux lisses
et elles deviennent horizontales lorsqu'elles entrent dans la zone de turbulence pleine.
Toutes ces courbes présentées dans la figure 5-4 ont été présentées par
Mood. C'est pourquoi la figure 5-4 est connue sous le nom de Diagramme de Moody, avec le
quel connu le nombre de Reynolds et la rugosité équivalente nous pouvons
trouver le facteur de frottement correspondant. Cependant, si cela n'est pas connu
le nombre de Reynolds, ayant pour inconnues le débit ou le diamètre, le facteur de
Le frottement s'obtient par des approximations successives.
233
Source : Mécanique des fluides appliquée
Robert L. Mott
Coefficients de rugosité pour certains matériaux
234
Source : Mécanique des Fluides Appliquée DIAGRAMME DE
Robert L. Mott
HUMEUR
(30)(0.0833)
NR 5,70 105
4,38 dix6
Par conséquent, le flux est turbulent. Maintenant, il faut calculer le
rugosité relative. Du tableau 5.1 nous avons Є = 8x10-4pieds. Alors,
la rugosité relative sera :
235
D0.0833
1,04 10 2 104
8 10 4
Observez que, afin que D/Є soit un quotient sans
Les dimensions, tant D que Є doivent être dans les mêmes unités.
Enfin, il faut localiser le nombre de Reynolds sur l'axe des abscisses.
du Diagramme de Moody (Fig5-4)
NR 5,70 105
nous projetons verticalement jusqu'à ce que la courbe soit atteinte
correspondant à D/Є=104, nous pouvons déjà utiliser la courbe numéro 100
que est proche. Nous projetons horizontalement vers la gauche, on
lisez la valeur qui dans ce cas est d'environ f=0.038.
f 64/NR
5.4.1
Cette relation est représentée dans le diagramme de Moody comme une ligne droite.
du côté gauche du diagramme. Bien sûr, pour les nombres de Reynolds
de 2000 à 4000, le flux se trouve dans la région critique et il est impossible
prédire la valeur de “f”.
Por encima del número de Reynolds de 4000, por lo general el flujo se conoce
comme turbulent. Cependant, en essence, il existe deux zones d'intérêt dans ce
236
point. Vers le côté droit du diagramme, le flux est dans la zone complète
turbulence. On peut observer que la valeur de « f » ne dépend pas du nombre de
Reynolds, mais seulement de rugosité relative. Dans cet intervalle, on applique la
ecuación:
1
2log10 3.7D/
f
5.4.2
La frontière de cette zone est la ligne pointillée qui court, en général, de la
partie supérieure gauche à la partie inférieure droite du diagramme, l'équation de
cette ligne est :
1 NR
f 200(D/ )
5.4.3
La troisième zone du diagramme de Moody, qui est connue sous le nom de zone de
transition, se trouve entre la zone de turbulences complètes et la ligne qui se
s'identifie comme des conduits lisses. La ligne des "conduits lisses" a l'équation :
1 NR f
2log10
f 2,51
5.4.4
5.4.5
237
L'équation suivante permet le calcul direct de la valeur du facteur de
friction, cela a été développé par P.K. Swamee et A.K. Jain,
0,25
f 2
1 5,74
journal
3.7(D/ ) N R0,9
5.4.6
L'équation 5.4.6 produit des valeurs pour “f” qui se situent entre 1.0%
Solution. Étant donné qu'il se trouve dans la zone turbulente, nous utilisons la
équation 5.4.6
0,25
f 2
1 5,74
journal
3.7(2000) (1 105) 0,9
f 0,0204
238
L'équation de Darcy pour la perte par friction peut être utilisée pour
sections transversales non circulaires si la géométrie est représentée par le
radio hydraulique au lieu du diamètre du conduit, comme utilisé dans les sections
circulaires. Après avoir calculé le rayon hydraulique, nous pouvons calculer le nombre
de Reynolds en remplaçant 4R par le diamètre, D. L'équation, alors,
reste comme :
(4R) (4R)
NR
v
5.5.1
Une substitution similaire peut être faite pour calculer la perte par friction.
Dans l'équation de Darcy :
l v2
hf f
4R2g
5.5.2
239
NR 2,96x104
240
Cette formule, établie tant pour les tuyaux que pour les canaux ouverts, a
a été beaucoup utilisée pour les tuyaux. La sélection des exposants a été faite avec l'idée
de parvenir à une variation minimale du coefficient pour tous les conduits de
même degré de rugosité.
5.6.1
Où :
Rh; radio hydraulique (m)
S; en attente de la ligne d'énergie
CHWcoefficient de Hazen-Williams, sans dimension
0,63 0,54
D hf
v 1.538C TH
4 l
5.6.2
Enlever
v1.85l
hf 2,271 1.85
C HWD1,16
5.6.3
Si nous multiplions et divisons le second membre par l/D, nous aurons :
241
v1,85 l
hf 2,271 1,85
C HWD1,16 D
5.6.4
2
De la même manière, nous multiplions et divisons le second membre par v 2g ,
nous avons :
44.56 l v2 44,56 l v2
hf 1,85
C HWD1.16v 0,15 D2g 1,85
C HWD 0,01 Dv 0,15 D2g
5.6.5
44,56 l v2
hf 1.85
C HWD Re 0,15v 0,15
0,01 D2g
5.6.6
44,56
f 1,85
C HWD Re 0,15v 0,15
0,01
5.6.7
242
Source: Introduction à l'Hydraulique
Texte auxiliaire de la matière
C'est une formule empirique, obtenue par expérimentation, et c'est une autre
alternative pour résoudre des problèmes d'écoulement dans des conduits fermés et des canaux
ouverts, cela est considéré pour des tuyaux d'un mètre de diamètre, étant très
fiable pour la gamme de diamètres compris entre 0,40 et 1,30 m.
1
v Rh 23/S 1/2
n
5.6.8
Où :
Rh; radio hydraulique
S; en attente de la ligne d'énergie
n ; coefficient de Manning, sans dimension
5.6.9
Despejando :
243
n 2 lv2
hf 6.34
D4/3
5.6.10
En multipliant et en divisant le membre de droite par D et par 2g, nous aurons :
n 2l v2
hf 124,4
D1/3 D2g
5.6.11
en égalant cette équation à l'équation de Darcy-Weisbach, nous aurons :
n2
f 124,4
D1/3
5.6.12
Le coefficient « n » dépend du matériau dont est fabriquée le tube, le la
Le tableau 5.3 présente quelques matériaux avec leurs valeurs respectives de « n ».
244
De la figure, il découle que le tronçon B-A alimente d'un seul côté, tandis que
que le tronçon D-C qui approvisionne des deux côtés devra entraîner plus de dépenses. Si se
tratara de una zona con la población uniformemente distribuida el tramo D-C
conduira le double de dépenses que le tronçon B-A. Lié à des dépenses avec
longitudes, il semblerait que le tronçon D-C mesure le double de la longueur que le tronçon
B-A, mais en réalité cela n'est pas vrai, par conséquent, on considère que sa longueur
virtuel est de 200 et sa longueur réelle de 100. Le segment B-A a une longueur
virtuel égal à sa longueur réelle de 100.
Pour les tuyaux d'alimentation LVirtuel=0
Pour les immeubles qui alimentent un seul côté LVirtuel= LRéel
Pour les propriétés qui alimentent les deux côtés LVirtuel= 2*LRéel
En additionnant ces longueurs virtuelles (de tout le système de réseau), on pourra
trouver, on obtient le coefficient de débit par mètre de tuyau :
Qmax h
q
LVirtuel [l/s*m]
Où :
q ; Coefficient de débit par mètre de tuyau
Qmax-h ; Caudal máximo horario (l/s)
LVirtuel Somme des longueurs virtuelles (m)
245
Cette méthode consiste à ajouter à la longueur réelle du tuyau une longueur
équivalent aux pertes de charge qui se produisent dans les pièces spéciales
existants dans le tuyau. Chaque pièce spéciale correspond à une longueur fictive et
adicional (1). Le tableau suivant montre la longueur équivalente pour différents
pièces.
246
Solution.
m
L V2V2
f Kx Z Un
Z factorisation
D2g2g x 1
B
m
V2 L
(f Kx) Z AZ B
2g D x 1
m
V2 8Q2 8Q2 L
comment (f Kx) Z AZ B
2g 2
D 4g D g
2 4
D x 1
Datos: g Z Un Z B
Q
8 L 1 m
ZA - ZB = 80 m f 5 ( Kx )
D D4 x 1
L = 1150 m
D N = 250 mm,
D I = 231,7 mm
g = 9,81 m/s²
π= 3.1416
Kcodo 22.5º=0.2 (voir annexes - pertes locales),
Kcodo 45º =0.4 ,
Kcodo 90º = 0,9
247
KVálve de comp. = 0,2
ΣKx = 3(0,2)+6(0,4)+2(0,9)+2(0,2)=5,2, ΣKx/D4 = 5,2/(0,2317)4 =1,804,26
L/D5 = 1150/(0.2317)5 = 1.7221 x106 ,
La formule pour calculer Q est :
De Reynolds
Re = (1.2606 x 106 ) (Q/0.2317 m) = 5.44 x 106 Q , avec Q en m3 /s
ε/D = 0,0015 mm / 231,7 mm = 6,474 x 10^6
Comme on peut l'observer, le nombre de Reynolds dépend du débit,
par conséquent, la f dépend également de la dépense, pour résoudre ce type de
des problèmes se font par tâtonnements (par essais et erreurs). Régulièrement, on en a entre
cinq et six itérations.
(m³/s)
1 0,100 5,4E+05 0,013006 0,200
2 0,200 1.1E+06 0.011614 0,211
3 0.211 1,1E+06 0,011520 0,212
4 0.212 1,2 million 0,011514 0,212
5 0,212 1,2E+06 0,011513 0,212
6 0,212 1,2 million 0,011513 0,212
Solution (a)
Q = 212 lps.
V= 5.02 m/s
V.9 PÉRDIDAS MENORES
248
Dans la plupart des systèmes de flux, la perte d'énergie primaire se
doit à la friction du conduit, comme cela a été vu dans le sous-titre précédent. Les autres
les types de pertes sont généralement petites en comparaison, et par conséquent
on les appelle des pertes mineures.
Les pertes mineures se produisent lorsqu'il y a un changement dans la section transversale de
la trajectoire d'écoulement ou dans la direction d'écoulement, ou quand la trajectoire d'écoulement
se encuentra obstruida, como sucede con una válvula. La energía se pierde bajo
ces conditions en raison de phénomènes physiques assez complexes. La prédiction
la théorie de l'ampleur de ces pertes est également complexe, et donc,
Normalement, on utilise les données expérimentales.
Dans certains cas, il peut y avoir plus d'une vitesse d'écoulement, comme avec les
dilatations ou dans les contractions. Il est de la plus haute importance que l'étudiant
sachez quelle vitesse utiliser avec chaque coefficient de résistance. Si la tête
de vitesse v 2/2g de l'équation 5.7.1 s'exprime dans les mêmes unités de
mètres, donc la perte d'énergie sera également en mètres ou N*m/N de
flux de fluide. Le coefficient de résistance n'a pas d'unités, car il représente
une constante de proportionnalité entre la perte d'énergie et la tête de
vitesse. La magnitude du coefficient de résistance dépend de la géométrie du
249
dispositif qui provoque la perte et parfois dépend de la vitesse de
flux.
En faisant passer un fluide d'un conduit de diamètre inférieur à un conduit de diamètre supérieur à travers une
hL K ( v 1 2/2g)
5.9.2
1 A/A 2 2
K 1 2 1 (D 1/D2) 2
5.9.3
250
Source : Mécanique des fluides appliquée
Robert L. Mott
251
EXEMPLE D'APPLICATION Déterminez la perte d'énergie que
cela se produira lorsque 10L/min d'eau s'écoulera à travers une dilatation soudain
provenant d'un tube en cuivre de 1 pouce à un tube de 3 pouces (type
K)
v1 2 3.322
m 0,56 m
2g2(9.81)
Pour trouver une valeur de k, la proportion du diamètre est requise,
alors nous aurons :
D2/D1 73,8/25,3 2,92 mm
Pendant l'écoulement d'un fluide d'un conduit vers un grand réservoir ou cuve
(Fig5-8), sa vitesse diminue presque à zéro. Dans le processus, l'énergie
cinétique que le fluide possédait dans le conduit, indiquée par la tête de vitesse,
se dissipe. Par conséquent, la perte d'énergie pour cette condition est :
252
hL 1.0(v1 2/2g)
5.9.4
Ceci est appelé la perte de sortie. La valeur de k=1,0 est utilisée indépendamment de la
forme de la sortie où le conduit se connecte au mur du réservoir.
Si la transition d'un conduit plus petit à un plus grand peut se faire moins
abrupta que la dilatación súbita de bordes carrés, la perte d'énergie se
réduire. Cela se fait généralement en plaçant une section entre les deux conduits
(Fig5-9). Les parois inclinées du cône tendent à guider le fluide pendant la
désaccélération et expansion du flux de courant.
253
hL k(v 12/2g)
5.9.5
254
Source: Mécanique
de Fluidos Appliquée - Robert
L. Mott
255
Source : Mécanique des fluides appliquée
Robert L. Mott
v12/2g 0,56 m
D2/D1 2,92
256
La perte d'énergie due à une contraction soudaine (Fig5-11) est calculée à
partir de :
hL k(v 22/2g)
5.7.6
257
Source : Mécanique des Fluides Appliquée
Robert L. Mott
258
EXEMPLE D'APPLICATION Déterminez la perte d'énergie que
il se produira un écoulement de 100l/min d'eau d'un tube en cuivre de 3 pouces (type
k) dans un tube en cuivre de 1 pouce (type K) à travers une contraction
soudain
v22/2g 0,56 m
De la figure 5-12 nous trouvons k=0.42, alors :
hL k(v 22/2g) (0,42)(0,56) 0,24 m
259
La perte d'énergie lors d'une contraction peut être considérablement réduite
faisant la contraction plus progressive. La figure 5-14 montre une contraction de
ce type, formé par une section conique entre les deux diamètres avec
changements brusques dans les joints. L'angle "θ" est appelé l'angle du cône.
260
cono mineurs, la transition entre les deux diamètres est très longue, ce qui
augmente les pertes de friction.
261
Source : Mécanique des fluides appliquée
Robert L. Mott
262
Dr> 0,15, aucune veine contractée ne se forme, la perte d'énergie est assez
petite, et nous utilisons K = 0.04.
263
Solution. Nous avons les données des exemples d'application précédents.
par conséquent
Pour :
a)
D2 25,3 mm
Un2 5,017x104 m 2
Q
v2 3,32 m/s
Un2
v22/2g 0,56 m
Pour une entrée de projection vers l'intérieur, K=1.0, alors :
hL k(v 22/2g) (1.0)(0.56) 0,56 m
b)
Pour une entrée bien arrondie, K=0.04, alors :
hL k(v 22/2g) (0,04)(0,56) 0,02 m
264
Source : Introduction à l'hydraulique - Texte
assistant de la matière
265
V.9.8.3 VANNES À PAPILLON
Selon les valeurs de la figure 5-20
V.9.9.1 BIFURCATIONS
Les bifurcations généralement désignées par leur forme sont souvent des "t" ou
oui, vous causez des pertes qui dépendent de l'angle formé par la tuyauterie
secondaire avec la principale ou la maîtresse, la relation entre les diamètres des deux
tuyaux et de la direction du courant. cette perte est plus grande là où elle se produit
union des deux tuyaux, c'est-à-dire là où se trouve la bifurcation.
La figure 5-21(a, b, c, d, e, f, g, h, i) montre plusieurs situations dans les
que se pourraient présenter ce type de bifurcations.
266
267
268
269
270
271
Source de graphiques 5-21a-5-21i : Introduction
à l'Hydraulique - Texte auxiliaire de la matière
Il est souvent plus pratique de courber un conduit ou un tuyau que d'installer un coude
commercialement fait. La résistance au flux d'un coude dépend de la
proportion du rayon "r" du coude avec le conduit à l'intérieur du diamètre D. Dans la
La figure 5-22 montre que la résistance minimale se produit lorsque le ratio r/D
c'est environ "3". La résistance se donne en termes de la proportion de
longueur équivalente Le/D, et par conséquent, l'équation :
k (Le/D)f T
5.9.7
Où :
Le; longueur équivalente et c'est la longueur du conduit droit du même diamètre
nominal comme l'accessoire qui aurait la même résistance que celui-ci.
fTfacteur de friction dans le conduit auquel l'accessoire est connecté.
272
Source : Mécanique des Fluides Appliquée
Robert L. Mott
Cette équation 5.7.7 doit être utilisée pour calculer le coefficient de résistance. La
La résistance montrée dans la figure 5-22 inclut à la fois la résistance du coude ainsi que la
résistance due à la longueur du conduit dans le coude. Lorsque nous calculons le
proportion r/D, "r" est défini comme le rayon à la ligne du centre du conduit ou
tuyau, appelé le rayon moyen (Fig5-23).
r Ri Faire/ 2
5.9.8
r Ro Faites/ 2
5.9.9
273
r (Ro Ri)/ 2
5.9.10
274
r Ri Faire/ 2
275
dv
p m
dt
5.10.1
5.10.2
Où :
hjeSurpression par coup de bélier
v ; vitesse de l'eau dans le tuyau
Eunmodule d'élasticité de l'eau
Et; modulo de elasticidad del material de la tubería
d ; Diamètre intérieur du tuyau
e ; Épaisseur du tuyau
E
MATERIAL
[Kg/cm2]
Acier 2,10*1006
Fonte 9,3*1005 Modules d'élasticité
Béton simple 1,25*100 para algunos materiales
Amiante
ciment 3,28*1005
PVC 3,14*1004
Polyéthylène 5,2*1003
Eau 2 067 * 1004
276
Source : Approvisionnement en eau
potable-Enrique Cesar Valdez
5.11.1
277
Où :
Qa; débits qui alimentent le réseau.
C : Consommations.
5.11.2
où pour un nœud donné, l'indice "j" indique le nœud vers lequel il va
le flux et l'indice "i" d'où il provient ; le terme C représente les
débits retirés.
à chaque nœud, il n'existe qu'un seul niveau d'énergie, par conséquent, les pertes de
un nœud à un autre, indépendamment du chemin de liaison sélectionné, sont
égaux. Cela conduit à la conclusion suivante : « si l'on considère la perte d'énergie
positive à celle qui se produit dans le sens contraire des aiguilles d'une montre et
viceversa, se podrá establecer para cada malla de red:
hij 0 K jej Qjenjij
5.11.3
V.11.1 SYSTÈME Q
278
Consiste à supposer des dépenses initiales quelconques dans les tuyaux de la
réseau qui respecte les équations de nœud. Chaque dépense initiale est affectée par
un accroissement correctif Q ce qui est considéré comme constant pour les tuyaux de
chaque maille. Les augmentations sont calculées par les équations de maille.
Le système d'équations comprendra alors "m" équations non linéaires
avec un nombre égal d'inconnues que sont les Q de chaque maille.
K ij(Qoij QJe) nij Je
0
……..
……..
K ij(Qoij QJe) nij m
0
Où :
Qo; dépense initiale supposée.
I, II,.. m; quantité de mailles
V.11.2 SYSTÈME H
Les inconnues dans ce système sont les niveaux ou hauteurs d'énergie "H" dans
chaque nœud. Les dépenses dans les tuyaux s'expriment en fonction de ces niveaux et à
la suite établit les équations de nœud correspondantes. Comme seul
il existe "n-1" équation indépendante de nœud, il est nécessaire de connaître le niveau
d'énergie dans un nœud quelconque. Il faut prendre en compte que :
1
n H i Hj nij
hij H je Hj K jej Q jej ij ; Qij
K ij
Le système H est :
n1ij
H je H j
C 0
K ij
1
1
nij
H je Hj
C 0
K ij
2
279
……
……
1
n
H je H j
ij
C 0
K ij
n1
Où :
1, 2,…., n-1; représentent les n-1 nœuds avec des équations indépendantes, les
équations de ce système est non linéaire.
V.11.3 SYSTÈME Q
Qij C 2
0
…..
Qij C n1
0
Équations de maillage :
K ijQjenjij Je 0
K ijQjenjijII 0
………
K ijQjenjijm 0
a) Méthode de Cross
280
Cette méthode consiste à rendre linéaires les équations du système Q
b) Méthode de Newton-Raphson
281
En chaque nœud du réseau, il y aura une fonction Fj, par conséquent, si les dépenses
Les alimentateurs de la tuyauterie sont connus, il y aura "n-1" fonctions et quand ce n'est pas le cas.
lo sean c'est-à-dire des mailles imaginaires le numéro sera « n ».
On commence à travailler avec des valeurs de H quelconques supposées. Sans
embargo, les suppositions doivent être raisonnables, sinon, la méthode
pourraient ne pas converger. Ces valeurs initiales sont corrigées en leur soustrayant un
increment fini H qui peut être calculé à l'aide de l'équation suivante,
exprimée sous forme matricielle :
F1 F 1 F1
....
H1 H2 Hj
H1 F1'
F2 F2 F2
.... H2 F '2
H1 H2 Hj
. .
. . . . * =
. .
. . . . . .
. . . . Hj
Fj Fj Fj F' j
....
H1 H 2 Hj
Où F1’, F2’, etc., ce sont les valeurs des fonctions Fjpour les valeurs
initiales de H. Pour les effets de placer des signes, il est opportun de rappeler que la
fonction Fjreprésente l'équation de continuité à chacun des nœuds,
cela implique que la dépense d'effluent est positive et que la dépense d'affluent est négative,
de plus, le sous-indice "i" correspond au nœud d'où provient le flux, et
«j» vers le nœud où il se dirige.
282