Next Article in Journal
Droplet Evaporation in a Gas-Droplet Mist Dilute Turbulent Flow behind a Backward-Facing Step
Next Article in Special Issue
Do the Volume-of-Fluid and the Two-Phase Euler Compete for Modeling a Spillway Aerator?
Previous Article in Journal
A Monitoring System of Sand Mining in Large Rivers and Its Application to the Ayeyarwady (Irrawaddy) River, Myanmar
Previous Article in Special Issue
Hydrodynamics and Free-Flow Characteristics of Piano Key Weirs with Different Plan Shapes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Apron and Cutoff Wall Scour Protection for Piano Key Weirs

1
Utah Water Research Laboratory, Department Civil and Environmental Engineering, Utah State University, Logan, UT 84321, USA
2
DESTEC-Department of Energy, Systems, Territory and Construction Engineering, University of Pisa, 56122 Pisa, Italy
*
Author to whom correspondence should be addressed.
Water 2021, 13(17), 2332; https://doi.org/10.3390/w13172332
Submission received: 16 June 2021 / Revised: 17 August 2021 / Accepted: 20 August 2021 / Published: 26 August 2021
(This article belongs to the Special Issue Advances in Spillway Hydraulics: From Theory to Practice)

Abstract

:
Piano key (PK) weirs are used in a variety of flow control structure applications, including spillway crests and open channel diversion structures. However, to the best of authors’ knowledge, structure-specific design guidance for scour mitigation is still needed. To fill this gap of knowledge, a systematic experimental campaign was conducted by testing different configurations of horizontal aprons with a cutoff wall. Protection structures were located at the toe of the PK weir. Namely, experiments were performed at large-scale to assess the effect of three apron lengths on downstream scour hole geometry under different hydraulic conditions. It was observed that a horizontal apron deflects the plunging jets originating from the PK weir, thus significantly reducing scour. Experimental evidence allowed corroboration that significant scour depth reduction occurs for an apron length 1.5 times the weir height, with longer aprons found to provide marginal benefits. Finally, also provided herein are tools to estimate the main scour characteristics and help practitioners in optimizing apron design.

1. Introduction

When considering current design standards, sustainability principles, and observed trends in flooding, the control of local scour represents a challenging task. Therefore, the analysis of the interaction between catchment-based processes (e.g., runoff hydrographs, river morphology, and transport of sediments, and debris) and hydraulic structure functioning is of fundamental importance to minimize the risk of failure [1,2,3,4,5]. Weirs are a common hydraulic structure found in rivers and incorporated in dams and levees. Often, localized scour phenomena are caused by falling nappe or oblique plunging jets originating from weirs. Depending on jet characteristics (e.g., submergence, geometry, and inclination), scour morphology at equilibrium could be either two- or three-dimensional (considering the vertical and streamwise directions with the third dimension being transverse) [4,5,6,7,8,9,10].
Despite the numerous studies dealing with drop structures, including linear weirs, the scour process downstream of piano key (PK) weirs has been only recently analyzed and results are still insufficient to correctly assess such a complex phenomenon. Although design manuals, e.g., [11,12,13,14,15,16,17] can represent a valid help for practitioners with upstream and downstream erosion at PK weirs discussed in [18], none of them provide guidance for PK weirs’ scour-induced processes and selection of a horizontal apron and cuttoff wall for scour protection.
Nevertheless, PK weirs have become quite popular. They have been recently employed in a number of dam rehabilitation projects in Asia, Australia, Europe, and North America, as well as in some river restoration projects [19,20,21] (e.g., the run-of-river Van Phong PK weir in Vietnam). It was found that PK weirs perform well as in-channel flow control structures due to their hydraulic efficiency, do not require operation (passive flow control), possess construction economy, and have the ability to pass floating debris [22,23,24,25]. Following their recent success in practical applications, there has been increased interest regarding scour morphology at PK weirs [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31].
However, to date there are very few studies providing specific tools to estimate scour characteristics and, apparently, only one dealing with scour mitigation. More specifically, from scaled hydraulic models of structures in Vietnam, [19] suggested that PK and labyrinth weirs in rivers with high submergence only need a short stilling basin and end sill to dissipate energy. The authors of [27] investigated the scour formation downstream of PK weirs in a 0.67 m-wide channel (i.e., relatively narrow channel), providing some interesting insights on the erosion mechanism and useful empirical equations to predict the main geometric characteristics of the equilibrium morphology. Following this study, [28] focused on the design of a pre-excavated rip-rap apron. However, they did not provide comprehensive guidelines to design effective countermeasures to stabilize the stilling basin downstream of PK weir.
To fill this gap of knowledge, this study analyzed the scour process occurring downstream of a PK weir in the presence of downstream apron. This protection structure was selected because it is a simple and cost-effective alternative for scour mitigation. In addition, cutoff walls are commonly included with concrete aprons to prevent a failure mode via undermining of the slabs [32,33,34,35]. To this end, tests were conducted varying both the length of the apron and the height of the cutoff wall under different hydraulic conditions and with two bed materials. A dedicated, large-scale model was built at Utah State University. Experimental evidence allowed the assessment of the effect of protection structure geometry on the induced scour mechanism. Furthermore, this study provides unprecedented design guidelines that can represent a valid tool for practitioners.

2. Materials and Methods

