Labyrinth weirs are commonly used to increase the capacity of existing spillways and provide more efficient spillways for new dams due to their high specific discharge capacity compared to the linear weir. In the present study, an experimental and numerical investigation was conducted to improve the rectangular labyrinth weir performance. In this context, four configurations were tested to evaluate the influence of the entrance shape and alveoli width on its discharge capacity. The experimental models, three models of rectangular labyrinth weir with a rounded entrance and one with a flat entrance, were tested in rectangular channel conditions for inlet width to outlet width ratios (a/b) equal to 0.67, 1 and 1.5. The results indicate that the rounded entrance increases the weir efficiency by up to 5%. A ratio a/b equal to 1.5 leads to an 8 and 18% increase in the discharge capacity compared to a/b ratio equal to 1 and 0.67, respectively. In addition, a numerical simulation was conducted using the open-source Computational Fluid Dynamics (CFD) OpenFOAM to analyze and provide more information about the flow behavior over the tested models. A comparison between the experimental and numerical discharge coefficient was performed and good agreement was found (mean absolute relative error of 4–6%).

  • This paper presents the efficiencies of rectangular labyrinth weirs with flat and rounded entrances for channel applications.

  • The rectangular labyrinth weir efficiency is mainly dependent on its inlet alveoli capacity.

  • Both numerical and experimental models provided close results.

  • An InterFoam solver produced more information about the effect of inlet shape and the alveoli widths.

Graphical Abstract

Graphical Abstract
Graphical Abstract

The labyrinth weir is a non-linear flood control structure providing an increase in the unit discharge capacity relative to linear weirs (Figure 1). The use of non-linear shape in plan view enables the labyrinth weir crest length to be increased within the existing spillway width. This option leads to an augmentation in the discharge capacity of the labyrinth weir by up to three to four times compared with a linear weir having the same width and upstream head (Tullis et al. 1995). Consequently, the adoption of a labyrinth weir can effectively reduce a direct and/or indirect cost of the projects of rehabilitation or construction of new dams. Contrary to a labyrinth weir, a linear weir requires a larger width or low crest elevation to increase the discharge efficiency due to its low specific discharge capacity. For example, the specific discharge for a standard Creager type is close to 2.15 h1.5, where h is the upstream head in meters (Hydrocoop 2013). A greater width of the crest of the spillway and downstream channel result in an increase in the volume of materials required for construction, and can directly affect the total cost of the structure. The construction of the crest of the dam with a lower level can reduce the normal elevation of the reservoir and thus the storage capacity of the dam, which may have an indirect negative impact on the project cost. As a result, the labyrinth weir is well suited for a dam or in a natural or artificial channel application to increase at low cost the spillway evacuation and/or storage capacity. The main geometric parameters of a labyrinth weir are the developed crest length L, the total width Wt, the cycle width w, the weir height P, the total upstream Ht, cycle number N, and the sidewall angle α. The first reported application of labyrinth weirs was conducted by Murphy (1909) to provide a high discharge capacity for the reinforced-concrete spillway of the Keno Canal, in the USA. It is highlighted that a linear weir of equivalent discharge capacity would need to be three times the length of the adopted labyrinth weir. However, the first study of the labyrinth weir performance was performed by Gentilini (1940). He tested a labyrinth weir designed with triangular shape in plan view and compared it with a linear weir capacity. Nevertheless, the significant development in the understanding of the effects of the labyrinth weir geometry on discharge capacity was conducted by Taylor (1968) and Hay & Taylor (1970), specifically for trapezoidal and triangular labyrinths. Later, an extensive investigation was mainly conducted on the trapezoidal labyrinth form by many researchers (e.g., Darvas 1971; Houston 1983; Hinchliff & Houston 1984; Lux & Hinchliff 1985; Tullis et al. 1995, 2007; Falvey 2003; Lopes et al. 2009; Crookston & Tullis 2012; Dabling et al. 2013). It was found in those works that the performance of trapezoidal labyrinth weirs is mainly governed by the total relative head ratio Ht/P, the developed crest length ratio, L/W, or the sidewall angle, α, the cycle width ratio w/P, and the crest shape.

Figure 1

Main geometrical parameters of a rectangular labyrinth weir.

Figure 1

Main geometrical parameters of a rectangular labyrinth weir.

Close modal

Recently, a numerical simulation method has been widely used to further analyze and design the labyrinth weir. This method has been applied on the trapezoidal labyrinth form by Savage et al. (2016) and it provided a good result compared to the experimental data (error of 3%). Carrillo et al. (2019) carried out a numerical simulation to analyze the trapezoidal labyrinth efficiency for free and submerged conditions. Ghaderi et al. (2020) performed numerical and experimental tests to analyze the energy dissipation and hydraulics of flow over trapezoidal and triangular labyrinth shape. They conclude that the numerical model is suitable to accurately predict the flow behavior of the labyrinth weir with a relative error of 3.05%. Safarzadeh & Noroozi (2017) numerically investigated the flow features over the rectangular labyrinth and designed Piano Key Weir with trapezoidal shape in plan view.

