An optimization of the Icosahedral grid modified by spring dynamics

We have investigated an optimum form of the modified icosahedral grid that is generated by applying the spring dynamics to the standard icosahedral grid system. The spring dynamics can generate a more homogeneous grid system than the standard icosahedral grid system by tuning the natural spring lenght: as the natural spring length becomes longer, the ratio of maximum grid interval to minimum one becomes closer to unit. When the natural spring length is larger than a critical value, however, the spring dynamic system does not have a stable equilibrium. By setting the natural spring length to be the marginally critical value, we can obtain the most homogeneous grid system, which is most efficient in terms of the CFL condition. We have analyzed eigenmodes involved in the initial error of the geostrophic balance problem [test case 2 of D. L. Williamson et al. (1992, J. Comput. Phys. 102, 211)]. Since the balance state in the discrete system differs slightly from the exact solution of the analytic system, the initial error field includes both the gravity wave mode and the Rossby wave mode. As the results of the analysis are based on Hough harmonics decompositions, we detected Rossby and gravity wave modes with zonal wavenumber 5, which are asymmetric against the equator. These errors are associated with icosahedral grid structure. The symmetric gravity wave mode with zonal wavenumber 0 also appears in the error field. To clarify the evolution of Rossby waves, we introduce divergence damping to reduce the gravity wave mode. From the simulated results of the geostrophic problem with various grid systems, we found that the spuriously generated Rossby wave mode is eliminated most effectively when the most homogeneously distributed grid system is used. It is therefore, concluded that the most homogeneous grid system is the best choice from the viewpoint of numerical accuracy as well as computational efficiency.

[1]  David L. Williamson,et al.  INTEGRATION OF THE PRIMITIVE BAROTROPIC MODEL OVER A SPHERICAL GEODESIC GRID , 1970 .

[2]  Yoshinobu Masuda,et al.  An Integration Scheme of the Primitive Equation Model with an Icosahedral-Hexagonal Grid System and its Application to the Shallow Water Equations , 1986 .

[3]  M. J. P. Cullen,et al.  Integrations of the primitive equations on a sphere using the finite element method , 1974 .

[4]  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 .

[5]  Application of double Fourier series to the shallow-water equations on a sphere , 2000 .

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

[7]  G. R. Stuhne,et al.  New icosahedral grid-point discretization of the shallow water equation on the sphere , 1999 .

[8]  Hyeong-Bin Cheong,et al.  Double Fourier Series on a Sphere , 2000 .

[9]  W. Skamarock,et al.  The stability of time-split numerical methods for the hydrostatic and the nonhydrostatic elastic equations , 1992 .

[10]  G. R. Stuhne,et al.  Vortex Erosion and Amalgamation in a New Model of Large Scale Flow on the Sphere , 1996 .

[11]  Michael Buchhold,et al.  The Operational Global Icosahedral-Hexagonal Gridpoint Model GME: Description and High-Resolution Tests , 2002 .

[12]  J. Hack,et al.  Spectral transform solutions to the shallow water test set , 1995 .

[13]  Hirofumi Tomita,et al.  Shallow water model on a modified icosahedral geodesic grid by using spring dynamics , 2001 .

[14]  Todd D. Ringler,et al.  Modeling the Atmospheric General Circulation Using a Spherical Geodesic Grid: A New Class of Dynamical Cores , 2000 .

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

[16]  M. J. P. Cullen,et al.  Forecasting and general circulation results from finite element models , 1979 .

[17]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

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

[19]  J. Thuburn.  A PV-Based Shallow-Water Model on a Hexagonal-Icosahedral Grid , 1997 .

[20]  Michael Selwyn Longuet-Higgins,et al.  The eigenfunctions of Laplace's tidal equation over a sphere , 1968, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[21]  P. Swarztrauber,et al.  The Vector Harmonic Analysis of Laplace's Tidal Equations , 1985 .