A rectangular flume (16-m long, 2-m wide, and 1.8-m deep) was constructed at the Utah Water Research Laboratory. A transparent acrylic wall section allowed monitoring of the scour evolution. The flume provided uniform flows in the headbox to a Type A [36] PK weir via a diffuser pipe and rock baffling. A four-key acrylic PK weir was located in the flume whose dimensions are summarized in Table 1. The movable bed section downstream of the PK weir was 5 m long, 2 m wide, and 1.09 m deep. Tailwater was controlled with a stop log structure located 9.68 m downstream of the PK weir. Figure 1 shows diagram sketches of the experimental apparatus along with the main geometric parameters, where W = total PK weir width, B = depth of the PK weir, Bb = streamwise length of the base of the PK weir, Bi = the length of the inlet key, Bo = the length of the outlet key, So = slope of the outlet key, Si = slope of the inlet key, Wu = cycle width, N = number of keys or cycles, P = weir height, Ts = sidewall thickness, Wi = width of the inlet key, Wo = width of the outlet key, Pd = height of the PK weir base relative to the invert of the channel, xM1, xM2, and xM3 are streamwise lengths from the apron to sphere measurement locations, La is the length of the apron, and Lc is the length or depth of the cutoff wall. Note that a false floor was not included upstream of the PK weir, as its presence would be negligible on the results of this study. We also note that the dimensions of the apron and cutoff wall are not indicative of actual thicknesses as may be structurally designed for specific applications.
Prior to each test, the movable bed was prepared to be level (zero bed slope) and at the elevation of the base of the PK weir. The various apron lengths and cutoff walls placed in the moveable bed were fabricated using dimensional lumber and coated for waterproofing. When the aprons were installed, they were also leveled at the elevation of the base of the PK weir, matching the bed surface.
Experimental tests involved two non-cohesive gravel substrates, i.e., Substrate 1 and 2. Substrate 1 was a coarse gravel with d90 = 20.00 mm, d50 = 13.00 mm, uniformity coefficient σ = 1.54, density ρs = 2604 kg/m3, and relative density G = ρs/ρ = 2.61, where dxx is the diameter for which xx% is finer and ρ is the water density. Substrate 2 was a fine gravel with d90 = 9.10 mm, d50 = 6.50 mm, σ = 1.30, ρs = 2604 kg/m3, and G = 2.61. A diagram sketch illustrating the main hydraulic and scour parameters is shown in Figure 2, with H indicating total hydraulic head (H = hu + V2/2g) where hu is the depth of flow relative to the weir crest (with total weir crest length = L), g is the standard acceleration due to gravity, and V is an average upstream flow velocity. Figure 2 also shows the PK weir foundation height Pd, the change in energy head, ΔH, hd the downstream flow depth, La the length of apron in the streamwise direction, Lmax is the maximum scour length, Xmax is the streamwise distance from the end of the apron to the maximum scour depth Zmax, and Zmax_ref is the maximum scour depth without an apron, referred to herein as the reference maximum scour depth.
A calibrated venturi meter (±0.25%) was used to measure flow discharge (Q). A point gauge (±0.75 mm) with a stilling well (i.e., point gauge with piezometer) was hydraulically connected two meters upstream of the weir (4.77P) to estimate the upstream piezometric head hu. Conversely, the downstream water level (and, consequently, the tailwater depth hd) was measured by using a Microsonic mic+130 ultra-sonic sensor (±1%) (Microsonic, Dortmund, Germany) [37,38] at 3.35 m downstream (x/P = 8). Initial investigations profiled the downstream water surface to select an appropriate location for computing hd.
Scour evolution and at equilibrium was monitored by using two techniques. Namely, low-buoyancy spheres (G = 0.53) (see [30] for more details) were spaced vertically in the substrate at 0.1 m in vertical columns (measurement accuracy is ±12 mm). Sphere column location was color coded, and the specific weight of the spheres was adjusted so that they would stay buried until approximately half of the sphere was exposed. Once the sphere was uncovered to this point, it quickly rose to the water surface. A GoPro Hero 7 Black high-resolution (4k) camera (GoPro, San Mateo, CA, USA) located above the channel allowed recording this process and estimate the time (t) at which each sphere was released from the bed material, i.e., the time when scour reached the depth at which the sphere was positioned. Overall, six monitoring locations were selected for the spheres and are indicated with numerical identifiers in Figure 1.
To determine equilibrium morphology, the substrate was scanned with an Intel® RealSenseTM D435 depth camera (accuracy of ±1 mm) (Intel, Santa Clara, CA, USA) [30,39,40] processed with a Utah State University (USU) custom MATLAB (vR2020) script. The camera data was cross-checked with point measurements taken on a grid layout with a point gage assembly (accuracy of ±1 mm). This camera is capable of capturing the surface of solids or roughened or opaque water surfaces in three dimensions; its specifications are summarized in Table 2. Post-processing included 2D and 3D scour morphology plots and allowed quantifying maximum scour location (Xmax) and its depth (Zmax), maximum length of the scour hole (Lmax), and other scour features.
Four apron lengths (La) were tested, i.e., La = 0P (reference configuration), 1P, 1.5P, and 2P, under three different combinations of Q and hd (Table 3). Note that all tests (see Table 4) were conducted with the lowest permissible tailwater depth in order to investigate the most critical hydraulic conditions. Each test began by slowly filling the headbox and substrate box until the weir was in a submerged condition. The target Q was set, the hd was lowered to the lowest permissible depth, and then the timer and video recording array were initiated. Note that for certain values of Q the horizontal bed prevented the target hd (~0.33P) from being achieved. Multiple experiments were conducted for more than 18 h each to confirm equilibrium conditions and necessary durations for the remainder of the test matrix. More specifically, to assess the time needed to reach the equilibrium configuration, this study conducted some selected longer tests. In particular, for Substrate 2 and Q = 600 L/s, the reference test 15 was conducted (lasting up 1170 min, i.e., 19.5 h, see Table 4). It was observed that the differences between the equilibrium morphologies at t = 6 h (i.e., the same time duration of Test 18, conducted under similar hydraulic conditions) and at 19.5 h were negligible. Thus, the duration of test 18 was long enough to reach the equilibrium configuration of the bed morphology. Similar considerations apply for all the tested discharges (e.g., for Q = 300 L/s test 14 was conducted with duration t = 1050 min).
Finally, it is worth mentioning that all the tests pertaining to Q ≤ 300 L/s were conducted under clear water conditions [41]. More specifically, the flow intensity U/Uc was less than 1, where U and Uc are the average flow velocity and the average threshold velocity, respectively, with Uc calculated using the methodology developed by [42]. As for tests conducted with Q = 600 L/s and Substrate 1, tests were characterized by U/Uc ≈ 1; Conversely, for tests with Q = 600 L/s and Substrate 2 (i.e., tests 15, 18, 21, and 24), U/Uc was slightly bigger than 1. Nevertheless, it was decided to also include these tests in these analyses, as their behavior in terms of the effect of apron geometry on scour characteristics was found to be essentially consistent with those of other tests. Therefore, they can provide useful information for practical applications.