Rectangular labyrinth weir performance

The rectangular labyrinth has a simple geometry, combines high performance and a low construction cost. However, far too little attention has been paid to the study of the performance of labyrinth weir having a rectangular shape in plan view. Kabiri-Samani et al. (2013) studied a rectangular labyrinth weir designed with a flat entrance. The influence of the free and submerged flow conditions on the discharge capacity of the rectangular labyrinth weir has been investigated by Azimi & Seyed Hakim (2019). They indicated that the rectangular labyrinth weir is more efficient than the sharp-crested weir and it is less sensitive to submergence and it behaves ideally within the upstream relative head (h0) range of 0.3≤h0/P≤0.4.

Lempérière & Ouamane (2003), Ouamane & Ben Said (2010) and (Ben Said & Ouamane 2011, 2018) studied several physical models of rectangular labyrinth weirs designed with entrances having a rounded shape and found that the rectangular labyrinth weir is more efficient than the trapezoidal weir mainly for moderate discharge. The rectangular labyrinth shape was found to provide the highest discharge capacity compared to the trapezoidal shape for L/W ratio equal to 4 and 5. Using a rectangular shape enables a reduction in the footprint of the weir, which reduces the construction cost and allows it to be installed in a limited area. In addition, through choosing a rectangular shape, the free height walls can be decreased by the partial filling of the downstream alveoli with ordinary concrete, allowing the steel-reinforced concrete to be reduced, which can result in a reduction in the amount of the construction cost without compromising its hydraulic efficiency (Ouamane & Ben Said 2010; Ben Said & Ouamane 2011). Figure 1 shows the main geometric parameters of the rectangular labyrinth weir including the total width of the weir Wt, the cycle width wu, the weir height P, the inlet alveoli width a, the outlet alveoli width b, and the lateral crest length B. In addition, the radius r of the rounded entrance was considered for the rectangular labyrinth having rounded shape (see Figure 4).

Table 1

Characteristics of models

ModelCycle numberL (cm)Wt (cm)P (cm)B (cm)Wu (cm)a (cm)b (cm)r (cm)L/WtW/Pa/bB/P
Flat entrance 351 90 15 25 58.5 7.5 7.5 – 3.9 1.67 
Rounded entrance 353 90.3 15 25 58.83 4.5 3.91 0.67 1.67 
Rounded entrance 355 90 15 25 58.5 7.5 7.5 3.75 3.9 1.67 
Rounded entrance 355 90 15 25 59.16 3.91 1.5 1.67 
ModelCycle numberL (cm)Wt (cm)P (cm)B (cm)Wu (cm)a (cm)b (cm)r (cm)L/WtW/Pa/bB/P
Flat entrance 351 90 15 25 58.5 7.5 7.5 – 3.9 1.67 
Rounded entrance 353 90.3 15 25 58.83 4.5 3.91 0.67 1.67 
Rounded entrance 355 90 15 25 58.5 7.5 7.5 3.75 3.9 1.67 
Rounded entrance 355 90 15 25 59.16 3.91 1.5 1.67 

The main objective of this study was to analyze the rectangular labyrinth weir discharge efficiency using experimental and numerical modeling. For that, the sensitivity of the rectangular labyrinth weir to the entrance shape and alveoli width was measured on four models of rectangular labyrinth weirs. Firstly, two models with the same a/b ratio = 1, were tested to investigate the influence of the entrance shape on the rectangular labyrinth. Then, the optimal model from the previous comparison was chosen as the reference model and compared to the two other models with ratios a/b = 0.67 and 1.5, to investigate the influence of the variation of the alveoli width on the rectangular labyrinth.

Labyrinth weir discharge capacity

The head discharge relationship for a linear sharp-crested weir is generally used to quantify the labyrinth weir capacity:
(1)
where
  • Q: discharge over weir (m3/s)

  • Cw: discharge coefficient

  • g: gravity acceleration (m/s2)

  • W: weir width (m)

  • H: total upstream head (m).

The height, the crest shape, and the crest length are the main parameters influencing the linear weir capacity (Darvas 1971). However, the principal geometric parameters affecting the performance of labyrinth weirs are: the total weir crest length L, the weir width W, the weir height P and the inlet and outlet alveoli widths a and b.

In this study, the discharge coefficient Cw is estimated by using the experimental H (H = h+U2/2 g) and Q data. Based on dimensional analysis, the discharge coefficient for rectangular labyrinth weir can be written as:
(2)
where
  • H: total upstream head

  • h: piezometric head

  • b: outlet alveoli width

  • Wc: channel width

  • a: inlet alveoli width

  • P: weir height

  • L: developed crest length

  • W: the weir width.

Experimental method

