A Discontinuous Galerkin Transport Scheme on the Cubed Sphere

Ramachandran D. Nair Scientific Computing Division, National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Ramachandran D. Nair in
Current site
Google Scholar
PubMed
Close
,
Stephen J. Thomas Scientific Computing Division, National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Stephen J. Thomas in
Current site
Google Scholar
PubMed
Close
, and
Richard D. Loft Scientific Computing Division, National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Richard D. Loft in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

A conservative transport scheme based on the discontinuous Galerkin (DG) method has been developed for the cubed sphere. Two different central projection methods, equidistant and equiangular, are employed for mapping between the inscribed cube and the sphere. These mappings divide the spherical surface into six identical subdomains, and the resulting grid is free from singularities. Two standard advection tests, solid-body rotation and deformational flow, were performed to evaluate the DG scheme. Time integration relies on a third-order total variation diminishing (TVD) Runge–Kutta scheme without a limiter. The numerical solutions are accurate and neither exhibit shocks nor discontinuities at cube-face edges and vertices. The numerical results are either comparable or better than a standard spectral element method. In particular, it was found that the standard relative error metrics are significantly smaller for the equiangular as opposed to the equidistant projection.

Corresponding author address: Ramachandran D. Nair, Scientific Computing Division, National Center for Atmospheric Research, 1850 Table Mesa Drive, Boulder, CO 80305. Email: rnair@ucar.edu

Abstract

A conservative transport scheme based on the discontinuous Galerkin (DG) method has been developed for the cubed sphere. Two different central projection methods, equidistant and equiangular, are employed for mapping between the inscribed cube and the sphere. These mappings divide the spherical surface into six identical subdomains, and the resulting grid is free from singularities. Two standard advection tests, solid-body rotation and deformational flow, were performed to evaluate the DG scheme. Time integration relies on a third-order total variation diminishing (TVD) Runge–Kutta scheme without a limiter. The numerical solutions are accurate and neither exhibit shocks nor discontinuities at cube-face edges and vertices. The numerical results are either comparable or better than a standard spectral element method. In particular, it was found that the standard relative error metrics are significantly smaller for the equiangular as opposed to the equidistant projection.

Corresponding author address: Ramachandran D. Nair, Scientific Computing Division, National Center for Atmospheric Research, 1850 Table Mesa Drive, Boulder, CO 80305. Email: rnair@ucar.edu

