In Mediterranean countries, users are often equipped with private tanks, which provide a temporary water storage capacity, able to compensate service interruptions due to either scarcity or irregularity of water supply. In the presence of private water storage, water supply is no longer linked to users' consumption and network-operating conditions can be off-design, therefore specific models have to be introduced in simulation models of water distribution networks. Here, a new mathematical model is proposed that is able to reproduce a tank's emptying/filling cycles. Specifically, by means of experimental analysis, a hyperbolic tangent law was tested to reproduce the filling process for private tanks. The flow rate is calculated by means of the classical Torricelli law, in which the float valve emitter coefficient and the valve area are calculated using a function that takes into account the water level within the private tank. The comparison obtained through the mathematical model and those observed from experiments confirmed the ability of the model to predict the flow rate balance within private tanks. The results show that the model is suitable for any length of float valve branch. The mathematical system can be easily used in a transient model to correctly estimate the supplied demand.
The reliability of water distribution systems is critical for urban activities, because continuous operation of urban water supply has to be ensured at sufficient pressure under different design and operational situations. Models developed for water distribution networks (WDNs) design and management generally take into consideration customer demands as fixed a priori outflows concentrated in nodes, independent of nodal heads (demand-driven simulation methods), so they ensure a complete balance between available discharge and required discharge in demand nodes. Under pressure-deficient conditions, head-driven simulation methods, which connect the outflow to the nodal head, are essential for a realistic hydraulic analysis of WDNs.
Head-driven analysis requires the introduction of a nodal outflow discharge model, which links the available discharge to the nodal pressure. Over the past three decades, several authors have derived various relationships, which can be classified into two categories: discontinuous and continuous formulas.
Bhave (1981) was probably the first to introduce the concept of minimum required nodal head value for normal working conditions and suggested a relationship which adopts a [0,1] concept for the head-outflow relationship and admits that full demand is available for heads higher than the minimum required value, whereas there are no discharge results for heads lower than the minimum required value. Germanopoulos (1985) proposed a formula which sets minimum and maximum threshold values such that a constant flow rate is assumed for a nodal head with a higher than maximum required value, zero value is obtained for pressure with a lower than the minimum threshold value and an exponential law is assumed for an intermediate level of pressure at the node.
Reddy & Elango (1989) suggested outflows wholly dependent on residual heads according to Qj=Sj (Hj–Hmin)0.5, where Sj is a node constant, Hj represents the available total head at demand node j and Hmin represents the pressure head, below which the service at demand node j is unavailable.
Gupta & Bhave (1996) addressed some weaknesses in Germanopoulos's equation and, more recently, several researchers proposed new expressions aimed at improving precision (see Tucciarelli et al. 1999; Ackley et al. 2001; Tanyimboh et al. 2001; Tanyimboh & Templeman 2004, 2010; Tabesh et al. 2014; Volpe et al. 2014).
However, all these pressure/discharge relationships fail when the service pipe fills a tank which flattens the water demand. Either located on the roof or beneath the ground, private tanks are often considered as very inefficient devices; however, users are often equipped with private tanks because they provide a temporary water storage capacity which is able to compensate possible service interruptions. In fact, either when the network operates under pressure deficient conditions or when water supplies are inadequate to meet consumers' continuous demands – due to either scarcity or irregularity of water resources – discontinuous water distribution and water resources rationing are often used as the main measures to cope with water shortage (De Marchis et al. 2011; Freni et al. 2014; Puleo et al. 2014; Giustolisi et al. 2017). In this operative scheme, which is very common all along the Mediterranean coast, users compensate for water service interruptions by collecting water during serviced periods, and delivering it when public water service is not available via gravity or pumping systems, according to the position of the water storage.
Due to filling/emptying cycles, the presence of the tank not only modifies the demand profile of a normal domestic user (with respect to design conditions (Giustolisi et al. 2012)) but also makes it impossible to employ classic head/discharge relationships because now the water level in the tank has to be taken into account. In fact, under the assumption that the required nodal pressure is sufficient to supply the tank, the filling process is governed by the action of a proportional float valve that opens partially or completely as a function of the water level in the tank. When the tank is almost full, the water level does not fall enough to induce the valve to open completely, but it dampens the instantaneous water demand and affects the pressure-consumption law. Therefore, in WDNs equipped with private tanks, specific models for the inflow/outflow process of the tank have to be considered; these relationships must be able to correctly estimate the flow rate in two different scenarios. In the former, the tank is not full and the valve is completely open. In this case the discharge is calculated by considering a constant value of the emitter coefficient. In the second case the specific pressure/discharge relationship has to take into account the closure of the orifice at the end of the filling process.
Recently, Giustolisi et al. (2014) considered a linear opening of the orifice with volume. However, the orifice feeding inline tanks are controlled by floating valves which do not follow linear behaviour. Therefore, a more accurate discharge law is needed in order to correctly analyse the effect of inline private tanks in WDN models. To this end, Criminisi et al. (2009) suggested a model which links the continuity equation for the tank with a float valve emitter law.
De Marchis et al. (2016) validated the tank model by Criminisi et al. (2009) against experimental data, showing that the model presents some weakness at the start of valve opening and closing.
The experimental research shows the existence of a new mathematical formulation which is able to model the entire experimental data set, obtained by considering three different valve branches. The model can be easily implemented in hydrodynamic models which have to take into account the presence of private tanks, regardless of tank and valve geometry.
MATHEMATICAL SUPPLY MODEL OF INLINE TANKS
This work introduces a pressure-discharge relationship that is able to take into account the presence of local storage. Water tanks are installed to supply water in the event of water shortages or intermittent distribution. The private tanks are usually located on the roof or beneath the ground. This proposed law models the action of the floating valve during the valve closure period when the valve discharge area progressively decreases as a function of the tank water level.
As recently observed by De Marchis et al. (2016), on the whole the formulation by Criminisi et al. (2009) is able to capture the valve-closing phase, but some deviations are observed at the beginning and end of the valve closure.
The formulation was obtained through the best fit of the numerical results to a set of laboratory measurements.
In Equation (7) the coefficient m is used to calibrate the slope of the hyperbolic tangent law. On the other hand, the coefficient n is the coefficient related to the shape of the float valve area during the closure phase.
Tests were conducted at the Environmental Hydraulic Laboratory of the University of Enna Kore, Italy, where a small WDN reproduced using high-density polyethylene (HDPE 100) pipes of nominal pressure equal to 16 bar (PN16). The WDN has three loops, nine nodes and eleven pipes. All the pipes have the characteristics DN 63 mm, thickness 5.8 mm and approximate length 45 m, and are arranged in horizontal concentric circles with bends, with radius variable in the range of 1.6–2.0 m. This way, a small WDN with three loops was reproduced in actual scale. Figure 1 plots a schematic representation of the WDN and a snapshot of the laboratory.
Four pumps (P) supply the needed discharge from the water tanks to the upstream air vessel (AV) which behaves as a constant head tank, keeping the pressure at the inlet node of the network constant and equal to a prescribed value by varying the speed of the pumps (total water head in the range of 10–60 m). The experimental setup aims at simulating the operational condition of an actual secondary pipe network.
The flow rate is monitored by electromagnetic flow meters, whereas pressure cells are used to measure the water head.
A private tank is located in the roof of the laboratory at approximately 17.5 m above the network level, and it is connected to the system by means of a high-density polyethylene (HDPE 100 PN16) pipe, 30 m long, with a diameter of 12.7 mm. The tank filling process is governed by a float ball valve. Details on the water distribution system can be found in De Marchis et al. (2013, 2015a, 2015b, 2016). Here, attention was focused on the filling process of the private tank and especially on the law governing the float valve.
In the proposed analysis, data collected by De Marchis et al. (2016) have been extended, increasing the number of experiments. Specifically, three sets of experiments were carried out by modifying the length of the float valve branch. In the first set of simulations we installed a branch approximately 20 cm long (referred to as test case 1, TC1); in the second set of experiments the float valve was equipped with a branch of 30 cm (referred to as test case 2, TC2). Finally, in the third set of experiments a branch 40 cm long (referred to as test case 3, TC3) was investigated. In Figure 2 a schematic representation of the water tank and float valve relative to the three cases is reported. As shown in the figure, the length of the valve branch affects the water level at which the valve starts to close. In the TC1 the undisturbed water level inside the tank reaches a depth of hmin= 0.60 cm, whereas when the water depth is hmax= 0.65 cm the valve can be considered completely closed. On the other hand, in the TC2 the water depth is hmin= 0.56 cm, and in the TC3 the depth is hmin= 0.53 cm. In all these cases, the experiments have been carried out considering that at the beginning of each test the water tank was completely empty and is full when the maximum water level is reached and the float valve is completely closed.
The experiments are carried out considering the tank was completely empty in order to reproduce a typical scenario of the intermittent water supply systems. During the event of absence of service, in fact, the users collect water from the private tank that is almost empty when the service starts. During the supply period, the users collect water from the network, thus guaranteeing the private tank filling process is complete.
For each test case, seven experiments were conducted by changing the pressure from 2 to 5 bar, with a step of 0.5 bar. This pressure range ensures that the plastic pipes are not affected by deformations.
The water volume which flows into the tank was evaluated through the electromagnetic flow meter installed in the WDN (see Figure 1). The hydraulic head H used in Equation (5) to estimate Qup was measured by means of a piezoresistive pressure transducer (with an accuracy of 0.1%). The data was collected with a sampling rate of 1 Hz. Each experiment was repeated twice and the average value was considered.
Experiments and mathematical laws: hyperbolic law for Cv and constant av
Figure 3 compares data collected from the TC1 experiments with results obtained from the hyperbolic law given by Equation (5). Data reported hereafter have been averaged, thereby achieving one datum per minute. The pressure varies from 2 to 5 bar with a step of 0.5. Overall, as expected, the more the hydraulic head grows the more the flow rate Q increases. On the other hand, the time of the filling process becomes increasingly shorter as the pressure increases. Finally, the front of the curve, which is representative of the closure phase of the valve, progressively increases its slope. For a pressure of 5 bar, the front is almost vertical due to the rapid closure of the valve. It can be seen from Figure 3(a) (P= 2 bar) that the filling process is completed in approximately 45 minutes and the float valve starts to close after approximately 25 minutes, resulting in a closure period of approximately 20 minutes. On the other hand, as shown in Figure 3(g), when the network pressure is increased (P= 5 bar), the tank is completely filled in approximately 30 minutes.
The valve begins the closure phase after 15 minutes and the closing period continues for 17 minutes. The results reported in Figure 3(b)–3(f) behave coherently with the above description.
The comparison between the experimental data and the proposed hyperbolic laws clearly shows that although the hyperbolic law is able to reproduce the experimental data, some deviations occur during the final phase of its closure. This deviation increases with pressure, due to the more rapid closure of the valve.
Similar results can be observed for the cases TC2 and TC3, plotted in Figures 4 and 5 respectively. The comparisons between the three cases show that increasing the length of the valve branch increases the time needed to complete the tank filling. In the TC3 case, with a pressure of 2 bar (Figure 4(a)) the tank is completely filled in approximately 1 hour, whereas the valve starts to close after 28 minutes. Coherently, Figure 4(g) shows that when the network pressure is increased (P= 5 bar), the tank is completely filled in approximately 34 minutes. The valve begins the closure phase after 10 minutes and the closing period continues for 24 minutes. The filling process time can be a useful indicator for the water utilities to estimate the minimum time of service required to guarantee the filling of the reservoirs. Similar results are achieved in the valve branch of 30 cm (case TC2). Interestingly, when the length of the float valve is increased, the agreement between the hyperbolic law and the collected data increases. In the TC3 an almost perfect overlap between the mathematical law and the data is observed.
The best agreement between the curve obtained from Equation (5) and the experimental data was obtained calibrating the coefficient m and . The best fitting was verified through the minimization of the root mean square error (RMSE) and the Nash–Sutcliffe model (N-S). Table 1 reports the values of the coefficients obtained for all the experiments.
|Test case .||Coefficient .||2 bar .||2.5 bar .||3 bar .||3.5 bar .||4 bar .||4.5 bar .||5 bar .|
|Test case .||Coefficient .||2 bar .||2.5 bar .||3 bar .||3.5 bar .||4 bar .||4.5 bar .||5 bar .|
In Figure 6 the variations of the coefficient m and versus the water head are plotted. It is clear that all the values of the coefficient m have a linear trend, as clearly shown in Figure 6(a).
Finding parametric equations that are able to reproduce physical phenomena is of fundamental importance for the development of numerical models with the ability to study physical phenomena of practical interest. In this framework, the results presented show that the hyperbolic tangent law can be used in numerical codes to reproduce the WDN, in which users address water scarcity conditions through the use of local reservoirs (see De Marchis et al. 2011).
Experiments and mathematical laws: hyperbolic law for Cv and av
Even though the results above show the ability of the mathematical model to reproduce the experimental data and to capture a fast closure of the float valve, some deviations were observed when the valve closure is nearly complete, at least for the short branch; furthermore, the model needs the calibration of several coefficients (see Equations (8) and (9)). In order to improve and simplify the float valve emitter law, the model was modified taking into account a hyperbolic law to model the closure of the discharge area av. Basically, the discharge is calculated according to Equation (7) above.
Figures 7–9 show the comparison between the experiments and the new formulation of the float valve emitter law. Interestingly, the agreement is clearly improved when the discharge area is calculated with a hyperbolic law. The deviations initially observed for the valve branch of 20 cm disappear. Moreover the coefficients m and n have a constant value. In all the test cases and for all the pressures considered, the results presented in Figures 7–9 have been obtained by assigning m= 2.5 and n= 4. This result is clearly attractive for the development of numerical code to reproduce WDN.
In order to show the ability of the mathematical model to reproduce the experimental data, in Table 2 the achieved values of N-S coefficients are reported. All the cases are characterized by values of N-S higher than 0.99.
|Test case .||Model .||2 bar .||2.5 bar .||3 bar .||3.5 bar .||4 bar .||4.5 bar .||5 bar .|
|Test case .||Model .||2 bar .||2.5 bar .||3 bar .||3.5 bar .||4 bar .||4.5 bar .||5 bar .|
Private tanks greatly affect the hydraulic behaviour of the network. In fact, the reservoirs create a high level of supply at the beginning of the service period, reducing the pressure level on the network. Under such conditions the formulation of the specific law accounting for the existence of inline tanks must be used. The existing formulations fail during the closure phase of the float valve. Experiments have been carried out to reproduce the filling process of private water tanks, directly connected to the WDN. A hyperbolic tangent law is proposed for the emitter coefficient Cv and for the discharge area av of the head-discharge relationship. The mathematical model has been compared with existing float valve law and with laboratory experiments, carried out for different values of pressure and for three branches of different lengths. The results demonstrate the ability of the hyperbolic tangent law to reproduce the filling process both when the valve is completely open and during the valve closure period. The mathematical system can be used in a transient model to correctly estimate the supplied demand in the presence of private tanks equipped with float valves. Further efforts are required to validate the ability of the proposed hyperbolic tangent law in pressure-driven steady state models. Moreover, new experiments are under consideration which would take into account the consumption pattern during the filling process. This way it will be possible to give more insights into the effect of the private tanks and validate the hyperbolic tangent law.