Discharge Coefficient of Rectangular Short-Crested Weir with Varying Slope Coefficients

Rectangular short-crested weirs are widely used for simple structure and high discharge capacity. As one of the most important and influential factors of discharge capacity, side slope can improve the hydraulic characteristics of weirs at special conditions. In order to systemically study the effects of upstream and downstream slope coefficients S1 and S2 on overflow discharge coefficient in a rectangular short-crested weir the Volume of Fluid (VOF) method and the Renormalization Group (RNG) κ-ε turbulence model are used. In this study, the slope coefficient ranges from V to 3H:1V and each model corresponds to five total energy heads of H0 ranging from 8.0 to 24.0 cm. Comparisons of discharge coefficients and free surface profiles between simulated and laboratory results display a good agreement. The simulated results show that the difference of discharge coefficients will decrease with upstream slopes and increase with downstream slopes as H0 increases. For a given H0, the discharge coefficient has a convex parabolic relation with S1 and a piecewise linearity relation with S2. The maximum discharge coefficient is always obtained at S2 = 0.8. There exists a difference between upstream and downstream slope coefficients in the influence range of free surface curvatures. Furthermore, a proposed discharge coefficient equation by nonlinear regression is a function of upstream and downstream slope coefficients.


Introduction
As common hydraulic engineering structures, low weirs are widely used to measure discharge in irrigation systems and to increase upstream water level in hydroelectric projects.Flat-topped weirs of finite crest length are classified into four types: long-crested weir (0 < H/δ ≤ 0.1), broad-crested weir (0.1 ≤ H/δ ≤ 0.4), short-crested weir (0.4 ≤ H/δ ≤ 1.5-1.9),and sharp-crested weir (1.5-1.9 ≤ H/δ), depending on the relative length of crest H/δ, where H is the crest depth over weir at free overflow condition and δ is the length of weir crest in the streamwise direction [1].The flat-topped low weirs generally include broad-crested weir and short-crested weir.Under the same inflow condition, the discharge coefficient of short-crested weir is approximately 0.33-0.46,while that of broad-crested weir is 0.32-0.385,hence the former is stronger than the latter in terms of discharge capacity.The short-crested weirs are designed into curvilinear profile and broken-line profile according to the longitudinal profile.The former is usually used as a high weir and the latter as a low weir.Considering as the small overflow weirs, although the discharge coefficient of former is larger, the latter is widely used in view of the convenience of construction.The broken-line short-crested weir is classified into rectangular profile and trapezoidal profile according to the longitudinal profile [2].Moreover, the rectangular short-crested weir is insensitive to the downstream submergence, and in special situations the geometrical configuration could be flexibly modified so as to increase discharge capacity and avoid cavitation damage of weir flow structures [3].
Extensive experimental and numerical studies on weir height and weir crest length of the short-crested weir were carried out [4][5][6][7], and the corresponding discharge coefficient formulae have been derived [1,[8][9][10].Although the existing research displays a close relation between slope coefficient and discharge coefficient of broad-crested weirs [10][11][12], recently there has been little research on the effect of slope coefficient on short-crested weirs [4,6].Azimi et al. [13] found that flow separated from the entrance of crest, and this part of flow would reattach to the weir crest with curvilinear flow existing over the entire weir before it left the crest again.Goodarzi et al. [14] showed that the size of flow separation space at the entrance of weir crest varied with upstream slope, which indirectly resulted in variation of discharge coefficients.Bos [4] found that a higher streamline curvature above the weir crest and a smaller size on separation region had a positive influence on the head-discharge relationship of weirs.Farhoudi et al. [15] concluded that decreasing upstream slope angle increased the discharge coefficient until it reached the maximum value at α = 25 • (2.1H:1V).Compared with a rectangular weir with a sloped upstream weir face, Sargison et al. [16] found that the weir with sloped downstream weir face could improve the discharge capacity.Recently, Tong et al. [10,17] have carried out much work on investigating the flow characteristics of short-crested weir and concluded that the maximum discharge coefficient existed at S 1 = 2.0-3.0.Chen et al. [18,19] analyzed the effect of upstream slope on the discharge coefficient by combining experiments and numerical simulations.Li et al. [20] pointed that the discharge coefficient of rectangular short-crested weir with S 1 = 0.5 is 2.0% larger than that of the rectangular short-crested weir, and the flow regime of the former was better than the latter.Haun et al. [21] applied computational fluid dynamics (CFD) software Flow-3D [22] to simulate the free overflow over the trapezoidal broad-crested weir, and the results were in good agreement with the experimental data [16].Paik et al. [23] used the RNG κ-ε turbulent model to calculate the two flow separation zones for free overflow over the rectangular broad-crested weir, which fit well with those of physical models.In conclusion, the approach of numerical simulation meets the requirement of studying the hydraulic characteristics of the free surface weir flow.
Almost all of the above conclusions are premised on the inflow discharge, whereas this study is premised on the total energy head H 0 over the crest to realize the flexible manipulation in engineering.Besides, the previous research on the discharge coefficient has not taken the slope coefficients into consideration [1,10].In order to meet the demand of high accuracy on applicability and construction design, performing research on the effect of slope coefficients on discharge capacity is necessary.The present paper uses the Flow-3D software to simulate the free overflow over the rectangular short-crested weir with varying upstream and downstream slope coefficients, and systematically studies the effects of slope coefficients on the discharge capacity of short-crested weir.Moreover, a calculation formula of discharge coefficient including the variables of upstream and downstream slope coefficients is derived by the nonlinear regression method, which provides references for the engineering design of short-crested weir.

Theoretical Analysis of Influential Factors of Discharge Coefficient
Figure 1 showed the definition sketch of short-crested weir constructed from three parts in the longitudinal cross-section: the upstream weir face U, the rectangular crest R and the downstream weir face D. The axis of x and y were along the longitudinal direction and vertical direction, respectively.The weir width B was equal to the flume width without considering lateral contraction.The common practice to determine discharge coefficient in most discharge measuring structures is based on the method of dimensional analysis.The discharge coefficient of the short-crested weir depends on hydraulic conditions, geometrical parameters and fluid properties [24].Considering these influential factors one could end up with the following relation: Water 2018, 10, 204 where Q is inflow discharge, m 3 /s; H 0 is total energy head over crest, H 0 = H + U 2 /2g, m; U is the approaching velocity, m/s; B is the width of weir, m; δ is the length of weir crest, m; P is the height of weir, m; S 1 and S 2 are up-and downstream slope coefficients, respectively; g is gravity acceleration, m/s 2 ; σ is surface tension of fluid, N/m; ρ is mass density of fluid, kg/m 3 ; µ is dynamic viscosity of fluid, Pa•s.Based on Buckingham's Π theorem, the three variables of H 0 , g and ρ are selected as the basic dimensions, and the relations between dimensionless numbers could be expressed as follows: = f δ/H 0 , P/H 0 , S 1 , S 2 , ρ 0.5 g 0.5 H 0 /σ 0.5 , ρg 0.5 H 0 1.5 /µ (2) Water 2018, 10, x FOR PEER REVIEW 3 of 18 where Q is inflow discharge, m 3 /s; H0 is total energy head over crest, H0 = H + U 2 /2g, m; U is the approaching velocity, m/s; B is the width of weir, m; δ is the length of weir crest, m; P is the height of weir, m; S1 and S2 are up-and downstream slope coefficients, respectively; g is gravity acceleration, m/s 2 ; σ is surface tension of fluid, N/m; ρ is mass density of fluid, kg/m 3 ; μ is dynamic viscosity of fluid, Pa•s.Based on Buckingham's Π theorem, the three variables of H0, g and ρ are selected as the basic dimensions, and the relations between dimensionless numbers could be expressed as follows: The left term represents the fundamental weir formulae (Equation ( 3)) where m is the discharge coefficient.
The first four terms on the right of Equation ( 2) describe the effect of the geometrical properties on the discharge capacity, while the last two terms respectively describe the effects of surface tension and viscosity on the discharge capacity.Isaacs [25] and Ranga Raju et al. [26] concluded that the influence of viscosity and surface tension on the water flowing through weirs could be ignored if the flow depth over the crest was greater than 5.0 cm.Regarding the above conclusion the experiments select the flow depths over the crest between 8.0 cm and 24.0 cm.Regardless of the influence of surface tension and viscosity, Equation ( 2) is simplified as: ( ) Four dimensionless numbers are included in Equation ( 4): π1 = δ/H0, π2 = P/H0, π3 = S1 and π4 = S2, which show the effects of the geometrical parameters and hydraulic conditions on the discharge coefficient of short-crested weir.A good deal of research has studied the influence of weir crest length and weir height on discharge coefficient in detail, so in this study the two elements are settled as constants and the research only focuses on the effect of varying upstream and downstream slope coefficients on discharge coefficient of short-crested weir.

Numerical Modeling
A brief introduction to the numerical models is given below, including governing equations, grids layout, and the calculation of the free water surface.
The first four terms on the right of Equation (2) describe the effect of the geometrical properties on the discharge capacity, while the last two terms respectively describe the effects of surface tension and viscosity on the discharge capacity.Isaacs [25] and Ranga Raju et al. [26] concluded that the influence of viscosity and surface tension on the water flowing through weirs could be ignored if the flow depth over the crest was greater than 5.0 cm.Regarding the above conclusion the experiments select the flow depths over the crest between 8.0 cm and 24.0 cm.Regardless of the influence of surface tension and viscosity, Equation ( 2) is simplified as: Four dimensionless numbers are included in Equation (4): π 1 = δ/H 0 , π 2 = P/H 0 , π 3 = S 1 and π 4 = S 2 , which show the effects of the geometrical parameters and hydraulic conditions on the discharge coefficient of short-crested weir.A good deal of research has studied the influence of weir crest length and weir height on discharge coefficient in detail, so in this study the two elements are settled as constants and the research only focuses on the effect of varying upstream and downstream slope coefficients on discharge coefficient of short-crested weir.

Numerical Modeling
A brief introduction to the numerical models is given below, including governing equations, grids layout, and the calculation of the free water surface.

Governing Equations
In the Cartesian coordinate system, the incompressible continuity (5) and Reynolds-Averaged Navier-Stokes (RANS) (6) equations are used to solve the water motion for turbulent flow.
where u i and u j are average velocity components, m/s; x i and x j are Cartesian coordinate axes; f i is body force, m/s 2 ; p is pressure, Pa; −ρu i u j is the term of Reynolds shear stress with i, j = 1, 2, 3.In this study, the turbulence is predicted by the renormalization group (RNG) κ-ε turbulent model (turbulent kinetic energy κ and its dissipation ratio ε) [27].This model could dispose of the high shear and curvature area with higher accuracy and showed better performance in the simulated area with flow separation [28,29].The wall function method was used for the flow close to the wall and the flow with the lower Reynolds number.

Disposal of Free Water Surface
In this study, the Volume of Fluid (VOF) method proposed by Hirth et al. [30] is used to accurately predict the variation of the free water surface.The transport Equation ( 7) of fluid fraction is expressed as follows: ∂F ∂t where F is the fraction of fluid in a cell.This is a two-phase approach where fluid and air are simulated in a structured grid of the finite difference algorithm.The method is based on the concept that every cell has a fraction of fluid (F), which is 0 when the element is full of air and 1 when the element is totally full of fluid.If the value is between 0 and 1, the element contains the free fluid surface.

Boundary Conditions
In this study, the laboratory experiments were conducted in order to validate the Computational Fluid Dynamics (CFD) model results.Thus, the boundary conditions in numerical models should consist with the real boundaries in laboratory experiments.Stagnation pressure is used as inflow and outflow boundary, and the corresponding flow depths are specified respectively.Moreover, keep the outflow depth as low as possible to ensure the tail water will not influence the upstream inflow.In view of the symmetry of the weir model and the existence of free surface, the symmetry boundary conditions are set at the symmetrical face and on the top of the flume, respectively.All the solid walls-including the side wall and the flume bed-are considered as the no-slip wall boundary conditions.Moreover, the initial fluid level is specified to save the computational time.

Grids Layout
The gridding of computational field was shown in Figure 2. The whole computational domain was gridded by nonuniform mesh blocks consisting of rectangular elements.In Table 1 the gridding characteristics of numerical models for free surface profiles were provided.To check the accuracy of numerical model the statistical variables of the Root Mean Square Error (RMSE) and the Mean Absolute Percent Error (MAPE) were calculated using the following equations: where R m and R s are the laboratory results and simulated results, respectively.
Water 2018, 10, x FOR PEER REVIEW 5 of 18 where Rm and Rs are the laboratory results and simulated results, respectively.1, the differences between the results of meshing 3 and meshing 4 were negligible and gridding 3 was chosen.Therefore, the whole computational domain was gridded by 632, 18, 50 rectangular elements in the direction X, Y and Z, respectively.Moreover, the number of grids in x and z direction will be adjusted properly according to the geometry of the weir model and the upstream water depth.

Numerical Simulation
In this study, the flume with a length of 10.0 m, width of 0.3 m, and height of 0.5 m was defined.The weir model was located at a distance of 4.0 m from the inlet.In Table 2, the schemes of numerical simulation were listed.The schemes were expressed as Mij, in which i = 1-7 correspond to seven simulated upstream slope coefficients S1 of V, 0.3H:1V, 0.5H:1V, 0.8H:1V, 1.0H:1V, 1.5H:1V and 2.0H:1V respectively, while j = 1-9 correspond to nine downstream slope coefficient S2 of V, 0.4H:1V, 0.8H:1V, 1.0H:1V, 1.3H:1V, 1.5H:1V, 1.8H:1V, 2.0H:1V and 3.0H:1V-all with effective square edges at the joins.Based on the definition of the short-crested weir by Tong et al. [10] the short-crested weir was limited to According to the above criteria, the height of weir P, the length of weir δ and the width of weir B were constant during this study and the values are 24.0 cm, 16.0 cm and 30.0 cm, respectively.There were 63 weir models in total, and each weir model was calculated corresponding to five total energy heads H0 of 8.0 cm, 12.0 cm, 16.0 cm, 20.0 cm, and 24.0 cm over the crest.The discharge coefficients of the numerical results could be obtained according to Equation (3).As shown in Table 1, the differences between the results of meshing 3 and meshing 4 were negligible and gridding 3 was chosen.Therefore, the whole computational domain was gridded by 632, 18, 50 rectangular elements in the direction X, Y and Z, respectively.Moreover, the number of grids in x and z direction will be adjusted properly according to the geometry of the weir model and the upstream water depth.

Numerical Simulation
In this study, the flume with a length of 10.0 m, width of 0.3 m, and height of 0.5 m was defined.The weir model was located at a distance of 4.0 m from the inlet.In Table 2, the schemes of numerical simulation were listed.The schemes were expressed as M ij , in which i = 1-7 correspond to seven simulated upstream slope coefficients S 1 of V, 0.3H:1V, 0.5H:1V, 0.8H:1V, 1.0H:1V, 1.5H:1V and 2.0H:1V respectively, while j = 1-9 correspond to nine downstream slope coefficient S 2 of V, 0.4H:1V, 0.8H:1V, 1.0H:1V, 1.3H:1V, 1.5H:1V, 1.8H:1V, 2.0H:1V and 3.0H:1V-all with effective square edges at the joins.Based on the definition of the short-crested weir by Tong et al. [10] the short-crested weir was limited to 0.67 < δ/H ≤ 1.5 ∼ 2.0,0.5 ≤ P/H < 3.0 (10) According to the above criteria, the height of weir P, the length of weir δ and the width of weir B were constant during this study and the values are 24.0 cm, 16.0 cm and 30.0 cm, respectively.There were 63 weir models in total, and each weir model was calculated corresponding to five total energy heads H 0 of 8.0 cm, 12.0 cm, 16.0 cm, 20.0 cm, and 24.0 cm over the crest.The discharge coefficients of the numerical results could be obtained according to Equation (3).

Validation of Numerical Models
To validate the numerical model results, the laboratory experiments were conducted.The physical model consisted of a self-circular, horizontal, rectangular flume with a short-crested weir.The side walls and weirs were made of Plexiglas.The upstream water level was measured with an accuracy of ±0.1 mm using a needle water level gauge and the volume flow rate with ±0.1 mm using 90 • V-notch weir.The tail water depth was selected in a way that would not affect the incoming flow.Table 3 listed the selected geometrical and hydraulic characteristic parameters used to validate the accuracy of the numerical models.The comparison of discharge coefficients between simulated and laboratory results at free overflow condition are shown in Figure 3.The discharge coefficients gradually increase as the total energy head H 0 increases, and the simulated results are slightly larger than those of laboratory results.Figure 4 shows the comparisons of free surface profiles between the simulated and laboratory results, in which the origin locates at the entrance of the weir crest and h is the fluid level along the streamwise direction and ξ = H 0 /(P + δ) is the relative total energy head.The deviation from free surface profiles might result from flow separation and the slightly less accurate simulation for the shapely varied free-surface profile.As shown in Figure 5, the changing trend of free surface elevation is nearly identical for simulated and laboratory results and the bottom roll backflow zone simultaneously appear in back of downstream weir face.Table 4 lists the RMSE and MAPE for the discharge coefficients and free surface profiles.In consideration of the comparative results it is obvious that the numerical model is able to accurately predict the water surface over the short-crested weir and the directly linked discharge coefficient.

Calculation Formula of Discharge Coefficient
Sixty-three weir models and 315 numerical simulations are conducted.From the dimensional analysis of Equation ( 4) a monotonically increasing linear relation between m and lnξ for different weir models is expressed as follows: where m′ is the fitted discharge coefficient and α and β are dimensionless parameters determined by S1 and S2.In Combining Equations ( 12)-( 15), the calculation formula of discharge coefficient of rectangular short-crested weir with varying upstream and downstream slope coefficients (URD) is derived as Equation (11), where ) The Equation ( 11) is valid for 0.0 ≤ S1 ≤ 2.0, 0.0 ≤ S2 ≤ 3.0, 0.67 ≤ δ/H ≤ 2.0, 1.0 ≤ P/H ≤ 3.0 and −1.61 ≤ lnξ ≤ −0.51.The range out of simulated results should be verified further.The deviation from m is plotted in Figure 7.The RMSE, MAPE and R 2 for discharge coefficient are computed with the values 0.36%, 0.02% and 0.977, respectively.Recently, it has not have a discharge coefficient formula

