{"title":"Optimal error estimates of discontinuous Galerkin methods with generalized fluxes for wave equations on unstructured meshes","authors":"Zheng Sun, Y. Xing","doi":"10.1090/mcom/3605","DOIUrl":null,"url":null,"abstract":"L2 stable discontinuous Galerkin method with a family of numerical fluxes was studied for the one-dimensional wave equation by Cheng, Chou, Li, and Xing in [Math. Comp. 86 (2017), pp. 121–155]. Although optimal convergence rates were numerically observed with wide choices of parameters in the numerical fluxes, their error estimates were only proved for a sub-family with the construction of a local projection. In this paper, we first complete the one-dimensional analysis by providing optimal error estimates that match all numerical observations in that paper. The key ingredient is to construct an optimal global projection with the characteristic decomposition. We then extend the analysis on optimal error estimate to multidimensions by constructing a global projection on unstructured meshes, which can be considered as a perturbation away from the local projection studied by Cockburn, Gopalakrishnan, and Sayas in [Math. Comp. 79 (2010), pp. 1351–1367] for hybridizable discontinuous Galerkin methods. As a main contribution, we use a novel energy argument to prove the optimal approximation property of the global projection. This technique does not require explicit assembly of the matrix for the perturbed terms and hence can be easily used for unstructured meshes in multidimensions. Finally, numerical tests in two dimensions are provided to validate our analysis is sharp and at least one of the unknowns will degenerate to suboptimal rates if the assumptions are not satisfied.","PeriodicalId":18301,"journal":{"name":"Math. Comput. Model.","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2020-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"8","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Math. Comput. Model.","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1090/mcom/3605","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 8
Abstract
L2 stable discontinuous Galerkin method with a family of numerical fluxes was studied for the one-dimensional wave equation by Cheng, Chou, Li, and Xing in [Math. Comp. 86 (2017), pp. 121–155]. Although optimal convergence rates were numerically observed with wide choices of parameters in the numerical fluxes, their error estimates were only proved for a sub-family with the construction of a local projection. In this paper, we first complete the one-dimensional analysis by providing optimal error estimates that match all numerical observations in that paper. The key ingredient is to construct an optimal global projection with the characteristic decomposition. We then extend the analysis on optimal error estimate to multidimensions by constructing a global projection on unstructured meshes, which can be considered as a perturbation away from the local projection studied by Cockburn, Gopalakrishnan, and Sayas in [Math. Comp. 79 (2010), pp. 1351–1367] for hybridizable discontinuous Galerkin methods. As a main contribution, we use a novel energy argument to prove the optimal approximation property of the global projection. This technique does not require explicit assembly of the matrix for the perturbed terms and hence can be easily used for unstructured meshes in multidimensions. Finally, numerical tests in two dimensions are provided to validate our analysis is sharp and at least one of the unknowns will degenerate to suboptimal rates if the assumptions are not satisfied.
《数学》杂志Cheng、Chou、Li和Xing研究了一维波动方程的具有一组数值通量的L2稳定不连续Galerkin方法。[p. 86 (2017), pp. 121-155]。虽然在数值通量参数选择范围广泛的情况下,在数值上观察到最优的收敛速率,但它们的误差估计仅通过构造局部投影来证明。在本文中,我们首先通过提供与文中所有数值观测相匹配的最佳误差估计来完成一维分析。关键是利用特征分解构造最优的全局投影。然后,我们通过在非结构化网格上构造一个全局投影,将最优误差估计的分析扩展到多维,该投影可以被认为是远离Cockburn, Gopalakrishnan和Sayas在[Math]中研究的局部投影的扰动。Comp. 79 (2010), pp. 1351-1367]杂交不连续Galerkin方法。作为主要贡献,我们使用了一个新的能量参数来证明全局投影的最优逼近性质。该技术不需要对扰动项进行矩阵的显式装配,因此可以很容易地用于多维非结构化网格。最后,给出了二维的数值实验来验证我们的分析是敏锐的,如果假设不满足,至少有一个未知数会退化到次优速率。