Next Article in Journal
Implementation of Quantitative Resilience Measurement Criteria in Irrigation Systems
Next Article in Special Issue
Flow Hydrodynamic in Open Channels: A Constantly Evolving Topic
Previous Article in Journal
Ecological Risk Evaluation of Baihetan Dam Based on Fuzzy Hazard Quotient Model
Previous Article in Special Issue
SWAN Modeling of Dredging Effect on the Oued Sebou Estuary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydraulic Characteristics of Lateral Deflectors with Different Geometries in Gentle-Slope Free-Surface Tunnels

1
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
2
Power China Guiyang Engineering Corporation Limited, Guiyang 550081, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(17), 2689; https://doi.org/10.3390/w14172689
Submission received: 16 July 2022 / Revised: 24 August 2022 / Accepted: 26 August 2022 / Published: 30 August 2022

Abstract

:
The gentle-slope tunnel has been adopted in many high dams, and aerators are usually required for high operating heads. For such tunnels, the lateral deflector is superior to the traditional bottom aerator, which loses its efficiency due to cavity blockage and fails to aerate the sidewalls. However, unfavorable flow patterns such as water-wings and shock waves are induced by the lateral deflectors. To address this problem, two novel lateral deflectors are proposed, and their hydraulic characteristics are comparatively investigated together with the triangular deflector by means of model test and numerical simulation. The triangular deflector was revealed to form a wide cavity that allows for the free rise up of the water-wings inside the cavity, leading to the development of a buddle-type shock wave, whereas the two-arc deflector yields a jet with a fluctuating surface, which induces water-wings and further develops into diamond-type shock waves. In contrast, the cavity formed behind the two-arc deflector with a straight downstream guiding line is stabler and shorter, thereby restricting the development of the rising flow and preventing the formation of water-wings and shock waves. Moreover, the two-arc deflector with a straight guiding line exhibits higher energy dissipation capacities and thus is recommended in practical engineering design.

1. Introduction

Due to complex topographic and geological conditions [1] as well as material transport difficulties, high dams usually adopt arch concrete or local material types [2,3], the dam-body flood-discharge capacity of which is significantly limited [4]. Consequently, complementary flood release structures are of great importance for a sound design of the entire project [4]. In the past century, the tunnel spillway has stood out from other alternatives [4,5], such as chute spillways and shaft spillways owing to its high flood-release capacity, good terrain compatibility, and decent construction cost [6]. To avoid potential cavitation and denudation damage [1,4,5,7], the flow pattern and velocity in the spillway tunnel have to be strictly controlled. The ideal flow pattern should at least meet the following two requirements: (i) the flow is free-surface flow without shock waves and with adequate space between the aerated surface and the tunnel soffit [8]; (ii) the flow velocity is within 25 m/s [8] and the structure surface is carefully smoothed [8,9]. Under such conditions, technically no aerator is required. In the past decades, the sagging dragon tail tunnel has been proven to effectively fulfill the above requirement and thus has gained popularity in many high dams (e.g., Jinping, Baihetan, Xiaowan [5], Wudongde [10], Shuangjiangkou [11]) in China. The underlying reason for this kind of tunnel’s satisfactory hydraulic behavior lies in its longitudinal layout, which is reflected in its name ‘sagging dragon tail’. It comprises an overflow ogee weir (dragon head) close to the reservoir, a long gentle-slope free-surface tunnel in the middle (dragon body) and a short open and steep chute connected to a ski-jump bucket (the sagging dragon tail) at the end, as shown in Figure 1. With such a design, the flow acceleration is negligible inside the gentle-slope tunnel [4,5] and the high-velocity flow only occurs at the short steep chute, where the flow is usually highly aerated and the potential local structure destructions can be easily found and repaired once it occurs.
However, as the desired maximum operation head further increased, the flow velocity in the gentle-slope tunnel could exceed the threshold value of some 25 m/s [8], which consequently increases the cavitation risk of the tunnel [5,7,12] and the energy dissipation burden of the downstream chute. Under such conditions, aerators have to be installed to reduce cavitation risk [9,13,14], and restrictions resulting from flow pattern control have to be adequately considered. The main challenges of designing bottom aerators for gentle-slope tunnels are the low aeration efficiency caused by cavity blockage [6,12,15,16] and the cavitation protection of tunnel sidewalls [17,18]. Moreover, Hager discovered that shock waves could also be introduced by local changes in the bottom profile [14], such as offset aerators and bottom deflectors [16,19]. Therefore, to improve the aeration efficiency, lateral deflectors are usually installed together with the bottom aerators [15,17,18,20,21,22,23], of which a typical application is the fall-expansion aerators that constructed at the inlet of the free-surface tunnel behind the gate chamber [1]. It has been discovered that the lateral deflectors could not only serve as effective aerators [17,18,24] independent of the tunnel slope, but also help dissipate a certain amount of kinetic energy [24,25,26] and reduce the energy dissipation pressure of downstream structures. From this point, the hydraulic characteristics of the standalone lateral deflector (SLD, i.e., without bottom aerators) deserve further investigation and it is expected that the SLD could achieve decent hydraulic performance with regard to flow pattern and aeration.
To date, the commonest lateral deflector design is the triangular deflector whose up- and downstream endpoints are connected to the tunnel sidewall abruptly, and without vertical variation. The most noteworthy two of the few exceptions are probably the right-angled tetrahedrons and vertical plate deflectors investigated by Hager [26]. With such designs, water-wings and shock waves are usually induced by the lateral deflectors and thus introducing additional denudation risks. For example, Nie et al. (2006) discovered the jet formed behind the deflector impacts the sidewall and then is deflected upstream, resulting in partial blockage of the lateral cavity and the unstable swirling flow inside the lateral cavity plays a vital role in causing surface erosion [27]. The water-wings [18,20,23,24,28,29,30] formation was observed downstream of the impact region, which further rises up to the tunnel soffit [31] and could possibly develop into air pocket flow [16]. Shock waves [14,16,20,24,26,29,32] were also found to be induced by lateral deflection, which could propagate far downstream, mainly in a diamond shape with water crowns colliding alternatively in the middle and on the sidewalls; thereby, causing sidewall erosion with repeating impact and pressure fluctuation. Under highly aerated situations, shockwaves could cause overtopping or trigger conduits choke [14], which dramatically impede the air transportation and thereby harm for the aeration protection [23]. Existing studies mostly focus on the optimization of the size [21,24,26,28,30] and layout [1,26,29,30,31] of the traditional triangular deflector, and few attempts can be found to improve the flow pattern and energy dissipation by means of modifying the deflector geometries.
In this paper, two new deflector geometries are proposed for the sake of both flow control and energy dissipation. The hydraulic characteristics of the newly proposed deflectors and the traditional triangular deflector are comparatively investigated using the hybrid approach of model test and numerical simulation. Special focus is put on the flow pattern improvement. The energy dissipation characteristics of these three deflectors are also compared based on the simulation results.

