Resistance in Open Channel Hydraulics

Similar documents
ROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)

Chapter 8 Laminar Flows with Dependence on One Dimension

Chapter 3 Water Flow in Pipes

RESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS

39.1 Gradually Varied Unsteady Flow

CHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

Aerodynamic Admittance Function of Tall Buildings

CEE 3310 Open Channel Flow,, Nov. 18,

12d Model. Civil and Surveying Software. Version 7. Drainage Analysis Module Hydraulics. Owen Thornton BE (Mech), 12d Model Programmer

CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a

5. Secondary Current and Spiral Flow

ERRATA. pp Appendix E.6 Add the following paragraph immediately before the last paragraph on p. 334:

When water (fluid) flows in a pipe, for example from point A to point B, pressure drop will occur due to the energy losses (major and minor losses).

What about water... What about open channel flow... Continuity Equation. HECRAS Basic Principles of Water Surface Profile Computations

Lesson 6: Apparent weight, Radial acceleration (sections 4:9-5.2)

Differential Equations

An Improved Expression for a Classical Type of Explicit Approximation of the Colebrook White Equation with Only One Internal Iteration

Transport Properties: Momentum Transport, Viscosity

12d Model. Civil and Surveying Software. Drainage Analysis Module Hydraulics. Owen Thornton BE (Mech), 12d Model Programmer.

Prof. B.S. Thandaveswara. A short horizontal reach of a prismatic channel is considered. Further, the external

Department of Hydro Sciences, Institute for Urban Water Management. Urban Water

Chapter 8: MULTIPLE CONTINUOUS RANDOM VARIABLES

NEWTONS LAWS OF MOTION AND FRICTIONS STRAIGHT LINES

Hydraulic validation of the LHC cold mass heat exchanger tube.

Controlling the Heat Flux Distribution by Changing the Thickness of Heated Wall

Chapter 11 Collision Theory

Exam 3 Review. F P av A. m V

Water Flow in Open Channels

NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION

4.1 & 4.2 Student Notes Using the First and Second Derivatives. for all x in D, where D is the domain of f. The number f()

Non-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings

Chapter 7: Boundary Layer Theory

FLUID MECHANICS. Lecture 7 Exact solutions

Chapter 7: Boundary Layer Theory

Chapter 14 PROBLEM SOLUTIONS Since vlight v sound, the time required for the flash of light to reach the observer is negligible in

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow

Fluid Physics 8.292J/12.330J

Eighteenth International Water Technology Conference, IWTC18 Sharm ElSheikh, March 2015 REVIEW OF FRICTION FORMULAE IN OPEN CHANNEL FLOW

8.3 Design of Base Plate for Thickness

7. TURBULENCE SPRING 2019

Unified Architecture Level Energy-Efficiency Metric

Watershed Sciences 6900 FLUVIAL HYDRAULICS & ECOHYDRAULICS

HYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.

A REVIEW OF FRICTION FORMULAE IN OPEN CHANNEL FLOW

( x) f = where P and Q are polynomials.

Two-dimensional analytical solution for compound channel flows with vegetated floodplains

NUMERICAL ANALYSES OF COLD-FORMED THIN-WALLED SECTIONS WITH CONSIDERATION OF IMPERFECTIONS DUE TO THE PRODUCTION PROCESS

Numerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder

Analysis of Non-Thermal Equilibrium in Porous Media

Phase Changes Heat must be added or removed to change a substance from one phase to another. Phases and Phase Changes. Evaporation

Buoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center

Submarine sand ripples formation in a viscous fluid: 2D and 3D linear stability analysis

OUTLINE DESIGN OF COLUMN BASE PLATES AND STEEL ANCHORAGE TO CONCRETE 12/21/ Introduction 2. Base plates. 3. Anchor Rods

Syllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.

A note on critical flow section in collector channels

LECTURE NOTES - VIII. Prof. Dr. Atıl BULU

Chapter 10 Flow in Conduits

Friction Factors and Drag Coefficients

Chapter (3) Motion. in One. Dimension

Fin efficiency of the newly developed Compartmented Coil of a Single Coil Twin Fan System

Figure 34: Coordinate system for the flow in open channels.

Incompressible Viscous Flows

COMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION

Increasing and Decreasing Functions and the First Derivative Test. Increasing and Decreasing Functions. Video

LINEAR SUPERPOSITION AND INTERFERENCE PHENOMENA

CRITICAL MASS FLOW RATE THROUGH CAPILLARY TUBES

Comments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium

Heat-fluid Coupling Simulation of Wet Friction Clutch

Kuldeep Rawat*, Ayushman Srivastav* *Assistant Professor, Shivalik College of Engineering, Dehradun.

Available online at ScienceDirect. Energy Procedia 83 (2015 ) Václav Dvo ák a *, Tomáš Vít a

NUMERICAL ANALYSIS OF FORCED CONVECTION HEAT TRANSFER FROM TWO TANDEM CIRCULAR CYLINDERS EMBEDDED IN A POROUS MEDIUM

Module 2. The Science of Surface and Ground Water. Version 2 CE IIT, Kharagpur

3. Several Random Variables

Second Order Slip Flow of Cu-Water Nanofluid Over a Stretching Sheet With Heat Transfer

Filtration. Praktikum Mechanical Engineering. Spring semester ML F16 Tel.:

Sediment transport and river bed evolution

EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018

2015 American Journal of Engineering Research (AJER)

Conservation of Mechanical Energy 8.01

FOCUS ON CONCEPTS Section 7.1 The Impulse Momentum Theorem

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

Chapter 4: Properties of Pure Substances. Pure Substance. Phases of a Pure Substance. Phase-Change Processes of Pure Substances

Uniform Channel Flow Basic Concepts Hydromechanics VVR090

EFFECTS OF HEAT AND MASS TRANSFER FLOW OF A JEFFREY FLUID THROUGH A VERTICAL DEFORMABLE POROUS STRATUM

7.6 Example von Kármán s Laminar Boundary Layer Problem

This is only a list of questions use a separate sheet to work out the problems. 1. (1.2 and 1.4) Use the given graph to answer each question.

Channel Structure Influence on the Thermal-Hydraulic Performance of. Zigzag PCHE

7. Basics of Turbulent Flow Figure 1.

HEADLOSS ESTIMATION. Mekanika Fluida 1 HST

15. Physics of Sediment Transport William Wilcock

Natural convection in a vertical strip immersed in a porous medium

Used to estimate energy loss due to friction in pipe. D = internal diameter of pipe (feet) L = length of pipe (feet) Penn State-Harrisburg

The wind-driven models of Stommel and Munk employed a linearization involving a small parameter, the Rossby number, which we need to reconsider.

Fs (30.0 N)(50.0 m) The magnitude of the force that the shopper exerts is f 48.0 N cos 29.0 cos 29.0 b. The work done by the pushing force F is

Optimum design of high frequency transformer for compact and light weight switch mode power supplies (SMPS)

Tangent Line Approximations

Definition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series.

EFFECTS OF VISCOUS DISSIPATION ON FREE CONVECTION BOUNDARY LAYER FLOW TOWARDS A HORIZONTAL CIRCULAR CYLINDER

Transcription:

16.3.1 Resistance in Open Channel Hdraulics I Manning and Chez equations are compared e1 2 3 1 1 1 2 0 = 2 2 0 1 R S CR S n 21 1-3 2 6 R R C= = n n 1 6 R C= n For laminar low: K R R e1 VR = VR K= 2 2 8gSR 8gVR S 8gR S 8g But K= = = R 2 2 2 V V V C 2 8g C = Re1 K VR 8g I R e1 = C 2 14 R e1 For triangular Smooth Channel (Reer: Chow) 24 R e1 For Rectangular Smooth Channel (Reer: Chow) e1 Sand Roughness Fixed to Flume Bed (Photograph - Thandaeswara)

16.3.2 Laminar Flow with Roughness 60 or a 90 V shape channel. Roughness 0.3023 mm R e1 33 R e1 1.0 0.8 0.6 0.4 Reerence: "Chow Ven Te- Open Channel Hdraulics", Mc Graw Hill Compan, International student edition, 1959, page - 10 0.2 0.1 0.08 0.06 14 R e1 24 R e1 0.04 0.02 0.01 0.008 0.006 Laminar Transitional Turbulent 0.004 10 2 4 6 8 10 2 2 4 6 8 10 3 2 4 6 8 10 4 2 4 6 810 5 2 4 6 810 6 R e1 Variation o riction coeicient with Renolds number Re1( = R ) in smooth channels

