The irrigation sector, being the largest user of water in India, requires 40% more water by 2050 to meet the food demand of the ever-increasing population. Hence, it becomes imperative to implement judicious water management technology to save water and enhance agricultural water productivity. Keeping this in view, the present study developed a portable, low-cost digital openchannel flow measuring device (a modified flume) to measure the discharge rate and volume of irrigation water. The modified flume was designed, fabricated, and tested at the hydraulic flume testing facility of the Water Technology Centre, ICAR-IARI, New Delhi. The device features two configurations: (1) a straight contraction with a curved diverging section and (2) a curved contraction with a straight diverging section. Both configurations were tested with three contraction ratios: 53.33%, 61%, and 70%. It was observed that the device with a 70% contraction ratio straight converging section (30° angles) and a curved diverging section was the best among all experimented designs. Subsequently, a digital sensing system with a waterproof ultrasonic sensor, microcontroller, data logger, and temperature sensor was mounted. The validated rating curve was embedded in the microcontroller through programed code to display the flow depth, discharge rate, and supply volume in a period of time. Nonetheless, the developed device can be tested under different channel shapes and sizes to ascertain its measurement accuracy and wider applicability.

  • A digital portable open channel flow measuring device was developed.

  • The JSN-SR04RT ultrasonic sensor was used for flow depth monitoring.

  • The straight-flow converging section performed better than the curved section.

  • The device measures the flow depth, discharge rate, and flow volume of irrigation water passing in an open channel.

  • Device suitable for earthen and masonry channels.

A good agricultural practice must include the knowledge of water usage by crop and techniques that permit efficient irrigation management. Therefore, the information on the quantity of water supplied to the field would assist in making decisions on when to stop and start irrigation events for different crops (Kumar & Sarangi 2022). When water is over-applied, the excess water percolates below the root zone of the crop without additional benefit to the crop (Abioye et al. 2022). On the other hand, a deficit amount of water, especially during critical growth periods of the crop, can result in the reduction of yield. The measurement of the quantity of water delivered and received by users will ensure that each user gets a fair share, establish a more equitable distribution of available water in the command area, and promote conservation of this precious resource. Monitoring how much water is being delivered in real-time with measuring devices can help farmers or canal authorities deliver the required amount of water to the fields according to crop growth stages and improve water productivity (Kumar & Sen 2020; Kumar et al. 2023). Moreover, for the development of an efficient on-farm water management plan, it is necessary to quantify how much water is being delivered to fields. Accurate water measurement provides the on-farm irrigation decision-maker with the information needed to achieve the best use of the irrigation water applied while minimizing negative environmental impacts. The amount of water applied to a field is a function of time, flow rate, and area. The area to be irrigated and the time of supply of irrigation water can be easily recorded, whereas the flow rate in an open channel requires appropriate instruments or methods for its measurement. Accurate flow measurement of an open channel can prove to be challenging because the flow is not under pressure, and pipe measuring devices such as venturi, electromagnetic, or strap-on transit-time flow meters are not viable.

Measuring water in an open channel is an important step toward water conservation (Samani et al. 1991; Kumar & Sarangi 2022; Kumar & Sen 2023). With the increasing demand for improved water management techniques around the world, there is an urgent need for low-cost and accurate flow measuring devices. Open channel flow measurement devices do not measure the flow directly. Instead, some devices measure the velocity of the flow, and others measure changes in head (flow depth) or pressure. Among the major types of measurement devices used in surface irrigation (open channels) are weirs, flumes, current meters, orifices, propeller meters, venturi meters, electromagnetic meters, turbine meters, ultrasonic meters, and pitot tubes (Santhosh & Roy 2012). Open channel flow measuring devices employ the principles of either the weir or the orifice, and each device is adapted to be used in specific locations (Kumar 2020; Kumar & Sarangi 2022). Moreover, the ideal measuring device would be inexpensive to construct, simple to operate and free from working parts, require little maintenance, provide accuracy in its measurement, not be affected by sand, silt, or floating trash and require minimal loss of head in the channel. Several design and development efforts for standardizing different dimensions of measuring devices in open channels have been made, viz. venturi flume (Cone 1917), Parshall flume (Parshall 1928), trapezoidal flume (Hyatt 1965), cut-throat flume (Skogerboe 1969), circular flume (Hager 1988), simple flume (Samani & Magallnez 2000). Ever since the development of the Parshall flume (Parshall 1928), attempts have been made to simplify its construction, improve the accuracy, and reduce the cost of open channel flow measuring devices. Parshall flume, cut-throat flume, trapezoidal flume are the most commonly used devices particularly in flat topography and with small channel gradients. The discharge becomes a function of only a single flow depth which is measured at some location upstream from where the critical depth/flow occurs. Critical flow in an open channel can be created through three general methods: raising the bottom of the channel (Replogle 1975; Bos et al. 1984), contracting the cross-sectional area of the flow (Skogerboe et al. 1967a, 1972; Hager 1988; Samani et al. 1991; Samani & Magallnez 1993, 2000), and reducing the bottom elevation to create a critical flow. Contracting the cross-section of the flow is one of the simplest methods of creating a critical flow for measuring the flow in an open channel (Samani et al. 1991; Hager 1988; Samani & Magallnez 2000). These attempts resulted in the development of the cut-throat flume (Robinson & Chamberlain 1960; Skogerboe et al. 1967b) and the Replogle, Bos, and Clemmens flume (Bos et al. 1984). Measuring flow by contracting the flow cross-section is often the simplest and cheapest method, as it does not require complex inflow and outflow transition (Samani 2017). The ideal condition for the measurement of discharge is that the throat section should be sufficiently constricted to produce critical depth. Therefore, most of the measurement structures utilize the principle of passing the flow through critical depth, leading to free flow conditions in which the discharge depends only on the upstream flow depth. However, when the downstream depth of flow increases to the point that the upstream depth of flow is affected, then free flow conditions cease, and the flow is said to be submerged. When submerged flow conditions exist, the stage and discharge relationship developed for free flow conditions is no longer valid (Hyatt 1965). The term free flow is often used when critical depth occurs in a flume. Flumes designed to operate under free flow conditions will be submerged sometimes, due to either unusual operating conditions downstream or the accumulation of mass and/or vegetation in the channel. Details of existing hydraulic measurement methods, such as weirs, orifices, calibrated flumes, and other indigenous flow measuring devices, along with their drawbacks, limitations, and future scope, are provided in Kumar & Sarangi (2022).

