English

Solving the Stokes Equations via a Least Squares Weak Galerkin Method

Numerical Analysis 2026-07-12 v1

Abstract

We present a least-squares weak Galerkin (LS-WG) finite element method for solving the Stokes equations on arbitrary polygonal and polyhedral meshes. By utilizing discrete weak derivatives on discontinuous polynomial spaces, the proposed framework naturally accommodates complex domain geometries and general partitions. Crucially, this least-squares formulation bypasses the traditional inf-sup (LBB) compatibility condition, transforming the standard indefinite saddle-point problem into an inherently symmetric and positive definite (SPD) discrete linear system. We establish the well-posedness of the numerical scheme and rigorously derive optimal-order error estimates in a custom discrete energy norm. Specifically, we prove convergence rates of O(hk)\mathcal{O}(h^k) for the discrete projection error and O(hk1)\mathcal{O}(h^{k-1}) for the global approximation error when employing polynomials of degree k1k \ge 1 for the velocity field and k1k-1 for the pressure. Extensive numerical experiments confirm these theoretical convergence rates, demonstrating the method's robustness, geometric flexibility, and overall efficiency.

Cite

@article{arxiv.2607.10831,
  title  = {Solving the Stokes Equations via a Least Squares Weak Galerkin Method},
  author = {Chunmei Wang and Shangyou Zhang},
  journal= {arXiv preprint arXiv:2607.10831},
  year   = {2026}
}

Comments

20 pages, 6 tables, 5 figures