pith. sign in

arxiv: 2604.13921 · v1 · submitted 2026-04-15 · 🧮 math.NA · cs.NA

A higher-order dual cell method for time-domain Maxwell equations

Lorenzo Codecasa , Bernard Kapidani , Joachim Sch\"oberl , Markus Wess This is my paper

Pith reviewed 2026-05-10 12:23 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Maxwell equationsdual cell methodhigher-order discretizationtime-domainbarycentric gridscurl-conforming spacesexplicit time integration
0
0 comments X

The pith

A higher-order dual cell method places electric and magnetic fields on mutually dual barycentric grids to reach arbitrary-order accuracy for time-domain Maxwell equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper extends the dual cell method to higher orders for three-dimensional time-domain Maxwell equations. It approximates the electric field on one barycentric grid and the magnetic field on its dual using curl-conforming polynomial spaces built from tensor-product Gauss-Radau interpolation. This construction produces block-diagonal mass matrices and sparse discrete curl operators that permit explicit time integration while exactly preserving a discrete energy identity. The spaces are made compatible across primal and dual meshes and keep tangential continuity after mapping from reference to physical elements. Numerical experiments on unstructured tetrahedral meshes confirm high-order convergence rates and the absence of spurious modes.

Core claim

Curl-conforming polynomial spaces can be constructed on mutually dual barycentric grids in three dimensions via tensor-product Gauss-Radau interpolation. These spaces ensure compatibility between the discrete curl operators on the primal and dual meshes, produce block-diagonal mass matrices, and support a discrete energy identity that holds exactly for the semi-discrete system. The resulting formulation therefore converges at arbitrary order, avoids spurious modes, and retains optimal sparsity on unstructured tetrahedral meshes.

What carries the argument

Compatible curl-conforming polynomial spaces on mutually dual barycentric grids, constructed via tensor-product Gauss-Radau interpolation that preserves tangential continuity under reference-to-physical mappings.

If this is right

  • The semi-discrete system converges at arbitrary order for the electric and magnetic fields.
  • The discrete spectrum contains no spurious modes beyond the expected null space.
  • Block-diagonal mass matrices allow explicit time integration without solving linear systems at each step.
  • The discrete energy identity holds exactly at the semi-discrete level.
  • Optimal sparsity is retained in the discrete curl operators for computational efficiency.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The energy-preserving property implies long-time stability for wave-propagation problems without added numerical dissipation.
  • The same dual-grid construction could be applied to other first-order hyperbolic systems that benefit from structure preservation.
  • Extension to general polyhedral meshes would widen the range of geometries that can be treated without remeshing.

Load-bearing premise

That compatible curl-conforming polynomial spaces can be constructed on mutually dual barycentric grids while preserving tangential continuity under reference-to-physical mappings.

What would settle it

A sequence of uniformly refined tetrahedral meshes for a closed cavity problem in which the observed L2 convergence rate of the electric field falls below the polynomial degree plus one or spurious zero-frequency modes appear in the discrete spectrum.

Figures

Figures reproduced from arXiv: 2604.13921 by Bernard Kapidani, Joachim Sch\"oberl, Lorenzo Codecasa, Markus Wess.

