Range of variation for the main parameters of the laboratory experiments.
The assessment of water discharge in open channel flow is one of the most crucial issues for hydraulic engineers in the fields of water resources management, river dynamics, eco-hydraulics, irrigation, hydraulic structure design, etc. Recent studies state that the entropy velocity law allows expeditive methodology for discharge estimation and rating curve development due to the simple mathematical formulation and implementation. A lot of works have been developed based on the entropy velocity profile supporting measurements in lab for rating curve assessment in regular ditch flows showing a good performance. The present work deals with the use of entropy velocity profile approach in order to give a general framework of threats and opportunities related to robust operational application of such laws in the field of rating curve assessment. The analysis has been carried on a laboratory flume with regular roughness under controlled boundary conditions and different stages generating an exhaustive dashboard for the better appraisal of the approaches. Finally, entropy model may represent a robust and useful tool for the water discharge assessment in rough ditches.
- entropy velocity ratio
- relative submergence
- aspect ratio
- water discharge
Water discharge assessment in open channel still represents a fundamental aspect for hydraulic engineer in several operative and technical fields like water resources management, ecological flow assessment and control, drainage and irrigation system as well as runoff and flood routing model calibration and implementation. Nevertheless, the water discharge evaluation in generic open channel is heavily affected by local fluid dynamics and geometric conditions, which well arise once flow velocity measurements and morphological boundaries are available at the same site. On the other hand, the drainage and irrigation channel present a regular cross section which might provide facilities in water discharge assessment and control, inducing also reduction in time and operative costs. That is, the implementation of operative procedures enabling operative charges simplifying the commitment of field activities, indeed, plays a fundamental role in channel monitoring for natural flow and manmade hydraulic structures. The main idea is related to the definition of expeditive procedures for flow field assessment and water discharge evaluation capable to optimize the surveying resources in time and efforts. Thus, the opportunity to manage with a simple and straightforward velocity law, different from the classical logarithm formulation but capable to provide suitable results is all the more technically fruitful. That is, an operative tool for expeditive velocity distribution assessment basing on simple and immediate parameters.
Recent theoretical and experimental studies endorse the informational content hold into the distributed velocity measurements following an entropy-probabilistic approach. That is, Chiu [1, 2] drew the correlation between the mean flow velocity and maximum flow velocity defining the entropy parameter, M introducing the velocity ratio Φ(M). Considering the important implication that this finding could have for monitoring of high flows in rivers, many authors investigated the reliability of this relationship using field data [3, 4, 5, 6, 7]. Overall, they found M as a river site depending and not influenced by the flood intensity both in terms of amount and duration. Thus, M should be considered a specific factor of the gauged cross section as outlined by Moramarco and Singh  exploring the dependence of M on the hydraulics and geometries of the river cross sections.
The study was able to explain the constancy of M value on the ground that M is not depending on the dynamic of flood, such as expressed by the energy or water surface slope, Sf and to identify a formula expressing M as a function of the hydraulic radius, Manning’s roughness and the location, y0, where the horizontal velocity is hypothetically equal to zero. For the latter, it was preliminarily found that if y0 was assessed by distinguishing low flows from high flows, then a better estimation of M would have been obtained across a gauged river site. However, considering that the y0 location is not of simple assessment and then might have high uncertainty, the assessment of M should be addressed, mainly for ungauged river sites, using hydraulic and geometric variables easy to acquire. Such a thought might be discussed introducing the relative submergence D/d (in which, D = average water depth and d = roughness dimension). That is, the velocity distribution in natural rivers depends on several variables like channel geometry, bed and bank roughness, and the vertical velocity distribution generally increases monotonically from 0 at the channel bed, to the maximum at the water surface and can be assumed 1-D flow dominant. Moreover, whenever the channel cannot be considered “wide”, that is the aspect ratio (B/D with B channel width and D water depth) is less than 6, besides the presence of the boundary, the velocity varies even transversely and a two-dimension distribution occurs, leading G as the 2D entropy parameter. The maximum velocity places below the water surface inducing dip-phenomenon and the position of maximum velocity is also influenced by the aspect ratio , which is of simple assessment once channel cross-section geometry is known. Thus, investigating the influence of bed roughness and cross section geometry on medium and maximum velocity ratio at the global scale assumes a relevant interest in the field of open channel flow.
Therefore, M might represent an intrinsic parameter of the gauged site and this insight led several authors to explore the dependence of M on hydraulic and geometric characteristics of the flow site [3, 7]. In the case of river flows, Greco  enlightened a different behavior of Φ(M) depending on the roughness dimension: the velocity ratio is heavily influenced by the magnitude of relative submergence if large or intermediate scale . Finally, the results support and validate a robust and fruitful operative chain to be implemented for expeditive water discharge assessment in rough and smooth irrigation ditch.
2. Entropy velocity profiles in open channels
The concept of informational entropy as a measure of uncertainty associated to a probability distribution was formulated for the first time in the field of hydraulics by Shannon . The principle of maximum entropy introduces the least-biased probability distribution of a random variable constrained by defined information system as well as the theorem of the concentration for hypothesis testing, introducing the informational entropy theory . A direct evaluation of uncertainty related to the probability distribution of a continuous random variable expressed in terms of entropy, H, is defined as follows
where, p(x) is the continuous probability density function of random variable x.
Using POME, entropy can be maximized through the method of Lagrange multiplier as follows:
in which, m > 0, gi(x) is the ith constraint function and λi is the constrain Lagrange multiplier as a weight in the maximization of entropy.
Chiu [1, 2] applied the concept of entropy to open-channel analysis to model velocity and shear stress distribution as well as sediment concentration. In such a way, the velocity distribution in the probability domain allows to obtain the cross-sectional mean velocity and the momentum and energy coefficients disregarding the geometrical shape of cross sections, which is generally complex in natural channels [2, 13].
Further, an assumption on the probability distribution in the space domain is needed to relate the entropy-based probability distribution to the spatial distribution. Therefore, defining u by the time-averaged velocity placed on an isovelocity curve with the assigned value ξ, the value of u is almost 0 at ξ0, which corresponds to the channel boundary, while u reaches Umax at ξmax, which generally occurs at or below the water surface, depending on the dip-phenomenon. Thus, the velocity u monotonically increases from ξ0 to ξmax and for each value of the spatial coordinate greater than ξ, the velocity is greater than u, and the cumulative distribution function can be written as
Thus, the Shannon entropy of velocity distribution can be written as:
Through a similar procedure, the probability density function of the velocity distribution is obtained by maximizing the Shannon entropy equation
in which, λ0 and λ1 are the Lagrange multipliers and the following constraint equations
Thus, Chiu’s 1D velocity distribution results as:
where M is the dimensionless entropy parameter introduced in the entropy-based derivation [14, 15]. Hence, M can be used as a measure of uniformity of probability and velocity distributions. The value of M can be determined by the mean, Um, and the maximum velocity values are derived from the following equation:
Φ(M) is a relevant parameter which contains relevant information about the flow field asset: the mean velocity value, the location of the mean velocity [14, 15, 16], and the energy coefficient [14, 16] can be obtained from
Eq. (10), in fact, represents the fundamental relationship, from an applied point of view, of the entropy velocity distribution and the assessment of the entropy parameter passing through the knowledge of the ratio between mean and maximum velocities,
In order to identify the dependence of
Even the maximum velocity plays an important role in the flow dynamics, and more than it magnitude, a relevant aspect is related to the position of the maximum velocity inside the flow domain. That is, the location of maximum velocity from the channel bottom,
In this context, Moramarco and Singh  identified the ratio between
with α = (
Eq. (14) highlights, indeed, a possible effect of bed roughness on the entropy velocity distribution in open channel flows, which depends on the roughness scale according to . The dependence between the ratio Φ(M) and the relative submergence, D/d, has been widely studied by Greco  using a wide volume of data collected in the field on several cross sections along different rivers and in the laboratory [20, 21, 22], showing values of Φ(M) ranging in the [0.5–0.9] interval.
3. Laboratory measurements in rectangular smooth and rough ditch
The experimental tests were carried out in the Hydraulics Laboratory of Basilicata University, on two free surface rectangular flumes of 9 m length and with a cross section of 0.5 × 0.5 and 1 × 1 m, whose slope can vary from 0 up to 1%. Figure 1 shows pictures about the flume, one of the bed configuration and the flow-meters.
The bed roughness (d) has been modulated between smooth surfaces, with 0.0005 m roughness height, and a rough bottom, obtained with both a sand bed, with a characteristic diameter of 0.002 m and standard deviation , and a set of wood spheres of 0.035 m in diameter.
The measurement reaches were placed at the distance of 4 m from the beginning of the flumes, in order to damp large-scale disturbances and allow a quasi-uniform water depth. In the end section of the flume, a grid was installed to regulate the water depth for each assigned discharge or rather to obtain a small longitudinal variation of the flow depth. The experiments were performed in steady flow conditions for different values of discharge (0.015–0.100 m3/s) and slope (0.05–1%). The measurement cross section was located in the middle of the rough reach in order to observe a fully developed flow, avoiding edge effects. The flow depth was measured by two hydrometers placed at both the beginning and the end of the measurement reach, and the water depth, D, was assumed as the average value. The velocity was acquired through a micro current-meter with a measuring head diameter of 0.01 m, while the water discharge was measured with a concentric orifice plate installed in the feed pipe and on a laboratory weir placed at the end of the flumes, and compared to the value calculated according to the velocity-area method , with a maximum error of around 1–2%. In particular, the adopted velocity-area method must be applied dividing the cross section into a fixed number of verticals and thus, on each vertical, a fixed measurement points are selected. In each point along the vertical, the velocity is acquired in order to compute the mean velocity of the flow along each vertical. Furthermore, the number of measures on each vertical was chosen with respect to the criterion that the difference in velocity between two consecutive points was less than 20%, of the higher measured velocity value, and the points close to the channel bottom and the water surface was fixed according to the size of the micro-current meter.
In such a way, two roughness configurations were enabled:
RRF: rough rectangular flume, with relative submergence ranging in between 1.89 and 6.43; and
SRF: smooth rectangular flume, with relative submergence greater than 50.
Table 1 synthetically reports the ranges of variation of the main parameters observed during the experiments for the RRF and SRF configuration, while Q is the water discharge, D is the water depth, D/d is the relative submergence, B/D is the aspect ratio, and
|Type||Q (mc/sec)||D (m)||D/d||B/D||Φ(M)|
For each configuration and for all the stages explored, a relevant bulk of velocity measurements was collected in order to provide a detailed reconstruction of the flow field allowing to obtain mean,
Figure 2 shows the linear relationship existing between the pairs (
From Figure 2, some useful issues arise. Even if the correlation among homogeneous data is very strong in both cases with R2 greater than 0.95, it is immediately realized a slight different behavior between rough and smooth channels. That is, for the smooth rectangular flow,
Figure 3 clearly outlines such an outcome, showing how the velocity ratio is austerely dependent on relative submergence in case of rough flows, while it is sufficiently uniform for values of D/d > 20. Furthermore, the same picture proposes several literature data collected by other authors during experimental laboratory campaigns carried on smooth and rough flumes [22, 24, 25, 26, 27], plotted and compared to those arising from the here presented research activity. The same Figure 4 immediately deals with the robust correspondence between data sets related to the low rough/smooth flow conditions for which the hypothesis of the constant value of mean-to-maximum velocities ratio might be assumed consistent, at least from an operative point of view for D/d > 20. At the same time, Eq. (14) still remains compelling for D/d < 20, but it needs to be recalibrated and the coefficients
Such a result can be immediately implemented in the operative chain of water discharge assessment, in order to derive the rating curve in a ditch or artificial channel. Furthermore, such knowledge allows us to assess the level of integrity of the channel in terms of sensitive changes in the bottom roughness, may be due to the local deposition of sediment or vegetation.
Furthermore, in case of
Finally, the use of the entropy velocity profile gives a robust feedback in terms of operative assessment of water discharge, due to the easy and immediate evaluation of the
4. Entropy velocity profile approach for rating curve assessment
The wide bulk of measurements obtained through the laboratory experiments allows us to perform a robust analysis in order to obtain suitable information for the use in the operative chain of water discharge assessment as well as in numerical flow dynamics modeling in regular open channel flow.
In Eq. (10), the mean velocity can be evaluated using Manning’s formula:
where n is the Manning’s roughness, R is the hydraulic radius, and Sf is the energy slope.
To determine the maximum velocity of the cross-section,
is the shear velocity (
The location of the maximum velocity, supporting the dip-phenomenon hypothesis, can be obtained by differentiating Eq. (16) and equating d
Experimental studies [2, 3, 4, 5, 6, 7, 8, 9] have shown that, for channels at different shapes of the cross-section, the velocity maximum is below the free surface around the 20–25% of the maximum depth. Thus, considering
From this latter equation, a new formulation of Manning’s roughness,
which takes into account the variation of a flow hydraulic and geometric characteristics following the change of the water discharge. Eq. (20) computes Manning’s roughness once the values of
Figure 4 shows the correspondence between
Defining the standard error, Se, as suggested by the ISO 1100-2 , through the following relationship:
where N is the number of available measures, the computed Se is permanently less than 5% for the rectangular rough flow (RRF), while increases up to 15%, with a generalized overestimation, in case of smooth rectangular flow. In both cases, the results support the use of this expeditive methodology in the chain of operative procedures leading a good assessment of the rating curve.
The use of a rating curve formulation derived from the entropy velocity theory complained to the assumption of a constant value of the dip velocity and taking into account the variables describing the geometric and hydraulic characteristics of a rectangular ditch, should allow us the improvement of water discharge assessment.
This approach was tested, in a first phase, on a suitable data set of water discharge measures collected in the laboratory on both rough and smooth rectangular cross section proposing practical and common flow conditions.
The rating curve evaluation, derived for the rough rectangular flow, underlines a standard error less than 5%, generally, favoring an expeditive assessment of the flow stage with a sufficient level of reliability, while such an error increase up to 15% in case of smooth cross section.
Chiu C-L. Entropy and probability concepts in hydraulics. Journal of Hydraulic Engineering, ASCE. 1987; 113(5):583-600
Chiu C-L. Velocity distribution in open channel flow. Journal of Hydraulic Engineering, ASCE. 1989; 115(5):576-594
Greco M, Mirauda D, Volpe Plantamura A. Manning’s roughness through the entropy parameter for steady open channel flows in low submergence. Procedia Engineering, ISSN 1877-7058, Published by Elsevier Ltd. 2014; 70:773-780. DOI: 10.1016/j.proeng.2014.02.084. 2-s2.0-84899680421
Greco M, Mirauda D. An entropy based velocity profile for steady flows with large-scale roughness. In: Lollino G, Arattano M, Rinaldi M, Giustolisi O, Marechal JC, Grant G, editors. Engineering Geology for Society and Territory. Vol. 3. Cham: Springer International Publishing; 2015. pp. 641-645. Print; ISBN: 9 978-331909054-2; 978-331909053-5, DOI: 10.1007/978-3-319-09054-2_128, WOS: 000358990300128, EID: 2-s2.0-84944599725
Mirauda D, Greco M, Moscarelli P. Practical method for flow velocity measurements in fluvial sections. WIT Transactions on Ecology and the Environment. 2011; 146:355-366. ISBN:978-1-84564-516-8 e ISSN 17433541. DOI: 10.2495/RM110301
Moramarco T, Saltalippi C, Singh VP. Estimation of mean velocity in natural channels based on Chiu’s velocity distribution equation. Journal of Hydrologic Engineering, ASCE. 2004; 9(1):42-50
Moramarco T, Singh VP. Formulation of the entropy parameter based on hydraulic and geometric characteristics of river cross sections. Journal of Hydrologic Engineering. 2010; 15(10):852-858
Yang SQ, Tan SK, Lim SY. Velocity distribution and dip-phenomenon in smooth uniform open channel flows. Journal of Hydrologic Engineering, ASCE. 2004; 130(12):1179-1186
Greco M. Effect of bed roughness on 1-D entropy velocity distribution in open channel flow. Hydrology Research. 2015; 46(1):1-10. ISSN 1998-9563, EID: 2-s2.0-84926637541. DOI: 10.2166/nh.2013.122
Bathurst JC, Li RM, Simons DB. Resistance equation for large-scale roughness. Journal of the Hydraulics Division. 1981; 107(HY12):1593-1613
Shannon CE. A mathematical theory of communication. The Bell System Technical Journal. 1948; 27:623-656
Tsallis C. Possible generalization of Boltzmann-Gibbs statistics. Journal of Statistical Physics. 1988; 52(1–2):479-487
Greco M, Moramarco T. Influence of bed roughness and cross section geometry on medium and maximum velocity ratio in open-channel flow. Journal of Hydraulic Engineering. 2016; 142(1):06015015. DOI: 10.1061/(ASCE)HY.1943-7900.0001064. Article number 0601501 EID: 2-s2.0-84950267443
Chiu C-L, Said CA. Maximum and mean velocity and entropy in open channel flow. Journal of Hydraulic Engineering, ASCE. 1995; 121(1):26-35
Luo H, Singh V. Entropy theory for two-dimensional velocity distribution. Journal of Hydrologic Engineering. 2011; 16(4):303-315
Chiu CL, Tung NC. Maximum velocity and regularities in open-channel flow. Journal of Hydraulic Engineering. 2002; 128(8):803-803
Xia R. Relation between mean and maximum velocities in a natural river. Journal of Hydrologic Engineering, ASCE. 1997; 123(8):720-723
Nezu I, Nakagawa H. Turbulence in open-channel flows. In: Balkema. The Netherlands: Rotterdam; 1993
Rouse H. Critical analysis of open-channel resistance. Journal of the Hydraulics Division. 1965; 91(4):1-23
Greco M, Mirauda D. Entropy parameter estimation in large-scale roughness open channel. Journal of Hydrologic Engineering. 2015; 20(1–2). DOI: 10.1061/(ASCE)HE.1943-5584.0001009. ISSN: 1084-0699 eISSN: 1943-5584, Article number 04014047, EID 2-s2.0-84921442470
Greco M. Entropy-based approach for rating curve assessment in rough and smooth irrigation ditches. Journal of Irrigation and Drainage Engineering. 2016; 142(3):04015062. DOI: 10.1061/(ASCE)IR.1943-4774.0000986. EID 2-s2.0-84958758546
Nikmehr S, Farhoudi J. Estimation of velocity profile based on Chiu’s equation in width of channels. Research Journal of Applied Sciences, Engineering and Technology. 2010; 2(5):476-479
ISO 748. 2007 Hydrometry—Measurement of liquid flow in open channels using current-meters or floats-ISO/TC 113/SC 1 Velocity area methods
Guy HP, Simons DB, Richardson EV. Summary of Alluvial Channel data from flume experiments, 1956–1961, Geological Survey Professional Paper, 1966, 462-I. US Government Printing Office, Washington, DC
Steffler PM, Rajaratnam N, Peterson AW. LDA Measurements of Mean Velocity and Turbulence Distribution in a Smooth Rectangular Open Channel. Edmonton: Department of Civil Engineering., University of Alberta; 1983
Bortz K. Parameter estimation for velocity distribution in open channel flow [M.S. Thesis]. Department of Civil Engineering, University of Pittsburgh: Pittsburgh, PA; 1989
Guo ZR. Personal communication in Chiu and Hsu, 2006. In: Southeast China Environmental Science Institute. China: Yuancun; 1990
ISO 1100-2: 2010 Hydrometry—Measurement of liquid flow in open channels—Part 2: Determination of the stage-discharge relation-ISO/TC 113/SC 1 Velocity area methods