Dispersion and stability analyses of the linearized two-dimensional shallow water equations in boundary-fitted co-ordinates

Document Type


Date of Original Version



In the present investigation, a Fourier analysis is used to study the phase and group speeds of a linearized, two-dimensional shallow water equations, in a non-orthogonal boundary-fitted co-ordinate system. The phase and group speeds for the spatially discretized equations, using the second-order scheme in an Arakawa C grid, are calculated for grids with varying degrees of non-orthogonality and compared with those obtained from the continuous case. The spatially discrete system is seen to be slightly dispersive, with the degree of dispersivity increasing with an decrease in the grid non-orthogonality angle or decrease in grid resolution and this is in agreement with the conclusions reached by Sankaranarayanan and Spaulding (J. Comput. Phys., 2003; 184:299-320). The stability condition for the non-orthogonal case is satisfied even when the grid non-orthogonality angle, is as low as 30° for the Crank Nicolson and three-time level schemes. A two-dimensional wave deformation analysis, based on complex propagation factor developed by Leendertse (Report RM-5294-PR, The Rand Corp., Santa Monica, CA, 1967), is used to estimate the amplitude and phase errors of the two-time level Crank-Nicolson scheme. There is no dissipation in the amplitude of the solution. However, the phase error is found to increase, as the grid angle decreases for a constant Courant number, and increases as Courant number increases. © 2003 John Wiley and Sons, Ltd.

Publication Title

International Journal for Numerical Methods in Fluids