Discontinuous Galerkin gradient discretisations for the approximation of second-order differential operators in divergence form

Robert Eymard (1) , Cindy Guichard (2, 3)

1 LAMA - Laboratoire d'Analyse et de Mathématiques Appliquées
2 LJLL - Laboratoire Jacques-Louis Lions
3 ANGE - Numerical Analysis, Geophysics and Ecology

We include in the Gradient Discretisation Method (GDM) framework two numerical schemes based on Discontinuous Galerkin approximations: the Symmetric Interior Penalty Galerkin (SIPG) method, and the scheme obtained by averaging the jumps in the SIPG method. We prove that these schemes meet the main mathematical gradient discretisation properties on any kind of polytopal mesh, by adapting discrete functional analysis properties to our precise geometrical hypotheses. Therefore, these schemes inherit the general convergence properties of the GDM, which hold for instance in the cases of the p−Laplace problem and of the anisotropic and heterogeneous diffusion problem. This is illustrated by simple 1D and 2D numerical examples.

Mots clés: Gradient Discretisation method; Discontinuous Galerkin method; Symmetric Interior Penalty Galerkin scheme; discrete functional analysis; polytopal meshes

Domaines: Analyse numérique [math.NA]

Robert Eymard, Cindy Guichard. Discontinuous Galerkin gradient discretisations for the approximation of second-order differential operators in divergence form. Computational & Applied Mathematics, 2018, 37 (4), pp.:4023-4054. ⟨10.1007/s40314-017-0558-2⟩. ⟨hal-01535147v2⟩