Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?
Since proposed in Zhang and Shu (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.
First ; Corresponding
|ESI Research Field|
Cited Times [WOS]:0
|Document Type||Journal Article|
|Department||Department of Mathematics|
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
|First Author Affilication||Department of Mathematics; SUSTech International Center for Mathematics|
|Corresponding Author Affilication||Department of Mathematics; SUSTech International Center for Mathematics|
|First Author's First Affilication||Department of Mathematics|
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).
|Files in This Item:||There are no files associated with this item.|
|Recommend this item|
|Export to Endnote|
|Export to Excel|
|Export to Csv|
|Similar articles in Google Scholar|
|Similar articles in Baidu Scholar|
|Similar articles in Bing Scholar|
Items in the repository are protected by copyright, with all rights reserved, unless otherwise indicated.