WST 086071681
WST 086071681
WST 086071681
287
Eyad Abushandi
Civil Engineering Stream, Faculty of Engineering, Sohar University, Sohar, Oman
E-mail: eabushandi@gmail.com
EA, 0000-0001-6160-6186
ABSTRACT
Open channel flow and velocity behaviour presents a dilemma for drainage systems designers where hydrologic parameters are fluctuating in
space and time. The experience of having extreme flash floods almost every year is flustering the need for understanding the flow behaviour
at different altitudes. In this study, open channel experimental efforts were made to model flow and velocity profiles. The three-dimensional
(3D) open channel flow and velocity profiles are generated at two types of roughness namely steel (smooth) and concrete (rough beds). The
experiments included different slope gradients and flow measurements at different distances. The channels slope ranges between 0 and
4.7% with an interval of 0.2%. Multiple linear regression (MLR) was applied to quantify the flow for longer distance while Kriging extrapolation
proxy was used to generate 3D surfaces of flow and velocity. The results showed that the flow in concrete channel is decreasing by moving to
the end of channel due to higher frictional resistance while it is rising up for the steel channel. In average, the velocity has been increased by
7.4% for steel and 6.1% for concrete channels at a changing slope.
Key words: Kriging extrapolation, material roughness, MLR, open channel, slope
HIGHLIGHTS
• Experimental validations of two open channel types are important to predict flow rates and velocities for different lengths with complex
slope and roughness conditions.
• The Pearson correlation shows a strong relation between the flow rates and the slope.
• The results showed the applicability of multiple linear regression model and 3D kriging model for analyzing real flow structure in a longer
distance.
INTRODUCTION
The status of urban drainage facilities as a part of the integrated infrastructure system is still under development in many
regions around the world. This, however, depends on the level of development and the society awareness of the drainage
system importance. As the accumulated flow reaches the catchment outlet, it may destroy the infrastructure of located city
or traffic roads, unless it is accommodated in a proper designed drainage system. The risk of damage caused by this flow
is basically based on water magnitudes, rainfall intensity, and storm duration. Therefore, a proper design of drainage network
at the catchment outlet section is important to avoid life losses and damages of infrastructure. Flow and velocity distributions
in drainage system are needed for a wide range of applications in hydrologic studies such as flood early warning, sediment
transport, urban drainage system design and maintenance. The research in this field started even before Darcy, Weisbach, and
Manning when ancient nations conveyed rainwater from the mountains into resident areas. The flow rate and velocity in
open channel flow have been studied physically and mathematically by numerous researchers, e.g. Alawdi & Prasad
(2018); Basu (2019); Benoumessad et al. (2014); Li Zeng & Bai (2020); Wang et al. (2019); Welderufael et al. (2019); and
Abushandi & Al Sarihi (2022). In addition to those studies, several experiments have been conducted to model the turbulent
flow and velocity in both two- and three-dimensional (2D and 3D) cases using Reynolds-averaged Navier-Stokes equations
(RANS) (Kang & Choi 2006; Wu et al. 2018; Mohotti Wijesooriya & Dias-da-Costa 2019; Welderufael et al. 2019). They
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/).
calculated the kinetic energy and turbulent dissipation at a fixed Reynolds number. Apart from some pilot studies at a catch-
ment scale, e.g. Abdel-Fattah et al. (2017); Abushandi & Merkel (2013); Aceves & Fuamba (2016); and Kordilla et al. (2013);
the flow magnitudes have been estimated based on empirical laboratory investigation. The flow attributes in the channels are
directly connected to channel slope and surface roughness. Therefore, several investigations in recent years have been con-
ducted to determine the flow and vertical distributions of velocities for turbulent fluids on both ‘smooth’ and ‘rough’ surfaces
(Kordilla et al. 2013; Singh Raushan & Debnath 2018; Bormashenko 2019). Computation findings for velocity distribution in
turbulent smooth-wall open channel flows were assessed and compared to experimental data (Welderufael et al. 2019). The
results concluded that the model based on the Reynolds-averaged Navier-Stokes (RANS) equation provides accurate assess-
ments. An ordinary differential equation (ODE) for velocity distribution in open channel flows is presented by Absi (2011)
based on an analysis of the Reynolds-averaged Navier–Stokes equations and a log-wake modified eddy viscosity distribution.
The method allows predicting the velocity-dip-phenomenon, i.e. the maximum velocity below the free surface. Pal et al.
(2016) attempt to model the impact of particles on a stream wise velocity profile at several flow conditions using mathemat-
ical modelling. Their method provided an accurate assessment of velocity vertical and horizontal profiles. In addition, the
assumption of steady-state, turbulent flow, or velocity distributions were the core for computational methods (Shi & Yu
2015; Govorukhin & Zhdanov 2018; Kimiaghalam et al. 2018; Luo et al. 2018; Zhang et al. 2018; Zhao et al. 2019). The
estimation of sediment transport rates in open channels is of great interest to engineers in addition to slope and vegetation
drag forces because of the practical importance of sediment transport in changing flow and velocity rates (Zhong et al. 2015;
Palucis et al. 2018; Elgueta-Astaburuaga & Hassan 2019; Tan & Yuan 2019).
The sampling interval or frequency has taken a place for accurate simulation (Ruonan et al. 2016). However, the state of art
technology in measurements and simulation were developed to improve the results such as unmanned aerial vehicle
(UAV)-based tracer tests using RGB (red, green, blue) images (Baek et al. 2019). Researchers have studied the velocity and
flow profiles intensively in open and close conduits. Kirkgoz (1989) measured the velocity using a laser-doppler anemometer
for a fully developed, rectangular, subcritical open channel flow on smooth and rough beds. The calculated velocities showed
an increasing tendency as the wall roughness increases. The representation of the overall data in terms of law-of-the-wall dis-
tribution seems reasonable; however, the velocity-defect distribution is not satisfactory because of small flow depth. A couple
of years later Kirkgoz & Ardlboglu (1997) approved better results, they found a strong linear relationship between the channel
length and flow depths. Castellarin (2014) used Kriging interpolation method to spatially estimate flow for ungagged basin
which dominated by same factors. It is a fact that Kriging method gives better result in a homogenous channel.
The primary objectives of this research are to provide a predictive tool for flow rates and velocities for complex real-world
conditions with different lengths and surface roughness. Furthermore, spatially describe flow and velocity distribution and
test slope and surface roughness affecting flow and velocity. For this purpose, two open channels were designed, and exper-
imental laboratory efforts have been conducted.
Experimental set-up
The steel and concrete channels were designed and manufactured at Sohar University/Hydraulics Laboratory with the
dimensions of 4 m long, 0.2 m wide and 0.3 m deep. Figure 1 shows the vertical and horizontal dimensions of the two
open channels and velocity measurement points. In addition, a 3D scanning was conducted to show the surface of a concrete
open channel. However, the dimensions of both channels are the same, in particular, the width of 20 cm, and will have the
same surface tension effects. Therefore, the effect on flow generated by scale is neglected. A small, electrically powered cen-
trifugal pump was used to provide water to a closed water circle with a flow rate of 0.5 m3/s. The roughness values used for
concrete is 0.2 mm while steel has a value of 0.02 mm (Native Dynamics 2020).
The flow rate tests of the two conveyance systems were measured at several lengths and slopes to ensure materials suit-
ability for the channel functionality. The flow rates were measured using float method and compared to ultrasonic
Figure 1 | Channels major elements, dimensions, and the points of measurements, and 3D scanning of the concrete channel using 3D
EinScan HX 3D-Scanner.
flowmeter to ensure the accuracy. The steel channel bed was thought to be suitable for a smooth surface wall experiment
(Figure 2), while the concrete is considered for the rough surface (Figure 3). The steel channel was manufactured by bending
the sheet without using any inside welding to avoid minor losses of energy. This, however, helps to consider the major losses
water flow area and avoid any affect in the velocity values accuracy; one long bolt was fixed as a jack to control the channel
slope. Furthermore, the most common used material to construct a drainage system nowadays is the concrete. Therefore, it
was an important aspect to test the flow rate through a designed concrete channel. The channels were designed to meet the
real-world adjustment. The models were performed for a variety of slope grades while always maintaining a constant flowrate
from the source. The system can be applied for agricultural and domestic purposes to collect harvest rainwater at different
slope gradients. It is a hand on the development of appropriate techniques for designing drainage system to accommodate
the unexpected flood magnitudes.
Figure 2 | Designed smooth steel channel to conduct flow rates and velocity rates.
The design is similar to the internal height and width of the steel channel in order to get a constant cross-sectional area
while conducting the flow tests and the regression model. Flow system was designed including a tank with 100 L capacity
and the pipes system was connected to get a water cycle. The slope interval for the two channels was between 0.0% and 4.7%.
Mathematical representation
There are several equations used for this research including Kriging interpolation method, Pearson correlation, multi-linear
regression model and Nash-Sutcliffe efficiency equation.
Kriging interpolation method was used to 3D interpolate the observed flow and velocity. Kriging is a method of interp-
olation named after the South African mining engineer D. G. Krige, who created the technique while investigating
correlation parameters in 1951 (Ryu et al. 2002). The method takes into considerations spatial autocorrelation of observed
records and the dimension of the channel.
Kriging is the same as inverse distance weighing method (IDW) that considers the value and location of an observed point
to find unmeasured values. The general formula of Kriging as follows:
X
N
^ 0) ¼
Z(s li Z(si )
i¼1
where:
Z(si ) is the measured value of flow or velocity at the location of 1–4 m
li is the unknown weight for flow or velocity at the location at i location
s0 is the predicted location
N is the number of measured values
To assess the efficiency of Kriging method results Pearson correlation was used in order to evaluate the relationship
between the two different flows based on the following formula:
P
(x x )(y y)
r ¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
(x x )2 (y y)2
The values of Pearson correlation are always between 1 and 1, and if x and y are not related, the correlation is equal to
zero. Furthermore, multiple linear regressions model was used to model the relationships between slope, friction, and flow or
velocity by fitting a linear equation to observed data. General formula of multiple linear regressions model is given by:
where
P
^ i )2
(qi q
Ef ¼ 1 P
)2
(qi q
It is necessary to apply Pearson’s correlation to show the relationship between discharge and slope. Although this obser-
vation has been known for nearly two centuries, there is a need to quantify such relationship and understand the flow
behaviour for different surfaces and slopes. Based on observations comparisons, Pearson correlation between the flow at
different distance and slope in concrete and steel showed a significant relationship between each two variables with an aver-
age value of r equal to 0.86 and 0.99 for concrete and steel, respectively. However, when the flow in the concrete channel is
declining, it has been rising in the steel channel. While r values between the slope and the concrete channel was 0.88, and for
the steel channel flow was 0.98. These values reflect how strong the influence of slope is on the flow. Tables 1 and 2 showed
detailed values of correlation between slope gradients and flow rates at different points of measurements. Generally, the cor-
relation is significant, thus, MLR can be successfully applied.
Table 1 | Pearson correlation for the concrete channel flow in response to slope
Variable Slope Q1 Q2 Q3 Q4
Variables Slope Q1 Q2 Q3 Q4
To simulate flow rates values at 1, 2, 3, and 4 m within different length of channels, an adjustment was performed based on
multiple linear regressions to model additional length of channels, Q5 represent the value of modelled flow for one additional
meter:
Concrete channel MLR:
Figure 8 | 3D wireframe of steel channel flow and velocity distributions (smooth surface) using Kriging Interpolation method at average slope
0.23%.
Figure 9 | 3D wireframe of concrete channel flow rate and velocity distributions (rough surface) using Kriging Interpolation method at slope
at average slope 0.23%.
However, the purpose of developing MLR equations is to simulate the flow in a further distance with slope and surface
roughness variability.
Based on the statistical analysis, the accuracy of the modelled values is around 97%. In addition, the Nash-Sutcliffe effi-
ciency (Ef ) was around 88.4 and 83.1 for concrete and steel channels, respectively.
In a different manner, 3D numerical models for calculating flow and velocity profiles in open channels were based on Kri-
ging formula. The models showed the changes with respect to the distance (Figures 8 and 9). The profile shows the values of
unmeasured points. Results indicated that Kriging method has the potential to spatially visualize the flow and velocity. In
addition, using MLR equation can further produce extrapolation for longer distances.
The performance of computed velocities and flows are rather good in comparison to observed records. Kriging method was
tested against several slopes for similar flow rate. All tests provided very good correspondence between modeled and
observed measured profiles. Generally, Kriging extrapolation method seems to be an appropriate tool for analysing open
channel spatial data. The accuracy of quantified estimation variance average reached 84.6 and 98.5% for concrete and
steel channels, respectively. Despite the interpolation method used, the advantage of velocity and flow is that the effects
of roughness and slope over the entire channel distance have been considered at same time, instead of temporal (points)
effect. However, this research was conducted at laboratory conditions. Real-world open channels have much more compli-
cated mechanisms due to sedimentation and slope fluctuation along the channel. In addition, natural open channels in arid
areas (wadis) are dominated by rocks, and limited vegetation cover. This, however, makes most modelling flow and velocity
profiles in vegetated channels based on a single physical concept (Nikora et al. 2013). Furthermore, modelling natural chan-
nels includes incident conditions including soil moisture content (Abushandi & Al Sarihi 2022).
CONCLUSION
Experimental validations of smooth and rough open channels are important to predict flow rates and velocities for different
lengths with complex slope and roughness conditions. As the slope is increasing flow rate and velocity of both type channels
are increasing. On average, the velocity has been increased by 7.4% for steel and 6.1% for concrete channels at each increas-
ing slope interval. While the velocity tendency over concrete surface was higher due to frictional resistance. However, further
studies are required to investigate real-world open channel hydrologic behavior including sediments transport and vegetation
cover influence. A hydraulic jump has been created in a steel channel while the flow rates and velocity are declining by
moving to the end of channel due to friction. MLR model can be applied as far as a linear relationship is found between differ-
ent variables. The measurements were conducted for 20 different slope gradients to understand numerically how this
parameter is influencing water flow and velocity in open channels. The three-dimensional representations provided a clear
image flow and velocity behavior. Further measurements will be conducted to understand vegetation cover and settlements
impacts on flow and velocity profiles.
ACKNOWLEDGEMENTS
The research leading to these results has received funding from the Research Council (TRC) of the Sultanate of Oman under
the Open Research Grant Program #BFP/RGP/EBR/19/164.
CONFLICT OF INTEREST
The authors declare there is no conflict.
REFERENCES
Abdel-Fattah, M., Saber, M., Kantoush, S., Khalil, M., Sumi, T. & Sefelnasr, A. 2017 A hydrological and geomorphometric approach to
understanding the generation of wadi flash floods. Water 9, 553.
Absi, R. 2011 An ordinary differential equation for velocity distribution and dip-phenomenon in open channel flows. Journal of Hydraulic
Research 49 (1), 82–89. doi:10.1080/00221686.2010.535700.
Abushandi, E. & Al Sarihi, M. 2022 Flash flood inundation assessment for an arid catchment, case study at Wadi Al Jizzi, Oman. Water
Practice and Technology wpt2022043. doi:10.2166/wpt.2022.043.
Abushandi, E. & Merkel, B. 2013 Modelling rainfall runoff relations using HEC-HMS and IHACRES for a single rain event in an arid region
of Jordan. Water Resources Management 27 (7), 2391–2409.
Aceves, M. & Fuamba, M. 2016 Methodology for selecting best management practices integrating multiple stakeholders and criteria. Part 2:
case study. Water 8, 56.
Ahadi, M., Bergstrom, D. J. & Mazurek, K. A. 2019 Application of the two-fluid model to prediction of sediment transport in turbulent open
channel flow. Physics and Chemistry of the Earth, Parts A/B/C 113, 73–82. https://doi.org/10.1016/j.pce.2019.06.001.
Alawdi, W. & Prasad, T. 2018 Comparison of performance of simplified RANS formulations for velocity distributions against full 3D RANS
model. International Journal of Hydraulic Engineering 7 (2), 33–42.
Baek, D., Seo, I. W., Kim, J. S. & Nelson, J. M. 2019 UAV-based measurements of spatio-temporal concentration distributions of fluorescent
tracers in open channel flows. Advances in Water Resources 127, 76–88. https://doi.org/10.1016/j.advwatres.2019.03.007.
Basu, S., 2019 Chapter III – Open-Channel Flow Measurement. In: Plant Flow Measurement and Control Handbook (Basu, S. ed.). Academic
Press, New York, NY, USA, pp. 257–331.
Benoumessad, K., Ilhem, K., Ali, F. & Abdelbaki, D. 2014 3D simulation of velocity profile of turbulent flow in open channel with complex
geometry. Eight International Conference on Material Sciences (CSM8-IMS5) 55, 119–128.
Bormashenko, E. 2019 Apparent contact angles for reactive wetting of smooth, rough, and heterogeneous surfaces calculated from the
variational principles. Journal of Colloid and Interface Science 537, 597–603. https://doi.org/10.1016/j.jcis.2018.11.068.
Castellarin, A. 2014 Regional prediction of flow-duration curves using a three-dimensional kriging. Journal of Hydrology 513, 179–191.
https://doi.org/10.1016/j.jhydrol.2014.03.050.
Elgueta-Astaburuaga, M. A. & Hassan, M. A. 2019 Sediment storage, partial transport, and the evolution of an experimental gravel bed under
changing sediment supply regimes. Geomorphology 330, 1–12. https://doi.org/10.1016/j.geomorph.2018.12.018.
Govorukhin, V. & Zhdanov, I. 2018 Steady-state flows of inviscid incompressible fluid and related particle dynamics in rectangular channels.
European Journal of Mechanics – B/Fluids 67, 280–290. https://doi.org/10.1016/j.euromechflu.2017.09.016.
Huai, W., Hu, Y., Zeng, Y. & Han, J. 2012 Velocity distribution for open channel flows with suspended vegetation. Advances in Water
Resources 49, 56–61. https://doi.org/10.1016/j.advwatres.2012.07.001.
Kang, H. & Choi, S.-U. 2006 Turbulence modeling of compound open-channel flows with and without vegetation on the floodplain using the
Reynolds stress model. Advances in Water Resources 29, 1650–1664. https://doi.org/10.1016/j.advwatres.2005.12.004.
Kimiaghalam, N., Clark, S. P. & Dow, K. 2018 Flow characteristics of a partially-covered trapezoidal channel. Cold Regions Science and
Technology 155, 280–288. https://doi.org/10.1016/j.coldregions.2018.08.016.
Kirkgoz, S. 1989 Turbulent velocity profiles for smooth and rough open channel flow. Journal of Hydraulic Engineering 115 (11), 1543–1561.
Kirkgoz, S. & Ardlboglu, M. 1997 Velocity profiles of developing and developed open channel flow. Journal of Hydraulic Engineering
123 (12), 1099–1105.
Kordilla, J., Tartakovsky, A. M. & Geyer, T. 2013 A smoothed particle hydrodynamics model for droplet and film flow on smooth and rough
fracture surfaces. Advances in Water Resources 59, 1–14. https://doi.org/10.1016/j.advwatres.2013.04.009.
Li, Q., Zeng, Y.-h. & Bai, Y. 2020 Mean flow and turbulence structure of open channel flow with suspended vegetation. Journal of
Hydrodynamics 32 (2), 314–325. doi:10.1007/s42241-020-0020-6.
Luo, H., Singh, V. & Schmidt, A. 2018 Comparative study of 1D entropy-based and conventional deterministic velocity distribution equations
for open channel flows. Journal of Hydrology 563, 679–693. https://doi.org/10.1016/j.jhydrol.2018.06.010.
Mohotti, D., Wijesooriya, K. & Dias-da-Costa, D. 2019 Comparison of Reynolds Averaging Navier-Stokes (RANS) turbulent models in
predicting wind pressure on tall buildings. Journal of Building Engineering 21, 1–17.
Native Dynamics. 2020 Absolute Roughness of Pipes Materials.
Nikora, N., Nikora, V. & O’Donoghue, T. 2013 Velocity profiles in vegetated open-channel flows: combined effects of multiple mechanisms.
Journal of Hydraulic Engineering 139 (10), 1021–1032.
Pal, D., Jha, S. K. & Ghoshal, K. 2016 Velocity lag between particle and liquid in sediment-laden open channel turbulent flow. European
Journal of Mechanics – B/Fluids 56, 130–142. https://doi.org/10.1016/j.euromechflu.2015.11.003.
Palucis, M. C., Ulizio, T. P., Fuller, B. & Lamb, M. P. 2018 Flow resistance, sediment transport, and bedform development in a steep gravel-
bedded river flume. Geomorphology 320, 111–126.
Ruonan, B., Liekai, C., Xingkui, W. & Danxun, L. 2016 Comparison of ADV and PIV measurements in open channel flows. Procedia
Engineering 154, 995–1001. https://doi.org/10.1016/j.proeng.2016.07.588.
Ryu, J. S., Kim, M. S., Cha, K. J., Lee, T. H. & Choi, D. H. 2002 Kriging interpolation methods in geostatistics and DACE model. KSME
International Journal 16 (5), 619–632.
Shi, H. & Yu, X. 2015 An effective Euler–Lagrange model for suspended sediment transport by open channel flows. International Journal of
Sediment Research 30, 361–370v.
Singh, S. K., Raushan, P. K. & Debnath, K. 2018 Combined effect of wave and current in rough bed free surface flow. Ocean Engineering 160,
20–32. https://doi.org/10.1016/j.oceaneng.2018.04.055.
Tan, W. & Yuan, J. 2019 Experimental study of sheet-flow sediment transport under nonlinear oscillatory flow over a sloping bed. Coastal
Engineering 147, 1–11. https://doi.org/10.1016/j.coastaleng.2019.01.002.
Wang, N., Lu, H., Xu, J., Guo, X. & Liu, H. 2019 Velocity profiles of granular flows down an inclined channel. International Journal of
Multiphase Flow 110, 96–107. https://doi.org/10.1016/j.ijmultiphaseflow.2018.09.002.
Welderufael, M., Absi, R. & Mélinge, Y. 2019 Assessment of models for velocity distribution in turbulent smooth-wall open channel flows.
ISH Journal of Hydraulic Engineering 1–11. doi:10.1080/09715010.2019.1677180.
Wu, Y., Liu, Z., Chen, Y. & Li, M. 2018 Investigation of velocity distribution and turbulence characteristics in subcritical circular open
channel flows using a modified Reynolds stress model. Journal of Hydro-Environment Research 19, 68–77. https://doi.org/10.1016/
j.jher.2018.02.002.
You, B.-H., Jeong, Y. H. & Addad, Y. 2019 Assessment of advanced RANS models ability to predict a turbulent swept liquid metal flow over a
wire in a channel. Nuclear Engineering and Design 353, 110206. https://doi.org/10.1016/j.nucengdes.2019.110206.
Zhang, J., Zhong, Y. & Huai, W. 2018 Transverse distribution of streamwise velocity in open-channel flow with artificial emergent vegetation.
Ecological Engineering 110, 78–86. https://doi.org/10.1016/j.ecoleng.2017.10.010.
Zhao, X., Bolognin, M., Liang, D., Rohe, A. & Vardon, P. J. 2019 Development of in/outflow boundary conditions for MPM simulation of
uniform and non-uniform open channel flows. Computers & Fluids 179, 27–33. https://doi.org/10.1016/j.compfluid.2018.10.007.
Zhong, D., Zhang, L., Wu, B. & Wang, Y. 2015 Velocity profile of turbulent sediment-laden flows in open-channels. International Journal of
Sediment Research 30, 285–296. https://doi.org/10.1016/j.compfluid.2018.10.007.
First received 16 May 2022; accepted in revised form 30 August 2022. Available online 3 September 2022