Editorial Type:
Article
Article Type:
Research Article

Quadrature-Free Implementation of a Discontinuous Galerkin Global Shallow-Water Model via Flux Correction Procedure

Ramachandran D. Nair Institute for Mathematics Applied to Geosciences, National Center for Atmospheric Research,* Boulder, Colorado

Search for other papers by Ramachandran D. Nair in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Apr 2015
DOI:
https://doi.org/10.1175/MWR-D-14-00174.1
Page(s):
1335–1346
Restricted access

Abstract

The discontinuous Galerkin (DG) discretization relies on an integral (weak) formulation of the hyperbolic conservation law, which leads to the evaluation of several surface and line integrals for multidimensional problems. An alternative formulation of the DG method is possible under the flux reconstruction (FR) framework, where the equations are solved in differential form and the discretization is free from quadrature rules, resulting in computationally efficient algorithms. The author has implemented a quadrature-free form of the nodal DG method based on the FR approach combined with spectral differencing (SD), in a shallow-water (SW) model employing cubed-sphere geometry. The performance of the SD model is compared with the regular nodal DG variant of the SW model using several benchmark tests, including a viscous test case. A positivity-preserving local filter is tested for SD advection that removes spurious oscillations while being conservative and accurate. In this implementation, the SD formulation is found to be 18% faster than the DG method for inviscid SW tests cases and 24% faster for the viscous case. The results obtained by the SD formulation are on par with the regular nodal DG formulation in terms of accuracy and convergence.

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, 1850 Table Mesa Dr., Boulder, CO 80305. E-mail: rnair@ucar.edu

Abstract

