A Theoretical Derivation and Comparison Method for the Optimal Location for Energy Dissipation Boxes

: In long-distance, high-elevation gravitational water supply systems, it is essential to incorporate energy dissipation to lower pipeline pressures. The energy dissipation box is a novel device for pressure reduction, extensively utilized in gravitational flow transition systems. Despite its appealing contribution, systematic selection methods still need to be used regarding its optimum location. Hence, this paper considers the sum-of-the-maximum hydrostatic pressure head (SMHPH) and derives an extreme location equation for the energy dissipation box (EDB) in the design stage, and then, it comprehensively accounts for the overcurrent capacity and proposes a theoretical comparison and selection method for the optimal location (OL) between the critical and extreme locations. A theoretical analysis with an engineering case study is conducted to compare the theoretical OL, and numerical simulations are carried out to compare the pressure protection effect of the same box cross-sectional area and initial water volume of the theoretical OL and other possible locations with different initial water depths. The results show that the location of the EDB in the design stage should comprehensively consider the overcurrent capacity and the SMHPH. On the basis of both conditions, the oscillation amplitudes of the pressure gradient beyond the box are significantly decreased as the OL approaches closer to the downstream pipe. If the initial water depth in the box is large, the EDB provides greater protection to the pressure head (PH) by decreasing the cross-sectional area and then maximizing the utilization of water depth in the box and decreasing the volume of the EDB.


Introduction
Constructing long-distance water supply systems allows for the transportation of water from water-rich areas to regions experiencing water shortages, effectively addressing severe water resource deficiencies [1].These systems typically employ either gravity flow or pumping station pressurization.Gravity flow systems are favored for their reliability and cost-effectiveness since they operate without the need for external power to maintain flow [2].In a long-distance gravitational water supply system (LDGWSS), when the hydrostatic pressure head (PH) approaches or exceeds the specified design pressure, sudden changes in flow rate can lead to overpressure events.Such hydraulic shocks can potentially deform or rupture pipelines and their components under high pressures, while low pressures may cause pipeline collapse [3,4].Consequently, pipelines are equipped with overpressure protection systems, including pressure relief valves (PRVs), surge tanks combined with PRVs [5,6], and valves paired with energy dissipation boxes (EDBs).
The PRV, an essential element in hydraulic systems, is typically positioned at the supply line inlet or at the highest-pressure location.It functions by automatically opening its valve port when the pressure exceeds a predetermined threshold, thereby maintaining a constant pressure within the controlled system or loop [7].The valve stability profoundly influences the performance of a pipe system, and thus, a nonlinear model was established to obtain the PRV system's actual dynamic characteristics and stability criteria [8].In addition, Izuchi [9] studied the stability of the PRV through experiments and identified that dynamic instability arose from the interaction between pressure wave propagation in the inlet pipe and movement of the valve disk.Li and co-workers [10], concentrating on a numerical analysis, observed that the axial force exerted by the fluid on the valve flap initially decreased as the valve opened, followed by an increase.
Several works have focused on designing and applying water hammer protective systems employing PRVs and surge suppressors.For instance, Foster [11] and Alireza [12] studied transient flow dynamics with a surge tank and PRV integrated into a pump valve reservoir system.These protective measures successfully mitigated pressure surges within the spiral case and prevented turbine overspeed.In addition, the flux split technology has led to an efficient numerical method with first-and second-order accuracy that simulates the water hammer problem in a comprehensive hydraulic network with PRVs and surge tanks [13].The numerical evidence showed that the protective devices had desirable properties.Various PRVs and surge tanks have been proposed to analyze the multi-level decompression in LDGWSSs [14,15], with the corresponding results demonstrating that the protective devices above could mitigate excess water head, minimize negative pressure, and assure the safe operation of the pipeline.
PRVs have been used successfully under various hydraulic conditions to mitigate harmful transient pressures in distribution systems.PRVs require low inertia or a highly responsive control system to promptly address detected overpressures [16].However, careful selection of the PRV is essential, as improper choices can compromise the system's ability to manage transients effectively [17].PRVs are often accompanied by intermittent significant energy loss and frequent vibration, significantly affecting their discharge capacity and safety pressure protection [8].Furthermore, PRVs function reactively, which limits their applicability in many hydraulic systems where broader or more cost-effective mitigation strategies may be preferred over localized, rapid protection provided by PRVs.
The EDB examined in this study is an innovative device designed for energy dissipation and reducing PHs.It integrates an energy dissipation valve (EDV) with the EDB, allowing water to pass through the valve structure to achieve multi-stage energy dissipation.This is because the EDB eliminates the requirement of a large water head and avoids the sustained functionality of the EDB with a narrow aperture.Meanwhile, EDBs have the function of surge regulation and storage, cutting off the transmission of the water hammer, and alleviating the pressure oscillation.Therefore, the LHGWSSs should employ a setup combining EDVs and EDBs.In [18,19], the authors investigated the regulation of valve operation for the hydraulic transition within the EDB by establishing a one-dimensional numerical model to simulate hydraulic transients.Most of the previous EDB locations are influenced by pipeline layouts and are blindly built in designated places; there is still a lack of systematic theoretical analyses for the optimal location.Therefore, motivated by the shortcomings above, in the condition of considering the overcurrent capacity following EDB installation, [20] theoretically derived the critical location equation and analyzed the magnitude of the water level (WL) with upstream location.The results indicated that the theoretical critical location offered better PH reduction protection for the identical valve response time.
However, the study above ignores the optimization of SMHPHs to reduce the initial PHs along the pipeline as much as possible in the design stage; the critical location of the EDB is not necessarily the optimal location.Therefore, considering the overcurrent capacity and the SMHPH, there still remains a need for comprehensive theoretical evaluation regarding the OL of the EDB.This paper derives the extreme location for the EDB on the basis of decreasing the SMHPH.Zhang [21] established a theoretical analysis model of the one-way surge tank location under the assumption that the centerline elevation was straight.Meanwhile, Zhang [22] also established a theoretical analysis model of the air valve location by dividing the actual water transmission system into sections; each section was assumed to be a straight line.Through the aforementioned method, the optimal Water 2024, 16, 2189 3 of 14 location of the one-way surge tank and air valve in water supply pipelines was presented, and the results were validated using a numerical model for one-dimensional hydraulic transients.At present, there is still a lack of theoretical equations for SMHPHs along pipes.Considering the intricate nature of pipelines in actual gravitational flow transition systems, which encompass both ascending and descending sections, the calculation amount is huge when the extreme location for the EDB is calculated repeatedly each time according to specific trends.Hence, this paper establishes the theoretical equation of SMHPHs by dividing the water transmission system into sections and ignoring the local topographic trend of each section.Then, this paper proposes a theoretical comparison and selection method for the optimal location (OL) between critical and extreme locations and analyzes the impact of the locations in protecting against the water hammer in an LHGWSS.
The major innovations of this paper are as follows: (1) Deriving the theoretical equations for determining the extreme location of the EDB under the condition of minimizing SMHPHs, (2) proposing a discrimination method for identifying the theoretical OL of the EDB considering the overcurrent capacity and SMHPHs, and (3) comparing and analyzing the effect of the EDB location on water hammer protection.

