Primitive variable regularization to derive novel Hyperbolic Shallow Water Moment Equations
Pith reviewed 2026-05-22 01:08 UTC · model grok-4.3
The pith
Regularizing shallow water moment equations in primitive variables produces hyperbolic models that preserve momentum and admit analytical steady states.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Performing hyperbolic regularization in the primitive variables rather than the convective variables allows derivation of new shallow water moment equation systems that are hyperbolic, preserve the momentum equation, and possess analytically computable steady states, as demonstrated through variable transformations and dam-break simulations.
What carries the argument
Hyperbolic regularization performed in primitive variables (velocity and height) of the shallow water moment equations, which enables the proofs of hyperbolicity and steady states via transformations to the convective form.
If this is right
- The new models maintain hyperbolicity for stable numerical solutions.
- Preserving the momentum equation ensures physical accuracy in simulations.
- Analytical steady states provide interpretable equilibrium solutions for verification.
- Dam-break tests demonstrate superior accuracy compared to previous models.
Where Pith is reading between the lines
- This regularization strategy could be applied to derive improved moment models for other types of flows with vertical structure.
- Such models might lead to better long-term simulations in hydrological applications by avoiding unphysical oscillations.
- Further analysis could reveal how the choice of variables affects conservation properties in general hyperbolic systems.
Load-bearing premise
The vertical velocity expansion combined with primitive-variable regularization maintains physical fidelity without adding unphysical artifacts.
What would settle it
A dam-break simulation where the new model fails to conserve momentum or produces results less accurate than existing models would falsify the accuracy claim.
Figures
read the original abstract
Shallow Water Moment Equations are reduced-order models for free-surface flows that employ a vertical velocity expansion and derive additional so-called moment equations for the expansion coefficients. Among desirable analytical properties for such systems of equations are hyperbolicity, accuracy, correct momentum equation, and interpretable steady states. In this paper, we show analytically that existing models fail at different of these properties and we derive new models overcoming the disadvantages. This is made possible by performing a hyperbolic regularization not in the convective variables (as done in the existing models) but in the primitive variables. Via analytical transformations between the convective and primitive system, we can prove hyperbolicity and compute analytical steady states of the new models. Simulating a dam-break test case, we demonstrate the accuracy of the new models and show that it is essential for accuracy to preserve the momentum equation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives novel hyperbolic shallow water moment equations by applying hyperbolic regularization in primitive variables rather than the conventional convective variables used in prior models. It analytically demonstrates that existing models suffer from failures in hyperbolicity, momentum equation preservation, or interpretable steady states. New models are obtained via this primitive-variable approach, with hyperbolicity and steady states established through analytical transformations between the convective and primitive formulations. A dam-break test case is simulated to illustrate improved accuracy, with the claim that preserving the momentum equation is essential.
Significance. If the central derivations and transformation equivalences hold, the work offers a systematic route to shallow-water moment models that simultaneously satisfy hyperbolicity, correct momentum balance, and explicit steady states. This could improve the reliability of reduced-order models for free-surface flows in numerical analysis and engineering applications, particularly where analytical structure is required for stability proofs or long-time behavior.
major comments (2)
- [Abstract and derivation of new models] Abstract and the section presenting the regularization and transformations: the claim that primitive-variable regularization followed by analytical transformation yields a hyperbolic convective system assumes the regularization operator commutes with the change-of-variables map. Regularization alters flux or source terms in the primitive system; the inverse map back to convective variables is not guaranteed to preserve the same Jacobian eigenvalues or steady-state balance without explicit verification. Please supply the regularized primitive equations, the transformed convective form, and the characteristic analysis (e.g., the flux Jacobian and its eigenvalues) to confirm the hyperbolicity proof survives the regularization step.
- [Numerical results] Dam-break test and accuracy claims: while the simulation is presented to demonstrate accuracy when the momentum equation is preserved, the manuscript does not report quantitative error norms, grid-convergence studies, or comparisons against the underlying full shallow-water system or other moment closures. This leaves open whether the primitive regularization introduces new unphysical artifacts that offset the claimed gains in hyperbolicity and steady-state fidelity.
minor comments (2)
- [Introduction] Clarify the precise vertical velocity expansion and the moment closure assumptions early in the manuscript, as these enter the weakest assumption underlying physical accuracy.
- [Derivation] Ensure all transformed equations are written out explicitly rather than described only by reference to prior work, to allow readers to follow the analytical steps without external lookup.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable feedback on our manuscript. We address the major comments point by point below, and indicate the revisions we plan to make.
read point-by-point responses
-
Referee: [Abstract and derivation of new models] Abstract and the section presenting the regularization and transformations: the claim that primitive-variable regularization followed by analytical transformation yields a hyperbolic convective system assumes the regularization operator commutes with the change-of-variables map. Regularization alters flux or source terms in the primitive system; the inverse map back to convective variables is not guaranteed to preserve the same Jacobian eigenvalues or steady-state balance without explicit verification. Please supply the regularized primitive equations, the transformed convective form, and the characteristic analysis (e.g., the flux Jacobian and its eigenvalues) to confirm the hyperbolicity proof survives the regularization step.
Authors: We appreciate the referee's careful scrutiny of the derivation. While the manuscript establishes hyperbolicity through analytical transformations between the primitive and convective formulations, we acknowledge that explicit verification of the regularization step and the resulting Jacobian would enhance clarity. In the revised manuscript, we will supply the regularized primitive equations, detail the transformation to the convective form, and present the characteristic analysis including the flux Jacobian and its eigenvalues to rigorously confirm that hyperbolicity is preserved. revision: yes
-
Referee: [Numerical results] Dam-break test and accuracy claims: while the simulation is presented to demonstrate accuracy when the momentum equation is preserved, the manuscript does not report quantitative error norms, grid-convergence studies, or comparisons against the underlying full shallow-water system or other moment closures. This leaves open whether the primitive regularization introduces new unphysical artifacts that offset the claimed gains in hyperbolicity and steady-state fidelity.
Authors: The referee makes a valid observation regarding the numerical section. Our current dam-break test provides a qualitative demonstration of improved accuracy. To address this, we will augment the numerical results with quantitative error norms, a grid-convergence study, and direct comparisons to the full shallow-water system as well as other existing moment models. This will provide stronger evidence for the accuracy claims and help identify any potential artifacts introduced by the regularization. revision: yes
Circularity Check
No circularity: derivation uses independent analytical transformations
full rationale
The paper's central derivation performs hyperbolic regularization directly in primitive variables, followed by explicit analytical transformations to the convective form. These steps are self-contained mathematical operations that do not reduce any claimed prediction or property (hyperbolicity, steady states, momentum preservation) to a fitted input or prior self-citation by construction. Dam-break simulations supply external numerical evidence. No load-bearing claim collapses to an input by definition or renaming.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard assumptions of shallow-water theory and vertical velocity expansion remain valid for the moment equations
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
performing a hyperbolic regularization not in the convective variables ... but in the primitive variables. Via analytical transformations ... prove hyperbolicity and compute analytical steady states
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
-
[1]
M. Bauerle, A. J. Christlieb, M. Ding, and J. Huang. On the Rotational Invariance and Hyperbolicity of Shallow Water Moment Equations in Two Dimensions.SIAM Journal on Mathematical Analysis, 57(1):1039–1085, 2025. 4, 8, 9
work page 2025
- [2]
-
[3]
Z. Cai, Y. Fan, and R. Li. Globally hyperbolic regularization of grad’s moment system in one dimensional space.Communications in Mathematical Sciences, 11(2):547–571,
-
[4]
Z. Cai, Y. Fan, and R. Li. Globally Hyperbolic Regularization of Grad’s Moment System.Communications on Pure and Applied Mathematics, 67(3):464–518, 2014. 7
work page 2014
-
[5]
Z. Cai, R. Li, and Y. Wang. Numerical Regularized Moment Method for High{M}ach Number Flow.Commun. Comput. Phys., 11(5):1415–1438, 2012. 7
work page 2012
-
[6]
Y. Cao, A. Kurganov, Y. Liu, and R. Xin. Flux Globalization Based Well-Balanced Path-Conservative Central-Upwind Schemes for Shallow Water Models.Journal of Scientific Computing, 92(2):69, 2022. 2
work page 2022
-
[7]
M. J. Castro and C. Par´ es. Well-balanced high-order finite volume methods for systems of balance laws.Journal of Scientific Computing, 82(2):48, 2 2020. 29
work page 2020
-
[8]
M. Ciallella and J. Koellermeier. High order global flux schemes for general steady state preservation of shallow water moment equations with non-conservative products. Computers & Fluids, 305:106887, 2026. 30
work page 2026
-
[9]
Y. Fan, J. Koellermeier, J. Li, R. Li, and M. Torrilhon. Model Reduction of Kinetic Equations by Operator Projection.Journal of Statistical Physics, 162(2):457–486,
-
[10]
J. Garres-D´ ıaz, M. J. C. D´ ıaz, J. Koellermeier, and T. M. de Luna. Shallow Wa- ter Moment models for bedload transport problems.Commun. in Comput. Phys., 30(3):903–941, 2021. 1
work page 2021
-
[11]
G. J. Gassner, A. R. Winters, and D. A. Kopriva. A well balanced and entropy conservative discontinuous Galerkin spectral element method for the shallow water equations.Applied Mathematics and Computation, 272:291–308, 2016. 31
work page 2016
-
[12]
W.-X. Huai, S. Li, G. G. Katul, M.-Y. Liu, and Z.-H. Yang. Flow dynamics and sed- iment transport in vegetated rivers: A review.Journal of Hydrodynamics, 33(3):400– 420, 2021. 1
work page 2021
- [13]
-
[14]
S. Jin. Asymptotic preserving (AP) schemes for multiscale kinetic and hyperbolic equations: a review.Rivista di Matematica della Universit` a di Parma. New Series, 2, 2012. 31 36
work page 2012
-
[15]
Koellermeier.Derivation and numerical solution of hyperbolic moment equations for rarefied gas flows
J. Koellermeier.Derivation and numerical solution of hyperbolic moment equations for rarefied gas flows. PhD thesis, RWTH Aachen University, 2017. 2
work page 2017
-
[16]
J. Koellermeier and E. Pimentel-Garc´ ıa. Steady states and well-balanced schemes for shallow water moment equations with topography.Applied Mathematics and Computation, 427:127166, 2022. 2, 7, 10, 20, 21, 22, 26, 27, 28, 29
work page 2022
-
[17]
J. Koellermeier and M. Rominger. Analysis and Numerical Simulation of Hyper- bolic Shallow Water Moment Equations.Communications in Computational Physics, 28(3):1038–1084, 2020. 2, 7, 8, 10, 12, 13, 14, 16, 17, 20, 27, 28
work page 2020
-
[18]
J. Koellermeier and M. Torrilhon. On new hyperbolic moment models for the Boltz- mann equation. InConference Proceedings of the YIC GACM 2015, 2015. 15, 22
work page 2015
-
[19]
J. Kowalski and M. Torrilhon. Moment approximations and model cascades for shal- low flow.Commun. Comput. Phys., 25(3):669–702, 2019. 1, 2, 4, 15, 27, 28
work page 2019
-
[20]
J. Mac´ ıas, M. J. Castro, S. Ortega, and J. M. Gonz´ alez-Vida. Performance assessment of Tsunami-HySEA model for NTHMP tsunami currents benchmarking. Field cases. Ocean Modelling, 152:101645, 2020. 1
work page 2020
-
[21]
J. McDonald and M. Torrilhon. Affordable robust moment closures for{CFD}based on the maximum-entropy hierarchy.Journal of Computational Physics, 251:500–523,
-
[22]
R. P. Schaerer and M. Torrilhon. On Singular Closures for the 5-Moment System in Kinetic Gas Theory.Communications in Computational Physics, 17(2):371–400,
-
[23]
R. Verbiest and J. Koellermeier. Capturing Vertical Information in Radially Sym- metric Flow Using Hyperbolic Shallow Water Moment Equations.Communications in Computational Physics, 37(3):810–848, 2025. 4
work page 2025
-
[24]
Q. Xu, Y. Shi, J. L. Bamber, C. Ouyang, and X. X. Zhu. Large-scale flood modeling and forecasting with FloodCast.Water Research, 264:122162, 2024. 1 37
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.