Smoothness-Increasing Accuracy-Conserving (SIAC) filters for derivative approximations of discontinuous Galerkin (DG) solutions over nonuniform meshes and near boundaries

详细信息    查看全文

• 作者：X. Li^a ; ^{X.Li-2@tudelft.nl" class="auth_mail" title="E-mail the corresponding author} ; J.K. Ryan^b ; ^{Jennifer.Ryan@uea.ac.uk" class="auth_mail" title="E-mail the corresponding author} ; R.M. Kirby^c ; ^{kirby@cs.utah.edu" class="auth_mail" title="E-mail the corresponding author} ; C. Vuik^a ; ^{C.Vuik@tudelft.nl" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Discontinuous Galerkin method ; Post-processing ; SIAC filtering ; Superconvergence ; Nonuniform meshes ; Boundaries
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2016
• 出版时间：1 March 2016
• 年：2016
• 卷：294
• 期：Complete
• 页码：275-296
• 全文大小：4188 K

文摘

Accurate approximations for the derivatives are usually required in many application areas such as biomechanics, chemistry and visualization applications. With the help of Smoothness-Increasing Accuracy-Conserving (SIAC) filtering, one can enhance the derivatives of a discontinuous Galerkin solution. However, current investigations of derivative filtering are limited to uniform meshes and periodic boundary conditions, which do not meet practical requirements. The purpose of this paper is twofold: to extend derivative filtering to nonuniform meshes and propose position-dependent derivative filters to handle filtering near the boundaries. Through analyzing the error estimates for SIAC filtering, we extend derivative filtering to nonuniform meshes by changing the scaling of the filter. For filtering near boundaries, we discuss the advantages and disadvantages of two existing position-dependent filters and then extend them to position-dependent derivative filters, respectively. Further, we prove that with the position-dependent derivative filters, the filtered solutions can obtain a better accuracy rate compared to the original discontinuous Galerkin approximation with arbitrary derivative orders over nonuniform meshes. One- and two-dimensional numerical results are provided to support the theoretical results and demonstrate that the position-dependent derivative filters, in general, enhance the accuracy of the solution for both uniform and nonuniform meshes.