Journal article
Linear and nonlinear Stability analysis for finite difference discretizations of higher order Boussinesq equations
Coastal, Maritime and Structural Engineering, Department of Mechanical Engineering, Technical University of Denmark1
Department of Mechanical Engineering, Technical University of Denmark2
Scientific Computing, Department of Informatics and Mathematical Modeling, Technical University of Denmark3
Department of Informatics and Mathematical Modeling, Technical University of Denmark4
This paper considers a method of lines stability analysis for finite difference discretizations of a recently published Boussinesq method for the study of highly nonlinear 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 nonlinear 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 nonlinearity 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 nonlinear analysis.
The various methods of analysis combine to provide significant insight into into the numerical behavior of this rather complicated system of nonlinear PDEs.
Language: | English |
---|---|
Year: | 2004 |
Pages: | 751-773 |
ISSN: | 10970363 and 02712091 |
Types: | Journal article |
ORCIDs: | Bingham, Harry B. and Madsen, Per A. |