Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models

Similar documents
Pressure Head: Pressure head is the height of a column of water that would exert a unit pressure equal to the pressure of the water.

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

Abstract. 1 Introduction

Lecture Note for Open Channel Hydraulics

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

CIE4491 Lecture. Hydraulic design

Hydraulics for Urban Storm Drainage

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

Design Data 17. Partial Flow Conditions of Box Culverts

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

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

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

D. MATHEMATICAL MODEL AND SIMULATION

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

Numerical Hydraulics

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

Application of Mathematical Modeling to Study Flood Wave Behavior in Natural Rivers as Function of Hydraulic and Hydrological Parameters of the Basin

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

88 MODERN SEWER DESIGN

Advanced /Surface Hydrology Dr. Jagadish Torlapati Fall 2017 MODULE 2 - ROUTING METHODS

FHWA Headloss in InfoWorks ICM. August 2013

Simulation of Transcritical Flow in Hydraulic structures

Development of Kanako, a wide use 1-D and 2-D debris flow simulator equipped with GUI

A NONLINEAR OPTIMIZATION MODEL FOR ESTIMATING MANNING S ROUGHNESS COEFFICIENT

FLOOD ROUTING FOR A SPECIFIC ORIENTATION OF PLANNED DEVELOPMENTS FOR AL-SHAMIYA RIVER IN IRAQ AS CASE STUDY

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.

Transactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X

Determination of Urban Runoff Using ILLUDAS and GIS

Stage Discharge Tabulation for Only Orifice Flow

IMPORTANCE OF NUMERICAL EFFICIENCY FOR REAL TIME CONTROL OF TRANSIENT GRAVITY FLOWS IN SEWERS

Abstract. 1 Introduction

Hydraulic Design Of Polyethylene Pipes

Experimental Study on Effect of Houses on Debris-Flow Flooding and Deposition in Debris Flow Fan Areas

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

CHANGES IN RIVER BED AROUND THE FUKAWA CONTRACTION AREA BY FLOODS AND CHANNEL IMPROVEMENT WORKS IN THE LOWER TONE RIVER

Open Channel Flow I - The Manning Equation and Uniform Flow COURSE CONTENT

MODELING FLUID FLOW IN OPEN CHANNEL WITH HORSESHOE CROSS SECTION

The most important equation to describe the water balance of a reservoir is the water balance: Equation 3.1

The most common methods to identify velocity of flow are pathlines, streaklines and streamlines.

Transactions on Ecology and the Environment vol 8, 1994 WIT Press, ISSN

1.060 Engineering Mechanics II Spring Problem Set 4

A finite volume method for solving the gravity wave-model equations

Flood Capacity of Shirakawa River at Tatsudajinnnai Area in Kumamoto Prefecture

Rock Sizing for Batter Chutes

Appendix E Guidance for Shallow Flooding Analyses and Mapping

Uniform Channel Flow Basic Concepts. Definition of Uniform Flow

UNIVERSITY OF BOLTON. ENGINEERING, SPORTS and SCIENCES BSC CIVIL ENGINEERING SEMESTER 1 EXAMINATION 2014/2015 WATER ENGINEERING MODULE NO: BLT3023

Notes: We all know that Toulmins Spring Branch is a sub-watershed of Three Mile Creek watershed. Some part of it is in Mobile area and rest of it is

Effects of Forming Conditions of Roll Offset Method on Sectional Shape at the Corner of Square Steel Pipe +

Experiment (4): Flow measurement

For example an empty bucket weighs 2.0kg. After 7 seconds of collecting water the bucket weighs 8.0kg, then:

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.

c 2007 by Arturo S. León. All rights reserved.

Flow in Open Channel Flow Conditions

Experimental Study for Investigating the Impact Force on a Wooden House by Driftwood in Steady and Unsteady Surge-type Flow

Block 3 Open channel flow

Basic Hydraulics June 2007

DYNAMICS OF FLOOD FLOWS AND BED VARIATIONS IN RIVER SECTIONS REPAIRED TO SHIP-BOTTOM SHAPED CHANNELS FROM COMPOUND CHANNLS

Two-dimensional urban flood simulation: Fukuoka flood disaster in 1999

