Stability Analysis of Geostrophic Adjustment on Hexagonal Grids for Regions with Variable Depth

Hexagonal grids have been used in a number of numerical studies, and especially in relation to atmospheric models. Recent studies have suggested that ocean circulation models may also benefit from the use of hexagonal grids. These grids tend to induce less systematic errors and have better horizontal isotropy properties than traditional square grid schemes. If hexagonal grids are to be applied in ocean models, a number of features that are characteristic of ocean circulation problems need to be attended to. The topography of the ocean basin is an important feature in most ocean models. Ocean modelers can experience instabilities due to depth variations. In the present paper, analysis of the propagation matrix for the spatially discretized system is used to explain unphysical growth of the numerical solutions of the linear shallow water equations when using hexagonal grids over domains with variable depth. It is shown that a suitable weighting of the Coriolis terms may give an energy-conserving and stable numerical scheme.

[1]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

[2]  S. Nickovic,et al.  Frequency of quasi-geostrophic modes on hexagonal grids , 1996 .

[3]  P. Morel,et al.  A FINITE-DIFFERENCE APPROXIMATION OF THE PRIMITIVE EQUATIONS FOR A HEXAGONAL GRID ON A PLANE , 1969 .

[4]  David L. Williamson,et al.  Integration of the barotropic vorticity equation on a spherical geodesic grid , 1968 .

[5]  Y. Sasaki.  Variational design of finite-difference schemes for initial value problems with an integral invariant , 1976 .

[6]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[7]  Slobodan Nickovic,et al.  Geostrophic Adjustment on Hexagonal Grids , 2002 .

[8]  B. Ådlandsvik,et al.  Production of monthly mean climatological archives for the Nordic Seas , 1998 .

[9]  U. Feudel,et al.  Grid geometry effects on convection in ocean climate models: a conceptual study , 2004 .

[10]  D. Randall,et al.  A Potential Enstrophy and Energy Conserving Numerical Scheme for Solution of the Shallow-Water Equations on a Geodesic Grid , 2002 .

[11]  J. M. Popovi.  Frequency of Quasi-Geostrophic Modes on Hexagonal Grids , .

[12]  Todd D. Ringler,et al.  Climate modeling with spherical geodesic grids , 2002, Comput. Sci. Eng..

[13]  Richard Asselin,et al.  Frequency Filter for Time Integrations , 1972 .

[14]  M. Cullen,et al.  Numerical Prediction and Dynamic Meteorology, 2nd Edn. By G. J. HALTINER and R. T. WILLIAMS. Wiley, 1980. 477 pp. £26.90. , 1984, Journal of Fluid Mechanics.

[15]  R. Sadourny,et al.  Compressible Model Flows on the Sphere. , 1975 .

[16]  Conservation of energy for schemes applied to the propagation of shallow‐water inertia‐gravity waves in regions with varying depth , 2000 .

[17]  David A. Randall,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid. Part II. A Detailed Description of the Grid and an Analysis of Numerical Accuracy , 1995 .

[18]  On the use of hexagonal grids for simulation of atmospheric processes , 1994 .

[19]  J. Berntsen,et al.  On a class of ocean model instabilities that may occur when applying small time steps, implicit methods, and low viscosities , 2004 .

[20]  Ionel Michael Navon,et al.  Implementation of A Posteriori Methods for Enforcing Conservation of Potential Enstrophy and Mass in Discretized Shallow-Water Equations Models , 1981 .

[21]  Dale B. Haidvogel,et al.  Numerical Ocean Circulation Modeling , 1999 .

[22]  A. Arakawa.  Computational design for long-term numerical integration of the equations of fluid motion: two-dimen , 1997 .

[23]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[24]  L. Kantha,et al.  Numerical models of oceans and oceanic processes , 2000 .

[25]  R. Sadourny.  The Dynamics of Finite-Difference Models of the Shallow-Water Equations , 1975 .