Interior penalty continuous and discontinuous finite element approximations of hyperbolic equations

In this paper we present in a unified setting the continuous and discontinuous Galerkin methods for the numerical approximation of the scalar hyperbolic equation. Both methods are stabilized by the interior penalty method, more precisely by the jump of the gradient across element faces in the continuous case whereas in the discontinuous case the stabilization of the jump of the solution and optionally of its gradient is required to achieve optimal convergence. We prove that the solution in the case of the continuous Galerkin approach can be considered as a limit of the discontinuous one when the stabilization parameter associated with the penalization of the solution jump tends to infinity. As a consequence, the limit of the numerical flux of the discontinuous method yields a numerical flux for the continuous method as well. Numerical results will highlight the theoretical results that are proven in this paper