Next Article in Journal
Investigation on Electrical and Thermal Performance of Glass Fiber Reinforced Epoxy–MgO Nanocomposites
Next Article in Special Issue
Fatigue Strength Analysis of a Prototype Francis Turbine in a Multilevel Lifetime Assessment Procedure Part I: Background, Theory and Assessment Procedure Development
Previous Article in Journal
Mineral Contents in Aboveground Biomass of Sedges (Carex L., Cyperaceae)
Previous Article in Special Issue
A Review of Different Aspects of Off-Earth Drilling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Hammer Control Using Additional Branched HDPE Pipe

1
Faculty of Building Services, Hydro and Environmental Engineering, Warsaw University of Technology, 00-653 Warsaw, Poland
2
Department of Mechanical Engineering and Mechatronics, West Pomeranian University of Technology, 70-310 Szczecin, Poland
3
Department of Technical Systems Operation and Maintenance, Faculty of Mechanical Engineering, Wroclaw University of Technology, 50-371 Wroclaw, Poland
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(23), 8008; https://doi.org/10.3390/en14238008
Submission received: 29 October 2021 / Revised: 19 November 2021 / Accepted: 22 November 2021 / Published: 30 November 2021

Abstract

:
In pressurised pipeline systems, various water hammer events commonly occur. This phenomenon can cause extensive damage or even lead to a failure of the pumping system. The aim of this work is to experimentally re-examine the possibility of using an additional polymeric pipe, installed at the downstream end of the main pipeline, to control water hammer. A previous study on this topic investigated additional polymeric pipes connected to the hydraulic system with a short joint section of the same diameter as the main pipeline. In the current research, a different method of including an additional pipe was considered which involved connecting it with a pipe of a smaller diameter than the main pipeline. Three additional HDPE pipes, with different volumes, were investigated. The performance of the devices was studied for hydraulic transients induced by both rapid and slow, manual valve closures. Experimental results show that the additional polymeric pipe can provide significant pressure surge damping during rapid water hammer events. As the valve closing time lengthens, the influence of the additional pipe on the maximum pressure increase is reduced. The additional HDPE pipe does not provide notable protection against hydraulic transients induced by slow valve closure in terms of reducing the first pressure peak. No relationship between the volume of the additional pipe and the damping properties was noticed. The observed pressure oscillations were used to evaluate a one-dimensional numerical model, in which an additional pipe is described as a lumped parameter of the system. The viscoelastic properties of the device were included using the one element Kelvin–Voigt model. Transient flow equations were solved with the implicit method of characteristics. Calculation results demonstrate that this approach allows one to reasonably reproduce unsteady flow oscillations registered during experiments in terms of the maximum pressure increase and pressure wave oscillation period.

1. Introduction

Water hammer is an undesirable pressure surge caused by a sudden change in flow velocity. Among engineers involved in designing pipelines, it is well-known that this phenomenon can cause extensive damage or even lead to a failure of the pumping system [1,2]. In practice, it is impossible to completely avoid the occurrence of severe hydraulic transients. Pressure pipeline systems are constantly subjected to flow disturbances and, thus, water hammer events commonly occur. For this reason, the efforts of many engineers and researchers are focused on developing effective methods to mitigate the destructive effects of this phenomenon.
One of the common strategies to manage extreme pressures involves optimising the operating procedures. Bazargan-Lari et al. [3] studied an optimal valve closing procedure for real time pressure control in pipes. Wan et al. [4] showed that it is significant to set a reasonable operational time difference in pump-valve systems in order to prevent backflow, reverse rotation, and overpressure. Optimal closure laws were also studied in pump-turbines systems by Yu et al. [5], Zhou et al. [6] and Pérez-Sánchez et al. [7]. This problem is also relevant for transient flow of compressible fluids. Subani et al. [8] analysed the flow of a hydrogen-natural gas mixture to determine relationship between pressure waves and different modes of closing and opening of valves.
Another method of protecting pipelines against water hammer events is to install surge control equipment at the transient-sensitive region of the hydraulic system [9]. The following devices are commonly employed to reduce dangerous pressure rises [10]: air chambers, surge tanks, flywheels, air inlet valves and pressure relief or pressure-regulating valves. Stephenson [11] developed a simple guide for the correct sizing of air vessels. The optimal selection of air chamber size was numerically and experimentally studied by Kim et al. [12]. The performance of modern water hammer control systems involving surge tanks was analysed by Wan et al. in [13,14]. The influence of surge tanks on hydropower plants was studied by Riasi et al. [15] and Vereide et al. [16]. Research on the behaviour of air valves during hydraulic transients by Balacco et al. [17] suggested that these devices should be designed using experimental investigations. The basic principles in designing pressure relief valves were described by Zhang et al. [18]. Kim and Kim [19] studied the performance of pressure relief valves in cavitation condition. All of the mentioned devices are usually costly. Several recent studies have been carried out on the influence of a combination of simple water hammer control devices, to mitigate hydraulic transients. Bettaieb and Haj Taieb [20] analysed the control of pressure oscillations by means of simple and economic devices. Mery et al. [21] conducted experiments revealing the effect of combining different water hammer control techniques on pressure waves.
Another tactic used to attenuate pressure waves is to modify the pipeline system in order to lower its effective bulk elasticity modulus. This can be performed by inserting a flexible tube into the pipeline [22], installing additional highly-deformable pipes in a pumping system, adding an in-line polymeric short section, branched penstock or by combining these techniques. Kubrak and Kodura [23] used a flexible gas-inflated tube inserted into the steel pipeline to control unsteady pressure oscillations. The experimental and numerical study performed by Pezzinga and Scandura [24] and Pezzinga [25] revealed that the ability of the additional polymeric pipe to dampen pressure waves depends mainly on the ratio between the volume of the device and that of the network system. The performance of redesigning existing steel pipelines to control severe pressure waves was investigated by Triki [26,27,28,29] who revealed that this technique may significantly enhance the reliability of hydraulic systems.
The aim of this research is to re-examine the possibility of using additional branched HDPE pipes to control water hammer in a steel pipeline. Previous papers on this topic [24,25] considered additional pipes connected to the hydraulic system with a short joint section of the same diameter as the main pipeline. This may not often be feasible for existing pipeline systems. To overcome this shortcoming, this paper presents a new approach to include additional polymeric pipe in the steel pipeline system, in order to mitigate high pressures. Furthermore, most experimental studies only focused on water hammer events initiated by rapid valve closure. Current research also examines the performance of the additional pipe during hydraulic transients induced with slow valve manoeuvring.
This paper is divided into six sections. Section 2 provides a brief mathematical description of unsteady flow in a pipeline system with the additional pipe installed. In Section 3, the experimental study is presented. The analysis of experimental results is Section 5. Section 5 examines numerical calculations simulating water hammer events in the pipeline system, with an additional pipe, induced by an instantaneous valve closure. Our conclusions are drawn in the final section. The flowchart of this study is presented in Figure 1.

