Auto-Stabilized Weak Galerkin Finite Element Methods for Stokes Equations on Non-Convex Polytopal Meshes

This paper introduces an auto-stabilized weak Galerkin (WG) finite element method for solving Stokes equations without relying on traditional stabilizers. The proposed WG method accommodates both convex and non-convex polytopal elements in finite element partitions, leveraging bubble functions as a key analytical tool. The simplified WG method is symmetric and positive definite, and optimal-order error estimates are derived for WG approximations in both the discrete $H^1$ norm and the $L^2$ norm.