Stability analysis of numerical methods for a 1.5-layer shallow-water ocean model (Q1789943)

Summary: A 1.5-layer reduced-gravity shallow-water ocean model in spherical coordinates is described and discretized in a staggered grid (standard Arakawa C-grid) with the forward-time central-space (FTCS) method and the Leap-frog finite difference scheme. The discrete Fourier analysis method combined with the Gershgorin circle theorem is used to study the stability of these two finite difference numerical models. A series of necessary conditions of selection criteria for the time-space step sizes and model parameters are obtained. It is showed that these stability conditions are more accurate than the Courant-Friedrichs-Lewy (CFL) condition and other two criterions [\textit{A. F. Blumberg} and \textit{G. L. Mellor}, ``A description of a three-dimensional coastal ocean model'', in: Three-Dimensional Coastal Ocean Models, Coastal and Estuarine Sciences 32. 1--16 (1987); \textit{V. Casulli}, J. Comput. Phys. 86, No. 1, 56--74 (1990; Zbl 0681.76022); \textit{V. Casulli} and \textit{R. T. Cheng} [Int. J. Numer. Methods Fluids 15, No. 6, 629--648 (1992; Zbl 0762.76068)]. Numerical experiments are proposed to test our stability results, and numerical model that is designed is also used to simulate the ocean current.