‘The Equator–Pole grid system’: an overset grid system for the sphere, with an optimal uniformity property

Göran Starius

doi:10.1080/16000870.2018.1541373

‘The Equator–Pole grid system’: an overset grid system for the sphere, with an optimal uniformity property

Tellus

Volume 70 (2018): Issue 1

By: Göran Starius

Open Access

|Jan 2018

Abstract

This article proposes a new grid system for the sphere, which consists of three orthogonal and almost uniform grids. The basic one is a latitude-longitude grid covering an annular band around the equator. For the rest of the sphere this grid is complemented by two grids covering the polar regions, and based on suitably modified stereographic coordinates. The rectangular structure of the grids makes highly efficient implementation on massively parallel computer systems possible. Numerical experiments with the new grid system are carried out for two advection examples, namely smooth deformational flow and rotation of the Cosine bell, and for the test problems 2, 3, and 6 from Williamson et al., concerning the non-linear shallow water equations. For problem 6, the Rossby–Haurwitz wave, we study conservation properties for mass and energy. The computational results compare favourably with results for other grids. Our focus is the definition of the grid system together with the connection between grids by overlapping and the use of centred interpolation formulas. The method of centred finite differences is used for the spatial discretisation of the differential equations, because it is well-known and easy to implement.

AMS Subject Classification: 35L45; 65M06

References

