Boundary conditions in shallow water models - An alternative implementation for finite element codes

Finite element solution of the shallow water wave equations has found increasing use by researchers and practitioners in the modelling of oceans and coastal areas. Wave equation models, most of which use equal-order C 0 interpolants for both the velocity and the surface elevation, do not introduce spurious oscillation modes, hence avoiding the need for artificial or numerical damping. An important question for both primitive equation and wave equation models is the interpretation of boundary conditions. Analysis of the characteristics of the governing equations shows that for most geophysical flows a single condition at each boundary is sufficient, yet there is not a consensus in the literature as to what that boundary condition must be or how it should be implemented in a finite element code. Traditionally (partly because of limited data), surface elevation is specified at open ocean boundaries while the normal flux is specified as zero at land boundaries. In most finite element wave equation models both of these boundary conditions are implemented as essential conditions. Our recent work focuses on alternative ways to numerically implement normal flow boundary conditions with an eye towards improving the mass-conserving properties of wave equation models. A unique finite element formulation using generalized functions demonstrates that boundary conditions should be implemented by treating normal fluxes as natural conditions with the flux interpreted as external to the computational domain. Results from extensive numerical experiments show that the scheme does conserve mass for all parameter values. Furthermore, convergence studies demonstrate that the algorithm is consistent, as residual errors at the boundary diminish as the grid is refined.

[1]  M. Foreman,et al.  An accuracy analysis of boundary conditions for the forced shallow water equations , 1986 .

[2]  W. Gray,et al.  NORMAL FLOW BOUNDARY CONDITIONS IN SHALLOW WATER MODELS - INFLUENCE ON MASS CONSERVATION AND ACCURACY , 1994 .

[3]  G. K. Verboom,et al.  Weakly-reflective boundary conditions for two-dimensional shallow water flow problems , 1984 .

[4]  William G. Gray,et al.  On the well posedness of some wave formulations of the shallow water equations , 1988 .

[5]  W. Gray,et al.  Quiet: A reduced noise finite element model for tidal circulation , 1983 .

[6]  Robert L. Lee,et al.  Don''t suppress the wiggles|they''re telling you something! Computers and Fluids , 1981 .

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

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

[9]  R. Luettich,et al.  THE INFLUENCE OF NORMAL FLOW BOUNDARY CONDITIONS ON SPURIOUS MODES IN FINITE ELEMENT SOLUTIONS TO THE SHALLOW WATER EQUATIONS , 1994 .

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

[11]  S. Sigurdsson.  Alternative Approaches to the Wave Equation in Shallow Water Flow Modelling , 1994 .

[12]  R. Luettich,et al.  Development of the direct stress solution technique for three‐dimensional hydrodynamic models using finite elements , 1994 .

[13]  Joannes J. Westerink,et al.  Tide and Storm Surge Predictions Using Finite Element Model , 1992 .

[14]  William G. Gray,et al.  Some inadequacies of finite element models as simulators of two-dimensional circulation. , 1982 .

[15]  Joannes J. Westerink,et al.  Aspects of nonlinear simulations using shallow-water models based on the wave continuity equation , 1994 .

[16]  A Comparison of Tidal Models for the Southwest Coast of Vancouver Island , 1988 .

[17]  Nicole Goutal,et al.  TELEMAC: A new numerical model for solving shallow water equations , 1991 .

[18]  Roy A. Walters,et al.  Accuracy of an estuarine hydrodynamic model using smooth elements , 1980 .

[19]  W. Gray.  A finite element study of tidal flow data for the North Sea and English Channel , 1989 .

[20]  Norman W. Scheffner,et al.  ADCIRC: An Advanced Three-Dimensional Circulation Model for Shelves, Coasts, and Estuaries. Report 1. Theory and Methodology of ADCIRC-2DDI and ADCIRC-3DL. , 1992 .

[21]  C. Taylor,et al.  Tidal and long wave propagation—a finite element approach , 1975 .

[22]  Joannes J. Westerink,et al.  A solution for the vertical variation of stress, rather than velocity, in a three-dimensional circulation model , 1991 .

[23]  Daniel R. Lynch.  Mass balance in shallow water simulation , 1985 .