Premium
A semi‐implicit lagrangian scheme for 3D shallow water flow with a two‐layer turbulence model
Author(s) -
Stansby Peter K.,
Lloyd Peter M.
Publication year - 1995
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1650200203
Subject(s) - mechanics , turbulence , boundary layer , turbulence modeling , flow (mathematics) , vortex , geometry , mixing (physics) , mathematics , reynolds stress , geology , physics , quantum mechanics
Abstract A semi‐implicit Lagrangian finite difference scheme for 3D shallow water flow has been developed to include an eddy viscosity model for turbulent mixing in the vertical direction. The α‐co‐ordinate system for the vertical direction has been introduced to give accurate definition of bed and surface boundary conditions. The simple two‐layer mixing length model for rough surfaces is used with the standard assumption that the shear stress across the wall region at a given horizontal location is constant. The bed condition is thus defined only by its roughness height (avoiding the need for a friction formula relating to depth‐averaged flow, e.g. Chezy, used previously). The method is shown to be efficient and stable with an explicit Lagrangian formulation for convective terms and terms for surface elevation and vertical mixing handled implicitly. The method is applied to current flow around a circular island with gently sloping sides which produce periodic recirculation zones (vortex shedding). Comparisons are made with experimental measurements of velocity using laser Doppler anemometry (time histories at specific points) and surface particle‐tracking velocimetry.