Study on the Characteristic Point Location of Depth Average Velocity in Smooth Open Channels: Applied to Channels with Flat or Concave Boundaries

Based on the flow partition theory, we derive a mathematical expression by using the log-law for the characteristic point location (CPL) of depth average velocity in channels with flat or concave boundaries. It can manifest the position of the characteristic points in the vertical direction relative to the channel side wall or bed. Taking rectangular and semi-circular channels as research objects, we put forward a method to calculate the discharge of channels with CPL. Additionally, we carried out some experiments on rectangular and semi-circular channel sections. CPL’s analytic expression is validated against experimental results through comparison of velocity and discharge. The proposed formulation of characteristic point location could be extensively employed in flow measurements of flat and concave boundary channels, which has practical application value in simplifying the flow measurement steps of open channels.


Introduction
The contradiction between the serious shortage of agricultural water resources and unreasonable utilization has been restricting the development of agricultural economy [1]. Precision agriculture provides the potential to enhance irrigation efficiency through specific methods of flow measurement in channels.
Flow measurement methods through weir-gate structures [2] and velocity-area methods are widely used in flow measurement. The principle of a sluice gate or weir is based on the concept of critical flow, making it possible to measure only the depth and calculate the flow rate, thus simplifying the continuous monitoring of flow rate [3]. However, critical flow in open channels can be formed through two general methods: raising the bottom of the channel or contracting the cross-sectional area of the flow [3,4]. Neither of these two methods is easy to implement. Impurities on the sluice gate or weir accumulate easily and need to be cleaned regularly, which is difficult to achieve during actual production. Apart from measuring with weir-gate structures, the velocity-area method for calculating discharge in an open channel is considered to be particularly reliable. However, a large number of measurements results in an associated time cost of measurement. Numerous researchers have proposed different forms of the velocity distribution law in order to predict flow discharge accurately.
The study of velocity distribution is critical to the calculation of the discharge. Relevant research can be traced back to the 1930s. In 1938, Keulegen [5] first proposed that the log-law could be used to describe the time-averaged velocity distribution of fully developed open channel uniform turbulence: where u is the velocity in the channel, u * is frictional velocity, ν is the coefficient of kinematic viscosity, κ is the von Kármán constant, and ω is constant. Nikuradse made an artificial sand rough pipe using pasting sand with a uniform particle size on the wall of a pipe [6,7]. The analysis of the experimental data of the smooth sand tube showed that the values of parameters κ and ω in log-law are 0.4 and 5.5. Different researchers have given different values of the von Kármán constant κ and constant ω according to their experimental data [8][9][10]. This paper uses the parameter values given by Nikuradse for further derivation.
Under the guide of velocity distribution, researchers have used different methods to measure velocity. Chiu [11][12][13] developed the entropy model, correlating the mean flow velocity to the maximum value using the linear relationship Um = ϕ(M)U max , which depends on the entropic parameter M. Maghrebi [14] used the similarity between the magnetic field of a current wire and the isovel contours in a channel cross-section to derive the isovel patterns. Then, one can easily obtain the discharge using a single point on the isovel patterns of velocity measurement. However, the measured points are only selected from the upper half of the water depth and away from the boundaries, which would be much closer to the measured discharge. Bretheim et al. [15] used the restricted nonlinear (RNL) model to simulate wall turbulence and obtained the real mean velocity distribution in the channel at a low Reynolds number. Hong et al. [16] employed a systematic measuring technology combining ground-penetrating radar and surface velocity radar and established the rating curves representing the relation of water surface velocity to the channel cross-sectional mean velocity and flow area. Then, stream discharge was deduced from the resulting mean velocity and flow area. Moramarco et al. [17] proposed a new method for estimating discharge from surface velocity monitoring (u surf ). Based on entropy theory and sampling u surf , it can identify the two-dimensional velocity distribution in the cross-sectional flow area. This method is more accurate for rivers with a lower aspect ratio where secondary currents are expected. Johnson and Cowen [18,19] predicted the mean streamwise velocity and the depth-averaged velocity by permitting remote determination of the velocity power-law exponent. Then, the volumetric discharge from surface measurements of currents can be determined. Khuntia et al. [20] presented a new methodology to predict the depth-averaged velocity. They used multi-variable regression analysis to develop five models to predict the point velocities in terms of non-dimensional geometric and flow parameters at any desired location. Through these attempts, researchers have discussed new thoughts of quantifying channel discharge. But few people determine the average velocity of a channel by looking for the location of characteristic points.
Overall, this study investigates a new concept of characteristic point location aimed at estimating the discharge in open channels, which is based on flow partitioning theory and the log-law. This method is applied to flat channels (e.g., a rectangular channel) and concave channels (e.g., a semi-circle channel). The derived theoretical expression of characteristic point location is also verified by the experimental data. We expect that the method for flow measurement conveniently proposed by this paper could be used widely for different varieties of irrigation channels.

