Communications on Applied Mathematics and Computation ›› 2020, Vol. 2 ›› Issue (3): 461-486.doi: 10.1007/s42967-019-00049-2

• ORIGINAL PAPER • 上一篇    下一篇

Superconvergence and the Numerical Flux: a Study Using the Upwind-Biased Flux in Discontinuous Galerkin Methods

Daniel J. Frean1, Jennifer K. Ryan1,2   

  1. 1 School of Mathematics, University of East Anglia, Norwich, UK;
    2 Applied Mathematics and Statistics, Colorado School of Mines, Golden, CO, USA
  • 收稿日期:2019-04-03 修回日期:2019-10-02 出版日期:2020-09-20 发布日期:2020-05-12
  • 通讯作者: Jennifer K. Ryan E-mail:jkryan@mines.edu

Superconvergence and the Numerical Flux: a Study Using the Upwind-Biased Flux in Discontinuous Galerkin Methods

Daniel J. Frean1, Jennifer K. Ryan1,2   

  1. 1 School of Mathematics, University of East Anglia, Norwich, UK;
    2 Applied Mathematics and Statistics, Colorado School of Mines, Golden, CO, USA
  • Received:2019-04-03 Revised:2019-10-02 Online:2020-09-20 Published:2020-05-12
  • Contact: Jennifer K. Ryan E-mail:jkryan@mines.edu

摘要: One of the benefcial properties of the discontinuous Galerkin method is the accurate wave propagation properties. That is, the semi-discrete error has dissipation errors of order 2k + 1 (≤ Ch2k+1) and order 2k + 2 for dispersion (≤ Ch2k+2). Previous studies have concentrated on the order of accuracy, and neglected the important role that the error constant, C, plays in these estimates. In this article, we show the important role of the error constant in the dispersion and dissipation error for discontinuous Galerkin approximation of polynomial degree k, where k = 0, 1, 2, 3. This gives insight into why one may want a more centred fux for a piecewise constant or quadratic approximation than for a piecewise linear or cubic approximation. We provide an explicit formula for these error constants. This is illustrated through one particular fux, the upwind-biased fux introduced by Meng et al., as it is a convex combination of the upwind and downwind fuxes. The studies of wave propagation are typically done through a Fourier ansatz. This higher order Fourier information can be extracted using the smoothness-increasing accuracy-conserving (SIAC) flter. The SIAC flter ties the higher order Fourier information to the negative-order norm in physical space. We show that both the proofs of the ability of the SIAC flter to extract extra accuracy and numerical results are unafected by the choice of fux.

关键词: Discontinuous Galerkin, Smoothness-increasing accuracy-conserving (SIAC) fltering, Superconvergence

Abstract: One of the benefcial properties of the discontinuous Galerkin method is the accurate wave propagation properties. That is, the semi-discrete error has dissipation errors of order 2k + 1 (≤ Ch2k+1) and order 2k + 2 for dispersion (≤ Ch2k+2). Previous studies have concentrated on the order of accuracy, and neglected the important role that the error constant, C, plays in these estimates. In this article, we show the important role of the error constant in the dispersion and dissipation error for discontinuous Galerkin approximation of polynomial degree k, where k = 0, 1, 2, 3. This gives insight into why one may want a more centred fux for a piecewise constant or quadratic approximation than for a piecewise linear or cubic approximation. We provide an explicit formula for these error constants. This is illustrated through one particular fux, the upwind-biased fux introduced by Meng et al., as it is a convex combination of the upwind and downwind fuxes. The studies of wave propagation are typically done through a Fourier ansatz. This higher order Fourier information can be extracted using the smoothness-increasing accuracy-conserving (SIAC) flter. The SIAC flter ties the higher order Fourier information to the negative-order norm in physical space. We show that both the proofs of the ability of the SIAC flter to extract extra accuracy and numerical results are unafected by the choice of fux.

Key words: Discontinuous Galerkin, Smoothness-increasing accuracy-conserving (SIAC) fltering, Superconvergence

中图分类号: