Dispersion and stability analyses of the linearized two-dimensional shallow water equations in Cartesian coordinates

Document Type


Date of Original Version



In the present study, a Fourier analysis is used to develop expressions for phase and group speeds for both continuous and discretized, linearized two-dimensional shallow water equations, in Cartesian coordinates. The phase and group speeds of the discrete equations, discretized using a three-point scheme of second order, five-point scheme of fourth order and a three-point compact scheme of fourth order in an Arakawa C grid, are calculated and compared with the corresponding values obtained for the continuous system. The three-point second-order scheme is found to be non-dispersive with grid resolutions greater than 30 grids per wavelength, while both the fourth-order schemes are non-dispersive with grid resolutions greater than six grids per wavelength. A von Neumann stability analysis of the two- and three-time-level temporal schemes showed that both schemes are stable. A wave deformation analysis of the two-time-level Crank-Nicolson scheme for one-dimensional and two-dimensional systems of shallow water equations shows that the scheme is non- dispersive, independent of the Courant number and grid resolution used. The phase error or the dispersion of the scheme decreases with a decrease in the time step or an increase in grid resolution. © 2003 Elsevier Ltd. All rights reserved.

Publication Title, e.g., Journal

Ocean Engineering