3. Results

3.1. PK Weir Hydraulics

The PK weir geometry creates a 3D flow field that can be characterized by near-vertical and oblique plunging jets exiting the inlet and outlet keys, respectively. The oblique jet had a larger unit discharge than the near-vertical jet, resulting in greater shear stress acting on the granular bed in correspondence with outlet keys. In the absence of any protection measures, both jets impinge on the gravel bed [10,30]. Because of the excess of shear stress, scour evolution is initially very rapid (developing phase); successively, the kinetics of the scour evolution reduces and scour depth asymptotically reaches its maximum value. More specifically, [4] showed that, for jet-driven scour processes, the rate of material removal scales with the excess of shear stress (τ/τc − 1)m, with τ indicating the shear stress, τc the critical shear stress, and m is a coefficient that is equal to 1.5 for non-cohesive materials according [43]. Therefore, the equilibrium configuration is achieved when τ = τc and, consequently, Z = Zmax, with hd being constant. The presence of a horizontal apron significantly modifies flow features in the downstream stilling basin. In particular, impinging jets are deflected horizontally and a more uniform distribution of the flow occurs downstream of the apron. In addition, a dissipative process may occur on the apron due to the formation of a submerged hydraulic jump. As a consequence, scour takes place downstream of the apron and is characterized by a significant reduction of the main lengths (Figure 3) in comparison with the reference case. However, depending on the tested conditions, the presence of scour may endanger the stability of the apron. Therefore, a cutoff wall is suggested to avoid undermining of the apron and a potential failure mode.

3.2. Analysis of Scour Processes

The analysis of the scour process was conducted by comparing the longitudinal scour profiles for tests with apron with that pertaining to the corresponding reference test (no apron). For the reference tests (e.g., Figure 3) some armoring occurred at equilibrium, but this was dependent upon σ and the angle of repose of the downstream slope of the scour hole, as not all smaller material was removed to result with armoring of only the largest grain sizes (e.g., d90). The overall observed trend was that scour depth and maximum scour dimensions significantly decrease with La under identical hydraulic conditions (Figure 4 and Table 4) (note that for Run 10 no scour occurred in the stilling basin). Results evidenced that the effect of an apron is always prominent, regardless the bed material characteristics (i.e., for both Substrates 1 and 2). As mentioned in the previous section, this behavior can be explained when considering that the flow dynamics are radically different for protected basins. Both flow features and sediment transport mode exhibit similarities with that occurring downstream of an apron caused by submerged horizontal jets. Namely, scour starts in correspondence with the edge of the structure because of the excess of shear stress. During the developing phase, bed material is mobilized and kept in suspension. Successively, when the scour hole enlarges (developed phase), the jet exiting from the apron is deflected downward. Consequently, the sediment transport mode changes, involving both suspended and bed loads. Furthermore, this study corroborates the findings of [35], which showed that the maximum scour depth is a monotonic decreasing function of the non-dimensional apron length and average diameter of the bed material. To this end are presented Figure 4a,c,e (Substrate 1) and Figure 4b,d,f (Substrate 2) where longitudinal scour profiles pertaining to tests are shown. Please note that the apron length and cuttoff wall for each configuration are color matched to the corresponding maximum scour profiles. For example, in Figure 4, A La/P = 1.5 is in red, but is overlapped by La/P = 1.0 that is blue. The red scour profile (dots with dashed line) corresponds to La/P = 1.5.
The percent reduction of the maximum scour depth was quantified as Z*max = 100(Zmax_refZmax)/Zmax_ref, which is caused by the different apron configurations with Zmax_ref indicating the maximum scour depth of the corresponding reference test (i.e., without an apron). Namely, tests were grouped according to the two substrate types and Z*max was plotted as function of the non-dimensional apron length La/P for different densimetric Froude numbers Frd90 = q/{hu[(G − 1)gd90]1/2} (Figure 5a,b). Note that q = Q/W and is the average unit discharge. A slight difference, depending on the inflow characteristics, can be pointed out between tests involving Substrate 1 and 2. Namely, for both the substrates Z*max is a monotonic increasing function of La/P. However, for tests conducted with Substrate 1 and Frd90 = 3.05 and 3.25, Z*max becomes almost constant for 1.5 ≤ La/P ≤ 2, i.e., it only depends upon Frd90 and the effect of La/P becomes negligible (Figure 5a). In addition, the presence of the apron prevents the scour formation for La/P = 2 and Frd90 = 2.44 (i.e., Z*max = 100%, as shown in Figure 5a). A similar behavior can be pointed out for tests with Substrate 2. Overall, Z*max increases with La/P, whereas it slightly decreases with Frd90 (Figure 5b). These results corroborate the findings of [35,44,45], who found that scour processes occurring in protected stilling basins are essentially influenced by the geometry of the protection structures and the inflow conditions. The mentioned analysis allowed the derivation of the following predicting equation, valid for 0 ≤ La/P ≤ 2:
Z m a x * = ( 2.8 F r d 90 6.3 ) ( L a P ) 2 + ( 2.8 F r d 90 + 64.9 ) ( L a P )
Figure 5c–f shows that Equation (1) provides an accurate estimation of the variable Z*max and represents a simple tool that practitioners can use to estimate scour characteristics in protected basins.
Furthermore, in this study the effect of the variable La was quantified regarding its effect on the maximum scour depth Zmax and length Xmax, along with volume of sediment V removed under different inflow conditions and for both the substrates. In particular, Figure 6 shows that, by adding a 1.0P = La horizontal apron, the average reduction of the scour depth is approximately equal to 57%. This percentage increases with La and ranges between 75% and 83% for La = 1.5P and 2.0P, respectively (see also Table 5). It is worth remarking that the change in scour depth for La = 1.5P and La = 2.0P was on average 8%. It implies that, for practical applications and Frd90 > 3, an additional protection (i.e., longer apron) may only minimize scour depth marginally. Likewise, the variable Lmax is also affected by La (Figure 7). However, in this case a different behavior can be pointed out depending on the sediment bed gradation and, in particular, on d50. Namely, for Substrate 1, a general reduction of the scour hole length occurs by increasing La, regardless of the inflow conditions. Conversely, for finer bed material, the maximum scour depth seems to be less affected by the apron length, especially for higher discharges. This behavior is similar to that observed by [45] downstream of block ramps, especially for low slopes of the ramp, for which the horizontal component of the total hydrodynamic force at the ramp toe is much larger than the vertical counterpart. Namely, in the presence of a horizontal, smooth apron, depending on the inflow conditions (i.e., especially for higher values of Frd90), a hydraulic jump may not occur. Thus, the scour downstream of the apron is essentially due to a horizontal, not submerged, jet. The shear stress due to jet generally exceeds the critical counterpart, resulting in a scour formation downstream of the apron. Nevertheless, for finer bed materials, the superficial layer can be more easily transported downstream, thus modifying the overall shape of the scour hole, which becomes more longitudinally extended (Figure 4b,d,f and Figure 7).
Finally, it does not escape attention that Equation (1) also represents a tool to estimate the maximum scour depth for protected basins. Namely, considering the definition of Z*max, Equation (1) can be re-written as follows:
100 ( Z m a x _ r e f Z m a x ) Z m a x _ r e f = ( 2.8 F r d 90 6.3 ) ( L a P ) 2 + ( 2.8 F r d 90 + 64.9 ) ( L a P )
Or, equivalently:
Z m a x = Z m a x _ r e f { 1 0.01 [ ( 2.8 F r d 90 6.3 ) ( L a P ) 2 + ( 2.8 F r d 90 + 64.9 ) ( L a P ) ] }
valid for 2.4 ≤ Frd90 ≤ 5.3, 0 ≤ La/P ≤ 2, 0.11 ≤ H/P ≤ 0.36, and 0.24 ≤ hd/P ≤ 0.59, but the results of this study trend to hd/P ≤1.03 (exceeds tested tailwater depth ratios presented herein). Note that the maximum scour depth in the corresponding reference test (Zmax_ref) can be estimated using the following Equation (4) proposed by [27]:
Z m a x _ r e f = 0.42 d 50 ( h c d 50 ) 1.7 ( Δ H h d ) 0.3
where hc = (q2/g)1/3 is the critical depth, with q indicating the unit discharge, g the acceleration due to gravity, ΔH the change in energy head, and hd is the tailwater depth. Figure 8 show the comparison between measured and calculated (using Equation (3)) values of the variable Zmax.