Methodology
This paper assumes that characteristic points exist on the channel cross-section whose velocities could be called characteristic velocities. We can measure these points to get the precise measurements of velocity and discharge. The positions of these points are called characteristic point locations of depth average velocity (CPL).

Existing Form of the Division Line in Rectangular Channels
Consider the case of steady, uniform flow in a smooth rectangular channel with an aspect ratio of b/h, where h = water depth and b = channel width. Due to existence of a secondary current, the division line is actually a zero Reynolds shear stress line in the flow area, and also a line without energy cross-transfer [21]. Yang et al. [22] have proven the physical existence of division lines by using experimental data recorded at the Hydraulics Laboratory, University of Wollongong, as well as the experimental data recorded previously by Melling and Whitelaw [23] and that of Tracy [24].
Einstein [25] proposed that the flow cross-sectional area can be divided into three sub-sections that correspond to the side wall and channel bed. However, he did not propose any method of determining the exact location of the division line. Chien and Wang [26] performed an in-depth study of Einstein's idea. They gave the total overall boundary shear stress, τ 0 , as: where ρ = fluid density; g = gravitational acceleration; R = hydraulic radius; and S e = energy slope. Then, they divided τ 0 into two parts corresponding to the bed, τ b , and side wall, τ w . They also defined two hydraulic radii, one for the wall R w , and one for the bed R b . Thus: where τ b = mean bed shear stress;τ w = mean wall shear stress; A b = flow area corresponding to channel bed; A w = flow area corresponding to channel side wall; R w = hydraulic radius corresponding to side wall; and R b = hydraulic radius corresponding to bed. According to Yang [27], the mechanical energy contained in any flow will be transmitted to the boundary nearest to the "relative distance" and dissipated. The dimensionless relative distance Φ is defined as the ratio of geometrical length between a point in the flow field and the characteristic boundary length, representing the energy dissipation capacity of the boundary (for a smooth boundary channel, this would be the viscous layer thickness D s ).
The side wall and bed will each have their own fair share of surplus energy transferred from the main flow depending on the minimum relative distance Φ. Accordingly, for any unit volume in the flow field, there are two possible ways to transfer the surplus energy: either toward the bed (Φ b ) or toward the side wall (Φ w ). If Φ w ≤ Φ b , the energy will be transferred toward the side wall. Conversely, if Φ w ≥ Φ b , then the energy will be transported to the bed. It follows that the condition where Φ w = Φ b will define a division line by which the flow region near the corner is divided into two sub-flow sections.
Therefore, the division line Figure 1a can be expressed as: where z and y denote the horizontal and vertical axes in the Cartesian coordinates system, respectively; y = vertical direction normal to bed; z = transverse flow direction;D sw = characteristic length for smooth side wall; and D sb = characteristic length for smooth channel bed. Substituting D sw = C s ν/u * w and D sb = C s ν/u * b into Equation (5) gives: where C s = constant; u * w = mean side wall shear velocity; and u * b = mean bed shear velocity. The geometric significance of k in Figure 1a is the reciprocal of the slope of the line OD. Point D is the intersection of the division line and the water surface line.
where s C = constant; *w u = mean side wall shear velocity; and *b u = mean bed shear velocity.
The geometric significance of k in Figure 1a is the reciprocal of the slope of the line OD. Point D is the intersection of the division line and the water surface line. Then, it can be obtained that: Using Equations (3)-(8), k can be evaluated, yielding: Or: Using a similar method, Yang and Lim [27] also gave the evaluation of k1 in Figure 1b: Then, it can be obtained that: Using Equations (3)-(8), k can be evaluated, yielding: Or: Using a similar method, Yang and Lim [27] also gave the evaluation of k 1 in Figure 1b: where k 1 represents the slope of the division line OD in Figure 1b. Highly inspired by the derivation of Yang and Lim [27], we also take the two cases of Figure 1 to analyze. It has been proven that the critical aspect ratio α is equal to 2 in smooth rectangular channels by Yang and Lim [27]. As the rectangle is axisymmetric, we take the vertical line of the channel section as the axis and only analyze half of it. The division line can be determined by Equation (10) when b/h ≥ 2 ( Figure 1a). Since 2h/b is a constant, solving Equation (10) will give the value of k.
Referring to the value of k, a comparison relationship between 1/k and 2h/b is as follows: Water 2020, 12, 430