Theoretical Equations for the OL of the EDB
This paper adopts the combined working modes of EDBs and EDVs in an LHG-WSS.The EDB comprises a rear riprap, energy dissipation bottom sill, energy dissipation partition, and folding stilling board.Figure 1 illustrates the corresponding structural diagram, highlighting that the water flow energy is primarily distributed through mixing/turbulence/breaking.Energy dissipation occurs through the exchange of mass, energy, and momentum between the box, water, gas, and internal friction.Examining the structural configuration of the EDB reveals that despite its initial small size, it absorbs the majority of flow energies due to hydraulic friction.The second part occupies a larger volume within the EDB, mainly adjusting the flow patterns and stabilizing flow velocities, functioning as a surge chamber.
straight.Meanwhile, Zhang [22] also established a theoretical analysis model of the air valve location by dividing the actual water transmission system into sections; each section was assumed to be a straight line.Through the aforementioned method, the optimal location of the one-way surge tank and air valve in water supply pipelines was presented, and the results were validated using a numerical model for one-dimensional hydraulic transients.At present, there is still a lack of theoretical equations for SMHPHs along pipes.Considering the intricate nature of pipelines in actual gravitational flow transition systems, which encompass both ascending and descending sections, the calculation amount is huge when the extreme location for the EDB is calculated repeatedly each time according to specific trends.Hence, this paper establishes the theoretical equation of SMHPHs by dividing the water transmission system into sections and ignoring the local topographic trend of each section.Then, this paper proposes a theoretical comparison and selection method for the optimal location (OL) between critical and extreme locations and analyzes the impact of the locations in protecting against the water hammer in an LHGWSS.
The major innovations of this paper are as follows: (1) Deriving the theoretical equations for determining the extreme location of the EDB under the condition of minimizing SMHPHs, (2) proposing a discrimination method for identifying the theoretical OL of the EDB considering the overcurrent capacity and SMHPHs, and (3) comparing and analyzing the effect of the EDB location on water hammer protection.

Theoretical Equations for the OL of the EDB
This paper adopts the combined working modes of EDBs and EDVs in an LHGWSS.The EDB comprises a rear riprap, energy dissipation bottom sill, energy dissipation partition, and folding stilling board.Figure 1 illustrates the corresponding structural diagram, highlighting that the water flow energy is primarily distributed through mixing/turbulence/breaking.Energy dissipation occurs through the exchange of mass, energy, and momentum between the box, water, gas, and internal friction.Examining the structural configuration of the EDB reveals that despite its initial small size, it absorbs the majority of flow energies due to hydraulic friction.The second part occupies a larger volume within the EDB, mainly adjusting the flow patterns and stabilizing flow velocities, functioning as a surge chamber.Assuming that the regulating valve before the box is non-operational and the friction resistance and elasticity of the water in the pipeline are negligible, we use the harmonic vibration theory and obtain the oscillation period T and the magnitude of the WL in the EDB ∆Hst [20]: where L2 is the pipeline length behind the EDB, Ast is the cross-sectional region of the EDB, v0 denotes the initial pipeline flow velocity, f represents the pipe cross-sectional region, Assuming that the regulating valve before the box is non-operational and the friction resistance and elasticity of the water in the pipeline are negligible, we use the harmonic vibration theory and obtain the oscillation period T and the magnitude of the WL in the EDB ∆H st [20]: where L 2 is the pipeline length behind the EDB, A st is the cross-sectional region of the EDB, v 0 denotes the initial pipeline flow velocity, f represents the pipe cross-sectional region, and g indicates the gravity acceleration.Equation (1) describes the downstream pipeline's pressure fluctuation intensity.When the EDB is near the downstream reservoir, the pressure fluctuations and oscillation periods along the downstream pipeline decrease for a given valve closure.Consequently, for LHGWSS, it is advisable to position the EDB as near as possible to the downstream reservoir.
After installing the EDB, both conditions must be met: (1) ensuring the absence of a negative PH within the pipeline during standard operational conditions and (2) maintaining the hydrostatic PH within permissible limits per pipeline PH standards.and g indicates the gravity acceleration.Equation (1) describes the downstream pipeline's pressure fluctuation intensity.When the EDB is near the downstream reservoir, the pressure fluctuations and oscillation periods along the downstream pipeline decrease for a given valve closure.Consequently, for LHGWSS, it is advisable to position the EDB as near as possible to the downstream reservoir.

