350 likes | 602 Views
AD aptive H ydraulics Modeling System Single Computational Engine for Navier -Stokes, Saturated Groundwater, and Shallow-water Equations Unstructured Platform Portability Multi-Physics Conservation Adaption. ADH Model.
E N D
ADaptiveHydraulics Modeling System • Single Computational Engine for Navier-Stokes, Saturated Groundwater, and Shallow-water Equations • Unstructured • Platform Portability • Multi-Physics • Conservation • Adaption ADH Model
Hydrodynamic model that engineers use routinely in rivers and estuaries. Development Philosophy • Navigation • Flood Control • Environmental Restoration
Structure Design GIWW Michoud Canal Structure Inset Bayou Bienvenue MRGO Proposed GIWW Structure Velocity Results (m/s)
Levee Breach Case Study New Orleans, USA
3D Hydraulics Case Study Lower Monumental Lock and Dam Washington, USA
3D Hydraulics Case Study 8100 ft of Snake River Reach 3300 ft wide @ widest section Structural Features: Lock Guard Wall, 8 Spillway Bays, & 6 Powerhouse Unit
Baroclinic Transport Mobile Bay Alabama, USA
Sediment Transport Bed Displacement Suspended Sediment
Research Versions: • 3D Hydrodynamics • Groundwater Interaction • Constituent Transport • Sediment Transport ADH Model Development Public Release v4.1: • 2D Hydrodynamics • Baroclinic Transport (Salinity & Temperature)
Variable time step: Time step determined by physics • Super and subcritical flow • Wetting and Drying: Shock Capturing • Adaption: Dynamic element refinement • Flow Control Structures: bridge decks, spillways, weirs, gates • Bendway correction: 3D effects of vorticity around bends ADH-2D Key Features
Mesh Construction and Mesh Quality • Defining material properties • Specifying Time and Model Control parameters • Assigning Boundary Conditions • Defining Control Structures • Setting Initial Conditions ADH Model Setup
Mesh Resolution 182 Elements 1200 Elements 4800 Elements 19200 Elements Initial Concentration
Mesh Resolution Computed Concentration @ 380 seconds
Adaption Computed Concentration Time Lapse with Adaption
Adaption Example Adaption around a moving vessel
Mesh−Element Quality Poor Element Quality Improved Element Quality
Element Quality Tips Less Desirable Better Preferred
Bed Roughness: • Manning’s n • Roughness height • Vegetation Eddy Viscosity: • Estimated • Kinematic Material Properties
Dynamic vs. Steady State • Automatic vs. Specified Time Steps • Output Interval Time Controls Initial time step guidance:
Subcritical • Inflow: Total Discharge, Unit Flow, or Vx & Vy • Outflow: Water Surface Elevation Supercritical • Inflow: Vx, Vy, and Depth • Outflow: Free flow Flow Boundary Conditions
Flow Control Structures Decks (Lid) Nodestrings & Nodes (Assign BC) Spillway Nodestrings (NB SPL Card) Weirs Nodestrings (WER & WRS Cards) Flap Gates Nodestrings (FLP & FGP Cards)
Rainfall and Evaporation Elements (Material Properties) Wind Friction Stations (Spatial Data Coverage) Ice Friction Elements (ICE & IRH cards) Sidewall Friction Nodestrings (Assign BC) Other Boundary Conditions
Constant WSE recommended for cold start • Include velocity only when hotstarting from an ADH solution • Major contributing factor for model success Initial Conditions Guidance
ADH Model Execution • Pre-ADH • ADH • PC LVL card for additional screen output Nonlinear Iteration Max Residual Node ID Failure Flag # Nodes Physics Time Step Percent Completion Linear Iteration Max Increment Norm Node ID Time (sec)
Build Mesh • Mesh Quality • Material Properties • Time Control Parameters • Boundary Conditions • Initial Conditions • Renumber • Save Project & Launch ADH ADH Checklist
The Modeling Process Data Collection Network Design Model Application Model Debugging Testing & Calibration F e e d b a c k L o o p
Mesh geometry • Voids? • Poorly formed elements? • Reasonable resolution? • Are the elements aligned with flow direction and wetting/drying boundary? Troubleshooting
Boundary conditions • Is the water level below elevation of nodes? • Are your initial conditions consistent with your boundary conditions? • Supercritical flow? • Is the initial depth zero for a significant portion of the domain? Troubleshooting
Model parameters • Adjust DTL (shock capturing parameter) – rule of thumb is 5% of average depth. • Adjust time steps • Try switching the tolerance used to determine convergence (NTL/ITL) • Increase the maximum number of iterations (don’t go higher than 10) Troubleshooting
Try running a steady state simulation (STD card) and using the resulting depths/velocities as a hotstart for the dynamic run • Check the solver (in .bc file) if you are using multiple processors or blocks • Try allowing ADH to automatically compute time steps (ATF card) Other Troubleshooting Tips