of 17
It can be shown from Inequality (13), so That is to say, when b/h ≥ 2, the slope of the division line is always greater than 2h/b. The intersection point of the division lines is on or above the water surface line.
The division line can be determined by Equation (11) when b/h ≤ 2 ( Figure 1b). Since b/2h is a constant, the value of k 1 can be obtained by solving Equation (11).
Referring to the value of k 1 , a comparison relationship between k 1 and 2h/b is as follows: It can be shown from Inequality (15), so That is to say, when b/h ≤ 2, the slope of the division line is always greater than 2h/b. The intersection point of the division lines is on or above the water surface line, the same as b/h ≥ 2.
By combining cases (a) and (b), we know that Figure 1b will not occur. Whether in a wide-shallow or in a narrow-deep channel, the existing form of the division line is always shown as Figure 1a.

CPL of Lines in Rectangular Channel
As shown in Figure  In region I, the discharge of the shadow rectangle area can be expressed using Equation (16): In region I, the discharge of the shadow rectangle area can be expressed using Equation (16): where dQ is the discharge of the shadow rectangle area; dA is the area of shadow rectangle; and u is the longitudinal velocity of any point in the area. By integrating y, the discharge Q can be obtained as follows: As the rectangular area can be expressed as A = h I dz, the average velocity of this flow area can be obtained as: According to the log-law of vertical velocity in channel section: Substituting Equation (19) to (18) gives: Combining Equations (19) and (20) gives: y is expressed as: where y I represents the location of depth average velocity in region I, and h I (= h) is the water depth of the channel. In region II, y II represents the CPL of depth average velocity, which is similar to region I: where y II represents the location of depth average velocity in region II, and h II is the distance perpendicular to the channel bed to the division line. In region III, as shown in Figure 2, the discharge of the shadow rectangle area can be expressed as: Integrating z along the direction of z axis, the discharge can be determined using Equation (25): where h III is the horizontal distance from the side wall of the channel to the division line. As the rectangular cross-sectional area can be expressed as A = h III dy, the average velocity of this flow area can be obtained as shown as Equation (26): According to the log-law of vertical velocity in channel section: Using Equations (26) and (27), z can be evaluated using: Giving: where z III represents the location of the depth average velocity in region III.

CPL of Regions in Rectangular Channel
For region I, according to Equation (22), since the water depth h in the channel is constant, measuring the velocity of the CPL of depth average velocity can provide the average velocity for the whole section of region I. Therefore, the average discharge can be obtained by multiplying the average cross-section velocity by the area of the region. Similarly, the flow in regions II and III can be obtained using the same method. Figure 3 shows that P I , P II , and P III are the characteristic points that represent the mean velocities in regions I, II, and III, respectively.  For region II, the discharge Q can be expressed using Equation (30) (10) where e is a constant; S is the slope of channel; and z0 is shown in Figure 3.
The mean velocity can be calculated from the value of Q/A, which can be expressed using For region II, the discharge Q can be expressed using Equation (30): where e is a constant; S is the slope of channel; and z 0 is shown in Figure 3. The mean velocity can be calculated from the value of Q/A, which can be expressed using Equation (31): Using Equations (6), (19) and (31), Equation (32) can be obtained: Therefore, the y and z value of the point P II can be calculated using Equation (33): (33) where z II and y II represent the coordinates of P II .
For region III, the analysis method is similar to that of region II. The z and y value of the point P III , which can be used to calculate the value of Q of region III, can be expressed using Equation (34): (34) where z III and y III represent the coordinates of P III .

