pith. sign in

arxiv: 2605.13048 · v3 · pith:HRA4DP66new · submitted 2026-05-13 · 🧮 math.AP · cs.NA· math.NA

Exact conservation and the Onsager threshold: a discrete exterior calculus theory for incompressible Navier-Stokes Equations

Pith reviewed 2026-05-22 10:38 UTC · model grok-4.3

classification 🧮 math.AP cs.NAmath.NA
keywords discrete exterior calculusincompressible Navier-Stokesexact conservationOnsager thresholdweak solutionsmeasure-valued solutionsstructure-preserving discretizationEuler equations
0
0 comments X

The pith

Exact algebraic conservation at the discrete level excludes entire classes of energy-dissipating weak solutions to the Euler equations.

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

The paper constructs a discretization of the incompressible Navier-Stokes and Euler equations on prismatic meshes that enforces exact algebraic conservation of energy and other invariants. This conservation property functions as a selection principle: it rules out weak solutions that standard discretizations can reach without restriction. The argument is carried through four regimes, from smooth solutions with explicit convergence rates to the inviscid limit where dissipative Euler solutions are structurally excluded. A sympathetic reader sees this as evidence that structure preservation can eliminate unphysical behaviors that arise in the absence of exact conservation. The work isolates the interval 1/3 < alpha < 1 as the remaining open region where conservation holds but uniqueness is unresolved.

Core claim

Exact algebraic conservation of the discrete energy prevents any subsequence of solutions from converging to an energy-dissipating weak solution of the Euler equations at any regularity. In the inviscid measure-valued regime the same conservation implies that limits are conservative measure-valued solutions whose concentration defect vanishes for Holder exponents alpha greater than 1/3 whenever a uniform C^{0,alpha} bound is available on the discrete solutions.

What carries the argument

Exact algebraic conservation of energy and circulation arising from the discrete exterior calculus formulation on prismatic Delaunay-Voronoi meshes.

Load-bearing premise

Discrete solutions must admit a uniform Holder bound C to the alpha in the inviscid regime so that the concentration defect can be shown to vanish above the Onsager threshold.

What would settle it

A concrete sequence of discrete solutions whose energy remains exactly conserved yet converges in L2 to a weak Euler solution that dissipates positive energy would falsify the exclusion result.

Figures

Figures reproduced from arXiv: 2605.13048 by Peter Korn.