The open channel flow measuring device serves as the backbone of water measurement networks for open channel irrigation to enhance surface irrigation efficiency (Goel et al. 2015). Different kinds of flumes, viz. Parshall flume, cutthroat, and long-throated flume, are being used in India and abroad for the measurement of flow in open channels. But these flumes are associated with many known and unknown errors. The major lacunae in the use of the Parshall flume are the difficulty in the configuration of the throat section, its sloping floor, and field installation. Furthermore, a major cause for error with Parshall flumes is the submergence due to the backwater effect. Besides this, the structural settlement of the Parshall flume causes errors in flow measurement (Skogerboe & Hyatt 1967; Abt & Staker 1990; Prasad 1991; Kruse 1992; Genovez et al. 1993; Abt et al. 1989, 1994; 1995; 1998). Furthermore, the incorrect entrance geometry and gauge location cause an error in flow measurement using the Parshall flume. Moreover, efforts have been made by researchers around the globe to modify the design of existing flumes and to develop a flume that would overcome the limitations of the Parshall flume.

A meticulous review of literature pertaining to flow-measuring devices and the digital sensing system for open field channels revealed the non-availability of a portable and digital flow-measuring device for surface irrigation. Keeping in view the above, an effort was made in this study to develop a low-cost digital flow-measuring device for open channels to obviate the majority of the limitations pertaining to the use of weirs, orifices, calibrated flumes, or other indigenous flow-measuring devices.

The experiment was undertaken at the hydraulic flume testing facility of Water Technology Centre, ICAR-IARI, New Delhi, using the developed flow measuring device to determine its best hydraulic dimensions. The detailed descriptions of the principle of discharge measurement in open channel flow, different components of the flow measuring device, and the digital sensing system, besides the calibration and testing of the developed device, are presented in the following sections.

Principle of discharge rate measurement in the open channel

The discharge rate in the open channel can be measured by multiplying the flow velocity with the wetted area of the flowing water in the channel. To accomplish this, the flow velocity at a given flow depth and the corresponding wetted area in the channel need to be measured as and when there is a change in the flow depth. Moreover, such an approach in an open channel becomes cumbersome and necessitates the periodic measurement of velocity and the corresponding wetted area of the channel. Besides this, the friction of the channel surface, surface tension, and aerodynamic pressure cause instability of the water surface and pose hindrances in quantifying the wetted area accurately on real-time basis. Therefore, the measurement of discharge at a particular section of an open channel is undertaken using a calibrated depth vs. discharge relationship through the development of a rating curve, which would give accurate discharge estimation until the shape of the channel remains intact and the free flow condition prevails (Parsaie et al. 2017; Kumar & Sen 2024). The field channels may be made up of brick- cement or earth. Limitations with these channels are that the masonry channel constructed in the field is of a constant shape without any section along the flow length to ensure a static (stable) depth of flow, which restricts its calibration, whereas the shape of earthen channels varies, which also requires frequent generation of the rating curve for accurate flow measurement. These limitations can be obviated by the installation of a calibrated permanent or portable structure with the rating curve to measure the discharge in the open channel. The basic principle in the design of the primary device (portable structure) was to create a critical flow condition by either contracting some part of the channel section,changing the bottom slope, raising the channel bottom or a combination of all these three approaches in the measuring devices (Skogerboe et al. 1967a, 1972; Replogle 1975; Bos et al. 1984; Hager 1988; Samani et al. 1991; Samani & Magallnez 1993, 2000). Hence, the occurrence of critical flow permits the application of hydraulic principles to determine the unique relationship between the head (depth of flow) and discharge rate with high precision and accuracy. In this invention, the primary flow measuring device was based on the contraction of the central part of the device, the level of the bottom floor, and was designed without a throat section to create the critical flow. Furthermore, a digital flow depth measuring ultrasonic sensor was integrated with the device, and the head vs. discharge equation (rating curve) was embedded in the microcontroller to measure the discharge rate instantaneously and the total volume of water passing through the device over a period.

Design the shape and dimensions of the open channel flow measuring device

The design of open channel flow measuring devices is primarily intended for use in field channels for surface irrigation. ‘Field channel’ is the open channel used for supplying water to farmers' fields in a command area and is fed by the water courses under the canal command network. Field channels in canal commands are designed to carry about 1 cusec (28.3 L s−1) of water, which is fed from water courses through adjustable proportionate modules. In India, a field channel is generally either an earthen or masonry (brick-cement) structure. Keeping this in view, the measuring device was designed to fit into the existing field channel in any canal command in India. However, the device can also be used in any open channel with minor modification of the channel shape or preparation of an approach section to fit the structure with minimal effort. The bottom width of the flow measuring device was determined based on the principle of the most economical rectangular channel to carry the maximum discharge rate in the field channel (28.3 L s−1). Based on this, the bottom width was designed to be 0.30 m (30 cm). The length of the device (85 cm) was decided from the experiment to ensure that the smooth transition of flow from the channel to the measuring the device and from the device to the downstream channel occurs.

The shape of the designed flow measuring device is based on hydraulic principles, i.e. the downstream water does not affect the upstream water depth up to a certain limit. This can be achieved by raising the bed, lowering the flow bed, reducing the cross-section of the channel, or a combination of these. Raising the flow bed causes sedimentation upstream leading to errors in flow measurement, and lowering induces structural complexity and installation problems. The designed device has a contracted cross-section at the central part of the device with a flatbed and without a throat section. The shape of contractions is designed in such a manner as to provide a smooth and parallel streamline without the effects of flow separation in the device. Besides this, it also ensures the occurrence of critical flow in the downstream section, thereby creating a tranquil flow with minimal celerity and stable depth of flow in the upstream section for the more accurate measurement of depth and subsequent discharge passing through the device. The contractions are chosen to avoid overtopping of water in the side walls of the flume or from the upstream channel. Two types of shapes converging and diverging sections with three contraction ratios are designed and tested to select the best shape and dimension of the primary flow-measuring device. The first one had a curve conversing and a straight diversion section (30° and 45°) (Figure 1(a)) and another straight conversing section (30° and 45°) and a curve diverging (Figure 1(b)). A detailed description of the different components and dimensions of the flow measuring device is presented in Figure 1 and Table 1. This device has four sections, viz. (i) an approach section, (ii) a converging section, (iii) a diverging section, and (iv) an outlet section. The bottom floor of the device is level. The rectangular approach and downstream section ensure a smooth transition of flow into the converging section and in the downstream channel, especially in the earthen channel. The shape of contractions is designed in such a manner as to provide a smooth and parallel streamline without the effects of flow separation in the device. The upstream approach and downstream segments are designed to ensure the smooth transition of water from the upstream channel to the contraction section and from the contraction section to the downstream channel.
Table 1

