Research on the Normal Use of a Plug Discharge Tunnel

,is study analyzed the normal use of an unusual flood-releasing tunnel with a plug dissipator. Firstly, normal physical model tests based on the Froude criterion (1 : 50) were finished. Secondly, depression physical model tests based on the Froude criterion (1 : 50) and cavitation similarity criterion were finished. ,irdly, 3-dimensional numerical simulation of flow field was finished, and free surface profile was captured, which was based on RNG k-ε two-equation turbulencemodel and VOFmethod.,e focus of this study is on the relationship between normal use and cavitation characteristics (e.g., pressure, turbulence kinetic energy, and cavitation number). ,e results show that lowering the reservoir water level, reduced by 20.41m at most, increases the risk of cavitation of a plug discharge tunnel, which means with the decrease of the flow cavitation number, the possibility of structural damage will increase dramatically, while reducing the outlet height can effectively raise the flow cavitation number, ensuring the safety of normal use. Under the conditions of free outflow, for the H1/e values of 4.45, 4.00, and 3.55, the conditions in which the tunnel meets the requirements of anticavitation are h/D≤ 0.42, h/D≤ 0.39, and h/D≤ 0.35, respectively. In addition, the discharge capacity of the tunnel is not significantly reduced with the lowering of outlet height, implying that operation under a low water head of the plug discharge tunnel, as low as 3.55 of H1/e in the test, is feasible. ,e results obtained in this study can serve as reference information in engineering design of the plug discharge tunnel.


Introduction
Rebuilding the diversion tunnel into the discharge tunnel is a cost-effective approach that can reasonably utilize waste resources and increase the operational flexibility of engineering operations. e plug is a new internal energy dissipator that realizes the goal of drastic energy dissipation through the sudden enlargement and contraction of the water-crossing section [1]. Unlike the orifice energy dissipator, the plug energy dissipator overcomes the drawback of easy cavitation at a sharp-edged orifice, enabling it to obtain an outstanding performance with regard to structural strength, energy dissipation effect, and antiresistance [2,3].
Flow recovering length is an important concern because it determines the distance of adjacent multistage plugs. Only when the plug spacing is greater than the flow recovering length can the energy dissipation rate be maximum. Many researchers have studied the shape design of the plug dissipator. Ai and Ding [4] researched reasonable distance between the upper and the lower plug in multistage plug discharge tunnel. Ai and Ding [5] also studied the effect of contraction ratio to the minimum wall pressure coefficient with RNG k-ε two-equation turbulence model and proposed a quantitative calculation formula. Li et al. [6] studied the way of adding vent holes in the recovering area of the plug to increase the possibility of anticavitation damage of the plug flood tunnel and obtained positive results. e results show that, by setting vent holes in the recovering area, air can be diffused to the two side walls, thereby effectively protecting the plug discharge tunnel. ese studies mainly focused on the influence of different plug geometries on the flow characteristics, such as wall pressure characteristics and aeration protection length.
At present, many hydraulic projects have adopted plug energy dissipator. e temporary weep hole of the Glen Canyon dam in the USA and the bottom outlet of the Mica dam [7] in Canada have both adopted the plug energy dissipator. e Mica dam rebuilt a diversion tunnel with a water head of 180 m into a permanent bottom outlet. Although the discharge is not large, its successful use has played a strong role in promoting the engineering application of the plug. In China, the first use of an orifice dissipator and its operational implementation was in the Xiaolangdi Project located in the Yellow River, which underwent a prototype test under low-water lever conditions with no observed cavitation [8,9]. However, strong cavitation noise was subsequently observed during normal operation. e Xiluodu Project located in the Jinsha River studied the flood plugging tunnel of the pressurized shaft, and it was found that its energy dissipation rate and anticavitation performance still meet the requirements under a flow rate as high as 2700 m 3 /s, indicating that the use of a plug energy dissipator for releasing larger discharge is feasible [10]. e Houziyan Hydropower Station located in the Dadu River used a multistage plug to transform the diversion tunnel into a flood discharge tunnel [11,12]. e atmospheric pressure test, numerical simulations, and decompression test together showed that the structure satisfies the anticavitation requirement, and the energy dissipation rate of the entire tunnel was also improved to a relatively high level.
ese engineering experiences show that plug dissipator has the advantages of high energy dissipation rate, safety, and simple structure geometry.
However, in view of the fact that the plug is still a new type of energy dissipator, it is often designed as an unusual flood discharge tunnel in engineering, which means it is only used in operation under extreme conditions. is causes the resources to be wasted and is not conducive to engineering practice and promotion of the plug energy dissipation. Since the plug energy dissipator must have sufficient wall pressure to avoid cavitation damage [13,14], the tunnel outlet elevation is generally designed to be relatively low so that the outlet could be submerged under the downstream water level, and the reservoir water level must be sufficiently high to ensure high pressure in the tunnel. Operation of the tunnel at a low water level may reduce the wall pressure or even give rise to cavitation. According to the analysis, reducing the height of the outlet can increase the pressure inside the tunnel and enhance the antiresistance of the flood discharge tunnel so that even relatively low upstream and downstream water levels can ensure the safety of the flood discharge of the tunnel. However, the decrease in the outlet height will affect the discharge capacity and energy dissipation rate of the tunnel. How to solve the contradiction between the normal operational safety and discharge capacity and the energy dissipation rate of the plug discharge tunnel is a problem worth exploring.
is paper attempts to answer this question. e numerical model is developed due to factors such as economy, portability, and time-saving by comparing it with physical experiments. e numerical model is very flexible for adjusting parametric parameters and provides complete information of flow field, such as velocity profile and pressure field. e main focus of the present study is to investigate the features of discharge capacity, pressure, energy dissipation rate, and cavitation number of a plug discharge tunnel.
is work can provide references to prevent or reduce the risk of cavitation in the plug discharge tunnel in hydraulic engineering.

