Answer Happy

Analysis of Multi-Pipe Systems Using Hardy Cross Method Write a Mathlab program that can solve for the water flow rate distribution in a network of pipes using Hardy Cross Method capable of handling a maximum of ten loops with a maximum of ten lines in each loop. The program has to use Hazen Williams Equation for major pipe losses and symbols given on page 2. Table 1 on page 3 outlines the input data to be read into your program for the test case. Table 2 outlines the final solution for the flow rate distribution as well as the head distribution with assumed total head of 120 ft. at node B. Delta Z across the entire system is assumed to be zero. Verify that your program works with this test case with C=100. Indicate how many iterations it took to arrive at the solution. What is the maximum percent error of your flow rates as compared to the solution given (obtained by Darcy Equation).
Computer Assignment: Write a computer program that will determine (a) the flow-rate distribution in a loop network system by the Hardy-Cross method and (b) the total head at each node. Assume that starting values for Qij will be read into the program. Allow for a maximum of ten loops with a maximum of ten lines in a loop. Assume that the lines in the network might be made from different materials. Write your program so that part (a) is first completed before it proceeds to part (b). Use the following symbols in your program: IW NL NJ (I) ZL (I,J) D(I,J) Q (I,J) EPSLN = No. of Iterations = Number of loops. = Number of lines in Ith loop DELQ (I) ID (I,J) = Length of the Jth line in Ith loop Diameter of the Jth line in Ith loop Volume-flow rate in Jth line in Ith loop. Maximum allowable difference in flow rates between iterations (use 0.0001) = = HF (I,J) = head loss in Jth line in Ith loop DHDQ (I,J) = Rate of change of the head loss with respect to flow rate for the Jth line of the Ith loop PH(I,K) Pressure head in the Kth node in Ith loop (K is the number of nodes in each of the loops) = Flow rate correction applied to each line in Ith loop. = Index for Jth line in Ith loop (Its value equals zero if line (I,J) is not a common line, and equals k if line (I,J) is in common with a line in loop k) = head loss coefficient for Jth line in Ith loop Minor head loss coefficient for Jth line in Ith loop ALPHA (I,J) = BETA (I,J) = F (I,J) Darcy-Friction factor for Jth line in Ith loop The determination of the total head at each node requires that the calculations proceed in consecutive order from node to node and from loop to loop.
4.0 10.8 cfs 3.8 2.4 cfs 1.6 Loop 1 3.0 Loop 2 2.1 cfs A typical three-loop network. 2.1 cfs 1.5 0.3 1.0 Loop 3 1.8 1.4 cfs 1.4 cfs 1.4 cfs Cast Iron Pipe € = 0.00085 ft ✓ = 1.082 X 10 $ = 8 22=10 S2+1=3 select: cfs ава QBC=+3.8 cfs Фен - 4.0 = -1.5 cfs
TABLE 1 Loop Line No. No. (ft) 1 2 17 3 123 4 1234 1234 4 Length Diameter (0) (in.) (cfs) Badd 10560 16 15840 12 10560 14 15840 16 15840 16 13200 14 10560 12 10560 16. 6426 15840 15840 15840 12 15840 12 12 10 NNON 3.0 1.5 -1.6 -4.0 3.8 1.7 0.3 -3.0 1.8 0.4 -1.0 -1.5 TABLE 2 Loop Node Line. Swamee-Jain No. No. 1 2 32 1234 BEHAB BUANAB В с E E ܪ 3 ܝܐ ܕܝܐ E 1234 1224 3 Q (cfs) 3.780 1.508 -1.640 -4.040 2.980 0.880 -0.520 -3.780 1.753 0.353 -1.047 -1.508 H (ft) 120.000 102.902 84.085 90.816 120.000 120.000 103.823 101.244 102.902 120.000 102.902 77.730 74.734 84.085 102.902