Discharge formula for flows over open-check dams

Discharge formula for flows over open-check dams

Journal Pre-proof Discharge formula for flows over open-check dams Hsun-Chuan Chan, Hsin-Kai Yang, Po-Wei Lin, Jung-Tai Lee PII: S0955-5986(19)30230-...

3MB Sizes 8 Downloads 59 Views

Journal Pre-proof Discharge formula for flows over open-check dams Hsun-Chuan Chan, Hsin-Kai Yang, Po-Wei Lin, Jung-Tai Lee PII:

S0955-5986(19)30230-4

DOI:

https://doi.org/10.1016/j.flowmeasinst.2019.101690

Reference:

JFMI 101690

To appear in:

Flow Measurement and Instrumentation

Received Date: 21 June 2019 Revised Date:

30 December 2019

Accepted Date: 30 December 2019

Please cite this article as: H.-C. Chan, H.-K. Yang, P.-W. Lin, J.-T. Lee, Discharge formula for flows over open-check dams, Flow Measurement and Instrumentation (2020), doi: https://doi.org/10.1016/ j.flowmeasinst.2019.101690. This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2019 Published by Elsevier Ltd.

1

Discharge formula for Flows over Open-Check Dams

2

Hsun-Chuan Chan 1, Hsin-Kai Yang 2, Po-Wei Lin 3, and Jung-Tai Lee 4,*

3 4 5 6 7 8 9 10 11 12

1

Professor, Department of Soil and Water Conservation, National Chung Hsing University, 145 Xingda Rd., Taichung 402, Taiwan (ROC). (Email: [email protected])

2

Ph.D. Student, Department of Soil and Water Conservation, National Chung Hsing University, 145 Xingda Rd., Taichung 402, Taiwan (ROC). (Email: [email protected])

3

Ph.D. Student, Department of Soil and Water Conservation, National Chung Hsing University, 145 Xingda Rd., Taichung 402, Taiwan (ROC). (Email: [email protected])

4

Assistant Professor, Department of Forestry and Natural Resources, National Chiayi University, 300 Syuefu Rd., Chiayi City 60004, Taiwan (ROC). (Email: [email protected])

* Correspondence: [email protected]; Tel: 886-5-2717482; Fax: 886-5-2717467.

Abstract

13

The present study proposed a discharge prediction model for open-check dams.

14

The open-check dam is composed of a trapezoidal spillway and a rectangular opening

15

located at the bottom of the spillway. The flow conditions, determined using the ratio

16

of water head and dam crest length, for flow over open-check dams were divided into

17

broad-crested and sharp-crested weir flow conditions. For the broad-crested weir flow

18

condition, the discharge equation for rectangular broad-crested weir was proposed to

19

estimate the discharge. For the sharp-crested weir flow condition, a linear

20

combination of discharge equations for rectangular and trapezoidal sharp-crested

21

weirs was proposed to describe the discharge. A series of laboratory experiments were

22

performed in order to calibrate the proposed discharge model. The results showed the

23

discharge coefficient and cross-sectional area ratio ( Ar ) were highly correlated under

24

two flow conditions. The averaged error between measured and predicted discharge

25

was less than 2%. The previous models displayed averaged errors between 4% and

26

25% under the present experimental conditions, whereas that of the present model

27

was considerably low. The present model demonstrated favorable accurate and

28

convenient estimation of discharge for flows over open-check dams.

29

Keywords: open-check dam; discharge; discharge coefficient; open channel flow 1

1

1. Introduction

2

Check dams are commonly constructed in upstream catchment areas to minimize

3

sediment disasters. However, because check dam design involves accumulating a

4

substantial amount of sediment, dam heights are often increased to expand sediment

5

storage space. Thus, check dam construction projects often have adverse effects on

6

river habitats and the continuity of sediment transport. In recent years, the increased

7

awareness of ecological conservation involved adding a rectangular opening

8

underneath the conventional spillway in check dams. Figs. 1(a) and (b) show a

9

conventional check dam with a trapezoidal spillway and that of an open-check dam

10

with a trapezoidal spillway and a rectangular opening respectively. The rectangular

11

opening in a check dam reduces the gap in elevation between the upstream and

12

downstream of the riverbed which maintains the continuity of sediment and aquatic

13

organism, thereby ameliorating the negative effects of the check dam on the natural

14

environments. In the design process of a conventional check dam (Piton & Recking,

15

2016), the discharge formula for sharp-crested weirs proposed by Zollinger (1983) is

16

commonly employed to decide the size of the spillway. However, open-check dams

17

differ considerably from conventional check dams in cross-sectional shape (the

18

addition of a rectangular opening). The use of the aforementioned discharge formula

19

to predict discharge may yield erroneous results. To effectively predict the discharge

20

of open-check dams, the present study superposed the discharge of the spillway and

21

opening of open-check dams. Considering the variations of water depth resulted in the

22

change in flow conditions, relevant data were incorporated into the discharge

23

coefficients to build a discharge prediction model of open-check dams under different

24

flow conditions. Five open-check dam models were tested in a laboratory flume for

25

large range of discharges. The predicted results are in agreement with the

26

