A Nonoscillatory Discontinuous Galerkin Transport Scheme on the Cubed Sphere

Yifan Zhang Division of Applied Mathematics, Brown University, Providence, Rhode Island

Search for other papers by Yifan Zhang in
Current site
Google Scholar
PubMed
Close
and
Ramachandran D. Nair National Center for Atmospheric Research,* Boulder, Colorado

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

Abstract

The discontinuous Galerkin (DG) method is high order, conservative, and offers excellent parallel efficiency. However, when there are discontinuities in the solution, the DG transport scheme generates spurious oscillations that reduce the quality of the numerical solution. For applications such as the atmospheric tracer transport modeling, a nonoscillatory, positivity-preserving solution is a basic requirement. To suppress the oscillations in the DG solution, a limiter based on the Hermite-Weighted Essentially Nonoscillatory (H-WENO) method has been implemented for a third-order DG transport scheme. However, the H-WENO limiter can still produce wiggles with small amplitudes in the solutions, but this issue has been addressed by combining the limiter with a bound-preserving (BP) filter. The BP filter is local and easy to implement and can be used for making the solution strictly positivity preserving. The DG scheme combined with the limiter and filter preserves the accuracy of the numerical solution in the smooth regions while effectively eliminating overshoots and undershoots. The resulting nonoscillatory DG scheme is third-order accurate (P2-DG) and based on the modal discretization. The 2D Cartesian scheme is further extended to the cubed-sphere geometry, which employs nonorthogonal, curvilinear coordinates. The accuracy and effectiveness of the limiter and filter are demonstrated with several benchmark tests on both the Cartesian and spherical geometries.

The National Center for Atmospheric Research is sponsored by the National Science Foundation.

Corresponding author address: R. D. Nair, Computational and Information System Laboratory, National Center for Atmospheric Research, Boulder, CO 80305. E-mail: rnair@ucar.edu

Abstract

The discontinuous Galerkin (DG) method is high order, conservative, and offers excellent parallel efficiency. However, when there are discontinuities in the solution, the DG transport scheme generates spurious oscillations that reduce the quality of the numerical solution. For applications such as the atmospheric tracer transport modeling, a nonoscillatory, positivity-preserving solution is a basic requirement. To suppress the oscillations in the DG solution, a limiter based on the Hermite-Weighted Essentially Nonoscillatory (H-WENO) method has been implemented for a third-order DG transport scheme. However, the H-WENO limiter can still produce wiggles with small amplitudes in the solutions, but this issue has been addressed by combining the limiter with a bound-preserving (BP) filter. The BP filter is local and easy to implement and can be used for making the solution strictly positivity preserving. The DG scheme combined with the limiter and filter preserves the accuracy of the numerical solution in the smooth regions while effectively eliminating overshoots and undershoots. The resulting nonoscillatory DG scheme is third-order accurate (P2-DG) and based on the modal discretization. The 2D Cartesian scheme is further extended to the cubed-sphere geometry, which employs nonorthogonal, curvilinear coordinates. The accuracy and effectiveness of the limiter and filter are demonstrated with several benchmark tests on both the Cartesian and spherical geometries.

The National Center for Atmospheric Research is sponsored by the National Science Foundation.