Reerence: "Chow Ven Te- Open Channel Hdraulics", Mc Graw Hill Compan, International student edition, 1959, page - 11 1.0 0.1 37.5 cm 0.8 0.08 Varwick 25 cm 0.6 0.06 0.4 0.04 0.2 0.1 0.08 14 R e1 33 Re1 60 Re1 0.02 10 3 2 4 6 8 104 2 4 6 810 5 2 4 6 8 10 6 2 4 6 810 7 Varwick 1 1 20 cm 0.06 0.04 0.02 Laminar Transitional Turbulent 0.01 0.008 0.006 10 2 4 6 810 2 2 4 6 810 3 2 4 6 8 104 2 4 6 810 5 2 4 6 8 10 6 2 4 6 810 7 Re1 Variation o riction coeicient with Renolds number Re1 in rough channels Rectangular Channel - Rough low (Roughness = 0.7188) ( = R ) Bazin conducted experiment using (500 measurements were made at greatest care) (1) Grael embedded in cement. (2) Unpolished wood roughened b transerse wooden strip (i) 27 mm long 10 mm high 10 mm spacing. (ii) 27 mm 10 mm at 50 mm spacing. 3) Cement lining 4) Unpolished wood I the behaior o n and C is to be inestigated then a number o basic deinitions regarding the tpes o hdrodnamic low must be recalled. Flow can be diided into

(i) Hdro dnamicall smooth turbulent low (ii) Hdro dnamicall Rough turbulent low (iii) Hdro dnamicall transition turbulent low. The boundar laer δ or low past a lat plate is gien b 1/ 2 δ Vx o = 5 Laminar x 1/ 5 δ Vx o 7 = 0. 38 turbulent Re 2 10 logarthmic elocit law holds x > V 99% V δ δ δ0 Pseudo boundar Turbulent Transitional region Viscous sub laer Velocit Velocit distribution

Smooth kc = 5 kc is critical roughness height k is roughness height k c = 100 or aerage condition δ k δ0 δ0 δ0 k k kc (a) Smooth kc (b) wa kc Dierent surace roughness (c) rough

Viscous sublaer ks (i) Hdrodnamicall smooth turbulent low (R e ) ks Viscous sublaer (ii) Hdrodnamicall transition low (R e, k s /) ks Viscous sublaer (iii) Hdrodnamicall rough turbulent low (k s /) For hdro dnamicall smooth condition, iscous sub laer submerges the roughness elements. For hdro dnamicall transitional case the roughness element are partl exposed with reerence to iscous sub laer. For hdro dnamicall rough turbulent low the roughness elements are completel exposed aboe the iscous sub laer. For hdro dnamicall rough turbulent low resistance is a unction o Renolds number and the roughness height. K s I we deine R e = shear Renolds number. ; and τo = grs =. ρ

The low is classiied as ollows: K s < 4 Hdrodnamicall smooth K s 4 < < 100 Hdrodnamicall transition K s > 100 Hdrodnamicall ull deeloped turbulent low Summar o Velocit-Proile Equations or Boundar laers with dp 0 dx = Zone Smooth Walls Rough Walls Law o the wall Uniersal equations Laminar sub laer 4 < = ( δ ) - Buer zone 4< < 30 to 70 - - Logarithmi c zone (also called turb ulent laer) > 30 to 70 < 0.15 δ Velocit-deect law Inner region 0.15 δ < (oerlaps with logarithmi c wall law) Outer region 0.15 δ < (approxim ate ormula) (3000 < R e < 70,000) - outer region = A log + B 5.6 log = + 4.9 Power Law = 8.74 k = A log + B k = 5.6 log + B B= (roughness size, shape and distribution) V = A log + B δ V = 5.6 log + 2.5 δ 1 7 V = A log δ V = 8.6 log δ - A and B are constants.

Table shows elocit distributions or dierent conditions Blasius equation or smooth low Smooth pipe low Nikurads e Rough pipe Nikurads e White and Colebroo k ormula VR VR Pipe low equation Re = Open channel low Re = 4 0.3164 5 upto R 0.25 e <10 C=18.755 R mks units or g = 9.806 m/sec e Re 0. 223 1 Re 0.25 =2log R e 251. 5 Re 8g Re > 10 C = 4 2g log 251. C ( ) 1 = 0.86 ln Re - 0.8 1 = 1.14-0.86 ln d 1 /d o 2.51 = 0.86 ln + 37. Re o 1/8 2 C = 17.72 log Re 8g 251. C 3.5294R C = 17.72 log C e R 8g C = 2 log e C 2. 51 8g C = 2 log 8g C = -2 log s 14.83R 12R ks k 2.52 8g + 8g Re Suggested modiication to equation is C ks 2.5 = -2 log + 8g 12R Re [ASCE Task Force Committee 1963]. R is hdraulic mean radius, 4R = Diameter o pipe. In open channel low ollowing aspects come into picture R (, K, C,N, F,U) e (1) (2) (3) In which R e is the Renolds number, K is the Relatie Roughness, C Shape actor o the cross-section, N is the Non- uniormit o the channel both in proile and in plan, F is the Froude number, U is the degree o unsteadiness. In the aboe equation, the irst term corresponds to, Surace Resistance (Friction), the second term corresponds to wae resistance and the third term corresponds to Non uniormit due to acceleration/ deceleration in low.