Dimension of flow measuring device curved entrance section and straight diverging section

Contraction ratio (%) = (B − 2R)/B × 10053.336170
Figure 1(a)  Diverging angle (θ, degree) 45° 45° 45° 
  30° 30° 30° 
Figure 1(b)  Conversing angle (θ, degree) 45° 45° 45° 
  30° 30° 30° 
Contraction ratio (%) = (B − 2R)/B × 10053.336170
Figure 1(a)  Diverging angle (θ, degree) 45° 45° 45° 
  30° 30° 30° 
Figure 1(b)  Conversing angle (θ, degree) 45° 45° 45° 
  30° 30° 30° 
Figure 1

The top view of the experimented device. la = length of the approach section, lb = length of the conversing section, lc = length of the diverging section, ld = length of the outlet section, le = location of the ultrasonic sensor in the approach section, B = width of the open channel (30 cm), R = radius of the curve section, θ = angle of straight converging or the diverging section.

Figure 1

The top view of the experimented device. la = length of the approach section, lb = length of the conversing section, lc = length of the diverging section, ld = length of the outlet section, le = location of the ultrasonic sensor in the approach section, B = width of the open channel (30 cm), R = radius of the curve section, θ = angle of straight converging or the diverging section.

Close modal

Digital flow sensing system

Measurement of flow in an irrigation channel either at a given time or on a continuous basis is an essential requirement not only to prevent wastage of water but also to increase water use efficiency at the farm level. A typical continuous measuring system consists of a primary flow device, a flow sensor, a transmitter, a flow recorder, and a totalizer. Instantaneous flow measurements can be obtained using the primary flow device, which is the heart of a typical continuous flow measurement system that basically produces predictable hydraulic responses about the depth of flow passing through it. A flow sensor is required to measure the hydraulic responses of the primary flow measurement device and transmit the responses to the recording system. In general, sensors include ultrasonic transmitters, floats, pressure transducers, capacitance probes, differential pressure cells, electromagnetic cells. A microcontroller-based flow depth sensing system with a data logger was fabricated for the measurement and recording of flow parameters (depth of flow, discharge rate, total amount delivered, required amount, and remaining amount) in open field channel irrigation. The components of the sensing system are Arduino, ultrasonic sensors (JSN-SR04T), a secure digital (SD) card module with an SD card, a real-time clock (RTC) module, and a liquid crystal display (LCD), as shown in Figure 2. The Arduino Mega is an open-source microcontroller programed with the Arduino integrated development environment software. The C-based programming language was used as code for the Arduino operation. The ultrasonic sensor is a non-contact distance measurement sensor. It is used for the instantaneous measurement of the depth of flow in the approach/upstream section and the corresponding discharge rate and quantification of the delivered amount of water. A micro-SD card reader module is used to read and write different data files. A RTC module is used to keep track of local dates and times. An LCD screen displays instantaneous flow parameters and other ranges of applications.
Figure 2

Circuit diagram of thedeveloped flow-measuring sensing system.

Figure 2

Circuit diagram of thedeveloped flow-measuring sensing system.

Close modal

Measurement of discharge corresponding to a given flow depth passing through the measuring device

The experiment was conducted with a flow rate varying from 1.5–20 L s−1, which was in line with the field channel in a canal command area. To experiment with the above-mentioned flow range, the discharge rates at various depths of flow were measured volumetrically. The depth of flow in the approach section was regulated by the gate valve to bypass the excess flow rates to the storage tank. At a given flow rate, a fixed volume of water was collected in the sample collecting tank, and the time taken for such collection of a fixed volume was recorded. Subsequently, the discharge rate (L s−1) through the flume was calculated by dividing the fixed volume (L) collected in the sample tank by the collection time (s). The same procedure was repeated three times for the measurement of each flow depth to avoid errors in measurement, and the average discharge was calculated for each flow depth passing through the measuring device. The occurrence of critical depth was measured after averaging the discharge for a particular depth of flow in the diverging section, and the location of critical depth was measured and noted.

Computation of critical depth

A theoretical criterion for critical flow is the state of flow in which the specific energy (H) is the minimum for a given discharge (Q).

From the continuity equation for rectangular channels (Chow 1988),
(1)
where Q is the volume rate of flow, Ls−1; V is the average velocity, ms−1; Vc is the critical velocity, m/s; A is the cross-sectional area of flow, m2; Ac is the cross-sectional area at critical flow, m2;
From Bernoulli's equation
(2)
where H is the total energy head, m; g is the acceleration due to gravity, ms−2; y is the depth of flow, m.
Substitution of (1) in (2) gives
(3)
The water area dA near the surface is equal to Bdy, then the equation becomes
At the critical state of flow, the above equation, therefore, gives
(4)
where Ac = BcYc, substituting in (4)
(5)
where the subscript ‘c’ relates to the critical condition. Froude number is at critical flow.

Development of flow depth and the discharge relationship of measuring devices

The flow depth vs. discharge relationship of different designed devices was developed. The relationship should be in the form of an equation, which should indicate no discharge under zero depth of flow to corroborate the real flow situation in any open channel. Therefore, the linear or quadratic form of the equation or equations with a constant that is not associated with an independent variable (flow depth) cannot be considered in the development of the rating curve. Hence, the following power equation, was used for the development of the flow depth versus discharge relationship.
(6)
where, Qis the discharge rate in L s−1; Ha is the depth of flow in the upstream section, cm; αis the coefficient; and β is an exponent.

The most appropriate values of a and b were determined by the nonlinear least square grid search (NLIN) procedure using the Statistical Analysis System (SAS, version-9.4) software Copyright (c) 2002–2012 by SAS Institute Inc., Cary, NC, USA. The PROC NLIN procedure performs univariate nonlinear regression using the least squares methods. Fitting nonlinear models can be a difficult task. There is no closed-form solution for the parameter estimates, and hence the process is iterative. This procedure solves the nonlinear regression problem by one of the following five algorithms: (i) the steepest-descent or gradient method, (ii) the Newton method, (iii) the modified Gauss–Newton method, (iv) the modified Gauss–Newton method, and (v) the Marquardt method. From the above methods, the Marquardt method used by the NLIN procedure is considered for the present investigation to estimate the coefficient and exponent having the smallest sum of squared errors (SSE). The initial values of ‘α’ and ‘β’ for the NLIN procedure were determined from the power function developed using the observed data of different flow depths passing through the device and corresponding discharge rates. Then, to validate the developed empirical models, performance indicators like mean absolute deviation (MAD), mean square error (MSE), root mean square error (RMSE), and mean absolute percentage error (MAPE) are used with actual discharge rate and estimated discharge rate.