Physical Structure of the Plug
Discharge Tunnel e object of this study is the plug type discharge tunnel of Houziyan Hydropower Station in China. e design of the plug discharge tunnel is shown in Figure 1. e dimensions of the inlet gate are 9 m × 11 m (width × height), and the altitude difference between the inlet and outlet floor is 93 m. e entire flood discharge energy dissipation system consists of 6 plugs. Two vertical plugs (plugs 1 and 2) with the heights of 15.0 m and 6.0 m, respectively, are arranged in the shaft section, the outlet diameters are 7.0 m and 8.5 m, respectively, and the spacing between the two vertical plugs is 24 m. e second vertical plug is connected to the top of a rectangular baffle chamber, and the dimensions of the chamber are 13.0 m × 15.0 m × 30.0 m (width × height × length). After the baffle chamber, a slope plug is placed (plug 3) with the length of 25.0 m and the outlet dimensions of 10.0 m × 4.5 m (width × height). ree horizontal plugs (plugs 4, 5, and 6) are located after the slope plug. e inlet diameter of all three plugs is 12.0 m, and their outlet diameters are 9.0 m, 9.2 m, and 9.8 m, respectively. e distance between plugs 3 and 4 is 60.0 m, and the spacing between the three horizontal plugs are both 40.0 m. e lower horizontal section is in the form of an arched section, with a straight wall height of 11.5 m and a top arch radius of 7.8 m. e outlet of the tunnel has a "curve pressure slope outlet" body shape, and the width of the pressure slope plane gradually narrows from 13.0 m to 11.0 m. e "curve pressure slope outlet" body shape is made up of a circular curve, which requires the horizontal distance limited to be 26.1 m, and the tangent line of curve end needs to be horizontal. e variation in the height h of the pressure outlet is 4.0∼7.5 m, and the top plate is connected by a single arc curve. e horizontal length of the curve is kept unchanged at 26.1 m, and the radius ensures that the angle between the endpoint tangent and the horizontal direction is zero.

Normal Test of Physical Model Experiments.
e scale of the normal model is 1 : 50 based on the Froude criterion for experimental investigations (Figure 2). e experimental operational conditions are shown in Table 1. e test numbers 1-8 in Table 1 are unusual operational conditions in which the upstream and downstream water levels are unchanged (H 1 /e � 5.40), and the outlet submergence degree H 2 /h increases from 2.76 to 5.18 with the decrease in the outlet height. e normal operational conditions 9-32 consider the downstream as the limiting condition, which means that the outlet is always free discharge. In the table, H 1 is the upstream hydraulic head of 49 m, 44 m, and 39 m, e is the height of gate, and H 2 is the downstream hydraulic head (test numbers 1-8 are 20.7 m, and test numbers 9-32 are equivalent to h, changing from 4.0 m to 7.5 m). e hydraulic diameter D of the lower horizontal section, as a constant, is 14.13 m, and D is defined as 4A 0 /C, in which A 0 (A 0 � 181.53 m 2 ) and C (C � 51.38 m 2 ) are the cross-sectional area and perimeter of a lower horizontal section, respectively.