The discontinuous Galerkin (DG) discretization relies on an integral (weak) formulation of the hyperbolic conservation law, which leads to the evaluation of several surface and line integrals for multidimensional problems. An alternative formulation of the DG method is possible under the flux reconstruction (FR) framework, where the equations are solved in differential form and the discretization is free from quadrature rules, resulting in computationally efficient algorithms. The author has implemented a quadrature-free form of the nodal DG method based on the FR approach combined with spectral differencing (SD), in a shallow-water (SW) model employing cubed-sphere geometry. The performance of the SD model is compared with the regular nodal DG variant of the SW model using several benchmark tests, including a viscous test case. A positivity-preserving local filter is tested for SD advection that removes spurious oscillations while being conservative and accurate. In this implementation, the SD formulation is found to be 18% faster than the DG method for inviscid SW tests cases and 24% faster for the viscous case. The results obtained by the SD formulation are on par with the regular nodal DG formulation in terms of accuracy and convergence.

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, 1850 Table Mesa Dr., Boulder, CO 80305. E-mail: rnair@ucar.edu
Save
Cover Monthly Weather Review
Monthly Weather Review

  • Allaneau, Y., and A. Jameson, 2011: Connections between the filtered discontinuous Galerkin method and the flux reconstruction approach to high-order discretizations. Comput. Methods Appl. Mech. Eng., 200, 3628–3636, doi:10.1016/j.cma.2011.08.019.

    • Search Google Scholar
    • Export Citation

  • Arnold, D. N., F. Brezzi, B. Cockburn, and L. D. Marini, 2002: Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal., 39, 1749–1779, doi:10.1137/S0036142901384162.

    • Search Google Scholar
    • Export Citation

  • Atkins, H., and C.-W. Shu, 1998: Quadrature-free implementation of the discontinuous Galerkin method for hyperbolic equations. AIAA J., 36, 775–782, doi:10.2514/2.436.

    • Search Google Scholar
    • Export Citation

  • Bao, L., R. D. Nair, and H. M. Tufo, 2014: A mass and momentum flux-form high-order discontinuous Galerkin shallow water model on the cubed-sphere. J. Comput. Phys., 271, 224–243, doi:10.1016/j.jcp.2013.11.033.

    • Search Google Scholar
    • Export Citation

  • Bassi, F., and S. Rebay, 1997: A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier–Stokes equations. J. Comput. Phys., 131, 267–279, doi:10.1006/jcph.1996.5572.

    • 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, 5019–5044, doi:10.1016/j.jcp.2008.01.033.

    • Search Google Scholar
    • Export Citation

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

  • Cockburn, B., and C.-W. Shu, 1998: The local discontinuous Galerkin for convection diffusion systems. SIAM J. Numer. Anal., 35, 2440–2463, doi:10.1137/S0036142997316712.

    • Search Google Scholar
    • Export Citation

  • De Grazia, D., G. Mengaldo, D. Moxey, P. E. Vincent, and S. J. Sherwin, 2014: Connections between the discontinuous Galerkin method and high-order flux reconstruction schemes. Int. J. Numer. Methods Fluids, 75, 860–877, doi:10.1002/fld.3915.

    • Search Google Scholar
    • Export Citation

  • Deville, M. O., P. F. Fisher, and E. M. Mund, 2002: High-Order Methods for Incompressible Flow. Cambridge University Press, 499 pp.

  • Gao, H., Z. J. Wang, and H. T. Huynh, 2013: Differential formulation of discontinuous Galerkin and related methods for the Navier–Stokes equations. Commun. Comput. Phys., 13, 1013–1044.

    • Search Google Scholar
    • Export Citation

  • Giraldo, F., J. Hesthaven, and T. Warburton, 2002: Nodal high-order discontinuous Galerkin methods for spherical shallow water equations. J. Comput. Phys., 181, 499–525, doi:10.1006/jcph.2002.7139.

    • 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, 89–112, doi:10.1137/S003614450036757X.

    • Search Google Scholar
    • Export Citation

  • Huynh, H. T., 2007: A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods. AIAA Paper 2007-4079, 18th AIAA Computational Fluid Dynamics Conf., Miami, FL, AIAA, 1–42.

  • Huynh, H. T., 2009: A reconstruction approach to high-order schemes including discontinuous Galerkin for diffusion. AIAA Paper 2009-403, 47th AIAA Aerospace Science Meeting, Orlando, FL, AIAA, 1–34.

  • Huynh, H. T., 2014: On formulation of discontinuous Galerkin and related methods for conservation laws. Tech. Note NASA/TM-2014-218135, NASA, 1–28.

  • Kopriva, D. A., and J. H. Kolias, 1996: A conservative staggered-grid Chebyshev multidomain method for compressible flows. J. Comput. Phys., 125, 244–261, doi:10.1006/jcph.1996.0091.

    • Search Google Scholar
    • Export Citation

  • Liang, C., C. Cox, and M. Plesniak, 2013: A comparison of computational efficiencies of spectral difference method and correction procedure via reconstruction. J. Comput. Phys., 239, 138–146, doi:10.1016/j.jcp.2013.01.001.

    • Search Google Scholar
    • Export Citation

  • May, G., 2011: On the connection between the spectral difference method and the discontinuous Galerkin method. Commun. Comput. Phys., 9, 1071–1080.

    • Search Google Scholar
    • Export Citation

  • Nair, R. D., 2009: Diffusion experiments with a global Discontinuous Galerkin shallow-water model. Mon. Wea. Rev., 137, 3339–3350, doi:10.1175/2009MWR2843.1.

    • Search Google Scholar
    • Export Citation

  • Nair, R. D., S. J. Thomas, and R. D. Loft, 2005: A discontinuous Galerkin global shallow water model. Mon. Wea. Rev., 133, 876–888, doi:10.1175/MWR2903.1.

    • Search Google Scholar
    • Export Citation

  • Nair, R. D., H.-W. Choi, and H. M. Tufo, 2009: Computational aspects of a scalable high-order discontinuous Galerkin atmospheric dynamical core. Comput. Fluids, 38, 309–319, doi:10.1016/j.compfluid.2008.04.006.

    • 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. H. Lauritzen et al., Eds., Vol. 80, Springer-Verlag, 189–250.

  • Sadourny, R., 1972: Conservative finite-difference approximations of the primitive equations on quasi-uniform spherical grids. Mon. Wea. Rev., 100, 136–144, doi:10.1175/1520-0493(1972)100<0136:CFAOTP>2.3.CO;2.

    • Search Google Scholar
    • Export Citation

  • Williamson, D., J. Drake, J. Hack, R. Jacob, and P. Swarztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys., 102, 211–224, doi:10.1016/S0021-9991(05)80016-6.

    • 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, 3091–3120, doi:10.1016/j.jcp.2009.12.030.

    • Search Google Scholar
    • Export Citation

  • Zhang, Y., and R. D. Nair, 2012: A nonoscillatory discontinuous Galerkin transport scheme on the cubed sphere. Mon. Wea. Rev., 140, 3106–3126, doi:10.1175/MWR-D-11-00287.1.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 315 97 11
PDF Downloads 236 84 12