Determination of Flow Resistance Coefficients Due to Shrubs and Woody Vegetation

Similar documents
Determination of Flow Resistance Coefficients Due to Shrubs and Woody Vegetation

"HIP Modeling Methodology Based on the Inherent Process Anisotropy

Social Studies 201 Notes for March 18, 2005

Relationship between surface velocity divergence and gas transfer in open-channel flows with submerged simulated vegetation

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis

Social Studies 201 Notes for November 14, 2003

Rupro, breach model used by Cemagref during Impact project

CONFIDENCE LIMITS ON PREDICTED VALUES

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE

MODELLING OF FRICTIONAL SOIL DAMPING IN FINITE ELEMENT ANALYSIS

Comparing Means: t-tests for Two Independent Samples

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES

ρ water = 1000 kg/m 3 = 1.94 slugs/ft 3 γ water = 9810 N/m 3 = 62.4 lbs/ft 3

MAE 101A. Homework 3 Solutions 2/5/2018

Subject: Useful formulas for sediment transport and morphology in rivers and channels

Interaction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Unified Correlation between SPT-N and Shear Wave Velocity for all Soil Types

Cake ltration analysis the eect of the relationship between the pore liquid pressure and the cake compressive stress

Nearshore Sediment Transport Modeling: Collaborative Studies with the U. S. Naval Research Laboratory

The Combined Effect of Wind and Rain on Interrill Erosion Processes

Seismic Loads Based on IBC 2015/ASCE 7-10

Estimating floor acceleration in nonlinear multi-story moment-resisting frames

CHAPTER 3 LITERATURE REVIEW ON LIQUEFACTION ANALYSIS OF GROUND REINFORCEMENT SYSTEM

Constant Force: Projectile Motion

APPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS

Study of a Freely Falling Ellipse with a Variety of Aspect Ratios and Initial Angles

THE RATIO OF DISPLACEMENT AMPLIFICATION FACTOR TO FORCE REDUCTION FACTOR

Calculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient

THEORY FOR HOPPER SEDIMENTATION.

STRAIN LIMITS FOR PLASTIC HINGE REGIONS OF CONCRETE REINFORCED COLUMNS

Bernoulli s equation may be developed as a special form of the momentum or energy equation.

A Comparison of Correlations for Heat Transfer from Inclined Pipes

Transitional behaviors in well-graded coarse granular soils. Associate professor, State Key Laboratory of Coal Mine Disaster Dynamics and Control,

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011

18.6 COMPTON SCATTERING OF PHOTONS FROM ELECTRONS IN MAGNETICALLY INSULATED TRANSMISSION LINES* Sandia Laboratories, Albuquerque, New Mexico 87185

PART I: AN EXPERIMENTAL STUDY INTO THE VISCOUS DAMPING RESPONSE OF PILE-CLAY INTERFACES

Bahram Noshad Department of Electrical Engineering, Bandar Deylam Branch, Islamic Azad University, Bandar Deylam, Iran.

Unified Design Method for Flexure and Debonding in FRP Retrofitted RC Beams

UNITS FOR THERMOMECHANICS

Soil water electrical conductivity determination based on the salinity index concept

Experimental Determination of Particle Sedimentation Velocity in Opaque Drilling Fluids

Research Article Reliability of Foundation Pile Based on Settlement and a Parameter Sensitivity Analysis

PARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES

Optimal Coordination of Samples in Business Surveys

Interaction Diagram - Tied Reinforced Concrete Column (Using CSA A )

Optimization model in Input output analysis and computable general. equilibrium by using multiple criteria non-linear programming.

Dynamic Behaviour of Timber Footbridges

Standard Guide for Conducting Ruggedness Tests 1

ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME

Sediment Transport in Shallow Overland Flow

Streaming Calculations using the Point-Kernel Code RANKERN

Inferences Based on Two Samples: Confidence Intervals and Tests of Hypothesis Chapter 7

External Flow: Flow over Bluff Objects (Cylinders, Spheres, Packed Beds) and Impinging Jets

PhysicsAndMathsTutor.com

1. Basic introduction to electromagnetic field. wave properties and particulate properties.

Experimental study of the heat transfer for a tube bundle in a transversally flowing air

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter

THEORETICAL CONSIDERATIONS AT CYLINDRICAL DRAWING AND FLANGING OUTSIDE OF EDGE ON THE DEFORMATION STATES

Fate of 2-Chloro Ethyl Ethyl Sulfide on 13X Molecular Sieve Adsorbent Implications for Regenerative Filtration

Midterm Review - Part 1

Alternate Dispersion Measures in Replicated Factorial Experiments

Solved problems 4 th exercise

CONFERENCE PROCEEDINGS VOLUME I

Finding the location of switched capacitor banks in distribution systems based on wavelet transform

Analysis of cavitating flow through a venturi

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Mucus Transport in the Human Lung Airways: Effect of Porosity Parameter and Air Velocity

BUBBLES RISING IN AN INCLINED TWO-DIMENSIONAL TUBE AND JETS FALLING ALONG A WALL

A PLC BASED MIMO PID CONTROLLER FOR MULTIVARIABLE INDUSTRIAL PROCESSES

pure VIF on bluff bodie, almot all the exiting parameter identification procee have baed on the diplacement repone of VIV, intead of the force hitory

Effect of river bed slope and particle size distribution on washing out condition of trees in rivers

Effects of vector attenuation on AVO of offshore reflections

Convective Heat Transfer

ASSESSING EXPECTED ACCURACY OF PROBE VEHICLE TRAVEL TIME REPORTS

Modeling of Transport and Reaction in a Catalytic Bed Using a Catalyst Particle Model.

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions

Christian Linde Olsen Griffith University, Faculty of Engineering and Information Technology, Gold Coast Campus.

DEVELOPMENT OF A STRUCTURED THERMOCLINE THERMAL ENERGY STORAGE SYSTEM

Computation of Velocity, Pressure and Temperature Profiles in a Cryogenic Turboexpander

Frames of Reference and Relative Velocity

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Critical Height of Slopes in Homogeneous Soil: the Variational Solution

SIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS

Module: 8 Lecture: 1

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

Horizontal Biaxial Loading Tests on Sliding Lead Rubber Bearing System

SEDIMENT TRANSPORT AT THE BOTTOM OF SEA WAVES. G. Vittori 1, P. Blondeaux 1

Extending MFM Function Ontology for Representing Separation and Conversion in Process Plant

Green-Kubo formulas with symmetrized correlation functions for quantum systems in steady states: the shear viscosity of a fluid in a steady shear flow

Analysis of Mesostructure Unit Cells Comprised of Octet-truss Structures

Clustering Methods without Given Number of Clusters

SIMULATING THE STRESS AND STRAIN BEHAVIOR OF LOESS VIA SCC MODEL

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject

Estimation of Peaked Densities Over the Interval [0,1] Using Two-Sided Power Distribution: Application to Lottery Experiments

Transcription:

December 000 Determination of Flow Reitance Coefficient Due to hrub and Woody Vegetation by Ronald R. Copeland PURPOE: The purpoe of thi Technical Note i to tranmit reult of an experimental invetigation into the effect of vegetation (particularly ground cover plant, mall tree, and hrub) on flow reitance. INTRODUCTION: An important conideration for determining the tage-dicharge relationhip in river and tream i the effect or influence of vegetation on the overall head lo along a channel and in the overbank. Plant in the floodplain and along the bank can increae or even decreae the effective flow reitance. The vegetation may be natural or it may have been planted to improve aethetic or habitat, to prevent eroion, or for other reaon. ydraulic loe and drag due to actual plant were meaured at the Utah tate Univerity Water Reearch Laboratory utilizing a large wide flume and a maller ectional flume. Reearch in the flume reulted in the collection of data from more than 0 experiment with 0 different plant pecie. Experiment were conducted with both homogeneou and mixed plant grouping. ingle-tem and multiple-tem plant were included in the plant type evaluated. Plant with and without leave were evaluated. Plant denity, pacing, and ize were varied in the experiment. Plant were evaluated over a range of velocitie and depth. A methodology wa developed from the laboratory data to predict head lo and reitance coefficient a a function of lope and depth. Input data can be collected from the field or etimated plant characteritic may be ued in the methodology. The evaluation of vegetative impact on propoed and exiting channel to determine flow capacity and water-urface elevation require proper hydraulic roughne value for hrub and other aethetically and environmentally deirable plant. Given the near complete lack of hydraulic roughne value for hrub and imilar vegetation, the accurate etimation of channel capacity and water-urface elevation ha previouly been difficult at bet. Detail of the tudy may be found in Freeman, Rahmeyer, and Copeland (in preparation). REITANCE COEFFICIENT: Reitance to flow i typically characterized by a roughne coefficient. The mot commonly ued equation for flow reitance i the Manning equation. The ratio of hear velocity to mean velocity, V * /V, i another form of reitance coefficient, and may be thought of a the ratio of hear tre to inertial force. All variable are defined in Appendix I. There are other reitance coefficient in ue including the Darcy-Weibach friction factor, f, and the Chezy C. Thee can all be converted eaily to Manning n. In thi tudy, reitance equation were developed for the hear velocity to average velocity ratio becaue it i dimenionle and ha a ound theoretical bai, and for the Manning coefficient, becaue it ue i widepread. The Manning reitance coefficient for vegetation i calculated in 1

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information i etimated to average 1 hour per repone, including the time for reviewing intruction, earching exiting data ource, gathering and maintaining the data needed, and completing and reviewing the collection of information. end comment regarding thi burden etimate or any other apect of thi collection of information, including uggetion for reducing thi burden, to Wahington eadquarter ervice, Directorate for Information Operation and Report, 115 Jefferon Davi ighway, uite 104, Arlington VA 0-430. Repondent hould be aware that notwithtanding any other proviion of law, no peron hall be ubject to a penalty for failing to comply with a collection of information if it doe not diplay a currently valid OMB control number. 1. REPORT DATE DEC 000. REPORT TYPE 3. DATE COVERED 00-00-000 to 00-00-000 4. TITLE AND UBTITLE Determination of Flow Reitance Coefficient Due to hrub and Woody Vegetation 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTOR() 5d. PROJECT NUMBER 5e. TAK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME() AND ADDRE(E) Army Engineer Reearch and Development Center,Coatal and ydraulic Laboratory,3909 all Ferry Road,Vickburg,M,39180 8. PERFORMING ORGANIZATION REPORT NUMBER 9. PONORING/MONITORING AGENCY NAME() AND ADDRE(E) 10. PONOR/MONITOR ACRONYM() 1. DITRIBUTION/AVAILABILITY TATEMENT Approved for public releae; ditribution unlimited 13. UPPLEMENTARY NOTE 14. ABTRACT 11. PONOR/MONITOR REPORT NUMBER() 15. UBJECT TERM 16. ECURITY CLAIFICATION OF: 17. LIMITATION OF ABTRACT a. REPORT unclaified b. ABTRACT unclaified c. TI PAGE unclaified ame a Report (AR) 18. NUMBER OF PAGE 8 19a. NAME OF REPONIBLE PERON tandard Form 98 (Rev. 8-98) Precribed by ANI td Z39-18

December 000 conformity with the Cowan (1956) method for additive reitance. Thi method conit of addition to roughne for variou urface irregularitie and vegetation. TIFFNE MODULU: The modulu of plant tiffne, E, i critical to the calculation of reitance becaue of the flexibility of the plant and the deformation of leaf mae due to the flow force. The modulu of plant tiffne i calculated by E F 45 F45 = = 6.791 4 3 I D (1) The data neceary to ue Equation 1 i obtained by meauring the force, F 45, neceary to bend the plant to an angle of 45 deg. The 45-deg angle i meaured from the initial vertical poition to the tem or leaf ma at the point where the force i meaured i.e., at / a hown in Figure 1. / / D F 45 Figure 1. Methodology for meauring plant tiffne for calculating E in the field, for plant with effective height of leaf ma approximately equal to the plant height The reearch performed in the laboratory and in the field indicated that the tiffne modulu can be etimated from the relationhip of E to the ratio of /D (Equation ). Thi equation give the modulu in pound per quare foot, while Equation 3 give the value in newton per quare meter. E lbf ft = 1.597E05 + 454 37. 8 D D D 3 () E ( N ) = 7.648E06 +.174E04 + 1.809E03 m D D D 3 (3)