3.3. Apron and Cutoff Wall Design Guidance

The design and construction of an apron and cutoff wall must first consider a practical perspective of the entire structure design and how it may be constructed. Specifically, the PK Weir would require a foundation if placed in a river where it cannot be founded on competent rock. Thus, an economic analysis coupled with corresponding risks would inform the selection of (a) a deeper PK weir foundation while allowing a scour hole versus (b) designing downstream protection as part of the structure. Thus, when considering scour protection, the proper design of an apron and cutoff wall can help protect run-of-river structures from detrimental scour events associated with hydraulic conditions during flood events. In the previous section, it was shown that if Z*max exceeds 75%, the effect of the apron length becomes negligible, i.e., a further longitudinal extension of the apron results in a very marginal additional reduction of the scour depth (Figure 5a,b). Therefore, for practical applications, it can be reasonably assumed that a design value of the apron length for which Z*max = 75. Consequently, Equation (2) can be rearranged as follows:
( 2.8 F r d 90 6.3 ) ( L a P ) 2 + ( 2.8 F r d 90 + 64.9 ) ( L a P ) 75 = 0
and, for 0 ≤ La/P ≤ 2:
L a P = ( 2.8 F r d 90 + 64.9 ) + ( 2.8 F r d 90 + 64.9 ) 2 + 300 ( 2.8 F r d 90 6.3 ) 2 ( 2.8 F r d 90 6.3 )
Equation (6) provides a tool to estimate the optimum relative length of the apron La/P in order to maximize the scour reduction and minimize the construction cost. Note that Equation (6) is valid in the same ranges of parameters of Equation (3). It is worth remarking that Equation (6) does not depend on the tailwater level, which is one of the most important parameters in scour processes as the shear stress acting on the movable bed scales with the inverse of the sum of the water depth and scour depth. However, as mentioned in one of the previous sections, the range of variation of this parameter was relatively limited. Consequently, under tested conditions, its effect was found to be negligible. Overall, Equation (6) represents a useful tool as it can provide good estimations of the scour characteristics under most critical conditions (i.e., design condition). However, further investigations are needed to assess the effect of the tailwater depth. Figure 9 shows the plot of Equation (6), i.e., La/P as function of the densimetric Froude number Frd90.
Finally, in practice for erodible foundations a cutoff wall is often included both for the PK weir and then again at the end of a concrete apron or slab to protect against undermining. Based on the results of these tests, the maximum scour depth occurs in correspondence or in proximity of the edge of the apron. Therefore, the length or depth Lc of a cutoff wall should be designed in such a way that it is larger than the maximum scour depth Zmax calculated using Equation (3). To this end, the following is proposed:
L c = F S × Z m a x
where FS is the factor of safety selected by the designer according to site conditions and project objectives; note that when considering the deviation between measured and calculated values of Zmax, it is less than 30% for nearly all cases (see Figure 8) and thus a minimum FS would be 1.3. It is worth mentioning that the length of the cutoff wall Lc accounts for both inflow conditions and geometry of the structure, as such variables are embedded in Zmax. Further considerations for selection of Lc might include depth to a nonerodable geologic feature to key into, any relevant field observations, and balancing risks and project economy (i.e., excavation, dewatering, control of water, etc.).
Finally, it is worth remarking upon some site-specific design details that are not discussed herein, such as the means and methods of construction (i.e., in the dry with diversion or in the wet with excavation, sheet piles, etc.), the selection of the apron elevation (herein set at the base of the PK weir), the structural design of the apron and cutoff wall (dynamic forces, uplift pressures and loading moments, structural joints and connections, abrasion, etc.), and any changes to the existing mobile bed elevation downstream of the structure such as by general bed degradation independent of the weir, lateral channel migration, or downstream-instigated headcuts that may migrate upstream to the PK weir apron, etc. In addition, all the proposed equations from this study should be validated using field data, as they can be subjected to scale effects, regardless of the model scale. However, for the present study, a large laboratory flume was used. Considering that the width of natural streams/rivers (where such structures are located) is generally less than 50 m, the model scale can be assumed to be bigger than 1:25. Furthermore, tests were conducted for a Reynolds number greater than 50,000. Under such conditions, [46] showed that scale effects are negligible.

