needed leads to improved accuracy for a given computational cost as compared to models that use
structured grid methods. However, both structured and unstructured grid method solutions to the
governing morphological equation can experience numerical robustness and accuracy problems
manifested in the form of spurious spatial oscillations, especially in the presence of steep
bathymetric gradients (see for example Johnson and Zyserman, 2002).
In this paper, we describe the development of a new unstructured grid morphodynamic
model system that uses a new class of highly accurate finite element methods for the solution of
the governing morphological equation. The hydrodynamic model component of our system is
provided by the well verified and validated unstructured grid model ADCIRC, developed by the
second author and a number of collaborators (Luettich and Westerink, 2004). ADCIRC is both a
two-dimensional, depth-integrated (2DDI) and three-dimensional (3D) free surface flow model.
In this paper, we focus specifically on the 2DDI ADCIRC model, which solves the shallow water
equations using the standard or continuous Galerkin (CG) finite element method in space. To
overcome well known problems in solving the shallow water equations using equal-order
interpolating spaces with the CG finite element method, the continuity equation is replaced by the
so-called generalized wave continuity equation (GWCE) (Lynch and Gray, 1979 and Kinnmark,
1986). The solution strategy used in ADCIRC has proven to be robust and computationally
efficient, and it has been validated in a large number of cases (see for example Blain, et al., 1994;
Westerink et al., 1994; Mukai, et al., 2002; Westerink et al., 2004).
Working with a well established hydrodynamic model then, the main focus of this paper
is the development and verification of a bed load sediment transport/morphological model
component to work in conjunction with ADCIRC. Mathematically, the morphological evolution
of the bed is defined by the so-called sediment continuity or Exner equation. This equation
simply states that the time rate of change of the bed elevation is equal to the divergence of the
sediment flux, which can be expressed in terms of the local flow properties through the use of an
empirical sediment transport formula. As is well known, solving advection dominated transport
2
Previous Page
