pith. sign in

arxiv: 2508.17147 · v3 · submitted 2025-08-23 · 🧮 math.NA · cs.NA

Some new properties of an Active flux type scheme: PamPa

R\'emi Abgrall , Philipp \"Offner , Yongle Liu This is my paper

Pith reviewed 2026-05-18 21:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords active fluxpampa schemediscontinuous galerkinbound preservingsummation by partsnumerical schemesconservation lawshigh-order methods
0
0 comments X

The pith

Active flux and PamPa schemes include the discontinuous Galerkin method as a building block while possessing intrinsic bound preservation and summation-by-parts properties.

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

The paper establishes that Active Flux and PamPa schemes admit a general interpretation in which the discontinuous Galerkin scheme functions as one of their constituent elements. It further shows that the present PamPa variant maintains solution bounds by design, without requiring separate limiters, and verifies this behavior through numerical examples. In one dimension the scheme satisfies the summation-by-parts identity, a property that supports discrete conservation and stability analysis. Readers working on high-order methods for hyperbolic conservation laws would care because these structural links suggest routes to combine existing techniques and to reduce reliance on ad-hoc fixes for maintaining physical solution ranges.

Core claim

In full generality the AF/pampa schemes can be interpreted such that the discontinuous Galerkin scheme is one of their building blocks; the current variant of pampa has intrinsic bound preserving properties; and at least in one dimension the pampa scheme has the summation by parts property.

What carries the argument

The general embedding of discontinuous Galerkin schemes inside Active Flux and PamPa formulations, together with the specific PamPa construction that directly enforces bound preservation and the one-dimensional summation-by-parts relation.

Load-bearing premise

The current specific variant of the pampa scheme admits the stated bound-preserving properties without additional limiting or post-processing steps, and the discontinuous Galerkin interpretation and summation-by-parts property extend directly from the scheme definition.

What would settle it

A direct calculation of the discrete inner-product matrices for the one-dimensional PamPa scheme that fails to satisfy the summation-by-parts identity, or a computed solution that exits the expected bounds in the absence of any limiter.

read the original abstract

In this paper, we provide a few new properties of Active Flux (AF)/Point-Average-Moment PolynomiAl-interpreted (\pampa) schemes. First, we show, in full generality, that the AF/pampa schemes can be interpreted in such a way that the discontinuous Galerkin (dG) scheme is one of their building blocks. Secondly we provide intrinsic bound preserving properties of the current variant of pampa. This is also illustrated numerically. Last, we show, at least in one dimension, that the pampa scheme has the summation by part (SBP) property.

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

0 major / 3 minor

Summary. The manuscript explores new properties of Active Flux (AF) and Point-Average-Moment PolynomiAl-interpreted (PamPa) schemes. It shows that AF/PamPa schemes admit an interpretation in which the discontinuous Galerkin (dG) scheme is one of the building blocks, establishes intrinsic bound-preserving properties for the current PamPa variant (with numerical illustrations), and proves that the PamPa scheme satisfies the summation-by-parts (SBP) property at least in one dimension.

Significance. If the derivations hold, the results supply a structural unification of AF/PamPa schemes with dG methods and confirm stability features (bound preservation via explicit algebraic cancellation under the stated CFL restriction and monotone flux, and SBP via direct summation of discrete operators) without extra limiting steps. These properties are load-bearing for the practical use of such schemes on conservation laws and constitute a clear contribution to the numerical analysis literature.

minor comments (3)
  1. [Introduction] The precise definition of the 'current variant' of PamPa (including the specific reconstruction and flux choices) should be stated explicitly at the beginning of the analysis sections so that the generality claims can be checked without cross-referencing earlier papers.
  2. [Numerical results] In the bound-preservation numerical examples, report the exact CFL numbers employed and include a short table or statement quantifying any observed maximum-principle violations (even if zero) to make the verification fully reproducible.
  3. [SBP property] The 1-D SBP identity is obtained by direct summation; a brief remark on whether the same cancellation pattern extends to the multi-dimensional case or requires additional assumptions would improve the discussion.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for recommending minor revision. The referee's summary accurately captures our main contributions: the general interpretation of AF/PamPa schemes with the discontinuous Galerkin method as a building block, the intrinsic bound-preserving properties of the current PamPa variant (with numerical illustrations), and the proof of the summation-by-parts property in one dimension. These results provide a structural unification and confirm stability features without additional limiting steps.

