A global finite-element shallow-water model supporting continuous and discontinuous elements

Abstract. This paper presents a novel nodal finite-element method for either continuous and discontinuous elements, as applied to the 2-D shallow-water equations on the cubed sphere. The cornerstone of this method is the construction of a robust derivative operator that can be applied to compute discrete derivatives even over a discontinuous function space. A key advantage of the robust derivative is that it can be applied to partial differential equations in either a conservative or a non-conservative form. However, it is also shown that discontinuous penalization is required to recover the correct order of accuracy for discontinuous elements. Two versions with discontinuous elements are examined, using either the g1 and g2 flux correction function for distribution of boundary fluxes and penalty across nodal points. Scalar and vector hyperviscosity (HV) operators valid for both continuous and discontinuous elements are also derived for stabilization and removal of grid-scale noise. This method is validated using four standard shallow-water test cases, including geostrophically balanced flow, a mountain-induced Rossby wave train, the Rossby–Haurwitz wave and a barotropic instability. The results show that although the discontinuous basis requires a smaller time step size than that required for continuous elements, the method exhibits better stability and accuracy properties in the absence of hyperviscosity.

[1]  M. Taylor.  The Spectral Element Method for the Shallow Water Equations on the Sphere , 1997 .

[2]  Mark A. Taylor,et al.  CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model , 2012, Int. J. High Perform. Comput. Appl..

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

[4]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

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

[6]  Saulo R. M. Barros,et al.  Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver , 1990 .

[7]  Bram van Leer,et al.  High-order finite-volume methods for the shallow-water equations on the sphere , 2010, J. Comput. Phys..

[8]  Feng Xiao,et al.  Global shallow water models based on multi-moment constrained finite volume method and three quasi-uniform spherical grids , 2014, J. Comput. Phys..

[9]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[10]  Kevin Hamilton,et al.  Explicit global simulation of the mesoscale spectrum of atmospheric motions , 2006 .

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

[12]  Antony Jameson,et al.  A New Class of High-Order Energy Stable Flux Reconstruction Schemes , 2011, J. Sci. Comput..

[13]  S. J. Thomas,et al.  The NCAR spectral element climate dynamical core: Semi-implicit eulerian formulation , 2005 .

[14]  Francis X. Giraldo,et al.  A discontinuous Galerkin method for the shallow water equations in spherical triangular coordinates , 2008, J. Comput. Phys..

[15]  John Thuburn,et al.  Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes , 2005 .

[16]  Chi-Wang Shu,et al.  High Order Strong Stability Preserving Time Discretizations , 2009, J. Sci. Comput..

[17]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[18]  Jan S. Hesthaven,et al.  Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations , 2002 .

[19]  Harold Ritchie,et al.  Application of the Semi-Lagrangian Method to a Spectral Model of the Shallow Water Equations , 1988 .

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

[21]  Wanai Li,et al.  The Discontinuous Galerkin Method , 2014 .

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

[23]  J. Boyd,et al.  The Erf1⁄2-Log Filter and the Asymptotics of the Euler and Vandeven Sequence Accelerations , 1999 .

[24]  James A. Rossmanith,et al.  A wave propagation method for hyperbolic systems on the sphere , 2006, J. Comput. Phys..

[25]  M. A. Tolstykh,et al.  Vorticity-Divergence Semi-Lagrangian Shallow-Water Model of the Sphere Based on Compact Finite Differences , 2002 .

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

[27]  Feng Xiao,et al.  Shallow water model on cubed-sphere by multi-moment finite volume method , 2008, J. Comput. Phys..

[28]  John Thuburn,et al.  Some conservation issues for the dynamical cores of NWP and climate models , 2008, J. Comput. Phys..

[29]  Jean Côté,et al.  Experiments with different discretizations for the shallow‐water equations on a sphere , 2012 .

[30]  Christiane Jablonowski,et al.  MCore: A non-hydrostatic atmospheric dynamical core utilizing high-order finite-volume methods , 2012, J. Comput. Phys..

[31]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Global Shallow Water Model , 2005, Monthly Weather Review.

[32]  Eric Deleersnijder,et al.  A finite element method for solving the shallow water equations on the sphere , 2009 .

[33]  Paul A. Ullrich,et al.  Understanding the treatment of waves in atmospheric models. Part 1: The shortest resolved waves of the 1D linearized shallow‐water equations , 2014 .

[34]  Mikhail A. Tolstykh,et al.  Vorticity-divergence mass-conserving semi-Lagrangian shallow-water model using the reduced grid on the sphere , 2012, J. Comput. Phys..

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

[36]  Todd D. Ringler,et al.  A multiresolution method for climate system modeling: application of spherical centroidal Voronoi tessellations , 2008 .

[37]  F. Xiao,et al.  A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid , 2008 .

[38]  John M. Dennis,et al.  A Comparison of Two Shallow-Water Models with Nonconforming Adaptive Grids , 2008 .

[39]  Antony Jameson,et al.  A New Class of High-Order Energy Stable Flux Reconstruction Schemes for Triangular Elements , 2012, J. Sci. Comput..

[40]  Jean Côté,et al.  An Accurate and Efficient Finite-Element Global Model of the Shallow-Water Equations , 1990 .

[41]  R. Sadourny.  Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

[42]  Todd D. Ringler,et al.  Exploring a Multiresolution Modeling Approach within the Shallow-Water Equations , 2011 .

[43]  Mark A. Taylor,et al.  A compatible and conservative spectral element method on unstructured grids , 2010, J. Comput. Phys..

[44]  Steven J. Ruuth.  Global optimization of explicit strong-stability-preserving Runge-Kutta methods , 2005, Math. Comput..

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

[46]  Lei Bao,et al.  A mass and momentum flux-form high-order discontinuous Galerkin shallow water model on the cubed-sphere , 2014, J. Comput. Phys..