4. Conclusions

Limited published guidance is available for designing a minimal apron and cutoff wall lengths to protect or mitigate local scour at the toe of a PK weir. This encouraged a large-scale physical model study on the effects of scour at PK weirs with and without aprons and cutoff walls. From this study the following conclusions have been made:
  • Scour depths in the two tested noncohesive gravel materials without an apron and cutoff wall were significant.
  • Jets ensuing from the PK weir structure diffused over the apron and cause horizontal jet scour downstream of the apron.
  • Adding a 1.0P long apron can reduce scour by an average of 61%, adding a 1.5P long apron can reduce scour by an average of 75%, and adding a 2.0P apron can reduce scour by an average of 83%. Thus, apron length La can be optimized considering PK weir hydraulics and geometry with geologic conditions at PK weir toe.
  • Equations (3) and (6) were created to estimate La, Zmax downstream of an apron along with Equation (4), and Equation (7) was created to estimate Lc to protect the structure from scour undermining. Note that Equations (3) and (6) are considered valid for 2.4 ≤ Frd90 ≤ 5.3, 0 ≤ La/P ≤ 2, 0.11 ≤ H/P ≤ 0.36, and 0.24 ≤ hd/P ≤ 1.03. All proposed equations should be validated using field data.
  • For this study, it was determined that there was an 8% difference in the reduction of scour from a 1.5P apron length to a 2.0P apron length, and a 1.5P apron may be an adequate and cost-effective length to minimize scour. Note that as the substrate diameter decreases the potential for scour depth and length downstream of the apron will increase.
  • Although the facility tested at USU was quite large and might be considered a 1:5 to 1:10 scaled model, quantification of any scale effects requires additional investigation.
As with other laboratory-scale scour studies, there are multiple limitations to this study, yet there are ways that practitioners can overcome these limitations. First, the steady-state nature of the testing: To overcome this limitation, practitioners can use Zmax values as conservative design values. Another limitation was that only three Q and one hd conditions were tested in this study. To overcome this limitation, practitioners could use graphical means to interpolate potential Zmax values based on the rating curves provided and might adjust their FS if tailwater depths are greater. Further, results may be scaled using the densimetric Froude number. Furthermore, only two relatively uniform substrate materials were studied. Practitioners could use methods similar to the Erodibility Index method of [47] to scale geologic material strengths in the field to a noncohesive laboratory-scale material, such as was selected and tested herein. Additionally, these results and equations only directly apply to horizontal aprons without a sill and very mild river slopes or energy gradients. Lastly, only one PK weir geometry was tested with q estimated as Q/W.

Author Contributions

Conceptualization, B.M.C., W.L., and M.P.; methodology, B.M.C., W.L., and M.P.; software, W.L.; validation, B.M.C., W.L., and M.P.; formal analysis, B.M.C., W.L., and M.P.; investigation, B.M.C., W.L., and M.P.; resources, B.M.C.; data curation, W.L.; writing—original draft preparation, B.M.C., W.L., and M.P.; writing—review and editing, B.M.C. and M.P.; visualization, W.L.; supervision, B.M.C.; project administration, B.M.C.; funding acquisition, B.M.C. All authors have read and agreed to the published version of the manuscript.

Funding

For this research, the second author (Brian Crookston) was funded by USGS Grant 104b through the Utah Water Research Laboratory and by a generous private donation by D. Campbell.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Some or all data, models or code that support the findings of this study are available at https://digitalcommons.usu.edu/etd/8090/ (accessed on 4 March 2021) or from the corresponding author upon reasonable request.

Acknowledgments