Circularity Check

0 steps flagged

No significant circularity; derivations are self-contained from scheme definitions

full rationale

The paper's central claims—the dG embedding via moment matching, intrinsic bound preservation via algebraic cancellation under the stated CFL and monotone flux, and the 1-D SBP identity via direct operator summation—are obtained by explicit algebraic and structural arguments applied to the PamPa scheme definition itself. No load-bearing step reduces to a fitted parameter renamed as prediction, a self-referential definition, or an unverified self-citation chain; each property is shown to follow directly from the discrete operators and update rules without hidden assumptions that presuppose the target result. The derivations remain internally verifiable from the given equations and hold independently of external fitted data.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, the work derives properties from the existing definition of Active Flux and PamPa schemes without introducing new free parameters, ad-hoc axioms, or invented entities.

pith-pipeline@v0.9.0 · 5628 in / 1114 out tokens · 46561 ms · 2026-05-18T21:19:50.085297+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    R. Abgrall , A combination of residual distribution and the active flux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d Euler equations , Commun. Appl. Math. Comput., 5 (2023), pp. 370– 402, https://doi.org/10.1007/s42967-021-00175-w

    work page doi:10.1007/s42967-021-00175-w 2023

  2. [2]

    R. Abgrall , A combination of residual distribution and the Active Flux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d Euler equations , Commun. Appl. Math. Comput., 5 (2023), pp. 370– 402

    work page 2023

  3. [3]

    Abgrall and W

    R. Abgrall and W. Barsukow , Extensions of active flux to arbitrary order of accuracy , ESAIM, Math. Model. Numer. Anal., 57 (2023), pp. 991–1027, https://doi.org/10.1051/ m2an/2023004

    work page arXiv 2023

  4. [4]

    Abgrall, M

    R. Abgrall, M. Jiao, Y. Liu, and K. Wu , Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme: one-dimensional case , Commun. Comput. Phys., (2025). In press

    work page 2025

  5. [5]

    Abgrall, J

    R. Abgrall, J. Lin, and Y. Liu , Active flux for triangular meshes for compressible flows problems, Beijing Journal of Pure and Applied Mathematics, (2025). in press, also Arxiv preprint 2312.11271

    work page arXiv 2025

  6. [6]

    Abgrall, J

    R. Abgrall, J. Lin, and Y. Liu , Active flux for triangular meshes for compressible flows problems, Beijing J. of Pure and Appl. Math., 2 (2025), pp. 1–33

    work page 2025

  7. [7]

    Abgrall and Y

    R. Abgrall and Y. Liu , A new approach for designing well-balanced schemes for the shallow water equations: a combination of conservative and primitive formulations , SIAM J. Sci. Comput., 46 (2024), pp. A3375–A3400

    work page 2024

  8. [8]

    Abgrall, Y

    R. Abgrall, Y. Liu, and W. Boscheri , Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) on polygonal meshes, 2024, https://arxiv.org/abs/2502. 10069, https://arxiv.org/abs/2502.10069

    work page arXiv 2024

  9. [9]

    Abgrall, J

    R. Abgrall, J. Nordstr ¨om, P. ¨Offner, and S. Tokareva, Analysis of the SBP-SAT stabi- lization for finite element methods. I: Linear problems , J. Sci. Comput., 85 (2020), p. 28, https://doi.org/10.1007/s10915-020-01349-z. Id/No 43

    work page doi:10.1007/s10915-020-01349-z 2020

  10. [10]

    Ahmed, A

    B. Ahmed, A. Alsaedi, F. Brezzi, L. Marini, and A. Russo , Projectors for Virtual Element Methods, Comput. Math. Appl., 66 (2013)

    work page 2013

  11. [11]

    Barsukow, The active flux scheme for nonlinear problems , J

    W. Barsukow, The active flux scheme for nonlinear problems , J. Sci. Comput., 86 (2021), pp. Paper No. 3, 34, https://doi.org/10.1007/s10915-020-01381-z, https://doi.org/10. 1007/s10915-020-01381-z

    work page doi:10.1007/s10915-020-01381-z 2021

  12. [12]

    Barsukow, P

    W. Barsukow, P. Chandrashekar, C. Klingenberg, and L. Lechner, A generalized active This manuscript is for review purposes only. 25 flux method of arbitrarily high order in two dimensions , 2025, https://arxiv.org/abs/2502. 05101, https://arxiv.org/abs/2502.05101

    work page arXiv 2025

  13. [13]

    Barsukow, J

    W. Barsukow, J. Hohm, C. Klingenberg, and P. L. Roe , The active flux scheme on Cartesian grids and its low Mach number limit , J. Sci. Comput., 81 (2019), pp. 594–622, https://doi.org/10.1007/s10915-019-01031-z

    work page doi:10.1007/s10915-019-01031-z 2019

  14. [14]

    Barsukow, C

    W. Barsukow, C. Klingenberg, L. Lechner, J. Nordstr¨om, S. Ortleb, and H. Ranocha, Stability of the active flux method in the framework of summation-by-parts operators, 2025, https://arxiv.org/abs/2507.11068, https://arxiv.org/abs/2507.11068

    work page arXiv 2025

  15. [15]

    Beir ˜ao da Veiga, F

    L. Beir ˜ao da Veiga, F. Brezzi, L. D. Marini, and A. Russo , The Hitchhiker’s guide to the virtual element method , Math. Models Methods Appl. Sci., 24 (2014), pp. 1541–1573, https://doi.org/10.1142/S021820251440003X

    work page doi:10.1142/s021820251440003x 2014

  16. [16]

    Calhoun, E

    D. Calhoun, E. Chudzik, and C. Helzel , The Cartesian grid Active Flux method with adaptive mesh refinement, J. Sci. Comput., 94 (2023), p. 54

    work page 2023

  17. [17]

    D. C. Del Rey Fern´andez, P. D. Boom, and D. W. Zingg, A generalized framework for nodal first derivative summation-by-parts operators, J. Comput. Phys., 266 (2014), pp. 214–239, https://doi.org/10.1016/j.jcp.2014.01.038

    work page doi:10.1016/j.jcp.2014.01.038 2014

  18. [18]

    D. C. Del Rey Fern ´andez, J. E. Hicken, and D. W. Zingg , Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations, Comput. Fluids, 95 (2014), pp. 171–196, https://doi.org/10.1016/j. compfluid.2014.02.016

    work page doi:10.1016/j 2014

  19. [19]

    Eyman, Active flux, PhD thesis, University of Michigan, 2013

    T. Eyman, Active flux, PhD thesis, University of Michigan, 2013

    work page 2013

  20. [20]

    Eyman and P

    T. Eyman and P. Roe , Active flux. 49th AIAA Aerospace Science Meeting, 2011

    work page 2011

  21. [21]

    Eyman and P

    T. Eyman and P. Roe , Active flux for systems . 20 th AIAA Computationa Fluid Dynamics Conference, 2011

    work page 2011

  22. [22]

    T. A. Eymann and P. L. Roe , Multidimensional Active Flux schemes , in 21st AIAA Compu- tational Fluid Dynamics Conference, A. I. of Aeronautics and Astronautics, eds., 2013

    work page 2013

  23. [23]

    Glaubitz, J

    J. Glaubitz, J. Nordstr¨om, and P. ¨Offner, Summation-by-parts operators for general func- tion spaces , SIAM J. Numer. Anal., 61 (2023), pp. 733–754, https://doi.org/10.1137/ 22M1470141

    work page 2023

  24. [24]

    Glaubitz, H

    J. Glaubitz, H. Ranocha, A. R. Winters, M. Schlottke-Lakemper, P. ¨Offner, and G. Gassner, Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods , J. Comput. Phys., 529 (2025), p. 26, https://doi. org/10.1016/j.jcp.2025.113841. Id/No 113841

    work page doi:10.1016/j.jcp.2025.113841 2025

  25. [25]

    F. He , Towards a New-generation Numerical Scheme for the Com- press- ible Navier-Stokes Equations with the Active Flux Method , PhD thesis, Applied and Interdisciplinary Mathematics, University of Michigan, 2021. https://deepblue.lib.umich.edu/handle/2027.42/169687

    work page 2021

  26. [26]

    Helzel, D

    C. Helzel, D. Kerkmann, and L. Scandurra , A new ADER method inspired by the active flux method, Journal of Scientific Computing, 80 (2019), pp. 35–61

    work page 2019

  27. [27]

    Liu and P

    X.-D. Liu and P. D. Lax , Solution of two-dimensional riemann problems of gas dynamics by positive schemes, SIAM J. Sci. Comput., 19 (1998), pp. 319–340

    work page 1998

  28. [28]

    J. Maeng , On the Advective Component of Active Flux Schemes for Nonlinear Hyperbolic Conservation Laws, PhD thesis, Applied and Interdisciplinary Mathematics, University of Michigan, 2017. https://deepblue.lib.umich.edu/handle/2027.42/138695

    work page 2017

  29. [29]

    Martinez, F

    A. Martinez, F. Piazzon, A. Sommariva, and M. Vianello , Quadrature-based polyno- mial optimization , Optim. Lett., 14 (2020), pp. 1027–1036, https://doi.org/10.1007/ s11590-019-01416-x, hdl.handle.net/11577/3292611

    work page 2020

  30. [30]

    Mattsson, Diagonal-norm upwind SBP operators , J

    K. Mattsson, Diagonal-norm upwind SBP operators , J. Comput. Phys., 335 (2017), pp. 283– 310, https://doi.org/10.1016/j.jcp.2017.01.042

    work page doi:10.1016/j.jcp.2017.01.042 2017

  31. [31]

    Olsson , Summation by parts, projections, and stability

    P. Olsson , Summation by parts, projections, and stability. I , Math. Comput., 64 (1995), pp. 1035–1065, s23–s26, https://doi.org/10.2307/2153482

    work page doi:10.2307/2153482 1995

  32. [32]

    Olsson , Summation by parts, projections, and stability

    P. Olsson , Summation by parts, projections, and stability. II , Math. Comput., 64 (1995), pp. 1473–1493, https://doi.org/10.2307/2153366

    work page doi:10.2307/2153366 1995

  33. [33]

    Sommariva and M

    A. Sommariva and M. Vianello , Gauss-green cubature and moment computation over arbi- trary geometries, J. Comput. Appl. Math., 231 (2009), pp. 886–896, https://doi.org/10. 1016/j.cam.2009.05.014

    work page 2009

  34. [34]

    Sommariva and M

    A. Sommariva and M. Vianello , Cheap and stable quadrature on polyhedral elements . Pre- print, arXiv:2502.03446 [math.NA] (2025), 2025, https://arxiv.org/abs/2502.03446

    work page arXiv 2025

  35. [35]

    Arzani, K

    M. Sv¨ard and J. Nordstr ¨om, Review of summation-by-parts schemes for initial-boundary- value problems , J. Comput. Phys., 268 (2014), pp. 17–38, https://doi.org/10.1016/j.jcp. 2014.02.031. Appendix A. VEM approximation: basic facts. This manuscript is for review purposes only. 26 R. ABGRALL, Y. LIU AND P. ¨OFFNER Any basis of VkpP q is virtual, meaning ...

    work page doi:10.1016/j.jcp 2014

  36. [36]

    wh is a polynomial of degree k on each edge e of P , that is pwhq|e P Pkpeq,

    work page

  37. [37]

    wh is continuous on BP ,

    work page

  38. [38]

    ş P whmα dx “ ş P π∇whmα dx for |α| “ k ´ 1, k

    ∆ wh P Pk´2pP q, 4. ş P whmα dx “ ş P π∇whmα dx for |α| “ k ´ 1, k. The degrees of freedom are the same as in VkpP q:

    work page

  39. [39]

    The value of wh on the vertices of P ,

    work page

  40. [40]

    On each edge of P , the value of vk at the k ´ 1 internal points of the k ` 1 Gauss–Lobatto points on this edge,

    work page

  41. [41]

    The L2 projection of wh is computable

    The moments up to order k ´ 2 of wh in P , mαpwhq :“ 1 |P | ż P whmα dx, |α| ď k ´ 2. The L2 projection of wh is computable. For any β, if π0pwhq “ ř α,|α|ďk sαmα, we have xπ0pwhq, mβy “ ÿ α,|α|ďk sαxmβ, mαy “ x wh, mβy. The left hand side is computable since the inner product xmβ, mαy only involves monomials. We need to look at the right hand side. If |α...

    work page