arxiv: 2605.08105 · v1 · submitted 2026-04-27 · ⚛️ physics.comp-ph

Recognition: no theorem link

On Reconstructing Conservative and Primitive Variables: An Eigenvector Analysis on Curvilinear Grids

Amareshwara Sainadh Chamarthi

Pith reviewed 2026-05-12 00:44 UTC · model grok-4.3

classification ⚛️ physics.comp-ph

keywords conservative eigenvectorsprimitive eigenvectorscurvilinear gridscharacteristic reconstructioncontact discontinuitiesentropy correctionhypersonic boundary layers

0 comments

The pith

Conservative curvilinear eigenvectors have exact metric-free zeros in the density row of their shear columns, isolating contact discontinuities to the entropy eigenvector alone.

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

The paper establishes that standard conservative-variable eigenvectors on curvilinear grids place exact zeros, independent of the metrics, in the density components of the shear-wave columns of the right-eigenvector matrix. This property means shear waves carry no density perturbation and that a contact discontinuity is represented solely by the conservative entropy eigenvector. The corresponding left eigenvectors yield a metric-independent projection that extracts the entropy amplitude with unit contact scaling while giving zero projection of total-energy perturbations onto the shear amplitudes. By contrast, the standard primitive-variable eigenvectors introduce metric-dependent density components in the shear right eigenvectors and metric-weighted tangential-velocity terms in the entropy left eigenvector. These algebraic facts supply the two requirements for an exact, sufficient, metric-invariant rank-one entropy correction in characteristic reconstruction.

Core claim

For the standard conservative curvilinear eigenvectors, the density row of the right-eigenvector matrix contains exact, metric-free zeros in the shear columns, so shear waves carry no density perturbation and a contact discontinuity is represented by the conservative entropy eigenvector alone. The conservative left eigenvectors provide the dual projection property: the entropy amplitude is obtained with a metric-independent left eigenvector and has unit contact scaling, while total-energy perturbations have zero projection onto the shear amplitudes. In the standard primitive curvilinear eigenvectors, by contrast, shear right eigenvectors contain metric-dependent density components and the 5.

What carries the argument

Standard conservative curvilinear right and left eigenvectors, whose density row of the right matrix has metric-free zeros in the shear columns and whose left eigenvectors give metric-independent entropy projection with unit scaling.

If this is right

The conservative formulation supplies the algebraic requirements for an exact, sufficient, metric-invariant, rank-one entropy correction.
A contact discontinuity is represented by the conservative entropy eigenvector alone, without contribution from shear waves.
The entropy amplitude can be obtained with a metric-independent left eigenvector that has unit contact scaling.
Total-energy perturbations have zero projection onto the shear amplitudes in the conservative left-eigenvector basis.
The same conservative state-space contact direction underlies WA-CR even on Cartesian grids.

Where Pith is reading between the lines

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

The metric independence may simplify implementation of characteristic-based limiters in codes that already use conservative variables.
Primitive-variable reconstructions could be modified by projecting onto the conservative contact direction to recover similar invariance.
The result suggests that any variable set whose eigenvectors lack metric dependence in the contact direction will inherit the same reconstruction advantage.

Load-bearing premise

The analysis assumes the standard forms of conservative and primitive curvilinear eigenvectors as defined in the referenced literature and that the rank-one entropy correction premise applies directly to the contact direction isolated by these eigenvectors.

What would settle it

Explicit numerical evaluation of the right-eigenvector matrix for any non-Cartesian curvilinear metric and direct verification that the two shear-column entries in the density row are exactly zero.

read the original abstract

