Drift-Free Conservative Dynamics from Quantized Interaction Rules

Myungjoo Kang; Park Junhu; Youngsoo Ha

arxiv: 2604.26383 · v2 · submitted 2026-04-29 · 🧮 math.NA · cs.NA· physics.comp-ph

Drift-Free Conservative Dynamics from Quantized Interaction Rules

Park Junhu , Youngsoo Ha , Myungjoo Kang This is my paper

Pith reviewed 2026-05-12 01:46 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-ph

keywords conservation lawsdiscrete dynamicsinteger transferquantized state spacedrift-freeentropy solutionsnumerical conservationmonotone schemes

0 comments

The pith

Conservation laws admit an exact discrete realization via antisymmetric integer-transfer operators on quantized states, enforcing invariance without round-off drift.

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

The paper replaces conventional floating-point flux evaluations for discretizing conservation laws with an exact interaction rule defined on a quantized state space. An antisymmetric integer-transfer operator performs the update, guaranteeing conservation at the arithmetic level rather than through approximate cancellation. For scalar laws the same operator uses monotone order-preserving transfers to select admissible shock structures inside the primitive step. Readers should care because this encodes both conservation and entropy selection directly at the operator level, removing a common source of numerical drift in long integrations and extending uniformly to systems and multiple dimensions.

Core claim

Conservation laws are realized as exact discrete interaction rules on a quantized state space. The update is performed by an antisymmetric integer-transfer operator that enforces conservation exactly at the arithmetic level and eliminates round-off drift. For scalar laws, monotone order-preserving transfers select admissible shock structures within the primitive update rather than through flux reconstruction. The construction extends to multidimensional problems and systems of conservation laws through oriented, vector-valued integer transfers.

What carries the argument

The antisymmetric integer-transfer operator, which exchanges integer quantities between neighboring states to enforce exact conservation and monotonicity at the arithmetic level.

If this is right

High-frequency transport is preserved near the Nyquist limit.
Sharply localized discontinuities are maintained in Burgers dynamics.
The same operator construction applies directly to multidimensional problems.
Systems of conservation laws are handled by oriented vector-valued integer transfers.
Entropy selection occurs inside the primitive update for scalar laws.

Where Pith is reading between the lines

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

The approach could be combined with existing finite-volume codes by swapping only the update kernel while retaining the same grid.
Long-time integrations of nearly incompressible flows might show reduced accumulation of conservation errors compared with floating-point methods.
The quantized representation might allow exact conservation even when using reduced-precision arithmetic hardware.

Load-bearing premise

An antisymmetric integer-transfer operator can be constructed that preserves monotonicity and selects admissible weak solutions for general conservation laws and systems without introducing new artifacts or losing accuracy.

What would settle it

Run a long-time integration of a linear advection or Burgers equation on a periodic domain and check whether the total integral remains exactly constant to machine precision in the new scheme while drifting in a standard floating-point flux scheme.

Figures

Figures reproduced from arXiv: 2604.26383 by Myungjoo Kang, Park Junhu, Youngsoo Ha.

**Figure 1.** Figure 1: FIG. 1. Relative error versus normalized frequency. FQNM view at source ↗

**Figure 2.** Figure 2: FIG. 2. Burgers shock comparison at identical grid resolution view at source ↗

**Figure 3.** Figure 3: FIG. 3. Magnified shock comparison at identical grid resolu view at source ↗

**Figure 4.** Figure 4: FIG. 4. Shu–Osher shock–entropy interaction. The integer view at source ↗

read the original abstract

Conservation laws are conventionally discretized through floating-point flux evaluation, with invariants obtained by cancellation of approximate interface contributions and admissible weak solutions selected by reconstruction and Riemann solvers. Here we introduce an operator-level formulation in which conservative dynamics is realized as an exact discrete interaction rule on a quantized state space. The update is defined by an antisymmetric integer-transfer operator, which enforces conservation exactly at the arithmetic level and eliminates round-off drift from the primitive evolution \cite{highamAccuracyStabilityNumerical2002}. For scalar laws, monotone order-preserving transfers select admissible shock structures within the primitive update, rather than through flux reconstruction. Numerical experiments show that the interaction rule preserves high-frequency transport near the Nyquist limit and maintains sharply localized discontinuities in Burgers dynamics. The same construction extends to multidimensional problems and systems of conservation laws through oriented, vector-valued integer transfers. These results indicate that conservative dynamics admits an exact discrete realization in which both invariance and entropy selection are encoded at the operator level, rather than arising from approximate flux cancellation.

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 gives a clean operator for exact arithmetic conservation via antisymmetric integer transfers on quantized states, with monotone selection working for scalars, but the systems extension is asserted without enough checks.

read the letter