experimental measurements. The model developed in the present study is expected as 2

1

the reference of the future researches in disaster prevention, sediment transport, and

2

habitat stability of open-check dams.

(a) Conventional check dam

(b) Open-check dam

Fig. 1. Different types of check dams. 3

2. Background

4

Check dams are transverse structures of rivers, which block water flow and cause

5

water levels to rise. Roushangar et al. (2014) asserted that water passing through

6

check dam spillways exhibits flow characteristics, including nappe flow, curved

7

streamlines, and air vent, similar to those of weir flows. Fig. 2 shows the flow patterns

8

of water flowing through a check dam. Regarding the flow conditions and the

9

geometry of the spillways, numerous studies have predicted the discharge of check

10

dams by using weir flow equations (Zollinger, 1983; Fang & Fang, 2012; Piton &

11

Recking, 2016).

Fig. 2. Flow patterns of water flowing through a check dam. 12 13

2.1. Flows over Weir According to Govinda Rao and Muralidhar (1963), weir flow conditions can be 3

1

divided into broad-crested and sharp-crested on the basis of weir water head and weir

2

crest length ratios ( h L ; where h and L are the water depth above crest weir and

3

the weir crest length, respectively).

4

2.1.1. Broad-crested Weir Flow Condition

5

When h L ≤ 1.5 , water flow makes contact with the weir crest, and the

6

friction created will cause energy loss in the water flow. This flow condition is

7

referred to as a broad-crested weir (Fig. 3(a)). The water touching the weir crest

8

upstream undergoes water flow separation, which, because of the effect of gravity,

9

causes the water to again come into contact with the weir crest. The critical water

10

depth occurs at the highest point of the water flow separation. Thus, discharge

11

predictions are often made using critical flow conditions. However, because

12

discharge prediction formulas ignore energy losses in the water flow due to

13

friction with the weir crest, the discharge predicted is generally underestimated

14

(Govinda Rao & Muralidhar, 1963; Azimi & Rajaratnam, 2009; Azimi et al.,

15

2014; Bijankhan et al., 2014).

16

2.1.2. Sharp-crested Weir Flow Condition

17

When h L > 1.5 , flows are considered to those of sharp-crested weir flows

18

(Fig. 3(b)). Because the water depth of a sharp-crested weir is greater than that of

19

a broad-crested weir, water flow in a sharp-crested weir possesses greater

20

potential energy. In a sharp-crested weir, when the water flow makes contact with

21

the weir crest upstream, the energy lost in the water flow is less than that lost in a

22

broad-crested weir. Because most of the potential energy converts into kinetic

23

energy, the water flows directly over the weir crest and lose contact with the crest.

24

In this situation, the streamlines are highly curved. Because discharge prediction

25

formulas for sharp-crested weirs ignore energy losses in the water flow due to 4

1

highly curved nappe flow, the discharge coefficients predicted are generally

2

overestimated (French, 1986; Swamee, 1988; Bos, 1989; Rehbock, 1929; Aydin et

3

al., 2011).

(a) Broad-crested weir flow condition

(b) Sharp-crested weir flow condition

Fig. 3. Weir flow condition classification. 4

2.2. Discharge Estimation of Weir Flows

5

Weirs can have a variety of cross sections, depending on the requirements. The

6

simplest type is the weir with a rectangular cross-section. The shapes of the single

7

cross-section opening are mostly rectangular and trapezoidal. By contrast, when the

8

opening shapes are a combination of different shapes, the weirs are referred to as the

9

compound cross-section weirs. Many stage-discharge relationships were used to

10 11

predict discharge of weir flows. 2.2.1. Weirs with a rectangular cross-section

12

For rectangular sharp-crested weirs, Rehbock (1929) conducted experimental

13

research and found that the discharge coefficients are determined by their water

14

head–height ratios

15

passing through a weir, French (1986) proposed using constriction factors to

16

identify the effect of cross section constriction on discharge coefficients, where

17

the discharge coefficients range from 0.59 to 0.64. In the rectangular

18

broad-crested weirs, Govinda Rao and Muralidhar (1963) considered the

(h P ) .

Because flow exhibits a constriction effect while

5

1

discharge coefficient to be a function of h L in order to take into account the

2

friction loss caused by weir crest. In addition, they revealed that the discharge

3

coefficient of a rectangular broad-crested weir ranged from 1.02 to 1.28. A

4

comprehensive review of the experimental investigations on finite crest length

5

weirs can be found in Azimi and Rajaratnam (2009). They also developed

6

empirical correlations for the discharge coefficient.

7

2.2.2 Weirs with Compound Cross-Section

8

Discharge over weirs with the compound cross-section have investigated by

9

Jan et al. (2005), Salmasi (2011), and Gogus (2016). Jan et al. (2005) discussed

10

four different compound weirs: weirs with rectangular compound cross-section

11

opening at both the top and bottom, weirs with a trapezoidal compound

12

cross-section opening at the top and rectangular compound cross-section opening

13

at the bottom, weirs with rectangular compound cross-section opening at the top

14

and triangular cross-section opening at the bottom, and weirs with trapezoidal

15

compound cross-section opening at the top and triangular cross-section opening at

16

the bottom. They calculated the discharges of the various cross sections and

17

linearly superposed them. However, Jan et al. (2005) only examined the high

18

discharge (i.e., when the flow flowed through both openings), and ignored

19

instances when the discharge was low (i.e., when the flow flowed through only the

20

bottom cross-section opening), consequently limiting the applicability of the study

21

findings. Salmasi (2011) and Gogus (2016) stated that the discharge coefficients

22

of compound weirs should change when water depths change. They investigated

23

the flow discharge for water flowing through only the bottom rectangular

24

compound cross-section opening and water flowing through both the top and

25

bottom openings. Salmasi (2011) showed that when 0.1 < h L < 0.6 , the 6

1

discharge coefficients of compound weirs change with h L ; in addition, the

2

trends for water flowing through only the bottom rectangular cross-section

3

opening differ from those for water flowing through both the top and bottom

4

rectangular cross-section openings, with discharge coefficients for these two

5

opening types ranging from 0.7 to 1.0. Gogus (2016) utilized flow coefficient Cv

6

to correct the negative effects created when the influence of velocity heads was

7

ignored during the theoretical derivation process, thereby elevating discharge

8

prediction accuracy. Gogus (2016) conducted experiments and revealed that when

9

h L < 1.7 , the discharge coefficients of compound weirs range from 0.74 to 1.10.

10

3. Theoretical Analysis

11

When water flows through check dams, the water depth ( hd ) and crest length

12

( Ld ) ratio ( hd Ld ) changes constantly due to changes in water level. Fig. 4 shows the

13

water flowing through a check dam for different discharges. Roushangar et al. (2014)

14

suggested that flow conditions through check dams are similar to that through weirs,

15

and that weir flow condition changes with the ratio of weir water depth and weir crest

16

length. Therefore, to ensure the prediction model effectively estimates the discharge

17

in different flow conditions, the present study used the weir flow classification

18

proposed by Govinda Rao and Muralidhar (1963) and divided the flows of check

19

dams into broad-crested weir and sharp-crested weir flow conditions according to the

20

ratios of dam water depth and dam crest length. Furthermore, the discharge prediction

21

model for open-check dams was developed under different flow conditions.

22 23 7

1 2

(a) Low discharge

(b) High discharge

Fig. 4. Water flows over open-check dam with different discharges. 3

3.1. Broad-Crested Weir Flow

4

When flowing through open-check dam with hd Ld ≤ 1.5 , the amount of

5

water passing through the dams is considered to have a low discharge. At this

6

time, most of the water flows through the rectangular opening below the

7

trapezoidal spillway. Fig. 5(a) shows the dimensions of an open-check dam, where

8

b0 and hd are the width of the rectangular opening and the water depth at the

9

head measurement section, respectively. Here, because the water in the upstream

10

of the dam contains scant potential energy, it converts into little kinetic energy

11

when the water flows through the weir crest, resulting in relatively short water

12

flow separation distances and causing the water to come in contact with the weir

13

crest. Regarding hydraulics, the water flow is classified as a broad-crested weir

14

flow. In this flow condition, the discharge of open-check dam was predicted using

15

the broad-crested weir water flow method.

8

(a) Front view

(b) Side view

Fig. 5. Definition sketch for broad-crested weir flow over open-check dam. 1

Fig. 5(b) presents the height, width energy head of a dam. Cross Section 1 shows

2

a stable water flow in the upstream of the dam, whereas Cross Section 2 reveals the

3

separation of water when the water comes in contact with the weir crest. Gogus

4

(2016) hypothesized that there are no energy losses from Cross Section 1 to Cross

5

Section 2 and that critical flow conditions occur at Cross Section 2. The energy

6

relationships between the two cross sections can be obtained using the following

7

energy equation:

8

H 1 = h2 +

9

where H 1 is the total energy head of Cross Section 1, h2 is the water depth of

10

Cross Section 2, V2 is the flow velocity of Cross Section 2, g is the gravitational

11

acceleration, yc is the critical water depth , and Vc is the critical velocity. By

12

incorporating the relationship between critical velocity and total energy head obtained

13

from Equation (3-1) into the flow rate formula (Q = A × V ) , the equation of discharge

14

can be obtained as:

15

Qb = Ac 2g(H1 − yc )

V22 V2 = yc + c 2g 2g

(1)

(2)

9

1

where Ac is the cross-sectional area where critical flow occurred and Qb is the

2

discharge of the check dam under broad-crested weir water flow condition. Because

3

hd Ld ≤ 1.5 , the cross-sectional shape of the water passing through the open-check

4

dam is rectangular. Assuming that the velocity head of Cross Section 1 can be

5

ignored, total energy head is therefore equal to the water depth at the head

6

measurement section (H1 = hd ) . Because the flow through the control section is

7

critical, the discharge equation of the flow passing through the rectangular opening at

8

the bottom of the open-check dam can be written as:

9

Qb =

10

2 2 1. 5 b0 g hd 3 3

(3)

3.2. Sharp-Crested Weir Flow

11

When hd Ld > 1.5 , water flows through both the bottom rectangular opening

