Application of EPANET 2.2 Software and Jal-Tantra Web System for Optimal Hydraulic Design of Water Distribution System for University of Kashmir †

: The availability of portable water in sufﬁcient quantity and standard quality is one of the basic requirements of a civilization. The study area of this present work is the University of Kashmir, where the existing WDS has served its function for more than 60 years. The UOK is a growing campus


Introduction
A well-designed, efficiently-performing WDS is a prerequisite for both urban and rural civilizations. Various software packages are available for modelling a WDS. The public domain software that is available free of charge include EPANET, Pipes and Loops, Branch, etc., and the paid ones available includes WaterCAD, Watergems, Synergy Water, etc. [1]. The paid versions are found to be more efficient in reference to functionality and compatibility with different computational systems, a graphical user interface, and searching and optimizing algorithms [2]. However, EPANET is freely available software that offers a fair degree of accuracy and ease of use when cost is a constraint [3]. EPANET 2.0 was used by various researchers to design the WDS for different areas [4,5]. The design includes the layout of the WDS, modelling the pressure head at nodes, and the velocity of flow in the links. However, no consideration was given to maximisation of the efficiency of performance or minimisation of the capital cost of the network.
The design of a WDS by EPANET must be such that it serves as an "optimal" solution, which can be obtained by setting out objective functions of maximum efficiency of performance and minimum capital cost of the network [6,7]. The decision variables and the constraints for the optimality of the network design have to be decided after conducting a study of the factors affecting the optimality of the design and consulting various codes and manuals laying down the standard range of values for those parameters [8,9]. The decision variables for optimal hydraulic design include the pressure head at nodes, the velocity of flow in the links, and the staging height of the storage tank. The check for cost optimality of the network formulated by EPANET can be conducted using the Jal-Tantra web system [10,11]. The optimality of the network can be further checked by carrying out a performance analysis using various indices [12][13][14].
The objective of this current study is to provide an optimally designed WDS for the University of Kashmir using EPANET 2.0 software and the Jal-Tantra web system. The existing distribution system has already served for more than 60 years and is unable to meet the requirements of the growing population of the campus. Moreover, the worn-out links and other components of the network lead to increased pipe and joint losses, reducing the efficiency of the existing distribution system. Therefore, there is an immediate need to redesign the water distribution network such that the design is an optimal solution meeting the objective functions of maximum efficiency of performance and minimum capital cost of the network, rather than just a simulation of the pressure head and velocity of flow.

Study Area Description
The study area for our current work is at the University of Kashmir, located in Hazratbal, Srinagar, Jammu and Kashmir, India. It is a residential campus comprising academic blocks, administrative blocks, hostels, a medical unit, a cafeteria, and staff quarters. It has a residential and floating population of 5324 and 1050, respectively, for the Naseem Bagh Campus and 8066 and 1720, respectively, for the Hazratbal Campus.

Methodological Approach
A quantitative research approach was employed to arrive at the optimal design of the WDS for the University of Kashmir by using various software. A literature review related to the work was carried out, and the related standard books and codes were consulted. A preliminary survey of the study area was followed by a division of the study area into two parts, with separate WDNs for each of them (the Hazratbal campus and the Naseem Bagh campus). The water source for the networks is the groundwater reserve within the campus, which is pumped by intermittently operating pumps into the overhead water tanks.

Methods of Data Collection
The quantitative data of the existing population at various academic blocks, hostels, family quarters, and administrative blocks was obtained. The maximum depth of water in the bore wells below ground level is 10 m. In addition, the campus plan of the study area was obtained. The data availability was facilitated by the engineering wing of the UOK. Elevation data was obtained using Google Earth Pro.

