Low Order Discontinuous Galerkin Methods for Second Order Elliptic Problems

We consider DG-methods for second order scalar elliptic problems using piecewise affine approximation in two or three space dimensions. We prove that both the symmetric and the nonsymmetric versions of the DG-method have regular system matrices without penalization of the interelement solution jumps provided boundary conditions are imposed in a certain weak manner. Optimal convergence is proved for sufficiently regular meshes and data. We then propose a DG-method using piecewise affine functions enriched with quadratic bubbles. Using this space we prove optimal convergence in the energy norm for both a symmetric and nonsymmetric DG-method without stabilization. All of these proposed methods share the feature that they conserve mass locally independent of the penalty parameter.