2
Mathematical Problems in Engineering

Depression Test of Physical Model Experiments.
e scale of the depression test model is 1 : 50, which is based on the Froude criterion and the cavitation similarity criterion for experimental investigations. All the shape parameters are just the same with normal test above. e time-averaged cavitation number [15] can be defined as where p u /c and v 2 /(2g) are the pressure head and velocity head of each section, respectively. p a is the atmospheric pressure, and p v is the saturated vapor pressure. g is the gravitational acceleration, 9.81 m/s 2 . c is the unit weight, 9810 N/m 3 . v is the velocity.
To ensure safety, the local atmospheric pressure shall be the atmospheric pressure corresponding to the elevation at which the upper horizontal section is located [15]: where p v /c is quite sensitive to temperature. e prototype value adopts the water temperature of 20°C, that is, 0.24 m water column pressure.  Under similar conditions, the cavitation number in the decompression test should be equal to the cavitation number in prototype, that is [15], e subscript m (model) in the formula represents the relevant hydraulic parameters in the decompression test, and the subscript p (prototype) represents the relevant hydraulic parameters in the prototype.
erefore, under the premise of gravity similarity criterion, if the cavitation number is similar, the atmospheric pressure in the depression tank should be as follows [15]: e vacuum degree to be controlled in the depression tank is as follows [15]: where p al /c is the atmospheric pressure at the location of the depression tank during the test and p am /c is the vacuum pressure to be controlled in the depression tank under similar conditions. Relative vacuum degree can be expressed as follows [15]: e maximum relative vacuum degree of the depression equipment in the laboratory is close to 99%. In addition to satisfying similar condition, there is still considerable surplus. If no cavitation occurs, we will continue to improve the relative vacuum degree until we find the incipient cavitation number. e related depression test equipment is shown in Figure 3. e outlet height of 7.5 m is verified by the 1 : 50 decompression test, and the sound pressure signals are shown in Figure 4. In the effective analysis frequency range of 35∼160 kHz, the maximum sound pressure difference of the section of plugs 3-6 under similar vacuum pressure does not exceed 5 dB. However, experimentally, an occasional weak cavitation was observed at the exit of plugs 5 and 6. Although the cavitation was not collapsed in the sidewall, the situation was almost critical, and if it was considered for normal use, it is very likely that cavitation damage would occur. When vacuuming to 98.5% by the decompression test, the water flow noise exceeds the environmental noise by 5 dB, and the cavitation number σ at this time was used as the incipient cavitation number σ i of the different parts of the plug discharge tunnel. e incipient cavitation numbers of X � 220 m, X � 280 m, X � 300 m, and X � 480 m are 15.73, 11.40, 8.51, and 6.61, respectively.
Note that the reason for the choice of effective analysis frequency range of 35∼160 kHz is that when the frequency is lower than 35 kHz, the experimental results are greatly affected by the ambient noise, while cavitation usually occurs at high frequency, which means low frequency signals do not have a substantial effect on experimental results. In addition, 160 kHz is the maximum sampling frequency of the hydrophone.