Evaluation of the EPANET 2.2 Input Parameters
A dead-end type of distribution system was provided for both parts of the campus. The population forecast for the design period of 30 years, with a constant increase of 5% per decade in both residential as well as the floating population was carried out [15,16]. The values of ADD for both parts of the network were evaluated using water demand as 135 lpcd for the residential population and 45 lpcd for the floating population [17]. The MDD was evaluated using a peak factor of '1.8' [18], and a peak factor of '3' was used to calculate the MHD from the AHD [18]. Provisions to meet the fire demand were kept as per Indian standards [19]. The base demand multiplier equals the coincident draft in litres per person per minute (Table 1). Two different demand patterns were used for the network nodes: one for the hostels and quarters (p1), and the other for academic buildings and administrative blocks (p2). The capacity of the storage tank was evaluated by the mass curve analysis for both parts of the network. The discharge to be provided by the pump was worked out by considering 12 h of pumping per day. The staging height of the tank above ground level was selected after iterations to successfully run the software and achieve standard pressure heads at the nodes. The economical diameter of the rising main was obtained using Dupit's expression. The pressure head developed by the pump has to be equal to the sum of the delivery head, suction head, and frictional head loss in the rising main. The input power required to run the pump equals (ωQH/n), and the efficiency of the pump (n) was taken as 75%. Energy required for a 12-hour per-day operation of the pump was evaluated to arrive at the cost of the electrical energy at Rs 5 per unit per day. The design flow rate of the pump was equal to the MHD in litres per minute ( Table 2). After entering the EPANET 2.2 input data [20] and successfully running the software, a hydraulic modeling of the networks was carried out. Simple controls (obtained after iterations to successfully run the software and achieve standard pressure heads at the nodes), as given in Table 3, were used to perform the hydraulic analysis.  The behavior of the nodes and pipes of the WDS can be evaluated using the TPIs (technical performance indices), which are indicative of the performance of the network with respect to the hydraulic parameters such as pressure head and flow. Two types of TPIs are used: TPI pressure [14] and TPI velocity [14]. The TPI values range from 0 to 1, '0' for poor service, and '1' for efficient service. The TPI values for both parts of the network were evaluated and compared with the standard values as follows: where Q i = nodal demand and (i) iterates over all the nodes of the network.
Here, P i = prevailing pressure head at the node, P min = 17 m, P max = 70 m, and TPI i = TPI. For the node i: where Q j = flow in pipe (j), and (j) iterates over all the pipes of the network.
Here, V i = prevailing flow velocity in pipe j, V min = 0.2 m/s, V max = 3 m/s, and TPI j = TPI for the pipe j, V mean = (V min + V max )/2.

Input Parameters for the Jal-Tantra Web System
The input data for nodes (elevations, base demand, and minimum pressures) and pipes (start nodes, end nodes, length, and roughness) is entered into the Jal-Tantra system. Commercially available pipe diameters and their cost per unit length according to the prevailing market rates are also provided. The input data is given in Tables 4 and 5.

Optimization of the Hydraulic Model
Objective functions: The first objective function is the maximization of the performance efficiency of the proposed WDS in terms of the hydraulics. The second objective function is the minimization of the capital cost of the network in terms of the pipe diameters and valves.
Decision variables: The decision variables are as follows: • P i = prevailing pressure head at the node. • V i = prevailing flow velocity in the pipe. • H = staging height of the storage reservoir.

Distribution Network Layout
A dead-end type of distribution network was drawn for both parts of the study area using EPANET 2.2, as shown in Figure 1a,b. The input data for nodes-elevation; base demand; and demand pattern-was provided for each node, and the input data for pipes-pipe diameter; roughness coefficient; and length-was provided. The pipe network was arranged for both WDS, keeping in mind the geometry of the study area; the dead-end type of distribution system was used for both parts of the study area. Moreover, the basic assumption of the Jal-Tantra web system of having no loops in the network has to be satisfied. 023, 5, x 6 of 11

