Conservative and skew-symmetric forms of the incompressible Navier-Stokes equations in sigma-coordinates

Jaeyoung Jung; Marco Giometto

arxiv: 2604.22060 · v1 · submitted 2026-04-23 · ⚛️ physics.flu-dyn

Conservative and skew-symmetric forms of the incompressible Navier-Stokes equations in sigma-coordinates

Jaeyoung Jung , Marco Giometto This is my paper

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

classification ⚛️ physics.flu-dyn

keywords equationscartesianconservativeformulationsskew-symmetricderivedformformulation

0 comments

The pith

New conservative and skew-symmetric formulations of incompressible Navier-Stokes equations in sigma-coordinates preserve energy properties and avoid metric-induced structural disruptions.

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

The work focuses on rewriting the equations for fluid flow in a special coordinate system called sigma-coordinates, which follows the shape of the terrain below. Standard ways of changing coordinates often break important mathematical properties that help keep simulations stable and physically accurate. The authors create a conservative version that follows general rules for conservation laws and then analyze how its wave-like behaviors differ from the usual flat-coordinate version. They also introduce a new set of variables to make a skew-symmetric form. This version keeps energy conserved when there is no friction and keeps energy from growing when friction is present. They also talk about special boundary conditions based on wave characteristics to maintain these good properties. Because only the abstract is available, the exact equations and proofs cannot be checked here, but the approach aims to fix issues that arise when terrain-following grids are used in computer models of air or water movement.

Core claim

A skew-symmetric formulation is derived by introducing a new set of variables, yielding a form that is energy-conserving for the Euler equations and energy-bounded for the Navier-Stokes equations.

Load-bearing premise

That the sigma-transformation and new variables can be chosen such that metric terms do not disrupt intrinsic structural properties, which may not hold for all terrain configurations or discretizations.

read the original abstract

This study derives conservative and skew-symmetric formulations of the incompressible flow equations in a terrain-following sigma-coordinate system that preserve key structural properties of the Cartesian formulation. Unlike conventional formulations based on the direct application of the sigma-transformation to Cartesian equations, in which metric-induced terms disrupt the intrinsic structure of the governing equations, the proposed formulations are designed to avoid these structural inconsistencies. A conservative form is derived in a manner consistent with general conservation laws, and its modified eigenstructure is analyzed relative to the Cartesian counterpart. A skew-symmetric formulation is then derived by introducing a new set of variables, yielding a form that is energy-conserving for the Euler equations and energy-bounded for the Navier-Stokes equations. Finally, we discuss characteristic-based boundary conditions to ensure energy boundedness of the system.

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 derives conservative and skew-symmetric forms of the incompressible NS equations in sigma-coordinates by introducing new variables that restore energy conservation and bounds despite the metric terms.

read the letter

The key point is that the authors derive conservative and skew-symmetric formulations of the incompressible Navier-Stokes equations specifically for sigma-coordinates. By introducing a new set of variables, they restore the energy-conserving property for the Euler equations and the energy-bounded property for the full NS system, which standard transformations lose due to metric terms. They do this by building the conservative form directly from conservation laws instead of transforming the Cartesian version. This keeps consistency with the underlying physics. The eigenstructure analysis shows how the wave speeds and such differ from the Cartesian case, which is helpful for understanding the numerics. The skew-symmetric part uses the new variables to make the convective term skew-symmetric in the appropriate inner product, so its contribution to the energy vanishes or stays controlled. They also cover characteristic-based boundary conditions to maintain the bounds at the domain edges. This addresses a real issue in simulations over varying terrain, where long-term stability is important for applications like weather modeling or engineering flows. The derivations appear careful and avoid the usual pitfalls. The main soft spot is the reliance on the specific choice of new variables to cancel all metric and Jacobian contributions exactly in the energy inner product. For arbitrary terrain slopes, this cancellation might not be automatic, and any residual could break the bound. The abstract claims it works, but the full algebra needs checking to confirm no hidden assumptions about the sigma transformation. If the paper verifies it only for particular cases, that would limit the generality. Overall the thinking is clear and engages honestly with the structural properties literature. No circular reasoning or invented entities that don't check out. This is aimed at researchers in computational fluid dynamics who care about preserving invariants in transformed coordinates. A reader working on stable discretizations for geophysical flows would find the explicit forms and analysis worthwhile. It is worth sending to peer review because the problem is well-posed and the solution is targeted, even if the broader applicability needs checking in the details.

Referee Report

2 major / 2 minor

Summary. The manuscript derives conservative and skew-symmetric formulations of the incompressible Navier-Stokes equations in terrain-following sigma coordinates. It claims that, unlike direct sigma transformation of the Cartesian equations, the new forms preserve structural properties: the conservative version is constructed from general conservation laws with an analyzed eigenstructure, while a new set of variables yields a skew-symmetric convective operator that is energy-conserving for the Euler equations and energy-bounded for the Navier-Stokes equations. Characteristic-based boundary conditions are discussed to maintain energy boundedness.