A fuzzy dynamic wave routing model

Continuing Education Course #101 Drainage Design with WinTR-55

39.1 Gradually Varied Unsteady Flow

SOLUTIONS TO SELECTED PROBLEMS

Presented by: Civil Engineering Academy

HEAD LOSSES IN SEWER JUNCTION

Prediction of landslide-induced debris flow hydrograph: the Atsumari debris flow disaster in Japan

Growing and decaying processes and resistance of sand waves in the vicinity of the Tone River mouth

MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS

9. Flood Routing. chapter Two

Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati

STABILITY OF OLD STONE GROINS IN THE KATSURA RIVER DURING FLOODS

Project Description. Project Options. End Analysis On... Apr 26, :00:00. Rainfall Details

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

Design Methodology for Hydraulic Ram Pump

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

Two-Fluid Model 41. Simple isothermal two-fluid two-phase models for stratified flow:

INFLOW DESIGN FLOOD CONTROL SYSTEM PLAN 40 C.F.R. PART PLANT YATES ASH POND 2 (AP-2) GEORGIA POWER COMPANY

ENGINEERING HYDROLOGY

Application of high-resolution (10 m) DEM on Flood Disaster in 3D-GIS

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

CASE STUDIES. Introduction

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

EFFECT OF TWO SUCCESIVE CHECK DAMS ON DEBRIS FLOW DEPOSITION

This file was downloaded from Telemark Open Research Archive TEORA -

The Influence of Piers on the Risk of Flooding Upstream from a Bridge

Minimum Specific Energy and Critical Flow Conditions in Open Channels

Surge Analysis Using Transient Pressure Theory

HYDRAULIC EFFICIENCY OF GRATE AND CURB INLETS FOR URBAN STORM DRAINAGE

Physical landscapes River landscapes in the UK

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

Stage Discharge Prediction in a Prismatic Compound Channel

Application of Fuzzy Set Theory to Flood Risk Analysis in Natural Rivers as a Function of Hydraulic Parameters

Storm Water Best Management Practice: Development of Debris Filtering Structure for Supercritical Flow

The Sensitivity Analysis of Runoff from Urban Catchment Based on the Nonlinear Reservoir Rainfall-Runoff Model

Continuing Education Associated with Maintaining CPESC and CESSWI Certification

Reservoir Oscillations with Through Flow

Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process

Sewer, pressurization, differential pressure monitoring, fully dynamic hydraulic modeling, air displacement modeling.

Beaver Creek Corridor Design and Analysis. By: Alex Previte

Simulation of flow discharge on Danube River

