pith. sign in

arxiv: 2604.21906 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

A structure-preserving semi-implicit finite volume scheme on vertex-staggered unstructured meshes

Pith reviewed 2026-05-09 20:34 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords structure-preservingsemi-implicit finite volumevertex-staggered meshunstructured gridsdivergence-freeincompressible flowsMHDGPR model
0
0 comments X

The pith

A semi-implicit finite volume scheme on vertex-staggered unstructured meshes preserves divergence-free and curl-free vector fields exactly.

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

The paper develops a structure-preserving semi-implicit finite volume method on meshes formed by Delaunay triangles and their dual polygons. Scalar quantities sit at the nodes and are advanced implicitly via a continuous finite-element pressure solve, while vector quantities sit at triangle barycenters and advance explicitly via a compatible finite-volume update. The staggered placement produces exact preservation of divergence-free velocity in the incompressible limit, divergence-free magnetic fields in MHD, and curl-free inverse deformation gradients in the GPR model, all without auxiliary corrections. Because the time step is no longer limited by the sound speed, the method remains efficient for low-Mach and incompressible regimes on general unstructured grids.

Core claim

By storing scalars at mesh nodes and vectors at triangle barycenters on a primal Delaunay triangulation paired with its dual star polygons, the semi-implicit discretization exactly preserves the divergence-free and curl-free properties of the vector fields for the incompressible Euler and Navier-Stokes equations, incompressible MHD, and the incompressible GPR model for continuum mechanics.

What carries the argument

Compatible vertex-staggered discretization on primal Delaunay triangles and dual polygons, with scalars at vertices and vectors at barycenters, that enforces exact structure preservation through the finite-volume and finite-element updates.

Load-bearing premise

The vertex-staggered placement of scalars at nodes and vectors at barycenters remains compatible and structure-preserving for all listed systems on arbitrary unstructured meshes without additional corrections.

What would settle it

A long-time incompressible MHD simulation on a highly distorted unstructured mesh in which the discrete divergence of the magnetic field grows above machine round-off error would falsify the exact-preservation claim.

Figures

Figures reproduced from arXiv: 2604.21906 by Elena Bernardelli, Elena Gaburro, Michael Dumbser.