Significance. If the derivations hold, the work provides a route to structure-preserving discretizations in sigma coordinates that retain exact energy conservation or boundedness without ad-hoc stabilization. This is relevant for CFD and geophysical flows over complex terrain, where conventional sigma formulations lose these properties due to metric terms. The explicit construction via new variables and the eigenstructure comparison are positive features that could support stable high-order schemes.

major comments (2)

[skew-symmetric formulation section] The energy conservation claim for the skew-symmetric form rests on the new variables ensuring that all metric and Jacobian contributions from the sigma transformation cancel exactly in the inner product with velocity. The manuscript must supply the explicit volume-integral calculation demonstrating this cancellation for arbitrary terrain slopes (i.e., non-constant metric coefficients), as any residual term would invalidate the energy estimate.
[conservative form derivation] The conservative form is stated to be 'consistent with general conservation laws,' yet the derivation steps that map the Cartesian divergence form to sigma coordinates while preserving the integral conservation property (including the modified eigenstructure) need to be shown in full detail; without this, it is unclear whether the form remains conservative for general sigma mappings.

minor comments (2)

Clarify the precise definition of the new set of variables introduced for the skew-symmetric form; a side-by-side comparison table with the original Cartesian variables would improve readability.
The abstract refers to 'incompressible flow equations' while the title specifies Navier-Stokes; ensure consistent terminology throughout regarding the treatment of the pressure term.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments on our manuscript. We address the major comments point by point below and will revise the manuscript to incorporate the suggested improvements.

read point-by-point responses

Referee: [skew-symmetric formulation section] The energy conservation claim for the skew-symmetric form rests on the new variables ensuring that all metric and Jacobian contributions from the sigma transformation cancel exactly in the inner product with velocity. The manuscript must supply the explicit volume-integral calculation demonstrating this cancellation for arbitrary terrain slopes (i.e., non-constant metric coefficients), as any residual term would invalidate the energy estimate.

Authors: We agree that providing the explicit calculation is necessary to fully substantiate the energy conservation property. Although the derivation in the manuscript is based on the new variables that lead to cancellation, the detailed volume-integral proof for non-constant metrics was not included. In the revised manuscript, we will add this explicit calculation, demonstrating the cancellation of all metric and Jacobian terms in the inner product for arbitrary terrain slopes. revision: yes
Referee: [conservative form derivation] The conservative form is stated to be 'consistent with general conservation laws,' yet the derivation steps that map the Cartesian divergence form to sigma coordinates while preserving the integral conservation property (including the modified eigenstructure) need to be shown in full detail; without this, it is unclear whether the form remains conservative for general sigma mappings.

Authors: We acknowledge that the derivation steps could be presented with greater detail. The conservative form is constructed by applying the transformation to the integral form of the conservation laws, which inherently preserves the conservation properties. We will revise the manuscript to include the complete derivation steps, explicitly showing the mapping from the Cartesian divergence form to the sigma-coordinate version and the analysis of the modified eigenstructure for general sigma mappings. revision: yes

Circularity Check

0 steps flagged

Derivation from conservation laws and variable redefinition is self-contained with no circular steps.

full rationale

The paper starts from general conservation laws to derive the conservative form in sigma-coordinates, then introduces a new set of variables to obtain the skew-symmetric convective operator. Energy conservation for Euler and boundedness for NS follow directly from the constructed inner-product properties after the transformation, without reducing to fitted inputs, self-definitions, or load-bearing self-citations. The abstract and derivation chain explicitly contrast the new forms against direct sigma-application, confirming the steps are independent constructions rather than tautological. This is the normal non-circular outcome for a structural derivation paper.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on standard assumptions of incompressible flow and coordinate transformation validity, plus introduction of a new variable set whose properties are asserted to deliver energy bounds.

axioms (2)

domain assumption The flow is incompressible
Standard assumption stated in the title and abstract for the Navier-Stokes equations.
domain assumption Sigma-coordinate transformation is valid and differentiable
Invoked implicitly when applying the transformation to derive the new forms.

invented entities (1)

New set of variables for the skew-symmetric formulation no independent evidence
purpose: To achieve energy-conserving and energy-bounded properties
Explicitly introduced in the abstract to derive the skew-symmetric form.

pith-pipeline@v0.9.0 · 5429 in / 1337 out tokens · 41113 ms · 2026-05-08T14:18:47.965123+00:00 · methodology

Conservative and skew-symmetric forms of the incompressible Navier-Stokes equations in sigma-coordinates

Core claim

Load-bearing premise

discussion (0)