“…In fact, the spurious oscillations may occur near the interfaces and wave fronts that are harmful to the robustness of the numerical algorithm. To overcome this difficulty, various schemes and approaches have been developed in the literature, such as interface tracking algorithms [20], diffusive kinetic schemes [3], relaxation schemes [10], finite difference/volume weighted essentially nonoscillatory (WENO) methods [2,24,27], entropy stable schemes with artificial viscosity [23], method of lines transpose (MOL T ) approach with nonlinear filters [12], discontinuous Galerkin (DG) methods with maximum-principle-satisfying limiters [36,41], local DG finite element methods [40], direct DG methods [30], etc. In this paper, we focus on the DG method and extend our previous work [31] to the nonlinear convection-diffusion problem (1.1).…”