Save
  • Biswas, R., K. Devine, and J. Flaherty, 1994: Parallel, advaptive finite element methods for conservation laws. Appl. Numer. Maths., 14 , 255283.

    • Search Google Scholar
    • Export Citation
  • Cockburn, B., and C. W. Shu, 1989: TVB Runge–Kutta local projection discontinuous Galerkin method for conservation laws. II: General framework. Math. Comput., 52 , 411435.

    • Search Google Scholar
    • Export Citation
  • Cockburn, B., and C. W. Shu, 1998: The Runge–Kutta discontinuous Galerkin finite element method for conservation laws. V: Multidimensional systems. J. Comput. Phys., 141 , 199224.

    • Search Google Scholar
    • Export Citation
  • Cockburn, B., and C. W. Shu, 2001: The Runge–Kutta discontinuous Galerkin method for convection-dominated problems. J. Sci. Comput., 16 , 173261.

    • Search Google Scholar
    • Export Citation
  • Cockburn, B., G. E. Karniadakis, and C. W. Shu, 2000: Discontinuous Galerkin Methods: Theory, Computation, and Applications. Lecture Notes in Computational Science and Engineering, Vol. 11, Springer-Verlag, 470 pp.

    • Search Google Scholar
    • Export Citation
  • Doswell III, C. A., 1984: A kinematic analysis of frontogenesis associated with a nondivergent vortex. J. Atmos. Sci., 41 , 12421248.

  • Dutton, J. A., 1986: Dynamics of Atmospheric Motion. Dover, 617 pp.

  • Fournier, A., M. A. Taylor, and J. J. Tribbia, 2004: A spectral element atmospheric model (SEAM): High-resolution parallel computation and localized resolution of regional dynamics. Mon. Wea. Rev., 132 , 726748.

    • Search Google Scholar
    • Export Citation
  • Giraldo, F. X., and T. E. Rosmond, 2004: A scalable spectral element Eulerian atmospheric model (SEE-AM) for NWP: Dynamical core test. Mon. Wea. Rev., 132 , 133153.

    • Search Google Scholar
    • Export Citation
  • Giraldo, F. X., J. S. Hesthaven, and T. Warburton, 2002: Nodal high-order discontinuous Galerkin methods for spherical shallow water equations. J. Comput. Phys., 181 , 499525.

    • Search Google Scholar
    • Export Citation
  • Gottlieb, S., C. W. Shu, and E. Tadmor, 2001: Strong stability-preserving high-order time discretization methods. SIAM Rev., 43 , 89112.

    • Search Google Scholar
    • Export Citation
  • Henderson, R. D., 1999: Adaptive spectral element method for turbulence and transition. High-Order Methods for Computational Physics, T. J. Barth and H. Deconinck, Eds., Lecture Notes in Computational Science and Engineering, Vol. 9, Springer-Verlag, 225–324.

    • Search Google Scholar
    • Export Citation
  • Karniadakis, G. E., and S. J. Sherwin, 1999: Spectral/hp Element Methods for CFD. Oxford University Press, 390 pp.

  • Lesaint, P., and P. A. Raviart, 1974: On finite-element methods for solving the neutron transport equations. Mathematical Aspects of Finite Elements in Partial Differential Equations, C. de Boor, Ed., Academic Press, 89–145.

    • Search Google Scholar
    • Export Citation
  • McGregor, J. L., 1996: Semi-Lagrangian advection on a conformal cubic grid. Mon. Wea. Rev., 124 , 13111322.

  • Nair, R. D., and B. Machenhauer, 2002: The mass conservative cell-integrated semi-Lagrangian advection scheme on the sphere. Mon. Wea. Rev., 130 , 649667.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., J. Côté, and A. Staniforth, 1999: Cascade interpolation for semi-Lagrangian advection over the sphere. Quart. J. Roy. Meteor. Soc., 125 , 14451468.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., J. S. Scroggs, and F. H. M. Semazzi, 2003: A forward trajectory semi-Lagrangian scheme over the sphere. J. Comput. Phys., 190 , 275294.

    • Search Google Scholar
    • Export Citation
  • Rančić, M. R., J. Purser, and F. Mesinger, 1996: A global-shallow water model using an expanded spherical cube. Quart. J. Roy. Meteor. Soc., 122 , 959982.

    • Search Google Scholar
    • Export Citation
  • Randall, D. A., T. D. Ringler, R. P. Heikes, P. Jones, and J. Baumgardner, 2002: Climate modeling with spherical geodesic grids. Comput. Sci. Eng., 4 , 3241.

    • Search Google Scholar
    • Export Citation
  • Ronchi, C., R. Iacono, and P. S. Paolucci, 1996: The “cubed sphere”: A new method for the solution of partial differential equations in spherical geometry. J. Comput. Phys., 124 , 93114.

    • Search Google Scholar
    • Export Citation
  • Sadourny, R., 1972: Conservative finite-difference approximations of the primitive equations on quasi-uniform spherical grids. Mon. Wea. Rev., 100 , 136144.

    • Search Google Scholar
    • Export Citation
  • Shu, C-W., 1988: Total variation diminishing time discretizations. SIAM J. Sci. Stat. Comput., 9 , 10731084.

  • Taylor, M., J. Tribbia, and M. Iskandrani, 1997: The spectral element method for the shallow water equations on the sphere. J. Comput. Phys., 130 , 92108.

    • Search Google Scholar
    • Export Citation
  • Thomas, S. J., and R. D. Loft, 2002: Semi-implicit spectral element model. J. Sci. Comput., 17 , 339350.

  • Williamson, D. L., J. B. Drake, J. Hack, R. Jacob, and P. N. Swartztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys., 102 , 211224.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1409 693 55
PDF Downloads 863 147 19