Figure 1
Figure 1. Figure 1: Regularity landscape and four convergence modes of the DEC scheme (Thm. [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
read the original abstract

We develop a rigorous theory for a structure-preserving discretisation of the incompressible Euler and Navier--Stokes equations, based on discrete exterior calculus on prismatic Delaunay--Voronoi meshes over closed Riemannian manifolds. The central result is a selection principle: exact algebraic conservation at the discrete level is not merely a fidelity property but rules out entire classes of weak solutions that other discretisations reach unconditionally. We establish this in four regimes. \emph{Smooth solutions}: convergence at rate $\mathcal{O}(h^{\min(r_{\rm rec},\,r_\star)}\,|\log h|^{\beta_d})$, uniformly in viscosity $\nu \ge 0$, with $\beta_3 = 0$ and $\beta_2 = 1$; first order on general meshes and second order on meshes with centroid proximity and reconstruction symmetry. \emph{Leray--Hopf weak regime}: subsequential $L^2$ limits are weak solutions of the viscous system. \emph{Inviscid measure-valued regime}: limits are conservative measure-valued Euler solutions; their concentration defect vanishes above the Onsager threshold $\alpha > 1/3$ \emph{provided the discrete solutions admit a uniform $C^{0,\alpha}$ bound there}. \emph{Dissipative regime}: no subsequence converges to an energy-dissipating Euler solution at any regularity, a structural exclusion that follows from exact discrete energy conservation and distinguishes the scheme. The gap $1/3 < \alpha < 1$, where energy conservation and defect-free convergence hold but uniqueness remains open, isolates the central open problem of inviscid fluid dynamics.

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

1 major / 2 minor

Summary. The manuscript develops a structure-preserving discretization of the incompressible Euler and Navier-Stokes equations via discrete exterior calculus on prismatic Delaunay-Voronoi meshes over closed Riemannian manifolds. It establishes convergence rates for smooth solutions uniform in viscosity, subsequential L2 convergence to Leray-Hopf weak solutions, and a selection principle: exact algebraic discrete conservation implies conservative measure-valued Euler limits whose concentration defect vanishes for Hölder exponents α > 1/3 (under a uniform C^{0,α} bound hypothesis on the discrete solutions) and, unconditionally, excludes energy-dissipating Euler solutions in the dissipative regime.

Significance. If the uniform Hölder bound can be established independently, the exact algebraic conservation would furnish a rigorous selection principle that rules out classes of weak solutions reached by other schemes, providing new insight into the Onsager threshold and the structure of inviscid limits. The algebraic (parameter-free) nature of the conservation and the regime-specific results are clear strengths.

major comments (1)
  1. [Abstract / Inviscid measure-valued regime] Abstract and inviscid measure-valued regime: the claim that concentration defect vanishes for α > 1/3 (yielding conservative measure-valued Euler solutions) is conditional on the discrete solutions admitting a uniform C^{0,α} bound. This bound is invoked as a hypothesis rather than derived from the discrete exterior calculus structure or the exact algebraic conservation; without an independent proof or uniform-in-h, ν=0 estimate on the Hölder seminorm, the structural exclusion of dissipative solutions does not follow unconditionally from conservation alone.
minor comments (2)
  1. [Smooth solutions] The convergence rate statement includes |log h|^{β_d} with β_3 = 0 and β_2 = 1; explicit values or a general expression for β_d in arbitrary dimension would improve clarity.
  2. [Mesh and discretization setup] Notation for the reconstruction and proximity conditions (centroid proximity, reconstruction symmetry) could be cross-referenced to the mesh construction section for readers unfamiliar with the discrete exterior calculus setting.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting this important distinction in the inviscid regime. We address the comment below and will revise the presentation to improve clarity.

read point-by-point responses
  1. Referee: Abstract / Inviscid measure-valued regime: the claim that concentration defect vanishes for α > 1/3 (yielding conservative measure-valued Euler solutions) is conditional on the discrete solutions admitting a uniform C^{0,α} bound. This bound is invoked as a hypothesis rather than derived from the discrete exterior calculus structure or the exact algebraic conservation; without an independent proof or uniform-in-h, ν=0 estimate on the Hölder seminorm, the structural exclusion of dissipative solutions does not follow unconditionally from conservation alone.

    Authors: We agree that the vanishing of the concentration defect for α > 1/3 is conditional on a uniform C^{0,α} bound on the discrete solutions; this is explicitly stated as a hypothesis in the abstract and the relevant theorem. Deriving such a bound independently from the discrete exterior calculus or from exact algebraic conservation alone is an open question and is not claimed in the manuscript. The exact conservation is nevertheless essential: it guarantees that the limiting objects are conservative measure-valued solutions, after which the additional Hölder hypothesis implies vanishing defect. In contrast, the unconditional exclusion of energy-dissipating Euler solutions (in the dissipative regime) follows directly from exact discrete energy conservation without any Hölder assumption. We will revise the abstract and the discussion of the inviscid measure-valued regime to more sharply separate the conditional defect-vanishing result from the unconditional dissipative exclusion, thereby addressing the concern about over-stating what follows from conservation alone. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central claims follow directly from discrete algebraic conservation without reduction to inputs

full rationale

The paper's selection principle rests on exact algebraic conservation as an independent discrete property that structurally excludes energy-dissipating Euler solutions in the dissipative regime. The inviscid measure-valued result is explicitly conditional on an assumed uniform C^{0,α} bound rather than derived from the conservation itself or any self-referential loop. No self-definitional equivalences, fitted inputs renamed as predictions, or load-bearing self-citations appear in the derivation chain; the smooth-regime convergence rates and Leray-Hopf limits are presented as standard consequences of the structure-preserving discretization. The overall argument remains self-contained against external benchmarks of discrete conservation and weak-solution theory.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The theory rests on the standard framework of discrete exterior calculus and geometric properties of prismatic Delaunay-Voronoi meshes over closed Riemannian manifolds; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Exact algebraic conservation holds for the discrete exterior calculus operators on the specified prismatic Delaunay-Voronoi meshes
    Invoked as the foundation for the selection principle across all regimes.
  • domain assumption Mesh properties including centroid proximity and reconstruction symmetry enable the stated convergence rates
    Referenced for first- and second-order convergence on general versus special meshes.

pith-pipeline@v0.9.0 · 5827 in / 1446 out tokens · 48732 ms · 2026-05-22T10:38:31.489486+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    Finiteelementexteriorcalculus, homological techniques, and applications.Acta Numer., 15:1–155, 2006

    D.N.Arnold, R.S.Falk, andR.Winther. Finiteelementexteriorcalculus, homological techniques, and applications.Acta Numer., 15:1–155, 2006

  2. [2]

    D. N. Arnold, R. S. Falk, and R. Winther. Finite element exterior calculus: from Hodge theory to numerical stability.Bull. Amer. Math. Soc. (N.S.), 47(2):281–354, 2010. 72

  3. [3]

    V. I. Arnold. Sur la géométrie différentielle des groupes de Lie de dimension infinie et ses applications à l’hydrodynamique des fluides parfaits.Annales de l’Institut Fourier, 16(1):319–361, 1966

  4. [4]

    L. C. Berselli, T. Iliescu, and W. J. Layton.Mathematics of Large Eddy Simulation of Turbulent Flows. Scientific Computation. Springer, Berlin, 2006

  5. [5]

    Bossavit

    A. Bossavit. Whitney forms: a class of finite elements for three-dimensional compu- tations in electromagnetism.IEE Proc. A, 135(8):493–500, 1988

  6. [6]

    Bossavit

    A. Bossavit. Computational electromagnetism: Variational formulations, complemen- tarity, edge elements.Academic Press, 1998. Monograph

  7. [7]

    Bossavit

    A. Bossavit. Generalized finite differences in computational electromagnetics.Progress in Electromagnetics Research, 32:45–64, 2001

  8. [8]

    Brenier, C

    Y. Brenier, C. De Lellis, and L. Székelyhidi Jr. Weak-strong uniqueness for measure- valued solutions.Comm. Math. Phys., 305(2):351–361, 2011

  9. [9]

    Brezzi, K

    F. Brezzi, K. Lipnikov, and M. Shashkov. Convergence of the mimetic finite differ- ence method for diffusion problems on polyhedral meshes.SIAM J. Numer. Anal., 43(5):1872–1896, 2005

  10. [10]

    Buckmaster, C

    T. Buckmaster, C. De Lellis, L. Székelyhidi Jr., and V. Vicol. Onsager’s conjecture for admissible weak solutions.Comm. Pure Appl. Math., 72(2):229–274, 2019

  11. [11]

    Caffarelli, R

    L. Caffarelli, R. Kohn, and L. Nirenberg. Partial regularity of suitable weak solutions of the Navier–Stokes equations.Comm. Pure Appl. Math., 35(6):771–831, 1982

  12. [12]

    Charnyi, T

    S. Charnyi, T. Heister, M. A. Olshanskii, and L. G. Rebholz. On conservation laws of Navier–Stokes Galerkin discretizations.J. Comput. Phys., 337:289–308, 2017

  13. [13]

    S. H. Christiansen. Stability of Hodge decompositions in finite element spaces of differential forms in arbitrary dimension.Numer. Math., 107(1):87–106, 2007

  14. [14]

    Constantin, W

    P. Constantin, W. E, and E. S. Titi. Onsager’s conjecture on the energy conservation for solutions of Euler’s equation.Comm. Math. Phys., 165(1):207–209, 1994

  15. [15]

    De Lellis and L

    C. De Lellis and L. Székelyhidi Jr. The Euler equations as a differential inclusion. Ann. of Math. (2), 170(3):1417–1436, 2009

  16. [16]

    De Lellis and L

    C. De Lellis and L. Székelyhidi Jr. Dissipative continuous Euler flows.Invent. Math., 193:377–407, 2013

  17. [17]

    Desbrun, A

    M. Desbrun, A. N. Hirani, M. Leok, and J. E. Marsden. Discrete exterior calculus,

  18. [18]

    Desbrun, E

    M. Desbrun, E. Kanso, and Y. Tong. Discrete differential forms for computational modeling. InDiscrete Differential Geometry, volume 38 ofOberwolfach Seminars, pages 287–324. Birkhäuser, 2008

  19. [19]

    Dodziuk and V

    J. Dodziuk and V. K. Patodi. Riemannian structures and triangulations of manifolds. J. Indian Math. Soc. (N.S.), 40(1–4):1–52, 1976. 73

  20. [20]

    J. Droniou. Finite volume schemes for diffusion equations: introduction to and review of modern methods.Math. Models Methods Appl. Sci., 24(8):1575–1619, 2014

  21. [21]

    Duchon and R

    J. Duchon and R. Robert. Inertial energy dissipation for weak solutions of incom- pressible Euler and Navier–Stokes equations.Nonlinearity, 13(1):249–255, 2000

  22. [22]

    Elcott, Y

    S. Elcott, Y. Tong, E. Kanso, P. Schröder, and M. Desbrun. Stable, circulation- preserving, simplicial fluids.ACM Trans. Graph., 26(1):4, 2007

  23. [23]

    Eymard, T

    R. Eymard, T. Gallouët, and R. Herbin. Finite volume methods. In P. G. Ciarlet and J. L. Lions, editors,Handbook of Numerical Analysis, Vol. VII, pages 713–1020. North-Holland, Amsterdam, 2000

  24. [24]

    Girault and P.-A

    V. Girault and P.-A. Raviart.Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms, volume 5 ofSpringer Series in Computational Mathematics. Springer-Verlag, Berlin, 1986

  25. [25]

    Guermond

    J.-L. Guermond. Finite-element-based Faedo–Galerkin weak solutions to the Navier– Stokes equations in the three-dimensional torus are suitable.J. Math. Pures Appl. (9), 85(3):451–464, 2006

  26. [26]

    Guermond

    J.-L. Guermond. Faedo–Galerkin weak solutions of the Navier–Stokes equations with Dirichlet boundary conditions are suitable.J. Math. Pures Appl. (9), 88(1):87–106, 2007

  27. [27]

    Anerrorestimateforafinitevolumeschemeforadiffusion-convectionprob- lem on a triangular mesh.Numer

    R.Herbin. Anerrorestimateforafinitevolumeschemeforadiffusion-convectionprob- lem on a triangular mesh.Numer. Methods Partial Differential Equations, 11(2):165– 173, 1995

  28. [28]

    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 J. Numer. Anal., 19(2):275–311, 1982

  29. [29]

    Hiptmair

    R. Hiptmair. Discrete Hodge operators.Numer. Math., 90(2):265–289, 2001

  30. [30]

    A. N. Hirani.Discrete Exterior Calculus. PhD thesis, California Institute of Technol- ogy, 2003

  31. [31]

    Hohenegger, P

    C. Hohenegger, P. Korn, L. Linardakis, R. Redler, R. Schnur, P. Adamidis, J. Bao, S. Bastin, M. Behravesh, M. Bergemann, et al. Icon-sapphire: simulating the compo- nents of the Earth system and their interactions at kilometer and subkilometer scales. Geoscientific Model Development, 16(2):779–811, 2023

  32. [32]

    E. Hopf. Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen. Math. Nachr., 4:213–231, 1951

  33. [33]

    P. Isett. A proof of Onsager’s conjecture.Ann. of Math. (2), 188(3):871–963, 2018

  34. [34]

    P. Korn. Formulation of an unstructured grid model for global ocean dynamics.J. Comput. Phys., 339:525–552, 2017. 74

  35. [35]

    P. Korn, N. Brüggemann, J. H. Jungclaus, S. J. Lorenz, O. Gutjahr, H. Haak, L. Linar- dakis, C. Mehlmann, U. Mikolajewicz, D. Notz, W. A. Müller, D. A. Putrasahan, V. Singh, J.-S. von Storch, X. Zhu, and J. Marotzke. Icon-o: The ocean component of the icon earth system model — global simulation characteristics and local telescoping capability.Journal of A...

  36. [36]

    in the large

    O. A. Ladyzhenskaya. Solution “in the large” of the nonstationary boundary value problem for the Navier–Stokes system with two space variables.Comm. Pure Appl. Math., 12:427–433, 1959

  37. [37]

    O. A. Ladyzhenskaya. New equations for the description of the motions of viscous incompressible fluids, and global solvability for their boundary value problems.Tr. Mat. Inst. Steklova, 102:85–104, 1967. English transl.: Amer. Math. Soc. Transl. Ser. 2, vol. 102, 1973

  38. [38]

    O. A. Ladyzhenskaya. On the uniqueness and smoothness of generalized solutions of the Navier–Stokes equations.Zap. Nauchn. Sem. LOMI, 5:169–185, 1967. English transl.: Sem. in Math., Steklov Inst., vol. 5, 1969

  39. [39]

    J. Leray. Sur le mouvement d’un liquide visqueux emplissant l’espace.Acta Math., 63:193–248, 1934

  40. [40]

    Lions.Quelques méthodes de résolution des problèmes aux limites non linéaires

    J.-L. Lions.Quelques méthodes de résolution des problèmes aux limites non linéaires. Dunod; Gauthier-Villars, Paris, 1969

  41. [41]

    J. E. Marsden and T. S. Ratiu.Introduction to Mechanics and Symmetry, volume 17 ofTexts in Applied Mathematics. Springer, 2nd edition, 1999

  42. [42]

    H. K. Moffatt. The degree of knottedness of tangled vortex lines.J. Fluid Mech., 35(1):117–129, 1969

  43. [43]

    Mullen, K

    P. Mullen, K. Crane, D. Pavlov, Y. Tong, and M. Desbrun. Energy-preserving inte- grators for fluid animation. InACM SIGGRAPH 2009 Papers, pages 1–8, 2009

  44. [44]

    Pavlov, P

    D. Pavlov, P. Mullen, Y. Tong, E. Kanso, J. E. Marsden, and M. Desbrun. Structure- preserving discretization of incompressible fluids.Phys. D, 240(6):443–458, 2011

  45. [45]

    Pedregal.Parametrized Measures and Variational Principles, volume 30 ofProgress in Nonlinear Differential Equations and their Applications

    P. Pedregal.Parametrized Measures and Variational Principles, volume 30 ofProgress in Nonlinear Differential Equations and their Applications. Birkhäuser, Basel, 1997

  46. [46]

    Conservationpropertiesofunstructuredstaggeredmeshschemes.J

    B.Perot. Conservationpropertiesofunstructuredstaggeredmeshschemes.J. Comput. Phys., 159(1):58–89, 2000

  47. [47]

    T. D. Ringler, J. Thuburn, J. B. Klemp, and W. C. Skamarock. A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C- grids.J. Comput. Phys., 229(9):3065–3090, 2010

  48. [48]

    Scheffer

    V. Scheffer. An inviscid flow with compact support in space-time.J. Geom. Anal., 3(4):343–401, 1993

  49. [49]

    Shnirelman

    A. Shnirelman. On the nonuniqueness of weak solution of the Euler equation.Comm. Pure Appl. Math., 50(12):1261–1286, 1997. 75

  50. [50]

    Smagorinsky

    J. Smagorinsky. General circulation experiments with the primitive equations. I. The basic experiment.Mon. Weather Rev., 91:99–164, 1963

  51. [51]

    Compensatedcompactnessandapplicationstopartialdifferentialequations

    L.Tartar. Compensatedcompactnessandapplicationstopartialdifferentialequations. In R. J. Knops, editor,Nonlinear Analysis and Mechanics: Heriot–Watt Symposium, Vol. IV, volume 39 ofResearch Notes in Mathematics, pages 136–212. Pitman, Boston, 1979

  52. [52]

    Vohralík

    M. Vohralík. Residual flux-based a posteriori error estimates for finite volume and related locally conservative methods.Numer. Math., 111(1):121–158, 2008

  53. [53]

    Whitney.Geometric Integration Theory

    H. Whitney.Geometric Integration Theory. Princeton University Press, Princeton, NJ, 1957

  54. [54]

    Wolibner

    W. Wolibner. Un théorème sur l’existence du mouvement plan d’un fluide parfait, homogène, incompressible, pendant un temps infiniment long.Math. Z., 37(1):698– 726, 1933

  55. [55]

    Zhang, A

    Y. Zhang, A. Palha, M. Gerritsma, and L. G. Rebholz. A mass-, kinetic energy- and helicity-conserving mimetic dual-field discretization for three-dimensional incompress- ible Navier–Stokes equations, Part I: Periodic domains.J. Comput. Phys., 451:110868, 2022. 76