A Systematic Study on Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

This article provides a systematic study for the weak Galerkin (WG) finite element method for second order elliptic problems by exploring polynomial approximations with various degrees for each local element. A typical local WG element is of the form , where is the degree of polynomials in the interior of the element T, is the degree of polynomials on the boundary of T, and is the degree of polynomials employed in the computation of weak gradients or weak first order partial derivatives. A general framework of stability and error estimate is developed for the corresponding numerical solutions. Numerical results are presented to confirm the theoretical results. The work reveals some previously undiscovered strengths of the WG method for second order elliptic problems, and the results are expected to be generalizable to other type of partial differential equations.