arxiv: 2605.13048 · v1 · submitted 2026-05-13 · 🧮 math.AP · cs.NA· math.NA

Recognition: unknown

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

Peter Korn

Authors on Pith no claims yet

Pith reviewed 2026-05-14 18:37 UTC · model grok-4.3

classification 🧮 math.AP cs.NAmath.NA

keywords discrete exterior calculusincompressible Navier-StokesOnsager thresholdstructure-preserving discretizationmeasure-valued solutionsenergy conservationweak solutionsEuler equations

0 comments

The pith

Exact algebraic conservation in a discrete scheme rules out dissipative weak solutions of 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 develops a structure-preserving discretization of the incompressible Navier-Stokes equations using discrete exterior calculus on prismatic Delaunay-Voronoi meshes over closed Riemannian manifolds. Exact algebraic conservation at the discrete level serves as a selection principle that excludes entire classes of weak solutions reachable by other discretizations. In the smooth regime, solutions converge at explicit rates uniform in viscosity. In the Leray-Hopf regime, limits satisfy the viscous weak equations. Above the Onsager threshold in the inviscid regime, limits become conservative measure-valued Euler solutions under a uniform Holder bound, while no subsequence can converge to an energy-dissipating Euler solution at any regularity.

Core claim

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. In the inviscid measure-valued regime, limits are conservative measure-valued Euler solutions whose concentration defect vanishes for alpha greater than 1/3 provided the discrete solutions admit a uniform C^{0,alpha} bound. In the dissipative regime, no subsequence converges to an energy-dissipating Euler solution at any regularity, a structural exclusion that follows from exact discrete energy conservation.

What carries the argument

The exact algebraic conservation of energy and other quantities enforced by the discrete exterior calculus structure on prismatic Delaunay-Voronoi meshes, which produces algebraic identities that prevent discrete energy dissipation.

Load-bearing premise

The discrete solutions admit a uniform Holder bound of order alpha greater than one third in the inviscid regime.

What would settle it

A sequence of discrete solutions that converges to a weak Euler solution dissipating positive kinetic energy at any regularity would falsify the structural exclusion of dissipative limits.

Figures

Figures reproduced from arXiv: 2605.13048 by Peter Korn.

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. % from all Galerkin and finite-volume methods. 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 paper 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 exact algebraic conservation of energy and other invariants at the discrete level, which is used to derive a selection principle: subsequential limits are conservative measure-valued Euler solutions whose concentration defect vanishes for Hölder regularity α > 1/3 (provided a uniform C^{0,α} bound on the discrete velocities), while energy-dissipating weak solutions are excluded in the dissipative regime. Additional results include convergence rates O(h^{min(r_rec, r_*)} |log h|^β_d) for smooth solutions (uniform in viscosity) and subsequential convergence to Leray-Hopf weak solutions in the viscous case.

Significance. If the central claims hold, the work is significant for demonstrating that exact discrete conservation provides a parameter-free mechanism to rule out entire classes of weak solutions (including dissipative ones) that are reachable by standard Galerkin or finite-volume schemes. The conditional result linking the Onsager threshold to vanishing concentration defect, combined with the uniform-in-viscosity convergence for smooth data, offers a rigorous framework for structure-preserving methods in fluid dynamics. The exact algebraic conservation is a clear technical strength.

major comments (1)

[Inviscid measure-valued regime] Inviscid measure-valued regime (as described in the abstract): the assertion that the concentration defect vanishes for α > 1/3 is conditional on the discrete velocity fields admitting a uniform C^{0,α} bound independent of h. No discrete Sobolev embedding, maximum principle, or a priori estimate is indicated that would deliver this bound, rendering the defect-free convergence claim dependent on an unverified assumption that is load-bearing for the selection principle in this regime.

minor comments (2)

[Smooth solutions] The convergence rate statement includes the factor |log h|^β_d with β_3 = 0 and β_2 = 1; the dependence of β_d on dimension and the precise definition of r_rec and r_* should be stated explicitly in the main theorem for smooth solutions.
[Introduction] The abstract notes that the scheme distinguishes itself from Galerkin and finite-volume methods, but the introduction or comparison section should include a brief table or paragraph contrasting the conservation properties and reachable weak solutions with at least one standard scheme.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful and constructive review. The central observation regarding the conditional nature of the defect-vanishing result is addressed below. We provide a point-by-point response and indicate where the manuscript will be revised for clarity.

read point-by-point responses

Referee: Inviscid measure-valued regime (as described in the abstract): the assertion that the concentration defect vanishes for α > 1/3 is conditional on the discrete velocity fields admitting a uniform C^{0,α} bound independent of h. No discrete Sobolev embedding, maximum principle, or a priori estimate is indicated that would deliver this bound, rendering the defect-free convergence claim dependent on an unverified assumption that is load-bearing for the selection principle in this regime.

Authors: We agree that the vanishing of the concentration defect is stated conditionally on a uniform C^{0,α} bound independent of h, as explicitly indicated in the abstract and in the statement of the main theorem for the inviscid measure-valued regime. No discrete Sobolev embedding, maximum principle, or a priori Hölder estimate is provided in the paper, and establishing such a bound remains an open question (analogous to the continuous setting). The exact algebraic conservation nevertheless yields an unconditional structural result: no subsequence can converge to an energy-dissipating Euler weak solution at any regularity. The conditional defect-free convergence above the Onsager threshold is presented as a selection principle that holds under the stated hypothesis. We will revise the introduction and the discussion of the measure-valued regime to emphasize the conditional character more explicitly and to frame the uniform bound as an important open problem for future work. revision: partial

Circularity Check

0 steps flagged

No circularity; selection principle follows from engineered exact conservation without definitional reduction or fitted predictions

full rationale

The paper's derivation begins with a discretization explicitly constructed via discrete exterior calculus to enforce exact algebraic conservation of energy and other quantities on the given meshes. This built-in property is then used to establish the selection principle excluding dissipative weak solutions in the inviscid regime, with the Onsager-threshold claim stated conditionally on an external uniform C^{0,α} bound rather than derived from the conservation alone. No step equates a claimed prediction to a fitted input by construction, invokes a self-citation as the sole justification for a uniqueness theorem, or renames a known result as a new unification. The chain remains self-contained under the paper's stated mesh and regularity assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The theory depends on the exact conservation property of the chosen discretization and on mesh and regularity assumptions; no free parameters are introduced in the abstract.

axioms (2)

domain assumption Discretization is performed on prismatic Delaunay-Voronoi meshes over closed Riemannian manifolds
The entire scheme and convergence statements are defined for these specific meshes.
ad hoc to paper Discrete solutions admit a uniform C^{0,α} bound in the inviscid regime
This bound is required for the claim that concentration defect vanishes above the Onsager threshold.

pith-pipeline@v0.9.0 · 5607 in / 1385 out tokens · 48719 ms · 2026-05-14T18:37:50.715637+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

[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

work page 2006
[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

work page 2010
[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

work page 1966
[4]

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

work page 2006
[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

work page 1988
[6]

Bossavit

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

work page 1998
[7]

Bossavit

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

work page 2001
[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

work page 2011
[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

work page 2005
[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

work page 2019
[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

work page 1982
[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

work page 2017
[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

work page 2007
[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

work page 1994
[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

work page 2009
[16]

De Lellis and L

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

work page 2013
[17]

Desbrun, A

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

work page
[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

work page 2008
[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

work page 1976
[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

work page 2014
[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

work page 2000
[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

work page 2007
[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

work page 2000
[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

work page 1986
[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

work page 2006
[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

work page 2007
[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

work page 1995
[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

work page 1982
[29]

Hiptmair

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

work page 2001
[30]

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

work page 2003
[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

work page 2023
[32]

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

work page 1951
[33]

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

work page 2018
[34]

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

work page 2017
[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...

work page 2022
[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

work page 1959
[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

work page 1967
[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

work page 1967
[39]

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

work page 1934
[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

work page 1969
[41]

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

work page 1999
[42]

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

work page 1969
[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

work page 2009
[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

work page 2011
[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

work page 1997
[46]

Conservationpropertiesofunstructuredstaggeredmeshschemes.J

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

work page 2000
[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

work page 2010
[48]

Scheffer

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

work page 1993
[49]

Shnirelman

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

work page 1997
[50]

Smagorinsky

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

work page 1963
[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

work page 1979
[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

work page 2008
[53]

Whitney.Geometric Integration Theory

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

work page 1957
[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

work page 1933
[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

work page 2022