Abstract: Hyperbolic systems of conservation laws in multiple spatial dimensions display features absent in the one-dimensional case, such as involutions and non-trivial stationary states. These features need to be captured by numerical methods without excessive grid refinement. The active flux method is an extension of the finite volume scheme with additional point values distributed along the cell boundary. For the equations of linear acoustics, an exact evolution operator can be used for the update of these point values. It incorporates all multi-dimensional information. The active flux method is stationarity preserving, i.e., it discretizes all the stationary states of the PDE. This paper demonstrates the experimental evidence for the discrete stationary states of the active flux method and shows the evolution of setups towards a discrete stationary state.

Key words: Structure preserving, Stationarity preserving, Active flux