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
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.
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
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.
Referee Report
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)
- [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)
- [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.
- [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
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
-
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
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
axioms (2)
- domain assumption Exact algebraic conservation holds for the discrete exterior calculus operators on the specified prismatic Delaunay-Voronoi meshes
- domain assumption Mesh properties including centroid proximity and reconstruction symmetry enable the stated convergence rates
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
exact algebraic conservation at the discrete level is not merely a fidelity property but rules out entire classes of weak solutions... Dissipative regime: no subsequence converges to an energy-dissipating Euler solution at any regularity
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Inviscid measure-valued regime: limits are conservative measure-valued Euler solutions; their concentration defect vanishes above the Onsager threshold α>1/3 provided the discrete solutions admit a uniform C^{0,α} bound
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
-
[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]
- [6]
- [7]
-
[8]
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]
-
[10]
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]
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]
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]
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]
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]
C. De Lellis and L. Székelyhidi Jr. Dissipative continuous Euler flows.Invent. Math., 193:377–407, 2013
work page 2013
-
[17]
M. Desbrun, A. N. Hirani, M. Leok, and J. E. Marsden. Discrete exterior calculus,
-
[18]
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]
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]
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]
- [23]
-
[24]
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]
- [26]
-
[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]
-
[30]
A. N. Hirani.Discrete Exterior Calculus. PhD thesis, California Institute of Technol- ogy, 2003
work page 2003
-
[31]
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]
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]
- [44]
-
[45]
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]
-
[49]
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]
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]
-
[53]
Whitney.Geometric Integration Theory
H. Whitney.Geometric Integration Theory. Princeton University Press, Princeton, NJ, 1957
work page 1957
- [54]
- [55]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.