Practical exercise 7. Surge tank

Similar documents
Applications of Bernoulli s theorem. Lecture - 7

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Flow in porous media

DIRECT CURRENT CIRCUITS

Conservation Law. Chapter Goal. 5.2 Theory

Math 8 Winter 2015 Applications of Integration

Measuring Electron Work Function in Metal

The Wave Equation I. MA 436 Kurt Bryan

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Question 1: Figure 1: Schematic

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.

Motion of Electrons in Electric and Magnetic Fields & Measurement of the Charge to Mass Ratio of Electrons

l 2 p2 n 4n 2, the total surface area of the

KINEMATICS OF RIGID BODIES

AB Calculus Review Sheet

( ) as a fraction. Determine location of the highest

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

( ) where f ( x ) is a. AB Calculus Exam Review Sheet. A. Precalculus Type problems. Find the zeros of f ( x).

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

CBE 291b - Computation And Optimization For Engineers

Math 1B, lecture 4: Error bounds for numerical methods

Sample Problems for the Final of Math 121, Fall, 2005

The Moving Center of Mass of a Leaking Bob

Department of Mechanical Engineering MECE 551 Final examination Winter 2008 April 16, 9:00 11:30. Question Value Mark

Name Class Date. Match each phrase with the correct term or terms. Terms may be used more than once.

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Scientific notation is a way of expressing really big numbers or really small numbers.

Solutions to Physics: Principles with Applications, 5/E, Giancoli Chapter 16 CHAPTER 16

ES.182A Topic 32 Notes Jeremy Orloff

Math 32B Discussion Session Session 7 Notes August 28, 2018

New data structures to reduce data size and search time

Part I: Basic Concepts of Thermodynamics

Objectives. Materials

Simple Harmonic Motion I Sem

#6A&B Magnetic Field Mapping

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

Terminal Velocity and Raindrop Growth

Forces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc

CAPACITORS AND DIELECTRICS

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

The momentum of a body of constant mass m moving with velocity u is, by definition, equal to the product of mass and velocity, that is

First Law of Thermodynamics. Control Mass (Closed System) Conservation of Mass. Conservation of Energy

Thermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report

13: Diffusion in 2 Energy Groups

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Vadose Zone Hydrology

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

13.4 Work done by Constant Forces

7.6 The Use of Definite Integrals in Physics and Engineering

THE IMPORTANCE OF INCLUDING ELASTIC PROPERTY OF PENSTOCK IN THE EVALUATION OF STABILITY OF HYDROPOWEWR PLANTS

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Trigonometric Functions

Lecture 13 - Linking E, ϕ, and ρ

Topic 1 Notes Jeremy Orloff

PHYSICS 211 MIDTERM I 21 April 2004

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Final Exam - Review MATH Spring 2017

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

Section 6: Area, Volume, and Average Value

1 Which of the following summarises the change in wave characteristics on going from infra-red to ultraviolet in the electromagnetic spectrum?

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

Log1 Contest Round 3 Theta Individual. 4 points each 1 What is the sum of the first 5 Fibonacci numbers if the first two are 1, 1?

Math 124A October 04, 2011

Math 0230 Calculus 2 Lectures

7.2 The Definite Integral

This lecture covers Chapter 8 of HMU: Properties of CFLs

Name Solutions to Test 3 November 8, 2017

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

G. MATEESCU 1 A. MATEESCU 1 C. SAMOILĂ 2 D. URSUŢIU 2

Math 113 Exam 1-Review

ME 309 Fluid Mechanics Fall 2006 Solutions to Exam3. (ME309_Fa2006_soln3 Solutions to Exam 3)

Section 7.2 Velocity. Solution

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

The Regulated and Riemann Integrals

Topics Covered AP Calculus AB

( ) Same as above but m = f x = f x - symmetric to y-axis. find where f ( x) Relative: Find where f ( x) x a + lim exists ( lim f exists.

University of Alabama Department of Physics and Astronomy. PH126: Exam 1

Operations with Polynomials

Calculus - Activity 1 Rate of change of a function at a point.

Table of Contents. 1. Limits The Formal Definition of a Limit The Squeeze Theorem Area of a Circle

Summary Information and Formulae MTH109 College Algebra

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

Main topics for the Second Midterm

MCR 3U Exam Review. 1. Determine which of the following equations represent functions. Explain. Include a graph. 2. y x

Algebra Readiness PLACEMENT 1 Fraction Basics 2 Percent Basics 3. Algebra Basics 9. CRS Algebra 1

The Form of Hanging Slinky

Electrical Drive 4 th Class

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

200 points 5 Problems on 4 Pages and 20 Multiple Choice/Short Answer Questions on 5 pages 1 hour, 48 minutes

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

Factors affecting the phonation threshold pressure and frequency

Mathematics Tutorial I: Fundamentals of Calculus

1. a) Describe the principle characteristics and uses of the following types of PV cell: Single Crystal Silicon Poly Crystal Silicon

