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

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.