Figure 1
Figure 1. Figure 1: On the left, a control volume ωc, with its barycenter and edge midpoints. In violet, the sub-cell ωpc which results by connecting the vertex with position xp, the cell barycenter of coordinates xbc = 1/(d+ 1)P p∈P(c) xp and the left and right midpoints of the edges impinging on node p, namely p −1/2 and p +1/2 . In the middle, the dual cell ωp (cyan) is obtained by the union of all sub-cells ωc sharing a g… view at source ↗
Figure 2
Figure 2. Figure 2: Numerical solution of the Taylor-Green vortex at the final time [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Numerical solution of the Taylor-Green vortex at the final time [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Numerical solution of the Sod shock tube problem obtained using our structure-preserving semi [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical solution of the Gresho vortex at the final time [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Numerical solution of the MHD vortex ad a final time [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Numerical solution of the Solid rotor problem at the final time [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
read the original abstract

We present a novel structure-preserving semi-implicit finite volume method on vertex-based staggered meshes for the compatible discretization of first order systems of time-dependent partial differential equations (PDEs). The method preserves divergence-free and curl-free vector fields exactly thanks to the compatible vertex-staggered discretization of the state variables on unstructured grids that are constituted by primal Delaunay triangles and their dual polygons. For the weakly compressible Euler equations, the scheme is asymptotic preserving, yielding a consistent discretization of the incompressible limit as the Mach number goes to zero. The new scheme applies to a broad spectrum of PDEs, including the weakly compressible and incompressible Euler and Navier-Stokes equations, the incompressible magnetohydrodynamics (MHD) system, and the incompressible version of the first-order hyperbolic Godunov-Peshkov-Romenski (GPR) model for continuum mechanics. The computational domain is covered by a primal triangular mesh and a dual tessellation made of so-called star polygons. Scalar quantities (pressure, density, viscous stress) are defined at nodes, with pressure updated implicitly in a continuous finite element fashion, yielding a symmetric and positive definite pressure system. Instead, vector fields (velocity, momentum, magnetic and distortion fields) are stored at triangle barycenters and evolved explicitly using a compatible finite volume scheme. Thanks to the semi-implicit discretization, the CFL condition is independent of the sound speed, allowing simulations at low Mach numbers. The fully compatible formulation ensures exactly divergence-free velocity field in the incompressible limit, exactly divergence-free magnetic field for MHD, and exactly curl-free inverse deformation gradient in solid mechanics. The method is validated through a wide set of test cases.

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

2 major / 2 minor

Summary. The paper presents a novel structure-preserving semi-implicit finite volume method on vertex-staggered unstructured meshes for first-order time-dependent PDE systems. Scalars (pressure, density) are placed at mesh nodes and updated implicitly via continuous finite elements, while vector fields (velocity, momentum, magnetic field, distortion) are stored at triangle barycenters and evolved explicitly via finite volumes on the dual star polygons. The central claims are exact preservation of divergence-free velocity (incompressible limit), divergence-free magnetic field (MHD), and curl-free inverse deformation gradient (GPR), together with asymptotic consistency to the incompressible Euler/Navier-Stokes equations as Mach number tends to zero and a CFL condition independent of sound speed.

Significance. If the claimed exact structure preservation is rigorously established, the scheme would constitute a useful advance for low-Mach and MHD computations on unstructured meshes, eliminating the need for artificial compressibility or separate projection steps while retaining a symmetric positive-definite pressure system. The vertex-staggered Delaunay-Voronoi construction is a natural choice for compatibility, and the breadth of target systems (Euler, NS, MHD, GPR) is attractive.

major comments (2)
  1. [§2 (discretization) and §3 (properties)] The exact preservation of div u = 0 (or div B = 0) after each time step is asserted to follow from the adjointness of the continuous FE gradient operator (nodal scalars) and the FV divergence operator (barycentric vectors). No discrete integration-by-parts identity, summation-by-parts formula, or explicit verification that the implicit pressure solve projects exactly onto the kernel of the discrete divergence is supplied for arbitrary Delaunay triangulations. This identity is load-bearing for the central claims in the abstract and §1; without it the “exactly divergence-free” property remains an unproven assertion rather than a demonstrated consequence of the discretization.
  2. [§4 (asymptotic analysis) and §5 (numerics)] The asymptotic-preserving property for the weakly compressible Euler equations as M → 0 is stated but not accompanied by a formal consistency analysis or uniform-in-M error estimates. The numerical examples in §5 show good behavior at low Mach, yet the absence of a rigorous limit analysis leaves open whether the scheme converges to a consistent discretization of the incompressible system on general unstructured meshes.
minor comments (2)
  1. [§2.1] Notation for the dual star polygons and the precise definition of the inner products used for the adjoint relation should be introduced earlier and used consistently throughout the operator definitions.
  2. [§5] Several test cases are mentioned in the abstract and §5, but the manuscript would benefit from a compact table summarizing the observed orders of accuracy and the measured divergence errors for each model.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below and will revise the paper to strengthen the theoretical foundations as suggested.

read point-by-point responses
  1. Referee: [§2 (discretization) and §3 (properties)] The exact preservation of div u = 0 (or div B = 0) after each time step is asserted to follow from the adjointness of the continuous FE gradient operator (nodal scalars) and the FV divergence operator (barycentric vectors). No discrete integration-by-parts identity, summation-by-parts formula, or explicit verification that the implicit pressure solve projects exactly onto the kernel of the discrete divergence is supplied for arbitrary Delaunay triangulations. This identity is load-bearing for the central claims in the abstract and §1; without it the “exactly divergence-free” property remains an unproven assertion rather than a demonstrated consequence of the discretization.

    Authors: We agree with the referee that an explicit discrete integration-by-parts identity is required to rigorously establish the exact preservation property on arbitrary Delaunay triangulations. The manuscript derives the property from the adjointness of the continuous finite-element gradient (nodal) and finite-volume divergence (barycentric) operators that arises naturally from the vertex-staggered Delaunay-Voronoi construction. However, we acknowledge that a self-contained summation-by-parts formula and direct verification of the projection onto the discrete divergence kernel were not supplied. In the revised manuscript we will add a dedicated paragraph in §3 that states and proves the relevant discrete integration-by-parts identity for general Delaunay triangulations, thereby confirming that the implicit pressure solve exactly enforces the divergence-free constraint. revision: yes

  2. Referee: [§4 (asymptotic analysis) and §5 (numerics)] The asymptotic-preserving property for the weakly compressible Euler equations as M → 0 is stated but not accompanied by a formal consistency analysis or uniform-in-M error estimates. The numerical examples in §5 show good behavior at low Mach, yet the absence of a rigorous limit analysis leaves open whether the scheme converges to a consistent discretization of the incompressible system on general unstructured meshes.

    Authors: We thank the referee for highlighting the need for a formal consistency analysis. The manuscript presents the semi-implicit scheme and demonstrates its low-Mach behavior through numerical experiments in §5, but we agree that a rigorous limit analysis as M → 0 was not included. In the revised version we will expand §4 to contain a consistency analysis showing that the scheme reduces to a consistent discretization of the incompressible Euler equations on unstructured meshes in the zero-Mach limit. While uniform-in-M error estimates lie beyond the scope of the present work, this addition will provide the missing theoretical justification for the asymptotic-preserving claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity; exact preservation follows from the staggered discretization choice rather than self-definition or fitted inputs.

full rationale

The paper presents a new semi-implicit finite-volume scheme whose structure-preserving properties (exact div-free and curl-free fields) are claimed to arise directly from the vertex-staggered placement of variables on Delaunay-Voronoi meshes, with nodal scalars updated via continuous FE and barycentric vectors via explicit FV. This is a constructive discretization choice, not a reduction of the output to the input by definition, fitting, or self-citation chain. No equations or claims in the provided text equate the preservation result to a parameter fit or rename a known result; the central claim remains an independent consequence of the mesh and operator compatibility that the authors must demonstrate separately. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the central claim rests on the assumption that a compatible discretization exists between primal Delaunay triangles and dual star polygons for the listed PDE systems. No free parameters are mentioned. No new physical entities are introduced.

axioms (1)
  • domain assumption The primal triangular mesh and its dual polygonal tessellation admit a compatible staggering of scalar and vector variables that exactly enforces divergence and curl constraints by construction.
    Invoked in the description of the vertex-staggered discretization and the exact preservation property.

pith-pipeline@v0.9.0 · 5599 in / 1313 out tokens · 23249 ms · 2026-05-09T20:34:00.649219+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

93 extracted references · 93 canonical work pages

  1. [1]

    Peshkov and E

    I. Peshkov and E. Romenski. A hyperbolic model for viscous Newtonian flows.Continuum Mechanics and Thermodynamics, 28:85–104, 2016

  2. [2]

    Dumbser, I

    M. Dumbser, I. Peshkov, E. Romenski, and O. Zanotti. High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat–conducting fluids and elastic solids.Journal of Computational Physics, 314:824–862, 2016

  3. [3]

    Dumbser, I

    M. Dumbser, I. Peshkov, E. Romenski, and O. Zanotti. High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro–dynamics. Journal of Computational Physics, 348:298–342, 2017

  4. [4]

    Klainerman and A

    S. Klainerman and A. Majda. Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit of compressible fluids.Communications on pure and applied Mathemat- ics, 34(4):481–524, 1981

  5. [5]

    Klainerman and A

    S. Klainerman and A. Majda. Compressible and incompressible fluids.Communications in Pure Applied Mathematics, 35:629–651, 1982

  6. [6]

    Harlow, J.E

    F.H. Harlow, J.E. Welch, et al. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface.Physics of fluids, 8(12):2182, 1965

  7. [7]

    Casulli and D

    V. Casulli and D. Greenspan. Pressure method for the numerical solution of transient, compressible fluid flows.International Journal for Numerical Methods in Fluids, 4(11):1001–1012, 1984. 20

  8. [8]

    A.J. Chorin. A numerical method for solving incompressible viscous flow problems.Journal of computational physics, 135(2):118–125, 1997

  9. [9]

    A.J. Chorin. Numerical solution of the Navier-Stokes equations.Mathematics of computation, 22(104):745–762, 1968

  10. [10]

    Patankar and D.B

    S.V. Patankar and D.B. Spalding. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. InNumerical prediction of flow, heat transfer, turbulence and combustion, pages 54–73. Elsevier, 1983

  11. [11]

    Patankar.Numerical Heat Transfer and Fluid Flow

    V.S. Patankar.Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing Corporation, 1980

  12. [12]

    J. Van Kan. A second-order accurate pressure-correction scheme for viscous incompressible flow. SIAM journal on scientific and statistical computing, 7(3):870–891, 1986

  13. [13]

    J.B. Bell, P. Colella, and H.M. Glaz. A second-order projection method for the incompressible Navier-Stokes equations.Journal of computational physics, 85(2):257–283, 1989

  14. [14]

    Hirt and B.D

    C.W. Hirt and B.D. Nichols. Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of computational physics, 39(1):201–225, 1981

  15. [15]

    Taylor and P

    C. Taylor and P. Hood. A numerical solution of the Navier-Stokes equations using the finite element technique.Computers & Fluids, 1(1):73–100, 1973

  16. [16]

    Brooks and T

    A.N. Brooks and T. Jr Hughes. Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations.Computer methods in applied mechanics and engineering, 32(1-3):199–259, 1982

  17. [17]

    Jr Hughes, M

    T. Jr Hughes, M. Mallet, and M. Akira. A new finite element formulation for computational fluid dynamics: II. Beyond SUPG.Computer methods in applied mechanics and engineering, 54(3):341– 355, 1986

  18. [18]

    M. Fortin. Old and new finite elements for incompressible flows.International Journal for numerical methods in fluids, 1(4):347–364, 1981

  19. [19]

    Verf¨ urth

    R. Verf¨ urth. Finite element approximation on incompressible Navier-Stokes equations with slip boundary condition.Numerische Mathematik, 50:697–721, 1986

  20. [20]

    Heywood and R

    J.G. Heywood and R. Rannacher. Finite element approximation of the nonstationary Navier-Stokes problem. I. Regularity of solutions and second-order error estimates for spatial discretization.SIAM Journal on Numerical Analysis, 19(2):275–311, 1982

  21. [21]

    Heywood and R

    J.G. Heywood and R. Rannacher. Finite element approximation of the nonstationary Navier-Stokes problem III. Smoothing property and higher order error estimates for spatial discretization.SIAM Journal on Numerical Analysis, 25(3):489–512, 1988

  22. [22]

    Lax and B

    P.D. Lax and B. Wendroff. Systems of conservation laws. InSelected Papers Volume I, pages 263–283. Springer, 2005

  23. [23]

    Finite difference methods for the computation of discontinuous solutions of the equa- tions of fluid dynamics.Mat

    S.K Godunov. Finite difference methods for the computation of discontinuous solutions of the equa- tions of fluid dynamics.Mat. Sb., 47:271–306, 1959

  24. [24]

    P.L. Roe. Approximate Riemann solvers, parameter vectors, and difference schemes.Journal of computational physics, 43(2):357–372, 1981. 21

  25. [25]

    Osher and F

    S. Osher and F. Solomon. Upwind difference schemes for hyperbolic systems of conservation laws. Mathematics of computation, 38(158):339–374, 1982

  26. [26]

    Harten, P.D

    A. Harten, P.D. Lax, and B. Leer. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws.SIAM review, 25(1):35–61, 1983

  27. [27]

    Einfeldt, C.-D

    B. Einfeldt, C.-D. Munz, P.L. Roe, and B. Sj¨ ogreen. On Godunov-type methods near low densities. Journal of computational physics, 92(2):273–295, 1991

  28. [28]

    C.-D. Munz. On Godunov-type schemes for Lagrangian gas dynamics.SIAM Journal on Numerical Analysis, 31(1):17–42, 1994

  29. [29]

    E.F. Toro, M. Spruce, and W. Speares. Restoration of the contact surface in the HLL-Riemann solver.Shock waves, 4(1):25–34, 1994

  30. [30]

    LeVeque.Finite volume methods for hyperbolic problems, volume 31

    R.J. LeVeque.Finite volume methods for hyperbolic problems, volume 31. Cambridge university press, 2002

  31. [31]

    Toro.Riemann solvers and numerical methods for fluid dynamics: a practical introduction

    E.F. Toro.Riemann solvers and numerical methods for fluid dynamics: a practical introduction. Springer Science & Business Media, 2013

  32. [32]

    A. Meister. Asymptotic single and multiple scale expansions in the low Mach number limit.SIAM Journal on Applied Mathematics, 60(1):256–271, 1999

  33. [33]

    C.-D. Munz, S. Roller, R. Klein, and K.J. Geratz. The extension of incompressible flow solvers to the weakly compressible regime.Computers & Fluids, 32(2):173–196, 2003

  34. [34]

    R. Klein. Semi-implicit extension of a Godunov-type scheme based on low Mach number asymptotics I: One-dimensional flow.Journal of Computational Physics, 121(2):213–237, 1995

  35. [35]

    Klein, N

    R. Klein, N. Botta, T. Schneider, C.-D. Munz, S. Roller, A. Meister, L. Hoffmann, and T. Sonar. Asymptotic adaptive methods for multi-scale problems in fluid mechanics.Journal of Engineering Mathematics, 39:261–343, 2001

  36. [36]

    Park and C.-D

    J.H. Park and C.-D. Munz. Multiple pressure variables methods for fluid flow at all Mach numbers. International journal for numerical methods in fluids, 49(8):905–931, 2005

  37. [37]

    Toro and M.E

    E.F. Toro and M.E. V´ azquez-Cend´ on. Flux splitting schemes for the Euler equations.Computers & Fluids, 70:1–12, 2012

  38. [38]

    Dumbser and V

    M. Dumbser and V. Casulli. A conservative, weakly nonlinear semi-implicit finite volume scheme for the compressible Navier- Stokes equations with general equation of state.Applied Mathematics and Computation, 272:479–497, 2016

  39. [39]

    Zampa and M

    E. Zampa and M. Dumbser. An asymptotic-preserving and exactly mass-conservative semi-implicit scheme for weakly compressible flows based on compatible finite elements.Journal of Computational Physics, 521:113551, 2025

  40. [40]

    Dumbser, D.S

    M. Dumbser, D.S. Balsara, M. Tavelli, and F. Fambri. A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics.Int. J. Numer. Meth. Fl., 89:16–42, 2019

  41. [41]

    K. 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, 1966. 22

  42. [42]

    Balsara and D.S

    D.S. Balsara and D.S. Spicer. A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations.Journal of Computational Physics, 149:270–292, 1999

  43. [43]

    D.S. Balsara. Second-Order Accurate Schemes for Magnetohydrodynamics with Divergence-Free Reconstruction.The Astrophysical Journal Supplement Series, 151:149–184, 2004

  44. [44]

    D.S. Balsara. Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrody- namic flows.J. Comput. Phys., 229:1970–1993, 2010

  45. [45]

    D.S. Balsara. Multidimensional Riemann problem with self-similar internal structure – Part I – Application to hyperbolic conservation laws on structured meshes.J. Comput. Phys., 277:163–200, 2014

  46. [46]

    Abgrall, S

    R. Abgrall, S. Busto, and M. Dumbser. A simple and general framework for the construction of thermodynamically compatible schemes for computational fluid and solid mechanics.Applied Math- ematics and Computation, 440:127629, 2023

  47. [47]

    Abgrall, M

    R. Abgrall, M. Dumbser, and P.-H. Maire. A simple and general framework for the construction of exactly div-curl-grad compatible discontinuous galerkin finite element schemes on unstructured simplex meshes.Journal of Computational Physics, 541:114340, 2025

  48. [48]

    Boscheri and A

    W. Boscheri and A. Thomann. A structure-preserving semi-implicit IMEX finite volume scheme for ideal magnetohydrodynamics at all Mach and Alfv´ en numbers.Journal of Scientific Computing, 100(3):67, 2024

  49. [49]

    Boscheri, M

    W. Boscheri, M. Dumbser, R. Loub` ere, and P.-H. Maire. A structure-preserving and thermodynam- ically compatible cell-centered Lagrangian finite volume scheme for continuum mechanics.SIAM Journal on Numerical Analysis, 2026. in press

  50. [50]

    R´ ıo-Mart´ ın, F

    L. R´ ıo-Mart´ ın, F. Dhaouadi, and M. Dumbser. An exactly curl-free finite-volume/finite-difference scheme for a hyperbolic compressible isentropic two-phase model.Journal of Scientific Computing, 102(1):13, 2025

  51. [51]

    Boscheri, M

    W. Boscheri, M. Dumbser, and P.-H. Maire. A new thermodynamically compatible finite volume scheme for lagrangian gas dynamics.SIAM Journal on Scientific Computing, 46(4):A2224–A2247, 2024

  52. [52]

    Peshkov, M

    I. Peshkov, M. Dumbser, W. Boscheri, E. Romenski, S. Chiocchetti, and M. Ioriatti. Simulation of non-Newtonian viscoplastic flows with a unified first order hyperbolic model and a structure- preserving semi-implicit scheme.Computers & Fluids, 224:104963, 2021

  53. [53]

    Hyman and M

    J.M. Hyman and M. Shashkov. Natural discretizations for the divergence, gradient, and curl on logically rectangular grids.Computers & Mathematics with Applications, 33(4):81–104, 1997

  54. [54]

    Margolin, M

    L.G. Margolin, M. Shashkov, and P.K. Smolarkiewicz. A discrete operator calculus for finite difference approximations.Computer methods in applied mechanics and engineering, 187(3-4):365–383, 2000

  55. [55]

    Lipnikov, G

    K. Lipnikov, G. Manzini, and M. Shashkov. Mimetic finite difference method.Journal of Computa- tional Physics, 257:1163–1227, 2014

  56. [56]

    Balsara, S

    D.S. Balsara, S. Garain, A. Taflove, and G. Montecinos. Computational electrodynamics in material media with constraint-preservation, multidimensional Riemann solvers and sub-cell resolution–Part II, higher order FVTD schemes.Journal of Computational Physics, 354:613–645, 2018. 23

  57. [57]

    Hazra, P

    A. Hazra, P. Chandrashekar, and D.S. Balsara. Globally constraint-preserving FR/DG scheme for Maxwell’s equations at all orders.Journal of computational physics, 394:298–328, 2019

  58. [58]

    N´ ed´ elec

    J.C. N´ ed´ elec. Mixed finite elements in R3.Numerische Mathematik, 35(3):315–341, 1980

  59. [59]

    S. Zhang. A new family of stable mixed finite elements for the 3D Stokes equations.Mathematics of computation, 74(250):543–554, 2005

  60. [60]

    Hiptmair

    R. Hiptmair. Finite elements in computational electromagnetism.Acta Numerica, 11:237–339, 2002

  61. [61]

    Arnold, R.S

    D.N. Arnold, R.S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications.Acta numerica, 15:1–155, 2006

  62. [62]

    Monk.Finite element methods for Maxwell’s equations

    P. Monk.Finite element methods for Maxwell’s equations. Oxford university press, 2003

  63. [63]

    Alonso Rodr´ ıguez and A

    A. Alonso Rodr´ ıguez and A. Valli. Finite element potentials.Applied Numerical Mathematics, 95:2– 14, 2015

  64. [64]

    P. M. Campos and E. Sonnendr¨ ucker. Gauss-compatible Galerkin schemes for time-dependent Maxwell equations.Mathematics of Computation, 85(302):2651–2685, 2016

  65. [65]

    Di Pietro and J

    D.A. Di Pietro and J. Droniou. An arbitrary-order discrete de Rham complex on polyhedral meshes: Exactness, Poincar´ e inequalities, and consistency.Foundations of Computational Mathe- matics, 23(1):85–164, 2023

  66. [66]

    Bruno and E

    L.B. Bruno and E. Zampa. Unisolvent and minimal physical degrees of freedom for the second family of polynomial differential forms.ESAIM: Mathematical Modelling and Numerical Analysis, 56(6):2239–2253, 2022

  67. [67]

    Zampa, A

    E. Zampa, A. Alonso Rodr´ ıguez, and F. Rapetti. Using the FES framework to derive new phys- ical degrees of freedom for finite element spaces of differential forms.Advances in Computational Mathematics, 49(2):17, 2023

  68. [68]

    Jeltsch and M

    R. Jeltsch and M. Torrilhon. On curl-preserving finite volume discretizations for shallow water equations.BIT Numerical Mathematics, 46(Suppl 1):35–53, 2006

  69. [69]

    Balsara, R

    D.S. Balsara, R. K¨ appeli, W. Boscheri, and M. Dumbser. Curl constraint-preserving reconstruction and the guidance it gives for mimetic scheme design.Communications on Applied Mathematics and Computation, 5(1):235–294, 2023

  70. [70]

    Boscheri, M

    W. Boscheri, M. Dumbser, M. Ioriatti, I. Peshkov, and E. Romenski. A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics.Journal of Computational Physics, 424:109866, 2021

  71. [71]

    Chiocchetti and M

    S. Chiocchetti and M. Dumbser. An exactly curl-free staggered semi-implicit finite volume scheme for a first order hyperbolic model of viscous two-phase flows with surface tension.Journal of Scientific Computing, 94(1):24, 2023

  72. [72]

    Jung and V

    J. Jung and V. Perrier. A curl preserving finite volume scheme by space velocity enrichment. Appli- cation to the low Mach number accuracy problem.Journal of Computational Physics, 515:113252, 2024

  73. [73]

    Dhaouadi and M

    F. Dhaouadi and M. Dumbser. A structure-preserving finite volume scheme for a hyperbolic refor- mulation of the Navier–Stokes–Korteweg equations.Mathematics, 11(4):876, 2023. 24

  74. [74]

    Barsukow, R

    W. Barsukow, R. Loub` ere, and P.-H. Maire. A node-conservative vorticity preserving finite volume method for linear acoustics on unstructured grids.Mathematics of Computation, 94:2299–2343, 2025

  75. [75]

    F. Fambri. A novel structure preserving semi-implicit finite volume method for viscous and resistive magnetohydrodynamics.International Journal for Numerical Methods in Fluids, 93(12):3447–3489, 2021

  76. [76]

    Dumbser, A

    M. Dumbser, A. Thomann, M. Tavelli, and W. Boscheri. A structure-preserving semi-implicit four- split scheme for continuum mechanics.Journal of Computational Physics, 559:114889, 2026

  77. [77]

    S.K. Godunov. Symmetric form of the equations of magnetohydrodynamics.Numerical Methods for Mechanics of Continuous Media, 3(1):26–31, 1972

  78. [78]

    Godunov and E.I

    S.K. Godunov and E.I. Romenski. Nonstationary equations of the nonlinear theory of elasticity in Euler coordinates.J. Appl. Mech. Tech. Phys., 13:868–885, 1972

  79. [79]

    Romenski

    E. Romenski. Hyperbolic systems of thermodynamically compatible conservation laws in continuum mechanics.Mathematical and computer modelling, 28(10):115–130, 1998

  80. [80]

    Maire, R

    P.-H. Maire, R. Abgrall, J. Breil, and J. Ovadia. A cell-centered lagrangian scheme for two- dimensional compressible flow problems.SIAM J. Sci. Comput., 29:1781–1824, 2007

Showing first 80 references.