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
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.
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
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.
Referee Report
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)
- [§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.
- [§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)
- [§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.
- [§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
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
-
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
-
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
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
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.
Reference graph
Works this paper leans on
-
[1]
I. Peshkov and E. Romenski. A hyperbolic model for viscous Newtonian flows.Continuum Mechanics and Thermodynamics, 28:85–104, 2016
work page 2016
-
[2]
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
work page 2016
-
[3]
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
work page 2017
-
[4]
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
work page 1981
-
[5]
S. Klainerman and A. Majda. Compressible and incompressible fluids.Communications in Pure Applied Mathematics, 35:629–651, 1982
work page 1982
-
[6]
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
work page 1965
-
[7]
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
work page 1984
-
[8]
A.J. Chorin. A numerical method for solving incompressible viscous flow problems.Journal of computational physics, 135(2):118–125, 1997
work page 1997
-
[9]
A.J. Chorin. Numerical solution of the Navier-Stokes equations.Mathematics of computation, 22(104):745–762, 1968
work page 1968
-
[10]
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
work page 1983
-
[11]
Patankar.Numerical Heat Transfer and Fluid Flow
V.S. Patankar.Numerical Heat Transfer and Fluid Flow. Hemisphere Publishing Corporation, 1980
work page 1980
-
[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
work page 1986
-
[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
work page 1989
-
[14]
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
work page 1981
-
[15]
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
work page 1973
-
[16]
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
work page 1982
-
[17]
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
work page 1986
-
[18]
M. Fortin. Old and new finite elements for incompressible flows.International Journal for numerical methods in fluids, 1(4):347–364, 1981
work page 1981
-
[19]
R. Verf¨ urth. Finite element approximation on incompressible Navier-Stokes equations with slip boundary condition.Numerische Mathematik, 50:697–721, 1986
work page 1986
-
[20]
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
work page 1982
-
[21]
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
work page 1988
- [22]
-
[23]
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
work page 1959
-
[24]
P.L. Roe. Approximate Riemann solvers, parameter vectors, and difference schemes.Journal of computational physics, 43(2):357–372, 1981. 21
work page 1981
-
[25]
S. Osher and F. Solomon. Upwind difference schemes for hyperbolic systems of conservation laws. Mathematics of computation, 38(158):339–374, 1982
work page 1982
-
[26]
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
work page 1983
-
[27]
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
work page 1991
-
[28]
C.-D. Munz. On Godunov-type schemes for Lagrangian gas dynamics.SIAM Journal on Numerical Analysis, 31(1):17–42, 1994
work page 1994
-
[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
work page 1994
-
[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
work page 2002
-
[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
work page 2013
-
[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
work page 1999
-
[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
work page 2003
-
[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
work page 1995
- [35]
-
[36]
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
work page 2005
-
[37]
E.F. Toro and M.E. V´ azquez-Cend´ on. Flux splitting schemes for the Euler equations.Computers & Fluids, 70:1–12, 2012
work page 2012
-
[38]
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
work page 2016
-
[39]
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
work page 2025
-
[40]
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
work page 2019
-
[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
work page 1966
-
[42]
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
work page 1999
-
[43]
D.S. Balsara. Second-Order Accurate Schemes for Magnetohydrodynamics with Divergence-Free Reconstruction.The Astrophysical Journal Supplement Series, 151:149–184, 2004
work page 2004
-
[44]
D.S. Balsara. Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrody- namic flows.J. Comput. Phys., 229:1970–1993, 2010
work page 1970
-
[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
work page 2014
-
[46]
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
work page 2023
-
[47]
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
work page 2025
-
[48]
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
work page 2024
-
[49]
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
work page 2026
-
[50]
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
work page 2025
-
[51]
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
work page 2024
-
[52]
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
work page 2021
-
[53]
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
work page 1997
-
[54]
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
work page 2000
-
[55]
K. Lipnikov, G. Manzini, and M. Shashkov. Mimetic finite difference method.Journal of Computa- tional Physics, 257:1163–1227, 2014
work page 2014
-
[56]
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
work page 2018
- [57]
-
[58]
J.C. N´ ed´ elec. Mixed finite elements in R3.Numerische Mathematik, 35(3):315–341, 1980
work page 1980
-
[59]
S. Zhang. A new family of stable mixed finite elements for the 3D Stokes equations.Mathematics of computation, 74(250):543–554, 2005
work page 2005
- [60]
-
[61]
D.N. Arnold, R.S. Falk, and R. Winther. Finite element exterior calculus, homological techniques, and applications.Acta numerica, 15:1–155, 2006
work page 2006
-
[62]
Monk.Finite element methods for Maxwell’s equations
P. Monk.Finite element methods for Maxwell’s equations. Oxford university press, 2003
work page 2003
-
[63]
A. Alonso Rodr´ ıguez and A. Valli. Finite element potentials.Applied Numerical Mathematics, 95:2– 14, 2015
work page 2015
-
[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
work page 2016
-
[65]
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
work page 2023
-
[66]
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
work page 2022
- [67]
-
[68]
R. Jeltsch and M. Torrilhon. On curl-preserving finite volume discretizations for shallow water equations.BIT Numerical Mathematics, 46(Suppl 1):35–53, 2006
work page 2006
-
[69]
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
work page 2023
-
[70]
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
work page 2021
-
[71]
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
work page 2023
-
[72]
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
work page 2024
-
[73]
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
work page 2023
-
[74]
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
work page 2025
-
[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
work page 2021
-
[76]
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
work page 2026
-
[77]
S.K. Godunov. Symmetric form of the equations of magnetohydrodynamics.Numerical Methods for Mechanics of Continuous Media, 3(1):26–31, 1972
work page 1972
-
[78]
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
work page 1972
- [79]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.