×

Accuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative. (English) Zbl 1349.76376

Summary: Many newly developed climate, weather and ocean global models are based on quasi-uniform spherical polygonal grids, aiming for high resolution and better scalability. J. Thuburn et al. [ibid. 228, No. 22, 8321–8335 (2009; Zbl 1173.86304)] and T. D. Ringler et al. [ibid. 229, No. 9, 3065–3090 (2010; Zbl 1307.76054)] developed a C staggered finite volume/difference method for arbitrary polygonal spherical grids suitable for these next generation dynamical cores. This method has many desirable mimetic properties and became popular, being adopted in some recent models, in spite of being known to possess low order of accuracy. In this work, we show that, for the nonlinear shallow water equations on non-uniform grids, the method has potentially 3 main sources of inconsistencies (local truncation errors not converging to zero as the grid is refined): (i) the divergence term of the continuity equation, (ii) the perpendicular velocity and (iii) the kinetic energy terms of the vector invariant form of the momentum equations. Although some of these inconsistencies have not impacted the convergence on some standard shallow water test cases up until now, they may constitute a potential problem for high resolution 3D models. Based on our analysis, we propose modifications for the method that will make it first order accurate in the maximum norm. It preserves many of the mimetic properties, albeit having non-steady geostrophic modes on the f-sphere. Experimental results show that the resulting model is a more accurate alternative to the existing formulations and should provide means of having a consistent, computationally cheap and scalable atmospheric or ocean model on C staggered Voronoi grids.

MSC:

76M12 Finite volume methods applied to problems in fluid mechanics
86-08 Computational methods for problems pertaining to geophysics
65M08 Finite volume methods for initial value and initial-boundary value problems involving PDEs
76B15 Water waves, gravity waves; dispersion and scattering, nonlinear interaction
86A05 Hydrology, hydrography, oceanography
86A10 Meteorology and atmospheric physics
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Staniforth, A.; Thuburn, J., Horizontal grids for global weather and climate prediction models: a review, Q. J. R. Meteorol. Soc., 138, 662, 1-26 (2012)
[2] Thuburn, J.; Ringler, T. D.; Skamarock, W. C.; Klemp, J. B., Numerical representation of geostrophic modes on arbitrarily structured C-grids, J. Comput. Phys., 228, 8321-8335 (2009) · Zbl 1173.86304
[3] Ringler, T. D.; Thuburn, J.; Klemp, J. B.; Skamarock, W. C., A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids, J. Comput. Phys., 229, 9, 3065-3090 (2010) · Zbl 1307.76054
[4] Thuburn, J.; Cotter, C., A framework for mimetic discretization of the rotating shallow-water equations on arbitrary polygonal grids, SIAM J. Sci. Comput., 34, B203-B225 (2012) · Zbl 1246.65155
[5] Weller, H., Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid, Geosci. Model Dev., 7, 3, 779-797 (2014)
[6] Skamarock, W. C.; Klemp, J. B.; Duda, M. G.; Fowler, L. D.; Park, S.-H.; Ringler, T. D., A multiscale nonhydrostatic atmospheric model using centroidal Voronoi tesselations and C-grid staggering, Mon. Weather Rev., 140, 3090-3105 (2012)
[7] Ringler, T.; Petersen, M.; Higdon, R. L.; Jacobsen, D.; Jones, P. W.; Maltrud, M., A multi-resolution approach to global ocean modeling, Ocean Model., 69, 0, 211-232 (2013)
[8] Dubos, T.; Dubey, S.; Tort, M.; Mittal, R.; Meurdesoif, Y.; Hourdin, F., DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility, Geosci. Model Dev., 8, 10, 3131-3150 (2015)
[9] Thuburn, J.; Cotter, C. J.; Dubos, T., A mimetic, semi-implicit, forward-in-time, finite volume shallow water model: comparison of hexagonal-icosahedral and cubed-sphere grids, Geosci. Model Dev., 7, 3, 909-929 (2014)
[10] Weller, H., Controlling the computational modes of the arbitrarily structured C grid, Mon. Weather Rev., 140, 10, 3220-3234 (2012)
[11] Weller, H.; Thuburn, J.; Cotter, C. J., Computational modes and grid imprinting on five quasi-uniform spherical C-grids, Mon. Weather Rev., 140, 8, 2734-2755 (2012)
[12] Dubos, T.; Kevlahan, N. K.-R., A conservative adaptive wavelet method for the shallow-water equations on staggered grids, Q. J. R. Meteorol. Soc., 139, 677, 1997-2020 (2013)
[13] Arakawa, A.; Lamb, V., Computational design of the basic dynamical processes of the UCLA general circulation model, (Methods in Computational Physics, vol. 17 (1977)), 173-265
[14] Cockburn, B.; Gremaud, P.-A., A priori error estimates for numerical methods for scalar conservation laws. Part II: flux-splitting monotone schemes on irregular Cartesian grids, Math. Comput., 66, 218, 547-572 (1997) · Zbl 0866.65061
[15] Peixoto, P. S.; Barros, S. R.M., Analysis of grid imprinting on geodesic spherical icosahedral grids, J. Comput. Phys., 237, 61-78 (2013) · Zbl 1286.65118
[16] Miura, H.; Kimoto, M., A comparison of grid quality of optimized spherical hexagonal-pentagonal geodesic grids, Mon. Weather Rev., 133, 10, 2817-2833 (2005)
[17] Heikes, R.; Randall, D. A., Numerical integration of the shallow-water equations on a twisted icosahedral grid. Part I: basic design and results of tests, Mon. Weather Rev., 123, 6, 1862-1880 (1995)
[18] Peixoto, P. S.; Barros, S. R., On vector field reconstructions for semi-Lagrangian transport methods on geodesic staggered grids, J. Comput. Phys., 273, 0, 185-211 (2014) · Zbl 1351.86005
[19] Du, Q.; Gunzburger, M. D.; Ju, L., Constrained centroidal Voronoi tessellations for surfaces, SIAM J. Sci. Comput., 24, 1488-1506 (2003) · Zbl 1036.65101
[20] Ju, L.; Ringler, T. D.; Gunzburger, M., Voronoi tessellations and their application to climate and global modeling, (Lauritzen, P.; Jablonowski, C.; Taylor, M.; Nair, R., Numerical Techniques for Global Atmospheric Models. Numerical Techniques for Global Atmospheric Models, Lecture Notes in Computational Science and Engineering, vol. 80 (2011), Springer: Springer Berlin, Heidelberg), 313-342
[21] Tomita, H.; Satoh, M.; Goto, K., An optimization of the icosahedral grid modified by spring dynamics, J. Comput. Phys., 183, 1, 307-331 (2002) · Zbl 1058.76591
[22] Williamson, D. L.; Drake, J. B.; Hack, J. J.; Jakob, R.; Swarztrauber, P. N., A standard test set for numerical approximations to the shallow water equations in spherical geometry, J. Comput. Phys., 102, 1, 211-224 (1992) · Zbl 0756.76060
[23] Gassmann, A., A global hexagonal C-grid non-hydrostatic dynamical core (ICON-IAP) designed for energetic consistency, Q. J. R. Meteorol. Soc., 139, 670, 152-175 (2013)
[24] Perot, B., Conservation properties of unstructured staggered mesh schemes, J. Comput. Phys., 159, 1, 58-89 (2000) · Zbl 0972.76068
[25] Bonaventura, L.; Ringler, T. D., Analysis of discrete shallow-water models on geodesic Delaunay grids with C-type staggering, Mon. Weather Rev., 133, 8, 2351-2373 (2005)
[26] Ham, D. A.; Kramer, S. C.; Stelling, G. S.; Pietrzak, J., The symmetry and stability of unstructured mesh C-grid shallow water models under the influence of Coriolis, Ocean Model., 16, 1-2, 47-60 (2007)
[27] Kleptsova, O.; Pietrzak, J.; Stelling, G., On the accurate and stable reconstruction of tangential velocities in C-grid ocean models, The Sixth International Workshop on Unstructured Mesh Numerical Modelling of Coastal, Shelf and Ocean Flows. The Sixth International Workshop on Unstructured Mesh Numerical Modelling of Coastal, Shelf and Ocean Flows, Ocean Model., 28, 1-3, 118-126 (2009)
[28] Nickovic, S.; Gavrilov, M. B.; Tosic, I. A., Geostrophic adjustment on hexagonal grids, Mon. Weather Rev., 130, 668-683 (2002)
[29] Galewsky, J.; Scott, R. K.; Polvani, L. M., An initial-value problem for testing numerical models of the global shallow-water equations, Tellus A, 56, 5, 429-440 (2004)
[30] Wan, H.; Giorgetta, M. A.; Zängl, G.; Restelli, M.; Majewski, D.; Bonaventura, L.; Fröhlich, K.; Reinert, D.; Rípodas, P.; Kornblueh, L.; Förstner, J., The ICON-1.2 hydrostatic atmospheric dynamical core on triangular grids. Part 1: formulation and performance of the baseline version, Geosci. Model Dev., 6, 3, 735-763 (2013)
[31] Gassmann, A., Inspection of hexagonal and triangular C-grid discretizations of the shallow water equations, J. Comput. Phys., 230, 7, 2706-2721 (2011) · Zbl 1316.76069
[32] Bell, M. J.; Peixoto, P. S.; Thuburn, J., Numerical instabilities of vector invariant momentum equations on rectangular C-grids, [Accessed: 15 Jan 2016]
[33] Thuburn, J.; Zerroukat, M.; Wood, N.; Staniforth, A., Coupling a mass-conserving semi-Lagrangian scheme (SLICE) to a semi-implicit discretization of the shallow-water equations: minimizing the dependence on a reference atmosphere, Q. J. R. Meteorol. Soc., 136, 646, 146-154 (2010)
[34] Ringler, T. D.; Jacobsen, D.; Gunzburger, M.; Ju, L.; Duda, M.; Skamarock, W., Exploring a multiresolution modeling approach within the shallow-water equations, Mon. Weather Rev., 139, 11, 3348-3368 (2011)
[35] Thuburn, J., Numerical wave propagation on the hexagonal C-grid, J. Comput. Phys., 227, 11, 5836-5858 (2008) · Zbl 1220.76018
[36] Thuburn, J.; Li, Y., Numerical simulations of Rossby-Haurwitz waves, Tellus A, 52, 2, 181-189 (2000)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.