CPL in Semi-Circular Channel
The concave boundary cross-section is one of the most commonly used cross-section forms in farmland canal systems, open-flow culverts, urban drainage pipelines and other projects. Compared with the rectangular section, it has a larger hydraulic radius under the same conditions, so the velocity of the cross section is larger and the water flow that can be transported is also larger [28]. In order to obtain the discharge of the concave boundary cross-section channel more quickly and accurately, it is necessary to master the velocity distribution of the cross-section. Figure 4 shows the semi-circle channel cross-section.
where III z and III y represent the coordinates of PIII.

CPL in Semi-Circular Channel
The concave boundary cross-section is one of the most commonly used cross-section forms in farmland canal systems, open-flow culverts, urban drainage pipelines and other projects. Compared with the rectangular section, it has a larger hydraulic radius under the same conditions, so the velocity of the cross section is larger and the water flow that can be transported is also larger [28]. In order to obtain the discharge of the concave boundary cross-section channel more quickly and accurately, it is necessary to master the velocity distribution of the cross-section. Figure 4 shows the semi-circle channel cross-section. As shown in Figure 4, a straight line, which crosses the center of the circle and is angled with the center line (θ), is taken from a semicircular channel with radius R and maximum water depth h.
The coordinate system is established by taking the straight line and the direction perpendicular to the straight line as the coordinate axis. According to Figure 4, the rectangle FGJK is Ln in length and dr in width, where Ln is the length of the underwater part of the radius. With the rectangle FGJK, As shown in Figure 4, a straight line, which crosses the center of the circle and is angled with the center line (θ), is taken from a semicircular channel with radius R and maximum water depth h. Water 2020, 12, 430 9 of 17 The coordinate system is established by taking the straight line and the direction perpendicular to the straight line as the coordinate axis. According to Figure 4, the rectangle FGJK is L n in length and dr in width, where Ln is the length of the underwater part of the radius. With the rectangle FGJK, the area is A FGJK = L n · dr and the discharge through FGJK can be expressed as: Then, taking a shadow rectangle whose edge lengths are dρ and dr (dρ and dr are infinitesimal), the discharge through the shadow rectangle is: By integrating dρ in the direction of y' and substituting the log-law, we get: Substituting the log-law into Equation (35) gives: For Equations (37) and (38), the coordinates of the characteristic points representing the average velocity of the straight line can be obtained from Equation (39): where ρ is the CPL of the depth average velocity of the line.

Discharge in a Rectangular Channel
According to Section 2.1.3, based on the CPLs of different regions, the discharge from the cross-section can be given by Equation (40): where A I , A II and A III represent the areas of three regions, and u P I , u P II and u P III represent the velocities of points P I , P II , and P III .

Discharge in a Semi-Circular Channel
Similar to a rectangular channel, in a semi-circle channel, after obtaining the CPL of every normal line, it is necessary to know how to calculate the discharge through simplified point distribution. As shown in Figure 5, two straight lines passing through the center of the circle have divided the angle into θ1, θ2 and θ3. The discharge can be expressed as Equation (41): where A 1 , A 2 and A 3 represent sub-areas, u CPL is the velocity of the normal line and α c is the velocity coefficient along the bank. α c can be determined according to the test conditions.
where A1, A2 and A3 represent sub-areas, CPL u is the velocity of the normal line and αc is the velocity coefficient along the bank. αc can be determined according to the test conditions.     [29]. Different groups were set up by changing the discharge (or water depth) of the channel. Using acoustic doppler velocimetry (ADV) and other velocity measurement tools to measure the vertical velocity of a cross-section, a large amount of velocity data can be obtained.