2. Experiment Setup

The experiments were carried out with the physical model of a sagging dragon tail tunnel constructed at the State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University. In the experiment, water was supplied by a circulating system composed of an underground reservoir, pumps, and pipelines. The non-pressured section of the tunnel was a B = 18.75 cm wide, h = 17.75 cm deep flume inclined with a i = 3% bottom slope. The corresponding scale fell into the range of 1/24~1/80 with regard to the typical width being 4.5~15 m of the high unit-width-discharge free-surface tunnels (e.g., Yele, Xiaowan, Xiluodu) in China [8]. Experimental studies adopting similar geometric scale or model tunnel width can be found in [33,34,35,36]. The scale effect is considered acceptable as the main focus of this paper are the flow pattern and other macroscopical flow characteristics such as energy dissipation behaviors [36,37], which are much less sensitive to the model scale compared to two-phase flow characteristics such as air concentration, bubble size and air vent discharge [19,23]. The lateral deflectors were installed symmetrically (i.e., with identical geometry and streamwise location) on the sidewalls with our bottom aerators. The flume and deflectors are both made of plexiglass to provide a transparent view of the flow pattern. Three types of deflectors (Figure 1c) were investigated, the width b of which were all 0.94 cm and the streamwise length of deflectors A and B was 3.75 cm, while deflector C was further extended 1.25 cm downstream with a straight guiding line parallel to the side walls. The detailed geometries of the three types of deflectors are sketched in Figure 1. Type A is the traditional triangular deflector, deflector B is composed of two arcs tangent to each other: one being negative and has a radius of 5.625 cm and the other being positive and features a radius of 2.35 cm. As for deflector C, it is the same as deflector B except for the aforementioned additional straight line. Thus, the deflectors feature an identical contraction ratio of 10% (2b/B) and a tangent value of 1/4, and the dimensionless tail extension of deflector C with regard to the streamwise length of the curved section l is 1/3.
In the experiment, the inflow conditions were controlled at the cross-section 33.75 cm upstream of the deflector. All the measurements in this paper were conducted under the situation of inflow depth hin = 11.25 cm and volume flow rate Qin = 57.65 L/s, featuring a Reynolds number R e = Q i n h i n R υ = 1.2 × 105 (R represents the hydraulic diameter calculated as b h i n b + 2 h i n   and υ = 1 × 10−6 m2/s is the kinematic viscosity of water). This R e value being larger than 1 × 105 indicates the scale effect arising from viscous stress can be neglected according to [38]. The corresponding depth-width ratio of 0.6 in the experiment is a representative value for the real-world high-discharge flood tunnels [15,22,23,33,34,35,36] in China. In this study, the flow depth was measured using a fluviograph (accuracy ± 0.18 mm). The water discharge was monitored using an electromagnetic flowmeter (IFM4080K, Jiangsu Runyi Instrument Co., Ltd., Huaian, China), featuring accuracy of 0.1 L/s. The flow velocity was measured with a propeller-type flow meter (LS300-A, Fuzhou Lesida Information Technology Co., Ltd., Fuzhou, China) featuring accuracy of 0.01 m/s.

3. Numerical Models and Simulation Setup

Numerical simulations were performed in this study to obtain full-field velocities and turbulence properties for the analysis of the underlying mechanisms of flow pattern improvement and energy dissipation characteristics. The simulations were conducted using the commercial software FLOW-3D, which is claimed to have advantages over other opponents for free-surface flows and has been widely used for spillway and tunnel flows [39,40,41].

3.1. Governing Equations

FLOW-3D utilizes the one-fluid framework for free-surface flow modeling. This was achieved by using the Tru-VOF [42] technique to dynamically track the interface and, in the meantime, impose proper boundary conditions at the free surface. In this way, the computational cost is significantly reduced. The continuity and momentum equations are:
u i x i = 0
u i t + u j u i x j = 1 ρ p x i + f i + υ 2 u i x j x j x j ( 2 3 k δ i j 2 υ t S i j )
Here, u is velocity, p is the pressure, ρ and υ are the fluid density and kinematic viscosity, respectively. f stands for the body force, k and μt stand for the turbulent kinetic energy, and turbulent viscosity. δij is the Kronecker delta, and Sij is the mean rate of strain tensor calculated as 1 2 ( u i x j + u j x i ) .
The VOF equation for interface tracking reads:
α t + u i α x i = 0
where α is the volume fraction of the simulated fluid.
The RNG k-ε turbulence model [43] was adopted to account for the turbulence contribution to the time-averaged momentum transport, in which a transport equation was solved for the turbulent kinetic energy k and the turbulent dissipation rate ε, respectively:
k t + x i ( k u i ) = x j ( ( υ + σ k υ t ) k x j ) + P k ε
ε t + x i ( ε u i ) = x j ( ( υ + σ ε υ t ) ε x j ) + C 1 ε P k ε k C 2 ε * ρ ε 2 k
where υ t = C μ k 2 ε , P k = 2 υ t S i j S i j , C 2 ε * = C 2 ε + C μ η 3 ( 1 η η 0 ) 1 + β η 3 , η = ( 2 S i j · S i j ) 1 / 2 k ε , and model parameters are: σ k = σ ε = 1.39 , C μ = 0.085 , C 1 ε = 1.42 , C 2 ε = 1.68 , η 0 = 4.38 , and β = 0.012 .

