A no-go theorem and its resolution for the discrete compressible barotropic Navier--Stokes equations

Peter Korn

arxiv: 2605.16554 · v3 · pith:QNCDMLV3new · submitted 2026-05-15 · 🧮 math.AP · cs.NA· math.NA

A no-go theorem and its resolution for the discrete compressible barotropic Navier--Stokes equations

Peter Korn This is my paper

Pith reviewed 2026-05-25 06:18 UTC · model grok-4.3

classification 🧮 math.AP cs.NAmath.NA

keywords discrete Navier-Stokes equationsenergy conservationvector-invariant formno-go theoremdensity-weighted mass matrixHollingsworth instabilityDelaunay-Voronoi meshbarotropic flow

0 comments

The pith

Density-independent mass matrices produce an unavoidable O(h²) energy residual in discrete compressible barotropic Navier-Stokes equations on Delaunay-Voronoi meshes.

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

The paper proves a no-go theorem showing that on Delaunay-Voronoi meshes, any density-independent mass matrix consistent with integration-by-parts for the divergence operator leads to a sharp second-order energy residual whose sign is indeterminate and cannot be eliminated by operator choices. This holds for various grid staggerings including A, B, C, D, and quasi-B grids. The resolution is the density-weighted mass matrix, which is the unique choice that achieves exact total energy conservation while preserving the vector-invariant momentum equation, Lamb antisymmetry, and topological conservation laws. This construction also enables global well-posedness, convergence to smooth solutions, and Lyapunov stability around specific equilibria, excluding the Hollingsworth instability.

Core claim

Every density-independent mass matrix with integration-by-parts-consistent divergence on a Delaunay-Voronoi mesh carries a sharp O(h²) energy residual of indeterminate sign that no operator choice can eliminate. The density-weighted mass matrix is the unique algebraic remedy that restores exact total energy while preserving the vector-invariant momentum equation, Lamb antisymmetry, and topological conservation laws, at the cost of an O(h^{r*}) Kelvin defect.

What carries the argument

The density-weighted mass matrix, which multiplies the mass matrix entries by local density values to cancel the energy residual term while maintaining consistency with the divergence operator.

If this is right

Exact discrete total energy conservation is restored without altering the vector-invariant form or topological properties.
The scheme achieves global well-posedness for nonnegative viscosity and converges uniformly to smooth solutions.
Asymptotic preservation holds in the low-Mach limit, unlike the density-free case where the residual diverges as O(M^{-1}).
Lyapunov stability is guaranteed around hydrostatic and constant-flow stratified states, and conditionally around sheared baroclinic states via a discrete Charney-Stern criterion.
Baroclinic instability appears only when the continuous system permits it.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The density-weighting principle could apply to other mesh families provided an integration-by-parts identity is available.
Structural elimination of the residual may improve statistical fidelity in long integrations of atmospheric and ocean models.
The same construction might extend to non-barotropic or compressible Euler cases while retaining the listed conservation properties.
Direct computation of the energy residual on concrete meshes would serve as an immediate check on the theorem's sharpness.

Load-bearing premise

The mesh is a Delaunay-Voronoi tessellation in two or three dimensions and the discrete operators satisfy the standard integration-by-parts identity for the chosen staggering.

What would settle it

A numerical test on a Delaunay-Voronoi mesh with a density-independent mass matrix that measures whether total energy shows an O(h²) deviation whose sign varies with different flows or initial conditions.

read the original abstract

The compressible barotropic Navier--Stokes equations in vector-invariant form preserve the vorticity structure of the system and underlie modern atmospheric and ocean dynamical cores, yet no PDE theory has been developed for the compressible discrete system in this form. On a Delaunay--Voronoi mesh we prove via discrete exterior calculus, that every density-independent mass matrix with integration-by-parts-consistent divergence carries a sharp $\OO(h^2)$ energy residual of indeterminate sign that no operator choice can eliminate. This no-go theorem covers A-, B-, C-, D-, and quasi-B-grid staggerings. The density-weighted mass matrix is the unique algebraic remedy: it restores exact total energy while preserving the vector-invariant momentum equation, Lamb antisymmetry, and the topological conservation laws, at the cost of an $\OO(h^{r_\star})$ Kelvin defect matching the convergence rate. The residual is the cause of the Hollingsworth instability that has shaped vector-invariant dynamical-core design; the density-weighted construction removes it structurally. For the density-weighted~(DW) scheme on closed oriented Riemannian manifolds in $d = 2, 3$ we establish global well-posedness for $\nu \ge 0$, convergence to smooth solutions uniformly in $\nu$, and asymptotic preservation in the low-Mach limit; the density-free residual diverges as $\OO(M^{-1})$. Via a discrete Arnold energy-Casimir construction, exact discrete conservation forces Lyapunov stability around three classes of equilibria, excluding Hollingsworth instability: unconditional stability around hydrostatic and constant-flow stratified states, and conditional stability around sheared baroclinic states under a discrete Charney--Stern criterion. The DW scheme admits genuine baroclinic instability only when the continuum equations themselves do.

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.

Desk Editor's Note private letter to a colleague

The paper proves a no-go for exact energy conservation with any density-independent mass matrix on Delaunay-Voronoi meshes in vector-invariant discrete barotropic NS, then shows the density-weighted version fixes it and yields well-posedness plus Arnold stability.

read the letter

Hey, the core result here is a no-go theorem: on Delaunay-Voronoi meshes, every density-independent mass matrix that respects integration-by-parts for the divergence produces an O(h²) energy residual of indeterminate sign, and no choice of operators removes it while keeping the vector-invariant form. The density-weighted mass matrix is presented as the unique algebraic fix that restores exact total energy, keeps Lamb antisymmetry and topological conservation, at the price of an O(h^r) Kelvin defect. They then prove global well-posedness for the resulting scheme when viscosity is nonnegative, uniform convergence to smooth solutions, low-Mach asymptotic preservation, and discrete Arnold energy-Casimir stability around hydrostatic, constant-flow, and certain sheared equilibria under a discrete Charney-Stern condition. This directly ties the residual to the Hollingsworth instability seen in operational models. The no-go and the explicit density-weighted construction look new relative to the cited literature on discrete schemes. The work is internally consistent on its own terms and connects a practical modeling headache to a structural algebraic cause. The main limitation is the reliance on exact integration-by-parts on Delaunay-Voronoi meshes for the specific staggerings listed; the stress-test note is right that relaxing the IBP identity or changing the mesh type could open other algebraic routes, so the uniqueness claim is conditional on those assumptions. The well-posedness and stability arguments are strong in outline but would need the full discrete exterior calculus derivations checked for gaps. This is for people working on structure-preserving discretizations or geophysical dynamical cores. It deserves a serious referee because it supplies a concrete remedy backed by stability theory rather than just another scheme. I'd send it out for review.

Referee Report

2 major / 2 minor

Summary. The paper proves via discrete exterior calculus on Delaunay-Voronoi meshes that every density-independent mass matrix consistent with integration-by-parts for the divergence operator produces a sharp O(h²) energy residual of indeterminate sign in the vector-invariant discrete compressible barotropic Navier-Stokes system; this covers A/B/C/D/quasi-B staggerings and explains Hollingsworth instability. It shows that the density-weighted mass matrix is the unique algebraic fix restoring exact total energy conservation while preserving the vector-invariant form, Lamb antisymmetry, and topological laws (at the cost of an O(h^{r*}) Kelvin defect). For the resulting DW scheme the authors establish global well-posedness for ν ≥ 0, convergence to smooth solutions uniformly in ν, low-Mach asymptotic preservation, and discrete Arnold energy-Casimir Lyapunov stability around hydrostatic, constant-flow, and (conditionally) sheared baroclinic equilibria.

Significance. If the no-go theorem and its resolution hold, the work supplies a structurally stable, energy-exact discretization for a class of models central to atmospheric and ocean dynamical cores. The discrete energy-Casimir construction yielding unconditional stability for hydrostatic and constant-flow states and conditional stability under a discrete Charney-Stern criterion is a concrete strength; the explicit identification of the residual as the source of Hollingsworth instability and its removal by density weighting is a clear advance over ad-hoc fixes.

major comments (2)

[discrete exterior calculus / no-go theorem] The no-go theorem (discrete exterior calculus section) derives the O(h²) residual of indeterminate sign from the exact integration-by-parts identity between the chosen mass matrix and divergence on Delaunay-Voronoi meshes. The manuscript must state explicitly whether this identity is verified for each listed staggering or whether it is an assumption; if the identity is relaxed to an approximate relation, the claim that no density-independent matrix can eliminate the residual becomes conditional and the uniqueness of the density-weighted remedy requires a separate argument.
[global well-posedness and convergence section] Global well-posedness and convergence statements for the DW scheme claim uniform-in-ν convergence to smooth solutions together with an O(h^{r*}) Kelvin defect. The estimates must clarify how the defect term is controlled in the energy estimates so that it does not prevent uniformity in ν or degrade the convergence rate; without this control the uniform-convergence claim is not yet load-bearing.

minor comments (2)

[abstract / low-Mach paragraph] The low-Mach statement that the density-free residual diverges as O(M^{-1}) should be derived from the discrete energy equation rather than asserted; the symbol M should be defined at first use.
[abstract] Notation r_star for the Kelvin defect order should be introduced with a reference to the underlying approximation order of the discrete operators.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on the no-go theorem and the well-posedness analysis. We address each major comment below and will incorporate the requested clarifications in the revised manuscript.

read point-by-point responses

Referee: [discrete exterior calculus / no-go theorem] The no-go theorem (discrete exterior calculus section) derives the O(h²) residual of indeterminate sign from the exact integration-by-parts identity between the chosen mass matrix and divergence on Delaunay-Voronoi meshes. The manuscript must state explicitly whether this identity is verified for each listed staggering or whether it is an assumption; if the identity is relaxed to an approximate relation, the claim that no density-independent matrix can eliminate the residual becomes conditional and the uniqueness of the density-weighted remedy requires a separate argument.

Authors: The integration-by-parts identity is a direct consequence of the discrete exterior calculus construction on Delaunay-Voronoi meshes (Section 2), not an assumption imposed separately on each staggering. All listed staggerings (A/B/C/D/quasi-B) are instances of this DEC framework, so the identity holds by construction for any density-independent mass matrix consistent with the DEC divergence. We will add an explicit statement in the revised text confirming that the identity is verified within the DEC setup and that the no-go result therefore applies uniformly to the listed staggerings. revision: yes
Referee: [global well-posedness and convergence section] Global well-posedness and convergence statements for the DW scheme claim uniform-in-ν convergence to smooth solutions together with an O(h^{r*}) Kelvin defect. The estimates must clarify how the defect term is controlled in the energy estimates so that it does not prevent uniformity in ν or degrade the convergence rate; without this control the uniform-convergence claim is not yet load-bearing.

Authors: The O(h^{r*}) Kelvin defect is controlled in the energy estimates of Section 4 by absorbing it into the discrete dissipation term using the a priori bounds furnished by the discrete energy-Casimir functional; because the defect order matches the spatial approximation rate, it does not introduce ν-dependent growth and preserves uniformity in ν as well as the convergence rate. We will insert a dedicated paragraph in the revised Section 4 making this absorption step fully explicit. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation is self-contained via discrete exterior calculus on stated mesh assumptions.

full rationale

The no-go theorem is derived directly from the exact integration-by-parts identity assumed for density-independent mass matrices on Delaunay-Voronoi meshes, using discrete exterior calculus identities for the listed staggerings. The density-weighted remedy is exhibited algebraically as restoring exact energy conservation while preserving other properties, without reducing to a fitted parameter, self-definition, or self-citation chain. The result is conditional on the IBP assumption and mesh type, but this is an explicit premise rather than a hidden reduction. No load-bearing step equates the target conclusion to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The argument relies on standard properties of discrete exterior calculus (integration-by-parts identity, topological conservation laws) and the geometric assumption that the mesh is Delaunay-Voronoi; no free parameters are introduced and no new entities are postulated.

axioms (2)

standard math Discrete exterior calculus operators on Delaunay-Voronoi meshes satisfy the integration-by-parts identity for the chosen staggering.
Invoked to derive the O(h²) energy residual for any density-independent mass matrix.
domain assumption The vector-invariant form of the compressible barotropic Navier-Stokes equations is the starting point.
The entire analysis is performed in this form; the no-go and resolution are specific to it.

pith-pipeline@v0.9.0 · 5845 in / 1568 out tokens · 32811 ms · 2026-05-25T06:18:36.465478+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean absolute_floor_iff_bare_distinguishability unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

every density-independent mass matrix with integration-by-parts-consistent divergence carries a sharp O(h²) energy residual of indeterminate sign that no operator choice can eliminate
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

density-weighted mass matrix Mρ₁ = Pᵀ diag(ρ) P restores exact total energy while preserving vector-invariant momentum equation, Lamb antisymmetry

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.