Discontinuous Galerkin methods based on weighted interior penalties for second order PDEs with non-smooth coefficients

13/2007

Discontinuous Galerkin methods based on weighted interior penalties for second order PDEs with non-smooth coefficients

Thursday 30th August 2007

Zunino Paolo

We develop and analyze a Discontinuous Galerkin (DG) method based on weighted interior penalties (WIP) applied to second order PDEs and in particular to advection-diffusion-reaction equations featuring non-smooth and possibly vanishing diffusivity. First of all, looking at the derivation of a DG scheme with a bias to domain decomposition methods, we carefully discuss the set up of the discretization scheme in a general framework putting into evidence the helpful role of the weights and the connection with the well known Local Discontinuous Galerkin schemes (LDG). Then, we address the a-priori and the a-posteriori error analysis of the method, recovering optimal error estimates in suitable norms. By virtue of the introduction of the weighted penalties, these results turn out to be robust with respect to the diffusion parameter. Furthermore, we discuss the derivation of an a-posteriori local error indicator suitable for advection-diffusion-reaction problems with higly variable, locally small diffusivity. Finally, all the theoretical results are illustrated and discussed by means of numerical experiments.