2. Hydraulic Transient Model

Transient flow equations result from the conditions of continuity and dynamic equilibrium and can be expressed in the following form:
H t + c 2 g A Q x = 0
Q t + g A H x + f Q | Q | 2 D A = 0
where f is the friction factor (-), c is the pressure wave velocity (m/s), H is the piezometric head (m), Q is the discharge (m3/s), g is the gravity acceleration (m/s2), x is the space coordinate (m), t is time (s) and A is the cross-sectional area of the pipe (m2).
In order to take into account the influence of the additional pipe in the unsteady flow equations, a simplified model was used. It involves considering an additional pipe as a lumped parameter [24]. This approach allows for significant simplifications as the unsteady flow parameters in the additional pipe (Q and H) are only calculated at one spatial point. In this model, frictional losses and the inertia term in Equation (1) are neglected and a constant value of pressure is assumed in the additional pipe. This is acceptable for cases in which the velocity of the pressure wave in the additional pipe is negligible, in comparison with the period of the water hammer phenomenon in the main pipeline. The presented model is valid for pressure oscillations induced by instantaneous valve manoeuvring. The boundary condition of the device must satisfy mass conservation and it expresses the connection between the discharge and the deformation of the additional pipe, as well as the water inside it. By assuming that the additional pipe is fabricated of the elastic material, the boundary condition is:
Q = g c A P 2 d H d t W A P
where cAP is the pressure wave velocity in the additional pipe (m/s) and WAP is the volume of the additional pipe (m3).
If the additional pipe is fabricated of polymers, its viscoelastic behaviour must be taken into account and the boundary condition takes the following form:
Q = ( 2 d ε d t + g c A P 2 d H d t ) W A P
where ε is the unit strain in the additional pipe wall (-).
The mathematical description of the viscoelastic behaviour of pipe-wall material during hydraulic transients was thoroughly derived in [30,31,32]. Here, only crucial equations are given, for brevity. The distinctive feature of viscoelastic materials is that their total strain can be expressed as the sum of instantaneous and retarded strains:
ε = ε 0 + ε r
where ε0 is the instantaneous (elastic) strain and εr is the retarded strain (-).
According to the Boltzmann superposition principle, the total stress can be expressed as:
ε ( t ) = J 0 σ ( t ) + 0 t σ ( t t ) J ( t ) t d t
where J0 is the instantaneous creep compliance (Pa−1), J(t’) is the creep function at t’ (Pa−1) time and σ is the stress (Pa).
It is assumed that the total circumferential strain is:
ε = D A P D A P 0 D A P 0
where DAP is the inner diameter of the additional pipe at time t (m) and DAP0 is the inner diameter of the additional pipe in steady (pre-transient) flow (m).
Assuming that the pipe material is homogenous, isotropic, has linear viscoelastic behaviour for small strains, and its Poisson’s ratio is constant, the total circumferential strain is given by:
ε ( t ) = α 0 2 S A P 0 [ p ( t ) p 0 ] J 0 + 0 t α ( t t ) D A P ( t t ) 2 s A P ( t t ) [ p ( t t ) p 0 ] J ( t ) t d t
where p(t) is the pressure at time t (Pa), p0 is the initial pressure, J(t) is the creep compliance at time t (Pa), SAP is the additional pipe wall thickness (m) and α is the dimensionless parameter (dependent on the additional pipe diameter and constraints).
If the material is viscoelastic, its response to the applied stress can be described as the behaviour of a combination of elastic springs and viscous dashpots, which constitutes the so-called modified Kelvin–Voigt model. The creep function describing the viscoelastic properties of the n Kelvin–Voigt elements is:
J ( t ) = i = 1 n J i ( 1 e t τ i )
where i is the number of a single Kelvin–Voigt element (-), and n is the total number of Kelvin-Voigt elements (-).
It should be noted that Equations (3) and (4) were related to the control of a water hammer negative-surge framework. In such cases, in the initial phase of unsteady pressure oscillations, water discharges out of the additional pipe. Hence, the right-hand side of Equations (3) and (4) have a negative sign. In this research, the performance of an additional pipe was tested during positively induced pressure surges, which caused a certain volume of water to initially flow into the device. Therefore, the boundary condition for the viscoelastic additional pipe is:
Q = ( 2 d ε d t + g c A P 2 d H d t ) W A P
The differential Equation (10) must be simultaneously considered with the retarded strain component, described by Equation (8).

3. Experimental Study