Surace Resistance: To be accounted based on Karman - Prandtl - elocit distribution. The constant in resistance equation is due to the numerical integration, and is a unction o shape o the cross-section. C 1 R = =A log +B 2g ' For circular section A = 2.0, B = -0.62 For rectangular section: A = 2, B = -0.79 (or large ratio o width/depth) It has remained customar to delineate roughness in terms o the equialent sand grain dimensions k s. For its proper description, howeer, a statistical characteristic such as surace texture requires a series o lengths or length deriaties, though the signiicance o successie terms in the series rapidl approach a minimum. Morris classiied the low into three categories namel (1) isolated roughness low, (2) Wake intererence low, and (3) Quasi smooth low. The igure proides the necessar details. k s s Isolated - roughness low (k/s) - Form drag dominates The wake and the ortex are dissipated beore the next element is reached. The ratio o (k/s) is a signiicant parameter or this tpe o low

k s Wake intererence low (/s) s s When the roughness elements are placed closer, the wake and the ortex at each element will interere with those deeloped b the ollowing element and results in complex orticit and turbulent mixing. The height o the roughness is not important, but the spacing becomes an important parameter. The depth '' controls the ertical extent o the surace region o high leel turbulence. (/s) is an important correlating parameter. k j j j j s s s k is surace roughness height s is the spacing o the elements j is the grooe width is the depth o low Quasi smooth low - k/s or j/s becomes signiicant acts as Pseudo wall Quasi smooth low is also known as skimming low. The roughness elements are so closed placed. The luid that ills in the grooe acts as a pseudo wall and hence low essentiall skims the surace o roughness elements. In such a low (k/s) or (j/s) pla a signiicant role. Concept o three basic tpes o rough surace low k, j, s should describe the characteristics o roughness in one dimensional situations is Areal concentration o or densit distribution o roughness elements. (ater Moris).

16.3.3 Areal concentration or Densit Distribution Roughness Elements Spheres Schlichting, 1936 1 2 Spatial distribution o roughness O'Loughlin and Mcdonald (1964) Cubes arranged as in (1) abd (2) also sand grains (2.5 mm dia)cemented to the bed. Koloseus (1958) and Koloseus and Daidian (1965) conducted experiments using Cubical Roughness Smmetrical diamond shaped pattern.

4 Spheres 3 Cubes Sand ks 2 1 0 0 0.2 0.4 0.6 0.8 1.0 Eectie roughness as a unction o orm pattern, and concentration o roughness elements. (Assuming high Renolds number) 10 Nikuradse λ Areal concentration Schlichting (1936) - Sphere spacing Koloseus (1958) Koloseus and Daidian (1965) Cubical Roughness Smmetrical diamond shaped pattern O'Loughlin and Mcdonald (1964) Cubes arranged as in 1 and in 2. Also sand grains centered to the sand grains (2.5 m diameter) ks 1 Sand 0.1 0.001 0.01 0.1 1.0 λ Areal concentration Logarithmic plot o data rom igure at low concentration Open channel resistance (ater H. Rouse, 1965)

b 1.5 F = 1.5 1.0 V 0.5 1.0 d = 3b 3b 0.5 0 0 0.5 1.0 1.5 2.0 2.5 Froude number, F Resistance o a bridge pier in a wide channel, ater Kobus and Newsham 1.5 S = 5 D 7.5 30 1.0 CD 0.5 D V d D d = 30 S 0 0 0.5 1.0 1.5 2.0 Froude number, F Variation o pier resistance with lateral spacing "S"