This paper's main contribution is an antisymmetric integer-transfer operator on a quantized state space that builds exact conservation into the update rule itself. No floating-point flux cancellation is needed, so round-off drift disappears at the arithmetic level. For scalar conservation laws the monotone order-preserving transfers are meant to select admissible shocks directly inside the primitive step rather than through separate reconstruction or Riemann solvers. The Burgers tests show the scheme keeps high-frequency content near the Nyquist limit and holds discontinuities sharp, which is a solid concrete check on the scalar case. That part of the work is new and worth looking at if you care about long-time invariant preservation in discrete settings. The construction is self-contained and does not lean on circular self-citation. The soft spot is the extension to systems and multidimensional problems. The abstract states that oriented vector-valued transfers carry the same properties over, but the numerical evidence stays limited to the scalar one-dimensional Burgers equation. There is no demonstration yet that the vector version selects the correct entropy solution when several admissible combinations exist, nor any check that it avoids introducing oscillations or violating hyperbolicity. That claim therefore rests more on the operator definition than on verified behavior. The paper is aimed at numerical analysts who build discretizations for hyperbolic conservation laws and at people running long simulations where accumulated drift matters. A reader already working on monotone schemes or exact conservation methods will see the clearest value. It deserves a serious referee because the scalar results are concrete and the operator idea is distinct from standard flux-based approaches, even though the systems section will need more evidence before the central claim can be taken as general. I would send it for review with a request to add at least one system example and a short analysis of entropy selection in the vector case.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes an operator-level discretization of conservation laws via an antisymmetric integer-transfer operator acting on a quantized state space. This construction is claimed to enforce exact conservation at the arithmetic level (eliminating round-off drift) and, for scalar laws, to select admissible weak solutions through monotone order-preserving transfers rather than flux reconstruction or Riemann solvers. The approach is asserted to extend to multidimensional problems and systems of conservation laws using oriented vector-valued transfers, with numerical support cited for high-frequency transport and sharp discontinuities in the Burgers equation.

Significance. If the central claims hold, the work offers a genuine alternative paradigm for conservative discretizations in which both invariance and entropy admissibility are encoded directly in the interaction rule. This could be significant for long-time integrations where floating-point drift is problematic and for applications requiring strict conservation. The quantized, integer-transfer formulation is a clear strength if it proves robust beyond the scalar 1-D case shown.

major comments (2)

[Extension to systems] § on extension to systems (following the scalar construction): the central claim that oriented vector-valued integer transfers select admissible weak solutions for systems without new artifacts or loss of hyperbolicity is load-bearing, yet the manuscript provides no numerical experiments or analysis for any system (e.g., Euler or shallow-water equations). Only the scalar Burgers case is demonstrated, leaving the weakest assumption untested.
[Scalar laws] § on monotone transfers for scalar laws: while order-preserving integer transfers are stated to select admissible shocks within the primitive update, the manuscript does not supply a general proof that this holds for arbitrary scalar conservation laws (beyond the specific Burgers numerics); the selection property appears to rest on the operator definition rather than a verified entropy condition.

minor comments (2)

[Introduction / Abstract] The abstract cites Higham (2002) for round-off issues but does not discuss how the integer-transfer rule interacts with existing literature on exactly conservative schemes (e.g., finite-volume methods with exact flux cancellation).
[Numerical experiments] Figure captions for the Burgers experiments should explicitly state the grid size, time-stepping details, and comparison baselines (e.g., standard Godunov or WENO) to allow reproducibility assessment.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments identify key areas where additional support would strengthen the claims. We respond point by point below and indicate the revisions we will make.

read point-by-point responses

Referee: [Extension to systems] § on extension to systems (following the scalar construction): the central claim that oriented vector-valued integer transfers select admissible weak solutions for systems without new artifacts or loss of hyperbolicity is load-bearing, yet the manuscript provides no numerical experiments or analysis for any system (e.g., Euler or shallow-water equations). Only the scalar Burgers case is demonstrated, leaving the weakest assumption untested.

Authors: We agree that the manuscript currently demonstrates the approach only for the scalar Burgers equation and lacks numerical experiments for systems. The extension is formulated at the operator level using oriented vector-valued integer transfers that inherit antisymmetry for exact conservation and use orientation to align with characteristic directions. To address the concern, we will add a numerical example for a simple system (e.g., linear hyperbolic system or 1D shallow-water equations) in the revised manuscript. This will illustrate that the transfers introduce no new artifacts and preserve the expected hyperbolic behavior in the tested cases. revision: yes
Referee: [Scalar laws] § on monotone transfers for scalar laws: while order-preserving integer transfers are stated to select admissible shocks within the primitive update, the manuscript does not supply a general proof that this holds for arbitrary scalar conservation laws (beyond the specific Burgers numerics); the selection property appears to rest on the operator definition rather than a verified entropy condition.

