Block 3 Open channel flow

Size: px
Start display at page:

Download "Block 3 Open channel flow"

Transcription

1 Numerical Hydraulics Block 3 Open channel flow Markus Holzner

2 Contents of the course Block 1 The equations Block Computation of pressure surges Block 3 Open channel flow (flow in rivers) Block 4 Numerical solution of open channel flow Block 5 Transport of solutes in rivers Block 6 Heat transport in rivers

3 Saint Venant equations in 1D b(h) l b A(h) continuity momentum equation z h 3

4 Saint Venant equations in 1D continuity (for section without inflow) Q A + = 0 x t Momentum equation from integration of Navier-Stokes/Reynolds equations over the channel cross-section: v v h P τ l ρ + ρα v = ρg t x x A 0 b 4

5 Saint Venant equations in 1D Friction force: Only the shear stress at the channel walls survives integration. Per unit volume it is f friction force τ0lbδx τ0 = = = volume AΔx R For the wall shear we can insert hy λ τ ρ λ λ 8 ( ) 0 = v = (Re, k/ R hy ) 5

6 Saint Venant equations in 1D The friction can be expressed as energy loss per flow distance: ΔE ΔE/ V τ 0 = freibung Δx = V Δx R Using friction slope and channel slope hy I E Δ τ = = = = Δ E/ V 0 λ 1 v z IS xρg ρgrhy 4 Rhy g x Alternative: Strickler/Manning equation 6

7 Saint Venant equations in 1D we finally obtain v v h + v = g( I ) S IE g t x x ( va( h)) A( h) + = 0 x t For a rectangular channel : A = bh h h v + v + h = 0 t x x 7

8 Approximations and solutions Steady state solution Kinematic wave Diffusive wave Full equations 8

9 Steady state solution (rectangular channel) dv v = g( I I ) g S E dx dh dv v + h = 0 dx dx dh dx 9

10 Steady state solution (rectangular channel) Uniform flow (no advective acceleration): ( I I ) = 0 S Full solution (insert second equation into first): E dh IS IE v = = 0 with Fr = dx 1 Fr gh yields water surface profiles 10

11 Classification of profiles h gr = h cr = water depth at critical flow h N = water depth at normal flow I 0 = slope of channel bottom I gr = I cr = critical slope h gr = h cr I gr = I cr Horizontal channel bottom I 0 = 0 H : h > h gr H 3 : h < h gr 11

12 Classification of profiles Mild slope: h N > h gr I 0 < I gr M1: h N <h > h gr M: h N > h > h gr M3: h N > h < h gr h gr = h cr I gr = I cr Steep slope: h N < h gr I 0 > I gr S1: h N <h > h gr S: h N < h < h gr S3: h N > h < h gr 1

13 Classification of profiles Critical slope h N = h gr I 0 = I gr C1: h N < h C3: h N > h h gr = h cr I gr = I cr Negative slope I 0 < 0 N: h > h gr N3: h < h gr 13

14 Examples for composite slopes and transitions h gr = h cr I gr = I cr 14

15 Numerical solution (explicit FD method) Subcritical flow: Computation in upstream direction hx ( +Δx) hx ( ) IS IE( x+δx) = = Δx 1 Fr( x+δx) 0 Solve for h(x) Supercritical flow: Computation in downstream direction hx ( +Δx) hx ( ) Δx IS IE( x) = = 1 Fr( x) 0 Solve for h(x+ x) 15

16 Approximations Kinematic wave Diffusive wave Dynamic wave (Full equations) v v h = v g + g IS I t x x h h v + v + h = 0 t x x ( ) E Rectangular channel 16

17 Approximations In the different approximations different terms in the equation of motion are neglected against the term gi S : 1. Approximation: Kinematic wave g I ( I ) = 0. Approximation: Diffusive wave h g + g IS IE = x Complete solution: Dynamic wave S E ( ) 0 v v h + v + g g IS IE = t x x ( ) 0 17

18 Kinematic wave Normal flow depth. Energy slope is equal to channel bottom slope. Therefore Q is only a function of water depth. E.g. using the Strickler/Manning equation: Q dq h = x dh x Q= Ak R I = Q( h) /3 1/ str hy S Inserting into the continuity equation yields h h + c ' = 0 t x This is the form of a wave equation (see pressure surge) with wave velocity c = v+c dq dh c ' = b Instead of using Q=Q(h) the equation can be derived using v=v(h) 18

19 Kinematic wave With the Strickler/Manning equation we get: /3 hb Q = hbkstr I h+ b 1/ S and /3 1/3 1/ 1/ str S str S dq dh bh bh b c' = = k I + k hi b b+ h 3 b+ h ( b+ h) For a broad channel the hydraulic radius is approximately equal to the water depth. The wave velocity then becomes c' = c = v + c 3 v = 5 3 k str I 1/ S h /3 = 5 3 v = v + 3 v 19