The main goal of this experimental study was to record pressure changes during the valve-induced water hammer phenomenon in a steel pipeline with an additional HDPE pipe installed. The experiments were conducted on a simple pressure tank-pipeline stand, the diagram of which is presented in Figure 2.
During the experiments, the pipeline (labelled 1 in Figure 2) was fed from a pressure tank (labelled 2 in Figure 2) with a Bourdon pressure gauge (labelled 3 in Figure 2) connected into the water supply system (labelled 4 in Figure 2). An additional HDPE pipe (labelled 5 in Figure 2) was installed at the downstream end of the main pipeline just before the valve which initiated the water hammer phenomenon (labelled 6 in Figure 2). The method of connecting an additional pipe (AP) to the steel main pipeline (MP) requires additional comment. In previous studies on this topic [24,25], the connecting pipe had the same diameter as the MP. However, for existing pipelines, this method of connection of the device seems impractical. In the current research, the MP is smaller in diameter than the AP. For this reason, in all experiments, the AP was connected to the MP with an L-shape bend copper pipe with an internal diameter of 12 mm and a total length of 0.61 m (labelled 7 in Figure 1).
The steel MP with an inner diameter of 0.053 m was 48 m long and consisted of 24 sections rigidly connected to each other and fixed to the floor. Pressure changes were measured with the piezo-resistive sensors located at the end of the steel pipeline (labelled PM in Figure 2) and in the AP (labelled PA in Figure 2). Pressure sensors were connected to a laptop through a DAQ device with the dedicated LabView program. The valve, which induced a pressure surge, was equipped with a system which was able to register its opening angle as a function of time. This allowed the measurement of time of the valve closure. In order to shorten the valve’s manoeuvring time, it was partially closed during steady flow. Transient flow was induced by manual valve closing. The volumetric flow rate was regulated with a valve (labelled 8 in Figure 2) and measured with an inductive flow meter (labelled 9 in Figure 2). In order to assess the influence of volume of the AP on the ability to dampen pressure waves, three additional HDPE pipes of different lengths (with the same internal diameters and wall thicknesses) were used in the experiments. The material parameters of the APs are summarised in Table 1.
All three additional pipes used in the experiments are shown in Figure 3a, and the individual elements of each device are presented on the example of AP no. 1 in Figure 3b.
Each of the devices consisted of an HDPE pipe (labelled 1 in Figure 3b) and steel lids (labelled 2 in Figure 3b) connected to the pipe with steel rings (labelled 3 in Figure 3b), a vent valve (labelled 4 in Figure 3b), pressure sensor socket (labelled 5 in Figure 3b), shut-off valve (labelled 6 in Figure 3b) and connection pipe (labelled 7 in Figure 3b).
Water hammer runs were conducted in two variants: with and without AP installed. For each experimental variant, two types of water hammer were induced: rapid and slow manual valve closure. For each variant of the experiment, eight runs were carried out, which included three rapid water hammer and five slow water hammer events. Overall, 32 water hammer runs were conducted. In order to compare individual experimental results with each other, the measurements were carried out with care, to maintain a constant discharge during a steady flow of approximately 183 m3/h, which corresponded to an initial velocity of 0.39 m/s.
The collected data enabled the creation of graphs of pressure oscillations and valve opening angles as a function of time. Their analysis was used to obtain the water hammer parameters: maximum pressure head, pressure wave velocity and time of valve closure. Examples of the measurement results for rapid and slow valve closures are presented in Figure 4 and Figure 5.
As shown in Figure 4, the maximum pressure increase in the MP without AP installed was equal to 50.84 m. The orange plot in Figure 4 and Figure 5 represents the valve opening angle and does not correspond to values of pressure head in the vertical axis. In order to determine the time of valve closure which initiated the water hammer phenomenon, the time difference between the start and the end of valve closing was read (marked with green circles in Figure 4 and Figure 5). Experimentally obtained pressure oscillations are a harmonic function with a period equal to the pipeline period of the reflection time. To determine pressure wave speed, the time difference between minimum and maximum pressure was considered (marked with blue circles in Figure 4 and Figure 5). The main characteristics of all experimental runs are reported in Table 2.
v0 is the initial flow velocity, TC is the valve closing time, ΔHMP is the maximum pressure increase in the main pipeline (obtained with sensor PM) and ΔHAP is the maximum pressure increase in the additional pipe (obtained with sensor PA).
Table 1 demonstrates that, for the runs with rapid valve manoeuvring, it was closed between 0.12 s and 0.15 s. For experiments with AP, the pressure wave velocity did not exceed 500 m/s, and the time period of the pressure wave reflection time was longer than 2 L/c = 2·48/500 = 0.192 s, so the valve was closed faster than the reflected pressure wave travel time. For cases without AP, the highest recorded pressure wave velocity was 1200 m/s and period of the wave reflection time was 0.080 s, which is shorter than the valve closing time. However, the observed maximum pressure increases are consistent with those calculated using the Joukowsky formula. Moreover, thorough analysis of the observed valve closing angle, along with pressure oscillations as a function of time (Figure 4) shows that the maximum pressure increase occurs before the valve is fully closed. For this reason, all the water hammer runs with rapid valve manoeuvring were considered as being induced with instantaneous valve closures. These results are consistent with observations from the similar high-frequency pipeline system presented in [33].

4. Analysis of Experimental Data