Practice Final. Name: Problem 1. Show all of your work, label your answers clearly, and do not use a calculator.

Fundamentals of Analytical Chemistry

Problem Set 3 Solutions

Review of basic calculus

Transcription:

Prcticl exercise 7. Surge tnk Introduction Surge tnk is used t hydro power plnts for reduction of wter hmmer tht occurs t closing of the turbine inlet vlve. In this exercise we will mesure mss osciltions t vlve closing on surge tnk model. We will do mesurements t two different vlve closing speeds nd t set flow rte nd hed. igure: culty of civil nd geodetic engineering, bortory for hydrulics. The surge tnk is cylinder mde of plexi glss, seen on the left imge. The right imge shows flow settling tnks with weirs (left nd center of the imge, blue color). Surge tnk Bsic fetures nd vlues of the surge tnk cn be summed in the following lines: - vlue of the surge tnk: t longer hed rce or til rce chnnels/tunnels, the hed rce tunnel cn be mde of thin wlls, which mens cheper, - the surge tnk is n expensive structure, - mss osciltions re present in the surge tnk, - the surge tnk is importnt when mss flow is vried (closing of vlves), - there is no wter hmmer in the hed rce tunnel, or the wter hmmer in the hed rce tunnel is within the rtio of cross sections of the hed rce tunnel nd surge tnk, - wter hmmer is still present in the penstock, - simple surge tnk is the surge tnk where the cross section vries continuously with its height. When the turbine vlve is closed, the wter level in the surge tnk rises due to the inflowing wter from the hed rce tunnel. When wter level in the surge tnk becomes higher thn wter level in the ccumultion lke, the wter strts to flow through the hed rce tunnel bck to the lke or t lest the flow rte in the hed rce tunnel is decresed. After few oscilltions the flow nd the height settle. 31

surge tnk penstock hed rce tunnel lke turbine igure: Surge tnk scheme. Surge tnks re used in power plnts with ccumultion lkes tht re distnt from the turbine (long hed rce tunnels). control vlve for flow pumps: 20 l/s 50 l/s 100 l/s slow settling tnk 1 (weir) flow settling tnk 2 (weir) mnul setting of height surge tnk vlve ngle meter closing vlve tnk in the bsement 60 m 3 pressure trnsducer pressure tp tringulr notch weir computer for dt cquisition view from the front side, dischrge to the bsement tnk igure: bortory experiment, digrm. low settling tnks serve for ccurte setting of hed, becuse wter level is exctly set by the weir. The second flow settling tnk enbles mnul setting of height, which llows simultion of lke height bove the surge tnk height. Assumptions, necessry for derivtion of equtions When deriving equtions for mss oscilltions in the surge tnk we use the following ssumptions: - cross section S S(z), mening tht the surge tnk cross section chnges continuously with respect to height, - incompressibility of fluid or wter, = const, - in ny instnt the flow rte Q in the hed rce tunnel does not vry long its length, - inerti of wter mss in the surge tnk nd penstock is neglected (smll mss compred to the mss in the hed rce tunnel), - kinetic energy of wter in the surge tnk is neglected, becuse the speed of wter level oscilltions in the surge tnk is smll nd we ssume hydrosttic pressure gh t surge tnk entry, - wter level in the lke is not chnging. We will write the dynmic nd continuity equtions for surge tnk. Derivtion nd nomenclture re prtilly tken from the book Hidrvlik nestlneg tok, Rjr R., Univerz v jubljni, kultet z grdbeništvo in geodezijo, jubljn, 1980. 32

z level S=S(z) cross-section hj lke -w,-q w,q -z g b Qt h hb hdin hstt turbine river igure: Derivtion of equtions for the surge tnk, symbols. Derivtion of dynmic eqution We proceed from the Euler eqution. Chnge of velocity in the hed rce tunnel is consequence of the pressure difference tht ccelertes the fluid flow (2nd term on the right) nd inerti due to mss forces m (1st term on the right, index m mens force per unit mss). Mss forces m consist of grvittionl force nd friction force. dw 1 p m dt x The sign of the pressure term (2nd term on the right) is negtive becuse flow ccelertes from the loction of smller pressure (loction ) to the loction of higher pressure (loction b). m g, m fr, m Sin() is written s rtio of sides (rtio of heights/), becuse for smll ngles sin(). Slope of hed rce tunnel is usully smll. g, m g sin m g h hb riction force is proportionl to drg, tht is squre of velocity or flow rte. E fr mg fr E fr fr 2 Q Q Q mge fr mgq Q Absolute vlue in the bove eqution serves for preservtion of friction force sign. riction force in the bove eqution hs positive sign, becuse the eqution for mss forces hs friction written with negtive sign. The friction force per unit mss is written with the following expression. 33

fr, m fr m mgq Q m gq Q We express prtil derivtive of pressure with respect to x s difference of pressures t the beginning nd the end of the hed rce tunnel, the pressure difference is further written with heights p p x p b p g h z h h h j b j g z h We insert ll the bove written terms into the Euler eqution. w h Q Q hb z h hb g g g t w Q Q z g g t We replce the chnge in velocity with the chnge in flow rte. Index hrt mens hed rce tunnel. 1 S hrt 1 Shrt gs hrt Q Q Q g g t z / g h Q Q Q z g g t Q z Q Q 0 dynmic eqution for surge tnk t b Derivtion of the continuity eqution We write the continuity eqution on the bsis of the following considertion: flow rte to the turbine equls flow rte through the hed rce tunnel minus flow rte to the surge tnk (if the surge tnk is emptying, the flow rte is negtive, which mens tht the flow rte to the turbine is higher thn the flow rte through the hed rce tunnel). Qs is flow rte through the surge tnk, Qhrt is flow rte through the hed rce tunnel nd Qt is flow rte through the turbine. Ss is the surge tnk cross-section nd Shrt is the tunnel cross-section. Symbol v mens velocity in the surge tnk nd w mens velocity in the tunnel. Q Q t Q S v S t dz v dt s hrt Q s hrt w Agin we write the flow rte Q in the hed rce tunnel without the index, the sme s in the dynmic eqution. We copy the dynmic eqution. 34

dz Qt Ss Q continuity eqution for surge tnk dt gshrt Q z Q Q 0 dynmic eqution for surge tnk t We vry the flow rte to the turbine by using vlve, therefore the flow rte to the turbine is known vrible. The unknown vribles re flow rte in the hed rce tunnel nd the wter level in the surge tnk. Q=Q(t) z=z(t) flow rte in the hed rce tunnel wter level in the surge tnk We hve two equtions with two unknowns if we don't tke into ccount the friction coefficient nd geometricl dimensions. Initil nd boundry conditions for solution of the system of equtions ) initil condition 2 z( t 0) Q 0 b) left side boundry condition, the lke wter level is not chnging hj=const c) right side boundry condition, we ssume closing of the vlve ccording to the following eqution 1 t Q t t Q0 T0 time of vlve closing t the turbine T0 Mesuring of flow rte with tringulr notch weir The tringulr notch (V-notch) weir is intended for mesuring of fluid flow rte in open chnnels. The wter height t the notch is used to determine the flow rte. The figure below shows two tringulr notch weirs with ngles of =90 nd =45. low rte through the weir is sum of flow through infinitesimlly thin surfce res (htched re) cross the whole wter height H. Since the wter height for ech surfce is different, the outflow velocity is lso different. We ssume n element of height h being on height h. We replce the tringulr weir with two right tringles. We write the element width b with the following eqution: b 2 2 H htn The surfce re A of the htched element on the figure equls the product of height nd width. 35

