Local time-integration for Friedrichs' systems

In this paper, we address the full discretization of Friedrichs' systems with a two-field structure, such as Maxwell's equations or the acoustic wave equation in div-grad form, cf. [14]. We focus on a discontinuous Galerkin space discretization applied to a locally refined mesh or a small region with high wave speed. This results in a stiff system of ordinary differential equations, where the stiffness is mainly caused by a small region of the spatial mesh. When using explicit time-integration schemes, the time step size is severely restricted by a few spatial elements, leading to a loss of efficiency. As a remedy, we propose and analyze a general leapfrog-based scheme which is motivated by [5]. The new, fully explicit, local time-integration method filters the stiff part of the system in such a way that its CFL condition is significantly weaker than that of the leapfrog scheme while its computational cost is only slightly larger. For this scheme, the filter function is a suitably scaled and shifted Chebyshev polynomial. While our main interest is in explicit local-time stepping schemes, the filter functions can be much more general, for instance, a certain rational function leads to the locally implicit method, proposed and analyzed in [24]. Our analysis provides sufficient conditions on the filter function to ensure full order of convergence in space and second order in time for the whole class of local time-integration schemes.