In order to compare the performance of a particular AP on the maximum pressure increase in the MP at the downstream valve section, for each of the experiments, values of ΔHMP are presented as a function of valve closing time in Figure 6.
Figure 5 clearly shows that installing the AP has a significant damping effect on the initial stage of pressure oscillations (first pressure peak) during rapid valve closure. In the pipeline without the device, for the shortest valve closing times, the maximum pressure increases exceed 50 m. The unsteady flow oscillations observed during the tests with the AP are approximately 30 m. The lowest pressure increase was recorded for AP no. 2 (ΔHMP is between 26.80 m and 27.18 m) and the highest for AP no. 3 (ΔHMP is between 30.30 m and 31.22 m). The lowest attenuation of pressure wave induced by instantaneous valve closure was obtained for AP no. 3 and was equal to 50.84 m − 31.22 m = 19.62 m. The best performance was observed during experiments with device no. 2 and the attenuation was equal to 50.84 m − 26.80 m = 24.04 m. In other words, experiments results reveal that the damping rates of the additional pipes were between 38% and 47%. Interestingly, no relationship between the volume of AP and the damping properties was observed. This is likely due to the properties of the copper pipe which connects AP to the MP. The small diameter of the connecting section probably limits the AP’s capacity to attenuate pressure waves. It should be noted that in experiments with a sudden valve closure, consistent values of maximum pressure increase were obtained for each of the APs, which proves the reliability of the results.
As the valve closing time lengthens, the influence of the AP on the maximum pressure oscillations is reduced. The results of slow water hammer experiments showed that the maximum pressure increases in the MP (with and without the device) are practically the same. By comparing the maximum pressure increases obtained during different runs for valve closing times of approximately 1.5 s (15.89 m for run without AP, 14.26 m for AP no. 1, 14.78 m for AP no.2 and 17.74 m for AP no. 3), it is apparent that differences between these cases do not exceed 3.5 m. Thus, the additional HDPE pipe does not provide notable protection against hydraulic transients induced by a slow valve closure in terms of reducing the first pressure peak. To illustrate the influence of the AP on the unsteady flow oscillations at the downstream end of the MP during a rapid water hammer event, data with and without the AP were compared (Figure 7).
Figure 7 shows that a highly deformable AP provides significant pressure surge damping during a rapid water hammer event. Much more severe transients were induced in the main pipeline without the device, both in terms of maximum pressure increase and the duration of unsteady flow oscillations. The observed maximum surge damping rate in the main pipeline system with an additional pipe is 40%. There were strong disturbances in the pressure fluctuations in the MP with AP installed, which are typical of pipeline systems with sudden cross section changes [34]. Figure 8 and Figure 9 illustrate experimental pressure oscillations induced by slow valve closures.
As can be seen in Figure 8 and Figure 9, in the case of slow rapid water hammer events, the AP has a much smaller impact on the maximum pressure increase. For valve closing times between 0.77 s and 0.78 s (Figure 8), the maximum pressure increase in the main pipeline with the device (ΔHMP = 17.91 m) is 73% of the maximum pressure increase in the pipeline without the device (ΔHMP = 24.68 m). Increasing the valve closure time to approximately 1.5 s (Figure 8) causes the maximum pressure increase to be practically the same in the MP, with and without the device. However, it is apparent from Figure 9 that the AP reduces the duration of the phenomenon, as the pressure wave is dissipated and time dispersed. Based on Figure 8 and Figure 9 and Table 2, it can be seen that installing the additional pipe spreads out the period of pressure oscillations. In the mentioned cases, the reduction in the observed pressure wave velocity is 51% (828 m/s in run 4 and 425 m/s in run 12) and 55% (727 m/s in run 6 and 400 m/s in run 22). This is due to viscoelastic mechanical behaviour of the additional pipe and its response to pressure load during transients. In addition, in the case of slow water hammer runs, pressure disturbances do not come from the connection pipe.
So far, the data obtained from the pressure sensor installed in the main pipeline (sensor PM) were analysed. Turning now to observed pressure changes in the AP, it should be noted that maximum pressure increases in the MP are not the same as in the AP (Table 2). In order to compare these values, ratios of ΔHMPHAP are presented in Figure 10 as a function of valve closing time.
As illustrated in Figure 10, for cases with rapid valve closure, the maximum pressure increase in the MP is actually higher than the pressure increase in the AP. For all cases, extending the valve closing time to approximately 0.7 s to 1.5 s causes the observed maximum pressure increase in MP to be lower than in AP. Interestingly, it can be noted that, as the valve closing time further increases, the values of the ΔHMPHAP ratios tend to 1. A possible explanation for this effect is presented on the examples of individual experimental results. The largest differences between pressure increases in the MP and AP were observed for AP no. 2. The comparison of pressure oscillations observed in the MP and the AP for a rapid valve closure are presented in Figure 11.
The unsteady flow oscillations in MP and AP, observed during the water hammer run induced by slow valve closure, are presented in Figure 12.
As can be seen in Figure 11, the maximum pressure increase in the MP is higher than in the AP due to the pressure disturbances. Transient pressure waves are partially reflected by any change in the physical structure of the pipe [35]. During rapid water hammer events, the AP connected to the MP acts as a branch with a different diameter, in which the transmission and reflection of pressure waves occurs [36,37]. This effect is potentially disadvantageous for damping pressure surges. For example, it is known that an insufficient connector discharge capacity of the surge tank may intensify water hammer elsewhere in the pipeline [14]. Figure 12 shows that during less severe hydraulic transients, the maximum pressure increase in the MP is lower than in the AP. This is likely due to the fact that the viscoelastic effect affects low frequencies more severely [38]. As the pressure wave attenuates, local losses (due to diameter change in the connection pipe) cause energy dissipation and the local maximum pressure recorded by the sensor in the AP (Figure 3) is lower than in the MP. Figure 11 and Figure 12 show that, as the pressure wave dampens, the differences between the unsteady flow pressure oscillations in MP and AP disappear, hence the values of ΔHMPHAP ratios presented in Figure 10 tend to 1.

5. Numerical Solution and Model Validation

The aim of the conducted numerical analysis was to simulate unsteady flow oscillations in the MP (with the AP installed), induced by instantaneous valve closure. The numerical analysis of water hammer is obtained by approximating the solutions of the continuity and momentum equations. Several numerical methods are available to solve transient flow equations. In particular, method of characteristics (MOC) schemes are well-known by both researchers and practitioners [39]. In this paper, a popular fixed grid MOC approach was used to simulate water hammer in the pipeline with the AP installed. The presence of the device in the pipeline system was included in the form of a downstream boundary condition. An illustration of the discretisation scheme is shown in Figure 13.
As can be seen in Figure 13, the AP is located downstream, i.e., in the junction ns. Therefore, the characteristic equation is positive and yields:
Q n s t + Δ t Q n s 1 t + g A c ( H n s t + Δ t H n s 1 t ) + f | Q n s 1 t | Q n s 1 t 2 D A Δ t = 0
By substituting differential ratios in Equation (10) with finite differences, one obtains:
Q n s t + Δ t = 2 ε t W A P + H n s t + Δ t g W A P c A P 2 Δ t H n s t g W A P c A P 2 Δ t
Combining Equations (11) and (12) yields:
H n s t + Δ t = 2 ε t W A P H n s t g W A P c A P 2 Δ t Q n s 1 t g A c H n s 1 t + f | Q n s 1 t | Q n s 1 t 2 D A g A c g W A P c A P 2 Δ t
The retarded strain in Equation (8) is computed as a sum of each Kelvin–Voigt element [30,40]:
ε r t Δ t t = i = 1 n ε r t Δ t t = k = 1 n [ α D A P 2 S A P J k τ k ρ g ( H n s t Δ t H n s t = 0 ) ε i t Δ t τ i ]
The numerical approximation of each retarded strain in each node is equal to:
ε r i t Δ t = J i F t Δ t J i e Δ t τ i F t 2 Δ t J i τ i ( 1 e Δ t τ i ) F t Δ t F t 2 Δ t Δ t + e Δ t τ i ε r i t 2 Δ t
With the function F at time tt described as:
F t Δ t = α D A P 2 S A P J i τ i ρ g ( H n s t Δ t H n s t = 0 )
After calculating the pressure at the downstream end from Equation (13), along with Equations (14)–(16), discharge can be calculated from:
Q n s t + Δ t = Q n s t + g A c ( H n s 1 t H n s t + Δ t )
As mentioned previously, the presented numerical model holds for water hammer events induced by rapid valve closure at the downstream end of the pipeline. In order to validate the approach of a lumped parameter in the context of the AP connected to the MP via a conduit with a small diameter, computed results were compared with the experimental data recorded during run no. 17 (Table 2). Calculations were carried out with the space step of Δx = 1 m and the time step of Δt = 0.00082 s which enabled the Courant number to be obtained as close to unity as possible (Ca = 0.99), eliminating numerical dissipation and dispersion. The dimensions of the main pipeline system used in the experimental tests and parameters of AP no. 2 are listed in Table 1 and these were used as an input. The pressure wave velocity was determined on the basis of measurements in the pipeline without the device (c = 1200 m/s). The steady-state friction factor was calculated for an assumed pipe roughness of 0.08 mm using the Colebrook-White formula. For simplicity, calibration of the numerical model was performed using only one Kelvin–Voigt element. The values of retardation time (τ = 0.0447 s), creep compliance (J = 3.33 × 10 9 Pa−1) and pressure wave velocity in the AP (cd = 190 m/s) were estimated by a trial and error method, to minimise the mean squared error between the calculated and experimental data. Figure 14 illustrates the comparison between calculated and observed H(t) functions.
As can be seen in Figure 14, the computation results agree well with the experimental data, in terms of pressure oscillation periods. Differences in calculated and observed maximum pressure increase in the main pipeline are due to the pressure wave reflections and transmissions discussed previously. Calculation results show that this model can be useful in practice, as it allows estimation of the pressure increase resulting from rapid valve closure. The main advantage of this approach is that it is easy to implement in the form of a boundary condition and a relatively small number of parameters characterising the device. Use of the viscoelastic model allows the reproduction of dissipation and dispersion of the pressure wave. However, the presented numerical model has obvious limitations. It does not take into account the method of connecting the AP into the MP and, thus, it does not reproduce pressure disturbances in the MP that were observed during experimental tests. Another weakness is that the choice of values of creep compliance, retardation time and the number of Kelvin–Voigt elements applied in simulations, determines the obtained results in a significant way. For the correct description of the viscoelastic response of the pipeline, it is necessary to determine the coefficients of the model of the flexible material used [41]. This is a well-known problem in modelling the viscoelastic response of pipe-wall material in pressurised pipeline systems [42,43]. Values of viscoelastic parameters, along with values of pressure wave velocities (both in the MP and the AP), are very difficult to assess and should be estimated on the basis of calibration procedures.

