An unstructured grid morphodynamic model with a discontinuous Galerkin method for bed evolution

A new unstructured grid two-dimensional, depth-integrated (2DDI), morphodynamic model is presented for the prediction of morphological evolutions in shallow water. This modelling system consists of two coupled model components: (i) a well-verified and validated continuous Galerkin (CG) finite element hydrodynamic model; and (ii) a new sediment transport/bed evolution model that uses a discontinuous Galerkin (DG) method for the solution of the sediment continuity equation. The DG method is a robust finite element method that is particularly well suited for this type of advection dominated transport equation. It incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep bathymetric gradients that may form in the solution, and it possesses a local conservation property that conserves sediment mass on an elemental level. In this paper, we focus specifically on the implementation and verification of the DG model. Details are given on the implementation of the method, and numerical results are presented for three idealized test cases which demonstrate the accuracy and robustness of the method and its applicability in predicting medium-term morphological changes in channels and coastal inlets. � 2005 Elsevier Ltd. All rights reserved.

[1]  Ingemar Kinnmark,et al.  The Shallow Water Wave Equations: Formulation, Analysis and Application , 1985 .

[2]  Clint Dawson,et al.  A New Generation Hurricane Storm Surge Model for Southern Louisiana , 2005 .

[3]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[4]  G. Karami.  Lecture Notes in Engineering , 1989 .

[5]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[6]  J. F. A. Sleath,et al.  Sea bed mechanics , 1984 .

[7]  Joannes J. Westerink,et al.  The influence of domain size on the response characteristics of a hurricane storm surge model , 1994 .

[8]  Mary F. Wheeler,et al.  Compatible algorithms for coupled flow and transport , 2004 .

[9]  Magnus Larson,et al.  A general formula for non-cohesive bed load sediment transport , 2005 .

[10]  R. Soulsby,et al.  Threshold of Sediment Motion in Coastal Environments , 1997 .

[11]  P. Raviart,et al.  Conforming and nonconforming finite element methods for solving the stationary Stokes equations I , 1973 .

[12]  M. Hayes.  General morphology and sediment patterns in tidal inlets , 1980 .

[13]  J. A. Roelvink,et al.  Medium-term 2DH coastal area modelling , 1993 .

[14]  R. Luettich,et al.  Modelling tides in the western North Atlantic using unstructured graded grids , 1994 .

[15]  W. Graf.  Hydraulics of Sediment Transport , 1984 .

[16]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[17]  F. Henderson.  Open channel flow , 1966 .

[18]  M. Iskandarani,et al.  Comparison of advection schemes for high-order h–p finite element and finite volume methods , 2005 .

[19]  Ida Brøker,et al.  Intercomparison of coastal area morphodynamic models , 1997 .

[20]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[21]  Ethan J. Kubatko,et al.  hp Discontinuous Galerkin methods for advection dominated problems in shallow water flow , 2006 .

[22]  William G. Gray,et al.  A wave equation model for finite element tidal computations , 1979 .

[23]  Bernardo Cockburn,et al.  The Runge-Kutta local projection P1-discontinuous-Galerkin finite element method for scalar conservation laws , 1988 .

[24]  R. Luettich,et al.  Eastcoast 2001, A Tidal Constituent Database for Western North Atlantic, Gulf of Mexico, and Caribbean Sea , 2002 .

[25]  Julio A. Zyserman,et al.  Controlling spatial oscillations in bed level update schemes , 2002 .