Figure 1
Figure 1. Figure 1: Entities of the primal mesh (red) and their dual counterparts (blue) for a mesh consisting of one [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: A primal (tetrahedral) and dual (polyhedral) cell. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Reference-to-physical trilinear mapping FK from the cube Kp to a generic polyhedral subcell K P K. Hpcurl, Trq it results in ÿ TrPTr ż Tr e ¨ εBtE dV “ ÿ TrPTr ˜ ´ ż BTrzBΩ e ˆ H ¨ n dS ` ż Tr curl e ¨ H dV ¸ , (1) where n is the unit vector function outward normal to BTr. Similarly for each h P Hpcurl, T q we obtain ÿ TPT ż T h ¨ µBtH dV “ ´ ÿ TPT ˆ ´ ż BT h ˆ E ¨ n dS ` ż T curl h ¨ E dV ˙ . (2) Let now … view at source ↗
Figure 4
Figure 4. Figure 4: Support of different types of basis functions on the dual mesh [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Support of different types of basis functions on the primal mesh [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Sparsity of the resulting mass matrices for a fixed mesh and increasing polynomial order. [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Convergence of the third, non-trivial eigenvalue. [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Convergence of the eigenfunctions corresponding to the second, non-trivial eigenvalue [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The CFL condition with respect to mesh-size and polynomial order. [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Typical meshes of the waveguide problems with perfectly matched layers. [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Snapshots of the plain waveguide problem. [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: L2pΩintq-error of the E field (cf [PITH_FULL_IMAGE:figures/full_fig_p015_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Convergence of the L2pr0, 1s; Ωintq-error of the time domain E field (cf [PITH_FULL_IMAGE:figures/full_fig_p015_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Waveguide with spherical inclusion 16 [PITH_FULL_IMAGE:figures/full_fig_p016_15.png] view at source ↗
read the original abstract

We present a higher-order extension of the dual cell method for the time-domain Maxwell equations in three spatial dimensions. The approach builds upon a variational reinterpretation of the Finite Integration Technique on dual meshes and generalises a previously developed two-dimensional high-order formulation. The electric and magnetic fields are discretised on mutually dual barycentric grids using curl-conforming polynomial spaces constructed via tensor-product Gauss--Radau interpolation. The resulting semi-discrete formulation yields block-diagonal mass matrices and sparse discrete curl operators, enabling explicit time integration while preserving a discrete energy identity. Special attention is devoted to the construction of compatible approximation spaces on the three-dimensional primal and dual meshes, the reference-to-physical element mappings, and the preservation of tangential continuity. We show that the method achieves arbitrary-order convergence, avoids spurious modes, and maintains optimal sparsity properties. Numerical experiments confirm spectral correctness, high-order accuracy, and computational efficiency on unstructured tetrahedral meshes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript presents a higher-order extension of the dual cell method for the time-domain Maxwell equations in three dimensions. It builds upon a variational reinterpretation of the Finite Integration Technique on dual meshes and generalizes a previously developed two-dimensional high-order formulation. The electric and magnetic fields are discretized on mutually dual barycentric grids using curl-conforming polynomial spaces constructed via tensor-product Gauss-Radau interpolation. The resulting semi-discrete formulation yields block-diagonal mass matrices and sparse discrete curl operators, enabling explicit time integration while preserving a discrete energy identity. Special attention is given to the construction of compatible approximation spaces on the three-dimensional primal and dual meshes, the reference-to-physical element mappings, and the preservation of tangential continuity. The authors claim that the method achieves arbitrary-order convergence, avoids spurious modes, and maintains optimal sparsity properties, with numerical experiments confirming spectral correctness, high-order accuracy, and computational efficiency on unstructured tetrahedral meshes.

Significance. If the central constructions and claims hold, this would represent a meaningful advance in compatible discretizations for Maxwell's equations. It combines the structural advantages of dual-grid methods (energy preservation, explicit time stepping, sparsity) with arbitrary-order accuracy on unstructured meshes, addressing a practical need in computational electromagnetics where high-order schemes that remain stable and efficient are valuable. The emphasis on preserving tangential continuity under mappings and forming a compatible discrete de Rham sequence aligns with established principles in discrete exterior calculus and finite element exterior calculus.

minor comments (2)
  1. The abstract and introduction would benefit from a brief forward reference to the specific sections containing the error analysis and convergence proofs, as the claims of arbitrary-order convergence are central to the contribution.
  2. Notation for the primal and dual barycentric grids (e.g., consistent use of subscripts or superscripts for primal/dual quantities) could be standardized earlier to aid readability across the construction and implementation sections.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive evaluation of the manuscript and the recommendation to accept. The review accurately captures the key contributions of the higher-order dual cell method, including the construction of compatible spaces, preservation of discrete energy, and numerical performance on unstructured meshes.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's core construction uses tensor-product Gauss-Radau interpolation to build curl-conforming spaces on mutually dual barycentric grids that form a compatible discrete de Rham sequence while preserving tangential continuity under mappings. This yields block-diagonal mass matrices and sparse curl operators by design of the spaces, enabling explicit time integration and a discrete energy identity. The claims of arbitrary-order convergence and absence of spurious modes are then verified through analysis and numerical experiments on tetrahedral meshes, without any step reducing a prediction or result to a fitted parameter, self-definition, or load-bearing self-citation by construction. The variational reinterpretation of the Finite Integration Technique draws on standard discrete exterior calculus ideas that remain independent of the present work's fitted values or prior author results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the existence of compatible polynomial spaces on dual meshes and the preservation of tangential continuity under mappings; these are standard domain assumptions in finite-element exterior calculus rather than new postulates.

axioms (2)
  • domain assumption Curl-conforming polynomial spaces on primal and dual barycentric grids admit block-diagonal mass matrices and sparse discrete curl operators
    Invoked to enable explicit time integration and energy conservation
  • domain assumption Reference-to-physical element mappings preserve tangential continuity of the discrete fields
    Required for well-posedness on unstructured tetrahedral meshes

pith-pipeline@v0.9.0 · 5460 in / 1314 out tokens · 25922 ms · 2026-05-10T12:23:42.600095+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    K.S. Yee. Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media.IEEE Transactions on Antennas and Propagation, 14(3):302–307, May 1966

    work page 1966

  2. [2]

    Taflove, A

    A. Taflove, A. Oskooi, and S. G. Johnson, editors.Advances in FDTD computational electrodynamics: photonics and nanotechnology. Artech House, Boston, 2013

    work page 2013

  3. [3]

    John Wiley & Sons Inc, Hoboken

    Jian-Ming Jin.The finite element method in electromagnetics. John Wiley & Sons Inc, Hoboken. New Jersey, third edition edition, 2014

    work page 2014

  4. [4]

    J. C. Nedelec. Mixed finite elements in$\mathbbR3$.Numerische Mathematik, 35(3):315–341, Septem- ber 1980. 19

    work page 1980

  5. [5]

    Sch¨ oberl and Sabine Zaglmayr

    J. Sch¨ oberl and Sabine Zaglmayr. High order N´ ed´ elec elements with local complete sequence properties. COMPEL - The international journal for computation and mathematics in electrical and electronic engineering, 24(2):374–384, January 2005

    work page 2005

  6. [6]

    Hesthaven and T

    J.S. Hesthaven and T. Warburton.Nodal Discontinuous Galerkin Methods, volume 54 ofTexts in Applied Mathematics. Springer New York, New York, NY, 2008

    work page 2008

  7. [7]

    Arnold, Franco Brezzi, B

    Douglas N. Arnold, Franco Brezzi, B. Cockburn, and L. Donatella Marini. Unified Analysis of Dis- continuous Galerkin Methods for Elliptic Problems.SIAM J. Numer. Anal., 39(5):1749–1779, January 2002

    work page 2002

  8. [8]

    T. Weiland. Time Domain Electromagnetic Field Computation with Finite Difference Methods.In- ternational Journal of Numerical Modelling: Electronic Networks, Devices and Fields, 9(4):295–319, 1996

    work page 1996

  9. [9]

    Codecasa and M

    L. Codecasa and M. Politi. Explicit, Consistent, and Conditionally Stable Extension of FD-TD to Tetrahedral Grids by FIT.IEEE Transactions on Magnetics, 44(6):1258–1261, June 2008

    work page 2008

  10. [10]

    Codecasa, B

    L. Codecasa, B. Kapidani, R. Specogna, and F. Trevisan. Novel FDTD Technique Over Tetrahedral Grids for Conductive Media.IEEE Transactions on Antennas and Propagation, 66(10):5387–5396, October 2018

    work page 2018

  11. [11]

    Kapidani, L

    B. Kapidani, L. Codecasa, and R. Specogna. The Time-Domain Cell Method Is a Coupling of Two Explicit Discontinuous Galerkin Schemes With Continuous Fluxes.IEEE Transactions on Magnetics, 56(1):1–4, January 2020

    work page 2020

  12. [12]

    Brezzi, M

    F. Brezzi, M. Fortin, and L. D. Marini. Error analysis of piecewise constant pressure approximations of Darcy’s law.Comput. Methods Appl. Mech. Engrg., 195(13-16):1547–1559, 2006

    work page 2006

  13. [13]

    Lee and Ragnar Winther

    Jeonghun J. Lee and Ragnar Winther. Local coderivatives and approximation of Hodge Laplace prob- lems.Math. Comp., 87(314):2709–2735, 2018

    work page 2018

  14. [14]

    A Yee-like finite-element scheme for Maxwell’s equations on unstruc- tured grids.IMA J

    Bogdan Radu and Herbert Egger. A Yee-like finite-element scheme for Maxwell’s equations on unstruc- tured grids.IMA J. Numer. Anal., 45(2):1028–1053, 2025

    work page 2025

  15. [15]

    Kapidani, L

    B. Kapidani, L. Codecasa, and J. Sch¨ oberl. An arbitrary-order Cell Method with block-diagonal mass- matrices for the time-dependent 2D Maxwell equations.Journal of Computational Physics, 433:110184, May 2021

    work page 2021

  16. [16]

    Mass lumping the dual cell method to arbitrary polynomial degree for acoustic and electromagnetic waves.J

    Markus Wess, Bernard Kapidani, Lorenzo Codecasa, and Joachim Sch¨ oberl. Mass lumping the dual cell method to arbitrary polynomial degree for acoustic and electromagnetic waves.J. Comput. Phys., 513:Paper No. 113196, 19, 2024

    work page 2024

  17. [17]

    PhD thesis, California Institute of Technology, June 2003

    Anil Nirmal Hirani.Discrete exterior calculus. PhD thesis, California Institute of Technology, June 2003

    work page 2003

  18. [18]

    Hirani, Kaushik Kalyanaraman, and Evan B

    Anil N. Hirani, Kaushik Kalyanaraman, and Evan B. VanderZee. Delaunay Hodge star.Computer-Aided Design, 45(2):540–544, February 2013

    work page 2013

  19. [19]

    Geevers, W

    S. Geevers, W. A. Mulder, and J. J. W. Van Der Vegt. New Higher-Order Mass-Lumped Tetrahedral Elements for Wave Propagation Modelling.SIAM J. Sci. Comput., 40(5):A2830–A2857, January 2018

    work page 2018

  20. [20]

    Egger and B

    H. Egger and B. Radu. A mass-lumped mixed finite element method for acoustic wave propagation. Numer. Math., 145(2):239–269, June 2020

    work page 2020

  21. [21]

    Egger and B

    H. Egger and B. Radu. A Second-Order Finite Element Method with Mass Lumping for Maxwell’s Equations on Tetrahedra.SIAM J. Numer. Anal., 59(2):864–885, January 2021. 20

    work page 2021

  22. [22]

    E. T. Chung, P. Ciarlet, and T.F. Yu. Convergence and superconvergence of staggered discontinu- ous Galerkin methods for the three-dimensional Maxwell’s equations on Cartesian grids.Journal of Computational Physics, 235:14 – 31, 2013

    work page 2013

  23. [23]

    Chung, B

    E. Chung, B. Cockburn, and G. Fu. The Staggered DG Method is the Limit of a Hybridizable DG Method.SIAM J. Numer. Anal., 52(2):915–932, January 2014

    work page 2014

  24. [24]

    Collino and P

    F. Collino and P. B. Monk. Optimizing the perfectly matched layer.Computer Methods in Applied Mechanics and Engineering, 164(1):157–171, October 1998

    work page 1998

  25. [25]

    F. L. Teixeira and W. C. Chew. Complex space approach to perfectly matched layers: a review and some new developments.International Journal of Numerical Modelling: Electronic Networks, Devices and Fields, 13(5):441–455, 2000

    work page 2000

  26. [26]

    Sch¨ oberl

    J. Sch¨ oberl. C++11 implementation of finite elements in ngsolve. Preprint 30/2014, Institute of Analysis and Scientific Computing, TU Wien, 2014

    work page 2014

  27. [27]

    The dual cell method in ngsolve

    Joachim Sch¨ oberl and Markus Wess. The dual cell method in ngsolve. https://ngsolve.github.io/dcm

    work page

  28. [28]

    A Krylov eigenvalue solver based on filtered time domain solutions

    Lothar Nannen and Markus Wess. A Krylov eigenvalue solver based on filtered time domain solutions. Comput. Math. Appl., 176:179–188, 2024. 21

    work page 2024