Is the Classic Convex Decomposition Optimal for Bound-Preserving Schemes in Multiple Dimensions?
From MaRDI portal
Publication:6405370
DOI10.1016/J.JCP.2022.111882arXiv2207.08849MaRDI QIDQ6405370FDOQ6405370
Authors: Shumo Cui, Shengrong Ding, Kailiang Wu 
Publication date: 18 July 2022
Abstract: Since proposed in [X. Zhang and C.-W. Shu, J. Comput. Phys., 229: 3091--3120, 2010], the Zhang--Shu framework has attracted extensive attention and motivated many bound-preserving (BP) high-order discontinuous Galerkin and finite volume schemes for various hyperbolic equations. A key ingredient in the framework is the decomposition of the cell averages of the numerical solution into a convex combination of the solution values at certain quadrature points, which helps to rewrite high-order schemes as convex combinations of formally first-order schemes. The classic convex decomposition originally proposed by Zhang and Shu has been widely used over the past decade. It was verified, only for the 1D quadratic and cubic polynomial spaces, that the classic decomposition is optimal in the sense of achieving the mildest BP CFL condition. Yet, it remained unclear whether the classic decomposition is optimal in multiple dimensions. In this paper, we find that the classic multidimensional decomposition based on the tensor product of Gauss--Lobatto and Gauss quadratures is generally not optimal, and we discover a novel alternative decomposition for the 2D and 3D polynomial spaces of total degree up to 2 and 3, respectively, on Cartesian meshes. Our new decomposition allows a larger BP time step-size than the classic one, and moreover, it is rigorously proved to be optimal to attain the mildest BP CFL condition, yet requires much fewer nodes. The discovery of such an optimal convex decomposition is highly nontrivial yet meaningful, as it may lead to an improvement of high-order BP schemes for a large class of hyperbolic or convection-dominated equations, at the cost of only a slight and local modification to the implementation code. Several numerical examples are provided to further validate the advantages of using our optimal decomposition over the classic one in terms of efficiency.
Hyperbolic equations and hyperbolic systems (35Lxx) Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems (65Mxx) Basic methods in fluid mechanics (76Mxx)
This page was built for publication: Is the Classic Convex Decomposition Optimal for Bound-Preserving Schemes in Multiple Dimensions?
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6405370)
