Abstract
This paper develops interior penalty discontinuous Galerkin (IP-DG) methods to approximate W2 , p strong solutions of second order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. The proposed IP-DG methods are closely related to the IP-DG methods for advection-diffusion equations, and they are easy to implement on existing standard IP-DG software platforms. It is proved that the proposed IP-DG methods have unique solutions and converge with optimal rate to the W2 , p strong solution in a discrete W2 , p-norm. The crux of the analysis is to establish a DG discrete counterpart of the Calderon–Zygmund estimate and to adapt a freezing coefficient technique used for the PDE analysis at the discrete level. To obtain such a crucial estimate, we need to establish broken W1 , p-norm error estimates for IP-DG approximations of constant coefficient elliptic PDEs, which is also of independent interest. Numerical experiments are provided to gauge the performance of the proposed IP-DG methods and to validate the theoretical convergence results.
Original language | English |
---|---|
Pages (from-to) | 1651-1676 |
Number of pages | 26 |
Journal | Journal of Scientific Computing |
Volume | 74 |
Issue number | 3 |
DOIs | |
State | Published - Mar 1 2018 |
Externally published | Yes |
Funding
The work of the Xiaobing Feng and Stefan Schnake was partial supported by the NSF through Grant DMS-1318486 and the work of the Michael Neilan was partially supported by the NSF Grant DMS-1417980 and the Alfred Sloan Foundation.
Keywords
- 35J25
- 65N12
- 65N30