The discontinuous Petrov–Galerkin method for elliptic problems

Abstract

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual–primal mixed (DPM) finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge–Kutta schemes of the collocation type exhibiting optimal behavior and we develop the stability analysis. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.