Bao, L., Nair, R. D. and Tufo, H. M. 2014. A mass and momentum flux-form higher-order discontinuous Galerkin shallow water model on the cubed-sphere. J. Comput. Phys. 271, 224–243. doi:10.1016/j.jcp.2013.11.033.
Open DOI Search in Google Scholar Back to article
Browning, B. L., Hack, J. J. and Swarztrauber, P. N. 1989. A comparison of three numerical methods for solving differential equations on the sphere. Mon. Wea. Rev. 117, 1058–1075. doi:10.1175/1520-0493(1989)117<;1058:ACOTNM>2.0.CO;2.
Open DOI Search in Google Scholar Back to article
Chen, C., Li, X., Shen, X. and Xiao, F. 2014. Global shallow water models based on multi-moment constrained finite volume method and three quasi-uniform spherical grids. J. Comput. Phys. 271, 191–223. doi:10.1016/j.jcp.2013.10.026.
Open DOI Search in Google Scholar Back to article
Clancy, C. and Pudykiewics, J. A. 2013. On the use of exponential time integration methods in atmospheric models. Tellus A 65, 20898. doi:10.3402/tellusa.v65i0.20898.
Open DOI Search in Google Scholar Back to article
Gates, W. L. and Riegel, A. 1963. A study of numerical errors in the integration of barotropic flow on a spherical grid. JGR 67, 406–423.
Search in Google Scholar Back to article
Gaudreault, S. and Pudykiewicz, J. A. 2016. An efficient exponential time integration method for the numerical solution if the shallow water equations on the sphere. J. Comput. Phys. 322, 827–848. doi:10.1016/j.jcp.2016.07.012.
Open DOI Search in Google Scholar Back to article
Giraldo, F. X. 2006. High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere. J. Comput. Phys. 214, 447–465. doi:10.1016/j.jcp.2005.09.029.
Open DOI Search in Google Scholar Back to article
Henshaw, W. D. and Schwendeman, D. W. 2006. Moving overlapping grids with adaptive mesh refinement for high-speed reactive and non-reactive flow. J. Comput. Phys. 216, 744–779. doi:10.1016/j.jcp.2006.01.005.
Open DOI Search in Google Scholar Back to article
Ii, S. and Xiao, F. 2010. A global shallow water model using high order constrained finite volume method and icosahedral grid. J. Comput. Phys. 229, 1774–1796. doi:10.1016/j.jcp.2009.11.008.
Open DOI Search in Google Scholar Back to article
Kageyama, A. and Sato, T. 2004. The Yin-Yang grid: An overset grid in spherical geometry. Geochem. Geophys. Geosyst. 5, Q09005.
Search in Google Scholar Back to article
Kurihara, Y. 1965. Integration of the primitive equations on a spherical grid. Mon. Wea. Rev. 93, 399–415. doi:10.1175/1520-0493(1965)093<;0399:NIOTPE>2.3.CO;2.
Open DOI Search in Google Scholar Back to article
Lanser, D., Blom, J. G. and Verwer, J. G. 2000. Spatial discretisation of the shallow water equations in spherical geometry using Osher’s scheme. J. Comput. Phys. 165, 542–565. doi:10.1006/jcph.2000.6632.
Open DOI Search in Google Scholar Back to article
Lauritzen, P. H., Nair, R. D. and Ullrich, P. A. 2010. A conservative semi-Lagrangian multi-tracer scheme (CSLAM) on a cubed-sphere grid. J. Comput. Phys. 229, 1401–1424. doi:10.1016/j.jcp.2009.10.036.
Open DOI Search in Google Scholar Back to article
Läuter, M., Giraldo, F. X., Handorf, D. and Dethloff, K. 2008. A discontinuous Galerkin method for the shallow water equations using spherical triangular coordinates. J. Comput. Phys. 227, 10226–10243. doi:10.1016/j.jcp.2008.08.019.
Open DOI Search in Google Scholar Back to article
Nair, R. D. and Jablonowski, C. 2008. Moving vortices on the sphere: A test case for horizontal advection problems. Mon. Wea. Rev. 136, 699–711. doi:10.1175/2007MWR2105.1.
Open DOI Search in Google Scholar Back to article
Nair, R. D., Thomas, S. J. and Loft, R. D. 2005. A discontinuous Galerkin transport scheme on the cubed sphere. Mon. Wea. Rev. 133, 814–828. doi:10.1175/MWR2890.1.
Open DOI Search in Google Scholar Back to article
Peixoto, P. S. and Barros, S. R. M. 2013. Analysis of grid imprinting on geodesic spherical icosahedral grids. J. Comput. Phys. 237, 61–78. doi:10.1016/j.jcp.2012.11.041.
Open DOI Search in Google Scholar Back to article
Phillips, N. A. 1959. Numerical integration of the primitive equations on the hemisphere. Mon. Wea. Rev. 87, 333–345. doi:10.1175/1520-0493(1959)087<;0333:NIOTPE>2.0.CO;2.
Open DOI Search in Google Scholar Back to article
Pudykiewics, J. A. 2011. On numerical solution of the shallow water equations with chemical reactions on icoshedral geodesic grid. J. Comput. Phys. 230, 1956–1991. doi:10.1016/j.jcp.2010.11.045.
Open DOI Search in Google Scholar Back to article
Putman, W. M. and Lin, S.-J. 2007. Finite-volume transport on various cubed-sphere grids. J. Comput. Phys. 227, 55–78. doi:10.1016/j.jcp.2007.07.022.
Open DOI Search in Google Scholar Back to article
Qaddouri, A. and Lee, V. 2011. The Canadian global environmental multiscale model on the Yin-Yang grid system. Qjr. Meteorol. Soc. 137, 1913–1926. doi:10.1002/qj.873.
Open DOI Search in Google Scholar Back to article
Robert, A. 1982. A semi-Lagrangian and semi-implicit numerical integration scheme for the primitive meteorological equations. J. Meteorol. Soc. Jpn. 60, 319–325. doi:10.2151/jmsj1965.60.1_319.
Open DOI Search in Google Scholar Back to article
Ronchi, C., Iacono, R. and Paolucci, P. S. 1996. The cubed sphere: A new method for the solution of partial differential equations in spherical geometry. J. Comput. Phys. 124, 93–114. doi:10.1006/jcph.1996.0047.
Open DOI Search in Google Scholar Back to article
Sadourny, R. 1972. Conservative finite-differencing of the primitive equations on quasi-uniform spherical grids. Mon. Wea. Rev. 100, 136–144. doi:10.1175/1520-0493(1972)100<;0136:CFAOTP>2.3.CO;2.
Open DOI Search in Google Scholar Back to article
Staniforth, A. and Thuburn, J. 2012. Horizontal grids for global weather and climate prediction models: a review. Qjr. Meteorol. Soc. 138, 1–26. doi:10.1002/qj.958.
Open DOI Search in Google Scholar Back to article
Starius, G. 1977a. Composite mesh difference methods for elliptic boundary value problems. Numer. Math. 28, 243–258. doi:10.1007/BF01394455.
Open DOI Search in Google Scholar Back to article
Starius, G. 1977b. Constructing orthogonal curvilinear meshes by solving initial value problems. Numer. Math. 28, 25–48. doi:10.1007/BF01403855.
Open DOI Search in Google Scholar Back to article
Starius, G. 1980. On composite mesh difference methods for hyperbolic differential equations. Numer. Math. 35, 241–255. doi:10.1007/BF01396411.
Open DOI Search in Google Scholar Back to article
Starius, G. 2014. A solution to the pole problem for the shallow water equations on a sphere. TWMS JPAM. 5, 152–170.
Search in Google Scholar Back to article
Ullrich, P. A., Jablonowski, C. and van Leer, B. 2010. High-order finite-volume methods for the shallow-water equations on the sphere. J. Comput. Phys. 229, 6104–6134. doi:10.1016/j.jcp.2010.04.044.
Open DOI Search in Google Scholar Back to article
Volkov, E. A. 1968. The method of composite meshes for finite and infinite regions with piecewise smooth boundaries. Proc. Steklov Inst. Math. 96, 145–185.
Search in Google Scholar Back to article
Williamson, D. L. 1968. Integration of the barotropic vorticity equation on a spherical geodesic grid. Tellus. 20, 642–653.
Search in Google Scholar Back to article
Williamson, D. L. 2007. The evolution of dynamical cores for global atmospheric model. JMSJ. 85B, 241–269. doi:10.2151/jmsj.85B.241.
Open DOI Search in Google Scholar Back to article
Williamson, D. L., Drake, J. B., Hack, J. J., Jakob, R. and Swarztrauber, P. N. 1992. A standard test set for numerical approximation to the shallow water equations in spherical geometry. J. Comput. Phys. 102, 211–224. doi:10.1016/S0021-9991(05)80016-6.
Open DOI Search in Google Scholar Back to article