12

and the top trapezoidal spillway of open-check dams. Fig. 6(a) shows the dimensions

13

of a check dam, where b0 is the width of the rectangular opening, b1 is the bottom

14

width of the trapezoidal spillway, b2 is the top width of the trapezoidal spillway, hrd

15

is the water depth of the rectangular opening, and htd is the water depth of the

16

trapezoidal spillway. Considering the kinetic energy of sharp-crested weir flow is

17

high, the water no longer coming into contact with the weir crest and creates a

18

sharp-crested weir water flow pattern. Therefore, the present study used sharp-crested

19

weir theories to predict discharge. The water passes both the rectangular opening and

20

trapezoidal spillway. The linear superposition method (Martínez et al., 2005;

21

Piratheepan et al., 2006; Jan et al., 2005) was adopted to superpose the discharges of

22

the rectangular opening and trapezoidal spillway.

10

1

In Fig. 6(b), Cross Section 1 presents stable water flow in the upstream of the

2

dam, whereas Cross Section 2 shows the water flowing over the weir crest. Assuming

3

that no energy losses occurred between Cross Section 1 to Cross Section 2, the energy

4

relationships between the two cross sections can be obtained as:

5

H1 = (hd − ∆h) +

V22 2g

(4)

6

Assuming that the approaching velocity of progressive flow in Cross Section 1 is

7

zero, then the velocity head in Cross Section 1 (V12 2 g ) can be ignored. The upstream

8

total energy head ( H1 ) is therefore equal to check dam water depth hd ( H1 = hd ),

9

and the flow velocity of Cross Section 2 ( V2 ) can be simplified as:

10

V2 = 2g∆h

(3–5)

11

Using the flow velocity ( V2 ) and cross-sectional areas of Cross Section 2, one

12

can obtain the discharge passing though the weir. Bos (1989) expressed the discharge

13

for flow over a rectangular sharp-crested weir in the following form:

14

2 1.5 Qr = b0 2g hrd 3

15

where Qr is the discharge of the rectangular opening of the open-check dam. For a

16

trapezoidal sharp-crested weir, Henderson (1966) expressed its discharge equation in

17

the form as:

18

Qt =

19

where Qt is the discharge of the trapezoidal spillway. Theoretical discharge equation

20

for an open-check dam with the sharp-crested weir Flow is developed by linearly

(6)

2 4   2 g htd1.5  b1 + mhtd  3 5  

(7)

11

1

superposing the discharge of the rectangular opening and the trapezoidal spillway. The

2

equation can be rearranged as follows:

3

Qc =

 1. 5 2 4   2 g b0 hrd + htd1.5 b1 + mhtd   3 5   

(8)

(a) Front view

(b) Side view

Fig. 6. Definition sketch for sharp-crested weir flow over open-check dam. 4

3.3. Discharge Coefficient ࡯ࢊ

5

In the previous developments of the discharge equations, a number of

6

assumptions are made for the sake of simplicity and therefore a discharge coefficient

7

was usually introduced in the real fluid conditions. (Swame, 1988; Bos, 1989;

8

Rehbock, 1929; Govinda Rao & Muralidhar, 1963; Azimi & Rajaratnam, 2009; Aydin

9

et al., 2011; Salmasi, 2011; Bijankhan et al., 2014; Gogus, 2016). Eqs. (9) and (10)

10

are the theoretical discharge prediction models containing the discharge coefficient

11

for open-check dams under sharp-crested and broad-crested weir flow conditions,

12

respectively.

13

Broad-crested weir flow conditions ( hd Ld ≤ 1.5 )

Q p1 = Cd 1 × Qb 14

2 2 1.5  = Cd 1 ×  b0 g hd  3 3  

(9)

12

1

Sharp-crested weir flow conditions ( hd Ld > 1.5 ) Q p 2 = C d 2 × Qc

2

= Cd 2 ×

 2 4   2 g b0 hrd1.5 + htd1.5 b1 + mhtd   3 5   

(10)

3

where Qp1 is the discharge through the rectangular opening, Qp 2 is the discharge

4

through both the rectangular opening and trapezoidal spillway, Cd1 is discharge

5

coefficient of broad-crested weir flow condition, and Cd 2 is the discharge coefficient

6

of sharp-crested weir flow condition. Aforementioned analysis revealed the discharge

7

coefficients of weir flows, applications of these discharge coefficients have

8

limitations and must be met for the applicable conditions. To effectively and

9

accurately predict the discharge of open-check dams, it is necessary to explore the

10

discharge coefficients under desired flow conditions.

11

4. Experimental Materials and Procedure

12

4.1. Experimental Setup

13

Experiments of this study were performed in a rectangular flume. The flume was

14

supplied by pumping water from the reservoir to the constant-head water tank. An

15

inflow control valve was used to adjust the discharge during the experiments, and the

16

water head of a V-shaped weir was used to determine the discharge. The check dam

17

model was built using 6-mm thick acrylic sheets. The inside of the model was filled

18

with quartz sand to eliminate the effect of buoyancy during the experiments. Silicone

19

was used to seal the contact areas between the flume sidewalls and the bed to prevent

