Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?

Cui，Shumo1,2,4; Ding，Shengrong1; Wu，Kailiang1,2,3

2023-03-01

10.1016/j.jcp.2022.111882

JOURNAL OF COMPUTATIONAL PHYSICS   影响因子和分区

0021-9991

1090-2716

Since proposed in Zhang and Shu (2010) [1], 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.

Bound-preserving schemes Discontinuous Galerkin methods Hyperbolic conservation laws Multidimensional Optimal convex decomposition

SCI ; EI

英语

第一 ; 通讯

National Natural Science Foundation of China[12171227]

Computer Science ; Physics

Computer Science, Interdisciplinary Applications ; Physics, Mathematical

WOS:000990051900001

ACADEMIC PRESS INC ELSEVIER SCIENCE

20230213380341

Galerkin methods

Numerical Methods:921.6 ; Systems Science:961

PHYSICS

2-s2.0-85146049575

Scopus

1.Department of Mathematics,Southern University of Science and Technology,Shenzhen,518055,China
2.SUSTech International Center for Mathematics,Southern University of Science and Technology,Shenzhen,518055,China
3.National Center for Applied Mathematics Shenzhen (NCAMS),Shenzhen,518055,China
4.Guangdong Provincial Key Laboratory of Computational Science and Material Design,Shenzhen,518055,China

Cui，Shumo,Ding，Shengrong,Wu，Kailiang. Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?[J]. JOURNAL OF COMPUTATIONAL PHYSICS,2023,476.

Cui，Shumo,Ding，Shengrong,&Wu，Kailiang.(2023).Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?.JOURNAL OF COMPUTATIONAL PHYSICS,476.

Cui，Shumo,et al."Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?".JOURNAL OF COMPUTATIONAL PHYSICS 476(2023).