Defra / Environment Agency FCERM Joint Science Programme. SC030218/SR Afflux Estimation System. August 2007 HYDRAULIC REFERENCE
|
|
- Sydney Hampton
- 5 years ago
- Views:
Transcription
1 Defra / Environment Agency FCERM Joint Science Programme SC030218/SR Afflux Estimation System August 2007 HYDRAULIC REFERENCE South Barn Broughton Hall SKIPTON North Yorkshire BD23 3AE UK t: +44 (0) f: +44 (0)
2
3 REVISION HISTORY Revision Ref./ Date Issued Technical notes: 5 December 2005 AES Hydraulic Reference Notes August 2007 Amendments Issued to Andrew Pepper, Mervyn Bramley, Technical Review Group CES/AES web site CONTRACT These research notes describe work commissioned by under Science Group contract Ref. SC s representative for the contract was Andrew Pepper of ATPEC. Prof. Peter Mantz, Dr. Serter Atabay, Dr. Rob Lamb, Steve Rose, Claire Packett and Jeremy Benn of carried out the work. Prepared by:... Peter Mantz, BSc, MSc, PhD, CEng, CPhys, MICE, MASCE = = = = = = mêáååáé~ä=^å~äóëí= Reviewed by:... Rob Lamb, MA PhD = = = = = = qéåüåáå~ä=aáêéåíçê= Approved by: Jeremy Benn, MA, MSc, CEng, FICE, FCIWEM, MASCE = = = = = = j~å~öáåö=aáêéåíçê= DISCLAIMER This document has been prepared solely as a series of research notes for. JBA Consulting accepts no responsibility or liability for any use that is made of this document or associated software other than by the Client for the purposes for which it was originally commissioned and prepared. i
4 ACKNOWLEDGMENTS We are grateful to the Defra/Environment Agency joint flood management research programme (Engineering Theme) for funding this study. In particular, we would like to acknowledge the support given to the project by Dr Mervyn Bramley OBE (formerly Theme Leader) and Andrew Pepper (External Theme Advisor). We are grateful for expert comment and review by Prof. Donald Knight, Dr Les Hamill, Prof. Nigel Wright, Dr Paul Samuels and Dr John Riddell. We are also grateful to the above for helping with the provision of data on afflux, and for having the opportunity to support a number of MSc laboratory projects at Birmingham University through this project. We would also like to thank Dr Galip Seckin for assistance with data. ii
5 EXECUTIVE SUMMARY This document This technical reference note describes the background and development of the Afflux Estimator. It has been compiled for circulation to the project Technical Review Group. The Afflux Estimation System (AES) Accurate estimation of flood water level underpins almost all sectors of flood risk management. A new support tool - the Afflux Estimation System (AES) for bridges and culverts - has been commissioned by the Environment Agency as part of the joint Defra/EA flood research programme to ensure that best available methods relevant to conditions in the UK are used for flood risk management. The Afflux Estimation System will comprise two main outputs: a simple Afflux Advisor providing quick, approximate calculations and guidance in an accessible spreadsheet form a more rigorous Afflux Estimator, produced as open source code, and available to be implemented as a software application or an additional component of existing river modelling packages. Afflux Estimator Afflux Estimator (AE) is being developed to model one dimensional water surface profiles at varied depths above and below a bridge or culvert located in a stream. The profiles are generated using the standard step method, which is based upon the principle of conservation of energy. Additionally, the code can model single bridge piers using pier drag coefficients together with the principle of conservation of momentum. To date, the bridge code is in place and sufficient testing has been undertaken to enable the development of AE with the Conveyance Estimation System software (CES, 2004). The AE model is designed for simplicity, therefore only three bridge reaches are used. These are the two transition reaches from the undisturbed stream sections upstream and downstream of the structure to the bridge (the flow contraction and expansion reaches respectively) and the bridge reach itself. The alongstream lengths of these transition reaches are estimated using a Transition Calculator. The latter computes algorithms derived from the downstream flow Froude number, the bridge Blockage ratio (the ratio of the area of blockage cause by the bridge and the total flow area at the bridge), and the relative floodplain width. The algorithms were derived from recent laboratory transition length measurements (Atabay and Knight, 2002), and field transition length data attained by computer simulation (HEC, 1995A). The transition calculator also estimates the associated, transition energy loss coefficients from the same data. There are four arch bridge types modelled in the code and two beam bridge types. The arch types are semicircular, parabolic, elliptical and user-defined. The beam types are continuous pier (for which the pier extends beneath the full bridge width) and single pier (for which the pier does not extend beneath the full bridge width). These model types have been tested by AE using the recent laboratory data from the University of Birmingham (Atabay and Knight, 2002; Seckin et al, 2004), and field data from large rivers surrounded by densely wooded floodplains (USGS, 1978). iii
6 This page is intentionally left blank. iv
7 CONTENTS Page REVISION HISTORY CONTRACT EXECUTIVE SUMMARY CONTENTS LIST OF FIGURES LIST OF TABLES NOTATION i i iii v vi vii viii 1 PRELIMINARY ANALYSIS FOR TRANSITION LENGTHS Summary The HRC method Transition length determinations for the UB (2002) data Transition length analyses for the UB (2002) data Transition length analyses for the HEC (1995) data Transition length analyses for combined UB (2002) and HEC (1995) data Discussion of the HEC (1995) transition length analyses Design tables for transition lengths Appendices to Chapter BACKWATER FLOW THROUGH A CONTRACTION Summary The standard step method Uniform flow for the UB(2002) Case 2 laboratory channel Uniform flow for the CES (2004) River Main Uniform flow for the USCS (1978) Yellow River Contracted flow for the UB(2002) Case 1 and Case 5 laboratory channel Contracted flow for the CES (2004) River Main Contracted flow for the USGS (1978) Yellow River FLOW PROFILES FOR A SINGLE ARCH AND BEAM BRIDGE Summary The Bridge Data Entry worksheet Conveyance estimation The Backwater code Rating Comparisons for a single parabolic arch and single beam bridge PHYSICAL CONSTRAINTS FOR BRIDGE DIMENSIONS Summary Combined Weir and Orifice modes Single Arch and Beam bridge of 20 m span and high road level Single Arch bridge of 10 m span and low road level Single Arch bridge of 10 m span and very low soffit level Single Arch bridge of eccentric 20 m span Summary of physical constraint data MULTI-SPAN BRIDGES Summary Numerically testing the Multi-arch code Numerically testing the Multi-beam code Comparison of AE and HEC-RAS discharge ratings v
8 5.5 Summary of Multi-span data BEAM BRIDGES WITH SINGLE PIERS Summary CMM used in HEC-RAS CMMR used in AE Pier drag coefficients (Cd) CMMR applied to a 20 m span beam bridge CMMR applied to a piered beam bridge compared with HEC-RAS modelling AE CALIBRATION FOR LABORATORY SCALE BRIDGES Summary The AE model AE calibration with UB data Comparison with HEC-RAS AE examples with UB data Appendix to Chapter 7 (summary of afflux data) AE CALIBRATION FOR FIELD SCALE BRIDGES Summary The AE model Calibration principles for the transition parameters Transition reach lengths Transition reach energy coefficients AE and HEC-RAS comparison for the Yellow River (USGS, 1978, Survey 610) DEVELOPMENT AND VERIFICATION OF AE FOR USE WITH FIELD SURVEYS Summary Arch bridge verification Arch bridge results Beam bridge verification Appendix 1: Summary of USGS (1978) data for modelling water surface profiles in AE COMPARISON OF AE CULVERT CODE WITH HEC-RAS MODELS Summary Design of the AE code for culverts Shallow culvert modelling Comparison with HEC-RAS models for the shallow culverts Deep culvert modelling Comparison with HEC-RAS models for deep culverts vi
9 LIST OF FIGURES Figure 1-1: Side elevation at a bridge contraction...2 Figure 1-2: Summary of Afflux data plotted on HRC dimensionless coordinates...2 Figure 1-3: Estimation of contraction lengths for the UB experiments...5 Figure 1-4: Estimation of contraction lengths for the UB Beam experiments...6 Figure 1-5: Estimation of expansion lengths for the PBEAM3 experiments...7 Figure 1-6: Estimation of expansion lengths for the AMOSC and BEAM1 experiments...8 Figure 1-7: Plot of Contraction length (Lc) against Afflux (dh) with linear trend...10 Figure 1-8: Plot of Expansion length (Le) against Afflux (dh) with possible linear trend...10 Figure 1-9: Variation of dimensionless contraction length with F3 and J3, and predicted trended length against measured length...11 Figure 1-10: Variation of dimensionless expansion length with F3 and J3, and predicted trended length against measured length...12 Figure 1-11: River cross section and data used in the HEC (1995) computational model length...13 Figure 1-12: Variation of dimensionless contraction length with F3 and J3 and predicted trended length against measured length...14 Figure 1-13: Variation of dimensionless expansion length with F3 and J3, and predicted trended length against measured length...15 Figure 1-14: Predicted against measured dimensionless contraction lengths (upper curve) and expansion lengths (lower curve) for the UB (open symbols) and HEC data...17 Figure 1-15: Predicted against measured dimensionless contraction lengths (upper curve) and expansion lengths (lower curve) for combined UB and HEC data...18 Figure 2-1: Plan and profile for flow through a river contraction...29 Figure 2-2: River reach for the standard step water surface computations...31 Figure 2-3: Uniform flow for the UB Case 2 channel...33 Figure 2-4: Uniform flow for the River Main...34 Figure 2-5: Uniform flow for the Yellow River...35 Figure 2-6: AE profile and Section worksheets for the UB Case 1 channel...37 Figure 2-7: AE profile and Section worksheets for the UB Case 5 channel...38 Figure 2-8: AE profile and Section worksheets for the CES (2004) River Main...39 Figure 2-9: AE profile and Section worksheets for the Yellow River...40 Figure 3-1: Subsoffit flow profile for a single arch bridge...42 Figure 3-2: Supersoffit flow profile for a single arch bridge...43 Figure 3-3: Bridge Data Entry Worksheet...43 Figure 3-4: Cross section data for Section 2D...44 Figure 3-5: Bridge and Culvert surface roughness coefficients...46 Figure 3-6: Bridge and Culver surface roughness coeffcients vii
10 Figure 3-7: Section 2D for a multi-arch bridge...47 Figure 3-8: Section 1 worksheet...48 Figure 3-9: Section 2D conveyance rating...49 Figure 3-10: Extrapolation methods for AE and HEC-RAS...50 Figure 3-11: Comparison ratings for a single parabolic arch bridge (AE is the Structure curve)...52 Figure 3-12: Comparison ratings for a single beam bridge (AE is the Structure curve)...53 Figure 4-1:The Affluc Estimator application...55 Figure 4-2: The four modes of flow modelled in AE...56 Figure 4-3: Combined Weir and Orifice mode...57 Figure 4-4: Rating curves at Section 4 for the unstructured (Rating 4) and structured flow...58 Figure 4-5: Rating curves at Section 4 for an arch bridge of 10m span (low road level)...59 Figure 4-6: Rating curves at Section 4 for an arch bridge of 10m span (low soffit level)...60 Figure 4-7: Rating curves at Section 4 for an asymmetrical arch bridge of 20m span...60 Figure 5-1: Multi-span bridge examples...63 Figure 5-2: Multi-arch bridge models...65 Figure 5-3: Multi-beam bridge models...66 Figure 5-4: Twin Arch Model and Ratings...67 Figure 5-5: Twin Beam Model and Ratings...68 Figure 5-6: 20 Arch model and ratings...69 Figure 5-7: 20 Beam model and ratings...70 Figure 5-8: HEC-RAS profiles for the 20 Beam bridge...70 Figure 6-1: Beam bridge with continuous and single piers...74 Figure 6-2: Cross section near and inside a bridge (from HEC-RAS Hydraulic Reference, Fig 5.3)...75 Figure 6-3: Discharge ratings for a 20m beam bridge with varied piers and openings...79 Figure 6-4: AE Bridge Data Entry worksheet for an 11 opening, 20m Beam bridge...80 Figure 6-5: AE and RAS 10 pier models for a 20m span Beam bridge...81 Figure 6-6: AE and RAS 19 pier models for an 80m span Beam bridge...81 Figure 7-1: Correlation of computed and measured afflux...84 Figure 7-3: Water surface profile for a multiple arch bridge (from UB, 2002)...85 Figure 7-4: Cross section data used in AE...85 Figure 7-5: Transitional lengths for the UB data...87 Figure 7-6: Afflux comparisons for the subcritical flow computations...88 Figure 7-7: Affflux comparisons for the critical flow computations...89 Figure 7-8: Afflux comparisons for all flow computations...90 Figure 7-9: HEC-Ras and AE afflux comparisons for ASOE data...90 Figure 7-10: Discharge ratings for the ASOE Case 2 experiment...91 Figure 8-2: Water surface profile for a multiple arch bridge (from UB, 2002) Figure 8-3: Submerged orifice flow modelled in Afflux Estimator Figure 8-4: River reach for the standard step water surface computations Figure 8-5: AE Transition calculator viii
11 Figure 8-6: Section 1 worksheet in AE showing friction slope calibrated for Q=6.6m 3 /s Figure 8-7: Bridge Data Entry for the USGS (1978) Survey Figure 8-8: Discharge ratings for the Q=56.6 m 3 /s flood event Figure 8-8: Discharge ratings for the Q=187.0 m 3 /s flood event LIST OF TABLES Table 1-1: Summary of the number of experimental discharges for the UB normal bridge flows with measured transition lengths....4 Table 1-2: HEC Contraction Ratios (from HEC, 1995, Table 15)...19 Table 1-3: Trended Contraction Ratios...19 Table 1-4: Difference between HEC and trended Contraction Ratios...20 Table 1-5: HEC Expansion Ratios (From HEC 1995, Table 15)...20 Table 1-6: Trended Expansion Ratios...20 Table 1-7: Difference between HEC and trended Expansion Ratios...20 Table 1-8: Simplified design for transition lengths...21 Table 1-9: Simplified design for transition lengths...22 Table 1-10: Detailed design for dimensionless expansion lengths (Le/D3)...23 Table 2-1: Comparison of cross section statistics for the UB Case 2 channel...33 Table 2-2: Comparison of cross section statistics for the River Main...34 Table 2-3: Comparison of cross section statistics for the Yellow River...35 Table 4-1: Summary of examples...61 Table 5-4: Multi-span examples...71 Table 6-1: Coefficients used for pier geometry (from HEC-RAS V3, 2001)...77 Table 6-2: Coefficients used for pier geometry and 6:1 aspect ratio (from Montes, 1998)...78 Table 6-4: Multi-pier examples...82 Table 7-4: Bridge data entry for the Case 3 examples...93 Table 8-1: AE Transition data examples for laboratory and field scales...99 Table 8-2: Mean and observed Water surface elevations modelled in HEC-RAs (metric units) Table 8-3: Water surface elevations modelled in HEC-RAS (from HEC, 1995B, feet units) Table 8-4: Bridge flow modes for different codes (Zd is the downstream water level at the bridge and Zu is the upstream water level) ix
12 NOTATION Mathematical symbols and notation have been used in the same way as in the original publications and are introduced in the text at each relevant equation. Since there is considerable variation in the notation used by different authors and organisations, this document uses the original notation in all places. Some symbols therefore have different meanings in different parts of the text. x
13 xi
14 1 PRELIMINARY ANALYSIS FOR TRANSITION LENGTHS 1.1 Summary A recent series of laboratory experiments which measured afflux for bridges that are normal to a compound channel flow has been conducted at the University of Birmingham, UB (Atabay and Knight, 2002). A total of 145 afflux measurements were made for bridge types of single semicircular openings (20 experiments, abbreviated ASOSC), multiple semi-circular openings (15 experiments, abbreviated AMOSC), single elliptical openings (15 experiments, abbreviated ASOE) and straight deck beam bridges with or without piers (95 experiments, abbreviated PBEAM and BEAM). A preliminary analysis of the total data was made using a dynamic similarity method first invoked at the University of Roorkee (Raju et al, 1983), and later by Hydraulics Research at Wallingford (Brown, 1988). A simplified version of this method involved expressing bridge afflux (dh) in dimensionless terms as follows: dh/d3 = f(f3, J3) where D3 is the normal flow depth without the presence of the structure, F3 is the Froude number of the flow without the presence of the structure, and J3 is a blockage ratio to the flow caused by the structure, and is given by the quotient: J3 = area of flow blocked by the downstream bridge section total flow area in the undisturbed channel These dimensionless variables were successfully applied to the UB (2002) laboratory data and to USGS field data (USGS, 1978) for use with the Afflux Advisor. Since the data were from compound channels, the variables are named herein as HRC dimensionless variables. They are now used for determining transition lengths for a structure located in a compound channel. The HEC (1995) model for a structural contraction in a subsoffit stream flow is illustrated for a bridge in Figure 1.1. Channel sections differ in nomenclature to that of the HRC model, and the downstream normal flow depth, y1 (D3 in HRC), is located at Section 1. There are sequentially 4 sections to the location of maximum backwater or afflux at Section 4. The alongstream distances between Sections 1 and 2 and Sections 3 and 4 are called the transition lengths. The former is a flow expansion length (Le) and the latter is a flow contraction length (Lc). The following notes analyse the transition lengths in terms of the HRC dimensionless variables. It is shown that Lc was almost a constant length for the UB (2002) data, and Le increased with afflux. Quadratic trending of the HRC variables were used to give predicted values as: Lc/D3 = f(j3, F3) ; Le/D3 = f(j3,f3) The HEC (1995) computational data were similarly analysed, and it was found that the effect of river scale was also a dominant variable. This scale variable was subjectively taken as the floodplain width, and a final predictive algorithm was derived for both narrow floodplain laboratory data (UB, 2002) and wide floodplain computational data (HEC, 1995). 1
15 Lc Bridge Le Backwater Flow Afflux y 4 y 3 y 2 y 1 Section 4 Section 3 Section 2 Section 1 Figure 1-1: Side elevation at a bridge contraction 1.2 The HRC method The HRC method is preferred as a similarity model for afflux hydraulics because the independent variables used are explicit and the dimensionless variables are physically bounded. The blockage ratio, J3 is dependent on the structure blockage area as a fraction of the total flow area. It is therefore specific to the structure. The Froude number (V3/ (gd3) 0.5 ) applies to the flow without the structure. It is therefore explicit to the undisturbed stream. The HRC method is simply expressing the fractional afflux as a function of the fractional structure blockage area and the Froude number of the undisturbed stream. 1.2 dh/d UB Low friction UB Medium friction USGS High friction HR Arch field data J3=0.9 J3=0.7 J3= J3=0.3 J3= F3 Figure 1-2: Summary of Afflux data plotted on HRC dimensionless coordinates 2
16 In contrast, the BPR similarity method uses a scale velocity that is the average velocity through the structure at normal depth (D3). This scale velocity is now an implicit variable which is also dependent on the other dependent variable J3 (which is a function of the BPR opening area). The BPR similarity method thus involves implicit variables which may lead to spurious correlations. The afflux data previously considered for the Afflux Advisor is summarised in terms of HRC variables in Figure 1.2. Both J3 and F3 have values from zero to unity. The plot is bounded at the abcissa by J3 = 0 and at the ordinate by J3 = 1. Hypothesised curves from about F3 = 0.8 bound the right side of the plot from the maximum J3 afflux conditions to F3 = 1. (For the latter condition of critical flow the afflux is theoretically zero, as indicated by the approximate HR Arch bridge data.) If the symbol hs represents the height of a structure soffit above bed level, then supersoffit flow (called the soffit threshold) begins approximately when: hs = dh +D3 or dh/d3 = (hs/d3) - 1 An approximate estimate for hs can be made using bridge design criteria. It is sometimes recommended that the soffit level should exceed the 200 year flood level (TAC, 2004). This flood level may be approximated as about twice the normal depth. The soffit threshold is thus about dh/d3 = 1, and this line provides the upper bound for Figure 1.2. The HRC plot is therefore applicable for subsoffit flows only. Note that with the present Afflux Advisor, the J3 curves are discontinuous at F3 = 0.1(being optimised for USGS field data) and F3 = 0.25 (being optimised for the medium friction UB data). These discontinuities may be smoothed in the future with an associated small loss in the accuracy of afflux prediction. Since the HRC method has successfully predicted afflux for the entire subsoffit flow through a structure, it is now tested on UB transition length data. 1.3 Transition length determinations for the UB (2002) data The available flow discharges for the determination of transition lengths are summarised in Table 1.1. For many experiments, supercritical flow conditions occurred through the structure and thus a flow control was located at the structure as well as downstream. For some cases, a hydraulic jump occurred and the flow was also unsteady. An expansion length could not be measured for these conditions (AMOSC-Case 1 abbreviated to AMOSC1, ASOE1, ASOSC1). For the remaining flows, the transition levels were measured at 0.1 m increments downstream for large water surface gradients (such as at standing waves) and at 1 m increments downstream for small gradients. All water level measurements were made along the channel centreline. Section 3 was located at the m mark from upstream and Section 2 at the m mark. Both sections were therefore taken as the 59.0 m mark for simplicity in transition length estimation. The complete Lc measurements are illustrated in Figures 1.3 and 1.4. For objectivity, Section 4 was taken as the first maximum elevation upstream of the 59.0 m bridge section position. This method contrasted with the method used in HEC (1995), which used maximum water surface contours over the entire cross section. Representative examples for the expansion length measurements are given in Figures 1.5 and 1.6. The use of a 2 point moving average trend for the water levels was found adequate to dampen the standing waves in this reach, and to establish an objective intersect with the normal water level profile. Figure 1.5 (piered beam experiments) illustrates the smallest Le of about 2 m encountered for a high friction (and low afflux) condition. Figure 1.6 (AMOSC2 experiments) illustrate a medium Le of about 4 m, and the largest Le of about 7 m is illustrated for the lowest friction (and highest afflux) BEAM1 experiments. Note that for the latter measurements only, a moving average trend of more than 2 points was sometimes necessary for an objective assessment. This method for Le determination differed considerably from the HEC (1995) method. For the latter, velocity vector plots were used to establish uniformity of the flow for Section 1. Furthermore, the Le estimated herein have expansion ratios (defined in HEC, 1995 as Le/Lobs, where Lobs is the average obstruction length equal to 0.5*(floodplain width bridge width)) that are much larger than the HEC (1995) values. It is however considered below that the differences can be explained in terms of the different scale and friction of the floodplain widths. Note that the objective estimates for the UB (2002) transition length data used herein are considered unambiguous. 3
17 Table 1-1: Summary of the number of experimental discharges for the UB normal bridge flows with measured transition lengths. = Case 1 Case 2 Case 3 Case 4 Case 5 oçìöüåéëë= j~áå=åü~ååéä=j~ååáåöûë=å= MKMN= MKMNR= MKMOU= MKMOV= MKMPO= cäçççéä~áå=j~ååáåöûë=å= MKMMV= MKMP= MKMQT= MKMSV= MKMUM= Bridge Type Bridge code Multiple semi-circular AMOSC 5 (Lc) 5 (Lc, Le) 5 (Lc, Le) arch Single elliptical arch ASOE 5 (Lc) 5 (Lc, Le) 5 (Lc, Le) Single semi-circular ASOSC 5 (Lc) 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) arch Beam (Span = 0.398) BEAM-4 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) Beam (Span = 0.498) BEAM-5 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) Beam (Span = 0.598) BEAM-6 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) Piered Beam (Span = PBEAM-4 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) 0.398) Piered Beam (Span = PBEAM-5 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) 0.498) Piered Beam (Span = PBEAM-6 5 (Lc, Le) 5 (Lc, Le) 5 (Lc, Le) 0.598) kçíéw== The total number of experiments with measured contraction length Lc is 145, and with measured expansion length Le is
18 Figure 1-3: Estimation of contraction lengths for the UB experiments (The abcissa is the centreline channel position of measurement from upstream in metres, the ordinate is the water level elevation in millimetres. The small black square denotes the contraction length upstream limit.) 5
19 Figure 1-4: Estimation of contraction lengths for the UB Beam experiments (The abcissa is the centreline channel position of measurement from upstream in metres, the ordinate is the water level elevation in millimetres for both normal and bridge flows.) 6
20 Figure 1-5: Estimation of expansion lengths for the PBEAM3 experiments (The abcissa is the centreline channel position of measurement from upstream in metres, the ordinate is the water level elevation in millimetres for both normal and bridge flows.) 7
21 Figure 1-6: Estimation of expansion lengths for the AMOSC and BEAM1 experiments (The abcissa is the centreline channel position of measurement from upstream in metres, the ordinate is the water level elevation in millimetres for both normal and bridge flows.) 1.4 Transition length analyses for the UB (2002) data The total data for transition lengths against HRC variables are summarised in Appendix 1.1. As a first appraisal, both Lc and Le were compared against afflux (dh) in Figures 1.7 and 1.8. It was seen that Lc was approximately constant with afflux thus: Lc = 0.50 ± 0.02 m at the 95% confidence level Since the average obstruction lengths, Lobs were mainly about 0.4 m for these data, the contraction ratio (Lc/Lobs) was about 1.4 (This value was similar to those predicted in HEC, 1995). In contrast, Le increased with afflux and a probable power law trend would give a higher correlation than the possible linear trend shown. 8
22 Physically, these results imply that the main flow control is at the structure. The afflux at Section 4 may be considered as an increased reservoir level at an approximately constant distance from the structure. The downstream flow expansion occurs further downstream with increased afflux, as the structure control increasingly dominates the downstream control. In accord with the afflux method, the dimensionless contraction lengths, Lc/D3 and Le/D3, were tested against the HRC variables in Figures 1.9 and Figure 1.9 first illustrates an approximate linear correlation of Lc/D3 increasing with F3 alone and decreasing with J3 alone. The F3 trend was meaningful since it was expected that Lc approached zero as F3 approached zero (i.e. still water). And the J3 trend indicated that Lc approached zero as J3 approached 1 (i.e. the flow was totally blocked). A quadratic surface trendline was applied to the 3 dimensionless variables, and the coefficient of determination was reasonably high at The parity between the measured and predicted results was of the same magnitude. The results were similar for the Le/D3 analysis. Figure 1.10 illustrates an approximate linear correlation of Le/D3 increasing with F3 alone and decreasing with J3 alone. The F3 trend was meaningful since it was expected that Lc approached zero as F3 approached zero (i.e. still water). And the J3 trend indicated that Lc approached zero as J3 approached 1 (i.e. the flow was totally blocked). A quadratic surface trendline was applied to the 3 dimensionless variables, and the coefficient of determination was high at The parity between the measured and predicted results was of the same magnitude. 9
23 Lc: m y = -11x R 2 = dh: mm Figure 1-7: Plot of Contraction length (Lc) against Afflux (dh) with linear trend Le: m y = x R 2 = dh: mm Figure 1-8: Plot of Expansion length (Le) against Afflux (dh) with possible linear trend 10
24 Lc/D3 Lc/D3 10 Case 3, 4, 5 Case 2 Case y = x R 2 = F3 10 y = x R 2 = J y = 0.608x R 2 = Predicted Measured Figure 1-9: Variation of dimensionless contraction length with F3 and J3, and predicted trended length against measured length 11
25 y = 174.2x R 2 = Le/D F y = x R 2 = Le/D J y = x R 2 = Predicted Measured Figure 1-10: Variation of dimensionless expansion length with F3 and J3, and predicted trended length against measured length 12
26 1.5 Transition length analyses for the HEC (1995) data A 2DH computational procedure was used in HEC (1995) for estimating transition lengths. Firstly, 5 flood events from the USGS (1978) river afflux data were calibrated for use with the 2DH model. The associated variables and dimensions were then used to simulate a single, symmetrically shaped river model for various conditions of bridge opening width (b), discharge (Q), overbank and main channel (n = 0.04 throughout) Mannings coefficients and bedslopes (Figure 1.11). A total of 76 computations were conducted using spill-through abutments and 11 computations for vertical abutments. The variables used for the spill through abutments were converted to SI units and are enumerated in Appendices 1.2, 1.3 and 1.4. Afflux Advisor was then used to rapidly compute the associated HRC dimensionless variables, and these are also listed in the Appendices. An HRC analysis was then conducted as for the UB (2002) data above. The Lc/D3 linear correlation against F3 and D3 is illustrated in Figure Note that the F3 range is smaller than that used for the UB laboratory data. However, Lc/D3 increased with F3 and approached zero as F3 approached zero (i.e. to the still water condition). The overall Lc/D3 correlation with J3 is more scattered than that for the UB data. This is because the data consisted almost equally of 3 different bridge opening widths (b). A further analysis showed that trendlines through each of the separate bridge width data approached zero as J3 approached zero (i.e. to complete blockage). A quadratic surface trendline was applied to the 3 dimensionless variables, and the coefficient of determination was moderate at The parity between the measured and predicted results was of the same magnitude. Similar comments apply to the Le/D3 analysis (Figure 1.13), however the quadratic trending coefficient of determination was significantly improved at Figure 1-11: River cross section and data used in the HEC (1995) computational model length 13
27 45 40 y = x R 2 = Lc/D F y = x R 2 = Lc/D J y = x R 2 = Predicted Measured Figure 1-12: Variation of dimensionless contraction length with F3 and J3 and predicted trended length against measured length 14
28 y = x R 2 = Le/D F y = 53.41x R 2 = Le/D J y = x R 2 = Predicted Measured Figure 1-13: Variation of dimensionless expansion length with F3 and J3, and predicted trended length against measured length 15
29 1.6 Transition length analyses for combined UB (2002) and HEC (1995) data Although the separate UB and HEC data used different methods for estimating transition lengths, an attempt was made to combine all data using the HRC variables. This gave coefficients of determination of 0.28 for Lc/D3 and 0.66 for Le/D3. It was clear that the data required a fourth dimensionless variable. The obvious additional variable was that of flow scale, since small laboratory flows were being combined with large field conditions. Since both data sets used a single floodplain width (B = 1.2 m and 305 m respectively), the scale variable B/Bmean (B/Bm) was considered appropriate. The combined transition data are first illustrated with their respective quadratic trends in Figure The floodplain width scale is immediately apparent for the Lc/D3 plot. Although a correspondence of combined Le/D3 predictions with the HRC variables occurs, the curve is significantly extended with the low friction, laboratory results at Le/D3 = 100. The algorithm chosen to represent the combined data sets used a simple linear interpolation between the quadratic trend coefficients of each data set since there was no other available data. For Lc/D3, the algorithm was: b0 = * B/Bm b1 = * B/Bm b2 = * B/Bm b3 = * B/Bm b4 = * B/Bm b5 = * B/Bm Lc_D3 = b0 + b1 * F3 + b2 * J3 + b3 * F3 ^ 2 + b4 * J3 ^ 2 + b5 * F3 * J3 And for Le/D3, the algorithm was: b0 = * B/Bm b1 = * B/Bm b2 = * B/Bm b3 = * B/Bm b4 = * B/Bm b5 = * B/Bm Le_D3 = b0 + b1 * F3 + b2 * J3 + b3 * F3 ^ 2 + b4 * J3 ^ 2 + b5 * F3 * J3 The combined data are illustrated in Figure 1.15, and the respective coefficients of determination were 0.89 and 0.80 for contraction and expansion lengths respectively. The high correlations imply that either dimensionless transition length may be computed to be within about 4% of its mean value (at the 95% confidence level). These algorithms may therefore be applied to all floodplain widths between about 1 and 1000 m. Inevitably, there is a need for intermediate field data for further verification and analyses. 16
30 45 40 y = x R 2 = Predicted y = 0.608x R 2 = Measured y = x R 2 = Predicted y = x R 2 = Measured Figure 1-14: Predicted against measured dimensionless contraction lengths (upper curve) and expansion lengths (lower curve) for the UB (open symbols) and HEC data 17
31 45 40 y = x R 2 = Predicted Measured y = x R 2 = Predicted Measured Figure 1-15: Predicted against measured dimensionless contraction lengths (upper curve) and expansion lengths (lower curve) for combined UB and HEC data 18
32 1.7 Discussion of the HEC (1995) transition length analyses Contraction lengths The procedure used in HEC (1995) to determine transition lengths involved a multiple regression analysis of dimensionless flow variables (Tables 14 and 15, HEC, 1995). Since the quadratic trending of 3 variables is a better method than multiple regression for interpolation and extrapolation, the HEC (1995) results are reproduced herein and discussed using quadratic trending. Initially, the independent variables were considered to be the flow discharge (Q), the bridge opening width (b), the average stream bedslope (S) and the overbank Mannings coefficient (nob). They are enumerated in Appendices 1.2, 1.3 and 1.4. The following relation (which included constants of the computations) was proposed for the contraction ratio (CR): CR = Lc/Lobs = f(q/qm, nob/nmc, S) = f(d, N, S) where Lobs is the average obstruction length equal to 0.5*(B-b) for a beam bridge, Qm is the average discharge, and nmc is the constant main channel friction coefficient (nmc = 0.04). Table 15 from HEC (1995) is reproduced in Table 1.3, and mean CR values are added to amplify the data. A quadratic trend analysis was first made on the variables for constant S, then the trend coefficients were varied with S to give the following algorithm for CR: b0 = * S ^ * S b1 = * S ^ * S b2 = * S ^ * S b3 = * S ^ * S b4 = * S ^ * S b5 = * S ^ * S CR = b0 + b1 * D + b2 * N + b3 * D ^ 2 + b4 * N ^ 2 + b5 * D * N This CR algorithm was used to regenerate Table 1.2 as Table 1.3, and the small differences are enumerated in Table 1.4. Note that the differences average about 0.1 for the total data. Although this dimensionless representation is recommended by HEC (1995) for non-iterative use, it has the deficiency that Q is a function of S, and these variables are therefore not independent. Table 1-2: HEC Contraction Ratios (from HEC, 1995, Table 15) N=1 N=2 N=4 S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D= Table 1-3: Trended Contraction Ratios N=1 N=2 N=4 S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=
33 Table 1-4: Difference between HEC and trended Contraction Ratios N=1 N=2 N=4 S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D= Expansion lengths A similar analysis by HEC (1995) was conducted for expansion ratios, except that the added variable b/b was used to regress the data (Table 1.5), thus: ER = Le/Lobs = f(q/qbar, nob/nmc, S, b/b) = f(d, N, S, b/b) Since Lobs = 0.5*(B-b) = 0.5*B(1 b/b) and B was constant for the data, the addition of the b/b variable was unnecessary. The data were therefore trended as for the CR data, and gave the following algorithm: b0 = * S b1 = * S ^ * S b2 = * S ^ * S b3 = * S ^ * S b4 = * S ^ * S b5 = * S ER = b0 + b1 * D + b2 * N + b3 * D ^ 2 + b4 * N ^ 2 + b5 * D * N This ER algorithm was used to regenerate Table 1.5 as Table 1.6, and the small differences are enumerated in table 1.7. Note that the differences average about 0.2 for the total data. Furthermore, Table 1.6 is simpler to use. Table 1-5: HEC Expansion Ratios (From HEC 1995, Table 15) N=1 N=2 N=4 b/b S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 b/b = b/b = b/b = Table 1-6: Trended Expansion Ratios N=1 N=2 N=4 S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D= Table 1-7: Difference between HEC and trended Expansion Ratios 20
34 N=1 N=2 N=4 b/b S D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 D=0.3 D=1.0 D=1.7 b/b = b/b = b/b = Design tables for transition lengths The HRC transition analyses In general, the HRC analyses accurately predict the transition lengths as: Lc/D3 or Le/D3 = f(j3, F3, B/Bm) However, the use of the flow scale (B/Bm) is considered herein as a mathematical convenience, since the general friction conditions for both data sets varied. For the UB (2002) data, the average main channel Mannings friction coefficient (nmc) was and the floodplain coefficient (nfp) was For the HEC (1995) data, the average nmc was and the average nfp was It was thus inevitable that the expansion ratio (Le/Lobs) was larger for the lower friction conditions of the UB data, since the decreased friction gave lower lateral turbulent viscosity coefficients. Further data coupled with computer analyses are required to improve the HRC physical model towards using a friction condition rather than the B/Bm variable. Meanwhile, the existing analyses are proposed for present design Simplified transition length design It was noted above that the UB (2002) contraction lengths (Lc) were near constant, with a 95% confidence limit (CL) of about 5% of the average magnitude. Similar average statistics were applied to both the UB and HEC data (Table 1.8). The CL for Lc were similar at 5%, and for Le they were about 10% of the absolute magnitudes. A simplified transition length, depending upon floodplain width (B) only, may thus be computed from these data using the prediction equations in Table 1.8. Note that these values are approximate, and a more detailed design is recommended below. Table 1-8: Simplified design for transition lengths B Lc Le Simple prediction Lc = 0.38B Le = 0.55B % CL 5% 10% Detailed transition length design The HRC algorithms for dimensionless transition length design (Lc/D3 and Le/D3) are reproduced for floodplain widths (B) of 1 m, 10 m, 100 m and 1000 m in Tables 1.10 and The scale of B = 1 m applies to laboratory flows, and that of B = 1000 m applies to flooded estuaries for the UK. Intermediate values may be obtained by linear interpolation, since a linear scale was used for the algorithm. For design cases where normal flow data are not readily available (such as that for the HEC, 1995 data), Afflux Advisor may be used to rapidly compute such flows from the Mannings friction conditions. The extrapolation of quadratic surface trends may sometimes lead to negative values at the extremities. Any tabulated negative data has simply been replaced herein with zeroes, since that is 21
35 the data trend. For the B = 1000 m expansion length table, most data for FR > 0.3 was negative, and this occurred because the HEC computations did not extend beyond this Froude number. Inevitably, high Froude number flows for such wide and deep floodplains lead to unrealistically extreme flows. It is finally noted that the HRC transition algorithms, or Tables 1.9 and 1.10, may be used without further detailed iteration such as required by the HEC (1995) method. Furthermore, the different laboratory and field transition lengths necessitate the use of a HEC-RAS model to determine expansion and contraction energy loss coefficients for all of the UB experiments. As a first approach, it is expected that the UB contraction energy loss coefficient (Cc) will be similar to that of the HEC (1995) value of Cc = 0.1, since Lc/Lobs is similar. Initial calibration may therefore begin using this value, and the expansion energy loss coefficient adjusted according to the HRC expansion length and associated afflux. Table 1-9: Simplified design for transition lengths B = 1.0 m B = 10 m J3 F J3 F B = 100 m B = 1000 m J3 F J3 F
36 Table 1-10: Detailed design for dimensionless expansion lengths (Le/D3) B = 1.0 m B = 10 m J3 F J3 F B = 100 m B = 1000 m J3 F J3 F
37 1.9 Appendices to Chapter 1 Appendix 1.1: Contraction and expansion lengths (Lc, Le) for the UB (2002) normal bridge flow data with HRC dimensionless variables C a s e Experiment Discharge F3 J3 dh/d3 Lc Le Lc/D3 Le/D3 Name (l/s) m m 1 AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC AMOSC ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOE ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC ASOSC D398121P b=398 D398124P with piers D398127P D398130P D398135P D b=398mm D no piers D D D D398321P b=398mm D398324P with piers D398327P D398330P D398335P D b=498mm D no piers D D D D498321P b=498mm D498324P with piers D498327P
BACKWATERRISE DUE TO FLOW CONSTRICTION BY BRIDGE PIERS
Thirteenth International Water Technology Conference, IWTC 1 009, Hurghada, Egypt BACKWATERRISE DUE TO FLOW CONSTRICTION BY BRIDGE PIERS Kassem Salah El-Alfy Prof. Dr., Irrigation &Hydraulics Dept., Faculty
EXPERIMENTAL STUDY OF BACKWATER RISE DUE TO BRIDGE PIERS AS FLOW OBSTRUCTIONS
Tenth International Water Technology Conference, IWTC1 6, Alexandria, Egypt 19 EXPERIMENTAL STUDY OF BACKWATER RISE DUE TO BRIDGE PIERS AS FLOW OBSTRUCTIONS Kassem Salah El-Alfy Associate Prof., Irrigation
Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
OPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK
OPEN CHANNEL FLOW Numerical Methods and Computer Applications Roland Jeppson TECHNISCHE INFORMATlONSBiBUOTHEK UNIVERSITATSB'BUOTHEK HANNOVER Si. i. CRC Press Taylor &.Francis Group Boca Raton London New
THE EFFECTS OF OBSTACLES ON SURFACE LEVELS AND BOUNDARY RESISTANCE IN OPEN CHANNELS
Manuscript submitted to 0th IAHR Congress, Thessaloniki, 4-9 August 00 THE EFFECTS OF OBSTACLES ON SURFACE LEVELS AND BOUNDARY RESISTANCE IN OPEN CHANNELS J. D. FENTON Department of Civil and Environmental
Closed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
UPPER COSUMNES RIVER FLOOD MAPPING
UPPER COSUMNES RIVER FLOOD MAPPING DRAFT BASIC DATA NARRATIVE FLOOD INSURANCE STUDY SACRAMENTO COUTY, CALIFORNIA Community No. 060262 November 2008 Prepared By: CIVIL ENGINEERING SOLUTIONS, INC. 1325 Howe
Evaluation of Different Software Packages in Flow Modeling under Bridge Structures
4*6 $ " $D )>2A " B*C )(=(@$1 )'2/ &35 &/( L!1' 2 (// QR4 /d/d G*) T)!MIKE11 54 /'# 4 G(%?,J # ( A# # /; ( 7;# F K u S E# 5#/ 7#%5. # 5"# " D#P 7"#B4
We will assume straight channels with simple geometries (prismatic channels) and steady state flow (in time).
56 Review Drag & Lift Laminar vs Turbulent Boundary Layer Turbulent boundary layers stay attached to bodies longer Narrower wake! Lower pressure drag! 8. Open-Channel Flow Pipe/duct flow closed, full,
Dealing with Zone A Flood Zones. Topics of Discussion. What is a Zone A Floodplain?
Dealing with Zone A Flood Zones Topics of Discussion Overview of Zone A Floodplains Permitting Development in Zone A Floodplains Estimating Flood Elevations in Zone A Flood Insurance Implications Letters
Hydromechanics: Course Summary
Hydromechanics: Course Summary Hydromechanics VVR090 Material Included; French: Chapters to 9 and 4 + Sample problems Vennard & Street: Chapters 8 + 3, and (part of it) Roberson & Crowe: Chapter Collection
Prediction of bed form height in straight and meandering compound channels
Water Resources Management III 311 Prediction of bed form height in straight and meandering compound channels R. D. Karamisheva, J. F. Lyness, W. R. C. Myers, J. O Sullivan & J. B. C. Cassells School of
Lateral Inflow into High-Velocity Channels
Lateral Inflow into High-Velocity Channels by Richard L. Stockstill PURPOSE: This Coastal and Hydraulics Engineering Technical Note (CHETN) investigates lateral flow discharging into a high-velocity channel.
CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - II GRADUALLY VARIED FLOW
CE 6403 APPLIED HYDRAULIC ENGINEERING UNIT - II GRADUALLY VARIED FLOW Dynamic equations of gradually varied and spatially varied flows - Water surface flow profile classifications: Hydraulic Slope, Hydraulic
Open Channel Flow Part 2. Ch 10 Young, notes, handouts
Open Channel Flow Part 2 Ch 10 Young, notes, handouts Uniform Channel Flow Many situations have a good approximation d(v,y,q)/dx=0 Uniform flow Look at extended Bernoulli equation Friction slope exactly
Hydraulics Part: Open Channel Flow
Hydraulics Part: Open Channel Flow Tutorial solutions -by Dr. K.N. Dulal Uniform flow 1. Show that discharge through a channel with steady flow is given by where A 1 and A 2 are the sectional areas of
HEC-RAS River Analysis System
US Army Corps of Engineers Hydrologic Engineering Center HEC-RAS River Analysis System Hydraulic Reference Manual Version 3.1 November 00 Approved for Public Release. Distribution Unlimited CPD-69 Table
y 2 = 1 + y 1 This is known as the broad-crested weir which is characterized by:
CEE 10 Open Channel Flow, Dec. 1, 010 18 8.16 Review Flow through a contraction Critical and choked flows The hydraulic jump conservation of linear momentum y = 1 + y 1 1 + 8Fr 1 8.17 Rapidly Varied Flows
LOMR SUBMITTAL LOWER NEHALEM RIVER TILLAMOOK COUNTY, OREGON
LOMR SUBMITTAL LOWER NEHALEM RIVER TILLAMOOK COUNTY, OREGON Prepared for: TILLAMOOK COUNTY DEPARTMENT OF COMMUNITY DEVELOPMENT 1510-B THIRD STREET TILLAMOOK, OR 97141 Prepared by: 10300 SW GREENBURG ROAD,
Two-Dimensional Simulation of Truckee River Hydrodynamics
Two-Dimensional Simulation of Truckee River Hydrodynamics by Stephen H. Scott PURPOSE: The purpose of this Coastal and Hydraulics Engineering Technical Note (CHETN) is to demonstrate the use of multidimensional
1.060 Engineering Mechanics II Spring Problem Set 8
1.060 Engineering Mechanics II Spring 2006 Due on Monday, May 1st Problem Set 8 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members
Accounting for increased flow resistance due to lateral momentum loss in restoration designs using 2-stage channels
Skamania 2005 Accounting for increased flow resistance due to lateral momentum loss in restoration designs using 2-stage channels Outline Aim and Objectives Definition Use of 2-stage channels in stream
CHAPTER 2- BACKGROUND. INVESTIGATIONS OF COMPOSITE ROUGHNESS COEFFICIENT IN A RIVER WITH LOW FLOW
2. Background 2.1 Introduction The estimation of resistant coefficient and hence discharge capacity in a channel or river is one of the fundamental problems facing river engineers. When applying Manning
MODELING FLUID FLOW IN OPEN CHANNEL WITH HORSESHOE CROSS SECTION
July. 2. Vol. 7. No. 2 MODELING FLUID FLOW IN OPEN CHANNEL WITH HORSESHOE CROSS SECTION 1 J. JOMBA, 2 D.M.THEURI, 2 E. MWENDA, 2 C. CHOMBA ABSTRACT Flow in a closed conduit is regarded as open channel
conservation of linear momentum 1+8Fr = 1+ Sufficiently short that energy loss due to channel friction is negligible h L = 0 Bernoulli s equation.
174 Review Flow through a contraction Critical and choked flows The hydraulic jump conservation of linear momentum y y 1 = 1+ 1+8Fr 1 8.1 Rapidly Varied Flows Weirs 8.1.1 Broad-Crested Weir Consider the
LOMR SUBMITTAL LOWER NESTUCCA RIVER TILLAMOOK COUNTY, OREGON
LOMR SUBMITTAL LOWER NESTUCCA RIVER TILLAMOOK COUNTY, OREGON Prepared for: TILLAMOOK COUNTY DEPARTMENT OF COMMUNITY DEVELOPMENT 1510-B THIRD STREET TILLAMOOK, OR 97141 Prepared by: 10300 SW GREENBURG ROAD,
CEE 3310 Open Channel Flow, Nov. 26,
CEE 3310 Open Channel Flow, Nov. 6, 018 175 8.10 Review Open Channel Flow Gravity friction balance. y Uniform Flow x = 0 z = S 0L = h f y Rapidly Varied Flow x 1 y Gradually Varied Flow x 1 In general
MODELING OF LOCAL SCOUR AROUND AL-KUFA BRIDGE PIERS Saleh I. Khassaf, Saja Sadeq Shakir
ISSN 2320-9100 11 International Journal of Advance Research, IJOAR.org Volume 1, Issue 8,August 2013, Online: ISSN 2320-9100 MODELING OF LOCAL SCOUR AROUND AL-KUFA BRIDGE PIERS Saleh I. Khassaf, Saja Sadeq
Evaluation of Scour Depth around Bridge Piers with Various Geometrical Shapes
Evaluation of Scour Depth around Bridge Piers with Various Geometrical Shapes Dr. P. D. Dahe * Department of Civil Engineering, SGGSIE&T, Vishnupuri, Nanded (Maharashtra) S. B. Kharode Department of Civil
Conveyance Estimation System (CES) Launch
(CES) Launch Conveyance Principles Caroline McGahey 10 June 2004 Why is there a new approach? Existing 1D software still based on historic hand-calculation methods - not based in rigorous physics - as
Discharge estimation in compound channels with fixed and mobile bed
Sādhanā Vol. 34, Part 6, December 2009, pp. 923 945. Indian Academy of Sciences Discharge estimation in compound channels with fixed and mobile bed GALIP SECKIN 1, MUSTAFA MAMAK 1, SERTER ATABAY 2 and
Dr. Muhammad Ali Shamim ; Internal 652
Dr. Muhammad Ali Shamim ali.shamim@uettaxila.edu.pk 051-904765; Internal 65 Channel Tranistions A channel transition is defined as change in channel cross section e.g. change in channel width and/or channel
Open Channel Hydraulics III - Sharpcrested
PDHonline Course H140 (2 PDH) Open Channel Hydraulics III - Sharpcrested Weirs Instructor: Harlan H. Bengtson, Ph.D., PE 2012 PDH Online PDH Center 5272 Meadow Estates Drive Fairfax, VA 22030-6658 Phone
NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS
TASK QUARTERLY 15 No 3 4, 271 282 NUMERICAL SIMULATION OF OPEN CHANNEL FLOW BETWEEN BRIDGE PIERS MICHAŁ SZYDŁOWSKI Faculty of Civil and Environmental Engineering, Gdansk University of Technology, Narutowicza
Lecture 10: River Channels
GEOG415 Lecture 10: River Channels 10-1 Importance of channel characteristics Prediction of flow was the sole purpose of hydrology, and still is a very important aspect of hydrology. - Water balance gives
Representing Bridges in or near TUFLOW
Representing Bridges in or near TUFLOW What to represent? Soffit Deck Piers Channel shape changes Associated structures (weirs, sluices, etc.) What are our options? Pure TUFLOW representations Energy Losses
APPENDIX B Hydraulic Considerations for Pipeline Crossings of Stream Channels
APPENDIX B Hydraulic Considerations for Pipeline Crossings of Stream Channels B-1 B-2 APPENDIX B HYDRAULIC CONSIDERATIONS FOR PIPELINE CROSSINGS OF STREAM CHANNELS Pipeline crossings of perennial, intermittent,
PENNSYLVANIA DEPARTMENT OF TRANSPORTATION ENGINEERING DISTRICT 3-0
PENNSYLVANIA DEPARTMENT OF TRANSPORTATION ENGINEERING DISTRICT 3-0 LYCOMING COUNTY S.R.15, SECTION C41 FINAL HYDROLOGIC AND HYDRAULIC REPORT STEAM VALLEY RUN STREAM RELOCATION DATE: June, 2006 REVISED:
9. Flood Routing. chapter Two
9. Flood Routing Flow routing is a mathematical procedure for predicting the changing magnitude, speed, and shape of a flood wave as a function of time at one or more points along a watercourse (waterway
Open Channel Flow I - The Manning Equation and Uniform Flow COURSE CONTENT
Open Channel Flow I - The Manning Equation and Uniform Flow Harlan H. Bengtson, PhD, P.E. COURSE CONTENT 1. Introduction Flow of a liquid may take place either as open channel flow or pressure flow. Pressure
EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018
EXAMPLES (SEDIMENT TRANSPORT) AUTUMN 2018 Q1. Using Cheng s formula estimate the settling velocity of a sand particle of diameter 1 mm in: (a) air; (b) water. Q2. Find the critical Shields parameter diameter
* Ho h h (3) D where H o is the water depth of undisturbed flow, D is the thickness of the bridge deck, and h is the distance from the channel floor t
The Seventh International Colloquium on Bluff Body Aerodynamics and Applications (BBAA7) Shanghai, China; September -6, 01 Numerical simulation of hydrodynamic loading on submerged rectangular bridge decks
Cavitation occurs whenever the pressure in the flow of water drops to the value of the pressure of the saturated water vapour, pv (at the prevailing
Cavitation occurs whenever the pressure in the flow of water drops to the value of the pressure of the saturated water vapour, pv (at the prevailing temperature); cavities filled by vapour, and partly
Section 4: Model Development and Application
Section 4: Model Development and Application The hydrologic model for the Wissahickon Act 167 study was built using GIS layers of land use, hydrologic soil groups, terrain and orthophotography. Within
Beaver Creek Corridor Design and Analysis. By: Alex Previte
Beaver Creek Corridor Design and Analysis By: Alex Previte Overview Introduction Key concepts Model Development Design Accuracy Conclusion Refresh v = Beaver Creek Site = Wittenberg Introduction Low head
Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati
Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati Module - 5 Channel Transitions Lecture - 1 Channel Transitions Part 1 Welcome back
The Influence of Piers on the Risk of Flooding Upstream from a Bridge
Archives of Hydro-Engineering and Environmental Mechanics Vol. 52 (2005), No. 4, pp. 287 301 IBW PAN, ISSN 1231 3726 The Influence of Piers on the Risk of Flooding Upstream from a Bridge Marek Sowiński
P10.5 Water flows down a rectangular channel that is 4 ft wide and 3 ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow.
P10.5 Water flows down a rectangular channel that is 4 ft wide and ft deep. The flow rate is 15,000 gal/min. Estimate the Froude number of the flow. Solution: Convert the flow rate from 15,000 gal/min
Q = α n AR2/3 h S1/2 0. bd 2d + b d if b d. 0 = ft1/3 /s
CEE 330 Open Channel Flow, Nov., 00 7 8.8 Review Open Channel Flow Gravity friction balance. In general we take an energy equation approach. y Uniform Flow x = 0 z = S 0L = h f where we could find h f
Presented by: Civil Engineering Academy
Presented by: Civil Engineering Academy Open-Channel Flow Uniform Flow (See CERM Ch. 19) Characterized by constant depth volume, and cross section. It can be steady or unsteady Non-uniform Flow *Not on
Opanuku Stream Benchmark Validation 1. Introduction. 2. The Opanuku Stream Model
Opanuku Stream Benchmark Validation 1. Introduction The model accuracy benchmark published by the Flood Risk Management Committee of the IAHR in http://members.iahr.org/imis/communitymanagement/communitylayouts/flood_risk_manageme
Chapter 8: Flow in Pipes
8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump
A note on critical flow section in collector channels
Sādhan ā, Vol. 26, Part 5, October 2001, pp. 439 445. Printed in India A note on critical flow section in collector channels 1. Introduction SUBHASISH DEY Department of Civil Engineering, Indian Institute
JOURNAL OF ENVIRONMENTAL HYDROLOGY The Electronic Journal of the International Association for Environmental Hydrology VOLUME
JOURNAL OF ENVIRONMENTAL HYDROLOGY The Electronic Journal of the International Association for Environmental Hydrology VOLUME 18 2010 REDUCED CHANNEL CONVEYANCE ON THE WICHITA RIVER AT WICHITA FALLS, TEXAS,
Gradually Varied Flow I+II. Hydromechanics VVR090
Gradually Varied Flow I+II Hydromechanics VVR090 Gradually Varied Flow Depth of flow varies with longitudinal distance. Occurs upstream and downstream control sections. Governing equation: dy dx So Sf
OPEN CHANNEL FLOW. One-dimensional - neglect vertical and lateral variations in velocity. In other words, Q v = (1) A. Figure 1. One-dimensional Flow
OPEN CHANNEL FLOW Page 1 OPEN CHANNEL FLOW Open Channel Flow (OCF) is flow with one boundary exposed to atmospheric pressure. The flow is not pressurized and occurs because of gravity. Flow Classification
NPTEL Quiz Hydraulics
Introduction NPTEL Quiz Hydraulics 1. An ideal fluid is a. One which obeys Newton s law of viscosity b. Frictionless and incompressible c. Very viscous d. Frictionless and compressible 2. The unit of kinematic
Summary of Hydraulic and Sediment-transport. Analysis of Residual Sediment: Alternatives for the San Clemente Dam Removal/Retrofit Project,
Appendix N SUMMARY OF HYDRAULIC AND SEDIMENT-TRANSPORT ANALYSIS OF RESIDUAL SEDIMENT: ALTERNATIVES FOR THE SAN CLEMENTE DAM REMOVAL/RETROFIT PROJECT, CALIFORNIA the San Clemente Dam Removal/Retrofit Project,
Engineering Hydrology (ECIV 4323) CHAPTER FOUR. Stream flow measurement. Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib
Engineering Hydrology (ECIV 4323) CHAPTER FOUR Stream flow measurement Instructors: Dr. Yunes Mogheir Dr. Ramadan Al Khatib -١ 4.1 Introduction - Surface water hydrology deals with the movement of water
Chapter 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
Lab 7: Nonuniform Flow and Open Channel Transitions
CE 3620: Water Resources Engineering Spring 2015 Lab 7: Nonuniform Flow and Open Channel Transitions BACKGROUND An open channel transition may be defined as a change either in the direction, slope, or
Geology 550 Spring 2005 LAB 3: HYDRAULICS OF PRAIRIE CREEK
Geology 550 Spring 2005 LAB 3: HYDRAULICS OF PRAIRIE CREEK Objectives: 1. To examine the distribution of velocity in a stream channel 2. To characterize the state of flow using dimensionless variables
Influence of Two-line Emergent Floodplain Vegetation on A Straight Compound Channel Flow
International Journal of Integrated Engineering, Vol. 5 No. 1 (2013) p. 58-63 Influence of Two-line Emergent Floodplain Vegetation on A Straight Compound Channel Flow Mazlin Jumain 1,*, Zulkiflee Ibrahim
Applying GIS to Hydraulic Analysis
Texas A&M University Department of Civil Engineering CVEN689 Applications of GIS to Civil Engineering Instructor: Francisco Olivera, Ph.D., P.E. Applying GIS to Hydraulic Analysis Lim, Chae Kwan April
FLOOD RISK MAPPING AND ANALYSIS OF THE M ZAB VALLEY, ALGERIA
River Basin Management IX 69 FLOOD RISK MAPPING AND ANALYSIS OF THE M ZAB VALLEY, ALGERIA AMEL OUCHERIF & SAADIA BENMAMAR National Polytechnic School of Algiers, Algeria ABSTRACT To contribute to flood
Pequabuck River Flooding Study and Flood Mitigation Plan The City of Bristol and Towns of Plainville and Plymouth, CT
Pequabuck River Flooding Study and Flood Mitigation Plan The City of Bristol and Towns of Plainville and Plymouth, CT Raymond Rogozinski and Maged Aboelata The City of Bristol and Towns of Plainville and
A STUDY OF LOCAL SCOUR AT BRIDGE PIERS OF EL-MINIA
A STUDY OF LOCAL SCOUR AT BRIDGE PIERS OF EL-MINIA Dr. Gamal A. Sallam 1 and Dr. Medhat Aziz 2 ABSTRACT Bridges are critical structures that require a substantial investment to construct and serve an important
YELLOWSTONE RIVER FLOOD STUDY REPORT TEXT
YELLOWSTONE RIVER FLOOD STUDY REPORT TEXT TECHNICAL REPORT Prepared for: City of Livingston 411 East Callender Livingston, MT 59047 Prepared by: Clear Creek Hydrology, Inc. 1627 West Main Street, #294
Design Hydraulic Study. Bridge 09C-0134, Blairsden-Graeagle Road over Middle Fork Feather River. Plumas County. Prepared for:
Design Hydraulic Study Bridge 09C-0134, Blairsden-Graeagle Road over Middle Fork Feather River Plumas County Prepared for: Quincy Engineering, Inc 3247 Ramos Circle Sacramento, CA 95827-2501 Prepared by:
Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow
Archives of Hydro-Engineering and Environmental Mechanics Vol. 61 (2014), No. 3 4, pp. 89 109 DOI: 10.1515/heem-2015-0006 IBW PAN, ISSN 1231 3726 Comparison of Average Energy Slope Estimation Formulas
Energy and momentum coefficients for wide compound channels
River Basin Management VII 87 Energy and momentum coefficients for wide compound channels P. K. Mohanty, S. S. Dash, K. K. Khatua & K. C. Patra Department of Civil Engineering, N.I.T. Rourkela, India Abstract
Basic Hydraulics. Rabi H. Mohtar ABE 325
Basic Hydraulics Rabi H. Mohtar ABE 35 The river continues on its way to the sea, broken the wheel of the mill or not. Khalil Gibran The forces on moving body of fluid mass are:. Inertial due to mass (ρ
4. BACKWATER EFFECTS OF PIERS IN SUBCRITICAL FLOW 5. Report Date October Performing Organization Code
Technical Report Documentation Page 1. Report No. FHWA/TX-0-1805-1. Government Accession No. 3. Recipient s Catalog No. 4. BACKWATER EFFECTS OF PIERS IN SUBCRITICAL FLOW 5. Report Date October 001 7. Author(s)
CIE4491 Lecture. Hydraulic design
CIE4491 Lecture. Hydraulic design Marie-claire ten Veldhuis 19-9-013 Delft University of Technology Challenge the future Hydraulic design of urban stormwater systems Focus on sewer pipes Pressurized and
HEC-RAS MODELING OF RAINBOW RIVER MFL TECHNICAL SUPPORT FRESHWATER STREAM FINAL REPORT SWFWMD TWA# 15TW
Appendix I HEC-RAS Modeling of Rainbow River, MFL Technical Support Freshwater Stream. Prepared by Environmental Consulting & Technology, Inc. for the Southwest Florida Water Management District, Brooksville,
Rating Curves: Part 1 Correction for Surface Slope
The Institution of Engineers, Australia Conference on Hydraulics in Civil Engineering Hobart 8 30 November 001, pp 309-317 Rating Curves: Part 1 Correction for Surface Slope J. D. Fenton Dip.C.E., B.E.,
Hydrodynamic model of St. Clair River with Telemac-2D Phase 2 report
Hydrodynamic model of St. Clair River with Telemac-2D Phase 2 report Controlled Technical Report CHC-CTR-084 revision 1 March 2009 NRC-CHC has prepared this report for the International Joint Commission
Comparing HEC-RAS v5.0 2-D Results with Verification Datasets
Comparing HEC-RAS v5.0 2-D Results with Verification Datasets Tom Molls 1, Gary Brunner 2, & Alejandro Sanchez 2 1. David Ford Consulting Engineers, Inc., Sacramento, CA 2. USACE Hydrologic Engineering
The most important equation to describe the water balance of a reservoir is the water balance: Equation 3.1
3 FLOOD PROPAGATION 3.1 Reservoir routing The most important equation to describe the water balance of a reservoir is the water balance: ds = I Q + A P E dt ( ) In finite differences form this equation
DRAFT Design Hydraulic Study. Bridge 04C-0055, Mattole Road Bridge over Mattole River at Honeydew. Humboldt County. Prepared for:
DRAFT Design Hydraulic Study Bridge 04C-0055, Mattole Road Bridge over Mattole River at Honeydew Humboldt County Prepared for: Morrison Structures Incorporated 1890 Park Marina Drive, Suite 104 Redding,
Flood Map. National Dataset User Guide
Flood Map National Dataset User Guide Version 1.1.5 20 th April 2006 Copyright Environment Agency 1 Contents 1.0 Record of amendment... 3 2.0 Introduction... 4 2.1 Description of the Flood Map datasets...4
Estimating of Manning s Roughness Coefficient for Hilla River through Calibration Using HEC-RAS Model
Estimating of Manning s Roughness Coefficient for Hilla River through Calibration Using HEC-RAS Model Luay Kadhim Hameed 1) and Salah Tawfeek Ali 2) 1) University of Kufa, Kufa, Iraq 2) University of Babylon,
Flood and runoff estimation on small catchments. Duncan Faulkner
Flood and runoff estimation on small catchments Duncan Faulkner Flood and runoff estimation on small catchments Duncan Faulkner and Oliver Francis, Rob Lamb (JBA) Thomas Kjeldsen, Lisa Stewart, John Packman
ESTIMATING JOINT FLOW PROBABILITIES AT STREAM CONFLUENCES USING COPULAS
ESTIMATING JOINT FLOW PROBABILITIES AT STREAM CONFLUENCES USING COPULAS Roger T. Kilgore, P.E., D. WRE* Principal Kilgore Consulting and Management 2963 Ash Street Denver, CO 80207 303-333-1408 David B.
B-1. Attachment B-1. Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling
Attachment B-1 Evaluation of AdH Model Simplifications in Conowingo Reservoir Sediment Transport Modeling 1 October 2012 Lower Susquehanna River Watershed Assessment Evaluation of AdH Model Simplifications
Keywords: flow characteristics, compound straight channel, bed morphology, floodplain
Flow Characteristics on Floodplain Vegetation in Compound Straight Channels Nur Amirah Nabilah Mohd Zamri 1, a, Zulhilmi Ismail 1,b,Zulkiflee Ibrahim 1,c 1 Faculty of Civil Engineering, Universiti Teknologi
Coastal and Hydraulics Laboratory
SAM Hydraulic Design Package for Channels Coastal and Hydraulics Laboratory William A. Thomas, Ronald R. Copeland, and Dinah N. McComas September 2002 Approved for public release; distribution is unlimited.
Chapter 4: Non uniform flow in open channels
Chapter 4: Non uniform flow in open channels Learning outcomes By the end of this lesson, students should be able to: Relate the concept of specific energy and momentum equations in the effect of change
Open Channel Hydraulics I - Uniform Flow
PDHonline Course H138 (2 PDH) Open Channel Hydraulics I - Uniform Flow Instructor: Harlan H. Bengtson, Ph.D., PE 2012 PDH Online PDH Center 5272 Meadow Estates Drive Fairfax, VA 22030-6658 Phone & Fax:
Calibration of Manning s Friction Factor for Rivers in Iraq Using Hydraulic Model (Al-Kufa River as Case study)
Calibration of Manning s Friction Factor for Rivers in Iraq Using Hydraulic Model (Al-Kufa River as Case study) Luay Kadhim Hameed, Civil Engineering Dept./ University of Kufa Hayder Sami Mohammed, Structure
WINFAP 4 QMED Linking equation
WINFAP 4 QMED Linking equation WINFAP 4 QMED Linking equation Wallingford HydroSolutions Ltd 2016. All rights reserved. This report has been produced in accordance with the WHS Quality & Environmental
Geomorphological Modelling in Coastal Waters
Abstract Geomorphological Modelling in Coastal Waters Morteza Kolahdoozan 1, Roger A. Falconer 2 (Fellow), Yiping Chen 3 Details are given herein of the development and application of a three dimensional
River Ice Modelling. Hydraulic Modelling of Mackenzie River at Ft. Province, NWT, Canada
River Ice Modelling Hydraulic Modelling of Mackenzie River at Ft. Province, NWT, Canada Submitted to F.E. Hicks, PhD, PEng, FCSCE Professor, Dept. of Civil and Environmental Engineering University of Alberta
Hydraulics for Urban Storm Drainage
Urban Hydraulics Hydraulics for Urban Storm Drainage Learning objectives: understanding of basic concepts of fluid flow and how to analyze conduit flows, free surface flows. to analyze, hydrostatic pressure
Kinetic energy and momentum correction factors in a stream
Kinetic energy and momentum correction factors in a stream Mehmet Ardıçlıoğlu, Mücella Vapur, 2Onur Genç, University of Erciyes, Department of Civil Engineering, Kayseri, Turkey 2 Melikşah University,
Analysis of the Cause for the Collapse of a Temporary Bridge Using Numerical Simulation
Engineering, 2013, 5, 997-1005 Published Online December 2013 (http://www.scirp.org/journal/eng) http://dx.doi.org/10.4236/eng.2013.512121 Analysis of the Cause for the Collapse of a Temporary Bridge Using
Review 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
Open-channel hydraulics
Open-channel hydraulics STEADY FLOW IN OPEN CHANNELS constant discharge, other geometric and flow characteristics depended only on position Uniform flow Non-uniform flow S; y; v const. i i 0 i E y 1 y
Fish Passage at Road Crossings
Fish Passage at Road Crossings 1 Crossing Design Workshop Outline 1:00 to 2:00 Intro, Design Overview, Channel Width 2:00 to 2:15 Break 2:15 to 3:15 No-Slope, Stream Simulation Design 3:15 to 3:30 Break
EQUATIONS FOR DISCHARGE CALCULATION IN COMPOUND CHANNELS HAVING HOMOGENEOUS ROUGHNESS * S. M. HOSSEINI **
Iranian Journal of Science & Technology, Transaction B, Vol. 28, No. B5 Printed in The Islamic Republic of Iran, 2004 Shiraz University EQUATIONS FOR DISCHARGE CALCULATION IN COMPOUND CHANNELS HAVING HOMOGENEOUS