The experimental tests on physical models of the rectangular labyrinth weir were conducted at the laboratory of Hydraulic Planning and Environment of Biskra University in a straight rectangular channel 12 m long, 1 m wide and 1 m deep (Figure 2). The experimental system was equipped with two pumps delivering up to 180 l/s into an upstream stilling basin containing a perforated wall to ensure uniform approach flow conditions. A portion of the sidewalls of the channel was designed with glass to allow the observation of flow behavior (Figure 3). An electromagnetic and ultrasonic flowmeter for measuring discharge and a series of ultrasonic level sensors to measure the flow depth were the instruments of measurement used in this work. The assembly operates within a closed circuit. It begins when the water exits the storage tank to enter the flume through the stilling basin. Then, for each discharge rate, the inflows were allowed to cross the crest of weir installed on the flume bottom for five minutes before recording the data. Finally, the evacuated water flows over the downstream bed of the flume back to the storage tank. Four rectangular labyrinth weir models with six cycles were tested (see Figure 4) for 0.67 ≤ a/b ≤ 1.5 and carefully inserted and leveled. All labyrinth models with a weir height of P = 15 cm were fabricated using 0.2 cm thick metal sheeting (Table 1). For all tests conducted in this study, free flow conditions were adopted, and around 15 discharges between 30 and 150 l/s were considered for each model. The tests were repeated three times to check the accuracy of the measurement, which gives a total of 180 runs.

Figure 2

View of the experimental flume.

Figure 2

View of the experimental flume.

Close modal
Figure 3

Channel flow over labyrinth weir with rounded entrance.

Figure 3

Channel flow over labyrinth weir with rounded entrance.

Close modal
Figure 4

Schematic of the labyrinth weir geometries considered in the analysis of the influence of the inlet entrance and alveoli width on discharge efficiency (see Table 1 for further details).

Figure 4

Schematic of the labyrinth weir geometries considered in the analysis of the influence of the inlet entrance and alveoli width on discharge efficiency (see Table 1 for further details).

Close modal

Numerical method

Governing equations

Numerical modeling of flow over rectangular labyrinth weir was performed using OpenFOAM software, which is an open-source Computational Fluid Dynamics (CFD). The InterFoam solver, distributed with OpenFOAM for treating two-phase flow cases, was adopted in this study. This solver uses the Finite Volume Method (FVM) for solving the Reynolds-Averaged Navier–Stokes (RANS) equations, where the pressure implicit splitting of operators (PISO) algorithm is used to define the flow pressure p and velocity fields U (Issa 1986). InterFoam is able to capture the interface between two fluids by determining the volume fraction α. The value of α can have a value of 0 when cell contains only air, or 1 when cell is fully occupied with water. An intermediate value represents the cell where the free surface should be located. The standard k-ε turbulence model is one among others used with InterFoam for resolving the turbulence parameters and was selected here for its acceptable ability for simulating free surface flow (Beg et al. 2017; Duguay et al. 2017). This turbulence model has two equations systems, in which the transport equations should be solved for estimating the turbulent kinetic energy (k) and the dissipation rate of energy (ε). The general form of the equations for continuity and momentum conservation for incompressible flows are shown as Equations (3) and (4) (Moukalled et al. 2015):
(3)
(4)
where U is the flow velocity, is the density, t is the time, p is the pressure, is the shear-rate tensor, fb is the body forces.

Computational mesh and boundary conditions

The geometry of the tested models of the rectangular labyrinth weir was generated by using the open-source parametric 3D modeler freeCAD. Then, the blockMesh and snappyHexMesh utilities provided in OpenFOAM have been used to discretize the computational domain. First, a background hex mesh defining the computational domain was created by the blockMesh utility. Then, the snappyHexMesh utility was used to fit the geometry based on the mesh previously generated by blockMesh. The upstream and downstream parts of the numerical channel were approximately 1.8 m (12P) and 0.35 m, respectively. However, the height of the numerical domain was a function of the upstream head. The applied boundary conditions on the model were the inlet at the entrance of the channel, the atmosphere at the top of the channel, the walls for sidewalls, the bottom and the weir, and the outlet for the end of the channel. The boundary conditions are similar for all cases with the exception that the discharge (Q) is defined at the inlet.

A variable HeightFlowRateInletVelocity boundary condition was considered for U at the entrance of the domain. At the top of the channel a pressure InletOutletVelocity condition for U was defined. At the walls, the boundary conditions were applied as no-slip and a zero gradient condition was used at the outlet. A structured mesh was adopted to discretize the computational domain of simulation. In order to reduce the required time of computation for the simulations, a grid refinement was used where the weir is located to accurately identify the weir geometry and the flow behavior along the different part of the weir.

To investigate the influence of the mesh resolution, the Grid Convergence Index (GCI) method was used to estimate the accuracy of the numerical results. Consequently, three meshes with different cell sizes, 11.2, 9.4 and 7.8 mm, corresponding to 324,305, 554,924 and 953,255 cells, respectively, were used. A refinement factor between the coarse and fine gride, r, equal to 1.2 was considered, which is greater than 10% as recommended by Roache (1997) and adopted from previous studies (e.g., Aydin & Emin Emiroglu 2013; Aydin & Emin Emiroglu 2016).

