Instabilities of Boussinesq models in non-uniform depth

The von Neumann method for stability analysis of linear waves in a uniform medium is a widely applied procedure. However, the method does not apply to stability of linear waves in a variable medium. Herein we describe instabilities due to variable depth for different Boussinesq equations, including the standard model by Peregrine and the popular generalization by Nwogu. Eigenmodes are first found for bathymetric features on the grid scale. For certain combinations of Boussinesq formulations and bottom profiles stability limits are found in closed form, otherwise numerical techniques are used for the eigenvalue problems. Naturally, the unstable modes in such settings must be considered to be as much a result of the difference method as of the underlying differential (Boussinesq) equations. Hence, modes are also computed for smooth depth profiles that are well resolved. Generally, the instabilities do not vanish with refined resolution. In some cases convergence is observed and we thus have indications of unstable solutions of the differential equations themselves. The stability properties differ strongly. While the standard Boussinesq equations seem perfectly stable, all the other formulations do display unstable modes. In most cases the instabilities are linked to steep bottom gradients and small grid increments. However, while a certain formulation, based on velocity potentials, is very prone to instability, the Boussinesq equations of Nwogu become unstable only under quite demanding conditions. Still, for the formulation of Nwogu, instabilities are probably inherent in the differential equations and are not a result of the numerical model. Copyright (C) 2008 John Wiley & Sons, Ltd.