Theoretical Equations for the EDB's Critical Location
After installing the EDB, both conditions must be met: (1) ensuring the absence of a negative PH within the pipeline during standard operational conditions and (2) maintaining the hydrostatic PH within permissible limits per pipeline PH standards.

𝛽 𝛼
(2) Equations ( 2) and (3) represent the critical location equations for determining EDBs' locations, which satisfies the pipeline's current overload capacity requirement.When β = α, the EDB's location is where a line parallel to the hydraulic slope intersects the upstream pipeline at the overcurrent point C.

Theoretical Equations for the EDB's Extreme Location
Considering the topographic trend of the gravity flow system, it is difficult to fully represent the SMHPH in a mathematical manner.Firstly, the position of the extreme point for EDBs based on either section of the pipe, which can be assumed to be straight, is analyzed theoretically.The water hammer protection measures should ensure that the maximum PH does not exceed 1.3-1.5 times the maximum working PH [23].Therefore, 1 + m times the maximum working PH is accepted as the pipeline PH standard, m∈(0.3,0.5).
Figure 3 presents a schematic diagram of the hydrostatic PH along Section I of the pipe for EDBs.The shaded region, AI1, is the SMHPH along the pipeline for the EDB (m 2 ), HI1 is the upstream pipe PH of Section I (m), HI2 is the downstream pipe PH of Section I (m), LI1 is the pipe length in front of the EDB (m), LI2 is the pipe length behind the EDB (m), δI is the slope angle of Section I, the angle is positive when upstream is higher than downstream, and hA is water depth, which is assumed to remain constant across different locations.
Equations ( 2) and (3) represent the critical location equations for determining EDBs' locations, which satisfies the pipeline's current overload capacity requirement.When β = α, the EDB's location is where a line parallel to the hydraulic slope intersects the upstream pipeline at the overcurrent point C.

Theoretical Equations for the EDB's Extreme Location
Considering the topographic trend of the gravity flow system, it is difficult to fully represent the SMHPH in a mathematical manner.Firstly, the position of the extreme point for EDBs based on either section of the pipe, which can be assumed to be straight, is analyzed theoretically.The water hammer protection measures should ensure that the maximum PH does not exceed 1.3-1.5 times the maximum working PH [23].Therefore, 1 + m times the maximum working PH is accepted as the pipeline PH standard, m ∈ (0.3, 0.5).
Figure 3 presents a schematic diagram of the hydrostatic PH along Section I of the pipe for EDBs.The shaded region, A I1 , is the SMHPH along the pipeline for the EDB (m 2 ), H I1 is the upstream pipe PH of Section I (m), H I2 is the downstream pipe PH of Section I (m), L I1 is the pipe length in front of the EDB (m), L I2 is the pipe length behind the EDB (m), δ I is the slope angle of Section I, the angle is positive when upstream is higher than downstream, and h A is water depth, which is assumed to remain constant across different locations.According to Figure 3, HI2= HI1 + LIsinδI, where LI is the pipe length (m).Therefore, AI1 can be obtained as follows: Equation ( 4) is the theoretical equation for SMHPHs.When LI1 = 0, the EDB is located upstream of the pipe; when LI1 = LI, the EDB is located downstream of the pipe, namely, where AI2 and AI3 are the SMHPHs for the EDB, located at the beginning and end of the pipe (m 2 ).For a given pipeline system, AI1 is a second-order equation concerning LI1.Therefore, the derivative of Equation ( 4) can be expressed as follows: when dAI1/dLI1 = 0, one can obtain the following: where LI1−ex is the pipe length before the EDB located at the extreme location.The SMHPH along the pipeline is optimized to achieve minimal value at LI1−ex.
As shown in Figure 2, the pipe is divided into three sections (Section I, II, and III), and the EDB is located at Section I.By ignoring the local topographic trend of each section, we theoretically derive the theoretical equation for the SMHPH along the overall pipe based on a straight line connecting the starting and ending centerline elevations in each section.Therefore, A I → III can be obtained as follows: where LII is the pipe length of Section II (m), LIII is the pipe length of Section III (m), and δII and δIII are the slope angle of the pipe of Section II and III (m), respectively.When d A I → III/dLI1 = 0, one can obtain the following: For extending the extreme location to the general pipe in gravity flow (the pipe can be divided to n sections, index 1, …, m, …, n), one can obtain the following: According to Figure 3, H I2 = H I1 + L I sinδ I , where L I is the pipe length (m).Therefore, A I1 can be obtained as follows: 4) is the theoretical equation for SMHPHs.When L I1 = 0, the EDB is located upstream of the pipe; when L I1 = L I , the EDB is located downstream of the pipe, namely, where A I2 and A I3 are the SMHPHs for the EDB, located at the beginning and end of the pipe (m 2 ).For a given pipeline system, A I1 is a second-order equation concerning L I1 .Therefore, the derivative of Equation ( 4) can be expressed as follows: when dA I1 /dL I1 = 0, one can obtain the following: where L I1−ex is the pipe length before the EDB located at the extreme location.The SMHPH along the pipeline is optimized to achieve minimal value at L I1−ex .
As shown in Figure 2, the pipe is divided into three sections (Section I, II, and III), and the EDB is located at Section I.By ignoring the local topographic trend of each section, we theoretically derive the theoretical equation for the SMHPH along the overall pipe based on a straight line connecting the starting and ending centerline elevations in each section.Therefore, A I→III can be obtained as follows: where L II is the pipe length of Section II (m), L III is the pipe length of Section III (m), and δ II and δ III are the slope angle of the pipe of Section II and III (m), respectively.
When d A I→III /dL I1 = 0, one can obtain the following: Water 2024, 16, 2189 6 of 14 For extending the extreme location to the general pipe in gravity flow (the pipe can be divided to n sections, index 1, . .., m, . .., n), one can obtain the following: where L m−ex is the pipe length in front of the EDB located at Section m (m) and H m1 is the upstream pipe PH of Section m (m).Equation ( 11) is the general theoretical equation for the extreme location under the assumption that each section of pipe can be approximated as a straight line.Equation ( 11) can basically represent the SMHPH in a complex pipe under the condition of refining the section of pipe.