The GCI was computed for the total upstream head H provided by the rectangular labyrinth weir with flat entrance using three rates of discharge equal to 37.1, 54.5 and 78.45 l/s.

Following this method which is outlined by Celik et al. (2008) based on Richardson extrapolation method, the representative cell size was defined as:
(5)
where ΔVi is volume of ith cell and N is number of cells. The apparent order P for the three-grid solution can be calculated by solving Equation (5):
(6)
where the value of the refinement factor r21 = and r32 = . The is the error between two adjacent meshes, where and . For this simulation, an average apparent order Pave equal to 6.39 was used to calculate the GCI. The GCI was estimated using the following relation:
(7)
where E is the approximate relative error between two meshes computed as:
(8)
f1, f2 and f3 are the key variables obtained from the numerical simulation on the three grids, where f is the total upstream head provided by the tested labyrinth weir.

Table 2 shows the GCI estimations and the relative change between two adjacent meshes are presented. GCI21 and GCI32 are the relative change from medium to coarse and from fine to medium mesh, respectively. Regarding these obtained results, the values of GCI32/rPGCI21 are close to 1, indicating that the results of the numerical model are in the asymptotic range. Consequently, a grid mesh with 554,924 cells was selected (Figure 5).

Table 2

Grid index convergence estimation

ModelDischarge Q (l/s)Total head Ht (m)
PaveE21E32GCI21GCI32GCI32/rp GCI21
f1f2f3
Rectangular labyrinth with flat entrance 37.1 0.0359 0.0358 0.0357 6.39 1.9973 1.9973 1.2244 1.2244 0.9973 
54.4 0.0558 0.0557 0.0561 0.1718 0.6886 0.1053 0.4221 1.0017 
78.45 0.0815 0.0814 0.0821 0.1151 0.8071 0.0706 0.4947 1.0012 
ModelDischarge Q (l/s)Total head Ht (m)
PaveE21E32GCI21GCI32GCI32/rp GCI21
f1f2f3
Rectangular labyrinth with flat entrance 37.1 0.0359 0.0358 0.0357 6.39 1.9973 1.9973 1.2244 1.2244 0.9973 
54.4 0.0558 0.0557 0.0561 0.1718 0.6886 0.1053 0.4221 1.0017 
78.45 0.0815 0.0814 0.0821 0.1151 0.8071 0.0706 0.4947 1.0012 
Figure 5

Numerical model mesh and boundary surface.

Figure 5

Numerical model mesh and boundary surface.

Close modal

Experimental results of the discharge coefficient

Entrance effect

Two models of rectangular labyrinth weirs were considered to test the influence of the entrance shape on its discharge capacity. These models have the same ratio L/W = 4 and the same alveoli width a/b = 1, with the exception of the entrance shape. The first model has a flat entrance and the second was designed with a rounded entrance. The obtained results have shown that the rectangular labyrinth weir with rounded entrance is slightly more efficient than the labyrinth with flat entrance. It provides a relatively constant gain in discharge efficiency around 5% (Figure 6(a)).

Figure 6

Experimental results of: (a) entrance effect for a/b = 1, (b) alveoli width effect for a/b: 0.67,1 and 1.5.

Figure 6

Experimental results of: (a) entrance effect for a/b = 1, (b) alveoli width effect for a/b: 0.67,1 and 1.5.

Close modal

Alveoli width effect a/b

In this investigation, the model of the rectangular labyrinth weir with a round entrance and symmetric alveoli width (a/b = 1) was chosen as the reference model to evaluate the influence of the alveoli width on the labyrinth weir efficiency. For that, two models of a rectangular labyrinth weir with a rounded entrance and varied alveoli width ratio a/b = 0.67 and 1.5, were compared to the refence model. The analysis of the experimental results showed that the rectangular labyrinth with ratio a/b = 1.5 gave the best discharge capacity mainly for moderate heads ratio (Ht/P). The model with a/b = 1.5 led to an increase by up to 8% compared to a ratio equal to 1 and 18% when the model has an a/b ratio = 0.67 as shown in Figure 6(b).

The labyrinth with alveoli width equal to 0.67 had the lowest weir capacity compared to the other models. This is mainly due to the reduction in the inlet alveoli width. Thus, as a/b values increase then the inlet alveoli capacity increases also, thus allowing an increase in the labyrinth weir efficiency. According to the studies that used a rectangular labyrinth weir with a rounded entrance, the optimal value of the ratio a/b is close to 1.5. However, increasing a/b further led to the reduction of the outlet alveoli capacity and thus the global weir efficiency (Ben Said & Ouamane 2011, 2018).

Numerical results

Discharge coefficient and validation