In wall-modelled large-eddy simulations of hypersonic boundary-layer transition, Hoffmann, Chamarthi and Frankel reported that characteristic reconstruction based on conservative-variable eigenvectors produced markedly better results than the corresponding primitive-variable implementation. The observation was empirical. A subsequent wave-appropriate conservative reconstruction (WA-CR) algorithm used a rank-one entropy correction based on the premise that contact-discontinuity error lies in a single conservative entropy/contact direction. This note gives the algebraic foundation for both observations. For the standard conservative curvilinear eigenvectors, the density row of the right-eigenvector matrix contains exact, metric-free zeros in the shear columns, so shear waves carry no density perturbation and a contact discontinuity is represented by the conservative entropy eigenvector alone. The conservative left eigenvectors provide the dual projection property: the entropy amplitude is obtained with a metric-independent left eigenvector and has unit contact scaling, while total-energy perturbations have zero projection onto the shear amplitudes. In the standard primitive curvilinear eigenvectors, by contrast, shear right eigenvectors contain metric-dependent density components and the primitive entropy left eigenvector contains metric-weighted tangential-velocity terms. Thus the conservative formulation supplies the two algebraic requirements for an exact, sufficient, metric-invariant, rank-one entropy correction: metric-independent entropy projection and a metric-independent entropy update direction. Curvilinear metrics make the distinction explicit, but the conservative state-space contact direction is already the natural direction underlying WA-CR even on Cartesian grids.

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

This note derives the exact algebraic reason conservative eigenvectors isolate contacts cleanly on curvilinear grids while primitive ones do not.

read the letter

The main thing to know is that the paper shows why conservative-variable characteristic reconstruction succeeds where primitive-variable versions fall short in curvilinear hypersonic simulations. For the standard conservative right-eigenvector matrix the density row has exact zeros in the shear columns with no metric dependence, so shear waves carry no density perturbation and contacts sit only in the entropy eigenvector. The left eigenvectors then give a metric-free projection that extracts the entropy amplitude with unit contact scaling and zero total-energy leakage onto shears. Primitive eigenvectors lack both properties because their shear columns pick up metric-dependent density terms and their entropy left vector mixes tangential velocities. This supplies the two requirements for an exact rank-one entropy correction without grid dependence. The algebra starts from the usual eigenvector definitions and extracts these features directly, which is new relative to the earlier empirical reports. The derivation is straightforward and internally consistent once the matrices are written out. It does not add fitting parameters or circular steps. The only soft spot is that the note stops at the algebra and does not include fresh numerical tests or code; it inherits whatever subtleties exist in the cited standard eigenvector forms. No load-bearing assumptions about grid orthogonality appear. This is for readers already working on high-order reconstruction schemes for compressible Navier-Stokes on curvilinear meshes, especially those doing wall-modeled LES for boundary-layer transition. A methods person who has seen the performance difference in practice will find the explanation useful. The work is clear and grounded enough to deserve referee time rather than a desk reject.

Referee Report

0 major / 2 minor

Summary. The manuscript provides an algebraic analysis of the right- and left-eigenvector matrices for the Euler equations in both conservative and primitive variables on curvilinear grids. It shows that the density row of the standard conservative right-eigenvector matrix contains exact, metric-independent zeros in the shear columns, implying that shear waves carry no density perturbation and that a contact discontinuity is represented solely by the conservative entropy eigenvector. The corresponding left eigenvectors yield a metric-independent entropy projection with unit contact scaling and zero projection of total-energy perturbations onto shear amplitudes. These properties are contrasted with the metric-dependent density components and tangential-velocity terms appearing in the primitive-variable eigenvectors. The analysis thereby supplies the algebraic justification for the empirical superiority of conservative-variable characteristic reconstruction and for the validity of the rank-one entropy correction used in the WA-CR algorithm, even on Cartesian grids.

Significance. If the algebraic claims hold, the note supplies a parameter-free derivation that directly explains the reported performance difference between conservative and primitive reconstructions in hypersonic wall-modelled LES and justifies the design of metric-invariant, rank-one entropy corrections. This strengthens the theoretical foundation for characteristic-based reconstruction schemes in curvilinear coordinates without introducing additional assumptions beyond the standard eigenvector forms.

minor comments (2)