Comparison Method for OL of the EDB
The analysis presented above reveals that the OL of the EDB is influenced by overcurrent capacity and the sum of the hydrostatic PH.As depicted in Figure 4, point A and point B are the critical and extreme locations for the EDB, respectively, L 1−ov is the pipeline length before the EDB situated at the critical location, A i→j represents the sum of the maximum hydrostatic PH from index i to j, and L hor 1−k represents the horizontal length of pipe before the EDB at point k (k is the index).∆H ov and ∆H ex denote the internal PH between the hydrostatic PH and centerline elevation at critical and extreme locations, respectively, and ∆H ′ ov and ∆H ′ ex are the internal water PH between the piezometric head and centerline elevation at the critical and extreme locations.
where Lm−ex is the pipe length in front of the EDB located at Section m (m) and Hm1 is the upstream pipe PH of Section m (m).Equation ( 11) is the general theoretical equation for the extreme location under the assumption that each section of pipe can be approximated as a straight line.Equation ( 11) can basically represent the SMHPH in a complex pipe under the condition of refining the section of pipe.

Comparison Method for OL of the EDB
The analysis presented above reveals that the OL of the EDB is influenced by overcurrent capacity and the sum of the hydrostatic PH.As depicted in Figure 4, point A and point B are the critical and extreme locations for the EDB, respectively,  is the pipeline length before the EDB situated at the critical location, Ai → j represents the sum of the maximum hydrostatic PH from index i to j, and  represents the horizontal length of pipe before the EDB at point k (k is the index).ΔHov and ΔHex denote the internal PH between the hydrostatic PH and centerline elevation at critical and extreme locations, respectively, and ∆ and ∆ are the internal water PH between the piezometric head and centerline elevation at the critical and extreme locations.
According to Equation (12), the installation of the EDB before the critical location can be performed at any elevation of the water.To maintain a defined margin for internal PH behind the EDB, the constant in Equation ( 12) is typically set to 3-5 m.
Based on Equation (1), it is preferable to position the EDB as proximate to the lower reservoir as feasible to alleviate pressure fluctuations.Therefore, when the difference in the SMHPH with an EDB located at critical and extreme locations is negligible, the OL is preferably closer to the downstream.However, when there is a significant difference between the two locations above, the SMHPH should be considered to determine the OL if the EDB corresponding to the extreme location is located prior to the critical location, and fluctuations of the WL inside the box and along the pipeline behind the box can be alleviated by reducing the valve operation time If the extreme location is located downstream of the critical location, theoretically the OL is at the extreme location.However, the general trend of the pipe centerline elevations To prevent negative PHs alongside the pipeline in the design stage, the minimum WL at the EDB can be determined as follows: According to Equation ( 12), the installation of the EDB before the critical location can be performed at any elevation of the water.To maintain a defined margin for internal PH behind the EDB, the constant in Equation ( 12) is typically set to 3-5 m.
Based on Equation (1), it is preferable to position the EDB as proximate to the lower reservoir as feasible to alleviate pressure fluctuations.Therefore, when the difference in the SMHPH with an EDB located at critical and extreme locations is negligible, the OL is preferably closer to the downstream.However, when there is a significant difference between the two locations above, the SMHPH should be considered to determine the OL if the EDB corresponding to the extreme location is located prior to the critical location, and fluctuations of the WL inside the box and along the pipeline behind the box can be alleviated by reducing the valve operation time If the extreme location is located downstream of the critical location, theoretically the OL is at the extreme location.However, the general trend of the pipe centerline elevations in gravity flow are typically from higher to lower, thus necessitating measures to avoid the significant water depth; the theoretical extreme location should be located near the critical location.Certain end segment pipelines are disregarded while retaining the extreme location downstream of the critical location.The position of the extreme point with respect to the EDB with lower water depth is determined through a comprehensive comparison.
The criterion of SMHPHs between the critical location and extreme location downstream of former location for certain terminal segment pipelines is not taken into consideration, as discussed below.
(1) When L 1−ex = L 1−ov , the extreme location coincides with the critical location.
(2) When L 1−ex < L 1−ov , the comparison of SMHPHs for EDBs located at the extreme location and critical location is between A up→A − A up→B − A B→A and A B→down − ex , the OL of the EDB is at the extreme point where the SMHPH reaches the minimum.
The dimension of the EDB and fluctuation of the WL inside the box should be also comprehensively considered to determine the OL during the operational control stage.The theoretical comparative methods presented above can offer direction and serve as a reference for the initial selection of EDB positions in practical engineering applications.Moreover, the methods can be flexibly utilized based on specific project requirements.It is worth noting that numerical simulations must be employed to determine an appropriate size for the EDB and establish valve operation protocols to prevent overflow or leakage.