The comparison of the discharge coefficient obtained from the CFD simulations and experimental tests are shown in Figure 7(a)–7(d) for ratio width a/b equal to 0.67,1,1.5 and 1 (for flat entrance), respectively. To quantify the difference between the computed and tested results, the estimation of the Mean Absolute Relative Error (MARE) was performed for the discharge coefficient. The MARE was computed using:
(9)
where Oi and Si are the experimental and simulated results, respectively, and N is the total number of data.
Figure 7

Comparison between experimental and numerical results for: (a) a/b = 0.67, (b) a/b = 1, (c) a/b = 1.5 and (d) a/b = 1 (flat entrance).

Figure 7

Comparison between experimental and numerical results for: (a) a/b = 0.67, (b) a/b = 1, (c) a/b = 1.5 and (d) a/b = 1 (flat entrance).

Close modal

In general, a good agreement was found between the experimental and simulated results especially for models having the ratio a/b = 1, 1.5 and 1 (flat entrance). However, the less agreement obtained from the model has a/b = 0.67, mainly for low upstream heads, may be attributed to the geometric accuracy of the experimental model. The evaluated value of MARE for a/b = 0.67 was 6.31% (ranging from 0.06 to 18%), a/b = 1 was 5.56% (ranging from 0.7 to 9.6%), a/b = 1.5 was 6.32% (ranging from 3.16 to 9.68) and for the model having flat entrance and ratio a/b = 1 the MARE was 4.161% (ranging from 0.47 to 6.9%).

Figure 8 presents the experimental and numerical flow over the rounded rectangular labyrinth weir having a ratio a/b = 1.5 for rate of discharge Q = 80 l/s. It shows that the CFD model provides flow behaviors over the weir close to the ones obtained from the experimentation.

Figure 8

Example of experimental and numerical flow over labyrinth weir for Q = 80 l/s: (a) view of the downstream, (b) view of upstream.

Figure 8

Example of experimental and numerical flow over labyrinth weir for Q = 80 l/s: (a) view of the downstream, (b) view of upstream.

Close modal

Entrance effect

In order to analyze the effect of the entrance on the labyrinth weir performance, discharge coefficients were numerically computed for two models with different inlet forms and having the same ratio a/b = 1.

As previously shown from the experimental results, the entrance shape has an effect on the discharge of the rectangular labyrinth weir, where the model designed with the rounded entrance is around 5% more efficient than the model with the flat entrance. The numerical results for the discharge coefficient are presented in Figure 9 for two configurations. This shows that the model designed with the rounded entrance is more efficient than the model designed with the flat entrance, which is in agreement with the results provided by the experimental model.

Figure 9

Discharge coefficient obtained from numerical simulation for the entrance effect.

Figure 9

Discharge coefficient obtained from numerical simulation for the entrance effect.

Close modal

To further analyze the effect of the entrance shape on the labyrinth weir capacity, simulated flow over two configurations was produced for Ht/P = 0.3. Figure 10 presents the flow streamlines for both flat and rounded entrances. It is clearly shown that the less flow contraction is provided by the model with the rounded entrance (Figure 10(b)) compared to the model with the flat entrance (Figure 10(a)). The use of the rounded entrance, providing a smoothly flow contraction, increases the rectangular labyrinth weir capacity by decreasing the energy losses at the inlet. Streamlines behavior at the inlet of the weir is presented in Figure 11(a) and 11(b) for flat and rounded shapes, respectively. The distribution of streamlines over the crest of two configurations is also presented in Figure 12. It was shown that the streamlines are distributed along the crest of the labyrinth designed with the rounded entrance (Figure 12(b)). However, the streamlines over the labyrinth with the flat entrance are mainly concentrated at the end portion of the crest as shown in Figure 12(a).

Figure 10

3D flow and streamlines over labyrinth weir, (a) with flat entrance and (b) with rounded entrance.

Figure 10

3D flow and streamlines over labyrinth weir, (a) with flat entrance and (b) with rounded entrance.

Close modal
Figure 11

Plan flow of streamlines over labyrinth weir, (a) with flat entrance and (b) with rounded entrance.

Figure 11

Plan flow of streamlines over labyrinth weir, (a) with flat entrance and (b) with rounded entrance.

Close modal
Figure 12

Streamlines distribution over the crest of labyrinth weir (a) flat entrance, (b) rounded entrance.

Figure 12

Streamlines distribution over the crest of labyrinth weir (a) flat entrance, (b) rounded entrance.

Close modal

Alveoli width effect

Figure 13 compares the numerically computed coefficient of discharge for three rectangular labyrinth weirs with varying ratios (a/b = 0.67, 1 and 1.5). It was observed that the numerical model enables prediction of the influence of the alveoli width variations and shows that the model with the a/b ratio = 1.5 is more efficient compared with the model with 0.65 and provides a slightly larger capacity than the model with ratio a/b = 1. Globally, this observation is in agreement with the observation provided by the experimental model (Figure 6(b)). Streamlines over the rectangular labyrinth weir of three models have been numerically generated for the same upstream head Ht/P = 0.3 (Figure 14). A high concentration of streamlines can be clearly seen within the inlet alveoli of three simulated models, indicating that the majority of flow passes through the inlet alveoli.

Figure 13

Discharge coefficient obtained from numerical simulation for the alveoli width.

Figure 13

Discharge coefficient obtained from numerical simulation for the alveoli width.

Close modal
Figure 14

Streamlines over the three tested models (a) a/b = 0.67, (b) a/b = 1 and (c) a/b = 1.5.

Figure 14

Streamlines over the three tested models (a) a/b = 0.67, (b) a/b = 1 and (c) a/b = 1.5.

Close modal

A 3D flow also presented for three models in Figure 15, this highlights that the inlet alveoli of the model with a/b = 0.67 is poorly supplied by the flow compared to the models with ratios a/b = 1 and 1.5. These observations confirmed that the use of inlet width larger than the outlet width leads to an increase in the discharge capacity.

Figure 15

Front and plan view of simulated flow over rectangular labyrinth weir with a ratio a/b equal 0.67, 1, and 1.5.

Figure 15

Front and plan view of simulated flow over rectangular labyrinth weir with a ratio a/b equal 0.67, 1, and 1.5.

Close modal

As previously discussed, the efficiency of the rectangular labyrinth weir designed with rounded entrance is larger than that obtained from the rectangular labyrinth weir having a flat entrance. The rounded shape enables an increase in the efficiency by 5% compared to the flat entrance. This can be attributed to the diminution in loss of energy at the entrance weir. Similar results were found by the experimental work carried out in reservoir application by Ouamane (2014) on two trapezoidal labyrinth weirs, one with a flat entrance and the second with a rounded entrance.

They confirmed that the trapezoidal labyrinth weir with rounded entrance provides higher discharge compared to the flat entrance. Their results indicate that the rounded shape provides a gain in efficiency of around 10%. Two reasons for this case were reported. First, adopting a rounded shape at the entrance facilitates the entrance of the flow in the inlet key, and thus the weir efficiency. Second, the rounded shape of the entrance enables elimination of the discontinuity points over the crest and thus to increase the effective crest length. Ouamane & Lempérière (2006) stated that installing a rounded or triangular nose under the upstream apex overhangs the Piano Key Weir, a new type of labyrinth weir, which increases the discharge capacity by 7%. This technic was used for designing Goulours and Saint-Marc dams (Laugier 2007, 2009). Anderson & Tullis (2013) also conducted tests to study the effect of the noses on the discharge capacity of Piano Key Weir. They indicated that noses with round shapes installed beneath the Piano Key Weir upstream apex overhangs increased the discharge capacity by up to 2.8%. They explain that the increase in efficiency by a decrease in inlet energy loss associated with the improved flow conditions at the entrance of the inlet.

The rectangular labyrinth weir designed with a larger inlet alveoli width produced the highest discharge capacity, as shown in the results of this study. The majority of the reported studies have been conducted to evaluate the effect of the alveoli width variation on the Piano Key Weir (Ouamane & Lempérière (2006), Machiels et al. (2011) and Andersonm & Tullis (2013). However, to the author's knowledge, the studies by Ouamane & Ben Said (2010) and Kabiri-Samani et al. (2013) are the only two studies carried out to investigate the influence of the alveoli width on the capacity of the rectangular labyrinth weir. These two studies were performed in two different cases of approach flow conditions. In the first study, the labyrinth weir designed with a rounded entrance was installed in the reservoir, and in the second study, the model with the flat entrance was tested in channel conditions. Both studies reported that the discharge capacity increased with increasing the alveoli width ratio a/b. For low upstream heads Ht/P, the rectangular labyrinth weir leads to an increase in the efficiency by up to 10%. However, it is limited to 2% for higher upstream heads (Ouamane & Ben Said 2010). Kabiri-Samani et al. (2013) stated that a value of the ratio of the upstream apex length to the total spillway width (b/B) equal or greater than 0.1, for Ht/P was a sufficient value to limit the effect of the lateral nappe interference on the labyrinth efficiency. The highest efficiency in the current study corresponded to the ratio a/b equal to 1.5, which is within the optimal range, a/b = 1.25–1.5, as recommended by Andersonm & Tullis (2013) for designing the Piano Key Weir. The gain in efficiency of the model with a/b = 1.5 is related to the increase in the inlet alveoli capacity due to the increase in the inlet cross-section, which allows the entrance of the greatest flow in the inlet alveoli and thus increases the weir efficiency.

Machiels et al. (2014) attest that increasing the inlet cross-section caused a decrease in the flow velocity over the lateral crest, and thus it increased the contribution of the lateral crest in the overall efficiency of the Piano Key Weir.

The purpose of this study was to investigate the influence of the geometric parameters on the rectangular labyrinth weir capacity. This study was conducted using an experimental and numerical model to investigate the entrance and the inlet and outlet width ratio a/b effect on the rectangular labyrinth weir efficiency. All models of the rectangular labyrinth weir have been tested in channel flow conditions for the same developed crest length ratio L/W = 4. The main conclusions of this study may be drawn as follows:

  • The weir entrance shape significantly affects the discharge efficiency. A rectangular labyrinth weir designed with a rounded entrance was found to provide the highest discharge capacity. The gain in discharge efficiency is close to 5% compared to the labyrinth designed with flat entrance shape, this is for the whole range of upstream heads considered in this study (Ht/P≤1).

  • For the rectangular labyrinth weir with a ratio L/W = 4 and designed with a rounded entrance, the results indicate that the increase of the ratio between the inlet and outlet alveoli width enables an increase in the discharge capacity of the labyrinth. The ratio a/b = 1.5 led to a 18% increase in the efficiency compared to the model with a/b = 0.67 and it is around 8% relative to the model with a/b = 1, mainly for upstream heads (Ht/P≤0.5).

  • Using an inlet width larger than the outlet width enables an increase in the efficiency of the rounded rectangular labyrinth weir at no additional cost.

  • Generally, a good agreement was found between the physical and numerical results. The obtained results indicated that the numerical model is able to compute the discharge capacity of the tested models of the rectangular labyrinth weir with mean average relative error of 4 and 6% for a/b = 1 (flat entrance) and a/b = 0.67, respectively.

  • The numerical simulation using OpenFOAM CFD helps to further analyze the rectangular labyrinth weir efficiency by providing important information about the flow behavior, such as free surface profile and streamlines.

  • The CFD simulations can be a useful tool to analyze and pre-design the labyrinth weir and help to reduce the cost of future weir designs.

All relevant data are included in the paper or its Supplementary Information.

Andersonm
R. M.
&
Tullis
B. P.
2013
Piano key weir hydraulics and labyrinth weir comparison
.
J. Hydraul. Eng.
139
(
3
),
246
253
.
Aydin
M. C.
&
Emin Emiroglu
M.
2013
Determination of capacity of labyrinth side weir by CFD
.
Flow Meas. Instrum.
29
,
1
8
.
Azimi
A. H.
&
Seyed Hakim
S.
2019
Hydraulics of flow over rectangular labyrinth weirs
.
Irrig. Sci.
37
(
2
),
183
193
.
Beg
M.
,
Carvalho
R. F.
,
Tait
S.
,
Brevis
W.
,
Rubinato
M.
,
Schellart
A.
&
Leandro
J.
2017
A comparative study of manhole hydraulics using stereoscopic PIV and different RANS models
.
Water Sci. Technol.
2017
(
1
),
87
98
.
Ben Said
M.
&
Ouamane
A.
2011
Study of the optimization of labyrinth weirs
. In:
Proceedings of the International Conference on Labyrinth and Piano Key Weirs – PKW 2011
,
Liège, Belgium
, pp.
67
74
.
Ben Said
M.
&
Ouamane
A.
2018
Remodelling labyrinth weir geometry to increase performance
.
Intl. J. Hydropower Dams
10
(
5
),
74
78
.
Carrillo
J. M.
,
Matos
J.
&
Lopes
R.
2019
Numerical modeling of free and submerged labyrinth weir for a large sidewall angle
.
Environ. Fluid Mech.
20
,
357
374
.
Celik
I. B.
,
Ghia
U.
,
Roache
P. J.
,
Freitas
C. J.
,
Coleman
H.
&
Raad
P. E.
2008
Procedure for estimation and reporting of uncertainty due to discretization in CFD applications
.
ASME J. Fluids Eng.
130
,
1
4
.
Crookston
B. M.
&
Tullis
B. P.
2012
Discharge efficiency of reservoir application specific labyrinth weirs
.
J. Irrig. Drain. Engr, ASCE
138
(
6
),
564
568
.
Dabling
M. R.
,
Tullis
B. P.
&
Crookston
B. M.
2013
Staged labyrinth weir hydraulics
.
J. Irrig. Drain. Engr, ASCE
139
(
11
),
955
960
.
Darvas
L. A.
1971
Discussion of performance and design of labyrinth weirs, by Hay and Taylor
.
J. Hydraul. Eng.
97
(
80
),
1246
1251
.
Duguay
J. M.
,
Lacey
R. W. J.
&
Gaucher
J.
2017
A case study of a pool and weir fishway modeled with OpenFOAM and FLOW-3D
.
Ecological Engineering
103
(
A
),
31
42
.
Falvey
H. T.
2003
Hydraulic Design of Labyrinth Weirs
.
ASCE Press
,
Reston
,
USA
.
Gentilini
B.
1940
Weirs with an Oblique Plant and Zigzag Crest
.
Memories and Studies of the Institute of Hydraulics and Hydraulic Constructions of the Royal Polytechnic of Milan
.
No. 48
.
Ghaderi
A.
,
Daneshfaraz
R.
,
Dasineh
M.
&
Di Francesco
S.
2020
Energy dissipation and hydraulics of flow over trapezoidal–triangular labyrinth weirs
.
Water
12
(
7
),
1992
.
Hay
N.
&
Taylor
G.
1970
Performance and design of labyrinth weirs
.
J. Hyd. Div. ASCE
96
(
2
),
2337
2357
.
Hinchliff
D. L.
&
Houston
K. L.
1984
Hydraulic design and application of labyrinth spillways
. In:
Proceedings of USCOLD Lecture Dam Safety and Rehabilitation, USCOLD
,
Denver, Colorado. USA
.
Houston
K. L.
1983
Hydraulic Model Study of Hyrum Auxiliary Labyrinth Spillway
.
Rep. No. GR-82-/3
,
Bureau of Reclamation
,
Denver, CO, USA
.
Hydrocoop
2013
New Labyrinth Weirs triple the spillways discharge, Posted on July 18, 2013 in Piano Keys Weirs. http://www.hydrocoop.org/new-labyrinth-weirs-triple-the-spillways-discharge/.
Kabiri-Samani
A.
,
Javaheri
A.
&
Borghei
S. M.
2013
Discharge coefficient of a rectangular labyrinth weir
. In:
Proceedings of the Institution of Civil Engineers, Water Management
166
(
8
),
443
451
.
Laugier
F.
2007
Design and construction of the first piano key weir spillway at Goulours dam
.
Int. J. Hydropower Dams.
14
(
5
),
94
100
.
Laugier
F.
2009
Design and construction of a labyrinth PK Weir spillway at Saint-Marc dam, France
.
Int. J. Hydropower Dams.
15
(
5
),
100
107
.
Lempérière
F.
&
Ouamane
A.
2003
The piano keys weir: a new cost-effective solution for spillways
.
Int. J. Hydropower Dams
10
(
5
),
144
149
.
Lopes
R.
,
Matos
J.
&
Melo
J. F.
2009
Discharge capacity for free-flow and submerged labyrinth weirs
. In:
Proceeding 33rd IAHR Congress, Water Engineering for A Sustainable 368 Environment
,
Vancouver, Canada
.
Lux
F.
&
Hinchliff
D.
1985
Design and construction of labyrinth spillways, Q95
. In:
15th ICOLD Congress
,
Lausanne, Switzerland
.
Machiels
O.
,
Erepicum
S.
,
Dewals
B. J.
,
Archambeau
P.
&
Pirotton
M.
2011
Experimental observations of flow characteristics over a piano key weir
.
J. Hydraul. Res.
49
(
3
),
359
366
.
Machiels
O.
,
Pirotton
M.
,
Archambeau
P.
,
Dewals
B. J.
&
Erepicum
S.
2014
Experimental parametric study and design of Piano Key Weirs
.
J. Hydraul. Res.
52
(
3
),
326
335
.
Moukalled
F.
,
Mangani
L.
&
Darwish
M.
2015
The Finite Volume Method in Computational Fluid Dynamics
.
Springer International Publishing
,
Switzerland
.
Murphy
D. W. A.
1909
Reinforced-concrete spillway with concentrated crest length
.
Eng. News
62
(
11
),
278
279
.
Ouamane
A.
2014
Improvement of labyrinth weirs shape. Proc. Int. Work. Labyrinth and Piano Key Weirs-PKW 2013, Paris, S. Erpicum, F. Laugier, M. Pfister, M. pp. 15-22
.
Ouamane
A.
&
Lempérière
F.
2006
Design of a new economic shape of weir
. In: Berga, L., Buil, J. M., Bofill, E. and De Cea, J. C. (eds)
Dams and reservoirs, societies and environment in the 21st century
.
Taylor & Francis Group
,
London, ISBN 0 415 404231
.
Ouamane
A.
&
Ben said
M.
2010
Economic Conceptions of the Labyrinth Weir Allowing the Improvement of the Management of Floods and the Increase of the Useful Capacity of the Reservoir
.
Rheinisch-Westfälische Technische Hochschule Aachen, Publication No. 158
,
Published by Shaker Verlag
,
Germany
.
Safarzadeh
A.
&
Noroozi
B.
2017
3D hydrodynamics of trapezoidal piano key spillways
.
Int. J Civ. Eng.
15
,
89
101
.
Savage
B.
,
Crookston
B. M.
&
Paxson
G.
2016
Physical and numerical modeling of large headwater ratios for a 15° labyrinth spillway
.
J. Hydrol. Eng.
142
(
11
),
04016046
.
Taylor
G.
1968
The Performance of Labyrinth Weirs
.
PhD thesis
,
University of Nottingham
,
UK
.
Tullis
J. P.
,
Amanian
N.
&
Waldron
D.
1995
Design of labyrinth spillways
.
J. Hydraul. Eng.
121
(
3
),
247
255
.
Tullis
B. P.
,
Young
J.
&
Chandler
M.
2007
Head-discharge relationships for submerged labyrinth weirs
.
J. Hydraul. Eng.
133
(
3
),
248
254
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).