Uniform-in-time convergence of numerical methods for non-linear degenerate parabolic equations

Jerome Droniou, Robert Eymard

Research output: Contribution to journalArticleResearchpeer-review

38 Citations (Scopus)


Gradient schemes is a framework that enables the unified convergence analysis of many numerical methods for elliptic and parabolic partial differential equations: conforming and non-conforming finite element, mixed finite element and finite volume methods. We show here that this framework can be applied to a family of degenerate non-linear parabolic equations (which contain in particular the Richards’, Stefan’s and Leray–Lions’ models), and we prove a uniform-in-time strong-in-space convergence result for the gradient scheme approximations of these equations. In order to establish this convergence, we develop several discrete compactness tools for numerical approximations of parabolic models, including a discontinuous Ascoli–Arzelà theorem and a uniform-in-time weak-in-space discrete Aubin–Simon theorem. The model’s degeneracies, which occur both in the time and space derivatives, also requires us to develop a discrete compensated compactness result.

Original languageEnglish
Pages (from-to)721-766
Number of pages46
JournalNumerische Mathematik
Issue number4
Publication statusPublished - 1 Apr 2016


  • 35K65
  • 46N40
  • 65M12

Cite this