A posteriori subcell finite volume limiter for general PNPM schemes: applications from gasdynamics to relativistic magnetohydrodynamics

TL;DR

This paper introduces an a posteriori subcell finite volume limiter for high-order PNPM schemes, enhancing their robustness near discontinuities while preserving accuracy across various hyperbolic PDE applications.

Contribution

A new subcell finite volume limiting strategy is developed for PNPM schemes, effectively handling shocks and discontinuities without sacrificing high-order accuracy.

Findings

Effective shock capturing in gasdynamics, MHD, and relativistic MHD.

Maintains high-order accuracy in smooth regions.

Works on adaptive Cartesian meshes with diverse PDE systems.

Abstract

In this work, we consider the general family of the so called ADER PNPM schemes for the numerical solution of hyperbolic partial differential equations with \textit{arbitrary} high order of accuracy in space and time. The family of one-step PNPM schemes was introduced in [Dumbser et al., JCP, 2008] and represents a unified framework for classical high order Finite Volume (FV) schemes (N=0), the usual Discontinuous Galerkin (DG) methods (N=M), as well as a new class of intermediate hybrid schemes for which a reconstruction operator of degree M is applied over piecewise polynomial data of degree N with M>N. In all cases with M >= N > 0 the PNPM schemes are linear in the sense of Godunov, thus when considering phenomena characterized by discontinuities, spurious oscillations may appear and even destroy the simulation. Therefore, in this paper we present a new simple, robust and accurate a posteriori subcell finite volume limiting strategy that is valid for the entire class of PNPM schemes. The subcell FV limiter is activated only where it is needed, i.e. in the neighborhood of shocks or other discontinuities, and is able to maintain the resolution of the underlying high order PNPM schemes, due to the use of a rather fine subgrid of 2N+1 subcells per space dimension. The paper contains a wide set of test cases for different hyperbolic PDE systems, solved on adaptive Cartesian meshes (AMR) that show the capabilities of the proposed method both on smooth and discontinuous problems, as well as the broad range of its applicability. The tests range from compressible gasdynamics over classical MHD to relativistic magnetohydrodynamics.