Impact of EDB Position on the Water Transition System
This section analyzes the LHGWSS illustrated in Figure 5, showcasing a pipeline with a total length of 48.09 km, a design flow rate of 2.1 m 3 /s, and a pipe diameter of 1.2 m.The WLs of the upstream and downstream reservoirs are 2017.50m and 1567.00 m, respectively, resulting in a head drop of 450.5 m between the two reservoirs.At the end of the pipe, two parallel-installed flow-regulating valves (DN800 and DN1200) are present.Employing the characteristic line method establishes a 1D hydraulic transient numerical model for the water transition system.Herein, to guarantee the stable and safe operation of the system, it is crucial to ensure that there are no negative PHs alongside the pipeline when closing the valves and that the maximum PH remains below the pipeline PH standard.Figure 5

Constant Flow Assessment of Gravity Flow System with No EDB
As illustrated in Figure 5, the piezometric head surpasses the pipe PH standard at L = 10,153 m (point D).From stake 10 + 153.00 to 11 + 753.00 and from stake 13 + 022.28 to stake 43 + 957.34, the water PH along the pipeline exceeds the maximum allowable positive PH (stake 0 + 00.00 to 6 + 033.00, 350 m; 6 + 033.00 to 11 + 753.00, 400 m; and stake 11

Constant Flow Assessment of Gravity Flow System with No EDB
As illustrated in Figure 5, the piezometric head surpasses the pipe PH standard at L = 10,153 m (point D).From stake 10 + 153.00 to 11 + 753.00 and from stake 13 + 022.28 to stake 43 + 957.34, the water PH along the pipeline exceeds the maximum allowable positive PH (stake 0 + 00.00 to 6 + 033.00, 350 m; 6 + 033.00 to 11 + 753.00, 400 m; and stake 11 + 753.00 to 48 + 103.18, 480 m).Notably, at stake 39 + 352.63, the disparity between the piezometric head and the maximum working PH achieves a peak of 133.19 m.Hence, it is crucial that an EDB is installed to lower the PH for LHGWSSs, thereby preventing extreme conditions such as pipe bursts and ensuring safe operation.

Theoretical Evaluation of the EDB's OL
Referring to the hydraulic slope angle shown in Figure 5, we determine that sinα = 0.002912 and cosα ≈ 1.Based on Equation ( 2), the critical EDB location satisfies the condition β = α, resulting in sinβ = 0.002912 and cosβ ≈ 1.To prevent the occurrence of negative PHs in the pipeline caused by the EDB, the overflow point C should be positioned at the highest point at the pipeline's end (stake 48 + 079.04 with the internal PH).Using Equation ( 2), a hypothetical straight line parallel to the hydraulic slope line is established between point C and the upstream pipeline.The intersection of this line with the centerline elevation (point A, stake 6 + 062.68) identifies the theoretical critical location of the EDB, with the internal PH at point A being 293.62 m.
According to the pipeline centerline elevation trend, the pipe is divided into three sections, as analyzed in Figure 5.The position of the extreme point for the EDB is theoretically analyzed based on a straight line connecting the starting and ending centerline elevations in each section.Considering the WL requirement in the EDB under overcurrent capacity, it is determined that the initial water depth of the EDB located in Section II exceeds 294.04 m, thus positioning the extreme location at Section I.In the design stage, the calculated water depth is initially selected to be identical to the initial water depth of the critical location (3.26 m).Table 1 compares the pipe parameters of each section.By introducing Equations ( 10) or (11), L I1−ex = 7782.68m, and the internal PH at point B (extreme location) is 268.93 m.From the theoretical analysis of the EDB location presented above, it can be observed that .44 m 2 , and ∆H ov L hor 1−ex − L hor 1−ov = 234,250.49m 2 .Therefore, according to the comparison method for assessing the OL of the EDB, the theoretical OL is at the extreme location.

Comparative Study of Protective Effects at Various EDB Locations
In this study, a DN1000 EDV is installed upstream of the box to ensure downstream overcurrent capacity.In all models, the reducing PHs of the EDVs are consistently 72.88 m to ensure the EDV remains fully open.At the critical locations (model a) and extreme locations (models b and c), considering the pipeline's capacity requirement, the water depths of the EDB are set to 3.26 m and 54.5 m, respectively, and are basically primarily situated along the same line parallel to the hydraulic slope.Figure 6 illustrates that, following EDB installation, the initial PH along the pipeline significantly decreases compared to the maximum allowable PH.Furthermore, the initial minimum PH along the pipe are 1.43 m, 2.99 m, and 2.99 m, respectively.There are no occurrences of negative PHs along the pipeline, satisfying the overcurrent capacity requirements in the design stage.Moreover, the protective effects of models a and b are compared under identical valve closure times and cross-sectional areas (25 m 2 ) but different initial water depths in the EDB.Additionally, for decreasing the dimension of the EDB, model a's protection effect is compared to model c's under the same valve operation and closure time, with identical initial water volume in the box (about 85.75 m 2 ) but differing cross-sectional areas (from 25 m 2 to 1.5 m 2 ).The valve closure protocols for the three models are detailed in Table 2, and the critical pressures along the pipeline for each model are documented in Table 3. Table 4 compares the EDB parameters and variations in water depth among the different models, while Figures 6-8 depict the protective effects against water hammer problems for models a, b, and c.According to Figure 5, the water surface lines of the EDB at the extreme and critical points are along a straight line parallel to the hydraulic slope line, the EDV and the EDB, resulting in a joint reduction of 287.41 m in heads for both three models.The EDV upstream of the box is completely open, reducing 72.88 m in the PH, and the PHs of EDBs are decreased by 214.53 m.As presented in Figure 6, model b has a significantly higher initial internal PH than model a (6062.68m versus 7782.68 m).In comparison, from 7782.68 m to the end, the initial internal PH in model b is less than in model a.The differences in the initial internal PH at the same stake behind the box remain at 0.84 m, namely, the hydraulic loss in the pipeline between the two locations.As presented in Table 2, when the EDB cross-sectional area and valve closing time are the same, the maximum PHs for models a and b are 304.06m and 344.9 m before the box and 442.29 m and 435.39 m behind the box, respectively.
In model b and c, the amplitude of maximum and minimum pressures after EDBs, particularly in regions with lower elevation at the centerline, is lower than in model a.This is because, according to Equation (1), the proximity of the EDB is lower to the reservoir in model b, resulting in a faster reflection of pressure waves induced by the closure of downstream valves, effectively reducing pressure oscillations along the pipeline.The highest internal water pressures in the models above do not surpass the pipeline's working pressure, and the lowest difference between the maximum value in model a and the pressure criteria is 37.71 m positioned at L = 39,426.45m at the point of minimum pipeline centerline elevation.Additionally, the minimum PH along the pipe in three models are 1.47 m, 2.26 m, and 2.24 m, respectively.Both models ensure there are no negative pressure heads and significant positive pressure heads because of the appropriate valve operation protocols.Table 3 highlights that the maximum PHs of models a, b, and c are 442.29,435.39, and 421.64 m, respectively.Additionally, the minimum PHs downstream of the EDB in models a, b, and c are 1.47, 3.81, and 3.93 m, respectively.Figure 7a reveals that in model c, the maximum pressure does not surpass its maximum working pressure.The model c exhibits lower maximum pressure compared to model b, while its minimum pressure is greater than in model b.In Figure 7b, the highest magnitudes of the maximum PH pointing downstream of the EDB in models a, b, and c are 115.14, 109.08, and 95.33 m, respectively.It is worth noting that reducing the cross-sectional region of the EDB from 25 m 2 to 1.5 m 2 results in model c having a maximum pressure amplitude 13.75 m lower than that of model b, significantly reducing pipeline pressure.This is because the WL magnitude velocity of the EDB varies inversely with its cross-sectional area.Therefore, a smaller crosssectional region for the EDB is needed to further alleviate the pressure oscillation.However, the smaller the cross-sectional area, the greater the WL fluctuation in the box, leading to the leakage or overflow of the EDB and consequently requiring higher EDB height requirements and increasing project costs.Table 3 highlights that the maximum PHs of models a, b, and c are 442.29,435.39, and 421.64 m, respectively.Additionally, the minimum PHs downstream of the EDB in models a, b, and c are 1.47, 3.81, and 3.93 m, respectively.Figure 7a reveals that in model c, the maximum pressure does not surpass its maximum working pressure.The model c exhibits lower maximum pressure compared to model b, while its minimum pressure is greater than in model b.In Figure 7b, the highest magnitudes of the maximum PH pointing downstream of the EDB in models a, b, and c are 115.14, 109.08, and 95.33 m, respectively.It is worth noting that reducing the cross-sectional region of the EDB from 25 m 2 to 1.5 m 2 results in model c having a maximum pressure amplitude 13.75 m lower than that of model b, significantly reducing pipeline pressure.This is because the WL magnitude velocity of the EDB varies inversely with its cross-sectional area.Therefore, a smaller crosssectional region for the EDB is needed to further alleviate the pressure oscillation.However, the smaller the cross-sectional area, the greater the WL fluctuation in the box, leading to the leakage or overflow of the EDB and consequently requiring higher EDB height requirements and increasing project costs.In model b and c, the amplitude of maximum and minimum pressures after EDBs, particularly in regions with lower elevation at the centerline, is lower than in model a.This is because, according to Equation (1), the proximity of the EDB is lower to the reservoir in model b, resulting in a faster reflection of pressure waves induced by the closure of downstream valves, effectively reducing pressure oscillations along the pipeline.The highest internal water pressures in the models above do not surpass the pipeline's working pressure, and the lowest difference between the maximum value in model a and the pressure criteria is 37.71 m positioned at L = 39,426.45m at the point of minimum pipeline centerline elevation.Additionally, the minimum PH along the pipe in three models are 1.47 m, 2.26 m, and 2.24 m, respectively.Both models ensure there are no negative pressure heads and significant positive pressure heads because of the appropriate valve operation protocols.
Table 3 highlights that the maximum PHs of models a, b, and c are 442.29,435.39, and 421.64 m, respectively.Additionally, the minimum PHs downstream of the EDB in models a, b, and c are 1.47, 3.81, and 3.93 m, respectively.Figure 7a reveals that in model c, the maximum pressure does not surpass its maximum working pressure.The model c exhibits lower maximum pressure compared to model b, while its minimum pressure is greater than in model b.In Figure 7b, the highest magnitudes of the maximum PH pointing downstream of the EDB in models a, b, and c are 115.14, 109.08, and 95.33 m, respectively.It is worth noting that reducing the cross-sectional region of the EDB from 25 m 2 to 1.5 m 2 results in model c having a maximum pressure amplitude 13.75 m lower than that of model b, significantly reducing pipeline pressure.This is because the WL magnitude velocity of the EDB varies inversely with its cross-sectional area.Therefore, a smaller cross-sectional region for the EDB is needed to further alleviate the pressure oscillation.However, the smaller the cross-sectional area, the greater the WL fluctuation in the box, leading to the leakage or overflow of the EDB and consequently requiring higher EDB height requirements and increasing project costs.
Figure 8 depicts the water depth variation across various EDB positions.Owing to the effect of valve timing, the EDB's WL undergoes an initial rise followed by a decline in all three scenarios.Figure 8a suggests that in the half-water hammer wave period of model b and c (t = L 1−down /c = 40,356.56/1000≈ 40.36 s), the water depth is still constant.Beyond t = 40.36s, the downstream positive wave reaches the EDB base, causing fluctuations in its water and altering the water depth ratio accordingly.According to Table 4, models a, b, and c show highest water depth amplitudes of 0.47 m, 0.43 m, and 4.54 m, respectively, with lowest depths of 1.25 m, 1.31 m, and 17.77 m.In Figure 7b, models a, b, and c show highest flow velocity in the EDB entering and exiting the pressurized pipeline of 0.084 m/s, 0.084 m/s, and 1.4 m/s.The flow velocity ratio undergoes an initial rise followed by a decline, and the initial, maximum, and minimum flow velocity ratio between models b and c are 0.6, 0.064, and 0.059, respectively.The flow velocity ratio remains nearly constant.
From Equation (1), it is noted that the closer the EDB is to the downstream reservoir, the smaller the WL fluctuation in the box is and the larger the oscillation period is.Model b places the EDB closer to the lower reservoir compared to model a, resulting in a smaller highest amplitude and longer fluctuation period of WL variation.Additionally, model a places the EDB closer to the upper reservoir compared to model b; this configuration allows for faster upstream replenishment of water into the EDB in model a when the water level drops under identical valve closure conditions, resulting in smaller lowest amplitude of WL variation (1.25 m versus 1.31 m).In addition, Equation (1) reveals that the WL magnitude of the EDB varies inversely with its cross-sectional area.After reducing the cross-sectional region of the EDB from 25 m 2 to 1.5 m 2 ; the magnitude of the EDB WL in model c is much more violent than that of model b (highest amplitudes of 0.43 m and 4.54 m, lowest amplitudes of 1.31 m and 17.77 m).
As depicted in Figure 8b, the flow velocity in models a, b, and c remain constant.After 120 s-the action time of EDVs-the flow velocity initially falls before stabilizing, and the flow velocity is measured at 0 m/s after the valves are completely closed.From the theory, it is noted that the flow velocity of the EDB is inversely proportional to its crosssectional area.Hence, after reducing the cross-sectional area, the flow velocity of the EDB in model c is much higher than that of model b until the valves are completely closed (highest velocities of 1.4 m/s and 0.084 m/s), and the flow velocities in models a and b remain the same, as does the cross-sectional area.It is worth noting that reducing the cross-sectional region of the EDB from 25 m 2 to 1.5 m 2 results in model c having a maximum pressure amplitude 13.75 m lower than that of model b, significantly reducing pipeline pressure.This is because the WL magnitude velocity of the EDB varies inversely with its cross-sectional area.The analysis above also elucidates that model c exhibits a more pronounced protective effect on the pressure amplitude along the downstream pipe compared to other models.

Figure 1 .
Figure 1.Structural diagram of a conventional EDB.

Figure 1 .
Figure 1.Structural diagram of a conventional EDB.

Figure 2
Figure 2 depicts the visual representation of the EDB position in an LHGWSS.∆H ′ A and ∆H ′ C are the internal water PHs between the piezometric head and centerline elevation at points A and C, respectively, L A→C denotes the pipe's length between points A and C, α and β represent the gradient angles of the hydraulic line and the angle formed by pipe points A and C, and h A is the water depth in the EDB measured in meters.

Figure 2
Figure 2 depicts the visual representation of the EDB position in an LHGWSS.∆ and ∆ are the internal water PHs between the piezometric head and centerline elevation at points A and C, respectively,  → denotes the pipe's length between points A and C,  and  represent the gradient angles of the hydraulic line and the angle formed by pipe points A and C, and hA is the water depth in the EDB measured in meters.

Figure 2 .
Figure 2. Visual representation of the theoretical, critical location of the EDB.

Figure 2
Figure 2 reveals that if the PH downstream of point D beyond the standard pipe PH, an EDB is installed at point A, situated between the upstream reservoir and point C (overcurrent point).The critical location of the EDB satisfies the condition outlined in [20]:

Figure 2 .
Figure 2. Visual representation of the theoretical, critical location of the EDB.

Figure 2
Figure 2 reveals that if the PH downstream of point D beyond the standard pipe PH, an EDB is installed at point A, situated between the upstream reservoir and point C (overcurrent point).The critical location of the EDB satisfies the condition outlined in [20]:

Figure 3 .
Figure 3. Schematic diagram of the theoretical extreme location of the EDB in Section I.

Figure 3 .
Figure 3. Schematic diagram of the theoretical extreme location of the EDB in Section I.

Figure 4 .
Figure 4. Schematic diagram of the comparative analysis for the EDB location in LHGWSS (L1−ex > L1−ov).To prevent negative PHs alongside the pipeline in the design stage, the minimum WL at the EDB can be determined as follows:

Figure 4 .
Figure 4. Schematic diagram of the comparative analysis for the EDB location in LHGWSS (L 1−ex > L 1−ov ).

2 Figure 6 .
Figure 6.Impact of EDB positions on pipeline pressures (models a and b).

Figure 6 .
Figure 6.Impact of EDB positions on pipeline pressures (models a and b).

Figure 7 .
Figure 7. (a) PH changes in models b and c, (b) PH changes at maximum PH point downstream of EDB in models a, b, and c.

Figure 8 .
Figure 8.(a) Impact of EDB position on water depth and water depth amplitude and on (b) flow velocity and velocity ratio in EDB.

Figure 8
Figure 8 depicts the water depth variation across various EDB positions.Owing to the effect of valve timing, the EDB's WL undergoes an initial rise followed by a decline in all three scenarios.Figure 8a suggests that in the half-water hammer wave period of model b and c (t = /c = 40,356.56/1000≈ 40.36 s), the water depth is still constant.Beyond t = 40.36s, the downstream positive wave reaches the EDB base, causing fluctuations in its water depth and altering the water depth ratio accordingly.According to Table 4, models a, b, and c show highest water depth amplitudes of 0.47 m, 0.43 m, and 4.54 m, respectively, with lowest depths of 1.25 m, 1.31 m, and 17.77 m.In Figure 7b, models a, b, and c show

Figure 7 .
Figure 7. (a) PH changes in models b and c, (b) PH changes at maximum PH point downstream of EDB in models a, b, and c.

Figure 7 .
Figure 7. (a) PH changes in models b and c, (b) PH changes at maximum PH point downstream of EDB in models a, b, and c.

Figure 8 .
Figure 8.(a) Impact of EDB position on water depth and water depth amplitude and on (b) flow velocity and velocity ratio in EDB.

Figure 8
Figure 8 depicts the water depth variation across various EDB positions.Owing to the effect of valve timing, the EDB's WL undergoes an initial rise followed by a decline in all three scenarios.Figure 8a suggests that in the half-water hammer wave period of model b and c (t = /c = 40,356.56/1000≈ 40.36 s), the water depth is still constant.Beyond t = 40.36s, the downstream positive wave reaches the EDB base, causing fluctuations in its water depth and altering the water depth ratio accordingly.According to Table 4, models a, b, and c show highest water depth amplitudes of 0.47 m, 0.43 m, and 4.54 m, respectively, with lowest depths of 1.25 m, 1.31 m, and 17.77 m.In Figure 7b, models a, b, and c show

Figure 8 .
Figure 8.(a) Impact of EDB position on water depth and water depth amplitude and on (b) flow velocity and velocity ratio in EDB.According to Figure 5, the water surface lines of the EDB at the extreme and critical points are along a straight line parallel to the hydraulic slope line, the EDV and the EDB, resulting in a joint reduction of 287.41 m in heads for both three models.The EDV upstream of the box is completely open, reducing 72.88 m in the PH, and the PHs of EDBs are decreased by 214.53 m.As presented in Figure 6, model b has a significantly higher initial internal PH than model a (6062.68m versus 7782.68 m).In comparison, from 7782.68 m to the end, the initial internal PH in model b is less than in model a.The differences in the initial internal PH at the same stake behind the box remain at 0.84 m, namely, the hydraulic loss in the pipeline between the two locations.As presented in Table 2, when the EDB cross-sectional area and valve closing time are the same, the maximum PHs for models a and b are 304.06m and 344.9 m before the box and 442.29 m and 435.39 m behind the box, respectively.In model b and c, the amplitude of maximum and minimum pressures after EDBs, particularly in regions with lower elevation at the centerline, is lower than in model a.This is because, according to Equation (1), the proximity of the EDB is lower to the reservoir in model b, resulting in a faster reflection of pressure waves induced by the closure of downstream valves, effectively reducing pressure oscillations along the pipeline.The highest internal water pressures in the models above do not surpass the pipeline's working pressure, and the lowest difference between the maximum value in model a and the pressure criteria is 37.71 m positioned at L = 39,426.45m at the point of minimum pipeline centerline elevation.Additionally, the minimum PH along the pipe in three models are 1.47 m, 2.26 m, and 2.24 m, respectively.Both models ensure there are no negative pressure heads and significant positive pressure heads because of the appropriate valve operation protocols.

Table 1 .
The pipe parameters of each section throughout the pipeline.

Table 4 .
Comparison of EDB parameters and variation in water depth for various models.

Table 2 .
Valve activation and deactivation guidelines for models a, b, and c.

Table 3 .
The PH across pipelines for various models.

Table 4 .
Comparison of EDB parameters and variation in water depth for various models.