20

experimental errors due to water leakage and dam movement. Fig. 7 shows definition

21

sketch of model dam and flume used in the experiment. Experimental analyses were

22

performed for the five check dams with different opening sizes. The geometric

23

dimensions of all the models are summarized in Table 1. 13

Fig. 7. Definition sketch of model dam and flume used in the experiment 1 Table 1. Parameters of the check dam model Distance between the Experimental

rectangular

group

opening and the riverbed

Rectangular opening width b0

Rectangular

Trapezoidal

Trapezoidal

opening

spillway

spillway

height hr

side slope

height ht

(m)

m

(m)

(m)

Pd (m)

2

Check dam crest length

Ld (m)

E1

0.15

0.05

0.10

0.5

0.15

0.07

E2

0.15

0.075

0.10

0.5

0.15

0.07

E3

0.15

0.10

0.10

0.5

0.15

0.07

E4

0.15

0.05

0.10

0.3

0.15

0.07

E5

0.15

0.05

0.10

1.0

0.15

0.07

4.2. Experimental Procedure

3

During the experiments, it was ensured that the water flow passed through the

4

check dam and became free flow with nappe flow characteristics. A needle gauge with

5

a precision of ±0.1 mm was used to measure the water depth of the flume.

6

Measurements were made from the nappe drop region to the upstream region until the

7

change in water depth was less than 1 mm (i.e., the change in water depth was 14

1

minimal), indicating a stabilized water depth. This water depth of the measurement

2

point was considered as the upstream water head ( h0 ). The check dam water depth

3

hd (= h0 − Pd ) was obtained by subtracting the weir height ( Pd ) from the upstream

4

water head ( h0 ). A total of 110 discharge experiments were performed, in which 22

5

discharge experiments were performed for five check dams. The experimental cases

6

and ranges of variables are summarized in Table 2, where Qm is the discharge, h0

7

is the upstream water head, and Fr is the Froude number based on the upstream

8

water head.

9 Table 2. Flow characteristics of the experimental cases. Experimental Number of

Discharge Qm

Check dam water depth h0

Froude number Fr

group

experiments

(m3/s)

E1

22

0.00023~0.03074

0.169~0.387

0.005~0.131

E2

22

0.00020~0.03416

0.168~0.383

0.005~0.149

E3

22

0.00018~0.02818

0.166~0.384

0.004~0.122

E4

22

0.00052~0.03340

0.171~0.383

0.011~0.145

E5

22

0.00038~0.03375

0.171~0.390

0.008~0.142

(m)

10 11

4.3. Deducing Reduction Factor for Discharge Coefficient

12

The discharge coefficient is generally described as a function of the relevant

13

dimensionless parameters based on experimental data. For a weir, the relevant

14

parameters of discharge ( Qm ) were fluid mass density ( ρ ), gravitational acceleration

15

( g ), dam water depth ( hd ), dynamic viscosity coefficient ( µ ), weir height ( Pd ), and

16

flume slope ( S ) (Aydin et al., 2011). For an open-check dam, the water flowing

17

through the dam may create a contraction effect. Therefore, the present study

15

1

considered the flow area of dam ( A ) as one of the factors that influenced discharge.

2

The discharge can be expressed by the following functional relationship:

3

Qm = f1 (ρ、g、h、 d µ、P、 d S、A)

(11)

4

Because the functional relationship in Eq. (11) represents a physical phenomenon

5

that does not depend on the choice of measurement units, using the Buckingham π

6

theorem, Eq. (11) can be expressed in a dimensionless form as follows:

Qm

 ρ g hd1.5 hd A  = f2  、 、 、S   µ Pd A0  

7

Cd =

8

where B is the width of the channel; A0 ( = B× h0 ) is the cross-sectional area

9

ρ g hd1.5 h upstream of the dam; is the Reynolds number of the flow, d is the ratio µ Pd

B g hd1.5

(12)

A is the ratio between the A0

10

of the check dam water depth and the weir height, and

11

flow area of the dam and that of the upstream. In the experiments flume slope is

12

settled as constant and the research only focuses on the effect of the ratio between the

13

flow area of the dam and that of the upstream on discharge coefficient of open-check

14

dam. The dimensionless parameter for an open-check dam then can be simplified as

15

follows:

16

Cd =

17

where Ar is the ratio between the flow area of the dam and that of the upstream. The

18

use of Ar enables the discharge coefficient to account for the contraction effect and

19

energy changes created by the differences between the cross-sectional areas of the

20

check dam opening and that of the channel.

 A Qm = f 3  1.5 B g hd  A0

  = f 3 ( Ar ) 

(13)

16

1

5. Results and Discussion

2

Figs. 8(a) and (b) show the variations of experimental discharge Qm and the

3

corresponding discharge coefficients Cd1 and Cd 2 at the broad-crested weir flow

4

and sharp-crested weir flow conditions, respectively. For the broad-crested weir flow

5

condition, discharge coefficient increases with increasing discharge or rectangular

6

opening width ( b0 ). For sharp-crested weir flow condition, the behavior of discharge

7

coefficient with discharge is reversed; i.e. discharge coefficient decreases with

