A Paper presentation on. Department of Hydrology, Indian Institute of Technology, Roorkee

Similar documents
SHOALING OF NONLINEAR INTERNAL WAVES ON A UNIFORMLY SLOPING BEACH

Multi-Fuel and Mixed-Mode IC Engine Combustion Simulation with a Detailed Chemistry Based Progress Variable Library Approach

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

P R = P 0. The system is shown on the next figure:

Numerical Simulation of the Dispersion of a Plume of Exhaust Gases from Diesel and Petrol Engine Vehicles

EVALUATION OF FORCE COEFFICIENTS FOR A 2-D ANGLE SECTION USING REALIZABLE k-ε TURBULENCE MODEL

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

Robustness Experiments with Two Variance Components

Variants of Pegasos. December 11, 2009

ELASTIC MODULUS ESTIMATION OF CHOPPED CARBON FIBER TAPE REINFORCED THERMOPLASTICS USING THE MONTE CARLO SIMULATION

Midterm Exam. Thursday, April hour, 15 minutes

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Linear Response Theory: The connection between QFT and experiments

12d Model. Civil and Surveying Software. Drainage Analysis Module Detention/Retention Basins. Owen Thornton BE (Mech), 12d Model Programmer

Cubic Bezier Homotopy Function for Solving Exponential Equations

FTCS Solution to the Heat Equation

Volatility Interpolation

Diffusion of Heptane in Polyethylene Vinyl Acetate: Modelisation and Experimentation

Solution in semi infinite diffusion couples (error function analysis)

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Research Article Cubic B-spline for the Numerical Solution of Parabolic Integro-differential Equation with a Weakly Singular Kernel

Data Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data

PHYS 1443 Section 001 Lecture #4

On One Analytic Method of. Constructing Program Controls

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Polymerization Technology Laboratory Course

Solving Parabolic Partial Delay Differential. Equations Using The Explicit Method And Higher. Order Differences

Outline. Energy-Efficient Target Coverage in Wireless Sensor Networks. Sensor Node. Introduction. Characteristics of WSN

Introduction to. Computer Animation

Numerical simulation of a solar chimney power plant in the southern region of Iran

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

WiH Wei He

TSS = SST + SSE An orthogonal partition of the total SS

January Examinations 2012

Accuracy of the Muskingum-Cunge method for constant-parameter diffusion-wave. channel routing with lateral inflow

Density Matrix Description of NMR BCMB/CHEM 8190

Comprehensive Integrated Simulation and Optimization of LPP for EUV Lithography Devices

NUMERICAL SIMULATION AND EXPERIMENTAL INVESTIGATION FOR INDOOR AIR ENVIRONMENT IN AN OFFICE ROOM

10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :

Anisotropic Behaviors and Its Application on Sheet Metal Stamping Processes

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

T q (heat generation) Figure 0-1: Slab heated with constant source 2 = q k

Relative controllability of nonlinear systems with delays in control

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

ON THE EFFICIENCY OF STORMWATER DETENTION TANKS IN POLLUTANT REMOVAL

Time-interval analysis of β decay. V. Horvat and J. C. Hardy

Density Matrix Description of NMR BCMB/CHEM 8190

2.1 Constitutive Theory

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

PERISHABLES INVENTORY CONTROL MODEL UNDER TIME- VARYING AND CONTINUOUS DEMAND

Stochastic Repair and Replacement with a single repair channel

Scattering at an Interface: Oblique Incidence

ID-1193 STUDY ON IMPACT BEHAVIOR OF COMPOSITES USING THE FINITE ELEMENT METHOD INTRODUCTION

Nonequilibrium models for a multi component reactive distillation column

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

DEEP UNFOLDING FOR MULTICHANNEL SOURCE SEPARATION SUPPLEMENTARY MATERIAL

Review of Numerical Schemes for Two Point Second Order Non-Linear Boundary Value Problems

Optimal Operation of the Cyclic Claus Process

Genetic Algorithm in Parameter Estimation of Nonlinear Dynamic Systems

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

Li An-Ping. Beijing , P.R.China

VEHICLE DYNAMIC MODELING & SIMULATION: COMPARING A FINITE- ELEMENT SOLUTION TO A MULTI-BODY DYNAMIC SOLUTION

Measurement of liquid holdup and axial dispersion in trickle bed reactors using radiotracer technique

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

Real time processing with low cost uncooled plane array IR camera-application to flash nondestructive

ON THE WEAK LIMITS OF SMOOTH MAPS FOR THE DIRICHLET ENERGY BETWEEN MANIFOLDS

Bayes rule for a classification problem INF Discriminant functions for the normal density. Euclidean distance. Mahalanobis distance

Chapter 2 Linear dynamic analysis of a structural system

A General Lagrangian Approach to Simulate Pollutant Dispersion in Atmosphere for Low-wind Condition

Computing Relevance, Similarity: The Vector Space Model

Transient Numerical of Piston Wind in Subway Station. Haitao Bao

IIR Band Pass and Band Stop Filter Design Employing Teaching-Learning based Optimization Technique

Single-loop System Reliability-Based Design & Topology Optimization (SRBDO/SRBTO): A Matrix-based System Reliability (MSR) Method