3.2. Simulation Setup

The computational domain adopted the same coordinate system as the experiment and extends from x = −33.75 to x = 135 cm in the streamwise direction. A grid convergence study involving three mesh schemes was carried out in advance to select a proper mesh for the simulation. The details of the grid convergence study were shown in Appendix A. One structured mesh block was used to discretize the domain into 10.25 million cuboid cells, the majority of which feature an average size of 0.34 cm × 0.125 cm × 0.2 cm (x × y × z). To ensure the spatial resolution of the deflector geometry using the FAVOR® technique [44,45], the computational domain was rotated by the flume slope of 3% to be aligned with the x-coordinate and the gravity vector g was adjusted accordingly. Moreover, the domain was locally refined to 0.188 cm × 0.094 cm × 0.188 cm in the region of −7.5 < x < 30 cm, −10 < y < −6.25 cm and 6.25 < y < −10 cm, and 0 < z < 12.5 cm.
The outlet boundary was configured with the outflow boundary condition (BC), implying a zero-gradient condition since the outflow is supercritical. The bottom and sidewalls were set as non-slip walls with equivalent roughness ks = 0.015 mm according to the plexiglass surface properties. The top boundary is specified with fixed relative pressure p = 0 (i.e., atmosphere pressure). As for the inlet boundary, a pre-simulation of the entire physical model (i.e., from the reservoir to the downstream cushion pool) was conducted first, and then all the flow parameters at x = −33.75 cm were mapped onto the inlet of the short domain using the grid overlay BC. The illustrative diagram of the simulation setup is shown in Figure 2.

4. Results and Discussion

4.1. Model Validation

The simulated and measured water surface profiles at two longitudinal slices are comparatively shown in Figure 3a. It can be seen that the simulation results are in good agreement with the experimental data. Only small deviations can be found near the deflectors (−5 < x < 10 cm at y = −8.375 cm). Moreover, vertical and horizontal profiles of the measured and calculated streamwise velocity at x = 112.5 cm are presented in Figure 3b and c respectively. The simulated vertical velocity profile overlaps well with the measured data in the middle of the flume (i.e., y = 0 cm), whereas the calculated velocities near the sidewall (y = −8.38 cm) are approximately 10% smaller than the experimental data, which can be attributed to the wall function effects. As for the horizontal velocity distribution, the simulated velocities are also generally consistent with the measured data, with the bottom values (z = 1 cm) noticeably lower than those in the middle (z = 6 cm) and near the surface (z = 11 cm). While the bottom and middle velocities exhibit typical U-shape horizontal profiles, the surface velocities show a remarkable attenuation from the middle to both sides due to influence of shock waves, which are also reflected in the fluctuant measured data. Nevertheless, the overall accuracy of the simulation is considered acceptable to support the analysis of the flow pattern and energy dissipation characteristics.

4.2. Flow Pattern

Figure 4 shows the side view of the flow pattern downstream of three different deflectors. Lateral jets can be found to form behind the deflectors, and the jets travel further at higher elevations, leading to a non-uniform lateral cavity that is consequently longer and wider at higher elevations. The jet impacting the sidewall is deflected immediately and keeps rising, clinging to the sidewall, associated with noticeable turbulence and aeration. In contrast, the water below the impacting region remains non-aerated. The three deflectors also form different boundary lines between the white and black water regions. With type A, the boundary line is smooth, above which the deflected flow develops into water-wings and continually rises over the sidewalls. With type B, the boundary line is fluctuant, and the deflected flow develops into water-wings that intermittently jump over the sidewall, and a few bubbles could be found in the lower clear water. As for type C, neither distinguished boundary line nor water-wings can be found since the aeration in the lower region is also noticeable. It appears that the deflected flow is submerged in the lateral jet and induces stronger turbulence and intensive aeration since a larger amount of air bubbles can be seen at the impacting region throughout the flow depth. This is probably caused by the rapid close of the lateral cavity and the stronger interaction between the jet and the deflected flow.
The jet trajectories behind three different deflectors are marked with white lines and shown in Figure 5. The detaching jet behind deflector A moves in the direction tangent to the hypotenuse, leading to the widest and largest lateral cavity, which further leaves enough space for the lower deflected flow to develop to the upper region without contacting the jet trajectory. However, the jet behind deflector B exhibits a fluctuating phenomenon and forms an unstable cavity smaller and narrower than that behind deflector A. The fluctuating phenomenon can be mainly ascribed to the continuously varying edge slope of the deflector since it could yield non-uniform velocity distribution that exacerbates jet surface fluctuations. A benefit from the additional straight line in comparison to deflector B, the flow at the tail of deflector C has developed to a more uniform condition before detaching away from the deflector. It exhibits a trajectory almost parallel to the sidewall right behind the deflector, forming the narrowest and smallest lateral cavity. The rapidly closed cavity blocks the rising passage of the underneath deflected flow, preventing the formation of water-wings.
The development of water-wings and the subsequent shock waves downstream of the three deflectors are shown in Figure 6. It can be discovered that deflectors A and B form shock waves in bundle shape and diamond shape, respectively, whereas no shock wave is observed downstream of deflector C. With deflector A, the water-wings from both sides quickly rise up and detach from the main flow, jetting to the middle and colliding with each other in the air, and then fall back to the main flow with a certain vertical velocity and the lateral momentum dissipated. The returned water further produces disturbance on the surface and thereby inducing a bundle-shaped shock wave. With deflector B, the lateral momentum of the water-wings partially counteracts the jet in the cavity, and thus the water-wings fail to contact each other before falling back to the main flow. The returned water still contains certain lateral momentum and moves together with the main flow, producing two repeated reverse oblique developing lines on the flow surface and finally developing into a diamond-shaped shock wave. As for deflector C, the deflected flow is restricted in the relatively small cavity and thus only forms small water crowns, which are restricted in the vicinity of the sidewall and disappear further downstream.
Based on the above observations, it can be concluded that a continuous variation of the lateral deflector surface at the tail with an additional flow guiding extension is the key to the elimination of the water-wings and shock waves.