8

increasing discharge. As the side slope of spillway m increases, moreover, there is a

9

decrease in discharge coefficient. It is clear from Figs. 8(a) and (b) that the discharge

10

coefficients of open-check dams change with discharge, the geometric shapes of the

11

check dams, and flow conditions change. The discharge coefficient is difficult to

12

represent as a simple trend. An attempt to directly fit the discharge to experimental

13

data with a simple mathematical expression involve using the results of the

14

dimensional analysis.

15

Figs. 9(a) and (b) plot the discharge coefficients

(i.e., Cd1 and C d 2 ) with

16

corresponding dimensionless flow area ratio Ar at different flow conditions,

17

respectively. The importance of flow area ratio Ar for contraction weir flow is

18

obvious. It is possible to develop best fit expressions of discharge coefficients. For the

19

broad-crested weir flow condition, C d 1 and Ar ranged from 1.0 to 1.19 and from

20

0.01 to 0.08, respectively, where Cd1 increases with Ar . The best fit coefficients are

21

obtained by a simple linear regression analysis and the relationship between Cd1 and

22

Ar is given as:

23

Cd1 = 3.278× Ar + 0.958

24

where the coefficient of determination ( R2 ) is 0.87, and discharge coefficient Cd1

25

and dimensionless parameter Ar are highly correlated. When water flow is

(14)

17

1

concentrated in the rectangular opening of the open-check dam, the flow experiences

2

a significant contraction effect and backwater occurs. It increases the potential energy

3

of the upstream water. The increased potential energy is subsequently converted into

4

kinetic energy while passing through the dam. As the potential energy is increased,

5

there is increasing upstream water level which increases the streamline curvature on

6

horizontal plane around the dam. This reduces the flow passing capacity of the weir

7

and corresponding an increase of discharge coefficient. For the sharp-crested weir

8

flow condition, Ar and Cd 2 ranged from 0.06 to 0.26 and from 0.71 to 0.79,

9

respectively, where Cd 2 decreases linearly when Ar increases. A simple linear

10

regression analysis also revealed a high correlation ( R2 = 0.848) between Ar and

11

Cd 2 . The best fit function for Cd 2 and Ar is as follows:

12

Cd 2 = −0.408× Ar + 0.817

(15)

13

When water flooded the spillways, because the geometry of trapezoid is narrow

14

at the bottom and wide at the top, the contraction effect partly vanished with

15

increasing flow depth and side slope (m). Thus, the discharge coefficient decreased

16

with increasing Ar . Figs. 10(a) and (b) show the comparisons of the experimental

17

discharge ( Qm ) and the calculated discharge ( Qp1 and Q p 2 ) obtained by using Eqs.

18

(14) and (15) for discharge coefficient under two flow conditions, respectively. For

19

the broad-crested weir flow condition, the error percentages between the experimental

20

data and predicted value ( (Qm − Q p 1 ) / Qm ) ranged between 0.14% to 2.88% with an

21

averaged value of 1.39%. For the sharp-crested weir flow condition, the error

22

percentages ( (Qm − Q p 2 ) / Qm ) ranged between 0.02% to 1.95% with an averaged

23

value of 0.99%. The prediction models successfully estimate the discharge passing

24

through open-check dams with error less than 3%. 18

1

To evaluate the performance of the discharge prediction model developed in the

2

present study, the results obtained using present model were compared with those

3

obtained using the previous models. For the broad-crested weir flow condition, the

4

previous models of Govinda Rao and Muralidhar (1963), Swamee (1988), Azimi and

5

Rajaratnam (2009), Salmasi (2011), and Gogus (2016) were used. Fig. 11 shows the

6

comparisons of the experimental discharge and the predicted discharge. It is seen that

7

previous models generally underestimated the discharge and the prediction errors

8

increased as discharge increased. The averaged error of different models is

9

summarized in Table 3. The error of the previous models ranges between 4.03% and

10

24.96% and three models have averaged error greater than 10%. From hydraulics

11

point of view, the approaching velocity must be low enough and unsteady phenomena

12

must be eliminated for the purpose of weir flow measurements. Therefore, the

13

contraction rate Ar of the broad-crested weir is sufficient to prevent possible

14

vorticity and turbulence produced at the upstream. In the present cases of weir flow,

15

the contraction effect was pronounced. It resulted in relatively large discharge errors

16

of the previous models than that of the present model. The model proposed by

17

Salmasi (2011) has a largest averaged error of approximately 25%. The discharge

18

prediction models presented by Gogus (2016) and Swamee (1988) displayes relatively

19

favorable prediction results among the previous models. The Gogus’s model shows

20

prediction errors ranging between 0.25% and 11.87% with an average of about 4%.

21

Nevertheless, Gogus’s model showed large deviations at high discharge values which

22

suggests further examinations. The errors of Swamee’s model range from 0.16% to

23

10.45% with an average of about 4%. However, this model involves a complex

24

calculation of parameters that gives a limitation for practically application. Among all

25

the models, the present model provides a lowest averaged error of 1.39%. It is seen

26

that the present model can successfully represent the discharge data with similar 19

1