A 2 2 H htn h igure: Tringulr notch weir intended for mesuring of fluid flow rte in open chnnels. Above is the weir with =90 ngle, below is the weir with =45 ngle. We write the ouflow velocity in the sme wy s the outflow velocity for fluids in tnk. The outflow velocity increses proportionl to the squre root of the wter height bove the outflow surfce. v 2gh low rte through the infinitesimlly thin element of the notch is product of its surfce re nd flow velocity. Q 2 2 2 H htn gh h If we integrte the bove expression between h=0 nd h=h, we get the following expression. Q 2 tn 2 8 15 tn 2 2g H 0 2g H Hh 5/ 2 1/ 2 h 3/ 2 dh 2 tn 2 2 2g H 3 5/ 2 2 5 H 5/ 2 The theoreticlly derived flow rte is not exctly the sme s the mesured one, therefore we introduce dischrge coefficient Cd, which mkes the ctul flow rte to be 8 Q Cd tn 2g H 15 2 5/ 2. 36

The dvntge of the tringulr notch over the rectngulr notch weir is tht the shpe of the nppe (the body of wter dischrging over the notch of the weir) is not chnging significntly. Tht mens tht the dischrge coefficient Cd is not chnging significntly for different flow rtes Q. The tringulr notch weir llows us to mesure wide rnge of flow rtes. Use the eqution bove when determining the operting point for this exercise. Description of mesurement loctions nd mesurement equipment Mesurements include the following vribles: - mesurement of wter height in the flow settling tnk, - mesurement of wter flow rte in the supply pipe, - mesurement of wter flow rte through the surge tnk t the tringulr notch weir, - mesurement of wter level in the surge tnk by using pressure trnsducer, - mesurement of vlve ngle, - determintion of integrl vribles or system dimensions, - recording of mesured vlues on computer. Wter height in the flow settling tnk is mesured by using tpe meter. Height in the second settling tnk is tken s reference. Wter flow rte in the supply pipe is mesured by the ABB Wtermster electromgnetic flowmeter. The flowmeter is wired to the DAQ crd nd the computer, which enbles disply nd recording of mesured vlues. Wter flow rte through the surge tnk is mesured by using 45 tringulr notch weir locted downstrem of the surge tnk nd vlve. It is not possible to mesure nd record the wter flow rte electroniclly. Use equtions tht were derived bove. Determine the dischrge coefficient Cd from digrm tht is locted in the lbortory. Wter level in the surge tnk is mesured by using the Endress Huser pressure trnsducer with rnge between 0 nd 4 br. Mesurement with this trnsducer is performed reltive to the tmospheric pressure. Tht mens tht only one pressure port hs to be connected, the other one is open to the tmosphere. Air hs to be removed from the pressure trnsducer nd connecting pipes before the mesurement. The pressure trnsducer is wired to the DAQ crd nd the computer. The voltge signl from the trnsducer hs to be clibrted on the computer to disply the ctul pressure. Since reltion between pressure nd wter height in the tnk is liner, the signl cn be directly clibrted to the wter level height in the surge tnk. Note: Electricl connection of the trnsducer is done by two wires. Supply voltge for the trnsducer is DC 24 V. or two-wire connection it is necessry to use voltge divider in order to cquire the signl by DAQ crd. or this purpose pproprite resistors re connected ccording to the figure below. The DAQ crd rnge is from 0 to 5 V. 37

pressure trnsducer power supply 24 V igure: Electricl connection of the pressure trnsducer to voltmeter or DAQ crd. Mesurement of vlve ngle will be performed by using precise rottionl resistor. Electricl resistnce of the rottionl resistor is dependent on the ngle of shft. The rottionl resistor hs three terminls for connection to the power supply nd to the DAQ crd. The rnge of the DAQ crd is from 0 to 5 V. Determine which ngle corresponds to individul voltges on DAQ crd nd mke clibrtion curve for the rottionl resistor. Integrl vribles re: cross sections of the supply pipe nd surge tnk nd time of vlve closing. Before or fter the mesurements you hve to determine the system dmping in stedy stte. Mesure the flow rte by using the tringulr notch weir. Assume tht the vlve closing is liner. Recording of flow rte, pressure nd vlve ngle to computer will be mde by using 16 bit DAQ crd NI 6036 nd NI bview softwre. Choose suitble time nd frequency of dt cquisition for recording to the hrddrive. The softwre will be prepred by the ssistnt. Assignment Assemble mesuring system tht will llow mesuring of the wter height in the surge tnk, the flow rte in the supply pipe nd the vlve ngle. Clibrte ll signls for the electronic meters by clculting equtions of line. Mesure flow rte nd wter height oscilltions in the surge tnk model. Mke mesurements t two different vlve closing speeds. In stedy stte compre flow rte indicted by the flowmeter to the flow rte mesured by the tringulr notch weir. 38