Structure Analysis and Optimization for Compacting System of Asphalt Paver

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Thermal fatigue estimation due to thermal stratification in the RCS branch line using one-way FSI scheme

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Additional Exercises for Chapter What is the slope-intercept form of the equation of the line given by 3x + 5y + 2 = 0?

Time Scale Evaluation of Economic Forecasts

Reactive Methods to Solve the Berth AllocationProblem with Stochastic Arrival and Handling Times

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

Gear System Time-varying Reliability Analysis Based on Elastomer Dynamics

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

Modeling and Solving of Multi-Product Inventory Lot-Sizing with Supplier Selection under Quantity Discounts

ESTIMATION OF HYDRAULIC JUMP LOCATION USING NUMERICAL SIMULATION

Sampling Procedure of the Sum of two Binary Markov Process Realizations

ISSN MIT Publications

02. MOTION. Questions and Answers

Investigation of Sway, Roll and Yaw Motions of a Ship with Forward Speed: Numerical modeling for flared up conditions

Comparison between the Discrete and Continuous Time Models

Chapters 2 Kinematics. Position, Distance, Displacement

Stability Analysis of Fuzzy Hopfield Neural Networks with Timevarying

RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA

(Time-dependent) Mean-field approaches to nuclear response and reaction

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

Fall 2010 Graduate Course on Dynamic Learning

Mohammad H. Al-Towaiq a & Hasan K. Al-Bzoor a a Department of Mathematics and Statistics, Jordan University of

Transcription:

A Paper presenaon on EXPERIMENTAL INVESTIGATION OF RAINFALL RUNOFF PROCESS by Ank Cakravar M.K.Jan Kapl Rola Deparmen of Hydrology, Indan Insue of Tecnology, Roorkee-247667

Inroducon Ranfall-runoff processes plays an mporan role n ydrologcal cycle. Te movemen of waer as manly arsen from e need o evaluae e amoun of avalable runoff waer a a parcular locaon o mee local demand as well as rsk of floodng due o excess waer. A number of models suc as pyscally based and concepual models ave been used o smulae e ranfall runoff process. However, due o s complexy and spao-emporal varaon, few models can accuraely smulae s gly non-lnear process.

Runoff s generaed by ransorms and s occurrence and quany are dependen on e caracerscs of e ranfall even,.e. nensy, duraon and dsrbuon. Te expermenal nvesgaons could be performed o undersand surface, sub-surface or ground waer runoff generaon mecansm. Te mecansm of surface runoff generaon s owever e focus of e presen nvesgaon. Te man objecve of s laboraory sudy s o quanfy e ranfall runoff process.

Descrpon of e laboraory se-up Te ranfall smulaor s desgned o smulae e ranfall even over e cacmen. Basc componens of e smulaor are- Tlng flume for slope adjusmen(recangular cross-secon avng.25 m wde and 5 m long), Spray Nozzles, Sorage Tank, Pressure Gauges, Elecrc Moor Slope Adjusmen Devce (0-5%) Runoff Recordng Sysem.

Fg. Advanced Hydrologc sysem S-2 MKII-50 (Ranfall Smulaor)

Meodology Te expermens s carred ou on a sand layer placed over an mpermeable plane surface (smoo meal see), Sedmen parcle dameer of 0.5-.0 mm s requred o fll over e flume. Smlar arrangemen for ranfall nensy from 30 mm/r o 90 mm/r. W e elp e spray nozzle arfcal ranfall s spread over e cacmen area n e flume flled w sand. Durng ranfall even waer flows over surface of sand reaces o e oule and e runoff waer flows roug collecng ank. Te Dep sensor s aaced below e collecng ank o gve e eg over wer.

Maemacal Modelng KW formulaon as been wdely used for modelng overland and sream flow. KW formulaon uses pysograpc parameers suc Overland rougness, slope, dranage area and leng, and sol caracerscs for compuaon of overland flow. Te KW approxmaon of San-Venan's equaon can be wren as:- A ( va ) x = q l.() m Q = αa.(2)

Were; A s e area of waersed (m 2 ), v s e flow velocy (m/s); Q s flow rae (m 3 /s), s e me (s), x s e dsance measured n e drecon of flow (m) α and m are parameers of e knemac wave model wc are closely relaed o e caracerscs of e flow. A α ma m A x = q l (3) α o m o m o x = r f (4) q o m = α o (5) o

Te governng dfferenal equaon were solved numercally usng explc numercal meod. From equaon () and (2) usng mplc fne dfference meod, = q x q ( r f ) (6) ( ) m q = α ( ) m q = α.(7). (8) ( ) m q = α. (9) Were r = ranfall nensy n mm/r, f = nflraon rae n mm/r, = 5 / 3, α o = s / 2 o n o m o

By solvng e above equaon usng Newon-Rapson erave meod, ( ) ( ) ( ) ( ) k k k k f f / = ( ) ( ) / = m m x f α Were Afer smplfyng e ese equaon, ( ) ( ) ( ) ( ) = m m f r x x f α α. (0).. (). (2)