Performance evaluation of designed flow measuring devices

Fabricated flow measuring devices with different dimensions of the contraction and diverging section besides contraction ratios (Table 1) were tested to select the best-designed flow measuring device. It was observed during the experiment that at the contraction section, a hydraulic drop in the water surface occurred due to side contractions causing appreciable curvature of the flow profile accompanied by an abrupt acceleration of flow. The flow profile in the contraction section of the device depends upon the discharge and side contractions. In general, the basic principle in the design of an open channel flow measuring device is to ensure the occurrence of critical flow in the downstream section, thereby creating a tranquil flow with minimal celerity and stable depth of flow in the upstream section for the more accurate measurement of depth and subsequent discharge passing through the device. Keeping this in view, the occurrence of critical depth for different designed dimensions was investigated. In general, the critical depth occured at different locations corresponding to different side contractions and different shapes of converging and diverging sections of designed dimensions. Also, the flow in the contraction and downstream sections of the device was more abrupt with an increase in the contraction ratio. The designed devices have curved converging and straight diverging sections with 30° and 45° angles, as shown in Figure 1(a) and vice versa with a straight converging section (i.e. with 30° and 45° angles), as in Figure 1(b) was experimented with under different contraction ratios and varying discharge rates. The occurrence and location of critical depth from the transition point of converging and diverging sections for both devices are presented in Table 2. It was observed from Table 2 that the contraction ratio played a major role in the occurrence of critical flow, and with lower contraction ratios, i.e. 53.33 and 60%, the critical flow occurred at a higher discharge rate of 4.5 and 2.5 L s−1 and above, respectively. Whereas with a contraction ratio of 70%, the critical flow occurred early at 1.5 L s−1. Besides this, the occurrence of critical depth was not affected much due to the different shapes of converging and diverging sections. However, with the curved converging section and straight diverging section, the location of the critical depth varied from 7.3 cm to 12.5 cm, corresponding to 70 and 53.33% contraction ratios, respectively. Moreover, for the straight converging section and curved diverging section, the location of critical depth was still lower, ranging from 6.3 to 10.2 cm apropos to 70 and 53.33% contraction ratios, respectively. The distance of critical flow from the transition point was still lower with a straight converging section of 45° angle (i.e. 6.3 cm) than that of 30° angle (i.e. 6.5 cm). It can be concluded that critical flow occurred for lower flow depths with a higher contraction ratio (70%) than lower contraction ratios (53.33 and 61%). Such flow behavior in the open channel may be the lower contraction ratio does not create sufficient head in the converging section to create a critical flow in the downstream section.

Table 2

Occurrence and location of critical depth with varying shapes of converging and diverging sections

Contraction ratio (%)Converging sectionDiverging sectionDischarge rate for the occurrence of critical flow (L s−1)Location of critical depth (cm) from transition point of converging and diverging sections
53.33 Curve 45° ≥4.5 9.4–12.3 
Curve 30° ≥ 4.5 10.3–12.5 
45° Curve ≥ 4.5 7.5–10.3 
30° Curve ≥ 4.5 7.3–10.2 
61 Curve 45° >2.5 7.8–10.2 
Curve 30° >2.5 8.6–10.4 
45° Curve >2.5 7–9.6 
30° Curve >2.5 6.8–9.5 
70 Curve 45° 1.5–20 7.3–9.5 
Curve 30° 1.5–20 8–9.7 
45° Curve 1.5–20 6.3–8.8 
30° Curve 1.5–20 6.5–9 
Contraction ratio (%)Converging sectionDiverging sectionDischarge rate for the occurrence of critical flow (L s−1)Location of critical depth (cm) from transition point of converging and diverging sections
53.33 Curve 45° ≥4.5 9.4–12.3 
Curve 30° ≥ 4.5 10.3–12.5 
45° Curve ≥ 4.5 7.5–10.3 
30° Curve ≥ 4.5 7.3–10.2 
61 Curve 45° >2.5 7.8–10.2 
Curve 30° >2.5 8.6–10.4 
45° Curve >2.5 7–9.6 
30° Curve >2.5 6.8–9.5 
70 Curve 45° 1.5–20 7.3–9.5 
Curve 30° 1.5–20 8–9.7 
45° Curve 1.5–20 6.3–8.8 
30° Curve 1.5–20 6.5–9 

Overall, based on the occurrence of critical flow and its location in the downstream section pertaining to different designed dimensions of the flow measuring device, it was indicated that the device with a 70% contraction ratio having a straight converging section with either a 30° or 45° angles and a curved diverging section was better than other experimented designs. The criterion of such selection was that the device could create a stable depth (sub-critical depth) in the upstream section under a wide discharge range besides the lowest discharge rate. Moreover, the critical flow occurred in the downstream section, and the length of such occurrence was minimal, which would lead to the fabrication of a smaller device, leading to savings in material, cost, and ease of portability. Further, the flow depth and discharge rate relationship of different designs without throat sections was developed to ascertain the best device that can measure the discharge more accurately.

Flow depth vs. discharge relationship of the device with a contraction ratio of 53.33%

Devices with different contraction ratios, besides varying shapes of converging and diverging sections, as mentioned in Table 1, were used in the experiment to acquire data on different flow depths and the corresponding discharge rate. These values were plotted with the dependent variable (i.e. discharge rate in L s−1) against the independent variable (i.e. flow depth in cm), and the regression model in the form of a power equation (Equation (6)) was developed. The fitted trend lines with power equations of different shapes for a contraction ratio of 53.33% are presented in Figure 3(a) and 3(b). The coefficient of determination (R2) of all four fitted equations ranged 0.9825–0.9981. The values of ‘α’ and ‘β’ in all design combinations as presented in Figures 3(a) and 3(b) were further used as the input values in the statistical procedure NLIN to estimate the optimal values of ‘α’ and ‘β’ with minimized residual sum of squares. The values of ‘α’ and ‘β’ generated using the NLIN procedure are presented in Table 3. It is reported (Hougaard 1985; Ratkowsky 1990) that the skewness of ‘α’ and ‘β’ parameters, i.e. γα and γβ should be within a range of 0.1–0.25 for better prediction accuracy of the developed model. Moreover, Table 3 shows that the skewness ranged from 0.1 to 0.25 for all parameters barring for ‘α’ under curved converging and straight diverging with 45° angle. Subsequently, ‘α’ and ‘β’ parameters obtained from the NLIN procedure and presented in Table 3 were used in the development of four equations, which are presented in the following:
Table 3