6. Conclusions

In this paper, the reduction in water hammer pressure waves in a steel pipeline using additional HDPE pipe was analysed. A previous study on this topic investigated additional polymeric pipes connected to the hydraulic system with a short joint section of the same diameter as the main pipeline. In the current research, a different method of including an additional pipe was considered which involved connecting it with a pipe of a smaller diameter than the main pipeline. This way of including an additional pipe in the hydraulic system makes it much more convenient to implement the device in the existing systems.
In general, the following conclusions can be made:
  • The experimental tests showed that the additional polymeric pipe can be successfully used to attenuate a pressure surge during a water hammer event induced by rapid valve closure. Comparing the pressure oscillations in the pipeline with and without the device, it is apparent that an additional pipe significantly reduces the maximum pressure increase and reduces the duration of the phenomenon.
  • As the valve closing time lengthens, the influence of the additional pipe on the maximum pressure increase is reduced. Results of experimental runs showed that, during less severe hydraulic transients, values of maximum pressure increase in the pipeline, with and without the device, are similar. Thus, additional HDPE pipe does not provide notable protection against hydraulic transients induced by slow valve closure in terms of reducing the first pressure peak.
  • In all analysed cases, both for rapid and slow valve closures, the presence of a polymeric additional pipe causes a reduction in pressure wave velocity and lengthens the period of pressure oscillations.
  • In this study, the performance of three devices with different volumes was experimentally analysed. The apparent lack of correlation between the volume of the additional polymeric pipe and its ability to attenuate pressure waves can be attributed to the small diameter of the connection pipe compared with the diameter of the device. This way of connecting additional pipes to the main pipeline, although justified by practical reasons, makes the device less effective. Further data collection is required to evaluate the influence of the connection pipe on the performance of the additional polymeric pipes.
  • In order to simulate the water hammer phenomenon, a numerical model in which an additional pipe is treated as a lumped parameter was used. Experimental data was used to calibrate the mechanical properties of the pipeline system. The calculation results showed that this approach allows the reasonable reproduction of unsteady flow oscillations registered during experiments, in terms of maximum pressure increase and pressure wave oscillation period.

Author Contributions

Conceptualisation, M.K. and A.K.; methodology, M.K., A.K. and A.M.; software, M.K., A.M. and K.U.; validation, M.K., M.S. and P.B.; formal analysis, M.K., K.U., A.M. and M.S.; investigation, M.K and K.U.; resources, M.K., A.K. and K.U.; writing—original draft preparation, M.K. and A.M.; writing—review and editing, A.K., K.U., P.B. and M.S.; visualisation, M.K.; supervision, M.K. and A.M.; project administration, M.K.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Faculty of Building Services, Hydro and Environmental Engineering, Warsaw University of Technology (Grant no. 504/04608).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The code generated during the study and experimental data are available from the corresponding author by request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Across sectional area of the main pipeline (m2)
cpressure wave velocity in the main pipeline (m/s)
cAPpressure wave velocity in the additional pipe (m/s)
Dmain pipeline internal diameter (m)
DAPadditional pipe internal diameter (m)
DAP0inner diameter of the additional pipe at time t = 0
ffriction factor (-)
ggravity acceleration (m/s2)
inumber of a single Kelvin–Voigt element (-)
Hpiezometric head (m)
J(t) creep compliance at time t (Pa−1)
J(t’) creep compliance at time t’ (Pa−1)
LAPlength of the additional pipe (m)
ntotal number of Kelvin–Voigt elements (-)
Qvolumetric flow rate (m3/s)
p0initial pressure (Pa)
p(t) pressure at time t (Pa)
SAPadditional pipe wall thickness (m)
TCvalve closing time (s)
ttime (s)
v0initial average flow velocity (m/s)
vi,Lvelocity on the left hand side of the connection node (m/s)
vi,Rvelocity on the right hand side of the connection node (m/s)
WAPvolume of the additional pipe (m3)
WMPvolume of the main pipeline (m3)
xspace coordinate (m)
ΔHAPmaximum pressure increase in the additional pipe (Pa)
ΔHMPmaximum pressure increase in the main pipeline (Pa)
Δttime step (s)
Δxspatial step (m)
εunit instantaneous strain (-)
ε0instantaneous (elastic) strain (-)
εrretarded strain (-)
σstress (Pa)
τithe retardation time of i-th Kelvin–Voigt element (s)
Acronyms
APadditional pipe
DAQdata acquisition system
MOCmethod of characteristics
HDPEhigh-density polyethylene
MPmain pipeline