20 Kinematic wave The wave velocity is not constant as v is a function of water depth h. Varying velocities for different water depth lead to self-sharpening of wave front Pressure propagates faster than the average flow. Advantage of approximation: PDE of first order, only one upstream boundary condition required. Disadvantage of approximation: Not applicable for bottom slope 0. No backwater feasible as there is no downstream boundary condition. 0

21 Diffusive wave Now Q is not only a function of h but also of Insertion into the continuity equation yields: h x = I S I E with I E = I E (Q/A) from Strickler or Darcy-Weisbach h / x h 1 Q + = 0 t b x h 1 Q( h, h/ x) h 1 Q( h, h/ x) h + + = t b h x b ( h/ x) x 0 1

22 Diffusive wave This equation has the form of an advectiondiffusion equation with a wave velocity w and a diffusion coefficient D: with 0 ' = + x h D x h c t h b h Q c = / ' b x h Q D ) / ( / =

23 Diffusive wave Using the Strickler/Manning equation and assuming a broad rectangular channel (h = R hy ) one obtains: and 5/3 Q= vhb= kstrbh IS h/ x Insertion into the continuity equation yields with = bv x 3 x I h/ x x Q 5 h vbh h h t c'= + c 5 3 v h x ' D and ( ) S h x D = = 0 vh / ( I h x) S 3

24 Diffusive wave D is always positive, as the energy slope I = I h/ x E is always positive in flow direction. S The wave moves downstream and flattens out diffusively. A lower boundary condition is necessary because of the second derivative. This allows the implementation of a backwater effect. 4

25 St. Venant equation as wave equation v v h + v = g + g I I S t x x h h v + v + h = t x x 0 ( ) E Linear combinations: Multiply second equation with ± l and add to first equation 5

26 St. Venant equation as wave equation v v v h h h + v + λh + λ + λv + g g IS IE = t x x t x x Write derivatives of h and v as total derivatives along a characteristic line: ( ) 0 v v h g h ( v λh) λ ( v ) g( I ) S IE = 0 t x t λ x Choosing λ =± g h the two characteristics have the same relative wave velocity (with respect to average water velocity). 6

27 St. Venant equation as wave equation v v g h h + ( v± gh) ± + ( v± gh) g( IS IE) = 0 t x gh t x The characteristics are therefore: dx = c' = v ± gh = v ± c dt and the relative wave velocity for shallow water waves is c= gh Difference to the surge in pipes: v cannot be neglected in comparison to c 7

28 General form of the St.-Venant equation General cross section: Width at water surface b(h), cross sectional area A(h) and da(h)/dh = b(h) v v h + v = g + g( I ) S IE t x x h h v bh ( ) + vbh ( ) + Ah ( ) = 0 t x x Rel. wave velocity: c= ga()/() h b h 8

29 Solution of the full St. Venant equation (Form with b = const. for clarity) v v h = v g + g I I S t x x h h v + v + h = t x x 0 ( ) E The solution is found numerically. E.g. with a difference method. In the following an explicit difference method is discussed first. 9

Numerical Hydraulics

Numerical Hydraulics ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow

More information

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

Department of Hydro Sciences, Institute for Urban Water Management. Urban Water Department of Hydro Sciences, Institute for Urban Water Management Urban Water 1 Global water aspects Introduction to urban water management 3 Basics for systems description 4 Water transport 5 Matter

More information

CIE4491 Lecture. Hydraulic design

CIE4491 Lecture. Hydraulic design CIE4491 Lecture. Hydraulic design Marie-claire ten Veldhuis 19-9-013 Delft University of Technology Challenge the future Hydraulic design of urban stormwater systems Focus on sewer pipes Pressurized and

More information

OPEN CHANNEL FLOW. One-dimensional - neglect vertical and lateral variations in velocity. In other words, Q v = (1) A. Figure 1. One-dimensional Flow

OPEN CHANNEL FLOW. One-dimensional - neglect vertical and lateral variations in velocity. In other words, Q v = (1) A. Figure 1. One-dimensional Flow OPEN CHANNEL FLOW Page 1 OPEN CHANNEL FLOW Open Channel Flow (OCF) is flow with one boundary exposed to atmospheric pressure. The flow is not pressurized and occurs because of gravity. Flow Classification

More information

UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW

UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW Derivation of uniform flow equation Dimensional analysis Computation of normal depth UNIFORM FLOW 1. Uniform flow is the flow condition obtained from a

More information

CEE 3310 Open Channel Flow, Nov. 26,

CEE 3310 Open Channel Flow, Nov. 26, CEE 3310 Open Channel Flow, Nov. 6, 018 175 8.10 Review Open Channel Flow Gravity friction balance. y Uniform Flow x = 0 z = S 0L = h f y Rapidly Varied Flow x 1 y Gradually Varied Flow x 1 In general

More information

Hydraulics for Urban Storm Drainage

Hydraulics for Urban Storm Drainage Urban Hydraulics Hydraulics for Urban Storm Drainage Learning objectives: understanding of basic concepts of fluid flow and how to analyze conduit flows, free surface flows. to analyze, hydrostatic pressure

More information

Chapter 4: Non uniform flow in open channels

Chapter 4: Non uniform flow in open channels Chapter 4: Non uniform flow in open channels Learning outcomes By the end of this lesson, students should be able to: Relate the concept of specific energy and momentum equations in the effect of change

More information

We will assume straight channels with simple geometries (prismatic channels) and steady state flow (in time).

We will assume straight channels with simple geometries (prismatic channels) and steady state flow (in time). 56 Review Drag & Lift Laminar vs Turbulent Boundary Layer Turbulent boundary layers stay attached to bodies longer Narrower wake! Lower pressure drag! 8. Open-Channel Flow Pipe/duct flow closed, full,

More information

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING Urban Drainage: Hydraulics. Solutions to problem sheet 2: Flows in open channels

DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING Urban Drainage: Hydraulics. Solutions to problem sheet 2: Flows in open channels DEPRTMENT OF CIVIL ND ENVIRONMENTL ENGINEERING Urban Drainage: Hydraulics Solutions to problem sheet 2: Flows in open channels 1. rectangular channel of 1 m width carries water at a rate 0.1 m 3 /s. Plot

More information

CIVL4120/7020 Advanced open channel hydraulics and design - Tutorial (1) Unsteady open channel flows

CIVL4120/7020 Advanced open channel hydraulics and design - Tutorial (1) Unsteady open channel flows School of Civil Engineering at the University of Queensland CIVL4120/7020 Advanced open channel hydraulics and design - Tutorial (1) Unsteady open channel flows Attendance to tutorials is very strongly

More information

Hydraulics Part: Open Channel Flow

Hydraulics Part: Open Channel Flow Hydraulics Part: Open Channel Flow Tutorial solutions -by Dr. K.N. Dulal Uniform flow 1. Show that discharge through a channel with steady flow is given by where A 1 and A 2 are the sectional areas of

More information

The University cannot take responsibility for any misprints or errors in the presented formulas. Please use them carefully and wisely.

The University cannot take responsibility for any misprints or errors in the presented formulas. Please use them carefully and wisely. Aide Mémoire Suject: Useful formulas for flow in rivers and channels The University cannot take responsiility for any misprints or errors in the presented formulas. Please use them carefully and wisely.

More information

HEC-RAS v5.0: 2-D applications

HEC-RAS v5.0: 2-D applications HEC-RAS v5.0: 2-D applications Tom Molls, Will Sicke, Holly Canada, Mike Konieczki, Ric McCallan David Ford Consulting Engineers, Inc. Sacramento, CA September 10, 2015: Palm Springs FMA conference What

More information

CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - II GRADUALLY VARIED FLOW

CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - II GRADUALLY VARIED FLOW CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - II GRADUALLY VARIED FLOW Dynamic equations of gradually varied and spatially varied flows - Water surface flow profile classifications: Hydraulic Slope, Hydraulic

More information

Dr. Muhammad Ali Shamim ; Internal 652

Dr. Muhammad Ali Shamim ; Internal 652 Dr. Muhammad Ali Shamim ali.shamim@uettaxila.edu.pk 051-904765; Internal 65 Channel Tranistions A channel transition is defined as change in channel cross section e.g. change in channel width and/or channel

More information

Gradually Varied Flow I+II. Hydromechanics VVR090

Gradually Varied Flow I+II. Hydromechanics VVR090 Gradually Varied Flow I+II Hydromechanics VVR090 Gradually Varied Flow Depth of flow varies with longitudinal distance. Occurs upstream and downstream control sections. Governing equation: dy dx So Sf

More information

Open Channel Flow Part 2. Ch 10 Young, notes, handouts

Open Channel Flow Part 2. Ch 10 Young, notes, handouts Open Channel Flow Part 2 Ch 10 Young, notes, handouts Uniform Channel Flow Many situations have a good approximation d(v,y,q)/dx=0 Uniform flow Look at extended Bernoulli equation Friction slope exactly

More information

Numerical Modelling and Hydraulics

Numerical Modelling and Hydraulics Numerical Modelling and Hydraulics Nils Reidar B. Olsen Department of Hydraulic and Environmental Engineering The Norwegian University of Science and Technology 3rd edition, 9. March 2012. ISBN 82-7598-074-7

More information

P = 2Rθ. The previous Manning formulas are used to predict V o and Q for uniform flow when the above expressions are substituted for A, P, and R h.

P = 2Rθ. The previous Manning formulas are used to predict V o and Q for uniform flow when the above expressions are substituted for A, P, and R h. Uniform Flow in a Partly Full, Circular Pipe Fig. 10.6 shows a partly full, circular pipe with uniform flow. Since frictional resistance increases with wetted perimeter, but volume flow rate increases

More information

VARIED FLOW IN OPEN CHANNELS

VARIED FLOW IN OPEN CHANNELS Chapter 15 Open Channels vs. Closed Conduits VARIED FLOW IN OPEN CHANNELS Fluid Mechanics, Spring Term 2011 In a closed conduit there can be a pressure gradient that drives the flow. An open channel has

More information

Uniform Channel Flow Basic Concepts Hydromechanics VVR090

Uniform Channel Flow Basic Concepts Hydromechanics VVR090 Uniform Channel Flow Basic Concepts Hydromechanics VVR090 ppt by Magnus Larson; revised by Rolf L Feb 2014 SYNOPSIS 1. Definition of Uniform Flow 2. Momentum Equation for Uniform Flow 3. Resistance equations

More information

Open-channel hydraulics

Open-channel hydraulics Open-channel hydraulics STEADY FLOW IN OPEN CHANNELS constant discharge, other geometric and flow characteristics depended only on position Uniform flow Non-uniform flow S; y; v const. i i 0 i E y 1 y

More information

Applicability of the local inertial approximation of the shallow water equations to flood modeling

Applicability of the local inertial approximation of the shallow water equations to flood modeling WATER RESOURCES RESEARCH, VOL. 49, 4833 4844, doi:10.1002/wrcr.20366, 2013 Applicability of the local inertial approximation of the shallow water equations to flood modeling Gustavo A. M. de Almeida 1

More information

39.1 Gradually Varied Unsteady Flow

39.1 Gradually Varied Unsteady Flow 39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in

More information

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

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004 OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie

More information

SOE2156: Fluids Lecture 7

SOE2156: Fluids Lecture 7 Weirs and SOE2156: Fluids Lecture 7 Vee Vee Last lecture { assumed the channel was uniform (constant depth, shape, slope etc.) { steady uniform Found that : location of free surface to be determined 2

More information

conservation of linear momentum 1+8Fr = 1+ Sufficiently short that energy loss due to channel friction is negligible h L = 0 Bernoulli s equation.

conservation of linear momentum 1+8Fr = 1+ Sufficiently short that energy loss due to channel friction is negligible h L = 0 Bernoulli s equation. 174 Review Flow through a contraction Critical and choked flows The hydraulic jump conservation of linear momentum y y 1 = 1+ 1+8Fr 1 8.1 Rapidly Varied Flows Weirs 8.1.1 Broad-Crested Weir Consider the

More information

Hydromechanics: Course Summary

Hydromechanics: Course Summary Hydromechanics: Course Summary Hydromechanics VVR090 Material Included; French: Chapters to 9 and 4 + Sample problems Vennard & Street: Chapters 8 + 3, and (part of it) Roberson & Crowe: Chapter Collection

More information

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow Uniform Channel Flow Basic Concepts Hydromechanics VVR090 Uniform occurs when: Definition of Uniform Flow 1. The depth, flow area, and velocity at every cross section is constant 2. The energy grade line,

More information

Q = α n AR2/3 h S1/2 0. bd 2d + b d if b d. 0 = ft1/3 /s

Q = α n AR2/3 h S1/2 0. bd 2d + b d if b d. 0 = ft1/3 /s CEE 330 Open Channel Flow, Nov., 00 7 8.8 Review Open Channel Flow Gravity friction balance. In general we take an energy equation approach. y Uniform Flow x = 0 z = S 0L = h f where we could find h f

More information

Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation

Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation Numerical Hydraulics Autumn semester 2016 Prof. Dr. Markus Holzner Author: Pascal

More information

EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018

EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 Q1. Using Cheng s formula estimate the settling velocity of a sand particle of diameter 1 mm in: (a) air; (b) water. Q2. Find the critical Shields parameter diameter

More information

Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow

Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 3 4, pp. 89 109 DOI: 10.1515/heem-2015-0006 IBW PAN, ISSN 1231 3726 Comparison of Average Energy Slope Estimation Formulas

More information

Numerical Hydraulics

Numerical Hydraulics ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction

More information

FE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)

FE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering) Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.

More information

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis. OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

Presented by: Civil Engineering Academy

Presented by: Civil Engineering Academy Presented by: Civil Engineering Academy Open-Channel Flow Uniform Flow (See CERM Ch. 19) Characterized by constant depth volume, and cross section. It can be steady or unsteady Non-uniform Flow *Not on

More information

Guo, James C.Y. (1999). "Critical Flow Section in a Collector Channel," ASCE J. of Hydraulic Engineering, Vol 125, No. 4, April.

Guo, James C.Y. (1999). Critical Flow Section in a Collector Channel, ASCE J. of Hydraulic Engineering, Vol 125, No. 4, April. Guo, James C.Y. (1999). "Critical Flow Section in a Collector Channel," ASCE J. of Hydraulic Engineering, Vol 15, No. 4, April. CRITICAL FLOW SECTION IN A COLLECTOR CHANNEL By James C.Y. Guo, PhD, P.E.

More information

Basic Hydraulics. Rabi H. Mohtar ABE 325

Basic Hydraulics. Rabi H. Mohtar ABE 325 Basic Hydraulics Rabi H. Mohtar ABE 35 The river continues on its way to the sea, broken the wheel of the mill or not. Khalil Gibran The forces on moving body of fluid mass are:. Inertial due to mass (ρ

More information

Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum

Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz

More information

1.060 Engineering Mechanics II Spring Problem Set 4

1.060 Engineering Mechanics II Spring Problem Set 4 1.060 Engineering Mechanics II Spring 2006 Due on Monday, March 20th Problem Set 4 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members

More information

Basic Fluid Mechanics

Basic Fluid Mechanics Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible

More information

3.2 CRITICAL DEPTH IN NONRECTANGULAR CHANNELS AND OCCUR- RENCE OF CRITICAL DEPTH

3.2 CRITICAL DEPTH IN NONRECTANGULAR CHANNELS AND OCCUR- RENCE OF CRITICAL DEPTH 3.2 CRITICAL DEPTH IN NONRECTANGULAR CHANNELS AND OCCUR- RENCE OF CRITICAL DEPTH Critical Depth in Non-Rectangular Channels Consider an irregular channel: da w dd dd d Specific energy is defined as: E

More information

The Shallow Water Equations

The Shallow Water Equations The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow

More information

Open Channel Flow - General. Hydromechanics VVR090

Open Channel Flow - General. Hydromechanics VVR090 Open Channel Flow - General Hydromechanics VVR090 ppt by Magnus Larson; revised by Rolf L Jan 2014, Feb 2015 SYNOPSIS 1. Introduction and Applications 2. The History of Open Channel Flow 3. Flow Classification

More information

Water Flow in Open Channels

Water Flow in Open Channels The Islamic Universit of Gaza Facult of Engineering Civil Engineering Department Hdraulics - ECIV 33 Chapter 6 Water Flow in Open Channels Introduction An open channel is a duct in which the liquid flows

More information

MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS

MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS M.L. Kavvas and L.Liang UCD J.Amorocho Hydraulics Laboratory University of California, Davis, CA 95616, USA Uncertainties

More information

Pressure in stationary and moving fluid Lab- Lab On- On Chip: Lecture 2

Pressure in stationary and moving fluid Lab- Lab On- On Chip: Lecture 2 Pressure in stationary and moving fluid Lab-On-Chip: Lecture Lecture plan what is pressure e and how it s distributed in static fluid water pressure in engineering problems buoyancy y and archimedes law;

More information

Open Channel Flow - General. Open Channel Flow

Open Channel Flow - General. Open Channel Flow Open Channel Flow - General Hydromechanics VVR090 Open Channel Flow Open channel: a conduit for flow which has a free surface Free surface: interface between two fluids of different density Characteristics

More information

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6 Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture

More information

NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS

NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS TASK QUARTERLY 15 No 3 4, 271 282 NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS MICHAŁ SZYDŁOWSKI Faculty of Civil and Environmental Engineering, Gdansk University of Technology, Narutowicza

More information

Compressible Duct Flow with Friction

Compressible Duct Flow with Friction Compressible Duct Flow with Friction We treat only the effect of friction, neglecting area change and heat transfer. The basic assumptions are 1. Steady one-dimensional adiabatic flow 2. Perfect gas with

More information

P10.5 Water flows down a rectangular channel that is 4 ft wide and 3 ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow.

P10.5 Water flows down a rectangular channel that is 4 ft wide and 3 ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow. P10.5 Water flows down a rectangular channel that is 4 ft wide and ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow. Solution: Convert the flow rate from 15,000 gal/min

More information

Flood routing. Prof. (Dr.) Rajib Kumar Bhattacharjya Indian Institute of Technology Guwahati

Flood routing. Prof. (Dr.) Rajib Kumar Bhattacharjya Indian Institute of Technology Guwahati Flood routing Prof. (Dr.) Rajib Kumar Bhattacharjya Indian Institute of Technology Guwahati Guwahati, Assam Email: rkbc@iitg.ernet.in Web: www.iitg.ernet.in/rkbc Visiting Faculty NIT Meghalaya Q (m 3 /sec)

More information

y 2 = 1 + y 1 This is known as the broad-crested weir which is characterized by:

y 2 = 1 + y 1 This is known as the broad-crested weir which is characterized by: CEE 10 Open Channel Flow, Dec. 1, 010 18 8.16 Review Flow through a contraction Critical and choked flows The hydraulic jump conservation of linear momentum y = 1 + y 1 1 + 8Fr 1 8.17 Rapidly Varied Flows

More information

Empirical versus Direct sensitivity computation

Empirical versus Direct sensitivity computation Empirical versus Direct sensitivity computation Application to the shallow water equations C. Delenne P. Finaud Guyot V. Guinot B. Cappelaere HydroSciences Montpellier UMR 5569 CNRS IRD UM1 UM2 Simhydro

More information

Modeling and simulation of bedload transport with viscous effects

Modeling and simulation of bedload transport with viscous effects Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR

More information

Fluid Mechanics II Viscosity and shear stresses

Fluid Mechanics II Viscosity and shear stresses Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small

More information

Linear Analysis of Coupled Equations for Sediment Transport

Linear Analysis of Coupled Equations for Sediment Transport Theme B of the XXVII IAHR Congress, San Francisco, 1-15 August, 1997, 156-161. Linear Analysis of Coupled Equations for Sediment Transport YANTAO CUI and GARY PARKER St. Anthony Falls Laboratory, University

More information

Beaver Creek Corridor Design and Analysis. By: Alex Previte

Beaver Creek Corridor Design and Analysis. By: Alex Previte Beaver Creek Corridor Design and Analysis By: Alex Previte Overview Introduction Key concepts Model Development Design Accuracy Conclusion Refresh v = Beaver Creek Site = Wittenberg Introduction Low head

More information

1.060 Engineering Mechanics II Spring Problem Set 8

1.060 Engineering Mechanics II Spring Problem Set 8 1.060 Engineering Mechanics II Spring 2006 Due on Monday, May 1st Problem Set 8 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members

More information

Chapter 8: Flow in Pipes

Chapter 8: Flow in Pipes 8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump

More information

OPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK

OPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK OPEN CHANNEL FLOW Numerical Methods and Computer Applications Roland Jeppson TECHNISCHE INFORMATlONSBiBUOTHEK UNIVERSITATSB'BUOTHEK HANNOVER Si. i. CRC Press Taylor &.Francis Group Boca Raton London New

More information

Lecture Notes - Sediment Transport The flow problem Overview

Lecture Notes - Sediment Transport The flow problem Overview Lecture Notes - Sediment Transport The flow problem Overview In the last lecture, we emerged with a transport model in which the rate of transport q s depends on the shear stress τ. Recall that, for the

More information

Sediment transport and river bed evolution

Sediment transport and river bed evolution 1 Chapter 1 Sediment transport and river bed evolution 1.1 What is the sediment transport? What is the river bed evolution? System of the interaction between flow and river beds Rivers transport a variety

More information

Lecture 9: Tidal Rectification, Stratification and Mixing

Lecture 9: Tidal Rectification, Stratification and Mixing Lecture 9: Tidal Rectification, Stratification and Mixing Chris Garrett 1 Additional Notes on Tidal Rectification This lecture continues the discussion of long-wavelength tidal flow over comparatively

More information

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.

Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis. OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric

More information

NOTES ON OPEN CHANNEL FLOW

NOTES ON OPEN CHANNEL FLOW NOTES ON OPEN CANNEL FLOW Prof. Marco Pilotti Facoltà di Ingegneria, Università degli Studi di Brescia Profili di moto permanente in un canale e in una serie di due canali - Boudine, 86 OPEN CANNEL FLOW:

More information

L. Pratt and J. Whitehead 6/25/ Parabolic Bottom

L. Pratt and J. Whitehead 6/25/ Parabolic Bottom 2.8 Parabolic Bottom Up to this point we have dealt strictly with channels with rectangular crosssections. The only allowable variation of bottom elevation has been in the longitudinal (y) direction. Although

More information

EFFECT OF SPATIAL AND TEMPORAL DISCRETIZATIONS ON THE SIMULATIONS USING CONSTANT-PARAMETER AND VARIABLE-PARAMETER MUSKINGUM METHODS

EFFECT OF SPATIAL AND TEMPORAL DISCRETIZATIONS ON THE SIMULATIONS USING CONSTANT-PARAMETER AND VARIABLE-PARAMETER MUSKINGUM METHODS INDIAN INSTITUTE OF TECHNOLOGY ROORKEE EFFECT OF SPATIAL AND TEMPORAL DISCRETIZATIONS ON THE SIMULATIONS USING CONSTANT-PARAMETER AND VARIABLE-PARAMETER MUSKINGUM METHODS Muthiah Perumal and C. Madhusudana

More information

Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws

Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws G. Bastin http://perso.uclouvain.be/georges.bastin/ Georges Bastin Jean Michel Coron Progress in Nonlinear

More information

NPTEL Quiz Hydraulics

NPTEL Quiz Hydraulics Introduction NPTEL Quiz Hydraulics 1. An ideal fluid is a. One which obeys Newton s law of viscosity b. Frictionless and incompressible c. Very viscous d. Frictionless and compressible 2. The unit of kinematic

More information

Lecture 12: Transcritical flow over an obstacle

Lecture 12: Transcritical flow over an obstacle Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem

More information

!y + v!d.!y = 0. (1.2.4)

!y + v!d.!y = 0. (1.2.4) L. Pratt J. Whitehead 5/23/06 1.2 The shallow water equations one-dimensional wave propagation. Traditional discussions of hydraulic effects such as those found in engineering text books are often based

More information

!y + v!d.!y = 0. (1.2.4)

!y + v!d.!y = 0. (1.2.4) L. Pratt J. Whitehead 9/13/05 1.2 The shallow water equations one-dimensional wave propagation. Traditional discussions of hydraulic effects such as those found in engineering text books are often based

More information

V/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0

V/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0 UNIT III FLOW THROUGH PIPES 1. List the types of fluid flow. Steady and unsteady flow Uniform and non-uniform flow Laminar and Turbulent flow Compressible and incompressible flow Rotational and ir-rotational

More information

Pressure in stationary and moving fluid. Lab-On-Chip: Lecture 2

Pressure in stationary and moving fluid. Lab-On-Chip: Lecture 2 Pressure in stationary and moving fluid Lab-On-Chip: Lecture Fluid Statics No shearing stress.no relative movement between adjacent fluid particles, i.e. static or moving as a single block Pressure at

More information

REE 307 Fluid Mechanics II. Lecture 1. Sep 27, Dr./ Ahmed Mohamed Nagib Elmekawy. Zewail City for Science and Technology

REE 307 Fluid Mechanics II. Lecture 1. Sep 27, Dr./ Ahmed Mohamed Nagib Elmekawy. Zewail City for Science and Technology REE 307 Fluid Mechanics II Lecture 1 Sep 27, 2017 Dr./ Ahmed Mohamed Nagib Elmekawy Zewail City for Science and Technology Course Materials drahmednagib.com 2 COURSE OUTLINE Fundamental of Flow in pipes

More information

Comparing HEC-RAS v5.0 2-D Results with Verification Datasets

Comparing HEC-RAS v5.0 2-D Results with Verification Datasets Comparing HEC-RAS v5.0 2-D Results with Verification Datasets Tom Molls 1, Gary Brunner 2, & Alejandro Sanchez 2 1. David Ford Consulting Engineers, Inc., Sacramento, CA 2. USACE Hydrologic Engineering

More information

Bernoulli and Pipe Flow

Bernoulli and Pipe Flow Civil Engineering Hydraulics Mechanics of Fluids Head Loss Calculations Bernoulli and The Bernoulli equation that we worked with was a bit simplistic in the way it looked at a fluid system All real systems

More information

FLOW IN CONDUITS. Shear stress distribution across a pipe section. Chapter 10

FLOW IN CONDUITS. Shear stress distribution across a pipe section. Chapter 10 Chapter 10 Shear stress distribution across a pipe section FLOW IN CONDUITS For steady, uniform flow, the momentum balance in s for the fluid cylinder yields Fluid Mechanics, Spring Term 2010 Velocity

More information

PIPE FLOWS: LECTURE /04/2017. Yesterday, for the example problem Δp = f(v, ρ, μ, L, D) We came up with the non dimensional relation

PIPE FLOWS: LECTURE /04/2017. Yesterday, for the example problem Δp = f(v, ρ, μ, L, D) We came up with the non dimensional relation /04/07 ECTURE 4 PIPE FOWS: Yesterday, for the example problem Δp = f(v, ρ, μ,, ) We came up with the non dimensional relation f (, ) 3 V or, p f(, ) You can plot π versus π with π 3 as a parameter. Or,

More information

FORMULA SHEET. General formulas:

FORMULA SHEET. General formulas: FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to

More information

FORMATION OF ROLL WAVES IN LAMINAR SHEET FLOW. by Pierre Y. Julien and David M. Hartley. January 1985

FORMATION OF ROLL WAVES IN LAMINAR SHEET FLOW. by Pierre Y. Julien and David M. Hartley. January 1985 FORMATION OF ROLL WAVES IN LAMINAR SHEET FLOW by Pierre Y. Julien and David M. Hartley January 1985 Civil Engineering Department Engineering Research Center Colorado State University Fort Collins, Colorado

More information

ScienceDirect. A pipe network simulation model with dynamic transition between free surface and pressurized flow

ScienceDirect. A pipe network simulation model with dynamic transition between free surface and pressurized flow Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 70 ( 2014 ) 641 650 12th International Conference on Computing and Control for the Water Industry, CCWI2013 A pipe network simulation

More information

NON UNIFORM FLOW IN OPEN CHANNEL

NON UNIFORM FLOW IN OPEN CHANNEL For updated version, please click on http://ocw.ump.edu.my HYDRAULICS NON - UNIFORM FLOW IN OPEN CHANNEL TOPIC 3.3 by Nadiatul Adilah Ahmad Abdul Ghani Faculty of Civil Engineering and Earth Resources

More information

CONVECTIVE HEAT TRANSFER

CONVECTIVE HEAT TRANSFER CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar

More information

1. Open Channel Hydraulics

1. Open Channel Hydraulics Open Channel Flow. Open Channel Hydraulics.... Definition and differences between pipe flow and open channel flow.... Types of flow.... Properties of open channels...4.4 Fundamental equations... 5.4. The

More information

WATER INJECTION DREDGING by L.C. van Rijn

WATER INJECTION DREDGING by L.C. van Rijn WATER INJECTION DREDGING by L.C. van Rijn (info@leovanrijn-sediment.com) Description of method Almost all harbour basins suffer from the problem of siltation of sediments. Usually, the deposited materials

More information

THE EFFECTS OF OBSTACLES ON SURFACE LEVELS AND BOUNDARY RESISTANCE IN OPEN CHANNELS

THE EFFECTS OF OBSTACLES ON SURFACE LEVELS AND BOUNDARY RESISTANCE IN OPEN CHANNELS Manuscript submitted to 0th IAHR Congress, Thessaloniki, 4-9 August 00 THE EFFECTS OF OBSTACLES ON SURFACE LEVELS AND BOUNDARY RESISTANCE IN OPEN CHANNELS J. D. FENTON Department of Civil and Environmental

More information

Q1 Give answers to all of the following questions (5 marks each):

Q1 Give answers to all of the following questions (5 marks each): FLUID MECHANICS First Year Exam Solutions 03 Q Give answers to all of the following questions (5 marks each): (a) A cylinder of m in diameter is made with material of relative density 0.5. It is moored

More information

Thelongwaveequations

Thelongwaveequations Alternative Hydraulics Paper 1, January 2010 Thelongwaveequations Institute of Hydraulic and Water Resources Engineering Vienna University of Technology, Karlsplatz 13/222, 1040 Vienna, Austria URL: http://johndfenton.com/

More information

Deriving one dimensional shallow water equations from mass and momentum balance laws

Deriving one dimensional shallow water equations from mass and momentum balance laws Deriving one dimensional shallow water equations from mass and momentum balance laws Amir Dehghanghadikolaei 1, Navid Namdari 2, Behrouz Mohammadian 3, Seyed Reza Ghoreishi 4 1Mechanical, Industrial and

More information

IMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2

IMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2 IMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2 1 Associated Professor, Irrigation and Hydraulic Department, College of Technical Engineering,

More information

Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations

Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations 1 Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations LAAS - CNRS Toulouse (France) Workshop on Irrigation Channels and Related Problems

More information

7. Basics of Turbulent Flow Figure 1.

7. Basics of Turbulent Flow Figure 1. 1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds

More information

A note on critical flow section in collector channels

A note on critical flow section in collector channels Sādhan ā, Vol. 26, Part 5, October 2001, pp. 439 445. Printed in India A note on critical flow section in collector channels 1. Introduction SUBHASISH DEY Department of Civil Engineering, Indian Institute

More information

CRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN

CRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN Criteria for the choice of flood routing methods in natural... CRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN NATURAL CHANNELS WITH OVERBANK FLOWS Roger Moussa 1 Abstract: The classification of river

More information