4.3. Velocity Distribution

The near-bottom (z = 1 cm) velocity distribution around three deflectors is comparatively shown in Figure 7. The jets behind deflectors A and B are firstly contracted and reach the maximum cavity width at approximately 1/3 of the cavity length. After that, the jets restart to spread to the sidewall. In comparison, the jet behind deflector C continuously moves close to the sidewall. Consequently, deflector A exhibits the largest cavity size, followed by deflectors B and C, which is in agreement with the experimental observation. It is worth emphasizing that the velocity vector distribution around deflector B is apparently more non-uniform compared to those around deflectors A and C, which is the reason for the fluctuating phenomenon of the jet surface. Moreover, the maximum velocity behind deflector C is about 0.2 m/s lower than those behind the other two deflectors, and the near-wall low-velocity region (i.e., the purple dashed rectangle in Figure 7) behind deflector C is obviously larger than those behind deflectors A and B. These phenomena all imply more effective energy dissipation of deflector C compared to A and B.
The near-wall (y = −9.3 cm) vertical distribution of velocity behind the deflectors is comparatively shown in Figure 8 to illustrate the kinematic characteristics of the deflected flow. Distinguished by the velocity magnitude and vector direction, the longitudinal flow behind deflectors A and B can be divided into three distinct regions: the main flow region that flows downstream, the water-wing region that rises up while flowing downstream and the impacting region that exhibits the highest velocities and separating the above two regions. The rising water-wing regions are in agreement with the experimentally observed boundary lines. As for deflector C, the flow direction is almost not affected by the jet since only a small region can be found that has upward velocities and therefore is the underlying reason for the fast recovery of the rising surface and the elimination of shock waves. Moreover, the velocities behind deflector C are lower than those behind deflectors A and B, especially dropping about 50% at the impacting region. This again confirms the more effective energy dissipation of deflector C.
Figure 9 shows the streamwise development of the cross-sectional velocity (CSV) with three different and without deflectors, which is computed from 17 flux surfaces with a spacing of 7.5 cm. The no-deflector scheme shows a linearly increased CSV, whereas for the other three schemes, the CSVs increase at the deflector region and decrease at the cavity region because of lateral flow contractions and expansions. With deflectors A and B, the CSVs first increase in the front of the cavity region and then decrease, whereas the CSVs of the flow behind deflector C exhibits a consistent decreasing trend due to continuous flow expansion. Affected by the shock waves, the CSVs of the downstream flow with deflectors A and B exhibit some fluctuations, and the fluctuating amplitude of deflector B is higher than that of deflector A. In contrast, the CSVs of deflector C feature a linear increase behavior similar to that of the no-deflector scheme in this region, although with lower values. Moreover, it can be found that the efficient energy dissipation mainly occurs at the impacting region (i.e., approximately 15 < x < 40 cm).

4.4. Energy Dissipation Characteristics

Figure 10 shows the streamwise development of the flux-averaged hydraulic head (FAHH) along the flume with three different deflectors and without deflectors. In FLOW-3D, the FAHH is recommended to evaluate the hydraulic head where significant changes occur. It is calculated using the following equation:
h = ( u 2 2 g + p ρ g + z ) d s φ d s
here, h is the FAHH, φ is the flux across a cell surface, and ds is the open area of the cell face.
The FAHH decreases linearly in the flume without lateral deflectors. The FAHH of the flow underwent some certain fluctuations around the no-deflector values when passing deflectors A and B. These fluctuations continue until they are twice the cavity length downstream (i.e., to approximately x = 60 cm). Further downstream, the FAHH values of the flow behind deflectors A and B turn out to be lower than the no-deflector values. In contrast, a rapid decrease in FAHH can be found at the deflector and cavity region (0 < x < 20) for the flume equipped with deflector C, and further downstream (x ≥ 45), the values of FAHH decrease steadily with a similar rate as those of the no-deflector scheme. A quantitative comparison of the energy dissipation effect of the three deflectors can be achieved using the local head loss coefficient ζ = ( FAHH 1 FAHH 2 ) 2 g v 2 between x = −17 cm and x = 45 cm with v being a reference CSV calculated at x = −17 cm. The local head loss coefficients for deflectors A, B and C are 2.72%, 2.90% and 4.2%, respectively. This behavior is in agreement with the streamwise development of the CSV and confirms the kinetic energy is mainly dissipated at the deflector and jet region.
To further analyze the energy dissipation mechanism of the lateral deflectors, the distribution of turbulent kinetic energy k and turbulence dissipation rate ε in the deflector and jet region are comparatively shown in Figure 11. The most noticeable difference in Figure 11 is the larger k and ε values at the impacting region behind deflector C compared to those behind deflectors A and B. This implies more mean kinetic energies are turned into turbulent energies and are then dissipated for the flow behind deflector C. This is consistent with the intensive aeration observed in the flow photos shown in Figure 6, as higher turbulent levels induce more intensive aeration [44,45,46].

