A mimetic, semi-implicit, forward-in-time, finite volume shallow water model: comparison of hexagonal-icosahedral and cubed-sphere grids

Abstract. A new algorithm is presented for the solution of the shallow water equations on quasi-uniform spherical grids. It combines a mimetic finite volume spatial discretization with a Crank–Nicolson time discretization of fast waves and an accurate and conservative forward-in-time advection scheme for mass and potential vorticity (PV). The algorithm is implemented and tested on two families of grids: hexagonal–icosahedral Voronoi grids, and modified equiangular cubed-sphere grids. Results of a variety of tests are presented, including convergence of the discrete scalar Laplacian and Coriolis operators, advection, solid body rotation, flow over an isolated mountain, and a barotropically unstable jet. The results confirm a number of desirable properties for which the scheme was designed: exact mass conservation, very good available energy and potential enstrophy conservation, consistent mass, PV and tracer transport, and good preservation of balance including vanishing ∇ × ∇, steady geostrophic modes, and accurate PV advection. The scheme is stable for large wave Courant numbers and advective Courant numbers up to about 1. In the most idealized tests the overall accuracy of the scheme appears to be limited by the accuracy of the Coriolis and other mimetic spatial operators, particularly on the cubed-sphere grid. On the hexagonal grid there is no evidence for damaging effects of computational Rossby modes, despite attempts to force them explicitly.

[1]  Rupert Ford,et al.  Balance and the Slow Quasimanifold: Some Explicit Results , 2000 .

[2]  M. Cullen.  On the accuracy of the semi‐geostrophic approximation , 2000 .

[3]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

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

[5]  Feng Xiao,et al.  A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid , 2010, J. Comput. Phys..

[6]  W. P. Crowley,et al.  NUMERICAL ADVECTION EXPERIMENTS1 , 1968 .

[7]  Yong Li,et al.  Numerical simulations of Rossby–Haurwitz waves , 2000 .

[8]  M. J. P. Cullen,et al.  A comparison of numerical solutions to the Eady frontogenesis problem , 2008 .

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

[10]  Christiane Jablonowski,et al.  Some considerations for high-order 'incremental remap'-based transport schemes: Edges, reconstructions, and area integration , 2013 .

[11]  R. K. Scott,et al.  An initial-value problem for testing numerical models of the global shallow-water equations , 2004 .

[12]  B. P. Leonard,et al.  The Flux-integral Method for Multidimensional Convection and Diffusion , 1994 .

[13]  James Kent,et al.  Cascades, backscatter and conservation in numerical models of two‐dimensional turbulence , 2014 .

[14]  J. Thuburn,et al.  Numerical wave propagation on the hexagonal C-grid , 2008, J. Comput. Phys..

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

[16]  M. J. P. Cullen,et al.  Modelling atmospheric flows , 2007, Acta Numerica.

[17]  Hiroaki Miura,et al.  An Upwind-Biased Conservative Transport Scheme for Multistage Temporal Integrations on Spherical Icosahedral Grids , 2013 .

[18]  David A. Randall,et al.  Optimized Icosahedral Grids: Performance of Finite-Difference Operators and Multigrid Solver , 2013 .

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

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

[21]  Wayne H. Schubert,et al.  Multigrid methods for elliptic problems: a review , 1986 .

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

[23]  J. Thuburn.  Multidimensional Flux-Limited Advection Schemes , 1996 .

[24]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[25]  R. Pielke,et al.  The forward-in-time upstream advection scheme:extension to higher orders , 1987 .

[26]  Qiang Du,et al.  Centroidal Voronoi Tessellations: Applications and Algorithms , 1999, SIAM Rev..

[27]  Hilary Weller,et al.  Controlling the Computational Modes of the Arbitrarily Structured C Grid , 2012 .

[28]  James Kent,et al.  Assessing implicit large eddy simulation for two-dimensional flow , 2012 .

[29]  Nigel Wood,et al.  An inherently mass‐conserving semi‐implicit semi‐Lagrangian discretisation of the shallow‐water equations on the sphere , 2009 .

[30]  William H. Lipscomb,et al.  An Incremental Remapping Transport Scheme on a Spherical Geodesic Grid , 2005 .

[31]  Colin J. Cotter,et al.  Mixed finite elements for numerical weather prediction , 2011, J. Comput. Phys..

[32]  William Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Reconstructions for Forward-in-Time Schemes , 2010 .

[33]  F. Mesinger,et al.  A global shallow‐water model using an expanded spherical cube: Gnomonic versus conformal coordinates , 1996 .

[34]  G. R. Stuhne,et al.  A higher order discontinuous Galerkin, global shallow water model: Global ocean tides and aquaplanet benchmarks , 2013 .

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

[36]  Bin Wang,et al.  Consistency problem with tracer advection in the Atmospheric Model GAMIL , 2008 .

[37]  P. Paolucci,et al.  The “Cubed Sphere” , 1996 .

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

[39]  B. P. Leonard,et al.  Positivity-preserving numerical schemes for multidimensional advection , 1993 .

[40]  Nigel Wood,et al.  On the relation between order of accuracy, convergence rate and spectral slope for linear numerical methods applied to multiscale problems , 2008 .

[41]  Shian-Jiann Lin,et al.  An explicit flux‐form semi‐lagrangian shallow‐water model on the sphere , 1997 .

[42]  Colin J. Cotter,et al.  A Framework for Mimetic Discretization of the Rotating Shallow-Water Equations on Arbitrary Polygonal Grids , 2012, SIAM J. Sci. Comput..

[43]  Robert Scheichl,et al.  Massively parallel solvers for elliptic partial differential equations in numerical weather and climate prediction , 2013, ArXiv.

[44]  John Thuburn,et al.  Horizontal grids for global weather and climate prediction models: a review , 2012 .