On Optimal Cell Average Decomposition for High-Order Bound-Preserving Schemes of Hyperbolic Conservation Laws

SIAM Journal on Numerical Analysis, Volume 62, Issue 2, Page 775-810, April 2024.
Abstract. Cell average decomposition (CAD) plays a critical role in constructing bound-preserving (BP) high-order discontinuous Galerkin and finite volume methods for hyperbolic conservation laws. Seeking optimal CAD (OCAD) that attains the mildest BP Courant–Friedrichs–Lewy (CFL) condition is a fundamentally important yet difficult problem. The classic CAD, proposed in 2010 by Zhang and Shu using the Gauss–Lobatto quadrature, has been widely used over the past decade. Zhang and Shu only checked for the 1D [math] and [math] spaces that their classic CAD is optimal. However, we recently discovered that the classic CAD is generally not optimal for the multidimensional [math] and [math] spaces. Yet, it remained unknown for a decade what CAD is optimal for higher-degree polynomial spaces, especially in multiple dimensions. This paper presents the first systematical analysis and establishes the general theory on the OCAD problem, which lays a foundation for designing more efficient BP schemes. The analysis is very nontrivial and involves novel techniques from several branches of mathematics, including Carathéodory’s theorem from convex geometry, and the invariant theory of symmetric group in abstract algebra. Most notably, we discover that the OCAD problem is closely related to polynomial optimization of a positive linear functional on the positive polynomial cone, thereby establishing four useful criteria for examining the optimality of a feasible CAD. Using the established theory, we rigorously prove that the classic CAD is optimal for general 1D [math] spaces and general 2D [math] spaces of an arbitrary [math]. For the widely used 2D [math] spaces, the classic CAD is, however, not optimal, and we develop a generic approach to find out the genuine OCAD and propose a more practical quasi-optimal CAD, both of which provide much milder BP CFL conditions than the classic CAD yet require much fewer nodes. These findings notably improve the efficiency of general high-order BP methods for a large class of hyperbolic equations while requiring only a minor adjustment of the implementation code. The notable advantages in efficiency are further confirmed by numerical results.