---

Created: 2026/10/21 10:30 Course: EWCE 461 - River Engineering and Flood Management Section: Section A (SA) - River Engineering Source: Lecture Materials (Lectures 1 to 12)

1. Thematic Curriculum Mapping
2. Detailed Hierarchical Outline
   1. Lecture 1: Rivers, Longitudinal Zones, Drainage Systems, and Channel Classification
      1. 1. Fundamentals of River Systems
      2. 2. Longitudinal Zonation of Rivers (Schumm's Basin Model)
      3. 3. Fluvial Morphology Terminology and Riparian Features
      4. 4. Drainage Network Patterns and Stream Topology
      5. 5. Multi-Criteria River Classification Schemes
      6. 6. Quantitative Planform Metrics
   2. Lecture 2: Channel Patterns, Governing Factors, Hydraulic Geometry, and River Longitudinal Profiles
      1. 1. Planform Patterns and Fluvial Dynamics
      2. 2. Governing Factors of Alluvial Channel Morphology
      3. 3. Hydrological Flood Response Classification
      4. 4. Theory of Hydraulic Geometry (Leopold & Maddock)
      5. 5. Channel-Forming and Dominant Discharges
      6. 6. River Longitudinal Profiles and Downstream Fining
   3. Lecture 3: Derived Downstream Hydraulic Relations, Meander Theories, and Bend Hydrodynamics
      1. 1. Derived Downstream Hydraulic Geometry Relationships
      2. 2. Planform Stability Boundaries and Thresholds
      3. 3. Hypotheses and Theories of Meander Formation
      4. 4. Morphometric Scaling of Meander Geometry
      5. 5. Bend Hydrodynamics and Helical Secondary Currents
   4. Lecture 4: Channel Stability Classification, Braiding Mechanics, and Sediment Load Modes
      1. 1. Mechanics and Genesis of River Braiding
      2. 2. Channel Stability Classification (Schumm's Fluvial Continuum)
      3. 3. Planform Stability (PS) Force Ratios
      4. 4. Comprehensive Sediment Load Classification Systems
   5. Lecture 5: River Hydrodynamics, Threshold of Movement, and the Shields Parameter
      1. 1. Hydrodynamic Forces and Boundary Shear Stress
      2. 2. Sediment Transport Concepts: Competence vs. Capacity
      3. 3. Mechanics of Incipient Motion (Threshold of Movement)
      4. 4. The Shields Parameter and Incipient Motion Criteria
      5. 5. Engineering Applications and Numerical Calculations
   6. Lecture 6: Flow Regimes, Bed Configurations, Flow Resistance, and Particle Settling Mechanics
      1. 1. Flow Regimes in Alluvial Sand-Bed Channels (Simons & Richardson)
      2. 2. Alluvial Flow Resistance and Roughness Partitioning
      3. 3. Particle Mechanics: Settling, Terminal, and Fall Velocity ($\omega$)
   7. Lecture 7: Quantitative Estimation of Sediment Transport Rates
      1. 1. Bed Shear Stress and Shear Velocity Foundations
      2. 2. Suspended Sediment Transport and the Rouse Profile
      3. 3. Bed Load Transport Formulations
      4. 4. Total Bed Material Load Prediction Methods
   8. Lecture 8: Dynamic Adjustment of Rivers, Lane’s Relationship, Aggradation, and Degradation
      1. 1. Emory Lane’s Fluvial Equilibrium Relationship
      2. 2. Quantitative Formulation of Lane’s Relation
      3. 3. Channel Aggradation and Degradation Processes
      4. 4. Mathematical Modeling & Field Case Studies
   9. Lecture 9: River Scour Fundamentals, Bridge Hydraulics, and Local Scour Mechanisms
      1. 1. Overview and Additive Components of Bridge Scour
      2. 2. Bridge Hydraulics and Contraction Scour
      3. 3. Mechanics of Local Pier Scour
      4. 4. Factors Influencing Local Pier Scour Depth
      5. 5. Abutment Scour Mechanics
   10. Lecture 10: Clear-Water vs. Live-Bed Scour, Complex Pier Scour Design, and Protection Measures
       1. 1. Clear-Water Scour vs. Live-Bed Scour
       2. 2. Complex Bridge Piers: Mechanics and Design Configurations
       3. 3. Effective Pier Diameter Formulation ($b_e$)
       4. 4. Empirical Scour Estimation Formulations
       5. 5. Bridge Pier Scour Countermeasures
       6. 6. Case Studies and Field Evaluations
   11. Lecture 11: River Training, Bank Protection, Guide Banks, and Nature-Based Solutions
       1. 1. Principles and Objectives of River Training
       2. 2. Spurs and Groynes (Transverse Training Works)
       3. 3. Guide Banks (Bell’s Bunds)
       4. 4. Permeable River Training Devices
       5. 5. Revetments and Longitudinal Bank Protection
       6. 6. Nature-Based Solutions (NBS) and Bio-Engineering
   12. Lecture 12: River Navigation, Channel Regulation, Dredging Engineering, and Environmental Impacts
       1. 1. Inland Waterway Transport (IWT) Fundamentals and Economics
       2. 2. Governing Factors of Channel Navigability
       3. 3. Engineering Methods for Navigability Improvement
       4. 4. Dredging Technology and Operational Phases
       5. 5. Navigation Infrastructure and Dredging in Bangladesh
       6. 6. Morphological and Environmental Impacts of Dredging

---

EWCE 461 / Section A: River Engineering — Comprehensive Course Outline

Thematic Curriculum Mapping

The content of EWCE 461 / Section A (River Engineering) spans twelve structured lectures, organized into five cohesive thematic modules:

EWCE 461 / Section A: River Engineering
│
├── Module I: Fluvial Geomorphology, Planform Patterns & Hydraulic Geometry
│   ├── Lecture 1: Rivers, Longitudinal Zones, Drainage Systems, and Channel Classification
│   ├── Lecture 2: Channel Patterns, Governing Factors, Hydraulic Geometry, and River Longitudinal Profiles
│   └── Lecture 3: Derived Hydraulic Relations, Meander Mechanics, and Bend Hydrodynamics
│
├── Module II: Channel Stability, Sediment Classification & Incipient Motion
│   ├── Lecture 4: Channel Stability Classification, Braiding Mechanics, and Sediment Load Modes
│   └── Lecture 5: River Hydrodynamics, Threshold of Movement, and the Shields Parameter
│
├── Module III: Flow Regimes, Bedforms, Settling Dynamics & Transport Estimation
│   ├── Lecture 6: Flow Regimes, Bed Configurations, Flow Resistance, and Particle Settling Mechanics
│   └── Lecture 7: Quantitative Estimation of Sediment Transport Rates (Bed, Suspended, and Total Load)
│
├── Module IV: Dynamic Equilibrium, Morphological Response & River Scour
│   ├── Lecture 8: Dynamic Adjustment of Rivers, Lane’s Relationship, Aggradation, and Degradation
│   ├── Lecture 9: River Scour Fundamentals, Bridge Hydraulics, and Local Scour Mechanisms
│   └── Lecture 10: Clear-Water vs. Live-Bed Scour, Complex Pier Scour Design, and Protection Measures
│
└── Module V: River Training, Bank Protection, Navigation & Dredging Engineering
    ├── Lecture 11: River Training, Bank Protection, Guide Banks, and Nature-Based Solutions
    └── Lecture 12: River Navigation, Channel Regulation, Dredging Engineering, and Environmental Impacts

---

Detailed Hierarchical Outline

---

Lecture 1: Rivers, Longitudinal Zones, Drainage Systems, and Channel Classification

1. Fundamentals of River Systems

• Definitions and Physical Hierarchy
  • Channelized flow vs. Overland / Sheet flow (depth, boundary confinement, flow concentration)
  • Surface runoff pathway: Rills / Rivulets $\to$ Gullies $\to$ Streams / Creeks $\to$ Rivers
• Energy Conversion Dynamics in Fluvial Systems
  • Conversion of potential energy (gravitational elevation head) into kinetic energy (velocity/motion)
  • Energy expenditure: Predominant dissipation via internal fluid and boundary friction
  • Effective work: $2\mathrm{\%}\text{to}4\mathrm{\%}$ of total potential energy expended in erosion and sediment transport
• The Three River Subsystems (Anatomy of a River Basin)
  • Collecting system (headwater branches): Tributary network funneling water and sediment
  • Transporting system (main trunk): Major channel conveyer where sediment routing, erosion, deposition, and pattern evolution occur
  • Dispersing system (distributary roots): Deltaic network dispersing water and sediment into sea, ocean, or lake

2. Longitudinal Zonation of Rivers (Schumm's Basin Model)

• Zone 1: Erosional Zone (Sediment Source)
  • Headwater reaches with steep channel gradients and high stream energy
  • Predominance of net bed degradation and downcutting
  • High frequency of confluences; bed material composed of boulders, cobbles, and coarse gravels
• Zone 2: Transfer / Transport Zone (Sediment Conveyance)
  • Intermediate reaches characterized by dynamic equilibrium between water/sediment inflow and outflow
  • Relatively stable average bed elevation; typical single-thread channels
  • Bed material dominated by gravels and sands of varying gradations
• Zone 3: Depositional Zone (Sediment Sink)
  • Lower reaches and deltaic plains characterized by net sedimentation and bed aggradation
  • Channel bifurcation, distributaries, and mild/flat hydraulic slopes
  • Bed material composed of fine sand, silt, and cohesive clay

3. Fluvial Morphology Terminology and Riparian Features

• Key Fluvial Definitions
  • Tributary vs. Distributary; Confluence vs. Outfall (River Mouth)
  • Upstream (headward) vs. Downstream (mouthward)
  • Thalweg: The continuous line linking the points of deepest water along the channel bed
• Riparian Topographic Features
  • Floodplain: Low-lying alluvial land adjacent to channels inundated at return periods of $1.5\text{to}2\text{years}$
  • Natural Levees, River Terraces (abandoned former floodplains), and Alluvial Ridges
• Comparative Analysis: River Outfalls (Delta vs. Estuary)
  • Dominant Forcing: Upland river flow (Delta) vs. Tidal energy and currents (Estuary)
  • Primary Geomorphic Process: Net deposition and progradation (Delta) vs. Net tidal scour/erosion (Estuary)
  • Physical Geometry: Triangular / fan-shaped multi-channel network (Delta) vs. Funnel-shaped wide opening (Estuary)
  • Land Building: Active seaward accretion of new land (Delta) vs. No seaward land building (Estuary)
  • Water and Soil Properties: Freshwater-dominated with high-fertility alluvium (Delta) vs. Brackish water with saline/marshy soils (Estuary)
  • Geographic Case Studies in Bangladesh: Lower Meghna Delta vs. Southwestern Coastal Estuaries

4. Drainage Network Patterns and Stream Topology

• Drainage Patterns (Collection System Typology)
  • Dendritic: Irregular branching tree-like pattern; develops on massive uniform rocks and flat-lying strata
  • Parallel: Sub-parallel channels formed on uniformly dipping or steep regional slopes
  • Radial: Drainage radiating outward from a central structural high (e.g., volcanic domes)
  • Trellis and Rectangular patterns: Geologically controlled by jointing and folding
• Stream Ordering (Horton-Strahler System)
  • First-order streams: Smallest fingertip unbranched tributaries
  • Higher-order confluence rules: Two streams of order $n$ join to form an order $n+1$ stream; a lower-order stream joining a higher-order stream leaves the order unchanged
  • Drainage basin hierarchy and Horton's morphometric laws

5. Multi-Criteria River Classification Schemes

• Classification by Nature of Substrate
  • Bedrock (Non-Alluvial) Rivers: Flow directly on solid rock; boundaries fixed and stable over geological time
  • Alluvial Rivers: Self-formed channels flowing through self-deposited sediment; mobile boundaries; continuous adjustments in cross-section, slope, and planform; morphologically unstable
• Classification by Groundwater Table Interaction
  • Influent (Losing) Streams: Water surface stands above local groundwater table; water seeps into banks/aquifer
  • Effluent (Gaining) Streams: Water surface lies below local groundwater table; receives baseflow discharge from aquifer
  • Seasonal and longitudinal reversals in perennial rivers
• Classification by Flow Permanency (Hydroperiod)
  • Perennial streams: Continuous flow year-round sustained by baseflow
  • Intermittent streams: Seasonal flow during wet periods, dry during dry seasons
  • Ephemeral streams: Flow exclusively during and immediately following storm events; dry beds for majority of the year
• Classification by Dominant Sediment Load Mode
  • Suspended-Load Channels: Carries fine sand, silt, and clay; low gradients; narrow and deep ($W/D<10$); high sinuosity; stable cohesive banks with minimal lateral shift
  • Mixed-Load Channels: Sinuous thalweg with alternate migrating bars; moderate width-to-depth ratio; shifting banklines and periodic cutoffs
  • Bed-Load Channels: High coarse sand/gravel transport; steep gradients; wide and shallow ($W/D>40$); multi-thread braided patterns; non-cohesive unstable banks; frequent bar migration and avulsion
• Classification by Planform / Morphology
  • Straight, Meandering, Braided, and Anastomosing rivers

6. Quantitative Planform Metrics

• Sinuosity Index ($SI$)
  • Mathematical definition: $$SI=\frac{\text{ChannelLength}(A\to B)}{\text{Valley}StraightDistance}\left(AB\right)\text{[PARSE ERROR: Undefined("RBrace")]}$$
  • Quantitative thresholds:
    • Straight channels: $SI\text{[PARSE ERROR: Undefined("Command("le")")]}1.05$ (rare in natural systems)
    • Sinuous channels: $1.05<SI\text{[PARSE ERROR: Undefined("Command("le")")]}1.50$
    • Meandering channels: $SI>1.50$ (can approach values of $4.0$)
• Planform Index ($PFI$) for Braided Rivers (Sharma, 2004)
  • Formulation: $$PFI=\left(\frac{\sum _{i=1}^N{T}_i}{B}\right)\times \frac{100}{N}$$ where $T_i$ = top width of individual channel $i$, $B$ = total braid belt flow width, $N$ = number of active channels
  • Classification thresholds:
    • Highly Braided: $PFI<3$
    • Moderately Braided: $3\text{[PARSE ERROR: Undefined("Command("le")")]}PFI\text{[PARSE ERROR: Undefined("Command("le")")]}20$
    • Low Braided: $PFI>20$
  • Case Study: Temporal analysis of the Jamuna River at RMJ11 (2005–2012 planform dynamics)

---

Lecture 2: Channel Patterns, Governing Factors, Hydraulic Geometry, and River Longitudinal Profiles

1. Planform Patterns and Fluvial Dynamics

• Primary Alluvial River Patterns
  • Straight rivers: High boundary resistance or young channels; artificial straightening creates dynamic instability
  • Meandering rivers: Single-thread sinuous curves; alternating erosion on outer concave banks and point bar accretion on inner convex banks
  • Braided rivers: Multi-thread channels separated by dynamic mid-channel bars and islands; unstable shifting thalwegs
• Empirical Planform Tracking
  • Long-term sinuosity tracking: Ganges River in Bangladesh ($1998\text{to}2025$)

2. Governing Factors of Alluvial Channel Morphology

• Discharge ($Q$)
  • Controls absolute channel dimensions: Width ($w$), Depth ($d$), meander wavelength ($M_L$), and amplitude
  • Independent variation of $Q$ does not dictate the fundamental planform type
• Valley Slope ($S$)
  • At a given discharge, meanders develop at lower slopes; braided channels develop at steeper slopes
  • Critical slope threshold separates single-thread meandering from multi-thread braided regimes
• Sediment Load and Grain Size
  • Abundant coarse bed load (gravel/sand) + steep slopes foster wide, shallow, braided patterns
  • Fine-grained wash/suspended load (silt/clay) fosters narrow, deep, stable meandering patterns
• Width-to-Depth Ratio ($W/D$)
  • Controlled by silt-clay percentage in perimeter sediment ($M$); high cohesion yields low $W/D$ ($<10$)

3. Hydrological Flood Response Classification

• Normal Rivers vs. Flashy Rivers
  • Normal rivers: Gradual rising and falling hydrograph limbs; broad flood attenuation capacity
  • Flashy rivers: Rapid hydrograph rise, extremely sharp flood peaks, and rapid recession
• Regional Application in Bangladesh
  • Flashy rivers in the Northeastern (Haor basin) and Southeastern hilly regions (e.g., Halda River at Narayanhat vs. Meghna River at Narsingdi)
  • Pre-monsoon early flash floods and catastrophic risks to Boro rice cultivation

4. Theory of Hydraulic Geometry (Leopold & Maddock)

• Core Parameters
  • River Stage ($h$) vs. Hydraulic Flow Depth ($d$); Discharge ($Q=V·A=w·d·u$)
• Power-Law Empirical Formulations
  • Width relation: $w=a{Q}^b$
  • Depth relation: $d=c{Q}^f$
  • Mean velocity relation: $u=k{Q}^m$
  • Mathematical continuity constraints: $$a·c·k=1.0\text{and}b+f+m=1.0$$
• Two Classical Formats of Hydraulic Geometry
  • At-a-Station Hydraulic Geometry: Describes temporal variations of $w,d,u$ as discharge fluctuates at a single fixed monitoring station over time
  • Downstream Hydraulic Geometry: Describes spatial scaling of $w,d,u$, and slope $S$ as discharge systematically increases downstream through tributary accumulation
• Practical Numerical Applications
  • Log-log regression fitting for width-discharge exponents ($b$) and calculation of corresponding velocity relationships ($m=1-b-f$)

5. Channel-Forming and Dominant Discharges

• Bankfull Discharge ($Q_{bf}$)
  • Definition: The discharge that completely fills the alluvial channel up to the active floodplain level without overtopping
  • Recurrence interval: Typically $1.5\text{to}2.0\text{years}$ on the annual maximum flood series
  • Quantitative identification criteria:
    • Depth corresponding to the absolute minimum width-to-depth ratio ($W/D$)
    • Inflection point where top width or cross-sectional area changes abruptly on stage-area curves
  • Numerical example: Determination of bankfull stage from river cross-section width-depth survey tables
• Dominant / Effective Discharge
  • The theoretical constant discharge that, if maintained indefinitely, would produce the same channel morphology and sediment transport capacity as the full hydrograph
  • Frequently equated to bankfull discharge or the $50\mathrm{\%}$ duration flow ($2year\text{[PARSE ERROR: Undefined("RBrace")]}$ return period)

6. River Longitudinal Profiles and Downstream Fining

• Concave-Upward Profile of Alluvial Rivers
  • Exponential longitudinal slope decline model: $$S=S_0{e}^{-\alpha x}$$
  • Integrated elevation profile: $$Y\left(x\right)=\frac{S_0}{\alpha }(e^{-\alpha x}-1)$$
• Physical Mechanisms Explaining Concave Profiles
  1. Particle size reduction via abrasion, grinding, and progressive hydraulic sorting requires a flatter slope downstream to maintain transport
  2. The upper watershed acts as the primary sediment factory; downstream reaches experience reduced sediment concentration due to floodplain deposition
  3. Tributary water inflows increase river discharge downstream faster than sediment load increases, reducing required equilibrium slope
• Sternberg’s Law of Downstream Sediment Fining
  • Exponential decrease in sediment particle weight: $$W=W_0{e}^{-\beta x}$$
  • Exponential decrease in median particle diameter: $$d=d_0{e}^{-cx}$$

---

Lecture 3: Derived Downstream Hydraulic Relations, Meander Theories, and Bend Hydrodynamics

1. Derived Downstream Hydraulic Geometry Relationships

• Downstream Resistance Variations
  • Manning’s roughness factor ($n$): Decreases downstream as bankfull discharge increases ($n\propto Q^{-0.16}$ to $Q^{-0.06}$)
  • Darcy-Weisbach friction factor ($f$): Decreases downstream with increasing discharge
  • Chezy’s resistance coefficient ($C$): Increases downstream ($C\propto Q^{0.39}$)
• Downstream Grain Size and Velocity Dynamics
  • Manning-Strickler linkage ($n\propto D_{{50}^{1/6}}$): Confirms median particle diameter $D_{50}$ decreases downstream with increasing $Q$
  • Downstream Velocity Trend: Mean flow velocity ($V$) exhibits a mild net increase or remains near-constant downstream ($V\propto d^{2/3}$ offsets slope flattening)

2. Planform Stability Boundaries and Thresholds

• Slope-Discharge Critical Boundaries (Lane / Leopold & Wolman)
  • Threshold equation separating meandering from braided channels: $$S_{cr}=0.012Q^{-0.44}\text{or}{S}_{cr}=0.06Q_{b^{-0.44}}$$
  • Meandering streams plot below $S_{cr}$; braided streams plot above $S_{cr}$
  • Critical slope declines systematically as discharge increases
• Boundary Strength and Multi-Channel Transitions
  • Relationship between normalized stream power / Froude parameter $(S/Fr)$ and inverse form ratio $(d/w)$
  • Bank strength index: High $d/w$ represents cohesive, resistant banks; low $d/w$ represents weak, non-cohesive banks
  • Multi-braid transition: As bank strength decreases, the number of active braids in a cross-section increases from single-thread to $2\mathrm{–}10$ braids

3. Hypotheses and Theories of Meander Formation

• 1. Earth's Rotation / Coriolis Force Theory
  • Deflection force normal to flow direction (deflects right in Northern Hemisphere, left in Southern Hemisphere); creates lateral flow imbalance
• 2. Disturbance Theory
  • Local obstructions, sediment shoals, or flood surges create asymmetric flow; induces cross-channel currents and lateral bank cutting
• 3. Excess Energy Dissipation Theory
  • Channels with excess energy dissipate stream power by increasing channel length (higher sinuosity), thereby reducing gradient toward minimum energy dissipation equilibrium
• 4. Secondary Flow / Transverse Circulation Theory
  • Turbulent boundary layer shear generates transverse secondary currents, causing alternating lateral scour and shoaling
• 5. General Wave Motion / Flow Instability Argument ("Flag Fluttering")
  • Dynamic instability in turbulent shear flows develops natural spatial oscillatory waveforms analogous to fluttering flags

4. Morphometric Scaling of Meander Geometry

• Universal Empirical Geometric Scaling (Leopold & Wolman Global Datasets)
  • Wavelength ($\mathrm{\Lambda }$ or $M_L$), Radius of curvature ($r_c$ or $R_M$), Channel width ($w$ or $B$), and Discharge ($Q$): $$\mathrm{\Lambda }=6.8w^{1.12}$$ $$\frac{\mathrm{\Lambda }}{4.53}=r_c⟹r_c\approx 1.5\text{to}2.5w$$ $$M_L=36Q^{0.5}$$ $$M_L=10.9B^{1.01}$$ $$M_L=4.7R_{M^{0.98}}$$ $$S=0.06Q_{b^{-0.44}}$$
• Pool-to-Pool Spacing Relationships
  • Riffle-to-riffle or pool-to-pool spacing typically equals $5\text{to}7$ channel widths ($SP\approx 5\mathrm{–}7w$)
  • Practical log-log regression problems: Deriving empirical meander length from spacing data

5. Bend Hydrodynamics and Helical Secondary Currents

• Superelevation Across Meander Bends
  • Dynamic balance between centrifugal acceleration ($F_c=\frac{m{V}^2}{r_c}$) and cross-stream hydrostatic pressure gradient ($\mathrm{\Delta }h$): $$\mathrm{\Delta }h=\frac{V^{2}w}{g{r}_c}$$
• Three-Dimensional Helical Flow Mechanism
  • Near the water surface: High longitudinal velocity causes centrifugal force to exceed the transverse pressure gradient $⟹$ radial flow directed toward the outer concave cut-bank
  • Near the riverbed: Boundary friction retards velocity; the inward hydrostatic pressure gradient exceeds centrifugal force $⟹$ radial flow directed toward the inner convex bank
  • Helical Spiral: Continuous three-dimensional spiral sweeping erodible bed material from outer cut-banks across the pool to deposit on inner point bars
• Cross-Sectional Asymmetry
  • Asymmetric deep scour pools at bend apexes vs. shallow, symmetric riffles/crossings at inflection points between bends

---

Lecture 4: Channel Stability Classification, Braiding Mechanics, and Sediment Load Modes

1. Mechanics and Genesis of River Braiding

• Primary Physical Causes of Braided Channels
  • Steep valley slopes providing high stream power
  • High sediment supply with coarse-textured non-cohesive bed and bank materials
  • Overloading: Sediment supply delivered to the reach exceeds the local transport capacity ($Q_s>\text{Capacity}$)
  • Deposition of coarse sediment into longitudinal and mid-channel bars; flow bifurcation into multi-thread channels
• Instability of Braided Systems
  • Highly dynamic, unpredictable thalweg shifts compared to the predictable migration cycles of meandering streams

2. Channel Stability Classification (Schumm's Fluvial Continuum)

• Classification Based on Sediment Load Type and Boundary Material
  • Pattern 1 & 3a (Suspended-Load Channels):
    • Bed material load $<3\mathrm{\%}$; wash load constitutes bulk of sediment
    • Narrow, deep cross-sections ($W/D<10$); high sinuosity ($SI>2.0$)
    • Highly stable cohesive silt-clay banks resistant to lateral shift
  • Pattern 3b (Mixed-Load Channels):
    • Bed material load $3\mathrm{\%}\text{to}11\mathrm{\%}$; intermediate width-to-depth ratio ($W/D\approx 10\mathrm{–}30$)
    • Sinuosity between $1.3\text{and}2.0$; active meander migration and alternate bar development
  • Patterns 4 & 5 (Bed-Load and Braided Channels):
    • Bed material load $>11\mathrm{\%}$; wide, shallow cross-sections ($W/D>40$)
    • Low sinuosity ($SI<1.3$); steep slopes; extensive submerged bars and vegetated islands
    • Non-cohesive sandy-gravel banks; rapid bank collapse and extreme planform instability

3. Planform Stability (PS) Force Ratios

• Planform Stability Ratio 1 ($P{S}_1$): Channel Boundary Resistance $$P{S}_1=\frac{\text{StreamErodingPower}}{\text{BoundaryMaterialStrength}}$$
  • Low $P{S}_1$: Boundary strength resists stream power $⟹$ stable single-thread meandering
  • High $P{S}_1$: Stream power easily overwhelms weak bank resistance $⟹$ channel widening and braiding
• Planform Stability Ratio 2 ($P{S}_2$): Sediment Conveyance Balance $$P{S}_2=\frac{\text{BedLoadSupply}}{\text{BedLoadTransportCapacity}}$$
  • Low $P{S}_2$: Stream has surplus capacity to evacuate incoming bed load $⟹$ single-thread equilibrium
  • High $P{S}_2$: Stream is overwhelmed by sediment delivery $⟹$ bed aggradation, mid-channel bar formation, and braiding

4. Comprehensive Sediment Load Classification Systems

• Classification Criteria
  • By physical nature and chemical state
  • By presence or absence in the channel bed material
  • By physical mode of transport through the water column
• Dissolved Load
  • Chemical solutes and ionic species ($HC{O}_{3^{-}},C{a}^{2+},S{O}_{4^{2-}},C{l}^{-}$, etc.)
  • Travels at the mean velocity of water; negligible direct mechanical impact on channel morphology
• Bed Load (Contact and Saltation Load)
  • Coarse sediment grains moving in continuous or intermittent contact with the bed via rolling, sliding, or saltating
  • Confined to a thin bed layer (approximately 2 grain diameters thick)
  • Purely bed-material load; strictly hydraulic capacity-limited
• Suspended Load
  • Particulate matter (fine sand, silt, clay) maintained in the main flow column by upward turbulent eddies
  • Upward turbulent velocity fluctuations must equal or exceed particle fall velocity ($\omega$)
  • Composed of suspended bed-material load plus wash load
• The Bed Load vs. Suspended Load Continuum
  • No sharp physical demarcation; dynamic, continuous exchange of grains between bed load and suspended load
  • A grain alternates between bed load and suspension depending on instantaneous turbulent velocity fluctuations
• Wash Load Characteristics
  • Ultra-fine particles (clay and fine silt, typically $<0.002\mathrm{–}0.0625\text{mm}$)
  • Maintained in permanent suspension by thermal molecular agitation (Brownian motion)
  • Supply-limited from watershed soil erosion; independent of local open-channel transport capacity
  • Uniform vertical concentration profile from riverbed to water surface
• Vertical Distributions of Concentration and Grain Size
  • Coarse sand exhibits pronounced vertical gradients (highest concentration and coarsest sizes near the bed)
  • Fine silt and clay exhibit nearly uniform vertical concentration profiles

---

Lecture 5: River Hydrodynamics, Threshold of Movement, and the Shields Parameter

1. Hydrodynamic Forces and Boundary Shear Stress

• Reach Force Balance in Open Channels
  • Downstream impelling gravity component: $F_g=\rho gAL\sin \alpha$
  • Opposing boundary friction force: $F_R={\tau }_{0}PL\cos \alpha$
  • Steady uniform flow equilibrium ($F_g=F_R$): $${\tau }_0=\rho gR\tan \alpha =\rho gRS$$
  • Wide-channel simplification ($R\approx d\text{or}h$): $${\tau }_0=\rho ghS$$
• Resistance Equations and Quadratic Shear Dependence
  • Chezy formula: $V=C\sqrt{hS}=C\sqrt{\frac{{\tau }_0}{\rho g}}$
  • Demonstrates that boundary shear stress is proportional to the square of flow velocity: $${\tau }_0\propto V^2$$
• Stream Power Formulations
  • Total stream power: $\mathrm{\Omega }=\gamma QS$ (rate of energy dissipation per unit channel length)
  • Unit stream power: $\omega =\frac{\mathrm{\Omega }}{w}={\tau }_{0}V\propto V^3$

2. Sediment Transport Concepts: Competence vs. Capacity

• Competence (Calibre)
  • The maximum particle size (diameter/weight) that a flow is capable of entraining and transporting
  • A strictly hydraulic limitation governed by critical shear stress and local velocity
• Capacity
  • The maximum total mass or volume rate of sediment of all sizes that a stream can carry for a given discharge
  • Governed by total stream power, channel geometry, and energy availability

3. Mechanics of Incipient Motion (Threshold of Movement)

• Forces Acting on an Individual Submerged Bed Particle
  • Submerged gravitational weight: $$F_g=\frac{4}{3}\pi {\left(\frac{d}{2}\right)}^3({\rho }_s-\rho )g=c_1{d}^3({\rho }_s-\rho )g$$
  • Hydrodynamic drag force ($F_D$) acting parallel to the bed due to surface shear and pressure differential
  • Hydrodynamic lift force ($F_L$) acting normal to the bed generated by vertical velocity gradients (Bernoulli effect)
• Equilibrium Limit for Initial Motion
  • Grain displacement occurs when disturbing forces overcome resisting gravity forces: $$F_D+F_L>F_g$$
• Theoretical Derivation of Critical Shear Stress (${\tau }_c$)
  • Moment balance at the point of rolling/sliding: $$a_1{c}_1{d}^{3}g({\rho }_s-\rho )\sin \alpha =a_2{c}_2{d}^2{\tau }_c\cos \alpha$$ $${\tau }_c=\left(\frac{a_1{c}_1}{a_2{c}_2}\right)dg({\rho }_s-\rho )\tan \alpha$$
  • Fundamental Hydrodynamic Principle: For coarse non-cohesive sediment, the critical shear stress required for entrainment depends directly and solely on particle diameter $d$ and submerged specific weight $({\rho }_s-\rho )$

4. The Shields Parameter and Incipient Motion Criteria

• Dimensionless Critical Shear Stress (Shields Parameter, ${\theta }_c$)
  • Definition: $${\theta }_c=\frac{{\tau }_c}{({\rho }_s-\rho )gd}$$
  • Plotted as a function of Shear / Boundary Reynolds Number: $$R{e}_{*}=\frac{u_{*}d}{\nu }\text{where}{u}_{*}=\sqrt{\frac{{\tau }_0}{\rho }}$$
• Key Regimes of the Shields Diagram
  1. Fine Cohesive Soil (Left side): High critical stress due to electrochemical interparticle cohesion and particles being sheltered within the laminar sublayer
  2. Sand Range Minimum ($0.06\mathrm{–}2.0\text{mm}$): Absolute lowest critical shear stress (${\theta }_c\approx 0.03\mathrm{–}0.035$); particles have minimal mass and lack cohesion, making sand the most easily erodible sediment
  3. Coarse Turbulent Regime (Right side): Grain size exceeds laminar sublayer thickness; ${\theta }_c$ levels off to a near-constant value (${\theta }_c\approx 0.045\mathrm{–}0.06$)
• Hjulström-Type Dynamics
  • Comparison of critical velocities required for erosion, transport, and deposition across sediment grain sizes
  • Cohesive clays require higher entrainment velocities than sand; once in suspension, wash load never settles under natural river velocities

5. Engineering Applications and Numerical Calculations

• Assessment of bed mobility: Calculating boundary shear stress ${\tau }_0=\rho ghS$ and comparing with critical shear stress ${\tau }_c={\theta }_c({\rho }_s-\rho )gd$
• Calculating threshold discharge and water depth to initiate motion in gravel-bed canals

---

Lecture 6: Flow Regimes, Bed Configurations, Flow Resistance, and Particle Settling Mechanics

1. Flow Regimes in Alluvial Sand-Bed Channels (Simons & Richardson)

• Lower Flow Regime ($Fr<1$, Bedforms and Water Surface Out of Phase)
  1. Plane bed without sediment motion: Flat bed prior to reaching critical shear stress
  2. Ripples ($Fr\ll 1$): Small, triangular bedforms formed only in sediment with $d_{50}<0.7\text{mm}$; wavelength $<0.4\text{m}$, height $<0.04\text{m}$; independent of flow depth
  3. Dunes with superposed ripples: Intermediate transitional stage with increasing flow intensity
  4. Dunes ($Fr<1$): Large, asymmetric triangular sand waves with gentle upstream luff slopes and steep downstream lee slopes (angle of repose $\approx {30}^{\circ }\mathrm{–}{35}^{\circ }$); flow separation and eddies in the lee trough; wavelength scales with flow depth
  5. Washed-out dunes: Progressive flattening of dune crests as flow intensity rises
• Transitional Flow Regime ($Fr\approx 1$)
  • Destruction of dunes; bed configuration rapidly transitions between dunes and plane bed; highly unstable flow resistance and sediment transport rates
• Upper Flow Regime ($Fr>1$, Bedforms and Water Surface In Phase)
  1. Plane bed with active sediment motion: Smooth bed with intense sheet transport; low Manning’s $n$ ($0.010\mathrm{–}0.013$)
  2. Standing waves ($Fr>1$): Symmetrical sinusoidal bed waves and water surface waves in perfect phase; no wave breaking
  3. Antidunes ($Fr>1$): Symmetrical bedforms that build up, move upstream or downstream, and break violently upstream against the current
  4. Chutes and pools ($Fr\gg 1$): High-energy supercritical flow characterized by alternating chutes of rapid flow and pools of turbulent deceleration
• Systematic Parameter Variations Across Bedform Sequences
  • Variation of Manning's $n$, Darcy-Weisbach $f$, stream power, and sediment transport rate
  • Manning's $n$ reaches a maximum during the dune stage (high form drag) and drops sharply during the upper plane bed stage

2. Alluvial Flow Resistance and Roughness Partitioning

• Roughness Decomposition Concept
  • Total boundary shear stress partitioned into grain and form components: $${\tau }_0={\tau }^{\prime }+{\tau }^{\prime }\text{'}\text{and}{k}_s=k_{s,\text{grain}}+k_{s,\text{form}}$$
• Grain (Skin) Roughness
  • Frictional resistance generated directly by individual sediment grains
  • Standard formulation: $k_s\approx 3d_{90}$
  • Dominates flow resistance in gravel-bed and boulder-bed streams
• Form Roughness
  • Form resistance generated by dynamic pressure differentials across bedforms
  • Acceleration on the upstream luff face (low pressure) vs. separation eddies on the downstream lee face (high pressure)
  • Dominates total alluvial resistance in sand-bed rivers
• Regime Flow Resistance Formulations
  • Lacey's silt factor: $f=1.76\sqrt{d_{50}}$
  • Lacey-Manning relationship: $n=0.022f^{0.2}$
  • Van Rijn equivalent roughness relationships ($k_{s,1},k_{s,2}$)

3. Particle Mechanics: Settling, Terminal, and Fall Velocity ($\omega$)

• Particle Geometric Descriptors
  • Sphericity: $$\text{Sphericity}=\frac{A_s}{A_p}$$ (Surface area of equivalent volume sphere divided by actual particle surface area)
  • Corely Shape Factor ($S.F.$): $$S.F.=\frac{c}{\sqrt{ab}}$$ (where $a$ = major axis, $b$ = intermediate axis, $c$ = minor axis)
• Hydrodynamic Force Balance for Settling Particles
  • Submerged particle weight: $$F_g=\frac{\pi }{6}{d}^3({\rho }_s-\rho )g$$
  • Opposing drag and viscous resisting forces ($F_v,F_d$)
• Particle Reynolds Number ($N_{Re}$) $$N_{Re}=\frac{\omega d}{\nu }$$
  • Governs whether settling behavior is dominated by viscous shear or inertial impact
• Stokes' Settling Law for Fine Particles ($N_{Re}<1$, Laminar Flow)
  • Balance of submerged weight and viscous resistance ($F_v=3\pi d\mu \omega$): $$\omega =\frac{1}{18}\frac{({\rho }_s-\rho )g{d}^2}{\mu }$$
  • Key Proportionality: Settling velocity $\omega$ is proportional to the square of diameter ($\omega \propto d^2$)
  • Inversely proportional to dynamic fluid viscosity $\mu$ (colder water reduces fall velocity)
  • Theoretical Limitations: Assumes smooth spherical particles; neglects flocculation; ignores boundary turbulence and inertial impact
• Inertial Settling Law for Coarse Particles ($N_{Re}>1000$, Turbulent Flow)
  • Balance of submerged weight and inertial impact drag: $$F_d=C_D\left(\frac{\pi }{4}{d}^2\right)\rho \frac{{\omega }^2}{2}$$ $$\omega =\sqrt{\frac{4g({\rho }_s-\rho )d}{3C_D\rho }}=\sqrt{\frac{4g(S_s-1)d}{3C_D}}$$
  • Key Proportionality: Settling velocity $\omega$ is proportional to the square root of diameter ($\omega \propto \sqrt{d}$)
  • Constant drag coefficient for spheres in turbulent flow: $C_D\approx 0.44$
• Transitional Settling Formulation ($1<N_{Re}<1000$)
  • Empirical drag coefficient formulation: $$C_D=\frac{24}{N_{Re}}+\frac{3}{\sqrt{N_{Re}}}+0.34$$
• Settling Duration Comparison (Through a $10m\text{[PARSE ERROR: Undefined("RBrace")]}$ Water Column)
  • Very coarse sand ($1.0\text{mm}$, $\omega =64\text{cm}s\text{[PARSE ERROR: Undefined("RBrace")]}$): $16\text{seconds}$
  • Fine sand ($0.1\text{mm}$, $\omega =0.64\text{cm}s\text{[PARSE ERROR: Undefined("RBrace")]}$): $26\text{minutes}$
  • Coarse silt ($0.01\text{mm}$, $\omega =0.0064\text{cm}s\text{[PARSE ERROR: Undefined("RBrace")]}$): $43\text{hours}$
  • Fine clay ($0.001\text{mm}$, $\omega =0.000064\text{cm}s\text{[PARSE ERROR: Undefined("RBrace")]}$): $180\text{days}$

---

Lecture 7: Quantitative Estimation of Sediment Transport Rates

1. Bed Shear Stress and Shear Velocity Foundations

• Boundary shear stress: ${\tau }_b={\rho }_{w}ghS$
• Bed shear velocity: $$u_{*}=\sqrt{\frac{{\tau }_b}{{\rho }_w}}=\sqrt{ghS}$$
• Dimensionless bed friction coefficient: $$C_f=\frac{{\tau }_b}{{\rho }_w{\omega }^2}={\left(\frac{u_{*}}{\omega }\right)}^2$$

2. Suspended Sediment Transport and the Rouse Profile

• Mathematical Definition of Suspended Transport
  • Depth-integrated unit suspended sediment flux: $$q_s={\int }_a^{h}c\left(y\right)·u\left(y\right)dy$$ where $a$ is the reference level near the bed, $C_a$ is reference concentration, and $h$ is total water depth
• The Rouse Concentration Profile Equation $$\frac{C\left(y\right)}{C_a}={\left(\frac{h-y}{y}·\frac{a}{h-a}\right)}^{R_n}$$ where the Rouse Number (Suspension Parameter) is defined as: $$R_n=\frac{\omega }{\kappa u_{*}}$$ (with von Kármán constant $\kappa =0.40$)
• Physical Interpretation of the Rouse Number ($R_n$)
  • $R_n>2.5$: Bed-load transport only; sediment grains cannot be held in suspension
  • $1.8<R_n\text{[PARSE ERROR: Undefined("Command("le")")]}2.5$: Incipient / partial suspension; particles attempt to enter suspension but drop back
  • $R_n<1.8$: Full, sustained sediment suspension threshold
  • $R_n<0.25$: Wash load; turbulence completely dominates settling, producing uniform vertical distribution
• Limitations of the Rouse Equation
  • Mathematical boundary condition forces $C\left(y\right)=0$ at the free surface ($y=h$), which deviates from real-world observations for fine silts
• Suspended Sediment Discharge Rating Models
  • Discharge rating curve: $Q_s=C_{0}Q$
  • Rubey’s theoretical concentration model: $C_{av}\propto R^{1/2}{S}^{2/3}$ combined with Chezy yields $Q_s\propto h^2{S}^{7/6}$
• Practical Numerical Applications
  • Calculating suspended sediment concentration at varying heights above the bed from a measured reference concentration $C_a$
  • Determining critical fall velocity threshold required for suspended particles to transition to bed load

3. Bed Load Transport Formulations

• Physical Characteristics of Bed Load
  • Confined to a thin layer $\sim 2d$ above bed; capacity-limited; typically represents $1\mathrm{\%}\text{to}10\mathrm{\%}$ of suspended load in large lowland rivers
• Dune Migration Method (Kinematic Method)
  • Computation of unit bed load transport based on measured dune geometry and migration velocity: $$q_b=\alpha v_b{\rho }_s(1-p)L\mathrm{\Delta }H$$ where $\alpha \approx 0.5\mathrm{–}0.7$ (shape factor), $v_b$ = dune migration speed, $p$ = bed porosity, $L$ = dune wavelength, $\mathrm{\Delta }H$ = dune height
  • Worked numerical calculation: Estimating total bed load discharge ($Q_b=q_b·W$)
• Meyer-Peter & Müller (M-P&M) Formula for Gravel-Bed Channels
  • Unit bed load transport rate per unit channel width: $$q_b=0.253({\tau }_0-{\tau }_{cr}{)}^{3/2}$$
  • Applied specifically where transport is dominated by coarse gravel and bed shear stress exceeds Shields critical shear stress

4. Total Bed Material Load Prediction Methods

• Practical Engineering Challenges
  • Yang (1996) demonstrated that total load estimates from different formulae may deviate by up to a factor of 4
• Recommended Total Load Formulations (Garde & Ranga Raju)
  • Ackers and White Formula (1973/1992):
    • Dimensionless grain size: $$D_{gr}=D_{35}{\left[\frac{g({\gamma }_r-1)}{{\nu }^2}\right]}^{1/3}$$
    • Sediment mobility parameter $F_{gr}$ and dimensionless transport function $G_{gr}=C{\left(\frac{F_{gr}}{A}-1\right)}^m$
    • Applies total boundary shear to fine suspended fractions and grain shear to coarse bed load fractions
  • Engelund and Hansen Formula (1967):
    • Stream power / shear approach: Simple and robust across sandy rivers, though neglects fluid temperature/viscosity
  • Other leading formulas: Brownlie (1981), Yang (1973), Karim-Kennedy (1990), Van Rijn (1984)
• Sediment Load Applications in Major Bangladesh Rivers
  • Methods applied to Ganges, Brahmaputra-Jamuna, Meghna, Teesta, and Surma-Kushiyara rivers

---

Lecture 8: Dynamic Adjustment of Rivers, Lane’s Relationship, Aggradation, and Degradation

1. Emory Lane’s Fluvial Equilibrium Relationship

• The Lane Balance Principle (1955)
  • Fundamental qualitative relationship governing alluvial river equilibrium: $$Q_s·d_{50}\propto Q_w·S$$ $$SedimentDischarge)\text{[PARSE ERROR: Undefined("RBrace")]}\times SedimentSize)\text{[PARSE ERROR: Undefined("RBrace")]}\propto WaterDischarge)\text{[PARSE ERROR: Undefined("RBrace")]}\times ChannelSlope)\text{[PARSE ERROR: Undefined("RBrace")]}$$
  • Proportional rearrangement for sediment discharge ($Q_s$): $$Q_s\propto \frac{Q_{w}S}{d_{50}}\propto \frac{\text{StreamPower}}{\text{Calibre}}$$
  • Proportional rearrangement for sediment concentration ($C_s=Q_s/Q_w$): $$C_s\propto \frac{S}{d_{50}}$$
• Equilibrium Adjustment Rules (Dynamic Balancing)
  • Increase in sediment load ($Q_s\uparrow$) $⟹$ Channel aggradation, Slope increases ($S\uparrow$)
  • Decrease in sediment load / clear water ($Q_s\downarrow$) $⟹$ Channel degradation, Slope decreases ($S\downarrow$)
  • Increase in water discharge ($Q_w\uparrow$) $⟹$ Channel degradation, Slope decreases ($S\downarrow$)
  • Decrease in water discharge / diversion ($Q_w\downarrow$) $⟹$ Channel aggradation, Slope increases ($S\uparrow$)
• Critical Limitations of Lane’s Qualitative Model
  • Non-quantitative: Cannot predict the exact numerical magnitude of morphological change
  • Lacks a time domain: Cannot predict the rate or timeframe required to achieve new equilibrium
  • Ignores planform and cross-sectional geometry adjustments (width and depth adjustments)

2. Quantitative Formulation of Lane’s Relation

• Analytical Derivation
  • Combining sediment transport function ($q_s={\rho }_w{k}_1{V}^3$), Darcy-Weisbach head loss ($S_0=\frac{f{V}^2}{gR}$), and Manning-Strickler equation ($n\propto d_{s^{1/6}}$): $$Q_s{\left(\frac{d_s}{y}\right)}^{1/3}\propto Q_w{S}_0$$
• Pre- and Post-Development Dimensionless Ratio Formulation
  • Dimensionless parameter ratios ($a=Q_{s2}/Q_{s1}$, $b=d_{s2}/d_{s1}$, $c=y_2/y_1$, $d=Q_{w2}/Q_{w1}$, $e=S_{02}/S_{01}$): $$a{\left(\frac{b}{c}\right)}^{1/3}=d·e$$
• Worked Engineering Example
  • Quantitative prediction of post-diversion bed slope adjustment and aggradation downstream of an irrigation intake

3. Channel Aggradation and Degradation Processes

• Definitions and Physical Dynamics
  • Aggradation: Long-term net elevation increase of the riverbed caused by sediment supply exceeding local transport capacity
  • Degradation: Long-term net lowering of the riverbed caused by sediment transport capacity exceeding sediment supply
• The Armouring Phenomenon
  • Selective winnowing of fine particles during channel degradation leaves an interlocking layer of coarse particles
  • Armour layer prevents further erosion once hydraulic shear stress can no longer mobilize the coarse armoring size
• Primary Driving Factors of Long-Term Bed Elevation Instability
  • Dam and reservoir impoundments (upstream aggradation vs. downstream degradation)
  • Watershed land-use changes (urbanization, deforestation)
  • Artificial or natural meander cutoffs
  • In-stream sand and gravel mining
  • Trans-boundary or upstream water diversions
  • Tectonic faulting and earthquakes (e.g., Deresuyu River faulting case study)
• Comparative Analysis: Aggradation vs. Degradation
  • Base level, kinetic energy state, sediment balance, and equilibrium timescales
• Hydraulic Structure Impacts (Dams and Barrages)
  • Upstream: Backwater formation reduces flow velocity $⟹$ deltaic aggradation
  • Downstream: Deposition in reservoir releases "clear water" (hungry water) with high turbulence $⟹$ severe downstream degradation and bank collapse

4. Mathematical Modeling & Field Case Studies

• One-Dimensional Sediment Continuity Equation (Exner Equation) $$\frac{\mathit{\partial }Z}{\mathit{\partial }t}+\frac{1}{{\gamma }_s(1-\lambda )}\frac{\mathit{\partial }{q}_t}{\mathit{\partial }x}=0$$ where $Z$ is bed elevation, $\lambda$ is bed porosity, and $q_t$ is total sediment transport rate
• Case Study: Aggradation and Degradation of Khowai River, Bangladesh (Hossain, 1997)
  • Observed narrowing of channel top width by $18.2\text{m}$ ($3\text{m}year\text{[PARSE ERROR: Undefined("RBrace")]}$) due to berm siltation between 1988 and 1994
  • Thalweg analysis: Upper reach experienced aggradation at $67\text{mm}year\text{[PARSE ERROR: Undefined("RBrace")]}$; lower reach experienced degradation at $83\text{mm}year\text{[PARSE ERROR: Undefined("RBrace")]}$
  • HEC-6 1D morphological modeling predictions of long-term cumulative bed aggradation

---

Lecture 9: River Scour Fundamentals, Bridge Hydraulics, and Local Scour Mechanisms

1. Overview and Additive Components of Bridge Scour

• Definition of Scour
  • The erosive removal of bed and bank material by flowing water around hydraulic structures
• The Three Additive Components of Total Scour
  1. General Scour (Long-Term Bed Elevation Changes): Natural aggradation or degradation across the entire river reach
  2. Contraction Scour: Channel-wide bed lowering caused by narrowing of the river flow area by bridge abutments, piers, and approach embankments
  3. Local Scour: Highly localized scour holes generated by vortex flow patterns around individual piers and abutments $$\text{TotalScourDepth}\left(d_{\text{total}}\right)=d_{\text{general}}+d_{\text{contraction}}+d_{\text{local}}$$

2. Bridge Hydraulics and Contraction Scour

• Hydraulic Response to Bridge Constriction
  • Reduction in waterway cross-sectional area produces upstream backwater rise and increased energy head
  • Flow acceleration through the bridge opening increases boundary shear stress above critical thresholds
• Quantitative Evaluation: With-Bridge vs. Without-Bridge Conditions
  • Numerical modeling comparison of flood levels, flow area, and mean approach velocities

3. Mechanics of Local Pier Scour

• Downflow and Vortex Dynamics at Bridge Piers
  • Stagnation pressure on the upstream pier face creates a strong downward vertical pressure gradient
  • Horseshoe Vortex System: Downward jet impinges on the bed and rolls up into a coherent horseshoe-shaped vortex wrapping around the pier base; excavates bed material from the front and sides of the pier
  • Wake Vortex System: Vertical-axis low-pressure vortices shed from the downstream separation points of the pier; suck up and transport entrained sediment downstream
• Scour Equilibrium Mechanism
  • As the scour hole deepens, the intensity of the horseshoe vortex diminishes until sediment transport out of the hole equals sediment transport into it

4. Factors Influencing Local Pier Scour Depth

• Flow Velocity ($V$)
  • Scour depth scales with velocity ratio: $d_s\propto (V/V_c{)}^n\times D$
• Flow Depth ($y$) and Pier Width ($b_p$) Regimes
  • Narrow Piers ($y/b_p>1.4$): Pier width controls scour depth ($d_s\propto b_p$); depth influence is negligible
  • Wide Piers ($y/b_p<0.2$): Flow depth controls scour depth ($d_s\propto y$)
  • Intermediate Piers ($0.2\text{[PARSE ERROR: Undefined("Command("le")")]}y/b_p\text{[PARSE ERROR: Undefined("Command("le")")]}1.4$): Scour depth depends on both pier width and flow depth
• Pier Nose Shape ($K_1$ Correction Factor)
  • Streamlining reduces horseshoe vortex strength
  • Square nose ($K_1=1.10$, $+20\mathrm{\%}$ higher scour than sharp nose), Round nose ($K_1=1.00$), Circular ($K_1=1.00$), Sharp nose ($K_1=0.90$)
• Angle of Attack of Flow ($\theta$ and $K_2$ Correction Factor)
  • Skewed flow exposes the entire pier length ($L$) to the current, drastically increasing effective obstruction width: $$K_2={\left(\cos \theta +\frac{L}{a}\sin \theta \right)}^{0.65}(\text{maximum}{K}_2=5.0)$$
  • For angles of attack $\theta >5^{\circ }$, the angle of attack dominates over pier nose shape ($K_1$ set to $1.0$)
• Bed Condition ($K_3$) and Bed Armoring ($K_4$) Correction Factors
  • $K_3$ accounts for bedforms (plane bed, small dunes, medium/large dunes)
  • $K_4$ accounts for coarsening and armoring by bed materials ($D_{50}\text{[PARSE ERROR: Undefined("Command("ge")")]}2\text{mm},D_{95}\text{[PARSE ERROR: Undefined("Command("ge")")]}20\text{mm}$)

5. Abutment Scour Mechanics

• Horizontal vortex formation along the embankment toe and downstream vertical wake vortex
• Comparative formulas for abutment scour (Laursen, Melville, Froehlich, Rahman & Haque)

---

Lecture 10: Clear-Water vs. Live-Bed Scour, Complex Pier Scour Design, and Protection Measures

1. Clear-Water Scour vs. Live-Bed Scour

• Fundamental Hydraulic Distinctions
  • Clear-Water Scour:
    • Occurs when approach flow velocity is less than or equal to critical velocity ($V\text{[PARSE ERROR: Undefined("Command("le")")]}{V}_c$)
    • Upstream bed is stationary; no bed material is transported into the scour hole
    • Maximum equilibrium scour depth is higher because there is no incoming sediment to replenish the hole
    • Scour hole develops slowly over an extended duration
  • Live-Bed Scour:
    • Occurs when approach flow velocity exceeds critical velocity ($V>V_c$)
    • Upstream bed is mobile; sediment is continuously transported into the scour hole
    • Equilibrium is reached rapidly; scour depth fluctuates cyclically as bedforms (dunes) pass through the bridge opening
    • Maximum scour depth is lower than clear-water scour due to continuous sediment replenishment and turbulence damping by suspended sediment
• Critical Approach Velocity Formulation (Laursen, 1963) $$V_c=10.95y_{1^{1/6}}{d}_{{50}^{1/3}}$$

2. Complex Bridge Piers: Mechanics and Design Configurations

• Structural Composition of Complex Piers
  • Multi-component assembly: Pier Column / Stem, Pile Cap / Footing, and Deep Pile Group
• Engineering Rationale for Complex Piers
  • Structural robustness: Resisting enormous vertical loads and transferring lateral forces (earthquake, wind, braking, ship collision)
  • Deep water foundations in alluvial floodplains: Bypassing weak, unconsolidated upper bed strata to anchor in deep bearing layers
• The Four Local Scour Cases Based on Pile Cap Elevation ($Y$)
  • Case I ($d_s<Y$): Top of pile cap remains buried deep below the base of the scour hole; behaves hydraulically as a simple narrow pier
  • Case II ($d_s>Y$, pile cap below bed level): Top of pile cap is exposed within the scour hole below bed level; intercepts the downward jet and reduces scour depth
  • Case III ($Y\text{[PARSE ERROR: Undefined("Command("le")")]}0$, pile cap at or above bed level): Pile cap sits at or above the original bed level; acts as a large obstruction, expanding the horseshoe vortex and dramatically increasing scour
  • Case IV: Pile cap extends partially or fully above the water surface level

3. Effective Pier Diameter Formulation ($b_e$)

• Equivalent Diameter Formula for Complex Piers $$b_e=b\left(\frac{d+Y^{*}}{d+b^{*}}\right)+b^{*}\left(\frac{b^{*}-Y^{*}}{b^{*}+d}\right)$$ where $b$ = pier column diameter, $b^{*}$ = pile cap width, $d$ = flow depth, and $Y^{*}$ = vertical distance from bed to top of pile cap
• Numerical Impact: Demonstrates that a complex pier with an exposed pile cap can increase scour depth by more than $2.2$ times compared to a simple pier

4. Empirical Scour Estimation Formulations

• HEC-18 CSU Equation (Richardson & Davis / FHWA) $$\frac{y_s}{y_1}=2.0K_1{K}_2{K}_3{K}_4{\left(\frac{a}{y_1}\right)}^{0.65}F{r}_{1^{0.43}}$$
• Indian Subcontinent Standard: Lacey-Inglis Method
  • Regime flow depth: $$D_{LQ}=0.47{\left(\frac{Q}{f}\right)}^{1/3}$$
  • Maximum scour depth below highest flood level (HFL): $$D_{se}=2.0D_{LQ}$$
• Comparative Formula Spectrum
  • Formulations by Breusers (1965), Laursen (1963), Neill (1987), Jain & Fischer (1980), Chitale (1988), Melville (1997)
• Numerical Worked Problems
  • Scour calculations for normal piers vs. complex piers
  • Comprehensive bridge scour design problem for the Kumar River Bridge

5. Bridge Pier Scour Countermeasures

• Armoring Countermeasures (Passive Resistance Layer)
  • Rock riprap, wire gabions, precast concrete blocks, and sand-filled geotextile bags
  • Provide a physical barrier that resists shear stresses and protects underlying fine sediments
• Flow-Altering Countermeasures (Active Vortex Mitigation)
  • Collars: Horizontal discs attached to the pier near bed level to deflect the downward jet and shield the bed from the horseshoe vortex
  • Sacrificial Piles: Upstream piles placed in front of the main pier to create wake zones, dissipate kinetic energy, and deflect approach velocity
  • Design Rules for Sacrificial Piles: Diameter $d=0.2\mathrm{–}0.5D$, Height $=0.6y$, Number $=3\mathrm{–}6$, Spacing $=2\mathrm{–}3d$, Distance upstream $=1\mathrm{–}3D$; provides $20\mathrm{\%}\text{to}60\mathrm{\%}$ scour reduction
• Application Constraints in Bangladesh
  • Sacrificial piles are not deployed in major Bangladesh rivers due to massive scour depths ($>60\text{m}$) and severe hazards to river navigation
  • Standard national practice: Deep pile foundations designed for maximum scour plus riprap, CC blocks, and geobags

6. Case Studies and Field Evaluations

• Field evaluation of local scour around complex piers of the Dhaleswari Bridge on the Kaliganga River (Islam, BUET IWFM Thesis, 2018)
• Nationwide analysis of scour design formulas across 239 bridges in Bangladesh (Mondal, 2022)

---

Lecture 11: River Training, Bank Protection, Guide Banks, and Nature-Based Solutions

1. Principles and Objectives of River Training

• Definition of River Training
  • Engineering structures and earthworks constructed in a river to guide and confine flow, regulate bed/bank configuration, and ensure safe movement of water and sediment
• Core Engineering Objectives
  • Safe and rapid conveyance of flood discharges without overtopping
  • Efficient routing of sediment load to prevent shoaling
  • Bank stabilization and prevention of lateral channel migration
  • Maintenance of adequate navigable depths and draft
  • Prevention of outflanking (bypassing) of bridges, barrages, and diversion structures
• Structural Classification
  • Transverse Structures (Active): Installed perpendicular to flow; divert flow and reduce boundary velocities (spurs, groynes, bandals, porcupines)
  • Longitudinal Structures (Passive): Installed parallel to river banks; absorb wave/current energy and armor boundaries (guide banks, revetments, riprap, dykes)
  • Hardcore Works (rigid concrete/rock structures) vs. Softcore Works (bio-engineering / Nature-Based Solutions)

2. Spurs and Groynes (Transverse Training Works)

• Functions of Spurs
  • Deflecting flow away from eroding banks
  • Creating slack water zones (slack flow) to induce silt deposition and establish a stable bankline
  • Channelizing wide braided rivers into defined navigation channels
• Classification by Permeability
  • Impermeable Spurs: Solid core of earth, sand, or gravel armored with rock pitching or concrete blocks; completely blocks and deflects flow
  • Permeable Spurs: Timber stakes, bamboo cages, or concrete frameworks; retards flow velocity, dissipates kinetic energy, and promotes rapid sedimentation
• Classification by Orientation / Alignment
  • Attracting Spurs: Point downstream; attract flow toward their head and maintain a deep channel along that bank
  • Deflecting Spurs: Set perpendicular (${90}^{\circ }$) to the bankline; divert flow without strong repulsion; create a slack zone behind
  • Repelling Spurs: Point upstream with an inclination of ${10}^{\circ }\text{to}{60}^{\circ }$ to the normal; repel current toward the opposite bank; promote heavy silt deposition upstream of the spur
• Classification by Head Geometry
  • Straight I-head, L-head, T-head, Hockey-head, and Circular heads
  • T-head and Hockey-head groynes provide superior sediment trapping and reduced scour at the nose
• Geometric Design Rules for Spurs
  • Length: Maximum length $\text{[PARSE ERROR: Undefined("Command("le")")]}1/5$ river width ($W$); minimum length $\text{[PARSE ERROR: Undefined("Command("ge")")]}2.5\times d_s$ (where Lacey scour depth $d_s=1.35(q^2/f{)}^{1/3}$)
  • Spacing: Typically spaced at $2.0\text{to}2.5$ times spur length
  • Layout: Array of 2 to 5 spurs in series along straight or convex banks

3. Guide Banks (Bell’s Bunds)

• Functions and Purpose
  • Confining wide alluvial braided flows to pass axially through bridge or barrage openings without outflanking
  • Preventing cross-flow and protecting approach embankments
• Geometric Design Standards
  • Upstream length: $1.0\text{to}1.5L$ (where $L$ is bridge waterway width); limits flow obliquity to $\text{[PARSE ERROR: Undefined("Command("le")")]}{30}^{\circ }$
  • Downstream length: $0.25\text{to}0.4L$
  • Upstream curved head: Radius $r_h=0.4\text{to}0.5L$; sweep angle $={120}^{\circ }\text{to}{145}^{\circ }$
  • Downstream curved tail: Radius $r_t=0.3\text{to}0.5r_h$; sweep angle $={45}^{\circ }\text{to}{60}^{\circ }$
  • Shank top width: $6\text{to}9\text{m}$ to accommodate vehicular maintenance access
• Launching Apron Mechanics
  • Heavy flexible stone layer placed horizontally on the bed; settles into developing scour holes to armor the slope against undermining
• Case Study: Guide bank design on the Arial Khan River at Madaripur Sadar

4. Permeable River Training Devices

• Concrete Porcupines
  • Permeable retarding devices constructed from reinforced concrete (RCC) or prestressed concrete (PSC) struts ($10\times 10\text{cm}$)
  • Tetrahedral (6 struts) and Prismatic (9 struts) configurations
  • Damps velocity, dissipates turbulence, and induces sediment accumulation along banklines
• Traditional Bamboo Bandaling
  • Traditional river training structures constructed from bamboo stakes and woven mats
  • Obstructs surface water flow (causing deposition near banks) while accelerating bottom currents (scouring a central navigation channel)
  • Case Study: Bandaling for bank protection and navigation on the Payra River at Lebukhali Cantonment

5. Revetments and Longitudinal Bank Protection

• Zonation of a River Bank
  • Embankment slope, Upper bank (between low water level and high flood level), and Lower underwater bank
• Revetment Structural Typology
  • Open-Joint (Flexible) Revetment: Unsealed joints allowing free movement of groundwater; relieves excess pore water pressure during rapid drawdown; ideal for alluvial rivers
  • Closed-Joint (Rigid) Revetment: Sealed joints with mortar or asphalt; susceptible to hydrostatic pressure buildup and blow-off failure; limited to upper banks above HWL
• Components of a Revetment System
  1. Cover / Armor Layer: Outer layer resisting shear forces, wave impact, and debris abrasion (rock riprap, CC blocks, geobags)
  2. Filter Layer (Granular or Geotextile): Prevents migration and leaching of fine subsoil particles while allowing groundwater drainage
  3. Toe Protection / Launching Apron: Trenching or flexible apron at the bank foot preventing toe undermining by scour
• Protective Armor Materials
  • Rock Riprap vs. Wire Gabions: Riprap requires mild slopes ($1:2$) and large stone ($>30\text{cm}$); gabions can be stacked on steep slopes using smaller stone ($8\mathrm{–}20\text{cm}$)
  • Precast Concrete Blocks: Cubical blocks ($300\text{mm}$ or $400\text{mm}$) manufactured on-site
  • Sand-Filled Geotextile Bags (Geobags): Highly flexible, self-adjusting revetments ($125\text{kg}$ and $250\text{kg}$ bags) widely deployed across Bangladesh

6. Nature-Based Solutions (NBS) and Bio-Engineering

• Principles of Green Infrastructure in River Engineering
  • Living, carbon-sequestering, biodegradable bank stabilization solutions
• Vegetative Species and Mechanics
  • Coastal/saline zones: Mangroves (Sundari, Kewra)
  • Freshwater zones: Vetiver grass (Chrysopogon zizanioides, deep fibrous roots up to $4\text{m}$) and Napier grass
  • Reinforcement mechanics: Roots bind soil matrix, increase apparent cohesion and shear strength, and attenuate wave energy
• Engineering Limitations of NBS
  • Vulnerable during initial growth phases before roots establish; unable to withstand immediate, severe, high-velocity hydraulic attacks

---

Lecture 12: River Navigation, Channel Regulation, Dredging Engineering, and Environmental Impacts

1. Inland Waterway Transport (IWT) Fundamentals and Economics

• Definition of Navigable Waters
  • Water bodies (rivers, canals, lakes) providing adequate depth, width, and flow calm for safe vessel transit in commerce
• Comparative Advantages of IWT
  • Lowest freight transport cost per ton-km
  • Highest energy efficiency and lowest fuel consumption
  • Lowest carbon emissions and air pollution
  • Land conservation (uses natural water bodies) and high safety record
  • Decongestion of overloaded road and rail networks
  • External environmental costs are 7 times lower than road transport
• Modal Transport Comparisons
  • Cost ratios (Water : Rail : Road): $1:2:5$ (international) / $1:3:8$ (China)
  • Bangladesh freight costs: Water (Tk 1.7) vs. Rail (Tk 5.6) vs. Road (Tk 6.5 per ton-km)
  • Hauling capability per horsepower: 1 hp pulls $4.0\text{tons}$ on water, $0.5\text{ton}$ on rail, and $0.15\text{ton}$ on road

2. Governing Factors of Channel Navigability

• Morphological Factors
  • Available Depth and Vessel Draft: Draft (keel to waterline); under-keel clearance; impact on fuel economy (fuel consumption drops by $25\mathrm{\%}$ when clearance increases from $0.6\text{m}$ to $1.5\text{m}$)
  • Channel Width: Fairway width for two-way vessel traffic; widening requirements at bends due to drift angles
  • Radius of Curvature at Bends: Minimum radius $r_c\text{[PARSE ERROR: Undefined("Command("ge")")]}3\times L_{\text{fleet}}$ (absolute minimum $2\times L_{\text{fleet}}$ to prevent speed loss)
  • Channel Cross-Section & Blockage Factor: Ratio of channel area to vessel submerged area ($A_c/A_v\text{[PARSE ERROR: Undefined("Command("ge")")]}6$) to prevent vessel squat, high resistance, and bank wash
• Flow Factors
  • Longitudinal current velocity (recommended maximum $\text{[PARSE ERROR: Undefined("Command("le")")]}3\text{m}s\text{[PARSE ERROR: Undefined("RBrace")]}$)
  • Transverse / helical currents (recommended maximum $\text{[PARSE ERROR: Undefined("Command("le")")]}0.3\text{m}s\text{[PARSE ERROR: Undefined("RBrace")]}$ to prevent lateral drift accidents)
  • Hazardous flow features: Shoals (submerged bars), rapids, eddies, vortices, and turbulent confluences

3. Engineering Methods for Navigability Improvement

• 1. Channel Regulation (Open River Training)
  • Use of spurs, groynes, training walls, and maintenance dredging to stabilize alignment and deepen the channel
• 2. Navigation Locks and Dams (Canalization of Rivers)
  • Construction of dams to transform rivers into a series of calm slack-water pools; locks transfer vessels between pools
  • Advantages (low velocities, wide waterways) vs. Disadvantages (high capital cost, lockage delays, rapid sedimentation in silt-laden rivers; e.g., Upper Mississippi River system)
• 3. Canalization (Bypass Navigation Canals)
  • Excavating artificial canals through land to bypass impassable river reaches or connect adjacent river basins
• 4. Artificial Meander Cut-offs
  • Dredging a pilot channel across meander necks to straighten rivers, shorten navigation paths, and improve flood discharge
  • Morphological Consequences: Local gradient steepening $⟹$ upstream nickpoint headcut erosion $⟹$ downstream sediment deposition/shoaling $⟹$ progressive relaxation toward a new equilibrium profile

4. Dredging Technology and Operational Phases

• Definition and the Dredging Cycle
  • Underwater excavation, hydraulic lifting, transport, and disposal of sediment
  • The four operational phases: Pretreatment, Extraction, Transportation, and Disposal / Placement (open-water, diked confinement, beneficial reuse)
• Dredger Classifications
  • Mechanical Dredgers: Bucket ladder, grab / clamshell, and backhoe / dipper dredgers
  • Hydraulic Dredgers: Plain suction dredger, Trailing Suction Hopper Dredger (TSHD), and Cutter Suction Dredger (CSD)
  • Cutter Suction Dredger (CSD): Stationary hydraulic dredger with rotating mechanical cutter head pivoting on a spud pole; excavates hard soil and pumps slurry through pipelines up to $10\text{km}$
• Operational Types
  • Capital / Development Dredging (initial deepening/widening) vs. Maintenance Dredging (routine removal of siltation)

5. Navigation Infrastructure and Dredging in Bangladesh

• National Waterway Network Overview
  • Total network of $24,000\text{km}$; monsoon navigable network: $5,968\text{km}$ ($25\mathrm{\%}$), shrinking to $3,865\text{km}$ ($16\mathrm{\%}$) in dry season
  • Traffic share: $12\mathrm{\%}$ of passenger traffic and $21\mathrm{\%}$ of national cargo traffic
• BIWTA Waterway Classification and Vertical Clearances
  • Class-I: Depth $3.66\mathrm{–}3.96\text{m}$, minimum vertical clearance $=18.30\text{m}$ from Standard High Water Level (SHWL)
    • Route 1: Chittagong – Chandpur – Narayanganj / Dhaka
    • Route 2: Shambhupura – Demra
    • Route 3: Shambhupura – Bhairab Bazar / Ashuganj
    • Route 4: Chowkighata – Barisal – Mongla – Khulna – Maheswarpasha
  • Class-II: Minimum vertical clearance $=12.20\text{m}$
  • Class-III: Minimum vertical clearance $=7.62\text{m}$
  • Class-IV: Minimum vertical clearance $=5.00\text{m}$
• Institutional Roles and Dredging Scale
  • Bangladesh Inland Water Transport Authority (BIWTA): Primarily maintenance dredging for navigation routes
  • Bangladesh Water Development Board (BWDB): Primarily capital dredging for flood control, drainage improvement, and river revival
  • Volume, costs, and fleet statistics

6. Morphological and Environmental Impacts of Dredging

• Channel Instability and the "Hungry Water" Effect
  • Deep dredging pits alter local equilibrium gradients and trap incoming sediment, starving downstream reaches
• Upstream Nickpoint Headcutting
  • Steeper gradient at the upstream lip of the dredge cut creates a nickpoint with high stream power, triggering headcut erosion propagating kilometers upstream and lowering river stages
• Downstream Incision and Bank Destabilization
  • Trapping of bedload inside the dredge cut causes clear-water downstream erosion, undermining riverbanks and bridge foundations
• The Sustainability Dilemma in Alluvial Rivers
  • In rivers with massive sediment transport, dredged cuts rapidly aggrade back to pre-dredge equilibrium profiles during monsoon floods
• Case Study: Gorai River Restoration Project (GRRP, Phase I, 1996–1999)
  • Massive pilot dredging ($200\text{m}$ base width, $13\text{m}$ depth) to augment dry-season flow; subsequent $1999$ flood completely silted up the dredged channel, illustrating that dredging alone without training structures is unsustainable

---