References

  1. Bergant, A.; Simpson, A.R.; Tijsseling, A.S. Water hammer with column separation: A historical review. J. Fluids Struct. 2006, 22, 135–171. [Google Scholar] [CrossRef] [Green Version]
  2. Ramos, H.; Covas, D.; Borga, A.; Loureiro, D. Surge damping analysis in pipe systems: Modelling and experiments. J. Hydraul. Res. 2004, 42, 413–425. [Google Scholar] [CrossRef]
  3. Bazargan-Lari, M.R.; Kerachian, R.; Afshar, H.; Bashi-Azghadi, S.N. Developing an optimal valve closing rule curve for real-time pressure control in pipes. J. Mech. Sci. Technol. 2013, 27, 215–225. [Google Scholar] [CrossRef]
  4. Wan, W.; Li, F. Sensitivity Analysis of Operational Time Differences for a Pump–Valve System on a Water Hammer Response. J. Press. Vessel. Technol. 2016, 138, 011303. [Google Scholar] [CrossRef]
  5. Yu, X.; Zhang, J.; Miao, D. Innovative Closure Law for Pump-Turbines and Field Test Verification. J. Hydraul. Eng. 2015, 141, 05014010. [Google Scholar] [CrossRef]
  6. Zhou, J.; Xu, Y.; Zheng, Y.; Zhang, Y. Optimization of Guide Vane Closing Schemes of Pumped Storage Hydro Unit Using an Enhanced Multi-Objective Gravitational Search Algorithm. Energies 2017, 10, 911. [Google Scholar] [CrossRef]
  7. Pérez-Sánchez, M.; López-Jiménez, P.; Ramos, H. PATs Operating in Water Networks under Unsteady Flow Conditions: Control Valve Manoeuvre and Overspeed Effect. Water 2018, 10, 529. [Google Scholar] [CrossRef] [Green Version]
  8. Subani, N.; Amin, N. Analysis of Water Hammer with Different Closing Valve Laws on Transient Flow of Hydrogen-Natural Gas Mixture. Abstr. Appl. Anal. 2015, 2015, 510675. [Google Scholar] [CrossRef] [Green Version]
  9. Boulos, P.F.; Karney, B.W.; Wood, D.J.; Lingireddy, S. Hydraulic Transient Guidelines for Protecting Water Distribution Systems. J. Am. Water Work. Assoc. 2005, 97, 111–124. [Google Scholar] [CrossRef]
  10. Chaudhry, M.H. Applied Hydraulic Transients, 3rd ed.; Springer: New York, NY, USA, 2014; ISBN 978-1-4614-8538-4. [Google Scholar]
  11. Stephenson, D. Simple Guide for Design of Air Vessels for Water Hammer Protection of Pumping Lines. J. Hydraul. Eng. 2002, 128, 792–797. [Google Scholar] [CrossRef]
  12. Kim, S.-G.; Lee, K.-B.; Kim, K.-Y. Water hammer in the pump-rising pipeline system with an air chamber. J. Hydrodyn. 2014, 26, 960–964. [Google Scholar] [CrossRef]
  13. Wan, W.; Zhang, B.; Chen, X.; Lian, J. Water Hammer Control Analysis of an Intelligent Surge Tank with Spring Self-Adaptive Auxiliary Control System. Energies 2019, 12, 2527. [Google Scholar] [CrossRef] [Green Version]
  14. Wan, W.; Zhang, B. Investigation of Water Hammer Protection in Water Supply Pipeline Systems Using an Intelligent Self-Controlled Surge Tank. Energies 2018, 11, 1450. [Google Scholar] [CrossRef] [Green Version]
  15. Riasi, A.; Nourbakhsh, A. Influence of Surge Tank and Relief Valve on Transient Flow Behaviour in Hydropower Stations. In Proceedings of the ASME-JSME-KSME 2011 Joint Fluids Engineering Conference, Hamamatsu, Japan, 24–29 July 2011; Symposia–Parts A, B, C, and D.; ASMEDC: Hamamatsu, Japan, 2011; Volume 1, pp. 1971–1977. [Google Scholar]
  16. Vereide, K.; Svingen, B.; Nielsen, T.K.; Lia, L. The Effect of Surge Tank Throttling on Governor Stability, Power Control, and Hydraulic Transients in Hydropower Plants. IEEE Trans. Energy Convers. 2017, 32, 91–98. [Google Scholar] [CrossRef] [Green Version]
  17. Balacco, G.; Apollonio, C.; Piccinni, A.F. Experimental analysis of air valve behaviour during hydraulic transients. J. Appl. Water Eng. Res. 2015, 3, 3–11. [Google Scholar] [CrossRef]
  18. Zhang, K.Q.; Karney, B.W.; McPherson, D.L. Pressure-relief valve selection and transient pressure control. J. Am. Water Work. Assoc. 2008, 100, 62–69. [Google Scholar] [CrossRef]
  19. Kim, H.; Kim, S. Optimization of pressure relief valve for pipeline system under transient induced cavitation condition. Urban Water J. 2019, 16, 718–726. [Google Scholar] [CrossRef]
  20. Bettaieb, N.; Taieb, E.H. Assessment of Failure Modes Caused by Water Hammer and Investigation of Convenient Control Measures. J. Pipeline Syst. Eng. Pract. 2020, 11, 04020006. [Google Scholar] [CrossRef]
  21. Mery, H.O.; Hassan, J.M.; EKaid, A.L. Water Hammer Mitigation by Air Vessel and Bypass Forward Configuration. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1094, 012052. [Google Scholar] [CrossRef]
  22. Wróbel, J.; Blaut, J. Influence of Pressure inside a Hydraulic Line on Its Natural Frequencies and Mode Shapes. In Advances in Hydraulic and Pneumatic Drives and Control 2020; Stryczek, J., Warzyńska, U., Eds.; Lecture Notes in Mechanical Engineering; Springer International Publishing: Cham, Switzerland, 2021; pp. 333–343. ISBN 978-3-030-59508-1. [Google Scholar]
  23. Kubrak, M.; Kodura, A. Water Hammer Phenomenon in Pipeline with Inserted Flexible Tube. J. Hydraul. Eng. 2020, 146, 04019054. [Google Scholar] [CrossRef]
  24. Pezzinga, G.; Scandura, P. Unsteady Flow in Installations with Polymeric Additional Pipe. J. Hydraul. Eng. 1995, 121, 802–811. [Google Scholar] [CrossRef]
  25. Pezzinga, G. Unsteady Flow in Hydraulic Networks with Polymeric Additional Pipe. J. Hydraul. Eng. 2002, 128, 238–244. [Google Scholar] [CrossRef]
  26. Triki, A. Water-hammer control in pressurized-pipe flow using an in-line polymeric short-section. Acta Mech. 2016, 227, 777–793. [Google Scholar] [CrossRef]
  27. Triki, A. Comparative assessment of the inline and branching design strategies based on the compound technique. J. Water Supply Res. Technol.–Aqua 2021, 70, 155–170. [Google Scholar] [CrossRef]
  28. Triki, A. Water-Hammer Control in Pressurized-Pipe Flow Using a Branched Polymeric Penstock. J. Pipeline Syst. Eng. Pract. 2017, 8, 04017024. [Google Scholar] [CrossRef]
  29. Triki, A.; Trabelsi, M. On the in-series and branching dual-technique–based water-hammer control strategy. Urban Water J. 2021, 18, 631–639. [Google Scholar] [CrossRef]
  30. Covas, D.; Stoianov, I.; Mano, J.F.; Ramos, H.; Graham, N.; Maksimovic, C. The dynamic effect of pipe-wall viscoelasticity in hydraulic transients. Part II—model development, calibration and verification. J. Hydraul. Res. 2005, 43, 56–70. [Google Scholar] [CrossRef]
  31. Keramat, A.; Tijsseling, A.S.; Hou, Q.; Ahmadi, A. Fluid–structure interaction with pipe-wall viscoelasticity during water hammer. J. Fluids Struct. 2012, 28, 434–455. [Google Scholar] [CrossRef] [Green Version]
  32. Keramat, A.; Kolahi, A.G.; Ahmadi, A. Waterhammer modelling of viscoelastic pipes with a time-dependent Poisson’s ratio. J. Fluids Struct. 2013, 43, 164–178. [Google Scholar] [CrossRef]
  33. Kodura, A. An Analysis of the Impact of Valve Closure Time on the Course of Water Hammer. Arch. Hydro-Eng. Environ. Mech. 2016, 63, 35–45. [Google Scholar] [CrossRef] [Green Version]
  34. Malesinska, A.; Kubrak, M.; Rogulski, M.; Puntorieri, P.; Fiamma, V.; Barbaro, G. Water Hammer Simulation in a Steel Pipeline System with a Sudden Cross-section Change. J. Fluids Eng. 2021, 143, 091204. [Google Scholar] [CrossRef]
  35. Brunone, B.; Meniconi, S.; Capponi, C. Numerical analysis of the transient pressure damping in a single polymeric pipe with a leak. Urban Water J. 2018, 15, 760–768. [Google Scholar] [CrossRef]
  36. Bergant, A.; Tijsseling, A.S.; Vítkovský, J.P.; Covas, D.I.C.; Simpson, A.R.; Lambert, M.F. Parameters affecting water-hammer wave attenuation, shape and timing—Part 1: Mathematical tools. J. Hydraul. Res. 2008, 46, 373–381. [Google Scholar] [CrossRef] [Green Version]
  37. Evangelista, S.; Leopardi, A.; Pignatelli, R.; de Marinis, G. Hydraulic Transients in Viscoelastic Branched Pipelines. J. Hydraul. Eng. 2015, 141, 04015016. [Google Scholar] [CrossRef]
  38. Duan, H.-F.; Ghidaoui, M.; Lee, P.J.; Tung, Y.-K. Unsteady friction and visco-elasticity in pipe fluid transients. J. Hydraul. Res. 2010, 48, 354–362. [Google Scholar] [CrossRef]
  39. Pal, S.; Hanmaiahgari, P.R.; Karney, B.W. An Overview of the Numerical Approaches to Water Hammer Modelling: The Ongoing Quest for Practical and Accurate Numerical Approaches. Water 2021, 13, 1597. [Google Scholar] [CrossRef]
  40. Bertaglia, G.; Ioriatti, M.; Valiani, A.; Dumbser, M.; Caleffi, V. Numerical methods for hydraulic transients in visco-elastic pipes. J. Fluids Struct. 2018, 81, 230–254. [Google Scholar] [CrossRef] [Green Version]
  41. Bogdevičius, M.; Karpenko, M.; Bogdevičius, P. Determination of rheological model coefficients of pipeline composite material layers based on spectrum analysis and optimization. J. Theor. Appl. Mech. 2021, 265–278. [Google Scholar] [CrossRef]
  42. Weinerowska-Bords, K. Viscoelastic Model of Waterhammer in Single Pipeline–Problems and Questions. Arch. Hydro-Eng. Environ. Mech. 2006, 53, 21. [Google Scholar]
  43. Weinerowska-Bords, K. Alternative Approach to Convolution Term of Viscoelasticity in Equations of Unsteady Pipe Flow. J. Fluids Eng. 2015, 137, 054501. [Google Scholar] [CrossRef]