Input and optimal values of ‘α’ and ‘β’ for the designed device

CR (%)CSDSInput values
Optimal values
SSESkewness
αβαΒγαγβ
53.33 Curve 45° 0.3436 1.4202 0.2499 1.5723 0.7275 0.2723 0.0410 
Curve 30° 0.2822 0.15036 0.2723 1.5193 0.2210 0.1707 0.0235 
45° Curve 0.1876 1.6464 0.1815 1.6615 0.3136 0.2351 0.0270 
30° Curve 0.1617 1.7077 0.1503 1.7416 0.3181 0.1508 0.0184 
61 Curve 45° 0.2414 1.5172 0.2327 1.5338 0.2929 0.1171 0.0217 
Curve 30° 0.2126 1.5581 0.2127 1.5581 0.4528 0.1304 0.0244 
45° Curve 0.1877 1.5922 0.1516 1.6852 0.5188 0.1644 .0250 
30° Curve 0.1637 1.6418 0.1712 1.6230 0.2527 0.1133 0.0076 
70 Curve 45° 0.1790 1.6123 0.1871 1.5945 0.6431 0.2329 0.0214 
Curve 30° 0.0991 1.8472 0.1131 1.7917 0.2985 0.1664 0.0147 
45° Curve 0.1912 1.5248 0.1614 1.5952 0.4032 0.1655 0.0244 
30° Curve 0.1121 1.7210 0.1119 1.7219 0.8247 0.1594 0.0162 
CR (%)CSDSInput values
Optimal values
SSESkewness
αβαΒγαγβ
53.33 Curve 45° 0.3436 1.4202 0.2499 1.5723 0.7275 0.2723 0.0410 
Curve 30° 0.2822 0.15036 0.2723 1.5193 0.2210 0.1707 0.0235 
45° Curve 0.1876 1.6464 0.1815 1.6615 0.3136 0.2351 0.0270 
30° Curve 0.1617 1.7077 0.1503 1.7416 0.3181 0.1508 0.0184 
61 Curve 45° 0.2414 1.5172 0.2327 1.5338 0.2929 0.1171 0.0217 
Curve 30° 0.2126 1.5581 0.2127 1.5581 0.4528 0.1304 0.0244 
45° Curve 0.1877 1.5922 0.1516 1.6852 0.5188 0.1644 .0250 
30° Curve 0.1637 1.6418 0.1712 1.6230 0.2527 0.1133 0.0076 
70 Curve 45° 0.1790 1.6123 0.1871 1.5945 0.6431 0.2329 0.0214 
Curve 30° 0.0991 1.8472 0.1131 1.7917 0.2985 0.1664 0.0147 
45° Curve 0.1912 1.5248 0.1614 1.5952 0.4032 0.1655 0.0244 
30° Curve 0.1121 1.7210 0.1119 1.7219 0.8247 0.1594 0.0162 

Note: CR, contraction ratio (%); CS, converging section; DS, diverging section; α, coefficient of rating curve; β, exponent of rating curve; SSE, sum of square error; γα, skewness of coefficient; γβ, skewness of exponent.

Figure 3

(a) Developed flow depth vs. discharge relationship for curved converging and straight diverging sections with a 45° (a) and 30° (b) angle. (b) Developed a flow depth vs. discharge relationship for the straight converging section with a 45° (a) and 30° (b) angle and a curved diverging section.

Figure 3

(a) Developed flow depth vs. discharge relationship for curved converging and straight diverging sections with a 45° (a) and 30° (b) angle. (b) Developed a flow depth vs. discharge relationship for the straight converging section with a 45° (a) and 30° (b) angle and a curved diverging section.

Close modal
Curved converging section with straight diverging (45° angle) section:
(7a)
Curved converging section with straight diverging (30° angle) section:
(7b)
Straight converging (45° angle) section with curved diverging section:
(7c)
Straight converging (30° angle) section with curved diverging section:
(7d)

Flow depth vs. discharge relationship of the device with a contraction ratio of 61%

Similarly, the flow depth versus discharge rate relationship was developed for all designed shapes of the device with a contraction ratio of 61%, as mentioned in Table 1 and presented in Figures 4(a) and 4(b). The R2 value of all fitted equations ranged from 0.9947 to 0.9987. The values of ‘α’ and ‘β’ pertaining to four design combinations, as presented in Figures 4(a) and 4(b), were subjected to NLIN analysis, and the optimized ‘α’ and ‘β’ values are presented in Table 3. It was observed from Table 3 that the skewness γα and γβ ranged from 0.0076 to a maximum of 0.1644 and were well within <0.25, which indicated that the developed model would be in line with the data acquired from the experiment. Subsequently, ‘α’ and ‘β’ parameters obtained from the NLIN procedure, as presented in Table 3, were used for the development of the flow depth versus discharge relationship, which are presented in the following:
Figure 4

(a) Developed flow depth vs. discharge relationship for curved converging and straight diverging sections a with divergence angles of 45° (a) and 30° (b). (b) Developed a flow depth vs. discharge relationship for the device with a straight converging section with angles of 45° (a) and 30° (b) and a curved diverging section.

Figure 4

(a) Developed flow depth vs. discharge relationship for curved converging and straight diverging sections a with divergence angles of 45° (a) and 30° (b). (b) Developed a flow depth vs. discharge relationship for the device with a straight converging section with angles of 45° (a) and 30° (b) and a curved diverging section.

Close modal
Curve converging section with straight (45° angle) diverging section:
(8a)
Curve converging section with straight (30° angle) diverging section:
(8b)
Straight (45° angle) converging section with curved diverging section:
(8c)
Straight (30° angle) converging section with curved diverging section:
(8d)

Flow depth vs. discharge relationship of the device with a contraction ratio of 70%

A similar procedure was adopted to develop the flow depth versus discharge relationship of the designed device with a contraction ratio of 70% having different shapes of converging and diverging sections. A fitted trend line with parameters ‘α’ and ‘β’ of the power equation is shown in Figures 5(a) and 5(b). It was observed from these figures that R2 of the fitted model ranged from 0.9941 to 0.9984. Subsequently, the iterated optimal values of ‘α’ and ‘β’ along with prediction error statistics corresponding to the NLIN procedure are presented in Table 3. It is observed from Table 3 that the skewness (γα and γβ) pertaining to parameters ‘α’ and ‘β’ ranged from 0.0147 to 0.2329, and hence all values were less than 0.25, corroborating that the developed model with parameters ‘α’ and ‘β’ would represent the real experimental situation with added accuracy. Subsequently, four equations were developed for a 70% contraction ratio using ‘α’ and ‘β’ parameters, which are given as follows:
Figure 5

(a) Developed flow depth vs. discharge relationship of the device with curved converging section and straight diverging sections with divergence angles of 45° (a) and 30° (b). (b) Developed a flow depth vs. discharge relationship of the device with a straight converging section with 45° (a) and 30° (b) angles and a curved diverging section.

Figure 5

(a) Developed flow depth vs. discharge relationship of the device with curved converging section and straight diverging sections with divergence angles of 45° (a) and 30° (b). (b) Developed a flow depth vs. discharge relationship of the device with a straight converging section with 45° (a) and 30° (b) angles and a curved diverging section.

Close modal
Curved converging section with straight diverging (45° angle) section:
(9a)
Curved converging section with straight diverging (30° angle) section:
(9b)
Straight converging (45° angle) section and curved diverging section:
(9c)
Straight converging (30° angle) section with curved diverging section:
(9d)

Overall, it was observed from the flow depth vs. discharge relationship of 12 sets of designed devices with three contraction ratios that 11 sets performed well as per the prediction error statistics. However, in the device with a contraction ratio of 53.33% and having curved converging and straight diverging sections with 45° angle, the coefficient ‘α’ exceeded the limiting skewness value of 0.25 and hence cannot be considered the best representative model simulating the real field situations. Further, the developed rating curves of all designed devices were validated under variable flow depth situations to arrive at the best-designed dimension of the flow measuring device.

Validation of developed flow depth vs. discharge relationship (rating curve) of designed flow measuring devices

Accurate measurement of flow depth in an open channel on a real-time basis is the prime factor leading to more accurate quantification of discharge by the flow measuring device. Moreover, the depth of flow in an open channel is a highly variable phenomenon, and properly calibrated sensors would capture an accurate flow depth in the channel. The developed flow measuring device, mounted with a digital sensing system, was validated under varying flow depths in the open channel. The flow depth was measured with MAD; 0.21 ± 0.01 cm, by the JSN-SR04T ultrasonic sensor (Kumar et al. 2023), and the discharge rate was displayed in the digital system as per the rating curve equations embedded in the microcontroller. However, the actual discharge was measured in the sample collection tank. The average discharge rate of five successive readings was taken as recorded in the flow sensing system. Prediction error statistics of validation results of all experimented devices were compared to select the best design of the flow measuring device.

Validation of models of contraction ratio 53.33% and different shapes of converging and diverging sections

Flow depth vs. discharge Equations (7a)–(7d) developed for different shapes of the converging and diverging sections, as detailed in section 3.1.1, were embedded in the digital sensing system pertaining to specific devices. These devices were experimented with in the hydraulic flume testing facility. The prediction error statistical parameters MAD, MSE, RMSE, and MAPE of all designed and tested devices with a contraction ratio of 53.33% are presented in Table 4. It is observed from Table 4 that the MAD, MSE, RMSE, and MAPE were lowest for curved-d converging sections with straight diverging sections with 45° angle and maximum for straight converging sections with 30° angle with curved diverging sections. However, the MAPE of all models was greater than 5, which indicated that the developed device with the rating equation will not perform at its best to simulate real field situations.

Table 4

Validation error statistics of flow measuring devices with varying shapes and contraction ratios

Contraction ratio (%)Converging sectionDiverging sectionMADMSERMSEMAPE
53.33 Curve 45° 0.31 0.10 0.31 5.95 
Curve 30° 0.49 0.27 0.52 8.27 
45° Curve 0.33 0.11 0.33 6.63 
30° Curve 0.44 0.20 0.45 8.02 
61 Curve 45° 0.50 0.26 0.52 8.52 
Curve 30° 0.41 0.17 0.41 6.99 
45° Curve 0.35 0.14 0.38 6.37 
30° Curve 0.31 0.11 0.33 5.27 
70 Curve 45° 0.32 0.11 0.33 5.92 
Curve 30° 0.30 0.11 0.34 4.57 
45° Curve 0.37 0.14 0.38 7.56 
30° Curve 0.29 0.09 0.30 4.33 
Contraction ratio (%)Converging sectionDiverging sectionMADMSERMSEMAPE
53.33 Curve 45° 0.31 0.10 0.31 5.95 
Curve 30° 0.49 0.27 0.52 8.27 
45° Curve 0.33 0.11 0.33 6.63 
30° Curve 0.44 0.20 0.45 8.02 
61 Curve 45° 0.50 0.26 0.52 8.52 
Curve 30° 0.41 0.17 0.41 6.99 
45° Curve 0.35 0.14 0.38 6.37 
30° Curve 0.31 0.11 0.33 5.27 
70 Curve 45° 0.32 0.11 0.33 5.92 
Curve 30° 0.30 0.11 0.34 4.57 
45° Curve 0.37 0.14 0.38 7.56 
30° Curve 0.29 0.09 0.30 4.33 

Validation of a model of a contraction ratio of 61% and different shapes of converging and diverging sections

Flow depth vs. discharge equations (i.e. Equations (8a)–(8d)) developed in section 3.1.2 were embedded in the microcontroller of the digital sensing system and used in the experimental setup for subsequent validation. The same procedure was followed, as mentioned in section 3.2.1. Referring to Table 4, the MAPE pertaining to the validation results of all four devices was greater than 5.0, which also raised an issue on its accuracy on subsequent use for flow measurement.

Validation of device with a contraction ratio of 70% and different shapes of converging and diverging sections

Similarly, the equations in section 3.1.3, i.e. Equations (9a)–(9d), under a 70% contraction ratio and different shapes of converging and diverging sections, were embedded in the microcontroller of the digital sensing device. The designed flow-measuring devices mounted with the sensing system were operated under varying flow depths, and the corresponding discharges were estimated and observed. It is observed from Table 4 that the MAD, MSE, RMSE, and MAPE values for the 70% contraction ratio were lower than the contraction ratios of 53.33 and 61%. Moreover, the MAPE value for the straight converging section with a 30° angle and curved- diverging section was the lowest (4.32) and was finally selected as the best-designed section for the measurement of open channel flow.

Selection of the hydraulically efficient and accurate flow measuring device

Based on the occurrence and location of critical depth, as detailed in section 3.1, and the validation results of the developed rating curves, the flow measuring device with a straight converging section (30° angle) and a curved diverging section with a 70% contraction ratio was finalized for subsequent testing under varying flow depths. The dimension of the device is presented in Figure 1(b). The isometric view and prototype of the finalized device are presented in Figure 6. It was observed that the device is a portable one with a total length (L) of 85 cm and a bottom width (B) of 30 cm. The device was developed for use in open field channels in canal command with a width of 30 cm, but its use can be tested on different channel widths by preparing a measuring section of 30 cm width along the main channel using locally available materials. A contraction ratio of 70% provided the stability of flow depth in the approach section for flow measurement, and the critical depth occurred for the entire discharge range of 1.5–20 L s−1. Also, the critical depth occurred at the shortest distance of 6.5 cm compared to devices with 53.33 and 60% contraction ratios. The flow depth vs. discharge rate relationship of the device with a 70% contraction ratio produced parameters ‘α’ and ‘β’ with a skewness of <0.25 (γα and γβ), indicating better accuracy of the developed equation than equations for devices with lower contraction ratios. The maximum depth of mounting based on the field of view of the sensor was estimated to be 56 cm, and the sensor was mounted at a 45 cm depth from the channel bead, and at a 15 cm distance from the starting point of the converging section in the approach section (Kumar et al. 2023).
Figure 6

Selected device for irrigation channel flow measurement: (a) isometric view, (b) prototype device with the flow sensing system.

Figure 6

Selected device for irrigation channel flow measurement: (a) isometric view, (b) prototype device with the flow sensing system.

Close modal

The output of the sensing system displays data and stores data in the data logger pertaining to discharge rate besides the total flow volume in a given time span, which is useful for irrigating specific farmland with a desired depth of water. A developed device with a straight converging section (30° angle) and a curved diverging section having a 70% contraction ratio estimated the discharge rate in close agreement with the observed discharge rate, which was corroborated with the lowest MAPE (4.3) among all tested devices. The finalized design of the new flow measuring device was 85 cm in length and 30 cm in width. The transition point in the contraction section was 54.5 cm from the upstream edge of the device (la + lb). The digital sensing system was mounted over it with a sensor depth of 45 cm from the channel bottom.

The development of a low-cost, portable, and digital open-channel flow measuring device marks a significant advancement in irrigation water management and a critical step toward improving water use efficiency in irrigation systems. However, multi-location trials of the device may be necessary to enhance its practical relevance and broader impact. Field validation under diverse channel shapes, sediment loads, and environmental factors is essential to ensure the versatility of the device for operation under real-world situations. Prioritizing user-friendliness and ease of maintenance is crucial for adoption by farmers with limited technical expertise. Additionally, assessing cost-effectiveness, production costs, scalability, and commercialization for mass deployment will establish its economic viability. Employing advanced data analysis techniques can enhance the interpretation and utility of collected data while ensuring compatibility with smart irrigation systems that can facilitate real-time monitoring and decision-making. Addressing environmental impacts, social and economic implications, and ethical concerns such as data privacy will further promote sustainable and responsible adoption. Therefore, refining these aspects will strengthen the device's contribution to achieving efficient and sustainable water management in agriculture.

Accurate measurement of irrigation water plays a significant role in not only saving water but also enhancing water productivity in agriculture. Indiscriminate use of water coupled with inaccurate scheduling of irrigation in agriculture poses hindrances to the attainment of sustainability in irrigated agriculture. To enhance the water productivity of canal commands and other regions under surface irrigation, there is a need for low-cost flow measuring devices in the field channels of canal commands to supply a measured quantity of water as per crop water requirements to farmers' fields. Keeping this in view, an attempt was made in this study to develop a portable, low-cost, and digital sensing device to measure the discharge rate and the volume of water passing through the irrigation channel in a given period of time. The flow measuring device was a modified flume without a throat section and a contraction (70%) of the central part to create the critical flow. Furthermore, a digital flow depth sensing system using an ultrasonic sensor was integrated with the device, and the head vs. discharge equation was embedded in the microcontroller to measure the discharge rate instantaneously and the total volume of water passing through the device over a given period of time. Nonetheless, the device will be useful for the canal commands of the country not only to enhance the agricultural water productivity but also to bridge the gap between the potential lag, i.e. irrigation potential created and utilized.

The authors acknowledge the financial assistance to undertake research under the sub-project of ‘Hub and Spoke model (DIC: 24-596)’ funded by MHRD, GoI, in collaboration with IIT, Kanpur, and operational at the Division of Agricultural Engineering, ICAR-Indian Agricultural Research Institute, New Delhi.

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

The authors declare there is no conflict.