December 000 Actual field meaurement of E are recommended where poible. ince the tiffne modulu varie depending on the plant ize, it wa determined that if the calculated modulu for a particular plant ize wa divided by (/D ) 1.5, the tiffne modulu became independent of plant ize and one value could be ued for all plant ize. Meaured tiffne moduli for plant ued in the experimental tudy are reported in Freeman, Rahmeyer, and Copeland (in preparation). REITANCE EQUATION FOR UBMERGED VEGETATION: Reult from large flume experiment were analyzed to determine the regreion of variable for ubmerged vegetation. The analyi found that log and polynomial relationhip gave a poor data fit while a power relationhip had very good reult. The parameter in the equation were modified to allow a direct olution for reitance (for a given depth) by combining the original parameter with Manning equation and the equation for hear velocity. Thi modification and combination of equation reulted in Equation 4 for hear velocity and Equation 5 for Manning n. In thee equation the reitance coefficient repreent the combined reitance of the bed and the plant. Reitance coefficient due only to vegetation mut be determined by ubtracting the bed reitance. In thee experiment the Manning bed reitance coefficient wa found to be 0.0 and V * /V for the bed wa found to be 0.069. V* V = g C E A = 0.183 ρa iv* 0.183 Y O 0.43 ( MA ) i 0.73 v V* R h 0.115 (4) n = K n E A 0.183 ρav i * 0.183 Y O 0.43 v V* R 0.115 1 V * 0.73 ( MA ) ( R ) / 3 ( ) 1/ i h h (5) It i important to note that the plant characteritic, A i, and A are the initial characteritic of the plant without the effect of flow ditortion. During the experiment, it wa oberved that ince the plant bent with flow, ubmergence occurred at flow depth le than 80 percent of the plant height. Equation 4 and 5 are to be applied only for ubmerged flow defined by Y o > 0.8. REITANCE EQUATION FOR PARTIALLY UBMERGED VEGETATION: Data for partially ubmerged vegetation were analyzed to determine the regreion of variable. The regreion analyi again found that a log relationhip gave a poor fit of data while a power relationhip produced very good reult. Equation 6 and 7 fit the data well and allow direct olution for reitance if the flow depth i known. ere again, in thee equation the reitance coefficient repreent the combined reitance of the bed and the plant. Reitance coefficient due only to vegetation mut be determined by ubtracting the bed reitance. In thee experiment the Manning bed reitance coefficient wa found to be 0.0 and V * /V for the bed wa found to be 0.069. V* V = g C = 0.150 0.6 E A * 0.166 V* Rh 3.487E 05 ( ) * ρ MAi Ai V* υ (6) 3

December 000 n = K n E A 3.487E 05 * ρ Ai V * 0.150 0.6 ( ) 3 * 0.166 V* R R / 1/ M A h h i ν V* (7) The blockage area in Equation 6 and 7 wa changed to an effective area, A i *, ince only a portion of the leaf ma produce blockage under partially ubmerged flow condition. CONCLUION: When plant were ubmerged, it wa oberved that the plant leaf ma tended to trail downtream forming a treamlined, almot teardrop-haped profile. The leaf ma hape changed with velocity and became more treamlined with increaing velocity. The effect of thi phenomenon wa a ignificant decreae in the drag coefficient and reitance coefficient with velocity. On the other hand, reitance increaed with depth for partially ubmerged plant a the blockage area increaed with depth until the plant were ubmerged. The tranition between ubmerged and partially ubmerged flow occurred at a depth of about 80 percent of the undeflected plant height. It wa alo oberved that the leaf ma or foliage canopy diverted flow beneath the canopy. The bottom flow reulted in ignificant velocitie along the channel bed cauing general cour and increaed ediment tranport. The bed velocitie were ufficient to tranport and move the larget ize of gravel found in the flume bed. The hydraulic roughne of a vegetated channel wa hown to be a function of the tiffne of the plant growing in the channel, the depth, velocity, and hydraulic radiu of the channel, plant denity, and frontal area of the plant obtructing the flow. It wa determined that the roughne can be calculated directly if the depth of flow i known. The modulu of plant tiffne, E, i critical to the calculation of reitance becaue of the flexibility of the plant and the deformation of leaf mae due to the flow force. The reearch performed in the laboratory and in the field indicated that the tiffne modulu can be etimated from the relationhip of E to the ratio of /D. Actual field meaurement of E are recommended where poible. The tiffne modulu can alo be etimated from meaured value of imilar plant. ince the tiffne modulu varie depending on the plant ize, it wa determined that if the calculated modulu for a particular plant ize wa divided by (/D ) 1.5, the tiffne modulu became independent of plant ize and one value could be ued for all plant ize. ADDITIONAL INFORMATION: For additional information contact Dr. Ronald R. Copeland, Coatal and ydraulic Laboratory, U.. Army Engineer Reearch and Development Center, 3909 all Ferry Road, Vickburg, M 39180, at 601-634-63 or e-mail Ronald.R.Copeland@erdc.uace.army.mil. The content of thi TN are not to be ued for advertiing, publication, or promotional purpoe. Citation of trade name doe not contitute an official endorement or approval of the ue of uch commercial product. 4