Distribution Network Layout
A dead-end type of distribution network was drawn for both parts of the study area using EPANET 2.2, as shown in Figure 1a,b. The input data for nodes-elevation; base demand; and demand pattern-was provided for each node, and the input data for pipes-pipe diameter; roughness coefficient; and length-was provided. The pipe network was arranged for both WDS, keeping in mind the geometry of the study area; the dead-end type of distribution system was used for both parts of the study area. Moreover, the basic assumption of the Jal-Tantra web system of having no loops in the network has to be satisfied.

Contour Plan of Pressure Head at Nodes at 9 a.m. (Peak Demand Hour)
The contour plans given in Figure 2a,b, show the distribution of the pressure head across the layout of the network, parts I and II. Different colours representing the range of the pressure head are depicted by the scale provided in each map. The analysis indicated that the pressure head at all the nodes is greater than 17 m, thus the water at each node can reach up to a minimum of three stories of the building [18]. There is no need to use the PSVS. Moreover, the pressure head is also less than 70 m at each of the nodes in both parts of the WDS, so there is no need to install the PRVs at any of the nodes.

Contour Plan of Pressure Head at Nodes at 9 a.m. (Peak Demand Hour)
The contour plans given in Figure 2a,b, show the distribution of the pressure head across the layout of the network, parts I and II. Different colours representing the range of the pressure head are depicted by the scale provided in each map. The analysis indicated that the pressure head at all the nodes is greater than 17 m, thus the water at each node can reach up to a minimum of three stories of the building [18]. There is no need to use the PSV S . Moreover, the pressure head is also less than 70 m at each of the nodes in both parts of the WDS, so there is no need to install the PRVs at any of the nodes.

Colour-Coded Diagram of Velocity in Links at 9 a.m. (Peak Demand Hour)
The colour-coded velocity diagram given in Figure 3a

Colour-Coded Diagram of Velocity in Links at 9 a.m. (Peak Demand Hour)
The colour-coded velocity diagram given in Figure 3a,b, shows the range of velocities of flow in each pipe of network parts I and II. The scale showing the range of velocity represented by each colour is given for each of the networks. The analysis indicated that the velocity of flow for all the pipes in both networks is in the range of 0.25 m/s to 3 m/s [18]. Thus, there is no danger of erosion and deposition in the pipes during peak flow.

Time Series Plots
The temporal variation of the pressure head at the peak demand nodes and the water level in the tank are shown in Figure 4a,b, and the velocity in the pipes feeding the peak demand nodes is shown in Figure 4c,d. The pressure head lies within the recommended ranges with respect to all the nodes shown throughout the day. However, the velocity of flow in the pipes at the hours of minimum demand is less than 0.25 m/s, which can be attributed to the demand pattern adopted and the fact that the diameters have been selected to provide sufficient pressure head and cater to the demands at the peak

Colour-Coded Diagram of Velocity in Links at 9 a.m. (Peak Demand Hour)
The colour-coded velocity diagram given in Figure 3a,b, shows the range of ve ties of flow in each pipe of network parts I and II. The scale showing the range of vel represented by each colour is given for each of the networks. The analysis indicated the velocity of flow for all the pipes in both networks is in the range of 0.25 m/s to 3 [18]. Thus, there is no danger of erosion and deposition in the pipes during peak flow

Time Series Plots
The temporal variation of the pressure head at the peak demand nodes and the ter level in the tank are shown in Figure 4a,b, and the velocity in the pipes feeding peak demand nodes is shown in Figure 4c,d. The pressure head lies within the rec mended ranges with respect to all the nodes shown throughout the day. However velocity of flow in the pipes at the hours of minimum demand is less than 0.25 m/s, w can be attributed to the demand pattern adopted and the fact that the diameters been selected to provide sufficient pressure head and cater to the demands at the hour. Moreover, the discharge in the pipes at these hours is less, and thus there is n rious hazard of deposition in the pipes. However, a velocity less than 0.25 m/s ma

