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.

Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme: one-dimensional case

Bound preserving Point-Average-Moment PolynomiAl-interpreted (PAMPA) scheme: one-dimensional case

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

Well-balanced Point-Average-Moment PolynomiAl-interpreted (PAMPA) Methods for Shallow Water Equations on Triangular Meshes

Conservative Piece-wise Parabolic Advective Remapping Method

Enhancement of the Material Point Method Using B‐spline Basis Functions

Bound-Preserving Framework for Central-Upwind Schemes for General Hyperbolic Conservation Laws

Arbitrary High-Order Non-Oscillatory Scheme on Hybrid Unstructured Grids Based on Multi-Moment Finite Volume Method.

Preliminary Effort in Developing the Smoothed Material Point Method for Impact

An Improved Quadrature Scheme in B-spline Material Point Method for Large-Deformation Problem Analysis

A high-order pseudo arc-length method with positivity-preserving flux limiter for compressible multi-medium flows

A Bound-Preserving Numerical Scheme for Space–Time Fractional Advection Equations

Circumventing volumetric locking in explicit material point methods: A simple, efficient, and general approach

A boundary point interpolation method for acoustic problems with particular emphasis on the calculation of Cauchy principal value integrals

A-ULMPM: An Arbitrary Updated Lagrangian Material Point Method for Efficient Simulation of Solids and Fluids

Virtual finite element and hyperbolic problems: the PAMPA algorithm

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws

Point-Based POMDP Algorithms: Improved Analysis and Implementation

High-Order Bound-Preserving Local Discontinuous Galerkin Methods for Incompressible and Immiscible Two-Phase Flows in Porous Media

A stable poro-mechanical formulation for Material Point Methods leveraging overlapping meshes and multi-field ghost penalisation

High-order -bound-preserving Fourier pseudo-spectral schemes for the Allen-Cahn equation