Towards an efficient and scalable discontinuous Galerkin atmospheric model

An efficient and scalable discontinuous Galerkin shallow water model on the cubed sphere is developed by extending the transport scheme of Nair et al. (2005). The continuous flux form nonlinear shallow water equations in curvilinear coordinates are developed. Spatial discretization is a nodal basis set of Legendre polynomials. Fluxes along internal element interfaces are approximated by a Lax-Friedrichs scheme. A third-order total variation diminishing Runge-Kutta scheme is applied for time integration, without any filter or limiter. The standard shallow-water test suite of Williamson et al. (1992) is used to validate the model. It is observed that the numerical solutions are accurate, the model conserves mass to machine precision, and there are no spurious oscillations in a test case where zonal flow impinges a mountain. Development time was substantially reduced by building the model in the high order method modeling environment (HOMME) developed at the National Center for Atmospheric Research (NCAR). Performance and scaling data for the steady state geostrophic flow problem is presented. Sustained performance in excess of 10% of peak is observed out to 64 processors on a Linux cluster.

[1]  Paul F. Fischer,et al.  Terascale Spectral Element Algorithms and Implementations , 1999, ACM/IEEE SC 1999 Conference (SC'99).

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

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

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

[5]  Stephen J. Thomas,et al.  Semi-Implicit Spectral Element Atmospheric Model , 2002, J. Sci. Comput..

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

[7]  R.D. Loft,et al.  Terascale Spectral Element Dynamical Core for Atmospheric General Circulation Models , 2001, ACM/IEEE SC 2001 Conference (SC'01).

[8]  John M. Dennis,et al.  Partitioning with space-filling curves on the cubed-sphere , 2003, Proceedings International Parallel and Distributed Processing Symposium.

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

[10]  H.M. Tufo,et al.  Terascale Spectral Element Algorithms and Implementations , 1999, ACM/IEEE SC 1999 Conference (SC'99).

[11]  Mark A. Taylor,et al.  The Spectral Element Atmosphere Model (SEAM): High-Resolution Parallel Computation and Localized Resolution of Regional Dynamics , 2004 .

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

[13]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

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

[15]  Christopher Edwards,et al.  Multi-scale geophysical modeling using the spectral element method , 2002, Comput. Sci. Eng..

[16]  Francis X. Giraldo,et al.  A Scalable Spectral Element Eulerian Atmospheric Model (SEE-AM) for NWP: Dynamical Core Tests , 2004 .

[17]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[18]  Stephen J. Thomas,et al.  A Discontinuous Galerkin Transport Scheme on the Cubed Sphere , 2005 .

[19]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[20]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

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

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