The abstract refers to the empirical observation in Hoffmann, Chamarthi and Frankel; the manuscript should ensure the full bibliographic entry appears in the reference list with the correct year and venue.
Although the derivations follow from standard matrix forms, an appendix or short section displaying the explicit conservative and primitive eigenvector matrices (with metric terms highlighted) would allow immediate verification of the claimed zero entries without requiring the reader to consult the referenced literature.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of the manuscript and the recommendation to accept.

Circularity Check

0 steps flagged

No significant circularity; algebraic properties derived directly from standard eigenvector definitions

full rationale

The manuscript performs an eigenvector analysis starting from the standard conservative and primitive curvilinear eigenvector matrices as defined in the referenced literature (Hoffmann et al.). The key results—exact metric-free zeros in the density row of the right-eigenvector shear columns, metric-independent left-eigenvector entropy projection with unit contact scaling, and zero projection of total-energy perturbations onto shear amplitudes—are obtained by direct algebraic inspection of those matrices and the conservative-to-primitive transformation. No parameters are fitted, no quantities are predicted from a subset of the same data, and no load-bearing premise reduces to a self-citation whose validity is assumed rather than re-derived. The self-citation to prior empirical observations by the same author group supplies only historical context for the WA-CR algorithm; the present derivation is self-contained and independent of that prior work. This is the normal, non-circular case of a paper that extracts structural properties from established definitions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper relies on the established mathematical definitions of eigenvectors for the compressible Euler equations in both conservative and primitive variable sets; no new free parameters, ad-hoc axioms, or invented physical entities are introduced.

axioms (1)

standard math Standard definitions of left and right eigenvectors for the Euler equations in conservative and primitive variables on curvilinear grids.
The analysis invokes the known analytic forms of these matrices as given in the referenced literature on characteristic reconstruction.

pith-pipeline@v0.9.0 · 5557 in / 1418 out tokens · 57815 ms · 2026-05-12T00:44:48.666363+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

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

[1]

Van Leer, Upwind and high-resolution methods for compressible flow: From donor cell to residual-distribution schemes, in: 16th AIAA Computational Fluid Dynamics Conference, 2003, p

B. Van Leer, Upwind and high-resolution methods for compressible flow: From donor cell to residual-distribution schemes, in: 16th AIAA Computational Fluid Dynamics Conference, 2003, p. 3559. 10

work page 2003
[2]

A. S. Chamarthi, N. Hoffmann, S. Frankel, A wave appropriate discontinuity sensor approach for compressible flows, Physics of Fluids 35 (6) (2023)

work page 2023
[3]

Hoffmann, A

N. Hoffmann, A. S. Chamarthi, S. H. Frankel, Centralized gradient-based reconstruction for wall modelled large eddy simulations of hypersonic boundary layer transition, Journal of Computational Physics (2024) 113128

work page 2024
[4]

A. S. Chamarthi, Wave-appropriate multidimensional upwinding approach for compressible multiphase flows, Journal of Computational Physics 538 (2025) 114157

work page 2025
[5]

A. S. Chamarthi, Physics appropriate interface capturing reconstruction approach for viscous compressible multicomponent flows, Computers & Fluids 303 (2025) 106858

work page 2025
[6]

A. S. Chamarthi, Wave-appropriate reconstruction of compressible flows: physics- constrained acoustic dissipation and rank-1 entropy wave correction, arXiv preprint arXiv:2604.02757 (2026)

work page internal anchor Pith review Pith/arXiv arXiv 2026
[7]

Haselbacher and B

A. Haselbacher and B. Landmann,NSCBC on Curvilinear Grids, 2019

work page 2019
[8]

G. K. Batchelor, An introduction to fluid dynamics, Cambridge university press, 1967

work page 1967
[9]

J. C. Meng, T. Colonius, Numerical simulation of the aerobreakup of a water droplet, Journal of Fluid Mechanics 835 (2018) 1108–1135. 11

work page 2018