Transcription:

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models T. UKON 1, N. SHIGETA, M. WATANABE 3 *, H. SHIRAISHI 4 and M. UOTANI 5 1 Department of Design and Development, NIHON KOGYO Co., Ltd., 61, Kaniwa, Miki, Kita, KAGAWA, 761-81, JAPAN. Graduate School of Civil and Environmental Engineering, Ehime University, 3, Bunkyo-cho, Matsuyama, EHIME, 79-8577, JAPAN. 3 Department of Civil and Environmental Engineering, Ehime University, 3 Bunkyo-cho, Matsuyama, EHIME, 79-8577, JAPAN. 4 Aratani Civil Engineering Consultants Co., Ltd., -1-, Yogo-naka, Matsuyama, EHIME, 79-45, JAPAN. 5 ZERO Co., Ltd., 199-1, Takasu-cho, Onomichi, HIROSHIMA, 79-141, JAPAN. *Corresponding author, e-mail nabemasa@dpc.ehime-u.ac.jp ABSTRACT The Preissmann slot model and the MOUSE slot model have been widely used to simulate a pressurized flow in urban drainage sewer pipe systems. In the runoff simulations utilizing theses two models, the slot width is often widened to make the computational time increment large and save the computing time. The resulting simulated water level, however, is likely to decrease to the extent that the simulated water level can not be considered for practical use. The main reason for the tendency mentioned above, in the case of the Preissmann slot model, is that the discharge through the slot, the friction loss is not taken into account, increases considerably according to the enlargement of the slot width, and the main reason, in the case of the MOUSE slot model, is that the excess storage of runoff water in the slot increases to a considerable extent according to the enlargement. In this paper, firstly, two correction methods that can improve the accuracy of the runoff simulations are proposed for the Preissmann slot model and the MOUSE slot model respectively. Secondly, the adaptabilities of the methods are examined through numerical runoff experiments. KEYWORDS Preissmann slot model, MOUSE slot model, Lateral model, manholes, pressurized flow, urban drainage sewer pipe system INTRODUCTION Recently, a localized torrential rain has happened frequently by global warming and heat island, and the occurrence of the flood damage in urban drainage sewer pipe systems has increased. The countermeasures against the flood damage, therefore, must be planned as soon as possible. In order to achieve this, it is necessary the simulation model that can simulate precisely the storm water runoff with a pressurized flow in the sewer pipe system has been developed. The Preissmann and MOUSE slot models have been widely used to simulate accurately the pressurized flow in urban drainage sewer pipe systems. In the runoff simulations utilizing these two models, the slot width is often widened without an assurance that the accuracy of T. Ukon, et al. 1

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 the simulation results does not decrease to make the computational time increment large and save the computing time. The resulting simulated water level, however, is likely to decrease to the extent that the simulated water level can not be considered for practical use, especially in the case of runoff simulations with flooding from the sewer pipe systems. In this paper, two correction methods that can improve the accuracy of the runoff simulations are proposed for the Preissmann slot model and the MOUSE slot model respectively, and the adaptabilities of the methods are examined through numerical runoff experiments. PREISSMANN SLOT MODEL, MOUSE SLOT MODEL AND LATERAL MODEL Preissmann slot model In the Preissmann slot model, a very narrow slot is started up in the top of the sewer pipe as shown in Figure 1, and then an actual pressurized flow can be handled as an open-channel flow. The following assumptions are made in the derivation of the flow equations of the model. (1) Water is incompressible. () The sewer pipe is rigid. (3) The slot wall does not take part in friction loss. (4) The slot width is very narrow. (5) The slot width is decided depending on the sound wave velocity of the pressure wave. Using the assumptions mentioned above, the flow equations of motion and continuity for the pressurized flow are derived as follows (Cunge and Wegner, 1964) : 1 V V V n V V + + I + = (1) 4 / 3 g t g R c + V + t g V = ga gb ga S p c = = a 1 + ( y D) ; a = (3) BS Ap BS where V = the mean cross-sectional velocity; y = the depth of flow; R = the hydraulic radius; I = the slope of the sewer pipe; n = Manning s roughness coefficient; c = the celerity of the small gravity wave; B S = the slot width; A = the flow cross-sectional area; A p = the crosssectional area of the sewer pipe; D = the diameter of the sewer pipe; a = the pressure wave celerity; g = the gravitational acceleration; t = the time; x = the distance along the sewer pipe. () (Surcharged Flow) Hypothetical Slot Free-surface (Open-channel Flow) Figure 1. Preissmann hypothetical slot. Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 MOUSE slot model In the MOUSE slot model, a narrow slot is started up in the top of the sewer pipe as shown in Figure, and then an actual pressurized flow can be handled as an open-channel flow. The following assumptions are made in the derivation of the flow equations of the model. (1) Water is compressible. () The sewer pipe wall is elastic. (3) The slot part acts as not a flow area but a storage area. (4) The slot width is very narrow. Using the assumptions mentioned above, the flow equations of motion and continuity for pressurized flow are derived as follows (Danish Hydraulic Institute, 3): Q Q α A + + ga = ga( I I f ) (4) t Q Q ρ gap + + = (5) ρ a t 1+ g( y D) 1 + g( y D) Er e a ρ = ρ, A = Ap, ar =, a = (6) a ar ρ D 1+ a / ar where Q = the discharge; α = velocity distribution coefficient; I f = the friction slope; ρ = the density of water for a free surface flow; a = the speed of sound in water; E r = the Young s modulus of elasticity; e = the pipe wall thickness. B = B S = ga p /a B > B S D Figure. Pipe with a fictitious slot. Lateral model In urban drainage sewer pipe systems, many lateral pipes are installed to transport storm water from inlets to sewer pipes. Considering the storage of runoff water in the lateral pipes as shown in Figure 3, the momentum and continuity equations for the pressurized flow in a sewer pipe are obtained (Watanabe et al. 199, 1991). 1 V V V n V V + + I + = (7) 4 / 3 g t g R T. Ukon, et al. 3

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 c + t g V L N = (8) ( sinθ ) A A P c = g (9) l where c = the pressure wave celerity; L = the length of the sewer pipe; A l, N, and θ = the cross-sectional area, number, and angle of the lateral pipes, respectively. Street Inlet Ground Surface Laterals House Inlet Ground Surface y D Laterals θ Combined Sewer Pipe Combined Sewer Pipe Figure 3. Sewer pipe and lateral pipes. CORRECTION METHOD FOR DROPPING OF SIMULATED WATER LEVEL In the runoff simulations utilizing these two models, the slot width is often widened without an assurance that the accuracy of the simulation results does not decrease to make the computational time increment large and save the computing time. The resulting simulated water level, however, is likely to decrease to the extent that the simulated water level can not be considered for practical use. The purpose of this paper is to develop the correction methods for dropping of simulated water level using Preissmann slot model and MOUSE slot model. Correction method for Preissmann slot model First, a correction factor, Г, is introduced in the momentum equation (1) of the Preissmann slot model and the factor is multiplied by the term of the friction loss. Second, the steady pressurized flow condition is assumed and Г is derived so that the hydraulic grade line of the Preissmann slot model coincides with that of the Lateral model. ( y D) B V V S full Γ = 1 + 1 F ± F A p c (1) c 4 / 3 R B S F = 1 + ( y D ) (11) 4 / 3 R A p where V full = the velocity of uniform pipe-full flow in the sewer pipe; R = the hydraulic radius except slot part; D = the length from invert to slot joint part; A p = the flow crosssectional area except slot part. 4 Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 The relation between the correction factor, Г, and y/d is shown in Figure 4. When the slot width is very narrow compared with the diameter of the sewer pipe (for example, B s =.1D), the simulated water level needs not being corrected because the discharge through the slot is very little compared with that through the sewer pipe. Therefore, Γ takes almost 1.. When the slot width is lager (for example, B s =.1D ~.5D), the discharge through the slot can not be ignored compared with that through the sewer pipe, and the simulated water level is likely to decrease because the slot wall does not have the friction loss. Therefore, in order to make the hydraulic grade line greater, Γ takes larger value than 1. according to the discharge through the slot. B S =.1D B S =.1D B S =.5D Γ 1. 1..8 1 3 y/d Figure 4. Correction factor, Г ~ y/d (D = 1. m). Correction method for MOUSE slot model The decrease of the simulated water level is due to the excess storage of runoff water, which is caused by widening the slot width, in the slot. In order to eliminate the excess storage, the same volume as the storage is reduced from the volume of the upper manhole as shown in Figure 5. The relation between the increment of the widened slot width and the reduced area of the manhole cross-sectional area is as follows: FM = α BS L (1) BS = B S BS (13) where FM = the reduced area of the manhole cross-sectional area; α = the correction coefficient (.5 α 1.); BS = the increment of the widened slot width; B S = the original slot width; B = the widened slot width. S V 1 = Excess storage of runoff water in slot Figure 5. Correction method for MOUSE slot model. V = Reduced volume of upper manhole V 1 = V T. Ukon, et al. 5

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 EXAMINATION OF ADDAPTABILITIES OF CORRECTION METHODS THROUGH NUMERICAL RUNOFF EXPERIMENTS The adaptabilities of the correction methods mentioned above were examined through the following numerical runoff experiments using the sewer pipe system as shown in Figure 6, Table 1, Table, and Table 3. The numerical runoff simulations under the following conditions were conducted: the inflow hydrographs at the upstream end of the pipe system are sinusoidal; the stage hydrographs at the downstream end are constant; the inundation on the ground surface never occurs; the periods of the inflow discharge hydrographs are 6 s and 1,8 s; the initial flow condition is the uniform pipe-full flow; the slot widths are assigned to be.1 m,.1 m, and.5 m. The simulated runoff hydrographs using the Preissmann slot model are shown in Figure 7 and Figure 8, and the ones using the MOUSE slot model are shown in Figure 9 and Figure 1. From these figures, it is apparent that the improvement of the dropping of the simulated water level as well as the accuracy of the simulated runoff hydrographs is achieved by making use of the correction methods. Q p Q Q full t T Level (m) 1 5 Q D Ground Surface L T G L T I Manholes 4 6 8 1 Distance (m) Figure 6. Sewer pipe system used in numerical experiments. Table 1. Dimensions of sewer pipe system used in numerical experiments. D F M L T G I N T (m) (m N ) (m) (m) ( ) (m -1/3 Q s) p (s) 6 1..33 6..5 3.8 17.13 1.3Q full 1,8 Table. Slot width and pressure wave celerity. (D = 1. m) B S (m) a (m/s).1 D.1 87.7.1 D.1 7.7.5 D.5 Table 3. Corrected manhole cross-sectional area, F' M. (MOUSE slot model) α B' S =.1 m B' S =.5 m F M F' M F M F' M.5.7.6 7.86.75.45 1.97.5.17 1..54 1.79.94 -.68 6 Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 9 8 B S =.1m B S =.1m B S =.5m. Water Level (m) 7 6 Discharge (m 3 /s) 1.8 5 6 1 18 6 1 18 6 1 18 6 1 18 (i) non- corrected (ii) corrected (i) non- corrected (ii) corrected (a) Water level hydrographs (b) Discharge hydrographs Figure 7. Simulated hydrographs (Preissmann slot model, T = 6 s). 9 8 B S =.1m B S =.1m B S =.5m. Water Level (m) 7 6 Discharge (m 3 /s) 1.8 5 1 4 36 1 4 36 1 4 36 1 4 36 (i) non- corrected (ii) corrected (i) non- corrected (ii) corrected (a) Water level hydrographs (b) Discharge hydrographs Figure 8. Simulated hydrographs (Preissmann slot model, T = 1,8 s). T. Ukon, et al. 7

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 9 8 B S =.1m B S =.1m B S =.5m. α =.75 α =.75 Water Level (m) 7 6 Discharge (m 3 /s) 1.8 5 6 1 18 6 1 18 6 1 18 6 1 18 (i) non- corrected (ii) corrected (i) non- corrected (ii) corrected (a) Water level hydrographs (b) Discharge hydrographs Figure 9. Simulated hydrographs (MOUSE slot mode, T = 6 s). 9 8 B S =.1m B S =.1m B S =.5m. α =.75 α =.75 Water Level (m) 7 6 Discharge (m 3 /s) 1.8 5 1 4 36 1 4 36 1 4 36 1 4 36 (i) non- corrected (ii) corrected (i) non- corrected (ii) corrected (a) Water level hydrographs (b) Discharge hydrographs Figure 1. Simulated hydrographs (MOUSE slot model, T = 1,8 s). 8 Correction methods for dropping of simulated water level utilising Preissmann and MOUSE slot models

11 th International Conference on Urban Drainage, Edinburgh, Scotland, UK, 8 CONCLUSIONS It was clarified that the dropping of the simulated water level caused by widening the slot width can be improved through the simulations making use of the correction methods proposed in this paper. REFERENCES Cunge, J. A. and Wegner, M. (1964). Numerical integration of Barré de Saint-Venant s flow equations by means of an implicit scheme of finite differences. La Houille Blanche., No.1, 33-39. Chaudhry, M. H. (1979). Applied Hydraulic Transients. VNR Company, New York, 43-448. Watanabe, M., Etoh, T. and Murota, A. (199). Runoff simulation model of sewer pipe systems considering pressure-relaxation effect by lateral pipes. Proc. 5 th Int. Conf. on Urban Storm Drainage..1, 75-8. Watanabe, M., Etoh, T. and Murota, A. (1991). Runoff simulation of sewer pipe systems with lateral pipes. Journal of Natural Disaster Science., 13(1), 19-8. MOUSE PIPE FLOW Reference Manual. (3). Danish Hydraulic Institute, Denmark, Chapter 3, 1-7. T. Ukon, et al. 9