pith. sign in

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

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

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

Pith reviewed 2026-05-17 21:02 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords PAMPA schemediscontinuous Galerkinbound preservationoscillation eliminationunstructured triangular mesheshyperbolic problemsnumerical methods
0
0 comments X

The pith

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.

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

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.

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.

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.

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.

Where Pith is reading between the lines

These are 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.

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.

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.

Figures

Figures reproduced from arXiv: 2511.16180 by R\'emi Abgrall, Yongle Liu.

Figure 1
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 φ = … view at source ↗
Figure 2
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. 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… view at source ↗
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. In the following, we select several representative and demanding examples, typical in the literature and sufficiently complex from our perspe… view at source ↗
Figure 5
Figure 5. Figure 5: Example 4: KT problem. Numerical solutions at t = 1 [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
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. Furthermore, [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
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 cell averages. However, the point values exhibit small remaining oscillations, indicating that the OE parameter is not sufficiently strong to f… view at source ↗
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.

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

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

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

specific steps
  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.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 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.

axioms (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.

pith-pipeline@v0.9.0 · 5496 in / 1340 out tokens · 45436 ms · 2026-05-17T21:02:52.677278+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

15 extracted references · 15 canonical work pages

  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 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

  9. [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,

    work page 2019

  10. [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–

    work page 2004

  11. [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,Viscous shock capturing in a time-explicit discontinuous Galerkin method, Math. Mode...

    work page arXiv 2020

  12. [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 inviscid and viscous flows, Comput. Fluids, 96 (2014), pp. 377–396. [59]J. S. Park a...

    work page 2012

  13. [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, arXiv preprint arXiv:2503.12473, (2025). [63]P.-O. Persson and J. Peraire,Sub-ce...

    work page arXiv 2025

  14. [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 dimensions, SIAM J. Sci. Comput., 47 (2025), pp. A3276–A3302. [66]P. Woodward and P. Cole...

    work page 1977

  15. [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 laws, J. Comput. Phys., 229 (2010), pp. 3091–3120. [73]X. Zhang and C.-W. Shu,On ...

    work page 2018