Numerical Simulations.
Since the research object is a closed pressurized structure, the conventional physical model test cannot effectively measure the internal characteristics of flow field, and numerical simulations provide a good solution for this problem. Considering that the sudden expansion and contraction of the plug will cause a sharp change in the streamline, the three-dimensional RNG k-ε two equation turbulence model was used to simulate the turbulent flow [16,17]. e control volume method was introduced to discretize governing equations. e SIMPLER method [18], which has good convergence characteristics, was used for pressure-velocity coupling [19]. e no-slip boundary condition was used on the wall [20]. e standard wall function was chosen to specify the near wall regions of the flow. Experience shows that, for reinforced concrete, roughness height, K s , equal to 0.0025 m is relatively appropriate. And this research considers that the roughness on the concrete surface is evenly distributed, so the value of roughness constant, C s , is equal to 0.5. e governing equations are follows [21].
Continuity equation: Momentum equation: k (turbulent kinetic energy) equation: ε (dissipation rate of turbulent kinetic energy) equation: where t is the time and ρ and μ are the average density of the volume fraction and molecular viscosity, respectively. p is the pressure, and G k represents the generation of turbulent kinetic energy due to the mean velocity gradients. μ t is the turbulent viscosity, which can be deduced for turbulence intensity k and energy dissipation rate ε: where S is the modulus of the mean strain rate tensor. General constants η 0 , β, C μ , C 2ε , σ k , and σ ε of each equation are 4.38, 0.012, 0.0845, 1.68, 0.7179, and 0.7179, respectively. e VOF method was introduced to simulate the free surface. e continuity equations of the volume fraction are as follows [22]: where u i is the velocity components and x i is the coordinates (i � 1, 2, 3). j is the sum suffix. e air-water interface is tracked by solving the continuity equations. e sum of the water and air volume fraction was one in the controlling body. ρ w and ρ a are the density of water and air, respectively. μ w and μ a are the viscosity of water and air, respectively. In the present study, the simulated results will be used for the analysis on hydraulic characteristics of the flow in the pressure tunnel. e inlet and outlet are both defined as pressure boundaries, where the inlet water level is given by a userdefined function (UDF) and the outlet is given by the boundary condition of the open channel flows. As shown in Figure 5, the grid spacing from the upper horizontal section to plug 1 is 1 m, and the grid spacing is 0.5 m at the energy dissipation section and lower horizontal section alongside the flow direction. e grid spacing of each section is 0.1∼0.5 m, and the grid spacing of the outlet channel is 0.3∼0.5 m. Except for the unstructured grid with a grid spacing of 0.2∼0.5 m in each connection part, the rest uses a structured grid. Figure 6 shows the distribution of the axial velocity of the energy dissipation section in the unusual operational conditions of h � 7.5 m.
e inlet position of each plug is squeezed by the water flow and forms a large flow velocity gradient, and a vortex of a certain intensity appears at the entrance angle. e flow of water between the different layers and the vortex area is strongly turbulent. At the end of the plug, due to the sudden enlargement of the cross section, a certain length of the vortex area is formed between the jet and the surrounding water. e jet boundary layer and the vortex area flow vigorously and exhibit turbulent mixing, and they are also the core area of the energy dissipation of the flow caused by the plug energy dissipator.

Results and Discussion
e grid independence was tested using the grid convergence index (GCI) method, which is a recommended and widely used method for the estimation of uncertainty due to the discretization in computational fluid dynamic (CFD) simulations [23], with grid numbers of 401,562, 562,268, and 755,624. Figure 7 presents the pressure with different grid numbers and uncertainty of the pressure results on the vertical plug section, where the GCI value is plotted in the form of error bars. From Figure 7, pressure results of the selected three groups of grids are relatively small, and the maximum relative error does not exceed 4%, and the maximum uncertainty was less than 3.2% (Table 2). erefore, the number of grids in this study was reasonably set to 562,268. Figure 8 shows the pressure comparison between experimental and calculated values at the top of the tunnel after plug 3. e maximum relative error in Figure 8 is not more than 13%, and the relatively large error is located in the backflow section after each plug, which is due to the complicated water flow in this area. Good agreement between the experimental and calculated values was obtained, indicating that the numerical simulation results are reliable. With the decrease in h, the pressure of the entire lower horizontal section increases, which is beneficial for the antiresistance of the tunnel. Due to the sudden contraction of the flow crosssection before each plug, the flow velocity increases and the pressure decreases. Due to the decreasing inner diameter, the plug pressure is gradually reduced until the outlet position of the plug is suddenly enlarged. e flow Mathematical Problems in Engineering forms a recirculating vortex around the jet, and the wall pressure continues to decrease. However, as the flow develops, the pressure begins to rise gradually. e pressure changing trend of each plug is almost the same, but the pressure generally decreases partially after each plug. e pressure tends to be stable after the backflow of the last plug.