Authors: The admissible-shock selection is a direct consequence of the monotone, order-preserving integer-transfer operator, which enforces the update without introducing non-physical oscillations or rarefactions. This property is verified numerically for the Burgers equation through preservation of sharp discontinuities and high-frequency transport. While the manuscript does not contain a general proof for arbitrary scalar laws, the operator construction mirrors the monotonicity that guarantees entropy satisfaction in the Kružkov sense. In the revision we will expand the discussion to clarify this link and supply additional justification tied to the operator definition. revision: partial

Circularity Check

0 steps flagged

No significant circularity; new operator introduced by definition

full rationale

The paper defines a new antisymmetric integer-transfer operator directly as the core of the discrete realization for conservation laws. The abstract and description present this as an original construction that enforces exact conservation and selects admissible solutions at the operator level, without reducing any central claim to a fitted parameter, prior self-citation, or renamed known result. The single external citation to Higham concerns round-off error and is not load-bearing for the derivation. No self-definitional, fitted-prediction, or uniqueness-imported patterns appear in the provided text.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The claim depends on the properties of this newly defined operator and the assumption that quantization does not alter the essential behavior.

axioms (1)

domain assumption The continuous conservation law dynamics can be faithfully represented on a quantized integer state space
Invoked to justify the discrete interaction rule as equivalent to the continuous case.

invented entities (1)

antisymmetric integer-transfer operator no independent evidence
purpose: Enforce exact conservation and select admissible solutions at the discrete level
Introduced as the core new construction in the paper.

pith-pipeline@v0.9.0 · 5473 in / 1185 out tokens · 66190 ms · 2026-05-12T01:46:31.761716+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

[1]

N. J. Higham,Accuracy and Stability of Numerical Algo- rithms, 2nd ed. (Society for Industrial and Applied Math- ematics, 2002)

work page 2002
[2]

E. F. Toro,Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction(Springer Sci- ence & Business Media, 2009)

work page 2009
[3]

R. J. LeVeque,Finite volume methods for hyperbolic problems, Vol. 31 (Cambridge university press, 2002)

work page 2002
[4]

Harten, Journal of Computational Physics135, 260 (1997)

A. Harten, Journal of Computational Physics135, 260 (1997)

work page 1997
[5]

Osher and S

S. Osher and S. Chakravarthy, SIAM Journal on Numer- ical Analysis21, 955 (1984)

work page 1984
[6]

Y. Ha, C. H. Kim, H. Yang, and J. Yoon, SIAM Journal on Scientific Computing38, A1987 (2016)

work page 2016
[7]

Y. Ha, C. H. Kim, H. Yang, and J. Yoon, SIAM Journal on Numerical Analysis59, 143 (2021)

work page 2021
[8]

C. E. Shannon, The Bell System Technical Journal27, 379 (1948)

work page 1948
[9]

T. M. Cover and J. A. Thomas,Elements of Information Theory(Wiley, 2006)

work page 2006
[10]

M. G. Crandall and L. Tartar, Proceedings of the Amer- ican Mathematical Society78, 385 (1980)

work page 1980
[11]

L. C. Evans,Partial Differential Equations, Vol. 19 (American Mathematical Society, 2010)

work page 2010

[1] [1]

N. J. Higham,Accuracy and Stability of Numerical Algo- rithms, 2nd ed. (Society for Industrial and Applied Math- ematics, 2002)

work page 2002

[2] [2]

E. F. Toro,Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction(Springer Sci- ence & Business Media, 2009)

work page 2009

[3] [3]

R. J. LeVeque,Finite volume methods for hyperbolic problems, Vol. 31 (Cambridge university press, 2002)

work page 2002

[4] [4]

Harten, Journal of Computational Physics135, 260 (1997)

A. Harten, Journal of Computational Physics135, 260 (1997)

work page 1997

[5] [5]

Osher and S

S. Osher and S. Chakravarthy, SIAM Journal on Numer- ical Analysis21, 955 (1984)

work page 1984

[6] [6]

Y. Ha, C. H. Kim, H. Yang, and J. Yoon, SIAM Journal on Scientific Computing38, A1987 (2016)

work page 2016

[7] [7]

Y. Ha, C. H. Kim, H. Yang, and J. Yoon, SIAM Journal on Numerical Analysis59, 143 (2021)

work page 2021

[8] [8]

C. E. Shannon, The Bell System Technical Journal27, 379 (1948)

work page 1948

[9] [9]

T. M. Cover and J. A. Thomas,Elements of Information Theory(Wiley, 2006)

work page 2006

[10] [10]

M. G. Crandall and L. Tartar, Proceedings of the Amer- ican Mathematical Society78, 385 (1980)

work page 1980

[11] [11]

L. C. Evans,Partial Differential Equations, Vol. 19 (American Mathematical Society, 2010)

work page 2010