Analysis with Rectangular Channels
The consistency of CPL of R1, R2 and R3 in rectangular channels is analyzed in Figure 6. These scatter points cover three regions of the rectangular cross-section. The abscissa of the scatter points in the figure represents the theoretical value of CPL and is calculated by Equations (22), (23) and (29). The longitudinal coordinates of points represent the measured value of CPL, which can be obtained by interpolating the experimental data. It can be seen that the data points are mostly distributed within the error line of ±15%. Therefore, it can be proven that the formulas of CPL in Section 2.1.2 can be applied to rectangular cross-section channels. scatter points cover three regions of the rectangular cross-section. The abscissa of the scatter points in the figure represents the theoretical value of CPL and is calculated by Equations (22), (23) and (29). The longitudinal coordinates of points represent the measured value of CPL, which can be obtained by interpolating the experimental data. It can be seen that the data points are mostly distributed within the error line of ±15%. Therefore, it can be proven that the formulas of CPL in Section 2.1.2 can be applied to rectangular cross-section channels. In order to further verify the correctness of the theoretical derivation, a series of experiments were conducted at different water depths at the CAU (China). We collated the experimental data of R4-R7 and selected the representative data points to draw in Figure 7. Since the perpendicular In order to further verify the correctness of the theoretical derivation, a series of experiments were conducted at different water depths at the CAU (China). We collated the experimental data of R4-R7 and selected the representative data points to draw in Figure 7. Since the perpendicular bisector is in region I, we calculated the point of CPL with Equation (22). Since the number of measuring points was less, and there were some test errors, the relative error between CPLt and CPLm still existed. With the increase in water depth, the influence of several measuring points on the results is apparent (R7). The velocities of CPLt and CPLm have a greater consistency when the number of measuring points is high. Hence, we have reason to think that the velocity of CPLt can represent the average velocity of the centerline in rectangular channels.
Water 2020, 12, x FOR PEER REVIEW 12 of 18 bisector is in region I, we calculated the point of CPL with Equation (22). Since the number of measuring points was less, and there were some test errors, the relative error between CPLt and CPLm still existed. With the increase in water depth, the influence of several measuring points on the results is apparent (R7). The velocities of CPLt and CPLm have a greater consistency when the number of measuring points is high. Hence, we have reason to think that the velocity of CPLt can represent the average velocity of the centerline in rectangular channels. The characteristic points of depth average velocity in each region are plotted in Figure 8 and recorded as PI, PII and PIII, which are calculated by Equations (22), (33) and (34), respectively. We can obtain the discharge by using the values of PI, PII and PIII in Equation (40). Curves in the graph are isovel patterns in the rectangular cross-section channel based on UOW test data, which is consistent with the contour line measured by Chiu and Chiou [30]. It can be seen from the figures that the velocity near the wall is small and the velocity far from the wall is large. This can be explained The characteristic points of depth average velocity in each region are plotted in Figure 8 and recorded as PI, PII and PIII, which are calculated by Equations (22), (33) and (34), respectively. We can obtain the discharge by using the values of PI, PII and PIII in Equation (40). Curves in the graph are isovel patterns in the rectangular cross-section channel based on UOW test data, which is consistent with the contour line measured by Chiu and Chiou [30]. It can be seen from the figures that the velocity near the wall is small and the velocity far from the wall is large. This can be explained by the log-law of velocity distribution. Increasing the monotone leads to the velocity increasing as the distance from the channel side wall or bed increases. Therefore, the maximum velocity point of the cross-section should theoretically be located at the intersection of the vertical line of the channel and the water surface.
Water 2020, 12, x FOR PEER REVIEW 13 of 18 to shift the maximum velocity from the free surface to the bed, its identification is still a complex task in hydraulics [34].