Corresponding author address: R. D. Nair, Computational and Information System Laboratory, National Center for Atmospheric Research, Boulder, CO 80305. E-mail: rnair@ucar.edu
Save
  • Balsara, D. S., C. Altman, C.-D. Munz, and M. Dumbser, 2007: Sub-cell based indicator for troubled zones in RKDG schemes and a novel class of hybrid RKDG+HWENO schemes. J. Comput. Phys., 226, 586620.

    • Search Google Scholar
    • Export Citation
  • Chen, C., and F. Xiao, 2008: Shallow water model on cubed-sphere by multi-moment finite volume method. J. Comput. Phys., 227, 50195044.

    • Search Google Scholar
    • Export Citation
  • Cockburn, B., 1997: An introduction to the Discontinuous-Galerkin method for convection-dominated problems. Advanced Numerical Approximation of Nonlinear Hyperbolic Equations, A. Quarteroni, Ed., Vol. 1697, Lecture Notes in Mathematics, Springer, 151–268.

  • Cockburn, B., and C.-W. Shu, 1989: TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: General framework. Math. Comput.,52, 411–435.

  • Cockburn, B., S. Hou, and C.-W. Shu, 1990: TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: The multidimensional case. Math. Comput., 54, 545581.

    • Search Google Scholar
    • Export Citation
  • Flyer, N., and G. B. Wright, 2007: Transport schemes on a sphere using radial basis functions. J. Comput. Phys., 226, 10591084.

  • Giraldo, F., J. 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 (1), 89112.

    • Search Google Scholar
    • Export Citation
  • Iskandarani, M., J. Levin, B.-J. Choi, and D. Haidgovel, 2005: Comparison of advection schemes for high-order h-p finite element and finite volume methods. Ocean Modell., 10, 233252.

    • Search Google Scholar
    • Export Citation
  • Krivodonova, L., 2007: Limiters for high-order discontinuous Galerkin methods. J. Comput. Phys., 226, 879896.

  • Läuter, M., F. X. Giraldo, D. Handorf, and K. Dethloff, 2008: A discontinuous Galerkin method for shallow water equations in spherical triangular coordinates. J. Comput. Phys., 227, 10 22610 242.

    • Search Google Scholar
    • Export Citation
  • LeVeque, R. J., 2002: Finite Volume Methods for Hyperbolic Problems. Cambridge University Press, 558 pp.

  • Levy, M. N., R. D. Nair, and H. M. Tufo, 2007: High-order Galerkin method for scalable global atmospheric models. Comput. Geosci., 33, 10221035.

    • Search Google Scholar
    • Export Citation
  • Liu, X.-D., and S. Osher, 1996: Nonoscillatory high order accurate self-similar maximum principle satisfying shock capturing schemes I. SIAM J. Numer. Anal., 33, 760779.

    • Search Google Scholar
    • Export Citation
  • Luo, H., J. D. Baum, and R. Lohner, 2007: A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids. J. Comput. Phys., 225, 686713.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., and C. Jablonowski, 2008: Moving vortices on the sphere: A test case for horizontal advection problems. Mon. Wea. Rev., 136, 699711.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., and P. H. Lauritzen, 2010: A class of deformational flow test cases for linear transport problems on the sphere. J. Comput. Phys., 229, 88688887.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., S. Thomas, and R. Loft, 2005a: A discontinuous Galerkin global shallow water model. Mon. Wea. Rev., 133, 876888.

  • Nair, R. D., S. Thomas, and R. Loft, 2005b: A discontinuous Galerkin transport scheme on the cubed sphere. Mon. Wea. Rev., 133, 814828.

    • Search Google Scholar
    • Export Citation
  • Nair, R. D., M. N. Levy, and P. H. Lauritzen, 2011: Emerging numerical methods for atmospheric modeling. Numerical Techniques for Global Atmospheric Models, P. Lauritzen et al., Eds., Vol. 80, Lecture Notes in Computational Science and Engineering, Springer, 251–311.

  • Pudykiewicz, J. A., 2006: Numerical solution of the reaction-advection-diffusion equation on the sphere. J. Comput. Phys., 213, 358390.

    • Search Google Scholar
    • Export Citation
  • Pudykiewicz, J. A., 2011: On numerical solution of the shallow water equations with chemical reactions on icosahedral geodesic grid. J. Comput. Phys., 230, 19561991.

    • Search Google Scholar
    • Export Citation
  • Qiu, J., and C.-W. Shu, 2005a: Hermite WENO schemes and their application as limiters for Runge–Kutta discontinuous Galerkin method II: Two dimensional case. Comput. Fluids, 34, 642663.

    • Search Google Scholar
    • Export Citation
  • Qiu, J., and C.-W. Shu, 2005b: Runge--Kutta discontinuous Galerkin method using WENO limiters. SIAM J. Sci. Comput., 26, 907929.

  • 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., 1997: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. Advanced Numerical Approximation of Nonlinear Hyperbolic Equations, A. Quarteroni, Ed., Vol. 1697, Lecture Notes in Mathematics, Springer, 325432.

  • Ullrich, P. A., C. Jablonowski, and B. van Leer, 2010: High-order finite-volume methods for the shallow-water equations on the sphere. J. Comput. Phys., 229, 61046134.

    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber, 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
  • Zhang, X., and C.-W. Shu, 2010: On maximum-principle-satisfying high order schemes for scalar conservation laws. J. Comput. Phys., 229, 30913120.

    • Search Google Scholar
    • Export Citation
  • Zhang, X., and C.-W. Shu, 2011: Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: Survey and new developments. Proc. Roy. Soc., A467, 27522776.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 298 215 7
PDF Downloads 109 51 4