Improved smoothness and homogeneity of icosahedral grids using the spring dynamics method

An icosahedral grid that has both high smoothness and homogeneity is proposed. The grid-generation method is based on the combination of the spring dynamics (SPR) method with zero natural spring length (SPR0) and transformation by a smooth analytic function around the 12 vertices of an icosahedron. As a preliminary step, we first showed that the grid interval of the grid generated by SPR0 was inversely proportional to a Lambert conformal conic projection map factor, with a map angle of 300^o around the vertices. Then, the transformation function was analytically determined, such that the resolution for the azimuthal direction became constant. In order to estimate cost-efficiency of numerical simulation with the newly proposed grid, we introduced an index defined as the ratio between the minimum grid interval and the squared maximum grid interval. It showed a 2.5% improvement from a recursive grid, and a 0.3-12% improvement from the best cases of the original SPR grid proposed by Tomita et al. (2002) [23] [hereafter, T02] dependent on global resolution. We also re-examined the original SPR method and found that the natural spring length proposed in T02 should be shortened to avoid instability when the global resolution is higher than grid-level 8. Finally, we examined the grids using advection/shallow water simulations.

[1]  William C. Skamarock,et al.  A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids , 2010, J. Comput. Phys..

[2]  P. Frederickson,et al.  Icosahedral Discretization of the Two-Sphere , 1985 .

[3]  Hirofumi Tomita,et al.  An optimization of the Icosahedral grid modified by spring dynamics , 2002 .

[4]  T. Ringler,et al.  Analysis of Discrete Shallow-Water Models on Geodesic Delaunay Grids with C-Type Staggering , 2005 .

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

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

[7]  Ots Price.  INTEGRAL AND SPHERICAL-HARMONIC ANALYSES OF THE GEOMAGNETIC FIELD FOR 1955.0 , 1963 .

[8]  Hirofumi Tomita,et al.  A new dynamical framework of nonhydrostatic global model using the icosahedral grid , 2004 .

[9]  Todd D. Ringler,et al.  A Multiscale Nonhydrostatic Atmospheric Model Using Centroidal Voronoi Tesselations and C-Grid Staggering , 2012 .

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

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

[12]  W. L. Sibley,et al.  Integral and spherical-harmonic analyses of the geomagnetic field for 1955.0, part 2 , 1963 .

[13]  Masaki Satoh,et al.  Nonhydrostatic icosahedral atmospheric model (NICAM) for global cloud resolving simulations , 2008, J. Comput. Phys..

[14]  Janusz A. Pudykiewicz.  On numerical solution of the shallow water equations with chemical reactions on icosahedral geodesic grid , 2011, J. Comput. Phys..

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

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

[17]  Piotr K. Smolarkiewicz,et al.  Predicting weather, climate and extreme events , 2008, J. Comput. Phys..

[18]  C. Jablonowski,et al.  Moving Vortices on the Sphere: A Test Case for Horizontal Advection Problems , 2008 .

[19]  Charles A. Doswell,et al.  A Kinematic Analysis of Frontogenesis Associated with a Nondivergent Vortex , 1984 .

[20]  J. Steppeler,et al.  Third-Order Finite-Difference Schemes on Icosahedral-Type Grids on the Sphere , 2008 .

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

[22]  William C. Skamarock,et al.  Numerical representation of geostrophic modes on arbitrarily structured C-grids , 2009, J. Comput. Phys..

[23]  Ning Wang,et al.  Geometric Properties of the Icosahedral-Hexagonal Grid on the Two-Sphere , 2011, SIAM J. Sci. Comput..

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

[25]  Alexander E. MacDonald,et al.  A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate , 2009 .

[26]  H. Miura.  An Upwind-Biased Conservative Advection Scheme for Spherical Hexagonal–Pentagonal Grids , 2007 .

[27]  Jean Côté,et al.  Cascade interpolation for semi‐Lagrangian advection over the sphere , 1999 .

[28]  Guoliang Xu,et al.  Discrete Laplace-Beltrami Operator on Sphere and Optimal Spherical Triangulations , 2006, Int. J. Comput. Geom. Appl..

[29]  Almut Gassmann,et al.  Towards a consistent numerical compressible non‐hydrostatic model using generalized Hamiltonian tools , 2008 .

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

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