The second author (Brian Crookston) thanks D. Campbell for his generous private donation to support this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Laursen, E.M. Observations on the Nature of Scour. In Proceedings 5th Hydraulics Conference Bulletin; University of Iowa: Iowa City, IA, USA, 1952; pp. 179–197. [Google Scholar]
  2. Ettema, R.; Yoon, B.; Nakato, T.; Muste, M. A review of scour conditions and scour-estimation difficulties for bridge abutments. Water Eng. 2004, 8, 643–650. [Google Scholar] [CrossRef]
  3. Zehe, E.; Becker, R.; Bárdossy, A.; Plate, E. Uncertainty of simulated catchment runoff response in the presence of threshold processes: Role of initial soil moisture and precipitation. J. Hydrol. 2005, 315, 183–202. [Google Scholar] [CrossRef]
  4. Bombardelli, F.A.; Palermo, M.; Pagliara, S. Temporal evolution of jet induced scour depth in cohesionless granular beds and the phenomenological theory of turbulence. Phys. Fluids 2018, 30, 085109. [Google Scholar] [CrossRef]
  5. Palermo, M.; Pagliara, S.; Bombardelli, F.A. Theoretical Approach for Shear-Stress Estimation at 2D Equilibrium Scour Holes in Granular Material due to Subvertical Plunging Jets. J. Hydraul. Eng. 2020, 146, 04020009. [Google Scholar] [CrossRef]
  6. Hoffmans, G.J.C.M.; Verheij, H.J. Scour Manual; Balkema, A.A., Ed.; Taylor and Francis Group: Rotterdam, The Netherlands, 1997. [Google Scholar]
  7. Adduce, C.; Sciortino, G. Scour due to a horizontal turbulent jet: Numerical and experimental investigation. J. Hydraul. Res. 2006, 44, 663–673. [Google Scholar] [CrossRef]
  8. Dey, S.; Raikar, R. Scour below a High Vertical Drop. J. Hydraul. Eng. 2007, 133, 564–568. [Google Scholar] [CrossRef]
  9. Ben Meftah, M.; Mossa, M. New Approach to Predicting Local Scour Downstream of Grade-Control Structure. J. Hydraul. Eng. 2020, 146, 04019058. [Google Scholar] [CrossRef]
  10. Lantz, W.; Crookston, B.M.; Palermo, M. Flood Infrastructure: Localized Scour at Piano Key Weirs. In Dam Safety 2020: Conference Proceedings; Association of State Dam Safety Officials: Lexington, KY, USA, 2020. [Google Scholar]
  11. United States Department of the Interior Bureau of Reclamation (USBR). Engineering Monograph No. 25: Hydraulic Design of Stilling Basins and Energy Dissipators; U.S Bureau of Reclamation: Denver, CO, USA, 1984.
  12. Thompson, P.; Kilgore, R. Hydraulic Design of Energy Dissipators for Culverts and Channels. In Federal Highway Administration; National Highway Institute: Washington, DC, USA, 2006. [Google Scholar]
  13. United States Department of Agriculture (USDA). The SAF Stilling Basin; U.S. Government Printing Office: Washington, DC, USA, 1959.
  14. Peterka, J. Hydraulic Design of Stilling Basins and Energy Dissipators; Bureau of Reclamation: Denver, CO, USA, 1958.
  15. Vischer, D.L.; Hager, W.H. Energy Dissipators: IAHR Hydraulic Structures Design Manuals 9; Taylor and Francis: Abingdon, UK, 1995. [Google Scholar]
  16. Novak, P.; Moffat, A.; Nalluri, C.; Narayanan, R. Hydraulic Structures, 4th ed.; E & FN Spon: Abingdon, UK, 1995. [Google Scholar]
  17. Khatsuria, R.M. Hydraulics of Spillways and Energy Dissipators; Marcel Decker: New York, NY, USA, 2005. [Google Scholar]
  18. Hager, W.H.; Schleiss, A.J.; Boes, R.M.; Pfister, M. Hydraulic Engineering of Dams; CRC Press: London, UK, 2020. [Google Scholar]
  19. Ho Ta Khanh, M.; Truong, C.H.; Nguyen, T.H. Main Results of the P.K Weir Model Tests in Vietnam (2004–2010); Labyrinth and Piano Key Weirs: PKW 2011; Taylor and Francis Group: London, UK, 2011; pp. 191–198. [Google Scholar]
  20. Crookston, B.M.; Erpicum, S.; Tullis, B.P.; Laugier, F. Hydraulics of Labyrinth and Piano Key Weirs: 100 Years of Prototype Structures, Advancements, and Future Research Needs. J. Hydraul. Eng. 2019, 145, 02519004. [Google Scholar] [CrossRef]
  21. Urban and Environmental Engineering (UEE). World Register of Piano Key Weirs Prototypes. Liege Universite. 2021. Available online: https://www.uee.uliege.be/cms/c_5026433/en/world-register-of-piano-key-weirs-prototypes (accessed on 11 March 2021).
  22. Schleiss, A.J. From Labyrinth to Piano Key Weirs—A Historical Review; Labyrinth and Piano Key Weirs: PKW 2011; Taylor and Francis Group: London, UK, 2011; pp. 3–16. [Google Scholar]
  23. Ribeiro, M.L.; Pfister, M.; Schleiss, A.J.; Boillat, J.L. Hydraulic design of A-Type Piano Key Weirs. J. Hydraul. Res. 2012, 50, 400–408. [Google Scholar] [CrossRef] [Green Version]
  24. Machiels, O.; Pirotton, M.; Pierre, A.; Dewals, B.; Erpicum, S. Experimental parametric study and design of Piano Key Weirs. J. Hydraul. Res. 2014, 52, 326–335. [Google Scholar] [CrossRef]
  25. López-Soto, J.; Wibowo, J.; Molina-Bas, O. Cost Reduction in Dam Infrastructure Using Arced Labyrinth Spillways. In Construction Research Congress 2016; San Juan, Puerto Rico, American Society of Civil Engineers (ASCE): Reston, VA, USA, 2016. [Google Scholar] [CrossRef]
  26. Laugier, F.; Vermeulen, J.; Lefebvre, V. Overview of Piano Key Weirs Experience Developed at EDF during the Past Few Years; Labyrinth and Piano Key Weirs II: PKW 2013; Taylor and Francis Group: London, UK, 2013; pp. 213–226. [Google Scholar]
  27. Jüstrich, S.; Pfister, M.; Schleiss, A.J. Mobile Riverbed Scour Downstream of a Piano Key Weir. J. Hydraul. Eng. 2016, 142, 04016043. [Google Scholar] [CrossRef]
  28. Pfister, M.; Jüstrich, S.; Schleiss, A.J. Toe-Scour Formation at Piano Key Weirs; Labyrinth and Piano Key Weirs III—PKW 2017; Taylor and Francis Group: London, UK, 2017; pp. 147–156. [Google Scholar]
  29. Palermo, M.; Crookston, B.; Pagliara, S. Analysis of Equilibrium Morphologies Downstream of a PK Weir Structure. In World Environmental and Water Resources Congress 2020; American Society of Civil Engineers (ASCE): Reston, VA, USA, 2020; pp. 43–51. [Google Scholar] [CrossRef]
  30. Lantz, W. A Laboratory Study on the Geometric Effects of Piano Key Weirs on Scour for Non-Cohesive Substrates and Simple Mitigation Techniques. Master’s Thesis, Utah State University, Logan, UT, USA, 2021. [Google Scholar] [CrossRef]
  31. Yazdi, A.M.; Hoseini, S.A.; Nazari, S.; Amanian, N. Effects of weir geometry on scour development in the downstream of Piano Key Weirs. Water Supply 2021, 21, 289–298. [Google Scholar] [CrossRef]
  32. Hassan, N.M.K.N.; Narayanan, R. Local Scour Downstream of an Apron. J. Hydraul. Eng. 1985, 111, 1371–1384. [Google Scholar] [CrossRef]
  33. Chatterjee, S.S.; Ghosh, S.N.; Chatterjee, M. Local Scour due to Submerged Horizontal Jet. J. Hydraul. Eng. 1994, 120, 973–992. [Google Scholar] [CrossRef]
  34. Sarkar, A.; Dey, S. Scour downstream of aprons caused by sluices. Proc. Inst. Civ. Eng. Water Manag. 2005, 158, 55–64. [Google Scholar] [CrossRef]
  35. Dey, S.; Sarkar, A. Scour Downstream of an Apron Due to Submerged Horizontal Jets. J. Hydraul. Eng. 2006, 132, 246–257. [Google Scholar] [CrossRef]
  36. Erpicum, S.; Laugier, F.; Boillat, J.L.; Pirotton, M.; Reverchon, B.; Schleiss, A.J. Labyrinth and Piano Key Weirs—PKW 2011; CRC Press: London, UK, 2011. [Google Scholar]
  37. Microsonic. “mic+ 130/IU/TC.” mic+ Ultrasonic Sensors. 2021. Available online: https://www.microsonic.de/en/distance-sensors/cylindrical/micplus/standard-sensors/standard-sensors/micplus130iutc.htm (accessed on 4 March 2021).
  38. Zhang, G.; Valero, D.; Bung, D.; Chanson, H. On the estimation of free-surface turbulence using ultrasonic sensors. Flow Meas. Instrum. 2018, 60, 171–184. [Google Scholar] [CrossRef] [Green Version]
  39. Intel. “Depth Camera D435.” Intel RealSense. 2021. Available online: https://www.intelrealsense.com/depth-camera-d435/ (accessed on 4 March 2021).
  40. Bung, D.B.; Crookston, B.M.; Valero, D. Turbulent free-surface monitoring with an RGB-D sensor: The hydraulic jump case. J. Hydraul. Res. 2020, 1–12. [Google Scholar] [CrossRef]
  41. Oliveto, G.; Hager, W.H. Temporal Evolution of Clear-Water Pier and Abutment Scour. J. Hydraul. Eng. 2002, 128, 811–820. [Google Scholar] [CrossRef]
  42. Neill, C.R. Note on Initial Movement Of Coarse Uniform Bed-Material. J. Hydraul. Res. 1968, 6, 173–176. [Google Scholar] [CrossRef]
  43. Julien, P.Y. Erosion and Sedimentation, 2nd ed.; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  44. Pagliara, S.; Palermo, M. Scour control and surface sediment distribution downstream of block ramps. J. Hydraul. Res. 2008, 46, 334–343. [Google Scholar] [CrossRef]
  45. Pagliara, S.; Roy, D.; Palermo, M. 3D plunge pool scour with protection measures. J. Hydro-Environ. Res. 2010, 4, 225–233. [Google Scholar] [CrossRef]
  46. Heller, V. Scale effects in physical hydraulic engineering models. J. Hydraul. Res. 2011, 49, 293–306. [Google Scholar] [CrossRef]
  47. Annandale, G.W. Erodibility. J. Hydraul. Res. 1995, 33, 471–494. [Google Scholar] [CrossRef]
