An arbitrary-order method for magnetostatics on polyhedral meshes based on a discrete de Rham sequence

Daniele A. Di Pietro, Jérôme Droniou

Research output: Contribution to journalArticleResearchpeer-review

11 Citations (Scopus)


In this work we develop a discretisation method for the mixed formulation of the magnetostatic problem supporting arbitrary orders and polyhedral meshes. The method is based on a global discrete de Rham (DDR) sequence, obtained by patching the local spaces constructed in [22] by enforcing the single-valuedness of the components attached to the boundary of each element. The first main contribution of this paper is a proof of exactness relations for this global DDR sequence, obtained leveraging the exactness of the corresponding local sequence and a topological assembly of the mesh valid for domains that do not enclose any void. The second main contribution is the formulation and well-posedness analysis of the method, which includes the proof of uniform Poincaré inequalities for the discrete divergence and curl operators. The convergence rate in the natural energy norm is numerically evaluated on standard and polyhedral meshes. When the DDR sequence of degree k≥0 is used, the error converges as hk+1, with h denoting the mesh size.

Original languageEnglish
Article number109991
Number of pages28
JournalJournal of Computational Physics
Publication statusPublished - 15 Mar 2021


  • Compatible discretisations
  • Discrete de Rham
  • Magnetostatics
  • Mixed methods
  • Polyhedral methods

Cite this