ABSTRACT 1 INTRODUCTION
|
|
- Milton Newman
- 5 years ago
- Views:
Transcription
1 Water distribution network steady state simulation AQUANET model H. Rahal Laboratory of Hydraulics, K. U. Leuven, de Croylaan 2, B-3001 Heverlee, Belgium ABSTRACT This article is devoted to the description of the algorithm upon which the steady state simulation model AQUANET has been developed. This algorithm solves the nonlinear system of the continuity and the loop equations formulated in terms of the link flows, using the Newton-Raphson iterative procedure. To ensure fast convergence, the initial solution is determined using an extended version of the linear method. The state of the art graph theory concepts were used to generate the fundamental circuits (or loops). To take into account, the pipes various types of material and ages, the model allows the use of the most used head loss equation for each pipe link. Among these head loss equations, the Darcy-Colebrook can be used without any trade-off, because of the efficient procedure developed to ensure the fastest convergence. A special algorithm for parallel pipe assembly can be used to reduce the size of the vector space spanned by the steady state equations. Beside the determination of the link discharges and the nodes heads, the algorithm checks the bounding constraints imposed on the link velocity and on the node head pressure. Such feature allows the reinforcement and expansion of existing networks. 1 INTRODUCTION Water distribution steady state simulation has been a major issue for the last half century. The oldest systematic solution was provided by Hardy- Cross [1]. Not only his method addresses the problem correctly, but also is adapted to hand calculations. It suffers however, from slow convergence and instabilities.
2 396 Hydraulic Engineering Software The first general solution of the full set of equations using the Newton- Raphson's numerical algorithm to solve the nodal head pressure formulation is reported to be addressed by Martin and Peters [2] and extended by Shamir and Howard [5] to take into account even variable roughness coefficients, beside flows and heads. The theoretical study achieved by Nielsen [3] proved the numerical superiority of the link flow method on the nodal head method, which lacks stability and fast convergence. The proposed algorithm solves the full steady state equations expressed in terms of the link flow discharges. In the following there is a detailed description of the AQUANET model algorithm and the subsequent issues raised such as: the nonlinear system solution, the use of nonlinear head loss equations, the algorithm convergence, the initial solution generation, the loop generation and the indirect optimization features. 2. BASIC HYDRAULIC BACKGROUND A hydraulic network is composed of n nodes interconnected by p pipe links. The set of nodes is composed of rifh fixed head nodes ( reservoirs,...) and rifd fixed demand ( consumption ) nodes. Any closed chain of links is called a loop. One of the chains of links connecting two fixed head nodes is called a fundamental pseudo-loop. The graph theory cyclomatic number stipulates that the number p of links is equal to the sum of the m? real loops, the rrip pseudo-loops and rifd fixed demand nodes. p = rn,. +?7%p + rt/d (1) The steady state governing equations major difficulty is due to the nonlinearity of the head loss equations. The most used equations can be put in the resistance form: Ay = -^- 9" = r?" (2) K D Where hj is the friction head loss, L is the pipe length, D is the pipe diameter, q is the pipe discharge, K is the roughness coefficient, r is the pipe resistance,and a, ra and n are various constants and exponents depending on the chosen head loss formula, such as Hazen-Williams, Manning- S trickier, or Darcy-Wiesbach. In contrast to the formulas of Hazen-Williams and Manning-S trickier which are explicit, the Darcy equation coupled with Colebrook formula form an implicit system of equations, which needs an iterative procedure. The proposed procedure to solve such system will be elaborated in Section 5.2. Beside the friction head loss occurring along the pipe, there are singular head losses due to various localfittingssuch as: valves, elbows, contractions, junctions,... Losses in fittings have been found, in general, to vary as: hs = ksq^, where k$ is the singularity head loss coefficients.
3 Hydraulic Engineering Software 397 In contrast to all apparatus inducing energy losses, the pumps provide boosting energy excess, according to a given characteristic curve. This curve expresses the pumping head hp versus the discharge q. The most common fitting curve is the quadratic expression: hp = a + bq + cq*. The total head loss h along a pipe is the sum of the friction and local head losses from which the eventual boosting head is subtracted. Its expression versus the pipe discharge is as follows: /,,(9) - /4W = r<t + 6,g2 - (a + 6% + c^) (3) 3. THE MODEL GOVERNING EQUATIONS The retained link flow formulation expresses the problem in terms of the link flow discharges vector Q = [^i, f/2,..., </p]*, which consists of a nonlinear system formed by the continuity and the loop equations. The governing system is formed by rifd fixed demand nodes continuity equations, rrip pseudo-loop equations and ra^ real loops equations. 3.1 The Continuity equations The continuity equation, with respect to any fixed demand node j, stipulates that the algebraic sum of discharges of the nl links connected to this node, is equal to the node demand discharge: dj. Where Sij is the link ij direction sign, with respect to node j. The conventional direction s^ = 1, when node j corresponds to the link downstream node. On the contrary, Sij 1, when the node j corresponds to the pipe upstream node. The choice of the upstream and downstream nodes of a pipe link is arbitrary. 3.2 The loop equations Beside the classical loop, consisting of a closed circuit of pipes, the pseudoloop ( or hydraulic path) is the chain of pipes connecting two fixed head nodes (reservoirs). The choice between the different paths connecting two fixed head nodes is irrelevant. The pseudo-loop equation stipulates, that the algebraic sum of the head losses along the path connecting the two fixed head end-nodes, is equal to the difference between the z heads of the path upstream and downstream nodes. For pseudo-loop?, the head equation is: k=l (5)
4 398 Hydraulic Engineering Software Where 5/^ = 1 when the kth link direction corresponds to the loop chain's direction, on the contrary s^ = 1. /is the number of links forming the pseudo-loop. The loops equations are identical to those of the pseudo-loops, except that the head difference between the end nodes is Az% = The resulting nonlinear system assembly In terms of the link flow discharges vector, The steady state nonlinear system of p equations can be written as: /i(9n 92?' ',%) = 0 /or = 0 4. THE ALGORITHM SOLUTION 4.1 Linearization The solution of a nonlinear system goes through the use of a Taylor's series expansion around a known vector Q^ = [q[, q%, -,<??]* f the ith iteration. The Newton- Raphson procedure assumes that the second order term of the expansion is negligible and fk(q^^) ~ 0. Therefore, the retained first order linearization turns to be: The above equation can be expressed in terms of Ag^ q^ q '. It can be extended, as well, to the whole system using the jacobian matrix J, as: J(Q(*'))AQ(''+i) = -F(Q^) (7) A linear system solver ( Gauss elimination ) determines at each iteration, the vector AQ^+*\ which allows the determination of the new discharge vector: ' i) - Q(*') (8)
5 Hydraulic Engineering Software The algorithm convergence Mathematically, the convergence is reached, when the nonlinear system F(Q) = 0. However, numerically the convergence is considered reached, when the residual error: e[f(q')] < e, where e is a fixed precision. The residual error is estimated using the Euclidean norm: To guarantee the best solution, the Newton-Raphson algorithm has been damped. In fact, at each iteration the residual error is computed. If the new residual error is higher (worse) than the subsequent one, the step size AQ(*+^ is divided by 2, as many times as possible, until the new iteration gives better residual error. The convergence rate depends on the closeness of the initial solution. The generation of this initial solution, not only grantees the convergence, but also allows its fastest rate. This has been confirmed by the determination of the jacobian matrix condition, which turns to be, virtually, never ill-conditioned. k=l 5. THE ALGORITHM OTHER FEATURES 5.1 The initial solution generation The overall Newton-Raphson algorithm convergence rate depends particularly on the closeness of the initial solution guess to the exact one. In this respect, an algorithm has been developed to approximate as good as possible the initial guess estimate. This algorithm is based on the linear theory method proposed by Wood [6]. It consists of assembling a linear system of p link flow equations. The n/d continuity equations are already linear. However, the loop equations are nonlinear, making their linearization necessary. The proposed linearization is based on the fact that practically the order of magnitude of the pipe velocity and its Darcy-Weisbach friction factor are respectively of about 1 m/s and Thus, the loop equation becomes: > - f E^ g^* = Az, (10) 5.2 Darcy-Colebrook formula use The Darcy-Weisbach head loss formula is the most indicated one, particularly, when it is coupled with the Colebrook- White formula for determining the friction coefficient. Such system of equations has the advantage of taking into account the fluid viscosity v and of covering a wide range of flow
6 400 Hydraulic Engineering Software regimes. However, the two equations are implicit and require an iterative procedure. The author solved the combined nonlinear function expressed in terms of the friction head loss: = 0 (11) The proposed algorithm uses the Net won- Raphson iterative procedure, which requires the determination of f(hf) and it's derivative f'(hf). The rate covergence depends upon the closeness of the initial guess. The best initial guess is achieved by the use of Zigrang-Sylvester [7] explicit friction coefficient formula. The implementation of such algorithm showed outstanding convergence. In fact, not more than two iterations are necessary to converge to even a very high fixed precision. 5.3 The loop generation Determining manually the set of loops of a large size network turns to be very constraining and some times a source of errors. Therefore, developing a routine that can generate automatically the loops is extremely important, particularly if a pipe break has to be simulated, since the loops might change at each simulation run. The AQUANET model loop generator routine starts by generating a short path tree (instead of a spanning tree, which is less efficient) to identify the tree branches from the co-tree chords. For each chord, there is an associated fundamental circuit. This chord forms thefirstlink of the circuit, the rest of the circuit links form a unique set of branches. The short path ( or skim ) tree is generated using Disjkstra's procedure combined with a binary search strategy. The circuit branches are determined by a minimum path routine. The pseudo-loops are generated the same way as for the real loops, except that the chord is replaced by the pseudo-chord joining the two fixed head nodes. 5.4 Parallel pipe assembly Existing water supply system end up, in most cases, with several types of pipe material of different ages, necessitating the use of different head loss formulas. In addition, to enhance the efficiency of the system to meet increasing demand, adding parallel pipes to reinforce the system is a common practice. An algorithm to assemble parallel pipes using different head loss formula developed by Rahal and Berlamont [6] was implemented and added to the AQUANET model to reduce the size of the nonlinear system of as many equations as the number of parallel links.
7 Hydraulic Engineering Software 401 This algorithm determines the equivalent pipe diameter of a set of p parallel pipes, using each one a particular head loss formula. The expression of the equivalent diameter is: A i=p AT'^r (12) Where the subscripts e and i refer respectively to the equivalent pipe and to the ith parallel pipe. The equivalent length is chosen as the weighted ( on the diameters) average length of the parallel pipes system. The parallel pipes coefficients and formulas are known. However, a decision has to made upon the choice of the formula for the equivalent system. The algorithm overcomes this difficulty by choosing the formula which has the highest sum of pseudo-resistance CY]- L^/Df). 5.5 Pseudo-optimization Beside computing the links flows and nodes head pressures, the AQUANET model allows an "indirect optimization". In fact, it provides the system total cost for a part or the whole system, constituted of new pipes, and indicates if there is any violated constraints. The imposed constraints are those of the bounded link velocities: Vmin < V < Vmax, and of the bounded node heads: z^,^ < z < z^^. This indirect optimization is remarkably useful, particularly for existing distribution systems, which needs to be expanded or reinforced. In most cases, only few links ought to be added or reinforced. Therefore, only few computer runs are necessary to optimize the system, without going through sophisticated optimization algorithms. 5.6 Final remark In this paper, the author omitted intentionally to discuss the dimensionality problems posed by the jacobian matrix for large size networks. In fact, a great research effort is still underway, to find the appropriate solutions for this issue. The author is studying whether, it is appropriate to add a sparse matrix routine solver or to reformulate the system of nonlinear equations through graph theory and various substitutions. 6. CONCLUSIONS The AQUANET model is developed, essentially to simulate the steady state in a water distribution network. This model has been under development
8 402 Hydraulic Engineering Software for the last few years during which the author tried to include as many efficient procedure as possible. The global algorithm convergence has been carefully studied by damping the Newton-Raphson procedure in a way the best solution is retained. Furthermore, this convergence is accelerated by generating an outstanding initial solution estimate based on an improvement of the linear method. The AQUANET model has not shown for the last few years, any cases of nonconvergence, except those due to input data errors. In addition the Jacobian matrix condition determination did not show any ill-conditioned cases. Finally, the author come up with the same conclusions as Nielsen [3] indicating that the best convergence is reached using the linear method for the first iteration and the Newton-Raphson for the subsequent iterations. Beside the algorithm convergence issue, many procedures have been developed with great success such as: the starting initial solution estimate, the use of comprehensive head loss formulas, the loops generation, the parallel pipes assembly and design constraints verifications which determines whether or not a specific design is feasible with respect to the link velocity and nodal head bound constraints. 7. REFERENCES [1]. Hardy Cross (1936). "Analysis of flow in networks of conduits or conductors." BWWm <?##, f&e C/muerszZ?/ o////nzozg Em/. EzperzmemZ ^Wzom, Urbana, Illinois. [2]. Martin, D. W., and Peters, G. (1963). "The application of Newton Raphson's method of network analysis by digital computers.",journal of the Institute of Water Engineers, 17, [3]. Nielsen, H. B. (1989). "Methods for analyzing pipe networks." Journal o/ffz/drwz'c #%#., 115 (2), [4]. Rahal, H. and Berlamont J.(1993). "Algorithm for parallel pipe assembly ",/%2e<7raW Compiler App/zmZzoms m MWer J)%pp/%/, l/b/.); Me2/to& &W Procedures /or ^z/sfems j^muwzmz &?W CWZro/, edited by Bryan Coulbek, Research Studies Press Lts and John Wiley Sons Inc. pp [5]. Shamir, U., and Howard, C. D. D. (1968). "Water distribution systems analysis." JowrW o/f/te Tf^raWzcs Dzm^orn, ASCE, 94 (HY1), [6]. Wood, D. J. and Charles, C.O.A (1972)." Hydraulic network analysis using linear theory", Jo%rW o/#%/<fr. DzmszW, ASCE, 98 (HY7), [7]. Zigrang, D.J. and Sylvester, N.D. (1982)." Explicit approximation to the solution of Colebrook's friction factor.", Journal of Amer. Inst. of 28 (3),
Fluids Engineering. Pipeline Systems 2. Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET
COURSE NUMBER: ME 423 Fluids Engineering Pipeline Systems 2 Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 SERIES PIPE FLOW WITH PUMP(S) 2 3 4 Colebrook-
More informationAn element by element algorithm for pipe network analysis
An element by element algorithm for pipe network analysis M.H. Afshar* and A. Afshar^ Imam Higher Education Center, Tehran, Iran * University ofscience and Technology, Tehran, Iran Email: mahabnet@dci.rian.com
More informationSTEADY FLOW THROUGH PIPES DARCY WEISBACH EQUATION FOR FLOW IN PIPES. HAZEN WILLIAM S FORMULA, LOSSES IN PIPELINES, HYDRAULIC GRADE LINES AND ENERGY
STEADY FLOW THROUGH PIPES DARCY WEISBACH EQUATION FOR FLOW IN PIPES. HAZEN WILLIAM S FORMULA, LOSSES IN PIPELINES, HYDRAULIC GRADE LINES AND ENERGY LINES 1 SIGNIFICANCE OF CONDUITS In considering the convenience
More informationHydraulics Prof Dr Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati
Hydraulics Prof Dr Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati Module No # 08 Pipe Flow Lecture No # 04 Pipe Network Analysis Friends, today we will be starting
More informationChapter (3) Water Flow in Pipes
Chapter (3) Water Flow in Pipes Water Flow in Pipes Bernoulli Equation Recall fluid mechanics course, the Bernoulli equation is: P 1 ρg + v 1 g + z 1 = P ρg + v g + z h P + h T + h L Here, we want to study
More informationThe Darcy-Weisbach Jacobian and avoiding zero flow failures in the Global Gradient Algorithm for the water network equations
The Darcy-Weisbach Jacobian and avoiding zero flow failures in the Global Gradient Algorithm for the water network equations by Elhay, S. and A.R. Simpson World Environmental & Water Resources Congress
More informationНаучный потенциал регионов на службу модернизации. Астрахань: АИСИ, с.
MODELING OF FLOWS IN PIPING TREES USING PROJECTION METHODS В.В. Войков, Астраханский инженерно-строительный институт, г. Астрахань, Россия Jason Mayes, Mihir Sen University of Notre Dame, Indiana, USA
More informationJournal of Hydraulic Engineering
Journal of Hydraulic Engineering Discussion of "Method to Cope with Zero Flows in Newton Solvers for Water Distribution Systems" by Nikolai B. Gorev, Inna F. Kodzhespirov, Yuriy Kovalenko, Eugenio Prokhorov
More informationA Spreadsheet Tool for the Analysis of Flows in Small-scale Water Piping Networks
A Spreadsheet Tool for the Analysis of Flows in Small-scale Water Piping Networks Kazeem B. Adedeji 1, Yskandar Hamam 1, Bolanle T. Abe 1 and Adnan M. Abu-Mahfouz 1,2 1 Department of Electrical Engineering,
More informationReview of pipe flow: Friction & Minor Losses
ENVE 204 Lecture -1 Review of pipe flow: Friction & Minor Losses Assist. Prof. Neslihan SEMERCİ Marmara University Department of Environmental Engineering Important Definitions Pressure Pipe Flow: Refers
More informationNONLINEAR DC ANALYSIS
ECE 552 Numerical Circuit Analysis Chapter Six NONLINEAR DC ANALYSIS OR: Solution of Nonlinear Algebraic Equations I. Hajj 2017 Nonlinear Algebraic Equations A system of linear equations Ax = b has a
More informationcha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO
cha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO ROOTS OF EQUATIONS PT2.1 MOTIVATION Years ago, you learned to use the quadratic formula x = b ± b 2 4ac 2a to solve f(x) = ax 2 + bx + c = 0 (PT2.1) (PT2.2)
More informationChapter 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 informationViscous Flow in Ducts
Dr. M. Siavashi Iran University of Science and Technology Spring 2014 Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate
More informationAnalysis of Complex Pipe Networks with Multiple Loops and Inlets and Outlets
Analysis of Complex Pipe Networks with Multiple Loops and Inlets and Outlets The techniques described previously for analysis of pipe flow are satisfactory if the pipe system is simple, consisting of one
More informationReynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment:
7 STEADY FLOW IN PIPES 7.1 Reynolds Number Reynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment: Laminar flow Turbulent flow Reynolds apparatus
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationA SPREADSHEET TEACHING TOOL FOR ANALYSIS OF PIPE NETWORKS
Engineering Journal of University of Qatar, Vol. 10, 1997, p. 33-50 A SPREADSHEET TEACHNG TOOL FOR ANALYSS OF PPE NETWORKS Aly Nabih El - Bahrawy Associate Professor, Civil Engineering Department University
More informationPiping Systems and Flow Analysis (Chapter 3)
Piping Systems and Flow Analysis (Chapter 3) 2 Learning Outcomes (Chapter 3) Losses in Piping Systems Major losses Minor losses Pipe Networks Pipes in series Pipes in parallel Manifolds and Distribution
More informationRelationship between Hazen-William coefficient and Colebrook-White friction factor: Application in water network analysis
European Water 58: 513-520, 2017. 2017 E.W. Publications Relationship between Hazen-William coefficient and Colebrook-White friction factor: Application in water network analysis M. Niazkar 1, N. Talebbeydokhti
More informationChapter 8: Flow in Pipes
Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate the major and minor losses associated with pipe flow in piping networks
More informationAn overview of the Hydraulics of Water Distribution Networks
An overview of the Hydraulics of Water Distribution Networks June 21, 2017 by, P.E. Senior Water Resources Specialist, Santa Clara Valley Water District Adjunct Faculty, San José State University 1 Outline
More information1-Reynold s Experiment
Lect.No.8 2 nd Semester Flow Dynamics in Closed Conduit (Pipe Flow) 1 of 21 The flow in closed conduit ( flow in pipe ) is differ from this occur in open channel where the flow in pipe is at a pressure
More informationCalculation of Pipe Friction Loss
Doc.No. 6122-F3T071 rev.2 Calculation of Pipe Friction Loss Engineering Management Group Development Planning Department Standard Pump Business Division EBARA corporation October 16th, 2013 1 / 33 2 /
More informationExperiment (4): Flow measurement
Experiment (4): Flow measurement Introduction: The flow measuring apparatus is used to familiarize the students with typical methods of flow measurement of an incompressible fluid and, at the same time
More informationMajor and Minor Losses
Abstract Major and Minor Losses Caitlyn Collazo, Team 2 (1:00 pm) A Technovate fluid circuit system was used to determine the pressure drop across a pipe section and across an orifice. These pressure drops
More informationUse of a smoothed model for pipe friction loss
Use of a smoothed model for pipe friction loss Bradley J. Eck, M.ASCE Please cite this as: Eck, B. (2016). Use of a Smoothed Model for Pipe Friction Loss. J. Hydraul. Eng., 10.1061/(ASCE)HY.1943-7900.0001239,
More informationChapter (3) Water Flow in Pipes
Chapter (3) Water Flow in Pipes Water Flow in Pipes Bernoulli Equation Recall fluid mechanics course, the Bernoulli equation is: P 1 ρg + v 1 g + z 1 = P ρg + v g + z h P + h T + h L Here, we want to study
More informationENGI9496 Modeling and Simulation of Dynamic Systems Bond Graphs
ENGI9496 Modeling and Simulation of Dynamic Systems Bond Graphs Topics covered so far: Analogies between mechanical (translation and rotation), fluid, and electrical systems o Review of domain-specific
More informationAnalysis Methods for Water Distribution Systems-3 rd Class )طرق تحليل أنظمة توزيع المياه( Dr. Sataa A. Al-Bayati (10-11)
تسم هللا الزحمه الزحيم Analysis Methods for Water Distribution Systems-3 rd Class )طرق تحليل أنظمة توزيع المياه( Dr. Sataa A. Al-Bayati (10-11) Methods of analysis are: )المقاطغ( Sectioning.1 )االوثىب
More information4 Pipelines and Pipe Network Hydraulics I
4 Pipelines and Pipe Network Hydraulics I The hydraulics of pipelines and pipe networks presented herein is limited to turbulent flow of water in closed conduits (pipes) flowing full. Closed conduits flowing
More informationAppendix F. + 1 Ma 1. 2 Ma Ma Ma ln + K = 0 (4-173)
5:39p.m. Page:949 Trimsize:8.5in 11in Appendix F F.1 MICROSOFT EXCEL SOLVER FOR NON-LINEAR EQUATIONS The Solver is an optimization package that finds a maximum, minimum, or specified value of a target
More informationCVE 372 HYDROMECHANICS EXERCISE PROBLEMS
VE 37 HYDROMEHNIS EXERISE PROLEMS 1. pump that has the characteristic curve shown in the accompanying graph is to be installed in the system shown. What will be the discharge of water in the system? Take
More informationENG3103 Engineering Problem Solving Computations Semester 2, 2013
Assessment: Assignment 2 Due: 16 September 2013 Marks: 100 Value: 10 % Question 1 (70 marks) Introduction You are designing a pipe network system that transfers water from the upper pipe to the lower pipe.
More informationCHAPTER 4 HYDRAULICS OF WATER DISTRIBUTION SYSTEMS
CHAPTER 4 HYDRAULICS OF WATER DISTRIBUTION SYSTEMS Kevin Lansey Department of Civil Engineering and Engineering Mechanics University of Arizona Tucson, AZ Larry W. Mays Department of Civil and Environmental
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationOPEN 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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationFACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4)
FACULTY OF CHEMICAL & ENERGY ENGINEERING FLUID MECHANICS LABORATORY TITLE OF EXPERIMENT: MINOR LOSSES IN PIPE (E4) 1 1.0 Objectives The objective of this experiment is to calculate loss coefficient (K
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationLOSSES DUE TO PIPE FITTINGS
LOSSES DUE TO PIPE FITTINGS Aim: To determine the losses across the fittings in a pipe network Theory: The resistance to flow in a pipe network causes loss in the pressure head along the flow. The overall
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
More informationFLOW FRICTION CHARACTERISTICS OF CONCRETE PRESSURE PIPE
11 ACPPA TECHNICAL SERIES FLOW FRICTION CHARACTERISTICS OF CONCRETE PRESSURE PIPE This paper presents formulas to assist in hydraulic design of concrete pressure pipe. There are many formulas to calculate
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationWATER DISTRIBUTION NETWORKS
WATER DISTRIBUTION NETWORKS CE 370 1 Components of Water Supply System 2 1 Water Distribution System Water distribution systems are designed to adequately satisfy the water requirements for a combinations
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationFlow in Pipes Simulation Using Successive Substitution Method
Flow in Pipes Simulation Using Successive Substitution Method Dr. Ismail Al-Rawi Arab Open University, (Kuwait Branch), P.O 830 Al Ardia, P.C 92400, Kuwait Abstract The aim of this paper is to present
More informationCOMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE
COMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE Rafael ÅMAN*, Heikki HANDROOS*, Pasi KORKEALAAKSO** and Asko ROUVINEN** * Laboratory
More informationLinear Second-Order Differential Equations LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS
11.11 LINEAR SECOND-ORDER DIFFERENTIAL EQUATIONS A Spring with Friction: Damped Oscillations The differential equation, which we used to describe the motion of a spring, disregards friction. But there
More informationMixed Integer Programming Solvers: from Where to Where. Andrea Lodi University of Bologna, Italy
Mixed Integer Programming Solvers: from Where to Where Andrea Lodi University of Bologna, Italy andrea.lodi@unibo.it November 30, 2011 @ Explanatory Workshop on Locational Analysis, Sevilla A. Lodi, MIP
More informationSaudi Journal of Engineering and Technology. DOI: /sjeat ISSN (Print)
DOI:10.21276/sjeat.2017.2.3.1 Saudi Journal of Engineering and Technology Scholars Middle East Publishers Dubai, United Arab Emirates Website: http://scholarsmepub.com/ ISSN 2415-6272 (Print) ISSN 2415-6264
More informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationspatial water demand (population density) spatial junctions distribution (building density) digital elevation map
SYSTEM ID: WDS-Designer NARRATIVE DESCRIPTION With the WDS Designer a tool for the algorithmic generation of synthetic water distribution systems (swds) based on GIS data was presented (Sitzenfrei et al.,
More informationComparative Analysis of Techniques for Solving the Hydraulics of Pressurized Irrigation Pipe Networks. Numan Mizyed
An-Najah Univ. J. Res., Vol. 11, (1997) 1-21 Comparative Analysis of Techniques for Solving the Hydraulics of Pressurized Irrigation Pipe Networks Numan Mizyed Plant Production Department, Faculty of Agriculture,
More informationChapter 3 Water Flow in Pipes
The Islamic University o Gaza Faculty o Engineering Civil Engineering Department Hydraulics - ECI 33 Chapter 3 Water Flow in Pipes 3. Description o A Pipe Flow Water pipes in our homes and the distribution
More informationComparison of MOC and Lax FDE for simulating transients in Pipe Flows
International Research Journal of Engineering and Technology (IRJET) e-issn: 395-0056 Volume: 04 Issue: 03 Mar -07 www.irjet.net p-issn: 395-007 Comparison of MOC and Lax FDE for simulating transients
More informationEnhancing Computer-Based Problem Solving Skills with a Combination of Software Packages
3420 Enhancing Computer-Based Problem Solving Sills with a Combination of Software Pacages Mordechai Shacham Dept. of Chemical Engineering, Ben Gurion University of the Negev Beer-Sheva 84105, Israel e-mail:
More informationME 305 Fluid Mechanics I. Part 8 Viscous Flow in Pipes and Ducts. Flow in Pipes and Ducts. Flow in Pipes and Ducts (cont d)
ME 305 Fluid Mechanics I Flow in Pipes and Ducts Flow in closed conduits (circular pipes and non-circular ducts) are very common. Part 8 Viscous Flow in Pipes and Ducts These presentations are prepared
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More information1. Method 1: bisection. The bisection methods starts from two points a 0 and b 0 such that
Chapter 4 Nonlinear equations 4.1 Root finding Consider the problem of solving any nonlinear relation g(x) = h(x) in the real variable x. We rephrase this problem as one of finding the zero (root) of a
More informationComparison of Explicit Relations of Darcy Friction Measurement with Colebrook-White Equation
Applied mathematics in Engineering, Management and Technology 2 (4) 24:57-578 www.amiemt-journal.com Comparison of Explicit lations of arcy Friction Measurement with Colebrook-White Equation Saeed Kazemi
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationApplying Solver Engine in Microsoft Excel to Analyze Water Distribution Network based on Linear Theory Method
J. Saf. Environ. Health Res. 1(2): XX XX, Winter 2017 DOI: 10.22053/jsehr.2017.79056.10** ORIGINAL RESEARCH PAPER Applying Solver Engine in Microsoft Excel to Analyze Water Distribution Network based on
More informationChapter 10 Flow in Conduits
Chapter 10 Flow in Conduits 10.1 Classifying Flow Laminar Flow and Turbulent Flow Laminar flow Unpredictable Turbulent flow Near entrance: undeveloped developing flow In developing flow, the wall shear
More informationPressure and Flow Characteristics
Pressure and Flow Characteristics Continuing Education from the American Society of Plumbing Engineers August 2015 ASPE.ORG/ReadLearnEarn CEU 226 READ, LEARN, EARN Note: In determining your answers to
More informationME 305 Fluid Mechanics I. Chapter 8 Viscous Flow in Pipes and Ducts
ME 305 Fluid Mechanics I Chapter 8 Viscous Flow in Pipes and Ducts These presentations are prepared by Dr. Cüneyt Sert Department of Mechanical Engineering Middle East Technical University Ankara, Turkey
More informationGuo, 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 informationA New Technique for Solving Sets of Coupled. Nonlinear Algebraic and Integro-Differential. Equations Encountered in Hydraulics
Int. J. Contemp. Math. Sciences, Vol. 3, 2008, no. 33, 1611-1617 A New Technique for Solving Sets of Coupled Nonlinear Algebraic and Integro-Differential Equations Encountered in Hydraulics Amir Nassirharand
More informationCS 450 Numerical Analysis. Chapter 5: Nonlinear Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationAn Equivalent Circuit Formulation of the Power Flow Problem with Current and Voltage State Variables
An Equivalent Circuit Formulation of the Power Flow Problem with Current and Voltage State Variables David M. Bromberg, Marko Jereminov, Xin Li, Gabriela Hug, Larry Pileggi Dept. of Electrical and Computer
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationHydraulics 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 informationEXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER
EXPERIMENT No.1 FLOW MEASUREMENT BY ORIFICEMETER 1.1 AIM: To determine the co-efficient of discharge of the orifice meter 1.2 EQUIPMENTS REQUIRED: Orifice meter test rig, Stopwatch 1.3 PREPARATION 1.3.1
More informationHydraulic (Fluid) Systems
Hydraulic (Fluid) Systems Basic Modeling Elements Resistance apacitance Inertance Pressure and Flow Sources Interconnection Relationships ompatibility Law ontinuity Law Derive Input/Output Models ME375
More informationFlood 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 informationDynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches
Dynamics of Multibody Systems: Conventional and Graph-Theoretic Approaches by Professor John J. McPhee, P.Eng. Systems Design Engineering (Cross-appointed to Mechanical Engineering) University of Waterloo,
More informationExperimental and Numerical Investigations of the Effect of Net Positive Suction Head on Water Hammer In Pipeline Systems
International Journal of Engineering and Advanced Technology (IJEAT) ISSN: 2249 8958, Volume-3, Issue-1, October 2013 Experimental and Numerical Investigations of the Effect of Net Positive Suction Head
More informationTwo-Layer Network Equivalent for Electromagnetic Transients
1328 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 18, NO. 4, OCTOBER 2003 Two-Layer Network Equivalent for Electromagnetic Transients Mohamed Abdel-Rahman, Member, IEEE, Adam Semlyen, Life Fellow, IEEE, and
More informationAn Expression for Obtaining Total Heads for Lift Pump Selection
American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-03, Issue-06, pp-169-176 www.ajer.org Research Paper Open Access An Expression for Obtaining Total Heads for
More informationSizing of Gas Pipelines
Sizing of Gas Pipelines Mavis Nyarko MSc. Gas Engineering and Management, BSc. Civil Engineering Kumasi - Ghana mariiooh@yahoo.com Abstract-In this study, an effective approach for calculating the size
More informationModule 10: Finite Difference Methods for Boundary Value Problems Lecture 42: Special Boundary Value Problems. The Lecture Contains:
The Lecture Contains: Finite Difference Method Boundary conditions of the second and third kind Solution of the difference scheme: Linear case Solution of the difference scheme: nonlinear case Problems
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationHydraulics and hydrology
Hydraulics and hydrology - project exercises - Class 4 and 5 Pipe flow Discharge (Q) (called also as the volume flow rate) is the volume of fluid that passes through an area per unit time. The discharge
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 3 3.4 Differential Algebraic Systems 3.5 Integration of Differential Equations 1 Outline 3.4 Differential Algebraic Systems 3.4.1 Constrained Dynamics 3.4.2 First and Second
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationPressure Head: Pressure head is the height of a column of water that would exert a unit pressure equal to the pressure of the water.
Design Manual Chapter - Stormwater D - Storm Sewer Design D- Storm Sewer Sizing A. Introduction The purpose of this section is to outline the basic hydraulic principles in order to determine the storm
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationImproved Method for Converting Equivalent Sand-grain Roughness to Hazen-Williams Coefficient
Proceedings of the 2 nd World Congress on Mechanical, Chemical, and Material Engineering (MCM'16) Budapest, Hungary August 22 23, 2016 Paper No. HTFF 119 OI: 10.11159/htff16.119 Improved Method for Converting
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationApplication of Mathematical Modeling to Study Flood Wave Behavior in Natural Rivers as Function of Hydraulic and Hydrological Parameters of the Basin
Hydrology Days 2010 Application of Mathematical Modeling to Study Flood Wave Behavior in Natural Rivers as Function of Hydraulic and Hydrological Parameters of the Basin Patrícia Freire Chagas 2, Silvia
More informationThe Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization
Physics Electricity & Magnetism fields Okayama University Year 2004 The Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization Koji Fujiwara
More informationLecture Note for Open Channel Hydraulics
Chapter -one Introduction to Open Channel Hydraulics 1.1 Definitions Simply stated, Open channel flow is a flow of liquid in a conduit with free space. Open channel flow is particularly applied to understand
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationFriction Factors and Drag Coefficients
Levicky 1 Friction Factors and Drag Coefficients Several equations that we have seen have included terms to represent dissipation of energy due to the viscous nature of fluid flow. For example, in the
More informationPresented at MASSEE International Conference on Mathematics MICOM-2009
M a t h e m a t i c a B a l k a n i c a New Series Vol. 26, 202, Fasc. 3 4 Lambert W Function in Hydraulic Problems Dejan Brkić Presented at MASSEE International Conference on Mathematics MICOM-2009 Darcy
More informationCHAPTER 11. A Revision. 1. The Computers and Numbers therein
CHAPTER A Revision. The Computers and Numbers therein Traditional computer science begins with a finite alphabet. By stringing elements of the alphabet one after another, one obtains strings. A set of
More informationTrust Regions. Charles J. Geyer. March 27, 2013
Trust Regions Charles J. Geyer March 27, 2013 1 Trust Region Theory We follow Nocedal and Wright (1999, Chapter 4), using their notation. Fletcher (1987, Section 5.1) discusses the same algorithm, but
More informationREE 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