December 000 REFERENCE Cowan, W. (1956). Etimating hydraulic roughne coefficient, Agricultural Engineering, 37(7), 473-475. Freeman, G. E., Rahmeyer, W.., and Copeland, R. R. (in preparation). Determination of reitance due to hrub and woody vegetation, Technical Report, U.. Army Engineer Reearch and Development Center, Vickburg, M. 5

December 000 Appendix I: Lit of Variable Numerou variable are ued throughout Freeman, Rahmeyer, and Copeland (in preparation). Thoe variable and their definition are preented here. Figure I1 and I how the meaurement that define the variable involving the leaf ma and plant dimenion for ubmerged and emergent (unubmerged or partially ubmerged) flow condition. A Cro ectional flow area, ft or m² A i Frontal area of an individual plant blocking flow, approximated by the equivalent rectangular area of blockage by W e, ft or m² * A i A Net ubmerged frontal area of a partially ubmerged plant, ft or m² Total cro-ectional area of all of the tem() of an individual plant, meaured at /4, ft² or m² b Width of channel flume, ft or m C Chezy reitance coefficient, ft 1/ / or m 1/ / C D Drag coefficient of vegetation, dimenionle D tem diameter, meaured at a height of /4, ft or m E Exponential cientific notation E Modulu of plant tiffne, lbf/ft² or N/m² f Darcy-Weibach friction factor, dimenionle f b Friction factor for the bed and plant, dimenionle f w Friction factor for the wall, dimenionle The horizontal force neceary to bend a plant tem 45 deg, lbf or N F 45 F D F r g ' Drag force, lbf or N Froude number, dimenionle Acceleration due to gravity = 3.17 ft/ or 9.806 m/² Average undeflected plant height, ft or m Undeflected height of the leaf ma of a plant, ft or m * Undeflected height of leaf ma that i below water urface for a partially ubmerged plant, ft or m (ee Figure I) I econd moment of inertia of cro ection of plant tem, ft 4 or m 4 K n Unit converion factor for Manning equation, 1.4861 ft 1/3 /ec or 1.0 m 1/3 /ec L Channel reach length, ft or m M Relative plant denity, number of plant per ft² or m² n Total Manning roughne coefficient, including idewall roughne 6

n b Manning reitance coefficient for vegetation and channel bed n veg Manning reitance coefficient for vegetation n o Manning reitance coefficient for the bed P Wetted perimeter, ft or m R e Reynold number, R e = V R h /ν R h ydraulic radiu, R h = flow area / wetted perimeter, ft or m R b ydraulic radiu for the bed and plant, ft or m R w ydraulic radiu for the wall, ft or m Bed or energy lope, dimenionle o Bed lope, dimenionle f Energy lope, dimenionle V Mean channel velocity, ft/ or m/ V P Local plant approach velocity in front of the leaf ma, ft/ or m/ V * hear velocity, V * = (g R h ) ½, ft/ or m/ V * /V Reitance coefficient, dimenionle Y o Flow depth, ft or m W e Equivalent average plant width, We = A i /, ft or m dy/dx Unit change in lope of the water urface γ pecific weight of water, lbf/ft 3 or N/m 3 ν Fluid dynamic vicoity, ft²/ or m²/ ρ Fluid denity, lug/ft 3 ( lbf-ec /ft 4 ) or kg/m 3 τ o hear tre on channel bottom, (τ o = γ R h ), lbf/ft or N/m² ERDC/CL CETN-VIII-3 December 000 7

December 000 W e Y O D Figure I1. Plant dimenion definition for ubmerged plant Figure I. Plant dimenion definition for partially ubmerged plant 8