0.6 0.4 /b = 1/16 ζ 0.2 /b = 1/8 2b 0.1 /b = 1/4 90 0 4b b 0.2 0.4 0.6 0.8 1.0 2.0 4.0 Froude number, F Loss at one o a series o channel bends ater Haet Some o the important Reerences: (i) Task orce on riction actors in open channels Proc. ASCE JI. o Hd. Dn. Vol. 89., No. H2, March 1963, pp 97-143. (ii) Rouse Hunter, "Critical analsis o open channel resistance", Proceedings o ASCE Journal o Hdraulic diision, Vol.91, Hd 4, pp 1-25, Jul 1965 and discussion pp 247-248, No. 1965, March 1966, pp 387 to 409. Schlichting, "Boundar laer theor", Mc Graw Hill Publication. 16.3.4 Open Channel Resistance There is an optimal area concentration 15% to 25% which produces greater relatie resistance. 1 R A log B DhS + h is the roughness height, S is the areal concentration (<15%), D is the constant which depends on shape and arrangement o the roughness elements. For sanded surace: D = 21 and B = 2.17 The existence o ree surace makes it diicult to assume logarthmic elocit distribution and to integrate oer the entire area o low or dierent cross-sectional shapes. The

lograthmic elocit distribution can be integrated onl or the wide rectangular and circular sections. Eect o boundar non-uniormit is normall ignored and particularl so or graduall aried low proile computation. The dependence on Froude number is clearl seen in case o pier. In case o unstead lows such as loods, it is assumed that the inertial eects are small in comparison with resistance. Hence, the resistance o stead uniorm low at the same depths and elocit is taken to be alid. Where the Froude number exceeds unit, the surace has instabilit in the orm o roll waes. Earlier ormulae or determining C (or details reer to Historical deelopment o Empirical relationships) 1. G.K. Formula (MKS) 2. Bazin s Formula 1897 (MKS) 3. Powell Formula (1950) FPS while using Powell ormula C must be multiplied b 0.5521 to get C in m 1/2 s -1 4. Paloskii Formula (1925) Manning equation is applicable to ull deeloped turbulent rough low. Slope o the straight line is 1:3

n 1/3 ks g g C R 2 C ks k s 1/3 ks 0.113 R I we replace k b diameter o the grain size (d) d 0.113 R s 1/3 8g 8g R C = = 0.113 d or MKS units g = 9.806 m/s 2 8 9.806 R R C = = 26.3482 0.113 d d R or C = 26.34 d R n = C 1 n = d = 0.0379 d 26.34 A number o empirical methods to relate n diameter o the particle are adanced. 1 Strickler (1923) 2 Henderson's interpretation o Strickler's ormula 3a Raudkii (1976) 3b Raudkii (1976) 3c Raudkii (1976) 4 Garde and Ranga Raju n = 0.02789 d [ d in m ] This is not applicable to mobile n = 0.034 d [ d in eet ] 50 n = 0.047 d [ d in m ] bed n = 0.013 d65 [ d in mm ] d 65 = 65 % o the material b n = 0.034 d [ d in eet ] 65 weight smaller. 1/ 6 n = 0.039 d50 [ d in eet ] 0. 039 ( 0. 3048) ( ) 5 Subramana n = 0.0475 d [ d in m ] 6 Meer and Peter and Muller 50 = 0. 039 0. 82036 = 0. 03199 ( 50 ) n = 0.03199 d, d is in 'm' [ ] 90 n = 0.038 d d in m (Signiicant proportion o coarse grained material)

7 Simons and Sentrrk (1976) 8 Lane and Carbon (1953) n = 0.047 d [ d in mm ] n= 0.026 d (d in inches and d 75 = 75% o the 75 material b weight is smaller) 8) Consider = g R S ks 4 < < 100 Transition low R R n = but C = 26.35 C d R d 1 n = = d = 0. 03795 d (d in m) R 2635 26. 35 (. ) Conditon or ull deeloped rough low k s n = 100 d = = 3. 3458 10 n 6 0.03795 n 1 g R S 0.03795 Assuming 6 6 ( ) 6 ( ) -6 2 2 = 1.01 10 m /s g = 9.806 m/s 6 9806. 1 n R S 100 6 6 101. 10 0.03795 n RS 9. 635 10 14 8 6

Blasius equation ( R e <10 5 ) 1 = 2.0 log ( Re ) 2.51 R e C = 4 2g log ( 8g ) 2.51C 180 150 140 130 120 110 100 90 80 70 60 Laminar low Smooth suraces 0.316 0.25 R e 1 8 (C = 15.746 Re, mks) ks = 100 do 2ks = 2R ks 507 252 126 60 30.6 15 10 50 40 Transition zone Full rough zone 1 = C = 2.0 log ( 12R ) or Manning 8g ks Commercial suraces Sand coated surace (Nikuradse) 30 10 3 2 4 6 10 4 10 5 10 6 10 7 10 8 Renolds number R e = 4 V R/ Modiied Mood Diagram showing the Behaior o the Chez C ater Henderson