COMPUTATIONAL MODELS FOR WEAKLY DISPERSIVE NONLINEAR WATER-WAVES

Numerical methods for the two- and three-dimensional Boussinesq equati ons governing weakly nonlinear and dispersive water waves are presente d and investigated. Convenient handling of grids adapted to the geomet ry or bottom topography is enabled by finite element discretization in space. Staggered finite difference schemes are used for the temporal discretization. resulting in only two linear systems to be solved duri ng each time step. Efficient iterative solution of linear systems is d iscussed. By introducing correction terms in the equations, a fourth-o rder, two-level temporal scheme can be obtained. Combined with (bi-) q uadratic finite elements, the truncation errors of this scheme can be made of the same order as the neglected perturbation terms in the anal ytical model, provided that the element size is of the same order as t he characteristic depth. We present analysis of the proposed schemes i n terms of numerical dispersion relations. Verification of the schemes and their implementations is performed for standing waves in a closed basin with constant depth. More challenging applications cover plane incoming waves on a curved beach and earthquake induced waves over a s hallow seamount. In the latter example we demonstrate a significantly increased computational efficiency when using higher-order schemes and bathymetry-adapted finite element grids. (C) 1998 Elsevier Science S. A. All rights reserved.