Analysis with Semi-Circular Channels
Mechanical energy is always transmitted to the closest boundary. As shown in Figure 9, this refers to the normal line direction. So, in a concave boundary channel, the analysis is carried out along the normal direction perpendicular to the tangent direction of the wall. In order to compare the accuracy of velocity distribution analysis along the vertical direction and along the normal wall direction, the following analysis has been made. The error values obtained by the experimental condition S2 (in Table 1) and log-law formulas are compared in tabular analysis in Table 2. The calculation formula of the average error value is given by: where E is the average error value in Table 2, Em denotes the measured value, and Ec denotes the calculated value. However, in Figure 8, we can see that the maximum velocity point is below the water surface, not on it. This phenomenon was first observed about a century ago [31,32], and further experimental studies showed that it was induced by the presence of secondary cross-sectional flow structures [33]. Due to the occurrence of anisotropic turbulence and cross-sectional secondary currents, which tend to shift the maximum velocity from the free surface to the bed, its identification is still a complex task in hydraulics [34].

Analysis with Semi-Circular Channels
Mechanical energy is always transmitted to the closest boundary. As shown in Figure 9, this refers to the normal line direction. So, in a concave boundary channel, the analysis is carried out along the normal direction perpendicular to the tangent direction of the wall. to shift the maximum velocity from the free surface to the bed, its identification is still a complex task in hydraulics [34].

Analysis with Semi-Circular Channels
Mechanical energy is always transmitted to the closest boundary. As shown in Figure 9, this refers to the normal line direction. So, in a concave boundary channel, the analysis is carried out along the normal direction perpendicular to the tangent direction of the wall. In order to compare the accuracy of velocity distribution analysis along the vertical direction and along the normal wall direction, the following analysis has been made. The error values obtained by the experimental condition S2 (in Table 1) and log-law formulas are compared in tabular analysis in Table 2. The calculation formula of the average error value is given by: Figure 9. Semi-circular channel cross-section and notation.
In order to compare the accuracy of velocity distribution analysis along the vertical direction and along the normal wall direction, the following analysis has been made. The error values obtained by the experimental condition S2 (in Table 1) and log-law formulas are compared in tabular analysis in Table 2. The calculation formula of the average error value is given by: where E is the average error value in Table 2, E m denotes the measured value, and E c denotes the calculated value. From Table 2, we can conclude that the average error value obtained by normal analysis along the wall is much lower than that obtained by vertical analysis when consistency analysis of the velocity measurements and the log-law is carried out. Hence, the correctness of the analysis of the concave boundary channel along the normal direction perpendicular to the wall is proved. Figure 10 shows the comparison of average velocity along the normal line under S1−S5. The horizontal ordinates in the figure are the average value of the velocity at each normal measuring line. The vertical coordinates are obtained by theoretical calculation. Specifically, the CPL is calculated by the formula from Section 2.2, and then the average velocity on the normal line is obtained by substituting it into the log-law.  From Table 2, we can conclude that the average error value obtained by normal analysis along the wall is much lower than that obtained by vertical analysis when consistency analysis of the velocity measurements and the log-law is carried out. Hence, the correctness of the analysis of the concave boundary channel along the normal direction perpendicular to the wall is proved. Figure 10 shows the comparison of average velocity along the normal line under S1−S5. The horizontal ordinates in the figure are the average value of the velocity at each normal measuring line. The vertical coordinates are obtained by theoretical calculation. Specifically, the CPL is calculated by the formula from Section 2.2, and then the average velocity on the normal line is obtained by substituting it into the log-law. Despite some fluctuations of a few test points in Figure 10, most of them are within the error range of 5%. The statistical rules of a large number of test points can still prove the accuracy of the theoretical deduction in Section 2.2. Therefore, the velocity corresponding to the CPL calculated by Equation (39) can represent the average velocity of the channel wall along the normal line to the water surface line. Most points, especially the measuring points close to the water surface, are located above the agreement line, which implies that Ut is bigger than Um. This is because the log-law appears to deviate near the water surface. Based on this, the log-wake law was proposed by Coles [35], which Despite some fluctuations of a few test points in Figure 10, most of them are within the error range of 5%. The statistical rules of a large number of test points can still prove the accuracy of the theoretical deduction in Section 2.2. Therefore, the velocity corresponding to the CPL calculated by Equation (39) can represent the average velocity of the channel wall along the normal line to the water surface line. Most points, especially the measuring points close to the water surface, are located above the agreement line, which implies that U t is bigger than U m . This is because the log-law appears to deviate near the water surface. Based on this, the log-wake law was proposed by Coles [35], which appears to be the most reasonable extension of the log-law. However, the value of Π in the log-wake law seems not to be universal [36]. So, more research is needed to correct the velocity near the water surface.
As shown in Figure 11, the theoretical and measured values of CPL on each normal line are very close, located in the middle and lower part of the normal line. The solid and dotted lines are smooth arc curves, which are consistent with the concave boundary. Therefore, the accuracy of the formula in Section 2.2 can be proved. It is advisable to use the velocity of CPL to represent the average velocity of the normal line.
Water 2020, 12, x FOR PEER REVIEW 15 of 18 As shown in Figure 11, the theoretical and measured values of CPL on each normal line are very close, located in the middle and lower part of the normal line. The solid and dotted lines are smooth arc curves, which are consistent with the concave boundary. Therefore, the accuracy of the formula in Section 2.2 can be proved. It is advisable to use the velocity of CPL to represent the average velocity of the normal line. When the S1 condition is applied to Equation (41), we can calculate the discharge of the semicircular channel. Here, αc = 0.8. Table 3 shows the calculated discharge of the channel with the combination of six kinds of angles. Only the first is a uniform partition; the others are all non-uniform. Relative error refers to the error value of the calculated discharge relative to the measured discharge. When the S1 condition is applied to Equation (41), we can calculate the discharge of the semi-circular channel. Here, α c = 0.8. Table 3 shows the calculated discharge of the channel with the combination of six kinds of angles. Only the first is a uniform partition; the others are all non-uniform. Relative error refers to the error value of the calculated discharge relative to the measured discharge. It can be found that the relative error value is the smallest only when the angle is divided uniformly, about 5%. Therefore, we can reasonably conclude that the accuracy of the discharge calculation is the highest only if the angle is divided uniformly while using the two-line method. Table 3. Analysis of calculated and measured discharge under different combinations of angles (S1).