5. Conclusions

In this paper, the hydraulic characteristics of lateral deflectors with three different geometries in gentle-slope free-surface tunnels were investigated using the hybrid approach of model test and numerical simulation. Special focus was placed on the flow pattern and energy dissipation features. The main findings are:
  • The cavity formed behind lateral deflectors usually features a right-angled trapezoid shape with a larger streamwise length at higher elevations because of non-uniform velocity distributions. This makes the deflected flow rise up along the impacting region inside the cavity and potentially induced shock waves depending on the interaction of the rising up water-wings and the jet surfaces.
  • The traditional triangular deflector forms an adequately wide cavity that allows for the free rising up of the water-wings inside the cavity, which further contributes to the development of the buddle-type shock wave, whereas the two-arc deflector yields a jet with fluctuating surface, resulting from the non-uniform planar velocity distribution caused by the continuously varying curvature of the arcs. Water-wings also develop inside the cavity and eventually produce a diamond-type shock wave downstream. In contrast, the jet behind the two-arc deflector with a straight guiding line at the tail is stabler and travels a shorter distance before impacting the side wall. The jet could thus restrict the development of the rising flow, and thereby eliminate the formation of water-wings and shock waves. Based on these observations, it is concluded that a continuous variation of the lateral deflector surface at the tail with an additional flow guiding extension is the key to the elimination of the water-wings and shock waves.
  • Compared to the triangular deflector and the two-arc deflector, the two-arc deflector with a straight line exhibits more effective energy dissipation, as reflected in the local energy loss coefficient. The underlying reason for its effective energy dissipation is the more intensive turbulence introduced by the stronger interaction between the deflected flow and the jet surface, which also leads to more intensive aeration.
Compared to real-world engineering problems, the present study has some limitations. The physical model is roughly 1/24~1/80 of real-world free-surface tunnels, and thus scale effect is expected, in particular for air–water two-phase flow characteristics [47,48,49]. However, the scale of the physical model is considered acceptable with regards to the main concern of the current paper, which is the comparative evaluation of the different lateral deflectors in terms of the flow pattern and the energy dissipation behaviors. Particularly, the jet trajectory behind the lateral deflector and its interaction with the deflected flow could be decently reproduced in the physical model. Moreover, considering that the numerical simulation approach remains immature for highly aerated flows but much more reliable for non-aerated flow [13,50,51], the hybrid approaches used in this paper could offer valuable insight into the underlying reason for the flow pattern improvement and higher energy dissipation observed with deflector type C. Aeration is another key concern for lateral deflector designs but is not investigated in detail in this paper due to facility reasons. Nonetheless, the findings from the study provide a novel and preferable lateral deflector design for gentle-slope free-surface tunnels, which could, if not resolve, significantly improve the unwanted flow patterns of water-wings and shock waves. Meanwhile, it could also achieve improved energy dissipation compared to traditional alternatives. The aeration characteristics of the novel deflector will be investigated in future research.

Author Contributions

The conceptualization of this research was directed by J.W. and the experimental data were from J.D. and S.D.; the numerical modeling was performed by J.D. under the supervision of J.W. and Z.D.; the validation is carried by J.D. and S.D.; the manuscript was mainly finished by J.D., Z.D. and J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2021YFC3090105), the China Postdoctoral Science Foundation (No. 2021M692754), and Guizhou Science and Technology Joint Support (2019) No. 2890.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

A grid convergence study involving three mesh schemes was carried out to evaluate the grid quality. The cell size (x × y × z) in the refined region of the three mesh schemes was 0.625 cm × 0.313 cm × 0.625 cm (coarse), 0.375 cm × 0.188 cm × 0.375 cm (medium), and 0.188 cm × 0.094 cm × 0.188 cm (fine), respectively.
The grid spacing of the 3D grids is defined as h k = h x , k · h y , k · h z , k 3 . The computational solution f for a grid k is defined as the computed maximum streamwise velocities umax in the middle longitudinal section at x = 112.5 cm. The GCI for each grid k is computed as G C I k = F s | E k | , where a value of 1.25 is chosen for the factor of safety F s , following the recommendations of Roache [52].
Table A1 shows the grid convergence index (GCI) calculated using the approach of Roache and Salas [52]. The observed order of accuracy is p = 1.84, and the GCI shows a convergence trend. The fine grid takes a very satisfactory value of 0.224%, which is an order of magnitude higher than the former and therefore is considered fine enough to resolve the flow field.
Table A1. Parameter of grid convergence calculation.
Table A1. Parameter of grid convergence calculation.
Mesh h /cm rp(umax)f(umax)/(m·s−1)ε(umax)/%GCI(umax)/%
coarse0.496-1.842.607--
medium0.2971.672.6330.9972.040
fine0.1491.992.6210.4560.224
The vertical distributions of streamwise velocity ux and turbulent kinetic energy k in the middle longitudinal section at x = 112.5 cm are comparatively shown in Figure A1. From Figure A1, it can be found that the main deviations of the ux profiles lie in the region close to the water surface and the bottom wall, whereas in the middle, all three meshes return almost identical velocity values. Moreover, the difference between the medium and fine meshes is less obvious compared to that between the coarse and the medium meshes. As for the k profile, the result from the coarse mesh is strikingly different from those from the medium and fine meshes, the difference between which is relatively negligible. Therefore, the fine grid was selected for the simulation.
Figure A1. Vertical distributions of streamwise velocity ux and turbulent kinetic energy k in the middle longitudinal section at x = 112.5 cm.
Figure A1. Vertical distributions of streamwise velocity ux and turbulent kinetic energy k in the middle longitudinal section at x = 112.5 cm.
Water 14 02689 g0a1

