Mathematics and Computers in Simulation, cilt.246, ss.491-508, 2026 (SCI-Expanded, Scopus)
In this paper, we analyze the error of the weak Galerkin finite element method (WG-FEM) for singularly perturbed reaction–diffusion equations using piecewise discontinuous bilinear polynomials on a 2D Shishkin mesh. To accurately capture boundary layer effects, we introduce a weighted and balanced norm that is stronger than the standard energy norm. Unlike traditional approaches, which rely on global or local L2-projections and face challenges in 2D settings, our balanced norm is defined using a weight function together with the bilinear interpolation operator. By employing integral identities, we establish supercloseness results on Shishkin meshes. Numerical experiments confirm the sharpness of the theoretical analysis.