On the Stability Barrier of Hermite Type Discretizations of Advection Equations

TL;DR

This paper investigates the stability limits of high-order Hermite-type discretizations for 1D advection equations, revealing a stability barrier related to stencil bias and extending the analysis to related Hermite methods.

Contribution

It establishes, both numerically and theoretically, a stability barrier for Hermite-type discretizations, a concept previously known for finite-difference schemes, and extends the analysis to Hermite WENO methods.

Findings

A stability barrier exists for Hermite-type discretizations with biased upwind stencils.

Semi-discretized HV methods become unstable when the stencil is sufficiently upwind-biased.

Tighter stability bounds are derived using combinatorial methods.

Abstract

In this paper we establish a stability barrier of a class of high-order Hermite-type discretization of 1D advection equations underlying the hybrid-variable (HV) and active flux (AF) methods. These methods seek numerical approximations to both cell-averages and nodal solutions and evolves them in time simultaneously. It was shown in earlier work that the HV methods are supraconvergent, providing that the discretization uses more unknowns in the upwind direction than the downwind one, similar to the "upwind condition" of classical finite-difference schemes. Although it is well known that the stencil of finite-difference methods could not be too biased towards the upwind direction for stability consideration, known as "stability barrier", such a barrier has not been established for Hermite-type methods. In this work, we first show by numerical evidence that a similar barrier exists for HV methods and make a conjecture on the sharp bound on the stencil. Next, we prove the existence of stability barrier by showing that the semi-discretized HV methods are unstable given a stencil sufficiently biased towards the upwind direction. Tighter barriers are then proved using combinatorical tools, and finally we extend the analysis to studying other Hermite-type methods built on approximating nodal solutions and derivatives, such as those widely used in Hermite WENO methods.