Pith. sign in

REVIEW 2 major objections 2 minor 1 cited by

Robust PAMPA Scheme in the DG Formulation on Unstructured Triangular Meshes: bound preservation, oscillation elimination, and boundary conditions

T0 review · 2 major / 2 minor · reviewed 2026-05-17 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗

Pith's one-line read The PAMPA scheme reinterpreted as a discontinuous Galerkin method on unstructured triangular meshes preserves bounds, eliminates oscillations, and achieves third-order accuracy for smooth solutions.

desk verdict This refines PAMPA into a globally continuous DG form on triangles, adds a non-oscillatory fix and cleaner BCs, but stays incremental on the authors' own prior results. read the letter →

arxiv 2511.16180 v3 submitted 2025-11-20 math.NA cs.NA

Rémi Abgrall , Yongle Liu This is my paper · ORCID
classification math.NAcs.NA
keywords PAMPAschemediscontinuousGalerkinboundpreservationoscillationeliminationunstructuredtriangularmesheshyperbolicproblemsnumericalmethods
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper develops an improved version of the PAMPA algorithm that seeks a globally continuous solution on unstructured triangular meshes. This formulation is locally conservative and does not require inverting a mass matrix. By reinterpreting PAMPA as a discontinuous Galerkin method for linear hyperbolic problems, the authors define a family of methods and handle boundary conditions rigorously. They complement prior bound-preserving techniques with a non-oscillatory approach. Truncation error analysis indicates third-order accuracy for smooth solutions, and this is verified through numerical experiments on a range of benchmarks where the scheme remains bound-preserving and free of oscillations.

What carries the argument

The reinterpretation of PAMPA as a discontinuous Galerkin method for the linear hyperbolic problem, which carries the definition of the method family and the rigorous implementation of boundary conditions.

What would settle it

A convergence test on a smooth solution where the measured order falls below three, or a benchmark run where the computed solution violates the preserved bounds or develops visible oscillations.

Watch

Extended reading notes

Core claim

The improved PAMPA scheme in the DG formulation seeks a globally continuous solution that is locally conservative without a mass matrix to invert. The reinterpretation as a discontinuous Galerkin method for linear hyperbolic problems defines a family of methods and enables rigorous boundary condition implementation. A complementary non-oscillatory method is introduced alongside bound preservation. Truncation error analysis shows the scheme is third-order accurate for smooth solutions, confirmed by numerical experiments demonstrating bound preservation and oscillation elimination across wide benchmarks on unstructured triangular meshes.

Load-bearing premise

The reinterpretation of PAMPA as a discontinuous Galerkin method for the linear hyperbolic problem allows defining the family of methods and implementing boundary conditions in a rigorous manner.

Editorial extensions

If this is right

  • The scheme remains locally conservative on unstructured triangular meshes.
  • No mass matrix inversion is needed during computation.
  • Third-order accuracy holds for smooth solutions as shown by truncation error analysis and tests.
  • The method stays bound-preserving and non-oscillatory across the presented numerical benchmarks.
  • Boundary conditions are applied consistently through the DG framework.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The non-oscillatory addition could be tested on problems with discontinuities such as contact waves or shocks.
  • The absence of a mass matrix may reduce computational cost relative to standard DG implementations.
  • The approach might be examined for extension to quadrilateral meshes or three-dimensional domains.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes an improved PAMPA scheme formulated as a globally continuous discontinuous Galerkin method on unstructured triangular meshes for hyperbolic problems. Leveraging the authors' prior reinterpretation of PAMPA as a DG method for linear hyperbolic problems, it defines a family of methods and implements boundary conditions rigorously. A non-oscillatory complement to the bound-preserving method is introduced, a truncation error analysis establishes third-order accuracy for smooth solutions, and numerical benchmarks demonstrate bound preservation and non-oscillatory behavior across a range of tests.

Significance. If the central claims hold, the work advances robust high-order schemes for hyperbolic conservation laws by providing a locally conservative, continuous formulation without mass-matrix inversion, with systematic boundary-condition treatment via the DG connection. The truncation error analysis and numerical validation are strengths that support the accuracy and robustness assertions on unstructured meshes; the combination of bound preservation with oscillation elimination addresses practical needs in simulations.

major comments (2)
  1. The reinterpretation as DG (building on Abgrall2025d) is used to claim rigorous boundary-condition implementation for the globally continuous variant, but the manuscript does not explicitly show how continuity is enforced (e.g., via zero-jump penalties or modified fluxes) while preserving the truncation-error order and bound properties without hidden mesh-quality or linearity assumptions; this is load-bearing for the central claim of rigorous BC treatment and third-order accuracy.
  2. Truncation error analysis section: the analysis indicates third-order accuracy for smooth solutions but does not address potential additional terms arising from global continuity enforcement or the non-oscillatory complement; explicit verification is needed to confirm the order is retained under the continuous formulation.
