Hydraulic Coeeficients
Hydraulic Coeeficients
Hydraulic Coeeficients
5-1991
S. C. Finkner
University of Nebraska-Lincoln
Gilley, John E. and Finkner, S. C., "Hydraulic Roughness Coefficients as Affected by Random Roughness" (1991). Biological Systems
Engineering: Papers and Publications. 122.
https://digitalcommons.unl.edu/biosysengfacpub/122
This Article is brought to you for free and open access by the Biological Systems Engineering at DigitalCommons@University of Nebraska - Lincoln. It
has been accepted for inclusion in Biological Systems Engineering: Papers and Publications by an authorized administrator of
DigitalCommons@University of Nebraska - Lincoln.
Gilley & Finkner in Transactions of the ASAE 34 (1991)
overland flow is less pronounced as the depth of overland assumed to be approximately equal to flow depth. For these
flow becomes greater. situations:
The objectiv巳of this investigation was to develop
Rn三旦 (6)
regression equations for estimating hydraulic roughness v
coefficients using random roughness and Reynolds number
as independent variables. Relationships were identified for where flow rate per unit width , q , is given as:
predicting both Darcy-Weisbach and Manning hydraulic
roughness coe征icients. The equations were derived using
random roughness values varying from 6 to 32 mm , and q=? (7)
Reynolds numbers ranging from 20 to 6000.
f= 8gRS
2
(1) y=-~
---
Vb
(9)
V
e--==--
whHugSVR In this study, water depth was determined indirectly using
'也略
cvovd.tsoc
d
aanh tvda
α 町
oke-en
nm-UB
u'au
晦 ed v
gh
.川
3 e
z
bUAUAU
, equation 9, and measurements of Q, V, and b.
巳ivρd
uv出
d
白
0
3
gIven as:
--nu)
,,
R=A
E飞
(2)
P
where A is cross-sectional flow area and P is wetted The Chezy hydraulic roughness coefficient can be
penmeter. determined directly from the Darcy-Weisbach hydraulic
roughness coe旺icient using the relationship:
The Manning hydraulic roughness coefficient , n , is
given as (Chow, 1959):
c=[制
l l)
/,‘、
2/LI/2
R-'·S
n= 一-一 (3)
V Information on existing flow characteristics is needed to
relate the Manning hydraulic roughness coefficient to
Manning and Darcy- Weisbach hydraulic roughness either Chezy or Darcy-Weisbach hydraulic roughness
coe证icients can be related using the following equation: coefficients , since:
,,,?&
I /6
F 4 1 l t t t t」
R一
句4d
、 ·1
g1- R
l l I Bι叶E E E E E 」
n= n=-一一 (1 2)
句
(4)
E
C
tt1
Rn= 立且 (5)
v
EXPERIMENTAL PROCEDURES
Field tests were conducted at the University of Nebraska
where v is kinematic viscosity. Kinematic viscosity can be Rogers Memorial Farm located in Lancaster County,
determined directly from water temperature. approximately 18 km east of Li ncoln , NE. The Sharpsburg
silty clay loam at the site (fine , montmorillonitic , mesic
For broad sheet flow conditions where flow width , b , is Typic Ar giudolls) formed on loess under prairi巳vegetation.
much greater than flow depth , y, hydraulic radius can be Average slope at the location was 6.4%.
The experimental design consisted of two randomized artificial tmf carpe t. Runoff was diverted into an HS flume
complete blocks , with the first block being located with a stage recorder for measurement of flow rate.
immediately upslope from the second. Each experimental Flow velocity was determined using dye tracing
block consisted of six tillage operations performed at techniques. Approximately 0.2 L of fluorescent dye was
random locations within the block. The tillage operations uniformly injected across the width of the plot , 0.76 m
included an anhydrous applicator, chisel plow , disk , field upslope from the lower boundary. A peristaltic pump was
cultivator, moldboard plow, and planter. These implements used to continuously withdraw flow at four points spaced
were chosen to provide a wide range of random roughness equally along the collection trough. Discharge was then
conditions. circulated through a flourometer and a scale on the
instrument provided a visual display of dye concentration.
SITE PREPARATION Average time of travel was calculated as the length of time
Existing wheat residue was first removed from the study required for the dye concentration peak to reach the lower
area by burning and hand raking. Selected tillage boundary.
operations were then performed along the contour at the Five measurements of travel time were obtained at each
study site. One-meter-square plots were established within inflow rate. The mean of the five readings was used to
each tillage treatment using galvanized sheet metal borders calculate flow velocity at a particular inflow rate.
for the top and both sides of the plots. A collection trough , Calculated flow velocities were used to determined flow
located at the bottom of the plots , was used to collect depths and corresponding hydraulic roughness coefficients.
runoff. When not in use , the plots were covered with Additional details concerning experimental procedures are
plywood which was placed several centimeters above the given by Finkner (1988).
soil surface. The plywood covering prevented weathering
of the soil surface.
Soil surface stabilization was required to prevent RESULTS AND DISCUSSION
destruction of soil form roughness during test procedures. Random roughness values were first identified for each
After measurements for random roughness were obtained , of six selected tillage operations. Darcy-Weisbach and
the plot surfaces were stabilized using a biodegradable, Manning hydraulic roughness coefficients were then
latex-base soil stabilizer. The stabilizer was sprayed over measured on each of the soil surfaces where random
the entire soil surface using a hand sprayer. The stabilizing roughness values were obtained. Independent variables
material penetrated the soil approximately 5 mm , influencing hydraulic roughness coefficients were
e旺巳ctively binding the soil particles together with a water determined. Finally, regression equations were developed
permeable layer. for estimating Darcy-Weisbach and Manning hydraulic
roughness coefficients from values of random roughness
RANDOM ROUGHNESS and Reynolds number.
Differences in soil surface height were recorded using a
mechanical profile meter. The surface profile meter , similar RANDOM ROUGHNESS
to the device described by Allmaras et al. (1967) , could be Random roughness was calculated using the procedure
easily rolled above the entire plot surface on a rectangular outlined by Allmaras et al. (1967). Table 1 presents random
support frame. The support frame was of variable height roughness measurements obtained in the present study, and
and was leveled in the horizontal plane. The rectangular values reported by Zobeck and Onstad (1987) in a review
frame was supported by four 250-mm steel stakes which of available literature. Random roughness values in the
were securely anchored into the soil to provide a horizontal present investigation ranged from 6 mm for the planter to
ref,巳renee. The upper left corner of each plot border as 32 mm for the moldboard plow treatmen t.
viewed from the bottom of the plots was used as a vertical
bench mark , creating a three-dimensional referencing TABLE 1. Random roughness values for selected
system. tillage operations
The profile meter consisted of a single row of equal Random
length, 3.2 mm diameter steel pins positioned at a spacing Random Roughness
of 6.4 mm. When lowered onto the soil surface , the top of Tillage Roughness* (Present Study)
the pins formed a nearly continuous line which was traced Operation mm 口1 m
onto a strip of paper located behind the pins. The profile Large offset disk 50
meter and frame were oriented so that surface elevations Moldboard plow 32 32
were measured parallel to the contour of the study area. Lister 25
Chisel plow 23 21
Transects were spaced every 50 mm along the slope and Disk 18 16
transect traces were later digitized at 25 mm spacings. A Field cultivator 15 14
total of 629 surface elevations were used for determination Row cultivator 15
of random roughness for each one meter square plo t. Rotary tillage 15
H征row 15
Anhydrous applicator 13 8
FLOW VELOCITY Rod weeder 10
Following surface stabilization with the latex-base soil Planter 10 6
stabilizer, flow was added to the top of each plot at twelve No-till 7
rates ranging from approximately 2 to 300 L/min. Flow Smooth surface 6
inlet energy was dissipated at the top of the plots using an * Zobeck, T. M. and C. A. Onstad (1987).
The random roughness values shown in Table 1 Random roughness , Reynolds number, and slope were
represent best estimates for a particular tillage operation. retained as possible predictors. In an analogy to pipe flow,
Differences in soil texture , water content at time of tillage , random roughness provided a measure of the physical
or tillage depth may affect surface conditions. In addition , roughness of the flow boundary and Reynolds number
variations in the physical characteristics of the tillage furnished a flow property. Bed slope was included because
implements may result in different random roughness it has previously been found to influence hydraulic
values. roughness coe旺icients (Issard , 1944; Emmett, 1970; and
The anhydrous applicator and planter caused little Yoon, 1970).
disturbance to the relatively smooth surface which existed A simple multiplicative relationship which included
at the study site. Random roughness values for these two random roughness , Reynolds number , and slope as
operations were less than those previously reported. For the independent variables was tested. The effect of adding each
other tillage operations , random roughness measurements of these three variables into the prediction equations was
obtained in the present study were in close agreement with evaluated using multiple linear regression analysis. Only
values reported by Zobeck and Onstad (1987). random roughness and Reynolds number were found to be
The addition of rainfall may serve to reduce random significant at the 0.10 probability leve l. Therefore , only
roughness. To quantify this reduction , a relative random these two variables were considered in subsequent
roughness term (RRR) was defined by Zobeck and Onstad analyses.
(1 987) as:
DARCy-WEISBACH HYDRAULIC ROUGHNESS COEFFICIENT百,
Darcy-Weisbach hydraulic roughness coefficients at
RRR= 旦旦 (13)
RR o varying Reynolds numbers for the moldboard plow and
planter treatments are presented in figure 1. The trends
presented for the moldboard plow and planter operations
where RR is the random roughness of a surface following are also characteristic of the other experimental treatments.
rainfall , and RR o is random roughness immediately after In gene时, hydraulic roughness coefficients can be seen to
tillage. From published data on relative random roughness , decrease with greater Reynolds number.
Zobeck and Onstad (1987) developed the following The moldboard plow and planter treatments produced
equa tIon: the largest and smallest random roughness values ,
respectively. The largest hydraulic roughness coefficients
-0.026 cumulative rainfall usually occurred on those plots with the greatest random
RRR=0.8ge (14) roughness. The planter treatments with relatively low
random roughness values produced the smallest hydraulic
where cumulative rainfall is given in centimeters. roughness coe旺icients.
Equations 13 and 14 can be used to estimate random Within the same tillage operation, substantial variations
roughness of a surface following rainfall from information in hydraulic roughness coefficients were found. Roughness
on cumulative rainfall since the last tillage operation. elements were sometimes larger than water depth. As
Reynolds number increased , variations in flow patterns
INDEPENDENT VARIABLES INFLUENCING HYDRAULIC sometimes occurred. In addition , transition from laminar to
ROUGHNESS COEFFICIENTS turbulent flow conditions may have resulted during a given
Three criteria were established for the model equations test series.
used to predict hydraulic roughness coefficients. The Information from the six tillage treatments was used to
equations should be simple and easily solved using the derive the following regression equation for estimating
fewest number of independent variables necessary to Darcy-Weisbach hydraulic roughness coefficients:
obtain reasonable results. The independent variables should
be generalized , and applicable to conditions beyond those
found in the present study. Finally , the independent 000
variables used in the relationships should be easily F
E -←一-- Moldboard Plow
identified at other locations. D
ι2
0-一一-0 Planter
Variables which could significantly affect hydraulic 』
Qj 100
D
roughness coefficients included random roughness , u
<J)
Reynolds number , slope , type of implement operation , and <J)
Q)
10
hydraulic radius. However , not all of these variable吕would C二
i二
口3
variable. 主 0. 1
>、
ιJ
Including hydraulic radius in the model equation would 」
m
Q
require an iterative procedure to solve the prediction
0.01
equations. An iterative solution would be more time 10 100 1.00e 10.000
consuming and difficult to solve , and could introduce Rey门口Ids Number
difficulties with convergence. Thus , hydraulic radius was Figure l-Darcy·Weisbach roughness coefficients vs. Reynolds
also eliminated from the analysis. number forselected tillage operations.
mmmm
如 V
C。-
O
D 20 40 60 80 0.742
100
Measured Hydraulic Roughness Coefficient - f O.I72RR
n=-一一--一一旦 (16)
Figure LPredtcted vs.measured Darey-weisbach hydraulic Rn 0.282
roughness coefficients.
o一一一。 P l a n te r
86420
。o
0
u-
-o
汇 @O U mm@ E L03O E O F
Em的白 C Z D 3
口00
飞、
1
。
庄0U一 -
币」3
UK工 E D F U -比
K CE币 飞4
30Ln
001 nu
10 100 1.00e 10 ,000 口 2 04 ~6 0.8
Reynolds Number Measured Hydraulic Roughness Coefficient - n
Figure 3- Manning roughness coefficients vs. Reynolds number for Figure 4-Predicted vs. measured Manning hydraulic roughness
目Iected tillage operations. coefficients.
Random roughness is calculated from surface elevation Engman, E.T. 1986. Roughness coefficients for routing
measurements. Information exists in the literature for surface runoff. J. Irrig. and Drain. Engr. , ASCE
relating random roughness values to single and multiple 112(1):39-53.
tillage operations. If cumulative rainfall since the last Finkner, S.C. 1988. Hydraulic roughness coefficients as
tillage operation is known , the reduction in random affected by random roughness. Unpublished M.S. thesis,
roughness caused by precipitation can also be estimated. Univ. ofNebraska, Lincoln.
A field study was conducted to identify random Gilley, J.E., E.R. Kottwitz and J.R. Simanton. 1990.
roughness and corresponding hydraulic roughness Hydraulic characteristics of rills. Transactions ofthe
coefficients over a wide range of conditions. Random ASAE 33(6):1900-1906.
roughness measurements were made following six tillage Issard, c.F. 1944. The surface profile of overland flow.
operations performed on initially smooth soil surfaces. Trans. AGU 4:959-968.
Random roughness measurements were found to be similar Li nden, D.R. and D.M. Van Doren , Jr. 1986. Parameters for
to previously reported values. characterizing tillage-induced soil surface roughness.
Following measurement of random roughness , plot Soil Sci. Soc. Am. J. 50:1560-1565.
surfaces were stabilized using a biodegradable , latex-base Liong, S.Y., S. Selvalingam and D.K. Brady. 1989.
materia l. Steady uniform flow conditions were then Roughness values for overland flow in subcatchments. J,
established for a wide variety of discharge rates. From Irrig. αnd Drain. Engr. , ASCE 1]5(2):203-214.
measurements of discharge rate and flow velocity, Darcy- Luttrell , D.H. 1963. The effect oftillage operations on bulk
Weisbach and Manning hydraulic roughness coefficients density and other physical properties of the soil.
were calculated. Unpublished Ph.D. diss. , Iowa State University , Ames.
Multiple linear regression analysis was used to identify Phelps , H.O. 1975. Shallow laminar flows over rough
the independent variables influencing hydraulic roughness granular surfaces. J. Hydr. Div. , ASCE 10](3):367-384.
coefficients. Hydraulic roughness coefficients were found Podmore , T.H. and L.F. Huggins. ]980. Surface roughne岱
to be significantly affected by random roughness and effects on overland flow. Transactions ofthe ASAE
Reynolds number. The field data were used to obtain 23(6):]43 4-]439 , ]445.
regression relationships which related Darcy-Weisbach and Potter, K.N., T.M. Zobeck and L. J. Hagan. 1990. A
Manning hydraulic roughness coefficients to random microrelief index to estimate soil erodibility by wind.
roughness and Reynolds number. The regression Transactions ofthe ASAE. 33(]):151-155.
relationships can be used for random roughness values Roels , J.M. 1984. Flow resistance in concentrated overland
varying from 6 to 32 mm , and Reynolds numbers ranging flow on rough slope surfaces. Earth Surface Processes
from 20 to 60∞. and Landforms 9:54]-55 1.
Several factors may contribute to total hydraulic Sadeghian , M.R. and J.K. Mitchell. 1990. Hydraulics of
resistance on a given upland site. Information is needed on micro锢braided channels: Resistance to flow on tilled
hydraulic roughness coe旺icients provided by each of these soils. Transactions ofthe ASAE 33(2):458 -468.
factors , their contribution to total hydraulic roughness , and Savat, J. 1980. Resistance to flow in rough supercritical
the effect of flow rate on hydraulic roughness coefficients. sheet flow. Earth SUI内ce Processes and Landforms
The accurate prediction of hydraulic roughness coe旺icients 5:102-122.
will improve our ability to understand and properly model Woo, D.C. and E.F. Brater. 196 1. Laminar flow in rough
upland flow hydraulics. rectangular channels. Journal ofGeophysical Research
66(12):4207-4217.
Yoon, N. ]970. The effect of rainfall on the mechanics of
REFERENCES steady spatially varied sheet flow on a hydraulically
Abrahams , A. D., A .J. Parsons and S.H. Luk. 1986. sm∞th surface. Ph.D. thesis , University of lllinois ,
Resistance to overland flow on desert hillslopes. J矿 Urbana.
Hydrol. 50:343-363. Zobeck , T.M. and C.A. Onstad. 1987. Tillage and rainfall
Allmaras , R.R., R.E. Burwell and R卫Holt. 1967. Plow- e仔ects on random roughness: A review. Soil and Til/age
layer porosity and surface roughness from tillage as Res.9:1-20.
affected by initial porosity and soil moisture at tillage
time. Soil Sci. Soc. Am. Proc. 31:550-556.
Chow, V.T. 1959. Open-Channel Hydraulics. New York: SYMBOLS
McGraw-Hill Book Co., Inc. A = cross-sectional flow area
Currence , D.H. and w' G. Lovely. 1970. The analysis of soil b = flow width
surface roughness. Transactions ofthe ASAE 13(6):710- c = Chezy roughness coe证icient
713. f = Darcy-Weisbach roughness coe征icient
Dunne , T. and w' E. Dietrich. 1980. Experimental study of g = acceleration due to gravity
Horton overland flow on tropical hillslopes. 2. Hydraulic n = Manning roughness coefficient
characteristics and hillslope hydrographs. Z. Geomphol. P = wetted perimeter
Suppl. Band 35:60-80. q = flow rate / unit width
Emmett , W,W. 1970. The hydraulics of overland flow on Q = flow rate
hillslopes. Professional Paper 662-A , U.S. Geological R = hydraulic radius
Survey. Rn = Reynolds number
ω曲 -
wwa IEl
vdn
00'm
v
eeω
-1-
····n
RR = random roughness of a surface following
dm
rainfall
m
oco-1V4
户l
山
6EE-
飞
RRo = random roughness immediately after tillage
v
丁A
RRR = relative random roughness
S = average slope