Abioye
E. A.
,
Hensel
O.
, Hensel, O., Esau, T. J., Elijah, O., Abidin, M. S. Z., Ayobami, A. S., Yerima, O. & Nasirahmadi, A.
(
2022
)
Precision irrigation management using machine learning and digital farming solutions
,
AgriEngineering
,
4
(
1
),
70
103
.
Abt
S. R.
&
Staker
K. J.
(
1990
)
Rating correction for lateral settlement of Parshall flumes
,
Journal of Irrigation and Drainage Engineering
,
116
(
6
),
797
803
.
Abt
S. R.
,
Thompson
K.
&
Staker
K. J.
(
1989
)
Discharge correction for longitudinal settlement of Parshall flumes
,
Transactions of the ASAE
,
32
(
5
),
1541
1544
.
Abt
S. R.
,
Genovez
A.
&
Florentin
B.
(
1994
)
Correction for settlement in submerged Parshall flumes
,
Journal of Irrigation and Drainage Engineering
,
120
(
3
),
676
682
.
Abt
S. R.
,
Florentin
C. B.
,
Genovez
A.
&
Ruth
B. C.
(
1995
)
Settlement and submergence adjustments for Parshall flume
,
Journal of Irrigation and Drainage Engineering
,
121
(
5
),
317
321
.
Abt
S. R.
,
Ruth
B. C.
&
Skowron
E. M.
(
1998
)
Rating adjustment for settlement of cutthroat flumes
,
Journal of Irrigation and Drainage Engineering
,
124
(
6
),
311
315
.
Bos
M. G.
,
Replogle
J. A.
&
Clemmens
A. J.
(
1984
)
Flow Measuring Flumes for Open Channel Systems
.
New York, NY, USA
:
John Wiley and Sons
, p.
321
.
Chow
V. T.
(
1988
)
Open-channel Hydraulics. Classical Textbook Reissue
.
New York, NY, USA: McGraw.
Cone
V. M.
(
1917
)
The Venturi Flume
.
Journal of Agriciltural Research, 9 (4), 115–129. https://www.biodiversitylibrary.org/item/261007#page/1/mode/1up
Genovez
A.
,
Abt
S.
,
Abt
B.
&
Garton
A.
(
1993
)
Correction for settlement in a Parshall flume
,
Journal of Irrigation and Drainage Engineering
,
119
(
6
),
1081
1091
.
Goel
A.
,
Verma
D. V. S.
&
Sangwan
S.
(
2015
)
Open channel flow measurement of water by using width contraction
,
World Academy of Science, Engineering and Technology
,
9
(
2
),
1557
1562
.
Hager
W. H.
(
1988
)
Mobile flume for circular channel
,
Journal of Irrigation and Drainage Engineering
,
114
(
3
),
520
534
.
Hougaard
P.
(
1985
)
The appropriateness of the asymptotic distribution in a nonlinear regression model in relation to curvature
,
Journal of the Royal Statistical Society: Series B (Methodological)
,
47
,
103
114
.
Hyatt
,
M. L.
(
1965
)
Design, Calibration, and Evaluation of A Trapezoidal Measuring Flume by Model Study
,
Utah Water Research Laboratory, Utah State University, Logan, UT, USA. https://digitalcommons.usu.edu/water_rep/386
Kruse
E. G.
(
1992
)
Discussion of rating correction for lateral settlement of Parshall flumes
,
Journal of Irrigation and Drainage Engineering
,
118
(
2
),
337
339
.
Kumar
V.
(
2020
)
A study of travel time for different open channels
,
Journal of The Institution of Engineers (India): Series A
,
101
(
2
),
399
407
.
Kumar, A. & Sarangi, A. (2022) Flow measuring devices in surface irrigation for enhancing agricultural water productivity. Current Science, 122 (10), 1135–1144.
Kumar
A.
,
Sarangi
A.
,
Singh
D. K.
,
Dash
S.
&
Mani
I.
(
2023
)
Evaluation of ultrasonic sensor for flow measurement in open channel
,
Journal of Scientific & Industrial Research
,
82
,
1091
1099
.
Parsaie
A.
,
Najafian
S.
,
Omid
M. H.
&
Yonesi
H.
(
2017
)
Stage discharge prediction in heterogeneous compound open channel roughness
,
ISH Journal of Hydraulic Engineering
,
23
(
1
),
49
56
.
Parshall
R. L.
(
1928
)
Improved Venturi Flume. Agricultural Experiment Station, Colorado Agricultural College, Bull. No. 336
. https://archives.colostate.edu/repositories/2/archival_objects/82673.
Prasad
S.
(
1991
)
Performance Evaluation of A Parshall Flume Under Different Simulated Settlements
.
Master's thesis
,
Division of Agricultural Engineering, ICAR-IARI
,
New Delhi-110012, India
.
Ratkowsky
D.
(
1990
)
Handbook of Nonlinear Regression Models
.
New York, NY, USA
:
Marcel Dekker
.
Replogle
J. A.
(
1975
). '
Critical flow flumes with complex cross section
',
Proc., Irrigation and Drainage Division Special Conf. ASCE
.
Reston, VA
, USA.
Robinson
A. R.
&
Chamberlain
A. R.
(
1960
)
Trapezoidal Flumes for Open-Channel Flow Measurement
. Transactions of the ASAE, 3 (2), 0120–0124. doi:10.13031/2013.41138
Samani
Z.
(
2017
)
Three simple flumes for flow measurement in open channels
,
Journal of Irrigation and Drainage Engineering
,
143
(
6
),
04017010
.
Samani
Z.
&
Magallnez
H.
(
1993
)
Measuring water in trapezoidal canals
,
Journal of Irrigation and Drainage Engineering
,
119
(
1
),
181
186
.
Samani
Z.
&
Magallnez
H.
(
2000
)
Simple flume for flow measurement in open channel
,
Journal of Irrigation and Drainage Engineering
,
126
(
2
),
127
129
.
Samani
Z.
,
Jorat
S.
&
Yousaf
M.
(
1991
)
Hydraulic characteristics of circular flume
,
Journal of Irrigation and Drainage Engineering
,
117
(
4
),
558
566
.
Santhosh
K. V.
&
Roy
B. K.
(
2012
)
An intelligent flow measurement technique using ultrasonic flow meter with optimized neural network
,
International Journal of Control and Automation
,
5
,
185
196
.
Skogerboe
G. V.
(
1969
)
Cut Throat Flow Measuring Flume for Flat Gradient Channels. Seventh Congress on Irrigation and Drainage. Mexico City ICID. Report 16
.
Utah Water Research Laboratory, Utah State University, Logan, UT, USA. https://digitalcommons.usu.edu/water_rep/96.
Skogerboe
G. V.
&
Hyatt
M. L.
(
1967
)
Rectangular Cutthroat Flow Measuring Flumes
.
Utah Water Research Laboratory, Utah State University, Logan, UT, USA. https://digitalcommons.usu.edu/water_rep/87.
Skogerboe
G. V.
,
Hyatt
M. L.
&
Eggleston
K. O.
(
1967a
)
Design and Calibration of Submerged Open Channel Flow Measurement Structures: Part 1-Submerged Flow
.
Utah Water Research Laboratory, Utah State University, Logan, UT, USA. https://digitalcommons.usu.edu/water_rep/93
Skogerboe
G. V.
,
Hyatt
M. L.
,
England
J. D.
&
Johnson
J. L.
(
1967b
)
Design and Calibration of Submerged Open Channel Flow Measurement Structures. Part 2: Parshall Flumes
.
Rep. WG31-3, Utah Water Research Laboratory, Utah State University, Logan, UT, USA. https://digitalcommons.usu.edu/water_rep/81.
Skogerboe
G. V.
,
Bennet
R. S.
&
Walker
W. R.
(
1972
)
Generalized discharge relations for cutthroat flumes
,
Journal of Irrigation and Drainage Engineering
,
98
(
4
),
569
583
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).