Summary
This study has investigated a fast flow measurement method based on the distribution of measuring points by conducting laboratory experiments. The principal achievements from the comparison of the theoretical derivation and experimental study are summarized below.
(1) Based on Yang et al.' s partitioning theory [21,22,27], this paper gives a re-description of the existing form of the division line of a rectangular cross-section channel. That is, whether the channel cross-section is wide-shallow or narrow-deep with the center line of the cross-section as the symmetrical axis, and whether the intersection points of the left and right division lines intersect on or above the water surface. (2) This paper analyzes characteristic points in flat channels (e.g., rectangular channel) and concave boundary channels (e.g., semi-circular channel). In the rectangular channel, the division line divides the section into three regions. In each region, the analysis is conducted in the direction perpendicular to the bottom or side wall of the channel. In the semi-circular channel, the analysis is conducted along the normal direction. Based on the log-law, the theoretical expressions for calculating the location of the average velocity characteristic points in flat and concave boundary channels are derived through the formula transformation. (3) The velocity data in different experimental sites are used to verify the validity of the CPL formulas applied to flat and concave boundary channels. Moreover, the discharge calculation formulas of channels are given through discussion with CPL.

Conclusions
The determination of CPL has provided a theoretical basis for simplifying and accurately measuring channel flow. It is precisely because of the realization of CPL that the average cross-section discharge can be calculated quickly so that the accuracy and stability of flow measurement in irrigation channels can be further improved.
The results of this paper can only serve for smooth open channels. However, if the roughness is non-uniform, the location of the division lines in the channel cross-section will be different [25]. It is also difficult to ensure that the velocity distribution in each sub-region is the same. With the comprehensive ecological treatment of channels becoming a hot topic, the transition from smooth channels to vegetation covered channels becomes relevant. So, it is of necessity to conduct experiments in open channels with roughness for future research.