OEDG: OSCILLATION-ELIMINATING DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC CONSERVATION LAWS

Peng, Manting1; Sun, Zheng2; Wu, Kailiang1,3,4

2024-07-01

10.1090/mcom/3998

MATHEMATICS OF COMPUTATION   影响因子和分区

0025-5718

1088-6842

. Suppressing spurious oscillations is crucial for designing reliable high-order numerical schemes for hyperbolic conservation laws, yet it has been a challenge actively investigated over the past several decades. This paper proposes a novel, robust, and efficient oscillation-eliminating discontinuous Galerkin (OEDG) method on general meshes, motivated by the damping technique (see J. Lu, Y. Liu, and C. W. Shu [SIAM J. Numer. Anal. 59 (2021), pp. 1299-1324]). The OEDG method incorporates an oscillation-eliminating (OE) procedure after each Runge-Kutta stage, and it is devised by alternately evolving the conventional semidiscrete discontinuous Galerkin (DG) scheme and a damping equation. A novel damping operator is carefully designed to possess both scale-invariant and evolution-invariant properties. We rigorously prove the optimal error estimates of the fully discrete OEDG method for smooth solutions of linear scalar conservation laws. This might be the first generic fully discrete error estimate for nonlinear DG schemes with an automatic oscillation control mechanism. The OEDG method exhibits many notable advantages. It effectively eliminates spurious oscillations for challenging problems spanning various scales and wave speeds, without necessitating problem-specific parameters for all the tested cases. It also obviates the need for characteristic decomposition in hyperbolic systems. Furthermore, it retains the key properties of the conventional DG method, such as local conservation, optimal convergence rates, and superconvergence. Moreover, the OEDG method maintains stability under the normal Courant-Friedrichs-Lewy (CFL) condition, even in the presence of strong shocks associated with highly stiff damping terms. The OE procedure is nonintrusive, facilitating seamless integration into existing DG codes as an independent module. Its implementation is straightforward and efficient, involving only simple multiplications of modal coefficients by scalars. The OEDG approach provides new insights into the damping mechanism for oscillation control. It reveals the role of the damping operator as a modal filter, establishing close relations between the damping technique and spectral viscosity techniques. Extensive numerical results validate the theoretical analysis and confirm the effectiveness and advantages of the OEDG method.

Hyperbolic conservation laws discontinuous Galerkin method oscil- lation elimination modal filter scale invariance optimal error estimates

SCI

英语

第一 ; 通讯

Shenzhen Science and Technology Program[RCJC20221008092757098] ; National Natural Science Foundation of China[12171227] ; NSF[DMS-2208391]

Mathematics

Mathematics, Applied

WOS:001280563300001

AMER MATHEMATICAL SOC

MATHEMATICS

Web of Science

1.Southern Univ Sci & Technol, Dept Math, Shenzhen 518055, Guangdong, Peoples R China
2.Univ Alabama, Dept Math, Tuscaloosa, AL 35487 USA
3.Southern Univ Sci & Technol, Shenzhen Int Ctr Math, Shenzhen 518055, Guangdong, Peoples R China
4.Natl Ctr Appl Math Shenzhen NCAMS, Shenzhen 518055, Guangdong, Peoples R China

Peng, Manting,Sun, Zheng,Wu, Kailiang. OEDG: OSCILLATION-ELIMINATING DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC CONSERVATION LAWS[J]. MATHEMATICS OF COMPUTATION,2024.

Peng, Manting,Sun, Zheng,&Wu, Kailiang.(2024).OEDG: OSCILLATION-ELIMINATING DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC CONSERVATION LAWS.MATHEMATICS OF COMPUTATION.

Peng, Manting,et al."OEDG: OSCILLATION-ELIMINATING DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC CONSERVATION LAWS".MATHEMATICS OF COMPUTATION (2024).