minor comments (2)
  1. Abstract: 'a wide range on numerical benchmarks' should read 'a wide range of numerical benchmarks'.
  2. Ensure all references to prior works (e.g., Abgrall2024a, Abgrall2025d) include complete bibliographic information and are clearly distinguished from the new contributions.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major comment below and have revised the manuscript to improve clarity on the points raised.

read point-by-point responses

  1. Referee: The reinterpretation as DG (building on Abgrall2025d) is used to claim rigorous boundary-condition implementation for the globally continuous variant, but the manuscript does not explicitly show how continuity is enforced (e.g., via zero-jump penalties or modified fluxes) while preserving the truncation-error order and bound properties without hidden mesh-quality or linearity assumptions; this is load-bearing for the central claim of rigorous BC treatment and third-order accuracy.

    Authors: We appreciate this observation. In the revised manuscript we have added an explicit description in Section 3 of how global continuity is enforced through the choice of a continuous finite-element space within the DG weak formulation. This choice sets the inter-element jumps to zero by construction, without requiring additional penalty terms or modified fluxes beyond those already present in the DG reinterpretation of Abgrall2025d. Boundary conditions are imposed via the standard numerical flux at the domain boundary. The truncation-error analysis and bound-preservation properties remain unchanged because the continuity constraint is satisfied exactly in the chosen space and does not introduce lower-order terms for smooth solutions. No hidden mesh-quality or linearity assumptions are added beyond those stated in the prior work. revision: yes

  2. Referee: Truncation error analysis section: the analysis indicates third-order accuracy for smooth solutions but does not address potential additional terms arising from global continuity enforcement or the non-oscillatory complement; explicit verification is needed to confirm the order is retained under the continuous formulation.

    Authors: We thank the referee for highlighting this gap. The truncation-error analysis in Section 4 is performed on the base scheme; because global continuity is enforced exactly by the continuous approximation space, it does not generate additional error terms for smooth solutions. The non-oscillatory complement is a limiter-type procedure that is inactive on smooth data and therefore does not affect the formal order. In the revised manuscript we have inserted a short paragraph in the analysis section that explicitly states these facts and cross-references the numerical convergence studies, which confirm third-order accuracy on the continuous formulation. revision: yes

Circularity Check

2 steps flagged · score 6.0 of 10

Core scheme definition, BC implementation, and bound preservation reduce to self-citations in Abgrall2025d and Abgrall2024a

  1. self citation load bearing [Abstract]
    "This method had been developed in a series of papers, see e.g. Abgrall2024a and the references therein. In Abgrall2025d, we had shown the connection between PAMPA and the discontinuous Galerkin method, for the linear hyperbolic problem. Taking advantage of this reinterpretation, we use it to define a family of methods, show how to implement the boundary conditions in a rigorous manner."

    The family of methods and rigorous boundary condition implementation are defined by directly adopting the DG reinterpretation from the authors' prior paper Abgrall2025d; the current work does not re-derive or independently justify this foundation but imports it via self-citation to establish the scheme.

  2. self citation load bearing [Abstract]
    "In addition, we propose a method that complements the bound preserving method developed in Abgrall2025d in the sense that it is non oscillatory."

    The bound-preserving property (central to the robustness claims of bound preservation on benchmarks) is imported wholesale from the authors' prior paper Abgrall2025d; this paper only adds a complementary non-oscillatory component rather than establishing bound preservation independently.

full rationale

The paper's central framework for the improved globally continuous PAMPA scheme is explicitly built by taking advantage of the DG reinterpretation and bound-preserving method from the authors' own prior work (Abgrall2025d), with the overall method history cited to Abgrall2024a. The truncation error analysis and non-oscillatory complement provide independent elements, but the load-bearing definitions of the method family, rigorous BCs, and bound preservation are imported without re-derivation here. This matches a moderate circularity level per the guidelines, as self-citations are load-bearing for the core claims while numerical benchmarks offer external checks.

Assumptions & free parameters 0 free parameters · 2 assumptions · 0 invented entities

The paper rests on the PAMPA-DG connection established in prior self-citations and standard assumptions for hyperbolic problems; no new free parameters or invented entities are explicitly introduced in the abstract.

assumptions (2)
  • domain assumption The connection between PAMPA and the discontinuous Galerkin method holds for the linear hyperbolic problem.
    Invoked to define the family of methods and boundary conditions.
  • standard math Standard truncation error analysis applies to the scheme on smooth solutions.
    Used to claim third-order accuracy.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Robust PAMPA Scheme in the DG Formulation on Unstructured Triangular Meshes: bound preservation, oscillation elimination, and boundary conditions." pith.science (2026). https://pith.science/paper/2511.16180

