Virtual finite element and hyperbolic problems: The PAMPA algorithm

In this paper, we explore the use of the Virtual Element Method (VEM) concepts to solve scalar and system hyperbolic problems on general polygonal grids. The new schemes stem from the Active Flux approach [1], which combines the usage of point values at the element boundaries with an additional degree of freedom representing the average of the solution within each control volume. Along the lines of the family of residual distribution schemes introduced in [2, 3] that integrate the Active Flux technique, we devise novel third order accurate methods that rely on the VEM technology to discretize gradients of the numerical solution by means of a polynomial-free approximation, by adopting a virtual basis that is locally defined for each element. The obtained discretization is globally continuous, and for nonlinear problems it needs a stabilization which is provided by a monolithic convex limiting strategy extended from [4]. This is applied to both point and average values of the discrete solution. We show applications to scalar problems, as well as to the acoustics and Euler equations in two dimension. The accuracy and the robustness of the proposed schemes are assessed against a suite of benchmarks involving smooth solutions, shock waves and other discontinuities.