Nas Suclffe Creron Were, NES(%) = n n ( Y ( Y o o Y Y c m ) ) 2 2 00 Y o Y c Y m = Observed flow (Expermenal ranfall flow) value a me = Compued flow (knemac flow) value a me = Mean of observed values. Error n Runoff Volume Compuaon Volumerc error, n % Y c = volume compued Y o = volume observed Y c 00 Y = o

Analyss of e daa Expermenal daa obaned from Advanced Hydrologc sysem (ranfall smulaor) s analyzed usng one dmensonal KW overland flow smulaon model usng FORTRAN programmng. Te smulaon was done usng me sep of 5 sec and spaal grd sze s aken as cm. Expermenal daa were smulaed usng developed model n order o sudy e overland flow rougness due o varaon of slope and ranfall nensy on e cacmen. Daa was observed for cacmen slope beween % o 4 % and ranfall nensy vares beween 30 o 90 mm/r. Developed KWM was calbraed for mannng s rougness coeffcen n o f no daa observed from expermen for 30 mm/r ranfall nensy for dfferen slopes.

Table Pernen caracerscs of observed and compued ydrograp for 30 mm/r nensy of ranfall. Volume Tme o Peak Slope NSE S.No Observed Compued Error Observed Compued Error (%) (l) (l) (%) (mn) (mn) (%) (%) 9.3 7.9 4.3 2.7 2.4. 98 2 2 6.3 5.9 6.3 2.3 2. 8.6 97.2 3 3 8.3 7.9 4.8 2. 2.0 4.7 9 4 4 7.6 7.3 3.9 2.0.8 5 93.2

Table 2 Pernen caracerscs of observed and compued ydrograp for 60 mm/r nensy of ranfall Volume Tme o Peak S.No Slope Observed Compued Error Observed Compued Error NSE (%) (l) (l) (%) (mn) (mn) (%) (%) 2.5 2.2 2.4 3.9 3.5 0.2 97.8 2 2 0.4 0. 2.8 3.2 2.9 9.3 96.6 3 3 6.2 5.8 2.4 2.8 2.6 7. 97.0 4 4 2.5 2. 3.2 2.5 2.4 4 92.6

Table 3 Pernen caracerscs of observed and compued ydrograp for 90 mm/r nensy of ranfall Volume Tme o Peak Slope NSE S.No Observed Compued Error Observed Compued Error (%) (l) (l) (%) (mn) (mn) (%) (%) 6. 5.7 2.4 3.6 3.3 8.3 96 2 2 8.3 7.9 2. 3. 2.9 6.4 97 3 3 5.4 4.9 3.2 3.2 3 6.2 96. 4 4 5. 4.7 2.6 2.6 2.5 3.8 96.7

Fg. 2 Comparson of observed and compued ydrograp for ran fall nensy 30 mm/r a % slope. Fg. 3 Comparson of observed and compued ydrograp for ran fall nensy 30 mm/r a 2 % slope.

Fg 4 Comparson of observed and compued ydrograp for ranfall nensy 60mm/r a 2% slope of e plane. Fg 5 Comparson of observed and compued ydrograp for ranfall nensy 60mm/r a 3% slope of e plane.

Fg 6 Varaon of Mannng s rougness w overland flow plane slope for ranfall nensy 30 mm/r.

Fg. 7 Varaon of Mannng s rougness w ranfall nensy for % slope of e plane.

Dscusson of resuls Te ploe beween observed and model compuaed ydrograp are sown n Fg. 2 o 5 for slope of %, 2%, and 3% respecvely I can be seen a, afer e sar of e ranfall on e cacmen, e waer flow rae a e oule ncreases rapdly w me, s poron of ydrograp s known as e rsng lmb. However a a ceran me e waer flow rae a e oule equals e ranfall nensy, s me s known as me o peak. Wen e ranfall nensy s sopped e waer flow rae a e cacmen oule sar reducng a a very slow rae, s s known as e fallng lmb of e ydrograp.

Te KWM could reproduce resonably well rsng lmb of e observed ydrograp as well as seady sae dscarge. In case of fallng lmb of e ydrograp, wo segmens are clearly vsble. Were only surface runoff domnaes wc s reproduced well by KWM. Te release of waer from sand bed, wc e KWM as no smulaed because e mecansm of release of waer from sand bed was no ncorporaed n e KWM. Form Fg 6 e value of n decreases lnearly w ncrease n overland flow plane slope. Ts may be due o rapd dranng of e waer due o ncreased surface slope.

Concluson Te overland runoff ydrograps presened n s sudy was conduced n e laboraory w equpmen conssng of a ranfall smulaor (nozzle spray), an mpermeable overland flow plane. Te resuls ndcae a ere s less dfference n e amoun of runoff volumes, peaks a dfferen ranfall nenses and slope of e cacmens. Bu consderable dfference n fallng lmb of ydrograp sapes. Laboraory and feld expermens wll es relaonsps for a wder range of condons a wll nclude e use of oer ranfall nensy, dfferen slopes and nflrang surfaces. Fuure work sould also nclude e comparson of e expermenal resuls w numercal resuls.e. knemac wave eory.

Tank You