Figure 1. (a) Plan and (b) side views of weir dimensions and flume setup.
Figure 1. (a) Plan and (b) side views of weir dimensions and flume setup.
Water 13 02332 g001
Figure 2. Profile of experimental setup with key parameters used in this study.
Figure 2. Profile of experimental setup with key parameters used in this study.
Water 13 02332 g002
Figure 3. Photographs from experiments under identical hydraulic conditions: (a) Reference test (Run 14); (b) Scour protection via horizontal apron and cutoff wall (Run 23).
Figure 3. Photographs from experiments under identical hydraulic conditions: (a) Reference test (Run 14); (b) Scour protection via horizontal apron and cutoff wall (Run 23).
Water 13 02332 g003
Figure 4. Comparison of longitudinal bed profiles for runs: (a) 1, 4, and 7; (b) 13, 16, 19 and 22; (c) 2, 5, 8 and 11; (d) 14, 17, 20 and 23; (e) 3, 6, 9, 9R and 12; and (f) 15, 18, 21 and 24.
Figure 4. Comparison of longitudinal bed profiles for runs: (a) 1, 4, and 7; (b) 13, 16, 19 and 22; (c) 2, 5, 8 and 11; (d) 14, 17, 20 and 23; (e) 3, 6, 9, 9R and 12; and (f) 15, 18, 21 and 24.
Water 13 02332 g004
Figure 5. Z*max as function of Frd90 for Substrate (a) 1 and (b) 2. Z*max as function of Frd90 along with the plot of Equation (1) for Frd90 = (c) 2.44, (d) 3.05, (e) 3.99, and (f) 4.98.
Figure 5. Z*max as function of Frd90 for Substrate (a) 1 and (b) 2. Z*max as function of Frd90 along with the plot of Equation (1) for Frd90 = (c) 2.44, (d) 3.05, (e) 3.99, and (f) 4.98.
Water 13 02332 g005
Figure 6. Effect of La on Zmax under different inflow conditions and Substrates (a) 1 and (b) 2. Note that 0P indicates corresponding reference tests.
Figure 6. Effect of La on Zmax under different inflow conditions and Substrates (a) 1 and (b) 2. Note that 0P indicates corresponding reference tests.
Water 13 02332 g006
Figure 7. Effect of La on Lmax under different inflow conditions and Substrates (a) 1 and (b) 2. Note that 0P indicates corresponding reference tests.
Figure 7. Effect of La on Lmax under different inflow conditions and Substrates (a) 1 and (b) 2. Note that 0P indicates corresponding reference tests.
Water 13 02332 g007
Figure 8. Comparison between measured and calculated (using Equation (3)) values of the maximum scour depth Zmax for protected basins, with 30% bounding lines referencing Zmax Observed.
Figure 8. Comparison between measured and calculated (using Equation (3)) values of the maximum scour depth Zmax for protected basins, with 30% bounding lines referencing Zmax Observed.
Water 13 02332 g008
Figure 9. Optimal design values of the variable La/P as function of Frd90.
Figure 9. Optimal design values of the variable La/P as function of Frd90.
Water 13 02332 g009
Table 1. PK weir geometric characteristics.
Table 1. PK weir geometric characteristics.
ParameterValue
B1.04 m
Bb0.52 m
Bi = Bo0.26 m
N4
P0.42 m
Pd1.09 m
Si = So0.55
Ts0.025 m
Wi0.248 m
Wo0.194 m
Wi/Wo1.280 m
Wu0.49 m
Table 2. Intel RealSense D435 specifications.
Table 2. Intel RealSense D435 specifications.
FeatureDetail
Global Shutter3 μm × 3 μm pixel size
IR Stereo: FOV86° × 57° (±3°)
IR Stereo: Resolution1280 × 720
RGB: FOV64°×41° × 77° (±3°)
RGB: Resolution1920 × 1080
Table 3. PK Weir Scour Protection Test Matrix.
Table 3. PK Weir Scour Protection Test Matrix.
Substrate TypeDischarge
Q
Headwater
Ratio
Tailwater †
hd
Apron Length
La
Substrate 1
d50 = 13 mm
d90 = 20 mm
σ = 1.54
Substrate 2
d50 = 6.5 mm
d90 = 9.1 mm
σ = 1.30
150 L/sH/P = 0.1116 cm = 0.38P2.0P = 0.84 m
1.5P = 0.63 m
1.0P = 0.42 m
0.0P = 0.00 m *
300 L/sH/P = 0.1817 cm = 0.41P2.0P = 0.84 m
1.5P = 0.63 m
1.0P = 0.42 m
0.0P = 0.00 m *
600 L/sH/P = 0.3522 cm = 0.52P2.0P = 0.84 m
1.5P = 0.63 m
1.0P = 0.42 m
0.0P = 0.00 m *
† reference test tailwater, see Table 5 for all tests; * reference test without apron, or apron length La = 0.
Table 4. Testing values for each experimental run.
Table 4. Testing values for each experimental run.
Run
(#)
Substrate
(-)
t
(min)
La/P
(-)
Q
(L/s)
1 *11800.0150
2 *12400.0300
3 *19000.0600
413601.0150
513601.0300
618701.0600
713601.5150
813601.5300
913601.5600
9R13601.5600
101no scour2.0150
1111202.0300
1213602.0600
13 *24500.0150
14 *210500.0300
15 *211700.0600
1628401.0150
1727801.0300
1827201.0600
1922401.5150
2028401.5300
2129901.5600
2224802.0150
2326002.0300
2427202.0600
R = repeat test; * reference test with no apron, or apron length La = 0.
Table 5. Experimental Results.
Table 5. Experimental Results.
Run
(#)
hd
(m)
Zmax
(m)
Xmax
(m)
Lmax
(m)
V
(m3)
1 *0.160.180.250.610.12
2 *0.170.330.421.070.36
3 *0.220.710.732.201.49
40.090.070.000.470.03
50.160.150.000.750.12
60.210.390.001.540.50
70.100.030.000.310.00
80.170.060.000.440.03
90.220.280.001.540.37
9R0.230.230.001.450.37
100.09No ScourNo ScourNo ScourNo Scour
110.150.060.000.360.02
120.220.190.001.180.27
13 *0.100.280.401.070.39
14 *0.140.560.521.811.03
15 *0.221.011.163.193.36
160.110.110.000.610.07
170.160.200.001.150.26
180.240.450.001.830.74
190.110.060.000.470.03
200.150.120.001.280.15
210.250.340.002.370.56
220.100.040.000.450.02
230.150.110.001.020.12
240.240.250.003.050.45
R = repeat test; * reference test with no apron, or apron length La = 0.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lantz, W.; Crookston, B.M.; Palermo, M. Apron and Cutoff Wall Scour Protection for Piano Key Weirs. Water 2021, 13, 2332. https://doi.org/10.3390/w13172332

AMA Style

Lantz W, Crookston BM, Palermo M. Apron and Cutoff Wall Scour Protection for Piano Key Weirs. Water. 2021; 13(17):2332. https://doi.org/10.3390/w13172332

Chicago/Turabian Style

Lantz, Wyatt, Brian Mark Crookston, and Michele Palermo. 2021. "Apron and Cutoff Wall Scour Protection for Piano Key Weirs" Water 13, no. 17: 2332. https://doi.org/10.3390/w13172332

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