250 likes | 255 Views
Transient Mixed Flow Modeling Capability in SRH-2D for Culverts and Bridges Yong G. Lai. Background Equations/Numerical Methods Verification/Demonstration Concluding Remarks. Outline. Transition between free surface open channel flows and pressurized conduit flows
E N D
Transient Mixed Flow Modeling Capability in SRH-2D for Culverts and Bridges Yong G. Lai
Background Equations/Numerical Methods Verification/Demonstration Concluding Remarks Outline
Transition between free surface open channel flows and pressurized conduit flows Examples: Culvert; Bridge; Sewer/Storm Water System & other closed conduit system What is Mixed Flow?
Large body of literature A recent review: Bousso et al. (2013), J. Hydraulic Eng 139(4) Challenges Two sets of equations with fast transient flows Discontinuous interface between free and pressurized flows Complex wave propagation Numerical Methods
Two Major Approaches One-Equation (Preissmann Slot) Two-Equation (Multi-Phase) Primarily 1D Models Suitable for most applications But 2D needed for, e.g., flooding with culverts/bridges, curved conduit etc. Existing Numerical Methods
A 2D Mixed Flow Model Few 2D Models: Maranzoni et al. (2015) (Advances in Water Resources) Built on SRH-2D (1D used at present) Capable of fast transient modeling? Preissmann Slot Method Extension to unstructured mesh Present Objective
Advantages: Pressurized zones reduce to open-channel equations “Slot” simulates fluid compressibility & conduit deformation Limitations: “False Volumes” added in the slot introduce solution errors! Slot width needs be minimized for accuracy But increased slot width reduces oscillations around shock waves Key Advantages & Limitations
Finite-Volume Collocated Method Pressure Head as Primary Variables (vs Density-Based) Implicit Time Integration SIMPLE-C Algorithm Sub-, Super-, and Trans-Critical Flows Arbitrary Shaped Mesh Cells Wetting-Drying Algorithm Numerical Algorithm
Left: Head=0.8 m Velocity= 2.0 m/s Right: Head=0.8 m Velocity=-2.0 m/s Ceiling: 1.0 m Slot_Ratio = 0.005 Test 1: Colliding Flows
Exact Solution: From Rankine-Hugoniot Relation and Exact Riemann Solver (Solid) Slower wave speed (due to slot introduced) Head Oscillation near shock is visible (symbols) Test 1: Comparison
Left: Head=3.0 m Velocity= 0.0 m/s Right: Head=0.5 m Velocity= 0.0 m/s Ceiling: 1.0 m What to expect? Right-traveling shock; Left-traveling rarefaction wave for depressurization Test 2: Dam-Break
Solid: Exact Solution at time=0.3 s Symbol: SRH-2D Slot_Ratio: 0.005 Test 2: Comparison
Left: Head=10 m Velocity= 0.0 m/s Right: Head= 1 m Velocity= 0.0 m/s Ceiling: 5.0 m Pressure Zone Radius: 11 m What to expect? Right-traveling shock; Left-traveling rarefaction wave for depressurization Model: Triangular Mesh Test 3: 2D Circular Dam-Break
Solid: SRH-2D Solution from 2D Triangular Mesh Dash: Reference “Exact” Solution from 1D Radial Solution Test 3: Comparison
Dam Water Depth: 15 m; Channel Water Depth: 1.5m Demonstration: Dam-Break Flood through a Bridge
Dam Water Depth: 15 m; Channel Water Depth: 1.5m Results at Time = 10 s
Preissmann Slot Concept is extended to 2D models with arbitrary mesh cells Equations are derived 2D transient mixed flow capability is successfully developed into SRH-2D Test cases demonstrate the accuracy of the model implementation Concluding Remarks
THANK YOU QUESTIONS ?