References

  1. Xie, Z.S.; Wu, Y.H.; Chen, W.X. New technology and innovation on flood discharge and energy dissipation of high dams in China. J. Hydraul. Eng. 2016, 47, 324–336. (In Chinese) [Google Scholar] [CrossRef]
  2. Ma, H.; Chi, F. Major technologies for safe construction of high earth-rockfill dams. Engineering 2016, 2, 498–509. [Google Scholar] [CrossRef]
  3. Ma, H.; Chi, F. Technical progress on researches for the safety of high concrete-faced rockfill dams. Engineering 2016, 2, 332–339. [Google Scholar] [CrossRef]
  4. Fan, L.; Zhang, H.W.; Liu, Z.P.; Zhang, D. Numerical study on hydraulic characteristic of free surface flow in spillway tunnel with different configuration. J. Hydroelectr. Eng. 2009, 28, 126–131. (In Chinese) [Google Scholar]
  5. Guo, J.; Zhang, D.; Liu, Z.P.; Fan, L. Achievements on hydraulic problems in large spillway tunnel with a high head and large discharge flow and its risk analysis. J. Hydraul. Eng. 2006, 37, 1193–1198. (In Chinese) [Google Scholar] [CrossRef]
  6. Pang, C.J.; Yuan, Y.Z. Research on type of aerator in large gentle-slope free-surface ‘dragon-head-up’ spillway tunnel. J. Hydraul. Eng. 1993, 6, 61–66. (In Chinese) [Google Scholar] [CrossRef]
  7. Bai, R.; Liu, S.; Wang, W.; Zhang, F. Experimental investigation of the dissipation rate in a chute aerator flow. Exp. Therm. Fluid Sci. 2019, 101, 201–208. [Google Scholar] [CrossRef]
  8. Suo, L.S.; Liu, N. Handbook of Hydraulic Structure Design, 2nd ed.; China Waterpower Press: Beijing, China, 2014. (In Chinese) [Google Scholar]
  9. Chanson, H. Study of air entrainment and aeration devices. J. Hydraul. Res. 1989, 27, 301–319. [Google Scholar] [CrossRef]
  10. Du, L.Q.; Tang, Y.Q.; Wang, J.S. Construction technology of slip form for bottom plate of turning point section of spillway tunnel of Wudongde Hydropower Station. Yangtze River 2020, 51, 193–196. (In Chinese) [Google Scholar] [CrossRef]
  11. Suo, H.M.; Zhou, Z.J.; He, L.; Zhang, J.M. Study on type of tunnel spillway in Shuangjiangkou Hydropower Station. Des. Hydroelectr. Power Stn. 2021, 37, 9–12. (In Chinese) [Google Scholar] [CrossRef]
  12. Wang, H.Y.; Dai, G.Q.; Yang, Q.; Liu, C. Experimental study on V-type aerator for spillway tunnel with inlet raised. J. Hydraul. Eng. 2005, 36, 102–105. (In Chinese) [Google Scholar] [CrossRef]
  13. Yang, J.; Teng, P.; Zhang, H. Experiments and cfd modeling of high-velocity two-phase flows in a large chute aerator facility. Eng. Appl. Comp. Fluid Mech. 2019, 13, 48–66. [Google Scholar] [CrossRef]
  14. Hager, W.H.; Boes, R.M. Hydraulic structures: A positive outlook into the future. J. Hydraul. Res. 2014, 52, 299–310. [Google Scholar] [CrossRef]
  15. Li, S.; Zhang, J.; Xu, W.; Chen, J.; Peng, Y.; Li, J.; He, X. Simulation and experiments of aerated flow in curve-connective tunnel with high head and large discharge. Int. J. Civ. Eng. 2016, 14, 23–33. [Google Scholar] [CrossRef]
  16. Juon, R.; Hager, W.H. Flip bucket without and with deflectors. J. Hydraul. Eng. 2000, 126, 837–845. [Google Scholar] [CrossRef]
  17. Wang, H.Y.; Yang, Q.; Dai, G.Q.; Liu, C. Study on hydraulic characteristic of aerator with side-walls sudden enlargement in free-flow spillway tunnel. Adv. Mater. Res. 2012, 614, 613–616. [Google Scholar] [CrossRef]
  18. Ma, X.; Dai, G.; Yang, Q.; Li, G.; Zhao, L. Analysis of influence factors of cavity length in the spillway tunnel downstream of middle gate chamber outlet with sudden lateral enlargement and vertical drop aerator. J. Hydrodyn. 2010, 22, 657–663. [Google Scholar] [CrossRef]
  19. Pfister, M.; Hager, W.H. Chute aerators. II: Hydraulic design. J. Hydraul. Eng. 2010, 136, 360–367. [Google Scholar] [CrossRef]
  20. Nie, M.; Wang, X.; Wu, G. Effect of lateral deflector on outlet cavity lengths. J. Hydraul. Eng. 2003, 129, 536–540. [Google Scholar] [CrossRef]
  21. Xu, J.; Wu, J.; Peng, Y.; Ma, F. Fin performance of 3-d aerator devices with backward lateral deflectors. J. Hydrodyn. 2020, 32, 410–413. [Google Scholar] [CrossRef]
  22. Siying, W.; Jinbo, Y.; Caihuan, W. Study on the aerator of spillway tunnel in one large-scale project. In Proceedings of the 35th World Congress of the International-Association-for-Hydro-Environment-Engineering-and-Research (IAHR), Chengdu, China, 8–13 September 2013. [Google Scholar]
  23. Ye, F.; Xu, W.; Wei, W. Experimental study of the air concentration diffusion in aerated chute flows downstream of lateral and bottom aerators. AIP Adv. 2022, 12, 25221. [Google Scholar] [CrossRef]
  24. Nie, M.X.; Wu, G.G. Effect of the lateral deflector on cavity design. J. Tsinghua Univ. 2002, 42, 546–550. (In Chinese) [Google Scholar] [CrossRef]
  25. Nie, M. Comparison of energy dissipation with and without aerators. Tsinghua Sci. Technol. 2000, 5, 222–225. [Google Scholar]
  26. Lucas, J.; Hager, W.H.; Boes, R.M. Deflector effect on chute flow. J. Hydraul. Eng. 2013, 139, 444–449. [Google Scholar] [CrossRef]
  27. Nie, M.X.; Duan, B.; Li, L.L. Hydraulic characteristics along the sidewall of a sluice downstream of a sudden lateral enlargement and vertical drop. J. Tsinghua Univ. 2006, 46, 1969–1972. (In Chinese) [Google Scholar] [CrossRef]
  28. Liu, C.; Zhang, G.K.; Li, N.W.; Zhang, G.B.; Wang, H.Y. The effect of base aerator on the length of lateral cavity. Adv. Eng. Sci. 2008, 40, 1–4. (In Chinese) [Google Scholar]
  29. Wu, G.L.; Ma, F.; Wu, J.H.; Xu, J.R.; Peng, Y.; Wang, Z. Experimental investigation of flow regime control for dragon-drop-tail section of the spillway tunnel. Chin. J. Hydrodyn. 2015, 30, 140–145. (In Chinese) [Google Scholar] [CrossRef]
  30. Wu, J.; Li, D.; Ma, F.; Qian, S. Fin characteristics of aerator devices with lateral deflectors. J. Hydrodyn. 2013, 25, 258–263. [Google Scholar] [CrossRef]
  31. Zhang, H.W.; Liu, Z.P.; Zhang, D.; Wu, Y.H. Study on hydraulic characteristics of sidewall aerators in high-head large-discharge spillway tunnel. J. Hydroelectr. Eng. 2015, 34, 111–116. (In Chinese) [Google Scholar] [CrossRef]
  32. Li, F.T.; Liu, P.Q.; Xu, W.L.; Tian, Z. Experimental study on effect of flaring piers on weir discharge capacity in high arch dam. J. Hydraul. Eng. 2003, 11, 43–47. (In Chinese) [Google Scholar] [CrossRef]
  33. Zheng, X.; Zhang, L.; Wu, S.; Song, K. Study on the shape of the aerator of high-head discharge tunnel with mild bottom slope. Water-Sui. 2021, 13, 2128. [Google Scholar] [CrossRef]
  34. Li, S.; Zhang, J.; Xu, W.; Chen, J.; Peng, Y. Evolution of pressure and cavitation on side walls affected by lateral divergence angle and opening of radial gate. J. Hydraul. Eng. 2016, 142, 05016003. [Google Scholar] [CrossRef]
  35. Yang, B.; Wang, Z.; Wu, M. Numerical simulation of naval ship’s roll damping based on cfd. Procedia Eng. 2012, 37, 287–293. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, H.; Zhang, L.; Wu, S.; Wang, F.; Zhan, Z.; Zheng, X.; Zhang, H.; Bao, W. Research on flood discharge and energy dissipation of a tunnel group layout for a super-high rockfill dam in a high-altitude region. Water 2021, 13, 3408. [Google Scholar] [CrossRef]
  37. Zhang, W.; Wang, J.; Zhou, C.; Dong, Z.; Zhou, Z. Numerical simulation of hydraulic characteristics in a vortex drop shaft. Water 2018, 10, 1393. [Google Scholar] [CrossRef]
  38. Pfister, M.; Hager, W.H. Chute aerators. I: Air transport characteristics. J. Hydraul. Eng. 2010, 136, 352–359. [Google Scholar] [CrossRef]
  39. Feurich, R.; Rutschmann, P. Efficient design optimization and investigation by combining numerical and physical models. In Proceedings of the Korea Water Resources Association Conference, Korea Water Resources Association, Seoul, Korea, 11–16 September 2005. [Google Scholar]
  40. Johnson, M.C.; Savage, B.M. Physical and numerical comparison of flow over ogee spillway in the presence of tailwater. J. Hydraul. Eng. 2006, 132, 1353–1357. [Google Scholar] [CrossRef]
  41. Stamou, A.I.; Chapsas, D.G.; Christodoulou, G.C. 3-d numerical modeling of supercritical flow in gradual expansions. J. Hydraul. Res. 2008, 46, 402–409. [Google Scholar] [CrossRef]
  42. Hirt, C.W.; Nichols, B.D. Volume of fluid (vof) method for the dynamics of free boundaries. J. Comput. Phys. 1981, 39, 201–225. [Google Scholar] [CrossRef]
  43. Yakhot, V.; Orszag, S.A.; Thangam, S.; Gatski, T.B.; Spezialea, C.G. Development of turbulence models for shear flows by a double expansion technique. Phys. Fluids A Fluid Dyn. 1992, 4, 1510–1520. [Google Scholar] [CrossRef]
  44. Dong, Z.; Bürgler, M.; Hohermuth, B.; Vetsch, D.F. Density-based turbulence damping at large-scale interface for reynolds-averaged two-fluid models. Chem. Eng. Sci. 2022, 247, 116975. [Google Scholar] [CrossRef]
  45. Hirt, C.W. Modeling Turbulent Entrainment of Air at a Free Surface; Technical Note FSI-01-12; Flow Science, Inc.: Santa Fe, NM, USA, 2003. [Google Scholar]
  46. Dong, Z.; Wang, J.; Vetsch, D.F.; Boes, R.M.; Tan, G. Numerical simulation of air–water two-phase flow on stepped spillways behind x-shaped flaring gate piers under very high unit discharge. Water 2019, 11, 1956. [Google Scholar] [CrossRef]
  47. Pfister, M.; Hager, W.H. Deflector-jets affected by pre-aerated approach flow. J. Hydraul. Res. 2012, 50, 181–191. [Google Scholar] [CrossRef] [Green Version]
  48. Pfister, M.; Lucas, J.; Hager, W.H. Chute aerators: Preaerated approach flow. J. Hydraul. Eng. 2011, 137, 1452–1461. [Google Scholar] [CrossRef]
  49. Pfister, M.; Chanson, H. Scale effects in physical hydraulic engineering models. J. Hydraul. Res. 2011, 49, 293–306. [Google Scholar] [CrossRef]
  50. Hohermuth, B.; Schmocker, L.; Boes, R.M.; Vetsch, D.F. Numerical simulation of air entrainment in uniform chute flow. J. Hydraul. Res. 2021, 59, 378–391. [Google Scholar] [CrossRef]
  51. Hohermuth, B.F. Aeration and Two-Phase Flow Characteristics of Low-Level Outlets. Ph.D. Thesis, ETH Zurich, Zürich, Switzerland, 2019. [Google Scholar]
  52. Roache, P.J. Quantification of uncertainty in computational fluid dynamics. Annu. Rev. Fluid Mech. 1997, 29, 123–160. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Schematic diagram of the experimental spillway tunnel and details of (b) the simulated region and (c) the three lateral deflectors (unit is cm).