Figure 1. Research flowchart.
Figure 1. Research flowchart.
Energies 14 08008 g001
Figure 2. Experimental setup: 1—main pipeline, 2—pressure tank, 3—Bourdon pressure gauge, 4—water supply system, 5—additional pipe, 6—valve initiating transient flow, 7—connecting pipe, 8—regulating valve, 9—inductive flow meter, PM—pressure sensor installed in main pipeline PA—pressure sensor installed in additional pipe.
Figure 2. Experimental setup: 1—main pipeline, 2—pressure tank, 3—Bourdon pressure gauge, 4—water supply system, 5—additional pipe, 6—valve initiating transient flow, 7—connecting pipe, 8—regulating valve, 9—inductive flow meter, PM—pressure sensor installed in main pipeline PA—pressure sensor installed in additional pipe.
Energies 14 08008 g002
Figure 3. Additional pipes used in experimental study. (a) Additional pipes no. 1,2 and 3, (b) Additional pipe no. 1: 1—HDPE pipe, 2—steel lids, 3—steel rings, 4—vent valve, 5—pressure sensor socket, 6—shut-off valve, 7—connection pipe.
Figure 3. Additional pipes used in experimental study. (a) Additional pipes no. 1,2 and 3, (b) Additional pipe no. 1: 1—HDPE pipe, 2—steel lids, 3—steel rings, 4—vent valve, 5—pressure sensor socket, 6—shut-off valve, 7—connection pipe.
Energies 14 08008 g003
Figure 4. Pressure oscillations recorded during run no. 3 (data obtained with sensor PM).
Figure 4. Pressure oscillations recorded during run no. 3 (data obtained with sensor PM).
Energies 14 08008 g004
Figure 5. Pressure oscillations recorded during run no. 5 (data obtained with sensor PM).
Figure 5. Pressure oscillations recorded during run no. 5 (data obtained with sensor PM).
Energies 14 08008 g005
Figure 6. Values of maximum pressure increase in main pipeline for each experiment as a function of valve time closure (data obtained with sensor PM).
Figure 6. Values of maximum pressure increase in main pipeline for each experiment as a function of valve time closure (data obtained with sensor PM).
Energies 14 08008 g006
Figure 7. Pressure oscillations recorded during runs no. 3 and no. 27.
Figure 7. Pressure oscillations recorded during runs no. 3 and no. 27.
Energies 14 08008 g007
Figure 8. Pressure oscillations recorded during runs no. 4 and no. 12.
Figure 8. Pressure oscillations recorded during runs no. 4 and no. 12.
Energies 14 08008 g008
Figure 9. Pressure oscillations recorded during runs no. 6 and no. 22.
Figure 9. Pressure oscillations recorded during runs no. 6 and no. 22.
Energies 14 08008 g009
Figure 10. Values of ΔHMPHAP ratios as a function valve closing time (data obtained with sensors PM and PA).
Figure 10. Values of ΔHMPHAP ratios as a function valve closing time (data obtained with sensors PM and PA).
Energies 14 08008 g010
Figure 11. Pressure oscillations observed during run no. 17.
Figure 11. Pressure oscillations observed during run no. 17.
Energies 14 08008 g011
Figure 12. Pressure oscillations observed during run no. 21.
Figure 12. Pressure oscillations observed during run no. 21.
Energies 14 08008 g012
Figure 13. Characteristic grid used in numerical calculations.
Figure 13. Characteristic grid used in numerical calculations.
Energies 14 08008 g013
Figure 14. Comparison between calculated and observed pressure oscillations.
Figure 14. Comparison between calculated and observed pressure oscillations.
Energies 14 08008 g014
Table 1. Parameters of the additional pipes used in experiments.
Table 1. Parameters of the additional pipes used in experiments.
AP NumberLAP
[m]
DAP
[m]
SAP
[m]
WAP
[m3]
WAP/WMP
[–]
10.400.0860.0120.00230.02
20.800.0860.0120.00460.04
31.500.0860.0120.00870.08
LAP—length of the additional pipe, DAP—inner diameter of the additional pipe, SAP—additional pipe wall thickness, WAP—volume of the additional pipe, and WMP—volume of the main pipeline.
Table 2. Data from experimental runs.
Table 2. Data from experimental runs.
RunVariant of ExperimentType of Valve Closure v0
[m/s]
TC
[s]
ΔHMP
[MPa]
ΔHAP
[MPa]
c
[m/s]
1no APrapid0.380.1449.79-1143
2no APrapid0.390.1450.33-1171
3no APrapid0.390.1350.84-1200
4no APslow0.400.7724.68-828
5no APslow0.390.8818.38-787
6no APslow0.411.5315.89-727
7no APslow0.392.9910.29-701
8no APslow0.404.874.85-238
9AP no. 1rapid0.380.1229.8326.58407
10AP no. 1rapid0.380.1429.2926.16407
11AP no. 1rapid0.390.1429.7426.63407
12AP no. 1slow0.380.7817.9123.10425
13AP no. 1slow0.391.4414.2617.55417
14AP no. 1slow0.391.9813.0715.18407
15AP no. 1slow0.393.598.659.39358
16AP no. 1slow0.384.439.189.99381
17AP no. 2rapid0.390.1327.1424.04393
18AP no. 2rapid0.390.1427.1824.17397
19AP no. 2rapid0.390.1426.8023.44393
20AP no. 2slow0.391.0317.4321.71410
21AP no. 2slow0.391.2416.7621.35403
22AP no. 2slow0.391.4814.7819.17400
23AP no. 2slow0.382.3513.8416.74414
24AP no. 2slow0.383.247.157.45333
25AP no. 3rapid0.380.1531.2231.66480
26AP no. 3rapid0.390.1530.3031.51480
27AP no. 3rapid0.390.1530.8731.11480
28AP no. 3slow0.381.4617.7421.99533
29AP no. 3slow0.381.9113.2014.71505
30AP no. 3slow0.382.468.739.10421
31AP no. 3slow0.393.466.326.32378
32AP no. 3slow0.393.966.116.13385
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kubrak, M.; Malesińska, A.; Kodura, A.; Urbanowicz, K.; Bury, P.; Stosiak, M. Water Hammer Control Using Additional Branched HDPE Pipe. Energies 2021, 14, 8008. https://doi.org/10.3390/en14238008

AMA Style

Kubrak M, Malesińska A, Kodura A, Urbanowicz K, Bury P, Stosiak M. Water Hammer Control Using Additional Branched HDPE Pipe. Energies. 2021; 14(23):8008. https://doi.org/10.3390/en14238008

Chicago/Turabian Style

Kubrak, Michał, Agnieszka Malesińska, Apoloniusz Kodura, Kamil Urbanowicz, Paweł Bury, and Michał Stosiak. 2021. "Water Hammer Control Using Additional Branched HDPE Pipe" Energies 14, no. 23: 8008. https://doi.org/10.3390/en14238008

APA Style

Kubrak, M., Malesińska, A., Kodura, A., Urbanowicz, K., Bury, P., & Stosiak, M. (2021). Water Hammer Control Using Additional Branched HDPE Pipe. Energies, 14(23), 8008. https://doi.org/10.3390/en14238008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop