Abstract
Large orifices are constructed for dams to release water and sediments from reservoirs. Such structures are called submerged spillways. Numerous studies have investigated discharge coefficient, velocity coefficient, and head loss coefficient of large orifices; however, the literature lacks data on the upper and lower nappes of the jets from these orifices. In the present experimental study, the upper and lower nappes are investigated up to 80 m head at different gate openings. The observed minor deviation between the lower nappe profile and trajectory profile equation suggests sensitivity to different factors. The significant role of the coefficient of velocity, averaging at 0.926, highlights its impact on minor deviation. Subsequently, the impact of the solid bottom profile on the discharge coefficient and upper nappe profile are also examined. The results show improvement in discharge coefficient of a sharp-edged large orifice, which increased from 0.69 to 0.74. The results also indicate that the upper nappe profiles and United States Bureau of Reclamation (USBR) profiles are similar. The improvement in the upper nappe profile indicates the significant role of the solid bottom profile, which consequently was found to be helpful in defining the roof profile of an orifice spillway. .
HIGHLIGHTS
Hydraulic analyses of lower and upper nappe profiles of sharp-edged large orifices.
Impact evaluation of a solid bottom profile on an upper nappe profile.
Hydraulic performance evaluation of solid bottom profile at various operational scenarios.
NOTATIONS
- D
gate opening/height of orifice
- b
width of orifice
- Re
Reynolds number
- Fr
Froude number
- We
Weber number
- H
operating head
- Cd
discharge coefficient
- Qth
theoretical discharge
- h
water level in stilling well of V-notch.
- hc
center line head over the orifice
- Cv
coefficient of velocity
- x, y
horizontal and vertical coordinates of lower nappe profile
- x1, y1
horizontal and vertical coordinates of upper nappe profile
INTRODUCTION
An orifice is an opening constructed in a dam structure to regulate the flow. The large rectangular orifices are generally constructed to pass the surplus flows and are used to dispose of the river sediments. Velocities more than 25 m/s over the spillway crest and discharge intensities of 200–300 m3/s/m make the design of the orifice spillway both complex and involved, requiring detailed investigation during planning besides studies on a scaled model (Bhosekar et al. 2014). However, the operation of these spillways also creates complicated issues, such as turbulence, sediment transport (Pu et al. 2014, 2016), natural bed forms (Pu et al. 2017), and cavitation (Sarwar et al. 2016). During spillway operation, energy dissipation of flowing water depends upon the Froude Number (Fr. No.). The higher the Fr. No., the greater the dissipation of energy. In vertical shafts, flow energy dissipation efficiency (ηs) decreases with increasing Fr. No. It was found to be between 10.80 and 62.29% (Mohmoudi-Rad & Najafzadehhe 2022, 2023). The size of an orifice is influenced by several factors such as water head, which is classified into three different categories i.e., Low Head, Medium Head, and High Head Orifice. A need for a physical and numerical model for modelling orifice spillway flows was identified. The discharge coefficient (Cd) was found in the range of 0.831–0.942 (Bos 1989; Swamee et al. 1998; Gadge et al. 2019; Mozaffari et al. 2022). Additionally, the actual size of the orifice may vary depending on the specific design requirements of the dam while other parameters that affect the size of the orifice are flow rates, size of reservoir, and expected floods (Amirkhani et al. 2017; Bhattarai & Sharma 2017).
Numerous studies investigated hydraulic coefficients for the flow through an orifice, but coefficient of discharge, coefficient of velocity, and head loss coefficient (Essien et al. 2019; Hussain et al. 2014) are the most significant ones and play vital roles to control, measure, and regulate the flow in the design and operation of orifices. Essien et al. 2019 conducted an experimental investigation on Cd. An empirical correlation between discharge coefficient and size of orifice was developed. Cd values obtained from experimental data and those from the empirical correlation were compared, and a mean standard deviation of 0.0231 was obtained. In addition, these parameters are also used for prediction of flow and energy loss (Babu et al. 2018; Manzano-Miura et al. 2022). Roof and bottom profiles of orifice spillways are important parts of orifice spillway geometry. These profiles are designed to optimize water discharge and control during spillway operation. The roof profile of an orifice spillway refers to the shape or contour of the structure's top surface, which directs water flow. Meanwhile, the bottom profile refers to the shape or contour of the spillway's lower surface, influencing how water exits the spillway. The design of roof profile and bottom profile has not been standardized (Khatsuria 2013). Vatankhah & Rafeifar (2020) conducted a series of laboratory runs (588 runs) for different values of orifice geometry. Using measurements obtained by laboratory runs, the proposed theoretical models of elliptical side orifices were calibrated under free outflow conditions. The model that includes the approach Froude number had an average error of 1.74%, while the other model that does not include the approach Froude number had an average error of about 2.43%. Haghbin & Sharafati (2022) provided a comprehensive review of the application of soft computing (SC) models for estimating Cd of different flow control structures such as ogee spillways, orifices, and side weirs. Besides, the ratio of orifice crest height to height of side orifice, the ratio of main channel width to length of side orifice, the ratio of main channel width to height of side orifice, and the ratio of the height of the side orifice to upstream flow depth were extensively employed to calculate Cd of orifice structures. Maher et al. (2019) investigated the effect of different inlet configurations on the discharge coefficient of large orifices and found that the discharge coefficient was highest for the circular inlet configuration. A study conducted by Bhattarai & Sharma (2017) noted a significant impact of orifice size on the spillway discharge capacity of the embankment dam. Conversely, Jothiprakash et al. (2015) noted that the velocity coefficient was affected by the upstream flow condition and the shape of the orifice. Similarly, Jithish & Ajay Kumar (2015) used Computational Fluid Dynamic (CFD) simulations to find out the effect of orifice size, its inlet and outlet configuration, and Reynolds number on the velocity coefficient. Hussain et al. (2016) developed relationships for coefficient of discharge for orifices under free and submerged flow conditions through analytical and experimental considerations. The computed discharges using developed relationships were within ±5% and ±10% of the observed ones for free and submerged orifices, respectively. A sensitivity analysis revealed that the discharge through the side orifice is more sensitive to the low head above the centre of the orifice. Apart from free flow and submerged flow conditions for orifice structure, Adam et al. (2016), Mali et al. (2020), Adam et al. (2019), and Wang et al. (2020) studied head loss coefficient for the large orifices. These studies focused on orifice size, submergence ratio, and inlet velocity on the head loss coefficient.
From the bibliographic analysis, due to the significance of spillways to discharge the floods safely, in addition to the discharge coefficient, velocity coefficient, and head loss, the study of spillway profiles, i.e., roof and bottom profiles are crucial in the design and operation of spillways. However, except for Gadge et al. (2019) who developed design guidelines and equations for the bottom and roof profiles (applicable to orifice spillways with head varying from 30 to 70 m), the literature lacks information on the design of roof profile and bottom profile for 40–80 m heads and 10–16 m orifice openings. Therefore, in the present study, a physical model is developed and operated under different water heads (40–80 m) and gate openings (10–16 m). This research study is focused on orifice spillways by investigating the flows through a sharp-edged large orifice. Understanding of lower and upper nappe profiles of the jet of large orifices may be helpful in determining the bottom and roof profiles of an orifice spillway.
MATERIAL AND METHODS
The present study has developed a physical model of 1:100 scale to investigate the roof and bottom profiles of a large orifice spillway. Before developing the physical model, dimensional analysis is carried out to identify the critical hydraulic parameters that affect the roof and bottom profiles of the studied spillway. The details of dimensionless numbers, physical model, and measurement techniques are explained in the proceeding sections 2.1, 2.2, and 2.3, respectively.
Scale modelling of flows through large orifice
The viscous and surface tension forces have a negligible effect on free surface flows whereas gravitational force has a considerable impact on free surface flows (Rezazadeh et al. 2020, Mahtabi & Arvanaghi 2018, and Daneshfaraz et al. 2023). Therefore, at present, Froude number similarities are employed to develop the physical model, whereas to minimize the viscous and surface tension forces, Reynolds number (Re > 105) and flow depth over the crest against the design head are employed for deciding the scale of the physical model (USBR 1980; Pfister & Chanson 2014; Pfister & Hager 2014; Gadge et al. 2016).
Development of physical model
The present physical model was developed in the hydraulic laboratory of the Centre of Excellence in Water Resources Engineering (CEWRE), Lahore, Pakistan. A scaled model of 1:100 ratio was constructed and placed in a 1-m wide and 15-m long flume.
Lower and upper nappe profile measurements
Discharge measurement
RESULT AND DISCUSSION
Coefficient of discharge for sharp-edged orifice
From Table 1 the Cd for a sharp-edged large orifice was found to range between 0.30 and 0.66, which agree with the values of Bansal 2010 and Prajakta et al. 2016. The variation in Cd values is the result of different factors such as head range, the size of the orifice, approach velocity, and orifice design and edge sharpness. In Table 1, H is the head above the crest of the orifice.
Gate opening . | Cd . | Average Value . | ||||
---|---|---|---|---|---|---|
H-40 m . | H-50 m . | H-60 m . | H-70 m . | H-80 m . | ||
10 m | 0.30 | 0.34 | 0.33 | 0.32 | 0.32 | 0.32 |
12 m | 0.42 | 0.42 | 0.38 | 0.39 | 0.37 | 0.40 |
14 m | 0.44 | 0.42 | 0.44 | 0.46 | 0.43 | 0.44 |
16 m | 0.54 | 0.55 | 0.55 | 0.63 | 0.66 | 0.59 |
Gate opening . | Cd . | Average Value . | ||||
---|---|---|---|---|---|---|
H-40 m . | H-50 m . | H-60 m . | H-70 m . | H-80 m . | ||
10 m | 0.30 | 0.34 | 0.33 | 0.32 | 0.32 | 0.32 |
12 m | 0.42 | 0.42 | 0.38 | 0.39 | 0.37 | 0.40 |
14 m | 0.44 | 0.42 | 0.44 | 0.46 | 0.43 | 0.44 |
16 m | 0.54 | 0.55 | 0.55 | 0.63 | 0.66 | 0.59 |
Nappe profiles of orifice flows
To develop the solid bottom and roof profiles for orifice spillways (Figure 2), upper and lower nappe profiles are drawn. These profiles are then compared with profiles obtained through the empirical relations 6 and 7.
Lower nappe profiles for sharp-edged large orifice
From the experiments, the values of the Cv were found to be ranged between 0.894 and 0.947 with an averaged value of 0.926. The variation in the Cv value was due to the change in k values. The results further indicated that any change in the values of Cv and k directly affected the lower profile of the orifice. For different heads and orifice gate openings, values of k and Cv are presented in Table 2.
Gate Opening (m) . | Head (m) . | Coefficient of velocity (Cv) . | Constant (k) . |
---|---|---|---|
10 | 40 | 0.894 | 3.278 |
50 | 0.904 | 3.446 | |
60 | 0.912 | 3.386 | |
70 | 0.925 | 3.468 | |
80 | 0.933 | 3.606 | |
16 | 40 | 0.890 | 3.497 |
50 | 0.909 | 3.373 | |
60 | 0.918 | 3.407 | |
70 | 0.930 | 3.513 | |
80 | 0.947 | 3.630 |
Gate Opening (m) . | Head (m) . | Coefficient of velocity (Cv) . | Constant (k) . |
---|---|---|---|
10 | 40 | 0.894 | 3.278 |
50 | 0.904 | 3.446 | |
60 | 0.912 | 3.386 | |
70 | 0.925 | 3.468 | |
80 | 0.933 | 3.606 | |
16 | 40 | 0.890 | 3.497 |
50 | 0.909 | 3.373 | |
60 | 0.918 | 3.407 | |
70 | 0.930 | 3.513 | |
80 | 0.947 | 3.630 |
From the results, it is found that the k value is ranged between 3.27 and 3.630, while the literature has indicated a maximum value of 4. In comparison to the conventional value of 4, the lower value of k significantly affected the curvature of the profile, which further resulted in a steeper curve. After employing the k values, the trends of the bottom profile of the orifice spillway are found to be in good agreement with the empirical data ().
Upper nappe profile for sharp-edged large orifice
Considering this significant difference in the upper nappe profiles of the present study and the existing literature representing the USBR approach, the observed upper nappe profiles in this study cannot be directly applicable for determining the roof profile of an orifice spillway. Therefore, in the following section, the impact of the solid bottom profile on the upper nappe profile in the form of experimental results is presented.
Significance of solid bottom profile for upper nappe of large orifice
Impact of solid bottom profile on coefficient of discharge
The solid bottom profile as shown in Figure 2 plays a crucial role in influencing the discharge coefficient of the orifice. Table 3 shows an improvement in the discharge coefficient on installation of the solid bottom profile as compared with the results of Table 1. The discharge coefficient varies from 0.39 to 0.74. This range of discharge coefficients is more than that without a solid bottom. The solid bottom profile was developed by analysing lower nappe profiles and using the equation . It helped in maintaining a consistent pressure distribution, reducing losses, and ultimately improving the discharge coefficient.
Gate opening . | Coefficient of Discharge . | Average value . | ||||
---|---|---|---|---|---|---|
H-40 m . | H-50 m . | H-60 m . | H-70 m . | H-80 m . | ||
10 m | 0.39 | 0.40 | 0.43 | 0.44 | 0.46 | 0.42 |
16 m | 0.63 | 0.64 | 0.66 | 0.69 | 0.74 | 0.67 |
Gate opening . | Coefficient of Discharge . | Average value . | ||||
---|---|---|---|---|---|---|
H-40 m . | H-50 m . | H-60 m . | H-70 m . | H-80 m . | ||
10 m | 0.39 | 0.40 | 0.43 | 0.44 | 0.46 | 0.42 |
16 m | 0.63 | 0.64 | 0.66 | 0.69 | 0.74 | 0.67 |
Pressure and velocity distribution along a solid bottom profile
At spillway chutes, the understanding of both positive and negative pressures is crucial from design and safety considerations. A balanced positive pressure ensures efficient water discharge. Conversely, negative pressures lead to cavitation. In the current study, the pressure head remained positive along the solid bottom profile, except at the last observation point as shown in Figure 7.
CONCLUSIONS
A physical model study was conducted to hydraulically investigate the flows through sharp-edged large orifices. The model was built at a 1:100 scale by considering discharge and space constraints. It was operated at various gate openings (10–16 m, in 2 m increments) and water heads (40–80 m, in 10 m increments). Upper and lower nappe profiles of a water jet and discharge coefficients were observed. A solid bottom profile (chute) was designed considering the lower nappe profiles. After installation of the chute, the upper nappe profiles, discharge coefficients, pressures, and velocities were measured under diverse operational scenarios.
The comparison between upper and lower profiles and the results of the trajectory profile equation revealed good agreement, underscoring the reliability of the model. However, small variations in trajectories signify the sensitivity of the profiles to different factors. Particularly, the Cv, averaging at 0.926, emerged as a critical factor. The variation in the coefficient of velocity significantly impacts the coefficient k of the trajectory profile equation. The maximum k value in this study was 3.63, slightly below the conventional threshold of 4. This deviation implies a steep curvature of the solid bottom profile (chute). Consequently, the curvature of the chute was designed using k value at 4.
The flows on the chute exhibited very high velocities (v > 35 m/s) and negative pressures near the discharging end, particularly at heads exceeding 40 m. This suggests that a k value exceeding 4 could yield a more suitable curvature for the 16 m orifice opening at heads above 40 m. In the absence of achieving k > 4, there's a recommendation to install an offset type of aerator to prevent cavitation damages, emphasizing the critical role of trajectory profile parameters in ensuring optimal performance of orifice spillways.
Upper nappe profiles clearly deviated from the USBR profiles at all models' operational scenarios before installation of the solid bottom profiles. The installation of solid bottom profiles significantly improved the alignment between upper nappe profiles and USBR profiles. This improvement in alignment highlights the importance of a solid bottom profile in determining upper nappe profiles. Despite the alignment, discrepancies persist, as the USBR equation lacks consideration for head variation. The solid bottom profile also contributed to enhancing the discharge coefficient, with a notable improvement from 0.66 to 0.74 at an 80-m head. These findings are crucial for refining roof profiles and optimizing discharge coefficients in orifice spillway design. Therefore, it is required to investigate the upper nappe profile for different combinations of heads and orifice openings to standardize the roof profile of orifice spillways.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.