Time Series Plots
The temporal variation of the pressure head at the peak demand nodes and the water level in the tank are shown in Figure 4a,b, and the velocity in the pipes feeding the peak demand nodes is shown in Figure 4c,d. The pressure head lies within the recommended ranges with respect to all the nodes shown throughout the day. However, the velocity of flow in the pipes at the hours of minimum demand is less than 0.25 m/s, which can be attributed to the demand pattern adopted and the fact that the diameters have been selected to provide sufficient pressure head and cater to the demands at the peak hour. Moreover, the discharge in the pipes at these hours is less, and thus there is no serious hazard of deposition in the pipes. However, a velocity less than 0.25 m/s may be adopted to meet the hydraulic standards only if there is a suitable method of scouring the pipes [18]. adopted to meet the hydraulic standards only if there is a suitable method of scouring the pipes [18].

Performance Evaluation by TPI
The TPI pressure and TPI velocity for WDS at the peak hour were evaluated for both parts of the study area. The TPI pressure for WDS part I = 1, indicating 100% efficiency of the Environ. Sci. Proc. 2023, 25, 83 9 of 11 network in terms of pressure head at the nodes at the peak hour. The TPI velocity for WDS part I = 0.615, indicating 61.5% efficiency of the network in terms of the velocity of flow in the pipes at the peak hour. The TPI pressure and TPI velocity for WDS part II are 1 and 0.645, respectively. Indicating the efficiency in terms of pressure head at nodes = 100% and in terms of velocity in the links = 64.5%. The lower values of TPI velocity can be attributed to the lower values of velocity during off-peak hours due to the demand pattern adopted.

Cost Optimization by Jal-Tantra Web System
The optimized value of the cost of the WDS in terms of pipe diameters was achieved by minimizing pipe diameters while maintaining the minimum specified pressure head for the nodes, as given in Table 6. However, the total cost before and after the optimization is approximately the same because the minimization of diameters can be seen in a smaller number of pipes and equals a maximum of 50 mm. Moreover, there is an increase in diameter of about 25 mm in some pipe lengths. The smaller difference in the costs indicates the cost-optimal design of the WDS for the study area by EPANET 2.2.

Conclusions
In this work, an "optimal" solution for the design of WDS for the University of Kashmir has been proposed to replace the existing and worn-out WDS, which has already served its function for more than 60 years. The design not only includes the modelling of hydraulic parameters such as pressure head and flow velocity, but the components are also designed in such a way that the objective functions of maximum efficiency of hydraulic performance and minimum capital cost in terms of pipe diameters and valves are satisfied. The main highlights of the work indicating the optimal hydraulic design include the following:

•
The WDS is designed to have a minimum pressure head of 17 m at each node throughout the 24 h of a day, so that the water can reach up to a minimum of the third floor level of all the buildings. The pressure head at all the nodes in both parts of the network is such that 17 ≤ P i ≤ 70 (m), indicating the absence of points of deficiency or excess in the pressure head in the network. As such, there is no need to install the pressure sustaining valves (PSV) or the pressure regulating valves (PRV), resulting in a further reduction in the capital cost of the network.

•
The velocity of flow in all the pipes during the hours of maximum discharge is within the standard range of 0.25 ≤ V i ≤ 3.0 (m/s), such that there is no danger of erosion or deposition during such hours. However, during the few hours of least discharge, the velocity of flow is less than 0.25 m/s in some pipes. Therefore, a scoring arrangement may be provided to remove the depositions, if any.

•
The TPIs were evaluated for both parts of the WDS. The performance at peak hour with respect to pressure at nodes is 100%, while the performance with respect to velocity in pipes is 61.5% for WDS part I and 64.5% for WDS part II. The lower values of TPI velocity can be attributed to the lower values of velocity at off-peak hours due to the demand pattern adopted. • Finally, the cost-optimality of the network was checked by using the Jal-Tantra web system, where the diameter of a few pipes was reduced by a maximum of 50 mm, while it increased by about 25 mm for a few pipes. The difference in the cost of the WDS as per EPANET modelling and that obtained by the Jal-Tantra is negligible, indicating the cost-optimal design of the WDS by EPANET.