Discharge Capacity and Discharge Coefficient.
e discharge capacity is one of the most important hydraulic factors in discharge structures. is study attempts to enhance the anticavitation ability of plug by reducing the height of the outlet to expand the use of the plug discharge tunnel, but the impact of outlet height on the discharge capacity is a factor that must be considered. e formula [15] for calculating the discharge capacity of a plug tunnel is where A is the outlet area, H up is the reservoir water level, H down is the downstream water level, ξ i is the local head loss coefficient of each part, λ is the resistance loss coefficient along the way, and L is the length of the resistance zone along the way. e formula (13) shows that when the upstream and downstream water levels are fixed, the discharge is mainly affected by the area of the outlet. In addition, Tong Haihong [24] collected a large amount of model test data and found that the decrease in the outlet height will increase the discharge coefficient of the pressured flood tunnel to a certain extent. Figure 9 shows the discharge capacity Q and discharge coefficient μ of the tunnel at different outlet heights. Under the same operational conditions, the discharge capacity does not have a linear relationship with the outlet height. With the decrease in h, the discharge capacity gradually decreases and the discharge coefficient gradually increases. At the same outlet height, the discharge coefficient of the tunnel basically does not change with the operational conditions; that is, the discharge coefficient of the plug discharge tunnel is only related to the body shape of the plug and is independent of the operational conditions. Under unusual operational conditions, the area of the outlet was reduced by 46.67% (from 82.5 m 2 to 44 m 2 ), and the discharge capacity of the tunnel was reduced by only 13.93% (from 1393 m 3 /s to 1199 m 3 /s). is indicates that the influence of the outlet height on the discharge capacity of the plug discharge tunnel is not significant; that is, it is feasible to improve the anticavitation ability of the tunnel by reducing the outlet height. In addition, the discharge capacity of the entire flood discharge system of the hydropower project is generally somewhat larger than the design value. e outlet height can be determined according to the excess amount of the entire flood discharge system, which means that the discharge reduced by the plug discharge tunnel can discharge through other flood discharge buildings.

Pressure Characteristics of Tunnel.
e pressure of the tunnel roof is shown in Figure 10. It can be concluded that, under the same outlet height condition, the pressure of each part of the flood discharge tunnel is reduced with the decrease in the upstream reservoir water level. is means the normal use of the plug discharge tunnel would increase the risk of cavitation erosion. e reason for this phenomenon is that, as the reservoir water level decreases, the piezometric head along the tunnel decreases. Under the same operational conditions, with the decrease in the outlet height, the pressure of the plug and subsequent sections increases.
is is because as the height of the outlet decreases, under the same upstream and downstream water level conditions, the discharge capacity of the tunnel will decrease, which reduces the flow velocity in the tunnel. According to the energy conservation equation, the pressure on the wall of tunnel increases. However, the change in the pressure of the upper horizontal section (X � 40 m) is relatively low; that is, the change in the pressure of the tunnel caused by the outlet height is mainly concentrated in the plug and its subsequent sections, so the design of the structural strength of the tunnel should focus on the plug and its subsequent sections. In fact, due to the decrease in the outlet height, the reduction in the discharge capacity of the flood discharge tunnel will result in a decrease in the flow velocity in the tunnel. Although the pressure in the upper horizontal section is hardly affected, the cavitation number in the tunnel is increased; that is, reducing the height of the outlet has a positive effect on the anticavitation of the entire flood discharge tunnel.

Energy Dissipation Rate.
e energy dissipation rate is one of the most important parameters of the hydraulic characteristics of energy dissipator, and it is also the most important factor for evaluating the performance of its energy dissipation characteristics. e energy dissipation rate of the plug discharge tunnel can be calculated according to the principle of energy conservation. e energy conservation equation between any two ends of the tunnel is given by the following equation [15]: where z i , p i , and v i are the potential energy, average pressure, and average flow velocity of section i, respectively, and h w is  the total head loss between the two sections. We define the energy dissipation rate between any two sections of the plug discharge tunnel as follows [15]: where H is the water head difference between upstream and downstream and is given by H � 93 m + H 1 − H 2 in this study. Figure 11 shows the variation in the energy dissipation rate of each plug. e variation rate of the energy dissipation rate of the plugs at different levels is basically the same under different operational conditions. Figure 11(a) shows the variation in the energy dissipation rate in the unusual operational conditions with the change in h. Figure 11(b) shows the change in the energy dissipation rate of the whole plug sections with the change in H 1 /e. With the decrease in the outlet height, the energy dissipation rate of the plugs showed a decreasing trend, which was most pronounced for plugs 1-3, while those for plugs 4-6 showed a little change. is is because the total energy dissipation of the system is mainly due to plugs 1-3, which can also be observed from the turbulent energy distribution shown in Figure 12. e most turbulent area of the entire flood discharge tunnel is the tail area near plugs 1-3. e energy dissipation rate of the entire plug system is also reduced accordingly, and the operational conditions of the flood discharge tunnel have a little effect on the total energy dissipation rate of the plug system; that is, the energy dissipation rate of the plug discharge tunnel is mainly affected by the body shape of the plug. In contrast, the energy dissipation rate is only weakly dependent on the working conditions.
At the end of each plug, due to the suddenly expanding boundary conditions and the inertia of the water flow, the sudden expanding water flow forms a strong turbulence, a strong shear jet, and a large-scale vortex. In the jet boundary layer and the large-scale vortex area, the flow field is turbulent and blended, resulting in significant energy dissipation. Large values of turbulence kinetic energy also appear on the floor and surrounding corners at the end of the #2 Mathematical Problems in Engineering plug. is is due to the high speed jet from the vertical hole plug impacting the floor. Large-scale vortices, strong turbulence in the impact point, and vortex area cause huge energy loss.