accuracy at the broad-crested weir flow condition.

2

For the sharp-crested weir flow condition, the present model was compared with

3

the discharge equation proposed by Jan et al. (2005) for the compound sharp-crested

4

weir with similar geometry as the open-check dam. A comparison between the

5

experimental discharge and the corresponding predicted discharge using the present

6

model and the model proposed by Jan et al. (2005) is presented in Fig. 12. The errors

7

of the model proposed by Jan et al. (2005) increases as discharge increases with the

8

averaged and largest errors of approximately 10% and 19%, respectively. Jan et al.

9

(2005) formulated the discharge by dividing into five parts as three rectangular and

10

two triangular sharp-crested weirs and linearly superposed discharge of each part.

11

However, the discharge coefficient of Jan et al. (2005) for the triangular sharp-crested

12

weir considered exclusively as a function of the geometric parameter. The most

13

important characteristic of the weir flow is the independence of the discharge

14

coefficient from the channel width, B, because of the completely contracted nature of

15

the flow over the weir. Errors of the present model fall far below the aforementioned

16

model due to incorporate contracted aspects of flow behavior into the formulation of

17

discharge coefficient. Because the compound cross sections were manufactured using

18

cross sections of varying shapes, the model proposed by Jan et al. (2005) first

19

predicted the discharge and discharge coefficients for each cross sections and linearly

20

superposed discharge of each cross section to obtain the total discharge of the

21

compound cross sections. The discharge coefficient used by Jan et al. (2005) could

22

not account for the changes in the cross sections of where water discharged as check

23

dam water level increased or decreased. Additionally, slit-check dams were exposed

24

to substantial contraction effects because of the changes in the cross sections where

25

water discharged. As a result, the discharge prediction errors created in the model

26

proposed by Jan et al. (2005) were greater than those developed in the present study 20

1

were.

2 3

(a)

(b) 21

Fig. 8. Discharge coefficient versus measured discharge: (a) broad-crested weir flow condition; (b) sharp-crested weir flow condition. 1

(a)

(b) 22

Fig. 9. Discharge coefficient versus flow area ratio: (a) broad-crested weir flow condition; (b) sharp-crested weir flow condition. 1

23

(a)

(b) Fig. 10. Comparison between calculated and measured discharge: (a) broad-crested weir flow condition; (b) sharp-crested weir flow condition. 1

24

Fig. 11. Comparison of predicted discharges by using the various models under broad-crested weir flow condition. 0.04 45 degrees exact line Present study Jan et al. (2005)

0.03

0.02

0.01

0 0

0.01

0.02

0.03

0.04

Measured Discharge (m3/sec) Fig. 12. Comparison of predicted discharges by using the various models under sharp-crested weir flow condition.

25

Table 3. Summary of the errors in the predicted discharges for different models. Error (%) Flow condition

Model Minimum

Maximum

Average

Standard deviation

0.14

2.88

1.39

0.77

4.76

18.15

11.13

3.01

0.16

10.45

4.05

2.49

5.77

19.98

12.48

3.12

17.65

31.20

24.96

3.02

Gogus (2016)

0.25

11.87

4.03

2.69

Sharp-crested Present study weir flow Jan et al. ( hd Ld > 1.5 ) (2005)

0.02

1.95

0.99

0.52

1.43

18.84

9.69

3.84

Present study Govinda Rao & Muralidhar (1963) Swamee Broad-crested (1988) weir flow Azimi & ( hd Ld ≤ 1.5 ) Rajaratnam (2009) Salmasi (2011)

1

6. Conclusions

2

In this study, the discharge equations for open-check dams with a trapezoidal

3

spillway at the top and a rectangular opening at the bottom were developed. The ratio

4

of check dam water depth ( hd ) and crest length ( Ld ) was used to determine the flow

5

conditions of the open-check dams. Theoretical discharge equations were developed

6

for different flow conditions. Considering the contraction effect created by dams, the

7

ratio between the cross-sectional areas of upstream and dam ( Ar ) was incorporated

8

into the discharge coefficient. Experiments were performed to investigate the effect of

9

Ar on discharge coefficients . The predicted discharges of the present study were

10

compared with those of other studies. From the analysis of these results, the following

11

conclusions can be stated: 26

1

1. Formulation of open-check dams depends on the ratio between the check dam

2

water depth and crest length. When hd Ld ≤ 1.5 , the water flow touched the weir

3

crest, creating flow condition that resembled broad-crested weir flow condition.

4

Thus, the broad-crested weir discharge formula can use to predict discharge; when

5

hd Ld > 1.5 , the flow condition resembled sharp-crested weir flow condition, and

6

the sharp-crested weir discharge formula can use to predict discharge.

7

2. Experimental observations indicate the discharge coefficients are significantly

8

affected by Ar and empirical equations of discharge coefficients for different

9

flow conditions were developed, respectively. Under broad-crested weir flow

10

conditions, the averaged error between the predicted and measured discharge was

11

within 2%; under sharp-crested weir flow conditions, the averaged error was

12

within 1%, indicating that the discharge prediction model of the present study

13

produced reasonable results.

14

