| Abstract: | This paper considers a method of lines stability analysis for finite difference discretizations of a recently published Boussinesq method for the study of highly non‐linear and extremely dispersive water waves. The analysis demonstrates the near‐equivalence of classical linear Fourier (von Neumann) techniques with matrix‐based methods for formulations in both one and two horizontal dimensions. The matrix‐based method is also extended to show the local de‐stabilizing effects of the non‐linear terms, as well as the stabilizing effects of numerical dissipation. A comparison of the relative stability of rotational and irrotational formulations in two horizontal dimensions provides evidence that the irrotational formulation has significantly better stability properties when the deep‐water non‐linearity is high, particularly on refined grids. Computation of matrix pseudospectra shows that the system is only moderately non‐normal, suggesting that the eigenvalues are likely suitable for analysis purposes. Numerical experiments demonstrate excellent agreement with the linear analysis, and good qualitative agreement with the local non‐linear analysis. The various methods of analysis combine to provide significant insight into the numerical behaviour of this rather complicated system of non‐linear PDEs. Copyright © 2004 John Wiley & Sons, Ltd. |
