Hydrodynamic Numerical Methods
Computational techniques for solving shallow water equations, Navier-Stokes systems, and coastal circulation models
Overview
Hydrodynamic numerical methods discretize the governing equations of fluid motion — continuity, momentum, and transport equations — to simulate ocean circulation, tidal dynamics, storm surge, and pollutant dispersion. The shallow water equations (SWE), a depth-integrated form of the Navier-Stokes equations, form the backbone of most coastal and estuarine models. Three primary spatial discretization schemes (finite difference, finite element, finite volume) offer different trade-offs for accuracy, stability, and geometric flexibility. Models like ADCIRC, FVCOM, ROMS, Delft3D, and SCHISM are workhorses of operational and research oceanography.
Governing Equations
∂η/∂t + ∂(Hu)/∂x + ∂(Hv)/∂y = 0 (Continuity)
∂u/∂t + u·∂u/∂x + v·∂u/∂y - fv = -g·∂η/∂x + τˢˣ/(ρH) - τᵇˣ/(ρH) + Aₕ∇²u (x-Momentum)
∂v/∂t + u·∂v/∂x + v·∂v/∂y + fu = -g·∂η/∂y + τˢʸ/(ρH) - τᵇʸ/(ρH) + Aₕ∇²v (y-Momentum)
where η = water surface elevation, H = total depth, u,v = velocity components, f = Coriolis parameter, g = gravity
Key Concepts
Finite Difference (FD)
FD methods discretize derivatives on structured grids using Taylor series approximations. Arakawa B/C grids stagger variables. Simple to implement but limited to rectangular domains. Used in ROMS (terrain-following, curvilinear).
Finite Element (FE)
FE methods use unstructured triangular/tetrahedral meshes with basis functions. Excellent for complex coastlines and variable resolution. ADCIRC uses linear triangular elements; SCHISM uses mixed triangular-quadrilateral elements.
Finite Volume (FV)
FV methods conserve fluxes across control volumes. Naturally conservative, well-suited for wetting/drying. FVCOM uses unstructured triangular FV; Delft3D-FM uses flexible mesh approach.
Time Integration
Explicit methods (Euler, RK4) are simple but restricted by CFL condition (Δt ≤ Δx/c). Implicit methods (Crank-Nicolson) allow larger time steps. Semi-implicit schemes balance accuracy and efficiency for baroclinic-barotropic mode splitting.
Vertical Coordinates
Sigma (terrain-following) coordinates resolve bottom boundary layers; z-level coordinates suit deep ocean; isopycnal coordinates follow density surfaces. Hybrid approaches (HYCOM, SCHISM) combine advantages.
Model Validation
Skill metrics include RMSE, correlation coefficient, Willmott skill score, and bias for water level, velocity, and salinity. Taylor diagrams summarize multi-variable performance. Tidal constituent comparison (M2, S2 amplitude/phase) is standard.
Common Models
| Model | Method | Grid Type | Application |
|---|---|---|---|
| ADCIRC | FE (CG) | Unstructured triangles | Storm surge, tides |
| FVCOM | FV | Unstructured triangles | Estuaries, coastal |
| ROMS | FD | Structured curvilinear | Regional ocean |
| Delft3D | FD/FV | Structured/Flexible | Morphodynamics |
| SCHISM | FE (DG/CG) | Mixed tri/quad | Cross-scale ocean |
| TELEMAC | FE | Unstructured triangles | Rivers, coastal |
Interactive Visualizations
CFL Stability Analysis — Time Step vs. Grid Spacing
Simulated 2D Shallow Water — Surface Elevation Field
Convergence Analysis — Error vs. Grid Resolution
Key References
- Luettich, R.A., Westerink, J.J. & Scheffner, N.W. (1992). ADCIRC: An advanced three-dimensional circulation model for shelves, coasts, and estuaries. Technical Report DRP-92-6, US Army Corps of Engineers.
- Chen, C., Liu, H. & Beardsley, R.C. (2003). An unstructured grid, finite-volume, three-dimensional, primitive equations ocean model. J. Atmospheric and Oceanic Technology, 20, 159–186.
- Shchepetkin, A.F. & McWilliams, J.C. (2005). The regional oceanic modeling system (ROMS). Ocean Modelling, 9(4), 347–404.
- Zhang, Y.J. et al. (2016). Seamless cross-scale modeling with SCHISM. Ocean Modelling, 102, 64–81.
- Kantha, L.H. & Clayson, C.A. (2000). Numerical Models of Oceans and Oceanic Processes. Academic Press.