Semi‐implicit finite difference methods for three‐dimensional shallow water flow

Abstract

SUMMARYA semi-implicit finite difference method for the numerical solution of three-dimensional shallow water flows is presented and discussed. The governing equations are the primitive three-dimensional turbulent mean flow equations where the pressure distribution in the vertical has been assumed to be hydrostatic. In the method of solution a minimal degree of implicitness has been adopted in such a fashion that the resulting algorithm is stable and gives a maximal computational efficiency at a minimal computational cost. At each time step the numerical method requires the solution of one large linear system which can be formally decomposed into a set of small three-diagonal systems coupled with one five-diagonal system. All these linear systems are symmetric and positive definite. Thus the existence and uniqueness of the numerical solution are assured. When only one vertical layer is specified, this method reduces as a special case to a semi-implicit scheme for solving the corresponding two-dimensional shallow water equations. The resulting two-and three-dimensional algorithm has been shown to be fast, accurate and mass-conservative and can also be applied to simulate flooding and drying of tidal mud-flats in conjunction with threedimensional flows. Furthermore, the resulting algorithm is fully vectorizable for an efficient implementation on modern vector computers.KEY WORDS Three-dimensional Semi-implicit Shallow water

The effects of river inflow and retention time on the spatial heterogeneity of chlorophyll and water–air CO<sub>2</sub> fluxes in a tropical hydropower reservoir

“…This model solves the 3-D hydrostatic, Boussinesq, Reynolds-averaged Navier-Stokes and scalar transport equations, separating mixing of scalars and momentum from advection. The hydrodynamic algorithms that are implemented in the ELCOM use a EulerLagrange approach for the advection of momentum adapted from the work of Casulli and Cheng (1992), whereas the advection of scalars (i.e., tracers, conductivity and temperature) is based on the ULTIMATE QUICKEST method proposed by Leonard (1991). The thermodynamics model considers the penetrative (i.e., shortwave radiation) and nonpenetrative components (i.e., longwave radiation, sensible and latent heat fluxes) (Hodges et al, 2000).…”

Challenges and opportunities for integrating lake ecosystem modelling approaches

“…The hydrodynamic model is based on the Navier-Stokes equations that describe dynamically the 3-D transport and mixing of water. An efficient numerical semiimplicit Eulerian-Lagrangian finite difference scheme was used in order to assure stability, convergence and accuracy (Casulli and Cheng 1992;Casulli and Cattani 1994). The chemical and biological dynamics in the open water and in the sediment are based on PCLake, with an adaptation to subtropical lakes by including omnivorous fish.…”

Physical changes within a large tropical hydroelectric reservoir induced by wintertime cold front activity

“…ELCOM solves unsteady, coupled, Reynolds-averaged Navier-Stokes equations and scalar transport equations using hydrostatic and Boussinesq approximations to model the velocity, salinity and temperature of water in 3-D space and time. The hydrodynamic algorithms that are implemented in ELCOM use an Euler-Lagrange approach for the advection of momentum adapted from Casulli and Cheng (1992), while the advection of scalars is based on the ULTIMATE QUICKEST method proposed by Leonard (1991).…”