Computational and Applied Mathematics, cilt.43, sa.1, 2024 (SCI-Expanded, Scopus)
In this paper, we design and develop a weak Galerkin finite-element numerical method for solving singularly perturbed convection–diffusion–reaction equations with a non-conservation convection term. Many advantages of the proposed method include support for the higher order of convergence and general polygonal meshes. The convergence study for the weak Galerkin algorithm is performed in both the triple-bar norm and the L2 norm. We achieve an optimal order of convergence of O(hk) in the triple-bar norm and O(hk+1) in the L2 norm. Several numerical experiments in a two-dimensional setting are carried out to demonstrate the convergence of our theories.