A Novel and Simple Invariant-Domain-Preserving Framework for PAMPA Scheme: 1D Case

TL;DR

This paper introduces a new invariant-domain-preserving PAMPA scheme for hyperbolic conservation laws, ensuring solutions stay within physical bounds without additional limiters, and demonstrates its effectiveness on various 1D problems.

Contribution

The paper proposes a simple, provably invariant-domain-preserving PAMPA scheme with an automatic reformulation inspired by machine learning functions, eliminating the need for post-processing limiters.

Findings

The scheme guarantees invariant domain preservation for cell averages.

It effectively captures strong shocks with suppressed spurious oscillations.

Numerical experiments confirm high accuracy and robustness across multiple 1D problems.

Abstract

The PAMPA (Point-Average-Moment PolynomiAl-interpreted) method, proposed in [R. Abgrall, Commun. Appl. Math. Comput., 5: 370-402, 2023], combines conservative and non-conservative formulations of hyperbolic conservation laws to evolve cell averages and point values. Solutions to hyperbolic conservation laws typically have an invariant domain, and ensuring numerical solutions stay within this domain is essential yet nontrivial. This paper presents a novel framework for designing efficient Invariant-Domain-Preserving (IDP) PAMPA schemes. We first analyze the IDP property for updated cell averages in the original PAMPA scheme, revealing the role of cell average decomposition and midpoint values in maintaining the invariant domain. This analysis highlights the difficulty of relying on continuous fluxes alone to preserve the invariant domain. Building on these insights, we introduce a simple IDP limiter for cell midpoint values, and propose a provably IDP PAMPA scheme that guarantees the preservation of the invariant domain for updated cell averages without requiring post-processing limiters. This approach contrasts with existing bound-preserving PAMPA schemes, which often require additional convex limiting to blend high-order and low-order solutions. Most notably, inspired by the Softplus and Clipped ReLU functions from machine learning, we propose an automatic IDP reformulation of the governing equations, resulting in an unconditionally limiter-free IDP scheme for evolving point values. We also introduce techniques to suppress spurious oscillations, enabling the scheme to capture strong shocks effectively. Numerical experiments on 1D problems, including the linear convection equation, Burgers equation, the compressible Euler equations, and MHD equations, demonstrate the accuracy and robustness of the proposed IDP PAMPA scheme.