Abstract:In this paper, we propose and analyze a numerically stable and convergent scheme for a convection-diffusion-reaction equation in the convection-dominated regime. Discontinuous Galerkin (DG) methods are considered since standard finite element methods for the convection-dominated equation cause spurious oscillations. We choose to follow a novel DG finite element differential calculus framework introduced in Feng et al. (2016) and approximate the infinite-dimensional operators in the equation with the finite-dimensional DG differential operators. Specifically, we construct the numerical method by using the dual-wind discontinuous Galerkin (DWDG) formulation for the diffusive term and the average discrete gradient operator for the convective term along with standard DG stabilization. We prove that the method converges optimally in the convection-dominated regime. Numerical results are provided to support the theoretical findings.
What problem does this paper attempt to address?
This paper explores a new method of using the Discontinuous Galerkin (DG) method to deal with convection-dominated problems. In this study, the authors propose a numerically stable and convergent scheme to solve the convection-diffusion-reaction equation under convection-dominant conditions. The standard finite element method often exhibits false oscillations when dealing with such problems, especially when the convection term dominates. To address this issue, the paper adopts the Dual-Wind Discontinuous Galerkin (DWDG) method to handle the diffusion term, and uses the averaged discrete gradient operator to handle the convection term, combined with standard DG stability techniques to avoid oscillations and achieve upwind methods.
The paper first introduces the framework of the Discontinuous Galerkin differential operator and defines the relevant discrete operators. Then, the authors analyze a simplified problem, where the diffusion coefficient ε approaches 0 in the convection-diffusion term, and propose a numerical approximation method, proving the consistency and optimal convergence of this method. Next, the paper extends these methods to the complete convection-diffusion-reaction equation and provides a fully discretized numerical method, proving sublinear convergence under convection-dominant conditions.
Through numerical experiments, the authors demonstrate the effectiveness of the method and support the theoretical findings. The experimental results show that the proposed method performs well in the convection-dominated region, capturing boundary layers and internal layers, and maintains optimal convergence rates under different conditions.
In summary, this paper proposes a novel DG method suitable for convection-dominated problems, which effectively addresses the challenges posed by the convection term and avoids oscillation issues encountered in standard finite element methods. This method holds potential application value for convection-dominated problems in numerical simulations of fluid dynamics and chemical reactions, among other fields.