3. In contrast with other relevant studies, the present model formulates discharge in

15

open-check dams in terms of the water depth and geometric dimensions of

16

open-check dams, which is easy to use practically. For the range of the

17

experiments, it is also demonstrated the present model provides favorable

18

prediction results.

19

7. References

20

1.

Finite Crest Length.” Journal of Hydraulic Engineering, 135(12), 120-125.

21 22

Azimi, A.H. and Rajaratnam, N. (2009). “Discharge Characteristics of Weirs of

2.

Aydin, I., Altan-Sakarya, A.B. and Sismanb, C. (2011). “Discharge formula for

23

rectangular

sharp-crested

weirs.”

24

Instrumentation, 22, 144-151.

Journal

27

of

Flow

Measurement

and

1

3.

Azimi, A.H., Rajaratnam, N. and Zhu, D. (2014). “Submerged Flows over

2

Rectangular Weirs of Finite Crest Length.” Journal of Irrigation and Drainage

3

Engineering, 140(5), 10.1061/(ASCE)IR.1943-4774.0000728, 06014001.

4

4.

Institute for Land Reclamation and Improvement, Wageningen. The Netherlands.

5 6

Bos, M.G. (1989). “Discharge measurement structures”, 2nd Ed., International

5.

Bijankhan, M., Stefano, C.D., Ferro, V. and Kouchakzadeh, S. (2014). “New

7

Stage-Discharge Relationship for Weirs of Finite Crest Length.” Journal of

8

Irrigation

9

10.1061/(ASCE)IR.1943-4774.0000670, 06013006.

and

Drainage

Engineering,

140(3),

10

6.

French, R.H. (1986). “Open-channel hydraulics.”, McGraw-Hill, New York.

11

7.

Fang, C. and Fang, K. (2012). “Derivation and Verification of New Common

12

Formulas for the Calculation of Instantaneous Dam-break Maximum Discharge.”

13

Journal of Procedia Engineering, 28, 648-656.

14

8.

Weirs of Finite-Crest Width.” Houille Blanche, 18(5), 537-545.

15 16

Govinda Rao, N.S. and Muralidhar, D. (1963). “Discharge Characteristics of

9.

Gogus, M., Al-Khatib, I.A., Atalay, A.E. and Khatib, J.I. (2016). “Discharge

17

prediction in flow measurement flumes with different downstream transition

18

slopes.” Journal of Flow Measurement and Instrumentation, 47, 28-34.

19 20

10. Henderson, F.M. (1966). “Open channel flow.” Prentice-Hall, Englewood Cliffs, N.J.

21

11. Jan, C.D., Chang, C.J. and Lee, M.C. (2005). “Discussion of Design and

22

Calibration of a Compound Sharp-Crested Weir.” Journal of Hydraulic

23

Engineering, 131(2), 112-116.

24

12. Martínez, J., Reca, K., Morillas, M.T. and López, J.G. (2005). “Design and

25

Calibration of a Compound Sharp-Crested Weir.” Journal of Hydraulic

26

Engineering, 13(2), 112-116. 28

1

13. Piratheepan, M., Winston, N.E.F. and Pathirana, K.P.P. (2006). “Discharge

2

Measurements in Open Channels Using Compound Sharp-Crested Weirs.”

3

Journal of the Institution of Engineers, 3, 31-38.

4

14. Piton, R. and Recking, A. (2016). “Design of Sediment Traps with Open Check

5

Dams. I: Hydraulic and Deposition Processes.” Journal of Hydraulic

6

Engineering, 142(2), 10.1061/(ASCE)HY.1943-7900.0001048, 04015045.

7 8

15. Rehbock, T. (1929). “Discussion of Precise Weir Measurements.” Journal of Transportation Engineering, 93, 1143-1162.

9

16. Roushangar, K., Akhgar, A., Salmasi, F. (2014). “Estimating Discharge

10

Coefficient of Stepped Spillways under Nappe and Skimming Flow Regime

11

Using Data Driven Approaches.” Journal of Flow Measurement and

12

Instrumentation, 59, 78-87.

13 14

17. Swamee, P.K. (1988). “Generalized Rectangular Weir Equations.” Journal of Hydraulic Engineering, 114(8), 945-949.

15

18. Salmasi, F., Sanaz, P., Dalir, A.H. and Zadeh, D.F. (2011). “Discharge Relations

16

for Rectangular Broad-Crested Weirs.” Journal of Agricultural Science, 17(4),

17

324-336.

18

19. Zollinger, F. (1983). “Die Vorgänge in einem Geschiebeabl Agerungsplatz ihre

19

Morphologie und die Möglichkeiten einer Steuerung.” Ph.D. thesis, ETH Zürich,

20

Zurich.

29

Highlights: •

Discharge prediction model for open-check dams was developed.



Flow conditions are divided as broad-crested weir flow and sharp-crested weir flow.



The results confirm the model was generalized with high accuracy.

Hsun-Chuan Chan: Conceptualization, Methodology. Hsin-Kai Yang: Data curation, Writing- Original draft preparation. Po-Wei Lin: Data curation. Jung-Tai Lee: Supervision, Writing- Reviewing and Editing.