Cavitation
Number. Transient cavitation number [15] can be defined as where p ′ is the transient pressure and ρ is the density of water. Based on a large amount of experimental data, Uberoi [25] found that ��� where u′ is the transient velocity. In the case of Gaussian distribution [15], p equals to 2 �� � p′ 2 most of the time, which en, where �� � u ′ 2 /v is the turbulence intensity. According to the predecessor's test results [15], the value of the turbulence intensity near the bottom of highly turbulent flow is close to 0.15, and then It can be seen that the transient cavitation number may be 0.06 lower than the time-averaged cavitation number.
For a plug discharge tunnel that is characterized as operating under extreme conditions, the transition from a submerged state to a free outflow will increase the discharge of the tunnel and increase the flow velocity in the tunnel. e decrease in the reservoir water level reduces the discharge  capacity and decreases the wall pressure ( Figure 10). ese two factors reduce the cavitation number and increase the possibility of cavitation damage in the tunnel. In this study, the outlet is considered as the limiting case; that is, the outlet is kept free. And the change in the cavitation number of different sections is obtained by reducing the water level of the reservoir, thus providing reference information for the normal use of the plug discharge tunnel. e transient cavitation number in each section is shown in Figure 13 for different operational conditions. Under the same operational conditions, as h decreases, the cavitation number gradually increases. is is because as the height of the outlet decreases, the flow velocity in the tunnel decreases and the pressure increases. According to the definition of the transient cavitation number, it can be known that the transient cavitation number of the flow increases. As the reservoir water level decreases, the decrease in the discharge capacity will result in the decrease in the velocity in the tunnel, while the wall pressure is also reduced. Generally, the transient cavitation number decreases as the reservoir water level decreases; that is, the normal use or reduction in the reservoir water level will increase the risk of cavitation erosion. e reason is that the reduction of the pressure head is lower than that of the velocity head, which results in a decrease in the transient cavitation number. When h/D is less than or equal to 0.53, the tunnel can meet the requirements of anticavitation under verification conditions. e requirements of antiresistance for H 1 /e of 4.45, 4.00, and 3.55, are h/D ≤ 0.42, h/D ≤ 0.39, and h/D ≤ 0.35 under the normal use conditions, respectively.

Conclusions
In this study, the normal use of the plug discharge tunnel is realized by changing the outlet height. e key hydraulic parameters are compared using model tests and numerical simulations. e results show the following: (i) e RNG k-ε two equation turbulence model can be used to perform numerical simulation of the complex plug discharge tunnel and obtain satisfactory results. (ii) e decrease in the outlet height does not significantly reduce the discharge capacity of the plug discharge tunnel. e discharge coefficient is mainly affected by the shape of plug geometry and is irrelevant with the operational conditions. (iii) e pressure of the upper horizontal section is mainly affected by the upstream reservoir water level. (iv) e energy dissipation rate of plugs 4-6 is relatively weakly affected by the outlet height, and the decrease in the outlet height mainly reduces the energy dissipation rate of plugs 1-3, leading to a decrease in the energy dissipation rate of the whole flood discharge tunnel. e energy dissipation rate of the entire flood discharge tunnel is also not related to the operational conditions.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interests.