@misc{pith2026251116180,
  author       = {Pith},
  title        = {Pith review of: Robust PAMPA Scheme in the DG Formulation on Unstructured Triangular Meshes: bound preservation, oscillation elimination, and boundary conditions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2511.16180}},
  note         = {Machine review of arXiv:2511.16180}
}
read the original abstract

We propose an improved version of the PAMPA algorithm where the solution is sought as globally continuous. The scheme is locally conservative, and there is no mass matrix to invert. This method had been developed in a series of papers, see e.g \cite{Abgrall2024a} and the references therein. In \cite{Abgrall2025d}, we had shown the connection between PAMPA and the discontinuous Galerkin method, for the linear hyperbolic problem. Taking advantage of this reinterpretation, we use it to define a family of methods, show how to implement the boundary conditions in a rigorous manner. In addition, we propose a method that complements the bound preserving method developed in \cite{Abgrall2025d} in the sense that it is non oscillatory. A truncation error analysis is provided, it shows that the scheme should be third order accurate for smooth solutions. This is confirmed by numerical experiments. Several numerical examples are presented to show that the scheme is indeed bound preserving and non oscillatory on a wide range on numerical benchmarks.

Figures

Figures reproduced from arXiv: 2511.16180 by the authors.

Figure 1
Figure 1. DoFs and normals of triangle K. 2.1 Third-order PAMPA Scheme and Its DG Formulation In the standard third-order PAMPA (or the so-called semi-discrete/generalized AF) scheme, the solution of (1) is approximated by a globally continuous finite element polynomial expansion uh within each element K: uh|K = X 6 i=1 uσiφσi + uKφ, (3) where the basis functions are given by ([9, 11]) • For the average, we set φ = 60λ1λ2λ3, … view at source ↗
Figure 2
Figure 2. Geometry for the first-order scheme. with uT K i being the arithmetic average of the u’s at three vertices of T K i (hence we use the average value here) and αT K i ≥ max σ∈T K i , n TK i σi normals of T K i ρ [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Example 2: Zalesak problem. Numerical solutions at t = 1. Example 3—KPP Problem In the third example, we consider the so-defined classical KPP (Kurganov–Petrova–Popov) problem, which admits infinitely many weak solutions [46]: ∂u ∂t + ∂ sin(u) ∂x + ∂ cos(u) ∂y = 0, prescribed in a domain [−2, 2]2 with the initial condition u(x, 0) =    7 2 π if ∥x − (0, 0.5)∥ ≤ 1, π 4 else. The unique entropy solution at the fi… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Example 3: KPP problem. Numerical solutions at t = 1. 4.2 Euler Equations of Gas Dynamics A wide variety of standard and challenging numerical tests for compressible Euler equations are available to evaluate the robustness and accuracy of a proposed numerical scheme. I…
Figure 5
Figure 5. Figure 5: Example 4: KT problem. Numerical solutions at t = 1 [PITH_FULL_IMAGE:figures/full_fig_p016_5.png]
Figure 6
Figure 6. Figure 6: –(a) and (b) display the cell averages and point values of density, respectively. Compared with the results reported in [9], the BP OE PAMPA method produces noticeably fewer oscillations near shock waves, while remains more accurate profile around the double Mach stem.…
Figure 7
Figure 7. Figure 7: displays the cell averages and point values of density computed on a mesh with 39,728 elements and 78,897 points at time t = 0.35. With the BP OE blending parameter enabled, the scheme successfully captures the physical solution and the complex wave patterns in the cel…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Construction of entropy satisfying Active Flux-type methods

    math.NA 2026-07 conditional novelty 6.0 of 10

    Entropy stability of Active Flux-type schemes reduces to an inequality on cell averages alone; a Tadmor-style blend of high- and low-order fluxes enforces it while retaining bound preservation and oscillation control.

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages · cited by 1 Pith paper

  1. [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. [2]R. Abgrall and W. Barsukow,Extensions of Active Flux to arbitrary order of accuracy, ...

    work page doi:10.1051/m2an/2023004 2023

  2. [2]

    Abgrall, W

    [6]R. Abgrall, W. Boscheri, and Y. Liu,Virtual finite element and hyperbolic prob- lems: the PAMPA algorithm, J. Comput. Phys., (2026). Accepted. Available at: https://doi.org/10.48550/arXiv.2412.01341. [7]R. Abgrall, M. Jiao, Y. Liu, and K. Wu,A novel and simple invariant-domain-preserving frame- work for PAMPA scheme: 1D case, SIAM J. Sci. Comput., (202...

    work page doi:10.48550/arxiv.2412.01341 2026

  3. [3]

    Barsukow,Semi-discrete Active Flux as a Petrov-Galerkin method, arXiv preprint arXiv:2508.15017, (2025)

    [15]W. Barsukow,Semi-discrete Active Flux as a Petrov-Galerkin method, arXiv preprint arXiv:2508.15017, (2025). [16]W. Barsukow and J. Berberich,A well-balanced Active Flux method for the shallow water equations with wetting and drying, Commun. Appl. Math. Comput., 6 (2024), pp. 2385–2430. [17]W. Barsukow, J. Kern, C. Klingenberg, and L. Lechner,Analysis ...

    work page Pith review arXiv 2025

  4. [4]

    Chudzik, C

    [22]E. Chudzik, C. Helzel, and A. Porfetye,A fully discrete truly multidimensional Active Flux method for the two-dimensional Euler equations, arXiv preprint arXiv:2508.06273, (2025). [23]B. Cockburn, S. Hou, and C.-W. Shu,The Runge–Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case, Ma...

    work page arXiv 2025

  5. [5]

    Eymann and P

    19 [31]T. Eymann and P. Roe,Active Flux schemes for systems, in 20th AIAA computational fluid dynamics conference, 2011, p

    work page 2011

  6. [6]

    [33]T. A. Eymann and P. L. Roe,Multidimensional Active Flux schemes, in 21st AIAA computational fluid dynamics conference, 2013, p

    work page 2013

  7. [7]

    Fan and K

    [34]C. Fan and K. Wu,High-order oscillation-eliminating Hermite WENO method for hyperbolic conser- vation laws, J. Comput. Phys., 519 (2024), p. 113435. [35]D. Fan,On the acoustic component of active flux schemes for nonlinear hyperbolic conservation laws, PhD thesis, University of Michigan,

    work page 2024

  8. [8]

    [36]H. Fu, J. Xia, and X. Ma,Shock capturing with the high-order flux reconstruction method on adaptive meshes based on p4est, Engineering Reports, 4 (2022), p. e12501. [37]C. Geuzaine and J.-F. Remacle,Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities, Int. J. Numer. Meth. Eng., 79 (2009), pp. 1309–1331. [38]F. He...

    work page 2022

Show all 15 references
  1. [9]

    Helzel, D

    [39]C. Helzel, D. Kerkmann, and L. Scandurra,A new ADER method inspired by the Active Flux method, J. Sci. Comput., 80 (2019), pp. 1463–1497. [40]J. S. Hesthaven and T. Warburton,Nodal discontinuous Galerkin methods: algorithms, analysis, and applications, Springer,

    2019

  2. [10]

    Hoteit, P

    [41]H. Hoteit, P. Ackerer, R. Mos ´e, J. Erhel, and B. Philippe,New two-dimensional slope limiters for discontinuous Galerkin methods on arbitrary meshes, Int. J. Numer. Meth. Eng., 61 (2004), pp. 2566–

    2004

  3. [11]

    Huang and Y

    [42]J. Huang and Y. Cheng,An adaptive multiresolution discontinuous Galerkin method with artificial viscosity for scalar hyperbolic conservation laws in multidimensions, SIAM J. Sci. Comput., 42 (2020), pp. A2943–A2973. [43]A. Kl ¨ockner, T. Warburton, and J. S. Hesthaven,Visc...

    2020

  4. [12]

    [57]J. S. Park and C. Kim,Multi-dimensional limiting process for finite volume methods on unstructured grids, Comput. Fluids, 65 (2012), pp. 8–24. [58]J. S. Park and C. Kim,Higher-order multi-dimensional limiting strategy for discontinuous Galerkin methods in compressible invi...

    2012

  5. [13]

    [61]M. Peng, Z. Sun, and K. Wu,OEDG: Oscillation-eliminating discontinuous Galerkin method for hyperbolic conservation laws, Math. Comput., 94 (2025), pp. 1147–1198. [62]M. Peng, K. Wu, and C. Yuan,Oscillation-eliminating central DG schemes for hyperbolic conser- vation laws, ...

    2025

  6. [14]

    Van Leer,Towards the ultimate conservative difference scheme

    [64]B. Van Leer,Towards the ultimate conservative difference scheme. IV. a new approach to numerical convection, J. Comput. Phys., 23 (1977), pp. 276–299. [65]G. Wissocq, Y. Liu, and R. Abgrall,A positive-and bound-preserving vectorial lattice Boltzmann method in two dimension...

    1977

  7. [15]

    You and C

    [71]H. You and C. Kim,High-order multi-dimensional limiting strategy with subcell resolution I. Two- dimensional mixed meshes, J. Comput. Phys., 375 (2018), pp. 1005–1032. [72]X. Zhang and C.-W. Shu,On maximum-principle-satisfying high order schemes for scalar conserva- tion l...

    2018

Pith tools

Reviewed May 17, 2026 · model on record in the stance chip above.