Calculation Formula of Discharge Coefficient
Sixty-three weir models and 315 numerical simulations are conducted.From the dimensional analysis of Equation ( 4) a monotonically increasing linear relation between m and lnξ for different weir models is expressed as follows: where m is the fitted discharge coefficient and α and β are dimensionless parameters determined by S 1 and S 2 .In Figure 6a-b the effects of S 2 on α and β are shown for all of the upstream slope coefficients S 1 , and the corresponding values are listed in Table 5.Both of α and β firstly increase until they reach the maximum of α max and β max at S 2 = 0.8, and then decrease until they reach the minimum of α min and β min at S 2 = 3.0.Moreover, a good linear relation between α max and S 1 and a good parabolic relation between β max and S 1 are shown in Figure 6c,d, and the relations could be expressed as follows: with R 2 = 0.989 and 0.979 respectively.The data are normalized by introducing the variables of α = α/α max and β = β/β max as shown in Figure 6e,f.Both of α and β have good piecewise linear relations with S 2 , which could end up with the expressions: Combining Equations ( 12)-( 15), the calculation formula of discharge coefficient of rectangular short-crested weir with varying upstream and downstream slope coefficients (URD) is derived as Equation (11), where α = (−1.3S The Equation ( 11) is valid for 0.0 ≤ S 1 ≤ 2.0, 0.0 ≤ S 2 ≤ 3.0, 0.67 ≤ δ/H ≤ 2.0, 1.0 ≤ P/H ≤ 3.0 and −1.61 ≤ lnξ ≤ −0.51.The range out of simulated results should be verified further.The deviation from m is plotted in Figure 7.The RMSE, MAPE and R 2 for discharge coefficient are computed with the values 0.36%, 0.02% and 0.977, respectively.Recently, it has not have a discharge coefficient formula including S 1 and S 2 .The discharge coefficient m 1 calculated by Equation ( 11) at S 1 = 0 and S 2 = 0 is compared with the Equation (18) proposed by Govinda Rao and Muralidhar [1], and comparative results are shown in the Figure 8 with the maximum relative error 1.92%.
The application range of Equation ( 18) is 0.4 ≤ H/δ ≤ 1.5-1.9.Analyzing Equation ( 11) concludes that the geometrical configuration corresponding to the maximum discharge coefficient relates closely to the value range of lnξ.In the range −1.06 < lnξ ≤ −0.51 (13.8 cm < H 0 ≤ 24 cm), m reaches m max = (7.43lnξ+ 49.53) × 10 −2 at S 1 = 0.459 − 0.75lnξ and S 2 = 0.8; in the range −1.61 ≤ lnξ ≤ −1.06 (8.0 cm ≤ H 0 ≤ 13.8 cm), m reaches m max = (5.46lnξ+ 47.49) × 10 −2 at S 1 = 2.0 and S 2 = 0.8.And the maximum discharge coefficient is always obtained at S 2 = 0.8.The existence of a maximum discharge coefficient at S 2 = 0.8 relates to the flow regime in back of downstream weir face.Comparisons of the velocity fields and free-surface elevations of weirs at S 1 = 1.0 with H/δ = 1.41 for S 2 = 0.0, 0.8 and 2.0 are shown in Figures 9 and 10, respectively.In Figure 9a, a bottom roll backflow zone forms between the downstream weir face and the nappe at S 2 = 0.0, and air is entrained into it simultaneously.The backflow will lift the nappe, resulting in the rise of downstream free surface profile, and reduce the discharge capacity.Moreover, the outline of backflow zone roughly forms an isosceles triangle on the side view, which is identical to the test observation.In Figure 9b, the backflow zone disappears and the downstream free surface profile falls after a rise at S 2 = 0.8 compared with the weir at S 2 = 0.0.In this situation the nappe attached on the weir face perfectly, which improves the discharge capacity.With increasing S 2 , the reattachment point produced by the combined action of inertial force and gravity is gradually covered by the downstream weir side, and the downstream free surface profile is lifted higher by the weir face than the weir at S 2 = 0.0 as shown in Figures 9c and 10.
Water 2018, 10, x FOR PEER REVIEW 9 of 18 including S1 and S2.The discharge coefficient m1 calculated by Equation ( 11) at S1 = 0 and S2 = 0 is compared with the Equation ( 18) proposed by Govinda Rao and Muralidhar [1], and comparative results are shown in the Figure 8 with the maximum relative error 1.92%.The application range of Equation ( 18) is 0.4 ≤ H/δ ≤ 1.5-1.9.Analyzing Equation ( 11) concludes that the geometrical configuration corresponding to the maximum discharge coefficient relates closely to the value range of lnξ.In the range −1.06 < lnξ ≤ −0.51 (13.8 cm < H0 ≤ 24 cm), m′ reaches m′max = ( 7  respectively.In Figure 9a, a bottom roll backflow zone forms between the downstream weir face and the nappe at S2 = 0.0, and air is entrained into it simultaneously.The backflow will lift the nappe, resulting in the rise of downstream free surface profile, and reduce the discharge capacity.Moreover, the outline of backflow zone roughly forms an isosceles triangle on the side view, which is identical to the test observation.In Figure 9b, the backflow zone disappears and the downstream free surface profile falls after a rise at S2 = 0.8 compared with the weir at S2 = 0.0.In this situation the nappe attached on the weir face perfectly, which improves the discharge capacity.With increasing S2, the reattachment point produced by the combined action of inertial force and gravity is gradually covered by the downstream weir side, and the downstream free surface profile is lifted higher by the weir face than the weir at S2 = 0.0 as shown in Figures 9c and 10

Effects of Slope Coefficients on Discharge Coefficient
When S2 = 0.0, Equation ( 11) would be simplified as: ( ) where mU′ is the fitted discharge coefficient for S2 = 0.0.Analyzing the fitted equation indicates that for a given lnξ, mU′ firstly increases until it reaches the maximum at S1 = 0.5822 − 0.7055lnξ, and then decreases with S1.
In Figure 11, the difference in discharge coefficients for different values of S1 shows a decreasing trend with increasing lnξ.In view of energy conservation law, the difference is produced by the frictional head loss of hf and local head loss of hj in the progress of free overflow.When lnξ holds constant, the gradually slowing upstream slope increases the flow run and reduces the streamline curvature of inflow, accordingly hf is increasing and hj is decreasing.In this progress, the increment of the former is less than the decrement of the latter, so mU′ continues increasing till S1 = 0.5822 − 0.7055lnξ.When the upstream slope continued slowing, the increment of hf exceeds the decrement of hj, so the discharge coefficient is gradually decreasing.With increasing total energy head of H0 over the crest the difference in the total head loss of hw where hw = hf + hj is gradually decreasing, so the difference in discharge coefficients is gradually decreasing for different values of S1.

Effects of Slope Coefficients on Discharge Coefficient
When S 2 = 0.0, Equation ( 11) would be simplified as: where m U is the fitted discharge coefficient for S 2 = 0.0.Analyzing the fitted equation indicates that for a given lnξ, m U firstly increases until it reaches the maximum at S 1 = 0.5822 − 0.7055lnξ, and then decreases with S 1 .
In Figure 11, the difference in discharge coefficients for different values of S 1 shows a decreasing trend with increasing lnξ.In view of energy conservation law, the difference is produced by the frictional head loss of h f and local head loss of h j in the progress of free overflow.When lnξ holds constant, the gradually slowing upstream slope increases the flow run and reduces the streamline curvature of inflow, accordingly h f is increasing and h j is decreasing.In this progress, the increment of the former is less than the decrement of the latter, so m U continues increasing till S 1 = 0.5822 − 0.7055lnξ.When the upstream slope continued slowing, the increment of h f exceeds the decrement of h j , so the discharge coefficient is gradually decreasing.With increasing total energy head of H 0 over the crest the difference in the total head loss of h w where h w = h f + h j is gradually decreasing, so the difference in discharge coefficients is gradually decreasing for different values of S 1 .
When S 1 = 0.0, Equation ( 11) would be simplified as: where m D is the fitted discharge coefficient for S 1 = 0.0.It reveals that the maximum of discharge coefficient m D max depends on both of S 2 and ξ.In the range −1.29 < lnξ ≤ −0.51, m D reaches the maximum of m D max = (7.85lnξ+ 49.17In Figure 12, the difference in discharge coefficients for different values of S2 shows an increasing trend with increasing lnξ.Similarly, the difference is caused by hf and hj, in which the local head loss is mainly produced by inflow dropping into the downstream with forming the bottom roll backflow zone between the downstream weir face and the nappe.Moreover, the size of the backflow zone is proportionate to the local head loss.In the range −1.29 < lnξ ≤ −0.51, increasing S2 will gradually decrease the size of the backflow zone until it disappears at S2 =0.8, and the local head loss plays the decisive role in this progress.Continuing increasing S2, hf tends to be dominant with increasing flow run.In the range −1.61 ≤ lnξ ≤ −1.29, although the backflow zone cannot be formed due to the lower crest depth, the subpressure zone is formed around the downstream weir face.Smaller the value of S2 is, larger the maximum of subpressure zone is.Although it increases the discharge coefficient, the instability of subpressure zone can lead to the instability of weir flow.Therefore, the design of short-crested weir should take the downstream slope into consideration.

Effects of Slope Coefficients on Free Surface Profiles
In this numerical study, the variations of free surface profiles and Froude numbers over the crest with various slope coefficients are investigated in overflow conditions.In Figure 13, the normalized free surface elevation of (h − P)/H is plotted against x/δ for different values of H/δ.As shown, the free surface profiles are self-similar.The range of (h − P)/H shows a decreasing trend with growing x/δ in all cases.The free surface profile can be divided into two parts due to the gradient: a slow descending segment and a sharp descending segment.For the former the curvature is nearly 0, but for the latter it descends quickly with hydraulic drop happening.The flow separation phenomenon happening around x/δ = 0.0 (at the entrance of weir crest) produces large In Figure 12, the difference in discharge coefficients for different values of S 2 shows an increasing trend with increasing lnξ.Similarly, the difference is caused by h f and h j , in which the local head loss is mainly produced by inflow dropping into the downstream with forming the bottom roll backflow zone between the downstream weir face and the nappe.Moreover, the size of the backflow zone is proportionate to the local head loss.In the range −1.29 < lnξ ≤ −0.51, increasing S 2 will gradually decrease the size of the backflow zone until it disappears at S 2 =0.8, and the local head loss plays the decisive role in this progress.Continuing increasing S 2 , h f tends to be dominant with increasing flow run.In the range −1.61 ≤ lnξ ≤ −1.29, although the backflow zone cannot be formed due to the lower crest depth, the subpressure zone is formed around the downstream weir face.Smaller the value of S 2 is, larger the maximum of subpressure zone is.Although it increases the discharge coefficient, the instability of subpressure zone can lead to the instability of weir flow.Therefore, the design of short-crested weir should take the downstream slope into consideration.In Figure 12, the difference in discharge coefficients for different values of S2 shows an increasing trend with increasing lnξ.Similarly, the difference is caused by hf and hj, in which the local head loss is mainly produced by inflow dropping into the downstream with forming the bottom roll backflow zone between the downstream weir face and the nappe.Moreover, the size of the backflow zone is proportionate to the local head loss.In the range −1.29 < lnξ ≤ −0.51, increasing S2 will gradually decrease the size of the backflow zone until it disappears at S2 =0.8, and the local head loss plays the decisive role in this progress.Continuing increasing S2, hf tends to be dominant with increasing flow run.In the range −1.61 ≤ lnξ ≤ −1.29, although the backflow zone cannot be formed due to the lower crest depth, the subpressure zone is formed around the downstream weir face.Smaller the value of S2 is, larger the maximum of subpressure zone is.Although it increases the discharge coefficient, the instability of subpressure zone can lead to the instability of weir flow.Therefore, the design of short-crested weir should take the downstream slope into consideration.

Effects of Slope Coefficients on Free Surface Profiles
In this numerical study, the variations of free surface profiles and Froude numbers over the crest with various slope coefficients are investigated in overflow conditions.In Figure 13, the normalized free surface elevation of (h − P)/H is plotted against x/δ for different values of H/δ.As shown, the free surface profiles are self-similar.The range of (h − P)/H shows a decreasing trend with growing x/δ in all cases.The free surface profile can be divided into two parts due to the gradient: a slow descending segment and a sharp descending segment.For the former the curvature is nearly 0, but for the latter it descends quickly with hydraulic drop happening.The flow separation phenomenon happening around x/δ = 0.0 (at the entrance of weir crest) produces large

Effects of Slope Coefficients on Free Surface Profiles
In this numerical study, the variations of free surface profiles and Froude numbers over the crest with various slope coefficients are investigated in overflow conditions.In Figure 13, the normalized free surface elevation of (h − P)/H is plotted against x/δ for different values of H/δ.As shown, the free surface profiles are self-similar.The range of (h − P)/H shows a decreasing trend with growing x/δ in all cases.The free surface profile can be divided into two parts due to the gradient: a slow descending segment and a sharp descending segment.For the former the curvature is nearly 0, but for the latter it descends quickly with hydraulic drop happening.The flow separation phenomenon happening around x/δ = 0.0 (at the entrance of weir crest) produces large energy losses and results in free surface profile descending sharply.Moreover, the range of the curvature of free surface profile shows a decreasing trend with growing H/δ in all cases.curved surface profile.When Fr > 1, the value of Fr for the same x/δ increases as H/δ decreases, indicating that for a smaller flow rate, the ratio of inertial force to gravity is so large it can keep the higher profile gradient.Moreover, for a given weir model, the position where supercritical flow takes place is relatively rearward as H/δ decreases.
Comparisons of free surface profiles and Fr over the weirs at S2 = 0.8 for S1 = 0.0, 1.0 and 2.0 are shown in Figure 13a-c and Figure 14a-c.The plots reveal that increasing S1 leads to the increase of surface curvature before the crest entrance and the decrease over the crest.When H/δ holds constant, increasing S1 will accelerate supercritical flow taking place.Analyzing the effect of different S2 on free surface profiles over crest at S1 = 1.0 in Figure 13b,d,e and the corresponding Froude numbers in Figure 14, it can be concluded that increasing S2 leads to the decrease of free surface curvature over the crest and it has no obvious influence on the curvature of free surface profile before the crest entrance.In Figure 14, Fr is plotted against x/δ for different values of H/δ.The range of Fr shows an increasing trend with growing x/δ in all cases.When Fr < 1, the value of F r for the same x/δ decreases Water 2018, 10, 204 15 of 18 as H/δ decreases, which indicates that for a smaller flow rate, the inertial force on the flow is smaller.Under the action of gravity, the free surface profile is easier to bend, resulting in a more curved surface profile.When Fr > 1, the value of Fr for the same x/δ increases as H/δ decreases, indicating that for a smaller flow rate, the ratio of inertial force to gravity is so large it can keep the higher profile gradient.Moreover, for a given weir model, the position where supercritical flow takes place is relatively rearward as H/δ decreases. (e)

Conclusions
The effects of varying upstream and downstream slope coefficients on discharge coefficients of rectangular short-crested weirs were studied systemically by the combination of laboratory experiments and numerical simulation.A calculation formula was proposed by using the dimensional analysis between hydraulic parameters and geometrical parameters: height of weir, length of weir, upstream overflow total energy head, and upstream and downstream slope Comparisons of free surface profiles and Fr over the weirs at S 2 = 0.8 for S 1 = 0.0, 1.0 and 2.0 are shown in Figures 13a-c and 14a-c.The plots reveal that increasing S 1 leads to the increase of surface curvature before the crest entrance and the decrease over the crest.When H/δ holds constant, increasing S 1 will accelerate supercritical flow taking place.Analyzing the effect of different S 2 on free surface profiles over crest at S 1 = 1.0 in Figure 13b,d,e and the corresponding Froude numbers in Figure 14, it can be concluded that increasing S 2 leads to the decrease of free surface curvature over the crest and it has no obvious influence on the curvature of free surface profile before the crest entrance.

Conclusions
The effects of varying upstream and downstream slope coefficients on discharge coefficients of rectangular short-crested weirs were studied systemically by the combination of laboratory experiments and numerical simulation.A calculation formula was proposed by using the dimensional analysis between hydraulic parameters and geometrical parameters: height of weir, length of weir, upstream overflow total energy head, and upstream and downstream slope coefficients.Moreover, the free surface profiles over the weirs were numerically investigated.
The upstream overflow total energy head played a decisive role on the discharge coefficient of rectangular short-crested weirs.For a certain configuration of a short-crested weir, the discharge coefficient is slightly increased as the total energy head H 0 increases.As H 0 increases, the difference in discharge coefficients for different values of S 1 shows a decreasing trend with increasing lnξ, while the trend is opposite for different values of S 1 .For a given total energy head over the crest, the rectangular short-crested weir with varying upstream slope has a convex parabolic relation between S 1 and the discharge coefficient, while the weir with varying downstream slope has a relation of piecewise linearity between the S 2 and the discharge coefficient.For a higher H 0 , the maximum discharge coefficient is always obtained at 0.8H:1V of downstream slope.
Increasing the upstream slope coefficients leads to the free surface curvature increasing before the crest entrance and decreasing over the crest, while increasing the downstream slope coefficients leads to the decrease of free surface curvature on the crest.In addition, the structural design for short-crested weirs should take the downstream slope into consideration.
The calculation formula of the discharge coefficient was derived by the nonlinear regression.It is a function of the upstream slope coefficient (S 1 ), the downstream slope coefficient (S 2 ), and the ratio of total energy head to the sum of crest length and crest height (H 0 /(P + δ)).The RMSE, MAPE and R 2 for the proposed equation are calculated, and the values are 0.36%, 0.02% and 0.977, respectively.

Figure 1 .
Figure 1.Definition sketch of free overflow over rectangular short-crested weir with sloped weir faces.

Figure 1 .
Figure 1.Definition sketch of free overflow over rectangular short-crested weir with sloped weir faces. m

Figure 2 .
Figure 2. Sketch of numerical model: (a) solving domain; (b) grids layout around the weir.

Figure 2 .
Figure 2. Sketch of numerical model: (a) solving domain; (b) grids layout around the weir.
Figure 6a-b the effects of S2 on α and β are shown for all of the upstream slope coefficients S1, and the corresponding values are listed in Table5.Both of α and β firstly increase until they reach the maximum of αmax and βmax at S2 = 0.8, and then decrease until they reach the minimum of αmin and βmin at S2 = 3.0.Moreover, a good linear relation between αmax and S1 and a good parabolic relation between βmax and S1 are shown in Figure6c,d, and the relations could be expressed as follows: R 2 = 0.989 and 0.979 respectively.The data are normalized by introducing the variables of α′ = α/αmax and β′ = β/βmax as shown in Figure6e,f.Both of α′ and β′ have good piecewise linear relations with S2, which could end up with the expressions:

Figure 7 .
Figure 7.Comparison between m′ and m.

Figure 7 .
Figure 7.Comparison between m′ and m.

Figure 7 .
Figure 7.Comparison between m′ and m.

Table 1 .
Gridding characteristics, RMSE and MAPE for the simulated profiles of free surface.

Table 1 .
Gridding characteristics, RMSE and MAPE for the simulated profiles of free surface.

Table 2 .
List of design schemes of numerical simulation.

Table 3 .
Range of variables for laboratory measurements.

Table 4 .
RMSE and MAPE for discharge coefficients and flow surface profiles.