Figure 1. (a) Schematic diagram of the experimental spillway tunnel and details of (b) the simulated region and (c) the three lateral deflectors (unit is cm).
Water 14 02689 g001
Figure 2. Schematic diagram of the simulation setup.
Figure 2. Schematic diagram of the simulation setup.
Water 14 02689 g002
Figure 3. Comparison of the experimental and numerical results: (a) water surface profile; (b) vertical; and (c) horizontal distribution of the streamwise velocity ux at x = 112.5 cm. The unit of the coordinates is cm.
Figure 3. Comparison of the experimental and numerical results: (a) water surface profile; (b) vertical; and (c) horizontal distribution of the streamwise velocity ux at x = 112.5 cm. The unit of the coordinates is cm.
Water 14 02689 g003
Figure 4. Side view of the flow pattern downstream of the deflectors.
Figure 4. Side view of the flow pattern downstream of the deflectors.
Water 14 02689 g004
Figure 5. Top view of the flow pattern near the lateral cavity.
Figure 5. Top view of the flow pattern near the lateral cavity.
Water 14 02689 g005
Figure 6. Top view of the flow pattern downstream of the deflectors.
Figure 6. Top view of the flow pattern downstream of the deflectors.
Water 14 02689 g006
Figure 7. Horizontal velocity distribution around the deflectors at z = 1.0 cm. The vectors are generated using the same density and scale for the three types of deflectors.
Figure 7. Horizontal velocity distribution around the deflectors at z = 1.0 cm. The vectors are generated using the same density and scale for the three types of deflectors.
Water 14 02689 g007
Figure 8. Velocity distribution near the sidewall (y = −9.3 cm). The vectors are generated using the same density and scale for the three types of deflectors.
Figure 8. Velocity distribution near the sidewall (y = −9.3 cm). The vectors are generated using the same density and scale for the three types of deflectors.
Water 14 02689 g008
Figure 9. Streamwise development of the cross-sectional velocity along the flume with three different deflectors and without deflectors. The brown color indicates regions where the cavity range of C (pink) overlaps those of A and B (green).
Figure 9. Streamwise development of the cross-sectional velocity along the flume with three different deflectors and without deflectors. The brown color indicates regions where the cavity range of C (pink) overlaps those of A and B (green).
Water 14 02689 g009
Figure 10. Streamwise development of the flux-averaged hydraulic head (FAHH) along the flume with three different and without deflectors. The brown color indicates regions where the cavity range of C (pink) overlaps those of A and B (green).
Figure 10. Streamwise development of the flux-averaged hydraulic head (FAHH) along the flume with three different and without deflectors. The brown color indicates regions where the cavity range of C (pink) overlaps those of A and B (green).
Water 14 02689 g010
Figure 11. Horizontal distribution of (a) turbulent kinetic energy k; (b) turbulence dissipation rate ε in deflector and cavity region at z = 1 cm.
Figure 11. Horizontal distribution of (a) turbulent kinetic energy k; (b) turbulence dissipation rate ε in deflector and cavity region at z = 1 cm.
Water 14 02689 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Da, J.; Wang, J.; Dong, Z.; Du, S. Hydraulic Characteristics of Lateral Deflectors with Different Geometries in Gentle-Slope Free-Surface Tunnels. Water 2022, 14, 2689. https://doi.org/10.3390/w14172689

AMA Style

Da J, Wang J, Dong Z, Du S. Hydraulic Characteristics of Lateral Deflectors with Different Geometries in Gentle-Slope Free-Surface Tunnels. Water. 2022; 14(17):2689. https://doi.org/10.3390/w14172689

Chicago/Turabian Style

Da, Jinrong, Junxing Wang, Zongshi Dong, and Shuaiqun Du. 2022. "Hydraulic Characteristics of Lateral Deflectors with Different Geometries in Gentle-Slope Free-Surface Tunnels" Water 14, no. 17: 2689. https://doi.org/10.3390/w14172689

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