Communications on Applied Mathematics and Computation ›› 2022, Vol. 4 ›› Issue (2): 634-656.doi: 10.1007/s42967-021-00133-6

• ORIGINAL PAPERS • 上一篇    下一篇

Discontinuous Galerkin Methods for a Class of Nonvariational Problems

Andreas Dedner1, Tristan Pryer2   

  1. 1 Mathematics Institute, University of Warwick, Coventry CV4 7 AL, UK;
    2 Department of Mathematical Sciences, University of Bath, Bath BA2 7 AY, UK
  • 收稿日期:2020-09-22 修回日期:2021-03-01 出版日期:2022-06-20 发布日期:2022-04-29
  • 通讯作者: Andreas Dedner, Tristan Pryer E-mail:A.S.Dedner@warwick.ac.uk;tmp38@bath.ac.uk

Discontinuous Galerkin Methods for a Class of Nonvariational Problems

Andreas Dedner1, Tristan Pryer2   

  1. 1 Mathematics Institute, University of Warwick, Coventry CV4 7 AL, UK;
    2 Department of Mathematical Sciences, University of Bath, Bath BA2 7 AY, UK
  • Received:2020-09-22 Revised:2021-03-01 Online:2022-06-20 Published:2022-04-29
  • Contact: Andreas Dedner, Tristan Pryer E-mail:A.S.Dedner@warwick.ac.uk;tmp38@bath.ac.uk

摘要: We extend the fnite element method introduced by Lakkis and Pryer (SIAM J. Sci. Comput. 33(2):786-801, 2011) to approximate the solution of second-order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the "fnite element Hessian" as an auxiliary variable in the formulation. Representing the fnite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble; thus, this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach. We conduct a stability and consistency analysis making use of the unifed frameworkset out in Arnold et al. (SIAM J. Numer. Anal. 39(5):1749-1779, 2001/2002). We also give an a posteriori analysis of the method in the case where the problem has a strong solution. The analysis applies to any consistent representation of the fnite element Hessian, and thus is applicable to the previous works making use of continuous Galerkin approximations. Numerical evidence is presented showing that the method works well also in a more general setting.

关键词: Nonvariational problems, Discontinuous Galerkin, Error estimates, Adaptivity

Abstract: We extend the fnite element method introduced by Lakkis and Pryer (SIAM J. Sci. Comput. 33(2):786-801, 2011) to approximate the solution of second-order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin (DG) framework. This is done by viewing the "fnite element Hessian" as an auxiliary variable in the formulation. Representing the fnite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems. Furthermore, the system matrix is very easy to assemble; thus, this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach. We conduct a stability and consistency analysis making use of the unifed frameworkset out in Arnold et al. (SIAM J. Numer. Anal. 39(5):1749-1779, 2001/2002). We also give an a posteriori analysis of the method in the case where the problem has a strong solution. The analysis applies to any consistent representation of the fnite element Hessian, and thus is applicable to the previous works making use of continuous Galerkin approximations. Numerical evidence is presented showing that the method works well also in a more general setting.

Key words: Nonvariational problems, Discontinuous Galerkin, Error estimates, Adaptivity

中图分类号: