pith. sign in

arxiv: 2605.16554 · v2 · 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
0 comments X

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.

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

  1. 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

  2. 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)

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

56 extracted references · 56 canonical work pages · 1 internal anchor

  1. [1]

    Amann.Ordinary Differential Equations: An Introduction to Nonlinear Analysis, volume 13 ofDe Gruyter Studies in Mathematics

    H. Amann.Ordinary Differential Equations: An Introduction to Nonlinear Analysis, volume 13 ofDe Gruyter Studies in Mathematics. Walter de Gruyter, Berlin, 1990. Translated from the German by Gerhard Metzen

    work page 1990

  2. [2]

    S. A. Anandan and G. Kanschat. Structure preserving schemes for the compressible Euler equations.Preprint, 2025

    work page 2025

  3. [3]

    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

  4. [4]

    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. 106

    work page 2010

  5. [5]

    V. I. Arnold. Conditions for nonlinear stability of stationary plane curvilinear flows of an ideal fluid.Dokl. Akad. Nauk SSSR, 162(5):975–978, 1965. English transl.: Soviet Math. Dokl. 6 (1965), 773–777

    work page 1965

  6. [6]

    Onanaprioriestimateinthetheoryofhydrodynamicstability.American Mathematical Society Translations, Series 2, 79:267–269, 1969

    V.I.Arnold. Onanaprioriestimateinthetheoryofhydrodynamicstability.American Mathematical Society Translations, Series 2, 79:267–269, 1969. Translation of Izv. Vyssh. Uchebn. Zaved. Matematika 54, no. 5 (1966), 3–5

    work page 1969

  7. [7]

    M. J. Bell, P. S. Peixoto, and J. Thuburn. Numerical instabilities of vector-invariant momentum equations on rectangular C-grids.Quarterly Journal of the Royal Meteo- rological Society, 143(702):563–581, 2017

    work page 2017

  8. [8]

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

    work page 2006

  9. [9]

    Bossavit.Computational Electromagnetism: Variational Formulations, Comple- mentarity, Edge Elements

    A. Bossavit.Computational Electromagnetism: Variational Formulations, Comple- mentarity, Edge Elements. Academic Press, San Diego, 1998

    work page 1998

  10. [10]

    Bossavit

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

    work page 2001

  11. [11]

    J. G. Charney and M. E. Stern. On the stability of internal baroclinic jets in a rotating atmosphere.Journal of the Atmospheric Sciences, 19(2):159–172, 1962

    work page 1962

  12. [12]

    F. H. Clarke.Optimization and Nonsmooth Analysis, volume 5 ofClassics in Applied Mathematics. SIAM, Philadelphia, 1990. Reprint of the 1983 Wiley original

    work page 1990

  13. [13]

    Discrete Exterior Calculus

    M. Desbrun, A. N. Hirani, M. Leok, and J. E. Marsden. Discrete exterior calculus. arXiv:math/0508341 [math.DG], 2005

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  14. [14]

    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

  15. [15]

    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

    work page 1976

  16. [16]

    Dubos and N

    T. Dubos and N. K.-R. Kevlahan. A conservative adaptive wavelet method for the shallow-water equations on staggered grids.Quart. J. Roy. Meteor. Soc., 141:1712– 1726, 2015

    work page 2015

  17. [17]

    Symmetric Instability of the Computational Kind

    N. Ducousso, J. Le Sommer, J.-M. Molines, and M. Bell. Impact of the “Symmetric Instability of the Computational Kind” at mesoscale- and submesoscale-permitting resolutions.Ocean Modelling, 120:18–26, 2017

    work page 2017

  18. [18]

    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

  19. [19]

    Feireisl, M

    E. Feireisl, M. Lukavcová-Medvivdová, H. Mizerová, and B. She.Numerical Anal- ysis of Compressible Fluid Flows, volume 20 ofMS&A – Modeling, Simulation and Applications. Springer, Cham, 2021. 107

    work page 2021

  20. [20]

    Gallouët, R

    T. Gallouët, R. Herbin, and J.-C. Latché. A convergent finite element–finite volume scheme for the compressible Stokes problem. I. The isothermal case.Math. Comp., 78(267):1333–1352, 2009. Published online 2008

    work page 2009

  21. [21]

    Gallouët, R

    T. Gallouët, R. Herbin, and J.-C. Latché. A convergent finite element–finite volume scheme for the compressible Stokes problem. Part I: The isothermal case.Math. Comp., 78:1333–1352, 2009

    work page 2009

  22. [22]

    Gallouët, R

    T. Gallouët, R. Herbin, and J.-C. Latché. A convergent finite element–finite volume scheme for the compressible Stokes problem. Part II: The isentropic case.Math. Comp., 79:649–675, 2010

    work page 2010

  23. [23]

    Gallouët, D

    T. Gallouët, D. Maltese, and A. Novotný. Error estimates for a numerical approxima- tion to the compressible barotropic Navier–Stokes equations.IMA J. Numer. Anal., 36:543–592, 2016

    work page 2016

  24. [24]

    Gassmann

    A. Gassmann. A global hexagonal C-grid non-hydrostatic dynamical core (ICON- IAP) designed for energetic consistency.Quarterly Journal of the Royal Meteorological Society, 139(670):152–175, 2013

    work page 2013

  25. [25]

    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, Berlin, 1986

    work page 1986

  26. [26]

    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

  27. [27]

    Hiptmair

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

    work page 2001

  28. [28]

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

    work page 2003

  29. [29]

    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

  30. [30]

    Hollingsworth, P

    A. Hollingsworth, P. Kållberg, V. Renner, and D. M. Burridge. An internal symmet- ric computational instability.Quarterly Journal of the Royal Meteorological Society, 109(460):417–428, 1983

    work page 1983

  31. [31]

    D. D. Holm, J. E. Marsden, T. Ratiu, and A. Weinstein. Nonlinear stability of fluid and plasma equilibria.Physics Reports, 123(1–2):1–116, 1985

    work page 1985

  32. [32]

    B. S. Jovanović and E. Süli.Analysis of Finite Difference Schemes for Linear Partial Differential Equations with Generalized Solutions, volume 46 ofSpringer Series in Computational Mathematics. Springer, London, 2014

    work page 2014

  33. [33]

    T. K. Karper. A convergent FEM–DG method for the compressible Navier–Stokes equations.Numer. Math., 125:441–510, 2013. 108

    work page 2013

  34. [34]

    Klainerman and A

    S. Klainerman and A. Majda. Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit of compressible fluids.Comm. Pure Appl. Math., 34(4):481–524, 1981

    work page 1981

  35. [35]

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

    work page 2017

  36. [36]

    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

  37. [37]

    Aconservativediscretizationoftheshallow-waterequations on triangular grids.Journal of Computational Physics, 375:871–900, 2018

    P.KornandL.Linardakis. Aconservativediscretizationoftheshallow-waterequations on triangular grids.Journal of Computational Physics, 375:871–900, 2018. Introduces the density-weighted mass matrix used in the companion papers on the barotropic and fully compressible systems

    work page 2018

  38. [38]

    O. A. Ladyzhenskaya.New equations for the description of the motions of viscous incompressible fluids and global solvability for their boundary value problems, volume

    work page

  39. [39]

    English translation in Proc

    1967. English translation in Proc. Steklov Inst. Math. 102 (1967), 95–118

    work page 1967

  40. [40]

    O. A. Ladyzhenskaya.On some nonlinear problems in the theory of continuous media. Mir, Moscow, 1968

    work page 1968

  41. [41]

    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

  42. [42]

    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

  43. [43]

    Melvin, B

    T. Melvin, B. Shipway, N. Wood, T. Benacchio, T. Bendall, I. Boutle, A. Brown, C. Johnson, J. Kent, S. Pring, C. Smith, M. Zerroukat, C. Cotter, and J. Thuburn. A mixed finite-element, finite-volume, semi-implicit discretisation for atmospheric dy- namics: Spherical geometry.Quarterly Journal of the Royal Meteorological Society, 150(764):4252–4269, 2024

    work page 2024

  44. [44]

    Métivier and S

    G. Métivier and S. Schochet. The incompressible limit of the non-isentropic Euler equations.Arch. Ration. Mech. Anal., 158(1):61–90, 2001

    work page 2001

  45. [45]

    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

  46. [46]

    Noelle, G

    S. Noelle, G. Bispen, K. R. Arun, M. Lukávcová-Medvivdová, and C.-D. Munz. A weakly asymptotic preserving low Mach number scheme for the Euler equations of gas dynamics.SIAM J. Sci. Comput., 36(6):B989–B1024, 2014

    work page 2014

  47. [47]

    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

  48. [48]

    Pedlosky.Geophysical Fluid Dynamics

    J. Pedlosky.Geophysical Fluid Dynamics. Springer, New York, 2nd edition, 1987. 109

    work page 1987

  49. [49]

    Rayleigh

    L. Rayleigh. On the stability, or instability, of certain fluid motions.Proceedings of the London Mathematical Society, s1-11(1):57–72, 1880

    work page

  50. [50]

    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

  51. [51]

    Schochet

    S. Schochet. The compressible Euler equations in a bounded domain: existence of solutions and the incompressible limit.Comm. Math. Phys., 104(1):49–75, 1986

    work page 1986

  52. [52]

    Schochet

    S. Schochet. Fast singular limits of hyperbolic PDEs.J. Differential Equations, 114(2):476–512, 1994

    work page 1994

  53. [53]

    W. C. Skamarock, J. B. Klemp, M. G. Duda, L. D. Fowler, S.-H. Park, and T. D. Ringler. A multiscale nonhydrostatic atmospheric model using centroidal Voronoi tesselations and C-grid staggering.Monthly Weather Review, 140(9):3090–3105, 2012

    work page 2012

  54. [54]

    Thuburn and C

    J. Thuburn and C. J. Cotter. A framework for mimetic discretization of the rotat- ing shallow-water equations on arbitrary polygonal grids.SIAM J. Sci. Comput., 34(3):B203–B225, 2012

    work page 2012

  55. [55]

    Whitney.Geometric Integration Theory

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

    work page 1957

  56. [56]

    Zängl, D

    G. Zängl, D. Reinert, P. Rípodas, and M. Baldauf. The ICON (ICOsahedral Non- hydrostatic) modelling framework of DWD and MPI-M: Description of the non- hydrostatic dynamical core.Quarterly Journal of the Royal Meteorological Society, 141(687):563–579, 2015. 110

    work page 2015