A high‐order Petrov–Galerkin finite element method for the classical Boussinesq wave model

Abstract

SUMMARYA high-order Petrov-Galerkin finite element scheme is presented to solve the one-dimensional depthintegrated classical Boussinesq equations for weakly non-linear and weakly dispersive waves. Finite elements are used both in the space and the time domains. The shape functions are bilinear in spacetime, whereas the weighting functions are linear in space and quadratic in time, with C 0 -continuity. Dispersion correction and a highly selective dissipation mechanism are introduced through additional streamline upwind terms in the weighting functions. An implicit, conditionally stable, one-step predictorcorrector time integration scheme results. The accuracy and stability of the non-linear discrete equations are investigated by means of a local Taylor series expansion. A linear spectral analysis is used for the full characterization of the predictor-corrector inner iterations. Based on the order of the analytical terms of the Boussinesq model and on the order of the numerical discretization, it is concluded that the scheme is fourth-order accurate in terms of phase velocity. The dissipation term is third order only affecting the shortest wavelengths. A numerical convergence analysis showed a second-order convergence rate in terms of both element size and time step. Four numerical experiments are addressed and their results are compared with analytical solutions or experimental data available in the literature: the propagation of a solitary wave, the oscillation of a flat bottom closed basin, the oscillation of a non-flat bottom closed basin, and the propagation of a periodic wave over a submerged bar.

Finite volume schemes for dispersive wave propagation and runup

“…Finite difference (FD) schemes [39,43,47], finite element methods [7,15,28,58] and spectral methods [60,62] have been proposed. More contemporary discontinuous Galerkin (DG) schemes have also been adapted with some success to dispersive wave equations [33,34,53,79] while the application of Finite Volumes (FV) or hybrid FV/FD methods remain most infrequent for this type of problems.…”

Numerical Modeling Tools Applied to Estuarine and Coastal Hydrodynamics: A User Perspective

“…And knowledge about the hydrodynamic patterns influenced by bottom morphological changes may allow assessment and forecasting of the effects of hazardous and extreme events, anthropogenic intervention, or climate change. Hydrodynamic modeling has therefore been the focus of a large number of previous works in estuarine environments [11,14,17,21,22,[40][41][42][43][44][45][46][47][48][49][50].…”

Finite volume and pseudo-spectral schemes for the fully nonlinear 1D Serre equations

“…3) 4) with φ being the velocity potential (by definition, the irrotational velocity field (u, v) = (∇φ, ∂ y φ), g the acceleration due to the gravity force and ∇ = (∂ x 1 , ∂ x 2 ) denotes the gradient operator in horizontal Cartesian coordinates and |∇φ| 2 ≡ (∇φ) · (∇φ). The incompressibility condition leads to the Laplace equation for φ.…”