Flow in pressurized conduits 65
Computer code 4.1
% Example 4.1 Pipe Network - Hardy Cross Method
% q = Volumetric discharge [m^3/s];
% d = Pipe diameter [m];
% r = 8*f*l/(g*pi^2d^5) Resistance coefficient [s^2/m^5];
% ctr = Convergence criterion;
% nb = Number of branches;
Table 4.2 Input data of the connectivity matrix
Loop
number
Branch number
1
2
3
4
5
6
7
8
9
10
11
12 13
1
–1
1
0
0
1
0
0
1
0
–1
0
0
0
2
0
0
1
1
–1
–1
0
0
0
0
0
0
0
3
0
0
0
0
0
1
1
–1
–1
0
0
0
0
4
0
0
0
0
0
0
0
0
1
1
–1
–1
1
6
5
4
3
2
1
0 1
2
3
4
5
6
7
Number of iterations
Max discharge difference between successive iterations
Figure 4.5 Convergence of the Hardy Cross method.